The extended Jacobi elliptic function expansion method is used for solving fractional differential equations
in the sense of Jumarie’s modified Riemann-Liouville derivative. By means of this approach, a few fractional differential
equations are successfully solved. As a result, some new Jacobi elliptic function solutions including solitary wave solutions
and trigonometric function solutions are established. The proposed method can also be applied to other fractional differential
equations.
1. Introduction
Fractional differential equations attracted attention in physics, biology, engineering, signal processing, systems identification, control theory, finance, and fractional dynamics [1–3]. Also, they are employed in social sciences such as food supplement, climate, finance, and economics.
Finding approximate and exact solutions to fractional differential equations is an important task. Various analytical and numerical methods have been introduced to obtain solutions of fractional differential equations, such as the Adomian decomposition method [4, 5], the variational iteration method [6–8], the homotopy analysis method [9–12], the homotopy perturbation method [13–15], the Lagrange characteristic method [16], the finite difference method [17], the finite element method [18], the differential transformation method [19], the fractional subequation method [20–24], the first integral method [25], the (G′/G)-expansion method [26–29], the fractional complex transform method [30], and the modified simple equation method [31–33].
In [34], Jumarie proposed a modified Riemann-Liouville derivative. With this kind of fractional derivative and some useful formulas, we can convert fractional differential equations into integer-order differential equations by variable transformation.
In this paper, we used extended Jacobi elliptic function expansion method [35–37] to establish exact solutions for three nonlinear space-time fractional differential equations in the sense of Jumarie’s modified Riemann-Liouville derivative, namely, the space-time fractional generalized reaction duffing equation, the space-time fractional bidirectional wave equations, and the space-time fractional symmetric regularized long wave (SRLW) equation. Also, we included figures to show the properties of some Jacobi elliptic function solutions of these fractional differential equations.
2. Jumarie’s Modified Riemann-Liouville Derivative and the Extended Jacobi Elliptic Function Expansion Method
In this section, we first give the definition and some properties of the modified Riemann-Liouville derivative which are used further in this paper.
The Jumarie modified Riemann-Liouville derivative of order α is defined by the expression [34](1)Dxαfx=1Γ-α∫0xx-ξ-α-1fξ-f0dξ,α<0,1Γ1-αddx∫0xx-ξ-αfξ-f0dξ,0<α<1,fnxα-n,n≤α<n+1,n≥1,where f:R→R, x→f(x) denote a continuous (but not necessarily differentiable) function.
Some properties of the fractional modified Riemann-Liouville derivative were summarized and three useful formulas of them are [34](2)Dxαxr=Γ1+rΓ1+r-αxr-α,Dxαfxgx=gxDxαfx+fxDxαgx,Dxαfgx=fg′gxDxαgx=Dgαfgxg′xα.
Next, let us consider nonlinear partial fractional differential equation(3)Pu,Dtαu,Dxβu,Dyγu,DtαDtαu,DtαDxβu,DxβDxβu,DxβDyγu,DyγDyγu,…=0,0<α,β,γ≤1,where u is an unknown function and P is a polynomial of u. In this equation, the partial fractional derivatives involving the highest order derivatives and the nonlinear terms are included.
Li and He [38] presented a fractional complex transform to convert fractional differential equations into ordinary differential equations (ODEs), so all analytical methods devoted to the advanced calculus can be easily applied to the fractional calculus. By using the traveling wave variable(4)ux,y,t=Uξ,ξ=δxβΓ1+β+ζyγΓ1+γ+λtαΓ1+α,where δ, ζ are nonzero arbitrary constants and λ is the wave speed, we can rewrite (3) as the following nonlinear ODE:(5)QU,U′,U′′,U′′′,…=0,where the prime denotes the derivation with respect to ξ. If possible, we should integrate (5) term by term one or more times.
Our main goal is to derive exact or at least approximate solutions, if possible, for this ODE. For this purpose, using the extended Jacobi elliptic function expansion method, U(ξ) can be expressed as a finite series of Jacobi elliptic functions, snξ, that is, the ansatz:(6)ux,y,t=Uξ=∑j=0najsnjξ+∑j=1nbjsn-jξ.The parameter n is determined by balancing the linear term(s) of highest order with the nonlinear one(s). And(7)cn2ξ=1-sn2ξ,dn2ξ=1-m2sn2ξ,ddξsnξ=cnξdnξ,ddξcnξ=-snξdnξ,ddξdnξ=-m2snξcnξ,where cnξ and dnξ are the Jacobi elliptic cosine function and the Jacobi elliptic function of the third kind, respectively, with the modulus m(0<m<1). Therefore, the highest degree of dpU/dξp is taken as(8)OdpUdξp=n+p,p=1,2,3,…,OUqdpUdξp=q+1n+p,q=0,1,2,…,p=1,2,3,….Substituting (6)–(8) into (5) and comparing the coefficients of each power of snξ in both sides, we get an overdetermined system of nonlinear algebraic equations with respect to λ, aj(j=0,1,…,n), and bj(j=1,2,…,n). Solving this system, with the aid of Mathematica, then λ, aj(j=0,1,…,n), and bj(j=1,2,…,n) can be determined. Substituting these results into (6), then some new Jacobi elliptic function solutions of (3) can be obtained. We can get other kinds of Jacobi doubly periodic wave solutions.
Since(9)limm→1snξ=tanhξ,limm→1cnξ=sechξ,limm→1dnξ=sechξ,limm→0snξ=sinξ,limm→0cnξ=cosξ,limm→0dnξ=1,u degenerates, respectively, as the following form.
Triangular function formal solution: (11)ux,y,t=∑j=0najsinjξ+∑j=1nbjcscjξ.
3. Applications of the Method
In this section, we present three examples to demonstrate the effectiveness of our approach to solve nonlinear fractional partial differential equations.
We have applied the extended Jacobi elliptic function expansion method to construct the exact solutions of space-time fractional generalized reaction duffing equation [39, 40] in the form(12)∂2αu∂t2α+p∂2αu∂x2α+qu+ru2+su3=0,0<α<1,where p, q, r, and s are all constants. Equation (12) reduces many well-known nonlinear fractional wave equations such as the following.
For our purpose, we introduce the following transformations:(18)ux,t=Uξ,ξ=lxαΓ1+α-λtαΓ1+α,where ξ is a wave variable and l and λ are constants; all of them are to be determined. Substituting (18) into (12), (12) is reduced into an ODE:(19)U′′ξ+qλ2+pl2Uξ+rλ2+pl2U2ξ+sλ2+pl2U3ξ=0,where U′=dU/dξ. Suppose that the solution of (19) can be expressed by(20)Uξ=∑j=0najsnjξ+∑j=1nbjsn-jξ.Considering the homogeneous balance between the highest order derivative U′′ and the highest order nonlinear term U3 in (19), we obtain n=1. So(21)Uξ=a0+a1snξ+b1sn-1ξ.Substituting (21) into (19) and comparing the coefficients of each power of snξ in both sides, we get an overdetermined system of nonlinear algebraic equations with respect to λ, a0, a1, and b1. Solving this system with Mathematica, we get the following results.
When the modulus m approaches to 1 in (25), (26), and (27), we can obtain solitary wave solutions of space-time fractional generalized reaction duffing equation, respectively:(28)u4=-r3s±r3stanhlxαΓ1+α∓ir2+18l2ps32sΓ1+αtα,u5=-r3s±r3scothlxαΓ1+α∓ir2+18l2ps32sΓ1+αtα,u6=-r3s±r3scoth2lxαΓ1+α+ir2+72l2ps38sΓ1+αtα.
3.1.2. Triangular Periodic Solutions
When the modulus m approaches to zero in (26), (27), we can obtain trigonometric function solutions of space-time fractional generalized reaction duffing equation, respectively: (29)u7=-r3s±2r3scsclxαΓ1+α∓ir2+9l2ps3sΓ1+αtα,u8=-r3s±2r3s+csclxαΓ1+α+ir2+9l2ps3sΓ1+αtα.
Let us apply our method to the space-time fractional bidirectional wave equations in the form [41, 42](30)Dtαv+Dxαu+uDxαv+vDxαu+aDxαDxαDxαu-bDxαDxαDtαv=0,Dtαu+Dxαv+uDxαu+cDxαDxαDxαv-dDxαDxαDtαu=0,0<α≤1,where x represents the distance along the channel, t is the elapsed time, the variable u(x,t) is the dimensionless horizontal velocity, v(x,t) is the dimensionless deviation of the water surface from its undisturbed position, and a, b, c, and d are real constants. When α=1, (30) is the generalization of bidirectional wave equations, which can be used as a model equation for the propagation of long waves on the surface of water with a small amplitude by Bona and Chen [43].
For our purpose, we use the following transformation:(31)ux,t=Uξ,vx,t=Vξ,ξ=RxαΓ1+α+StαΓ1+α,where R and S are nonzero constants. Substituting (31) into (30), we obtain(32)SV′+RU′+RUV′+RVU′+aR3U′′′-bR2SV′′′=0,SU′+RV′+RUU′+cR3V′′′-dR2SU′′′=0,where U′=dU/dξ. Suppose that the solutions of (32) can be expressed by(33)Uξ=∑j=0n1ajsnjξ+∑j=1n1bjsn-jξ,Vξ=∑j=0n2cjsnjξ+∑j=1n2djsn-jξ.Balancing the highest order derivative terms and nonlinear terms in (32), we can obtain n1=n2=2. So we have(34)Uξ=a0+a1snξ+a2sn2ξ+b1sn-1ξ+b2sn-2ξ,Vξ=c0+c1snξ+c2sn2ξ+d1sn-1ξ+d2sn-2ξ.Proceeding as in the previous case, we get the following results.
When the modulus m approaches to 1 in (37), (38), we can obtain solitary wave solutions of the space-time fractional bidirectional wave equations, respectively:(40)u3=i2c+dc2+6aR2-bc21-8cR22R3bcbc2-2dc2+6aR2-i·23cbc2Rbbc2-2dc2+6aR2·tanh2RxαΓ1+α+SΓ1+αtα,v3=-bc2+4cR23+2c2-2dc2+6aR212bcR2+c2·tanh2RxαΓ1+α+SΓ1+αtα,u4=i2c+dc2+6aR2-bc21-8cR22R3bcbc2-2dc2+6aR2-i·43cbc2Rbbc2-2dc2+6aR21+2csch22RxαΓ1+α+SΓ1+αtα,v4=-bc2+4cR23+2c2-2dc2+6aR212bcm2R2+2c21+2csch22RxαΓ1+α+SΓ1+αtα,where(41)S=-icc2+6aR23bbc2-2dc2+6aR2.
3.3. The Space-Time Nonlinear Fractional SRLW Equation
We consider the space-time nonlinear fractional SRLW equation [44, 45](42)Dt2αu+Dx2αu+uDtαDxαu+DtαuDxαu+Dt2αDx2αu=0,0<α≤1,which arises in several physical applications including ion sound waves in plasma. For our purpose, we use the following transformation:(43)ux,t=Uξ,ξ=kxαΓ1+α+ctαΓ1+α+ξ0,where k, c, and ξ0 are constants with k, c≠0. Substituting (43) into (42), we obtain(44)2k2c2U′′+2k2+c2U+kcU2=0.Suppose that the solutions of (44) can be expressed by(45)Uξ=∑j=0najsnjξ+∑j=1nbjsn-jξ.Considering the homogeneous balance between the highest order derivative U′′ and the highest order nonlinear term U2 in (44), we obtain n=2. So we have(46)Uξ=a0+a1snξ+a2sn2ξ+b1sn-1ξ+b2sn-2ξ.Proceeding as in the previous cases, we get the following results.
When the modulus m approaches to 1 in (53), (54), and (55), we can obtain solitary wave solutions of the space-time nonlinear fractional SRLW equation, respectively:(56)u4=12cksech2kxαΓ1+α+iΓ1+αk21+4k2tα+ξ0,u5=-8ck-12ckcsch2kxαΓ1+α+iΓ1+αk21-4k2tα+ξ0,u6=-48ckcsch22kxαΓ1+α+iΓ1+αk21+16k2tα+ξ0.
3.3.2. Triangular Periodic Solutions
We can obtain trigonometric function solutions of the space-time nonlinear fractional SRLW equation, when the modulus m approaches to zero; for example, (54), (55) give the same solution:(57)u7=-12ck·csc2kxαΓ1+α+iΓ1+αk21-4k2tα+ξ0,u8=8ck-12ck·csc2kxαΓ1+α+iΓ1+αk21+4k2tα+ξ0.
4. Conclusion
In this paper, we used the extended Jacobi elliptic function expansion method for solving fractional differential equations and applied it to find exact solutions of the space-time fractional generalized reaction duffing equation, the space-time fractional bidirectional wave equations, and the space-time fractional symmetric regularized long wave (SRLW) equation. With the aid of Mathematica, we successfully obtained some new Jacobi elliptic function solutions including solitary wave solutions and trigonometric function solutions for these equations. This method is effective and can also be applied to other fractional differential equations.
Conflict of Interests
The author declares that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The author wishes to thank the referees for their valuable suggestions. He thanks Mr. Mehmet Ekici from the Department of Mathematics, Bozok University, Yozgat, Turkey. This paper is supported by the Scientific and Technological Research Council of Turkey (TUBITAK).
MillerK. S.RossB.1993New York, NY, USAJohn Wiley & SonsMR1219954KilbasA. A.SrivastavaH. M.TrujilloJ. J.2006San Diego, Calif, USAElsevierPodlubnyI.1999198San Diego, Calif, USAAcademic PressMathematics in Science and EngineeringMR1658022El-SayedA. M. A.RidaS. Z.ArafaA. A. M.Exact solutions of fractional-order biological population model200952699299610.1088/0253-6102/52/6/04MR26830182-s2.0-73849139817SafariM.GanjiD. D.MoslemiM.Application of He's variational iteration method and Adomian's decomposition method to the fractional KdV-Burgers-KURamoto equation20095811-122091209710.1016/j.camwa.2009.03.043MR25573322-s2.0-70350566604IncM.The approximate and exact solutions of the space- and time-fractional Burgers equations with initial conditions by variational iteration method2008345147648410.1016/j.jmaa.2008.04.007MR2422665ZBL1146.353042-s2.0-43949121726WuG.-C.LeeE. W.Fractional variational iteration method and its application2010374252506250910.1016/j.physleta.2010.04.034MR26400232-s2.0-77953478991FouladiF.HosseinzadehE.BarariA.DomairryG.Highly nonlinear temperature-dependent fin analysis by variational iteration method201041215516510.1615/heattransres.v41.i2.402-s2.0-77951790468SongL.ZhangH.Solving the fractional BBM-Burgers equation using the homotopy analysis method20094041616162210.1016/j.chaos.2007.09.042MR25312522-s2.0-65449123574AbbasbandyS.ShirzadiA.Homotopy analysis method for multiple solutions of the fractional Sturm-Liouville problems201054452153210.1007/s11075-009-9351-7MR2670667ZBL1197.650962-s2.0-77954459409BararniaH.DomairryG.GorjiM.RezaniaA.An approximation of the analytic solution of some nonlinear heat transfer in fin and 3D diffusion equations using HAM201026111310.1002/num.20404MR25889022-s2.0-73449135842RashidiM. M.DomairryG.DoostHosseiniA.DinarvandS.Explicit approximate solution of the coupled KdV equations by using the homotopy analysis method200829–12581589MR2482743GanjiZ. Z.GanjiD. D.GanjiA. D.RostamianM.Analytical solution of time-fractional Navier-Stokes equation in polar coordinate by homotopy perturbation method201026111712410.1002/num.20420MR25889102-s2.0-73449083810GepreelK. A.The homotopy perturbation method applied to the nonlinear fractional Kolmogorov-Petrovskii-Piskunov equations20112481428143410.1016/j.aml.2011.03.025MR27936462-s2.0-79955470464GuptaP. K.SinghM.Homotopy perturbation method for fractional Fornberg-Whitham equation201161225025410.1016/j.camwa.2010.10.045MR27541332-s2.0-78651257332GuyJ.Lagrange characteristic method for solving a class of nonlinear partial differential equations of fractional order200619987388010.1016/j.aml.2005.10.016MR22404772-s2.0-33745108563CuiM.Compact finite difference method for the fractional diffusion equation2009228207792780410.1016/j.jcp.2009.07.021MR2561843ZBL1179.651072-s2.0-69049086472HuangQ.HuangG.ZhanH.A finite element solution for the fractional advection-dispersion equation200831121578158910.1016/j.advwatres.2008.07.0022-s2.0-55649099424OdibatZ.MomaniS.Fractional Green function for linear time-fractional equations of fractional order200821194199ZhangS.ZhangH.-Q.Fractional sub-equation method and its applications to nonlinear fractional PDEs201137571069107310.1016/j.physleta.2011.01.029MR27650132-s2.0-79251635229GuoS.MeiL.LiY.SunY.The improved fractional sub-equation method and its applications to the space-time fractional differential equations in fluid mechanics2012376440741110.1016/j.physleta.2011.10.056MR28777502-s2.0-84855193408LuB.Bäcklund transformation of fractional Riccati equation and its applications to nonlinear fractional partial differential equations201237628-292045204810.1016/j.physleta.2012.05.013MR29293922-s2.0-84861576550TangB.HeY.WeiL.ZhangX.A generalized fractional sub-equation method for fractional differential equations with variable coefficients201237638-392588259010.1016/j.physleta.2012.07.018MR29611212-s2.0-84864671185ZhengB.WenC.Exact solutions for fractional partial differential equations by a new fractional sub-equation method20132013, article 19910.1186/1687-1847-2013-199MR30841942-s2.0-84889055780LuB.The first integral method for some time fractional differential equations2012395268469310.1016/j.jmaa.2012.05.066MR29482592-s2.0-84864025292ZhengB.(G′/G)-expansion method for solving fractional partial differential equations in the theory of mathematical physics201258562363010.1088/0253-6102/58/5/02MR30890912-s2.0-84870265361GepreelK. A.OmranS.Exact solutions for nonlinear partial fractional differential equations2012211111020410.1088/1674-1056/21/11/1102042-s2.0-84869126980BekirA.GünerÖ.Exact solutions of nonlinear fractional differential equations by (G′/G)-expansion method2013221111020210.1088/1674-1056/22/11/1102022-s2.0-84888632803ZhengB.Exact solutions for some fractional partial differential equations by the (G′/G)
method201320131382636910.1155/2013/826369MR3123716SuW.-H.YangX.-J.JafariH.BaleanuD.Fractional complex transform method for wave equations on cantor sets within local fractional differential operator20132013, article 9710.1186/1687-1847-2013-972-s2.0-84879315907JawadA. J. M.PetkovićM. D.BiswasA.Modified simple equation method for nonlinear evolution equations2010217286987710.1016/j.amc.2010.06.030MR26786032-s2.0-77955664465ZayedE. M. E.A note on the modified simple equation method applied to Sharma-Tasso-Olver equation201121873962396410.1016/j.amc.2011.09.025MR28514952-s2.0-80054983838ZayedE. M. E.IbrahimS. A. H.Exact solutions of nonlinear evolution equations in mathematical physics using the modified simple equation method20122964060201JumarieG.Modified Riemann-Liouville derivative and fractional Taylor series of nondifferentiable functions further results2006519-101367137610.1016/j.camwa.2006.02.001MR22376342-s2.0-33745742268ZhangH.Extended Jacobi elliptic function expansion method and its applications200712562763510.1016/j.cnsns.2005.08.003MR22875042-s2.0-33845601704AbdouM. A.ElhanbalyA.Construction of periodic and solitary wave solutions by the extended Jacobi elliptic function expansion method20071271229124110.1016/j.cnsns.2006.01.013MR23215492-s2.0-34047154416BhrawyA. H.AbdelkawyM. A.BiswasA.Cnoidal and snoidal wave solutions to coupled nonlinear wave equations by the extended Jacobi's elliptic function method201318491592510.1016/j.cnsns.2012.08.034MR2996603ZBL1261.350442-s2.0-84868193363LiZ.-B.HeJ.-H.Fractional complex transform for fractional differential equations2010155970973MR2777702EslamiM.Fathi VajargahB.MirzazadehM.BiswasA.Application of first integral method to fractional partial differential equations201488217718410.1007/s12648-013-0401-62-s2.0-84898009585ArnousA. H.MirzazadehM.Bäcklund transformation of fractional Riccati equation and its applications to the space-time FDEs201410.1002/mma.3371LuB.Bäcklund transformation of fractional Riccati equation and infinite sequence solutions of nonlinear fractional PDEs20142014657205210.1155/2014/572052MR31666302-s2.0-84893159743LiW.YangH.HeB.Exact solutions of the space-time fractional bidirectional wave equations using the (G′/G)-expansion method20142014915370610.1155/2014/153706MR3228124BonaJ. L.ChenM.A Boussinesq system for two-way propagation of nonlinear dispersive waves19981161-219122410.1016/s0167-2789(97)00249-2MR1621920XuF.Application of Exp-function method to symmetric regularized long wave (SRLW) equation2008372325225710.1016/j.physleta.2007.07.035MR23800842-s2.0-37249080286AlzaidyJ. F.The fractional sub-equation method and exact analytical solutions for some nonlinear fractional PDEs201311141910.12691/ajma-1-1-3