Original article

Explicit solutions to the time-fractional generalized dissipative Kawahara equation

  • Marwan Alquran ,
  • Mohammed Ali ,
  • Omar Alshboul
Expand
  • Department of Mathematics and Statistics, Jordan University of Science and Technology, Irbid 22110, Jordan

Received 31 August 2021, Revised 28 January 2022, Accepted 16 February 2022, Available online 19 February 2022, Version of Record 22 July 2024.E-mail:marwan04@just.edu.jo;myali@just.edu.jo;oaalshboul19@sci.just.edu.jo

Received date: 2021-08-31

  Revised date: 2022-01-28

  Accepted date: 2022-02-16

  Online published: 2024-08-13

Abstract

In this paper, the generalized dissipative Kawahara equation in the sense of conformable fractional derivative is presented and solved by applying the tanh-coth-expansion and sine-cosine function techniques. The quadratic-case and cubic-case are investigated for the proposed model. Expected solutions are obtained with highlighting to the effect of the presence of the alternative fractional-derivative and the effect of the added dissipation term to the generalized Kawahara equation. Some graphical analysis is presented to support the findings of the paper. Finally, we believe that the obtained results in this work will be important and valuable in nonlinear sciences and ocean engineering.

Highlights

● The conformable-time-fractional generalized dissipative Kawahara equation is presented and solved by using the tanh-coth-expansion and sine-cosine-function methods.

● The quadratic-case and cubic-case are investigated for the proposed model.

● The new solutions highlighted the effect of the presence of the alternative fractional-derivative and the effect of the added dissipation term to the generalized Kawahara equation.

● Some graphical analysis is presented to support the findings of the current work.

Cite this article

Marwan Alquran , Mohammed Ali , Omar Alshboul . Explicit solutions to the time-fractional generalized dissipative Kawahara equation[J]. Journal of Ocean Engineering and Science, 2024 , 9(4) : 348 -352 . DOI: 10.1016/j.joes.2022.02.013

1. Introduction

Finding exact solutions to nonlinear models is very important part in applied mathematical sciences due to their significant role in investigating and explaining the mechanisms of propagation nonlinear waves in different complex media such as in oceans, water-waves, physics, optics, and much more [1], [2], [3], [4]. In order to obtain a clear view of the dynamic predictions of many natural phenomena and link them to their historical records, the mathematical representations of these nonlinear models were revisited and the time-varying rates were considered as fractional derivatives [5], [6], [7], [8], [9], [10]. The fractional-derivative is a generalization to the integer-derivative, and it possesses inherited memory. As a motivation by this trend, we study a time-fractional version to a model that arising in the field of ocean engineering and science. The suggested model is the generalized Kawahara equation (KE).
The KE is the modified version of the fifth-order KdV, and it describes the dynamics of shallow water waves and capillary gravity water waves. This equation takes the form Kawahara [11]
ψt(x,t)+bm(ψm(x,t)x)+Aψxxx(x,t)Bψxxxxx(x,t)=0.
In the literature, the cosh-function method is suggested in [12] as a solution to (1) for the case m=1,2 . In [13] the tanh-method and exp-method are used to solve KE for the case of m=1,2 and b=1 . Our aim is to study the solution of KE with dissipative term and conformable-time derivative. In fact, the proposed model is
Dtαψ(x,t)+aψx(x,t)+bm(ψm(x,t))x+Aψxxx(x,t)Bψxxxxx(x,t)=0,0<α1,
where A,B are the dispersion parameters, a is the dissipative coefficient, b is the nonlinearity-strength, and Dtα is the conformable-derivative which is defined by Khalil et al. [14] as
Dtαψ(x,t)=limϵ0ψ(x,t+ϵt1α)ψ(x,t)ϵ.
We refer to (3) as the conformable-time dissipative KE.
Since the fractional derivative plays the role of inheriting some historical properties of a natural phenomenon or a physical system, researchers revisited many well-known nonlinear models with the presence of conformable-derivative. The dynamics of these fractional problems have been investigated through extracting their explicit solutions. For this purposes, many effective techniques are used, such as, different versions of the (G/G) -expansion scheme [15], [16], [17], modified exp-expansion scheme [18], extended Kudryashov method [19], [20], [21], Polynomial function scheme [22], [23], [24], Jacobi elliptic function method [25], [26], Sine-Gordon equation [27], [28], F-expansion and simple equation algorithms [29], and many others.
The contribution of this research is that the model given in (2) is studied for the first time. In fact, we get explicit solutions to it by using two effective methods valid for this type of equations. Moreover, we study the propagation of the obtained solitary solutions to the fractional dissipative Kawahara upon varying the values of the conformable-derivative.
The paper is organised as follows: In Section 2, we suggested the tanh-coth expansion technique [30], [31] and the sine-cosine function method [32], [33] to study the solution of (2) when the nonlinearity term is of order m=2 . In Section 3, the same format is given to solve (2) for m=3 . Finally, we summarized our findings in Section 4.

2. Exact solutions to the fractional dissipative KE: m=2

The conformable-time dissipative KE with m=2 takes the following form
Dtαψ(x,t)+aψx(x,t)+bψ(x,t)ψx(x,t)+Aψxxx(x,t)Bψxxxxx(x,t)=0.
By using the following relation [14]
Dtα(ψ(ϕ(t)))=t1αDtϕ(t)Dtψ(ϕ(t)),
the new wave-variable ζ=xctαα transform (4) into the following differential equation
(ca)ψ(ζ)+b2ψ2(ζ)+Aψ(ζ)Bψ(4)=0.
We aim to implement two methods to retrieve explicit solutions to (6).

2.1. Tanh-coth expansion method

The tanh-coth expansion method suggests that the solution of (6) is
ψ=ψ(ζ)=i=0nλiZi,
where Z=Z(ζ) is the solution of the auxiliary differential equation
Z=μ(1Z2),
whose solution is
Z=tanh(μζ),
or
Z=coth(μζ).
We apply the balance procedure on the terms ψ2 against ψ(4) to get the relation 2n=n+4 . Thus, the index n given in (7) has the value of 4. Accordingly, the solution of (6) is
w(ζ)=λ0+λ1Z+λ2Z2+λ3Z3++λ4Z4.
Next, we substitute (11) and its necessary derivatives in (6) to get a polynomial of eighth-order in the variable Z with the following coefficients:
Z0:λ0(ac)+2λ2μ2(A+8Bμ2)+bλ02224Bλ4μ4,Z1:λ1(a2Aμ2+bλ016Bμ4c)+6λ3μ2(A+20Bμ2),Z2:λ2(a8Aμ2+bλ0136Bμ4c)+12λ4μ2(A+40Bμ2)+bλ122,Z7:λ3(bλ4360Bμ4),Z8:bλ422840Bλ4μ4.
Equating the coefficient of Zi:i=0,1,,8 to zero, results in a non algebraic system with the unknowns λ0,λ1,λ2,λ3,λ4,μ and c. Solving the obtained system gives
λ0=105A2169bB,λ1=λ3=0,λ2=210A2169bB,λ4=1680Bμ4b,μ=A213B,c=169aB+36A2169B.
Accordingly, two solutions to (4) are obtained and given by the following:
ψ1(x,t)=105A2169bB210A2tanh2(A(xtα(169aB+36A2)169αB)213B)169bB+105A2tanh4(A(xtα(169aB+36A2)169αB)213B)169bB,
and
ψ2(x,t)=105A2169bB210A2coth2(A(xtα(169aB+36A2)169αB)213B)169bB+105A2coth4(A(xtα(169aB+36A2)169αB)213B)169bB.
Figure 1, presents the 3D solutions of (14) and (15) for the integer case α=1 . In Figs. 2 and 3, we study the effect of the conformable-derivative acting on the propagation of ψ1(x,t) and ψ2(x,t) . It can be observed that for α>0.2 , the solutions are mapping continuously to each other with the same physical shape. However, for α<0.2 , we report a chaotic behaviors of the solutions where we can tell nothing [34]. Based on these observations, we may say that the fractional derivative can be considered as a memory index in a way that after a specific value, correct heritage-information can be retrieved regarding the physical structures of the given model.
Fig.1 3D soliton and singular-soliton solutions of (14) and (15), respectively, for the integer case α=1 and the assigned parameters a=b=A=B=1 .
Fig.2 The effect of the conformable-derivative acting on the propagation of ψ1(x,t) (14) when t=0.5 and a=b=A=B=1 .
Fig.3 The effect of the conformable-derivative acting on the propagation of ψ2(x,t) (15) when t=0.5 and a=b=A=B=1 .

2.2. Sine-cosine function method

The sine-cosine function method suggests that the solution of (6) is
ψ(ζ)=λsinn(μζ),
or
ψ(ζ)=λcosn(μζ).
Then, we substitute (16) in (6) to obtain
0=2Bμ4n(n36n2+11n6)+2μ2(n1)nsin2(μz)(A+2Bμ2(n22n+2))+2sin4(μz)(a+μ2(n2)(A+Bμ2n2)c)+bλsinn+4(μz).
Now, we equate the exponent n+4 to zero, and we set each coefficient of sini(μz):i=0,2,4 to zero. Then, we solving the resulting system to have the following:
λ=105A2169bB,μ=iA213B,c=169aB+36A2169B.
Accordingly, the solution of (4) is
ψ3(x,t)=105A2csch4(A(xtα(169aB+36A2)169αB)213B)169bB.
On the other side, if (17) is used, one more solution will be obtained and given by
ψ4(x,t)=105A2sech4(A(xtα(169aB+36A2)169αB)213B)169bB.
Remark: By using the identities csch2(μz)=coth2(μz)1 and sech2(μz)=1tanh2(μz) , the solutions of (4) given in (20) and (21) are identical, respectively, to (15) and (14).

3. Exact solutions to the fractional dissipative KE: m=3

The conformable-time dissipative KE with m=3 takes the following form:
Dtαψ(x,t)+aψx(x,t)+bψ(x,t)2ψx(x,t)+Aψxxx(x,t)Bψxxxxx(x,t)=0.
The wave-variable ζ=xctαα transforms (22) into the following differential equation
(ca)ψ(ζ)+b3ψ3(ζ)+Aψ(ζ)Bψ(4)=0.

3.1. Tanh-coth expansion method

Balancing the order of the term ψ3(ζ) against ψ(4)(ζ) gives n=2 . Thus, the tanh-coth solution of (23) is
ψ(ζ)=a0+λ1Z+λ2Z2.
Implicit differentiations of (24) leads to
ψ(ζ)=2μ2(Z21)(λ2(3Z21)+λ1Z),ψ(4)(ζ)=8μ4(Z21)(λ1Z(3Z22)+λ2(15Z415Z2+2)).
By substitution of (24)-(25) in (23), we get a polynomial of sixth-order in the variable Z with the following coefficients:
Z0:λ0(ac)+2λ2μ2(A+8Bμ2)+bλ033,Z1:λ1(a2Aμ2+bλ0216Bμ4c),Z2:λ2(a8μ2(A+17Bμ2)c)+bλ2λ02+bλ12λ0,Z5:λ1(bλ2224Bμ4),Z6:bλ233120Bλ2μ4.
Equating the coefficient of Zi:i=0,1,,6 to zero, results in a non algebraic system with the unknowns λ0,λ1,λ2,μ , and c. Solving the obtained system gives
λ0=±3A10bB,λ1=0,λ2=610Bμ2b,μ=A25B,c=25aB+4A225B.
Accordingly, two solutions to the quadratic Kawahara equation are obtained and given by the following:
ψ5(x,t)=(3Atanh2(A(xtα(25aB+4A2)25αB)25B)10bB3A10bB),
and
ψ6(x,t)=(3Acoth2(A(xtα(25aB+4A2)25αB)25B)10bB3A10bB).
Figure 4, presents the 3D solutions of (28) and (29) for the integer case α=1 . In Figs. 5 and 6, we study the effect of the conformable-derivative acting on the propagation of ψ5(x,t) and ψ6(x,t). We observe the same findings as depicted in Figs. 2 and 3.
Fig.4 3D soliton and singular-soliton solutions of (28) and (29), respectively, for the integer case α=1 and the assigned parameters a=b=A=B=1 .
Fig.5 The effect of the conformable-derivative acting on the propagation of ψ5(x,t) (28) when t=0.5 and a=b=A=B=1 .
Fig.6 The effect of the conformable-derivative acting on the propagation of ψ6(x,t) (29) when t=0.5 and a=b=A=B=1 .

3.2. Sine-cosine function method

In this section, we skip the details of implementing the sine-cosine function method in solving (22), and we only present its findings:
n=2,λ=3A10bB,μ=±iA25B,c=25aB+4A225B.
As a result, the following two solutions are reported to the conformable-time cKE:
ψ7(x,t)=(3Acsch2(A(xtα(25aB+4A2)25αB)25B)10bB),
and
ψ8(x,t)=(3Asech2(A(xtα(25aB+4A2)25αB)25B)10bB).
We should point here that ψ7(x,t) and ψ8(x,t) are equivalent, respectively, to ψ6(x,t) and ψ5(x,t) .

4. Conclusions

The conformable-fractional dissipative generalized Kawahara equation is introduced for the first time in the paper. We applied two numerical schemes to extract explicit solitary solutions to the proposed model for two states; the quadratic and the cubic cases. Moreover, we provided a graphical analysis that visualizes the effect of the value of the fractional derivative on the propagation of the obtained solutions. In addition, we observed that the solutions admit the bifurcation or the chaotic behaviours around a specific value of the fractional order. The findings of the current work confirm the fact that the fractional derivative acts as a tool for retrieving some properties inherited in a particular scientific model.
As a future work, we recommend the following actions: (1) Include other types of fractional derivatives and compare it with the obtained results of this work. (2) Apply analytical schemes and recover approximate solutions. (3) Study other fractional nonlinear models arise in ocean engineering and science.
Declaration of Competing Interest
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

[1]
M. Alquran. Opt. Quant. Electron., 53 (588) (2021)

[2]
R.W. Ibrahim, C. Meshram, S.B. Hadid, S. Momani. J. Ocean Eng. Sci., 5 (2) (2020), pp. 186-195

[3]
A. Yusuf, T.A. Sulaiman, A. Abdeljabbar, M. Alquran. J. Ocean Eng. Sci. (2021),

[4]
L. Akinyemi, M. Senol, U. Akpan, K. Oluwasegun. Opt. Quant. Electron., 53 (394) (2021)

[5]
M. Parto-Haghighi, J. Manafian. J. Ocean Eng. Sci., 5 (3) (2020), pp. 197-204

[6]
M.A. Shallal, K.K. Ali, K.R. Raslan, H. Rezazadeh, A. Bekir. J. Ocean Eng. Sci., 5 (3) (2020), pp. 223-229

[7]
P. Veeresha, D.G. Prakasha, N. Magesh, A.J. Christopher, D.U. Sarwe. J. Ocean Eng. Sci., 6 (3) (2021), pp. 265-275

[8]
M. Alquran, F. Yousef, F. Alquran, T.A. Sulaiman, A. Yusuf. Math. Comput. Simul., 185 (2021), pp. 62-76

[9]
M. Ali, M. Alquran, I. Jaradat. Int. J. Appl. Comput. Math, 7 (90) (2021)

[10]
E. Atilgan, M. Senol, A. Kurt,O. Tasbozan. China Ocean Eng., 33 (2019), pp. 477-483 DOI: 10.1007/s13344-019-0045-1

[11]
T. Kawahara. J. Phys. Soc. Jpn., 33 (1972), pp. 260-264

[12]
A. Biswas. Appl. Math. Lett., 22 (2) (2009), pp. 208-210

[13]
E. Yusufoǧlu, A. Bekir. Comput. Math. Appl., 55 (6) (2008), pp. 1113-1121

[14]
R. Khalil, M. Al-Horani, A. Yousef, M. Sababheh. J. Comput. Appl. Math., 264 (2014), pp. 65-70

[15]
M. Alquran,O. Yassin. Opt. Quantum Electron., 50 (2018), p. 354

[16]
M.M. Miah, A.R. Seadawy, H.M.S. Ali, M.A. Akbar. J. Ocean Eng. Sci., 5 (3) (2020), pp. 269-278

[17]
M.A. Shallal, K.K. Ali, K.R. Raslan, H. Rezazadeh, A. Bekir. J. Ocean Eng. Sci., 5 (3) (2020), pp. 223-229

[18]
I. Jaradat, M. Alquran. J. Ocean Eng. Sci. (2021)

[19]
I. Jaradat,M. Alquran. Mathematics, 8 (7) (2020), p. 1127 DOI: 10.3390/math8071127

[20]
I. Jaradat, M. Alquran, M. Ali, N. Al-Ali, S. Momani. Heliyon, 6 (6) (2020), p. e04057

[21]
M. Alquran. Results Phys., 28 (2021), p. 104577

[22]
M. Alquran, M. Ali,H. Jadallah. J. Ocean Eng. Sci. (2021),

[23]
M. Alquran, T.A. Sulaiman,A. Yusuf. Opt. Quantum Electron., 53 (2021), p. 227

[24]
M. Alquran, I. Jaradat, T.A. Sulaiman,A. Yusuf. Opt. Quantum Electron., 53 (2021), p. 26

[25]
E.A. Az-Zo'bi, W.A. Alzoubi, M. Senol, B.S. Masaedeh. Mod. Phys. Lett. B, 35 (15) ( 2021), p. 2150254 DOI: 10.1142/s0217984921502547

[26]
M. Alquran, A. Jarrah. J. King Saud Univ. -Sci., 31 (4) (2019), pp. 485-489

[27]
A. Ali, A.R. Seadawy. J. Ocean Eng. Sci., 6 (1) (2021), pp. 85-98

[28]
M. Mirzazadeh, L. Akinyemi, M. Senol,K. Hosseini. Optik, 241 (2021), p. 166318

[29]
R.M. El-Shiekh, M. Gaballah. J. Ocean Eng. Sci., 5 (2) (2020), pp. 180-185

[30]
W. Malfliet, W. Hereman. Phys. Scr.,54 (1996), pp. 563-568 DOI: 10.1088/0031-8949/54/6/003

[31]
A.M. Wazwaz. Commun. Nonlinear Sci. Numer. Simul., 11 (2006), pp. 148-160

[32]
M. Alquran. Appl. Math. Inf. Sci., 6 (1) (2012), pp. 85-88

[33]
M. Alquran, A. Qawasmeh. Nonlinear Stud., 20 (2) (2013), pp. 263-272

[34]
M. Alquran, K. Al-Khaled, M. Ali, O.A. Arqub. Waves Wavelets Fractals Adv. Anal., 3 ( 2017), pp. 31-39 DOI: 10.1515/wwfaa-2017-0003

Options
Outlines

/