SciELO - Scientific Electronic Library Online

 
vol.65 número5Asymptotically equivalent functions and ultrafilters applied to noncommutative quantum cosmologyNew reflections on higher dimensional linearized gravity índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Revista

Articulo

Indicadores

Links relacionados

  • No hay artículos similaresSimilares en SciELO

Compartir


Revista mexicana de física

versión impresa ISSN 0035-001X

Rev. mex. fis. vol.65 no.5 México sep./oct. 2019  Epub 23-Abr-2020

https://doi.org/10.31349/revmexfis.65.529 

Research

in Gravitation, Mathematical Physics and Field Theory

Symmetry properties and exact solutions of the time fractional Kolmogorov-Petrovskii-Piskunov equation

M.S. Hashemi a  

M. Inc b  

M. Bayram c  

a Department of Mathematics, Basic Science Faculty, University of Bonab, Bonab 55517, Iran.

b Department of Mathematics, Science Faculty, Firat University, 23119 Elazig, Turkey.

c Department of Computer Engineering, Istanbul Gelisim University, Istanbul, Turkey.


Abstract

In this paper, the time fractional Kolmogorov-Petrovskii-Piskunov (TFKPP) equation is analyzed by means of Lie symmetry approach. The TFKPP is reduced to ordinary differential equation of fractional order via the attained point symmetries. Moreover, the simplest equation method is used in construct the exact solutions of underlying equation with recently introduced conformable fractional derivative.

Keywords: Time fractional Kolmogorov-Petrovskii-Piskunov equation; Lie symmetry analysis; Erdélyi-Kober fractional derivative; Riemann-Liouville derivative; conformable fractional derivative; simplest equation method

PACS: 02.30.Uu; 04.20.Jb; 05.40.Fb; 05.60.-k

1. Introduction

The fractional calculus (FC) began to wind up exceptionally famous in a few parts of science and engineering. Numerous important event, that is, acoustics, anomalous diffusion, chemistry, control processing, electro-magnetics, and visco-elasticity have been expressed by FC. It is known that a systematic method for extracting the analytical solution of both ordinary differential equations (ODEs) and partial differential equations (PDEs) was first proposed by the Norwegian mathematician Sophus Lie in the early 19th century. The fundamental overview of this strategy is the estimation of variable changes that can leave differential condition unchanged. Therefore, a vital role in the field of FC is to attain the Lie symmetries and the solutions of the equations with the FC derivatives. There have been some properties of the fractional sense that could not be found in classical sense, owing to this we feel motivated to establish the symmetries of TFKPP equation. This equation has the following generalized form [1-4]

tαu=uxx+λu+μu2+γu3,λ+μ+γ=0,ψ2=μ2-4λγ0, (1)

where tαu:=Dtαu stands for Riemann-Liouville of order α, expressed as [5]

Dtαu(x,t)=1Γn-αntn0tt-ξn-α-1ux,ξdξ,                                                                             n-1<α<nnutn,  α=nN (2)

The TFKPP Eq. (1), has a large application and includes as particular cases the time fractional Fitzhugh-Nagumo equation (λ= -c, μ = c + 1, γ = -1; 0 < c < 1), which is used in population genetics, the time fractional Newell-Whitehead equation (λ = 1, μ = 0, γ = -1). Recently, the homotopy perturbation method and homotopy analysis method have utilized to consider the TFKPP equation by Gepreel [1] and Hariharan [2], respectly with λ = μ = 0 and γ = -2.

In FC, there are large amount of differential derivatives were defined e.g. [6-9]. In the calculus, the chain rule is a useful and an applicable. It is also hold for conformable fractional derivatives.

As far as we know, every proposed fractional derivative has some disadvantages. Therefore, Khalil et al., [9], proposed a new definitions:

Definition 1.1. Surmise that 𝑓: [a, b] × (0, ∞) ℝ, then the conformable fractional derivative of f is given by

tTαfx,t=limϵ0fx,t+ϵt1-α-fx,tϵ,α0,1, (3)

for all t > 0.

Theorem 1.1 [9] Suppose that a, b ∈ ℝ and α ∈ (0, 1], then

i tTα(au+bv)=atTα(u)+btTα(v),iitTαtλ=λtλ-α,λRiiitTα(uv)=utTα(v)+vtTα(u)vitTαuv=utTαv-vtTαuv2vtTαut=t1-αu't, uC1.

More than that, the chain rule is valid for conformable fractional derivatives, shown by Abdeljawad [10].

Theorem 1.2. Surmise that f: (0, ∞) is a real differentiable, α-differentiable function. Assume that g is a function defined in the range of f and also differentiable; then, one has the following rule:

tTαfogt=t1-αg'tf'gt. (4)

There are many investigation about conformable fractional derivatives [11-14] and also some physical interpretations of this newly introduced fractional derivative are described in [15].

The organization of the manuscript is given below: In Sec. 2, we provide some preliminaries. Section 3, is devoted to the description of Lie symmetry analysis of TFKPP Eq. (1). General similarity forms and symmetry reductions are established. In Sec. 4, exact solutions to the TFKPP equation with conformable fractional derivative are investigated. Finally, the last section is devoted to conclusions.

2. Lie symmetry analysis of fractional partial differential equations

Here, some description for solving fractional partial differential equations (FPDEs) via Lie symmetry analysis will be provided. Surmise that FPDE having as in [16-26]

tαu=Fx,t,u,ux,uxx,0<α<1. (5)

If (5) is invariant under a one parameter Lie group of point transformations

t¯=t¯x,t,u;ϵ, x¯=x¯x,t,u;ϵ, u¯=u¯x,t,u;ϵ, (6)

the vector field of an evolution type of equation is as follows:

V=ξtx,t,ut+ξxx,t,ux+ϕx,t,uu,

where the coefficients ξ t , ξ x and φ of the vector field are to be determined. When V satisfy the Lie symmetry condition, the vector field (7) generates a symmetry of (5),

prα,2VΔ|Δ=0=0,Δ=tαu-F

Thus the extension operator take the form

pr(α,2)V=V+ϕα0tαu+ϕxux+ϕxxuxx,

where

ϕx=Dxϕ-uxDxξx-utDxξt, ϕxx=Dxϕx-uxtDxξt-uxxDxξx, ϕα0=Dtαϕ+ξxDtαux-Dtαξxux+DtαDtξtu-Dtα+1ξtu+ξtDtα+1u.

The condition of invariance

ξtx,t,u|t=0=0,

is inevitable for the (6), due to the (2).

The α th extended infinitesimal is presented as:

ϕα0=Dtαϕ+ξxDtαux-Dtαξxux+DtαDtξtu-Dtα+1ξtu+ξtDtα+1u, (8)

where Dtα exhibits the total fractional derivative operator. The fractional generalized Leibnitz rule is expressed as

Dtαutvt=n=0αnDtα-nutDtnvt, α>0, (9)

here

αn=-1n-1αΓn-αΓ1-αΓn+1.

Therefore using (9) one can represent (8) as

ϕα0=Dtαϕ-αDtξtαutα-n=1αnDtnξxDtα-nux-n=1αn+1Dtn+1ξtDtα-nu. (10)

Using chain rule

dmfgtdtm=k=0m r=0kkr×1k!-gtrdmdtmgtk-rdkfgdgk,

and setting 𝑓(t) = 1, one can get

Dtαϕ=αϕtα+ϕuαutα-uαϕutα+n=1αnnϕutnDtα-nu+ϑ,

where

ϑ=n=2 m=2n k=2m r=0k-1αnnmkr1k!×tn-αΓ(n+1-α)-urmtmuk-rn-m+kϕtn-muk.

Therefore

ϕα0=αϕtα+(ϕu-αDt(ξt))αutα-uαϕutα+ϑ+n=1αnαϕutα-αn+1Dtn+1ξtDtα-nu-n=1αnDtnξxDtα-nux.

3. Symmetry representation of TFKPP equation

In view of the Lie theory, we have:

ϕα0=ϕxx+λϕ+2μϕu+3γϕu2. (11)

Substituting (10) into (11), the determining equations for Eq. (1) is attained, consequently, we have

ξt=4tc3, ξx=c1+2αxc3,ϕ=c2u+3α-2uc3+Cx,t,

where c 1 ; c 2 and c 3 are constants and C (x, t) is a solution of Eq. (1). Therefore, the algebra g of Eq. (1) can be written as

V1=x, V2=uu,V3=4tt+2αxx+(3α-2)uu,V4=Cx,tu.

For V 3, one can write

dt4t=dx2αx=du3α-2u,

and this give

ζ=xt-α2, u(x,t)=t3α-24F(ζ). (12)

Theorem 3.1. The transformation (12) reduces (1) to the following:

P2α-α4+12,αFζ=F+λF+μF2+γF3, (13)

with the Erdélyi-Kober (EK) fractional differential operatorPβτ,αdefined by

Pβτ,αF:=j=0n-1τ+j-1βζddζKβτ+α,n-αFζ, n=α+1, αα,    α

where

Kβτ,αFζ:=1Γα1u-1α-1u-τ+αFζu1βdu,F(ζ), α=0,

is the EK fractional integral operator.

Proof: Let n −1 < α < n, n = 1,2,3,.... By means of Reimann-Liouville, one reaches

αutα=ntn1Γn-α×0tt-sn-α-1s3α-24Fxs-α2ds. (14)

Letting ρ = t/s, one can get ds = -(t/ρ 2), therefore (14) can be written as

αutα=ntntn-α4-12K2α3α+24,n-αFζ.

Taking into account the relation (ζ = xt -α/2 ), we can obtain

ttϕζ=tζtdϕζdζ=-α2ζdϕζdζ.

Therefore one can get

ntntn-α4-12K2α3α+24,n-αFζ=n-1tn-1ttn-α4-12K2α3α+24,n-αFζ=n-1tn-1tn-α4-32n-α4-12-α2ζddζ×K2α3α+24,n-αFζ==t-α4-12j=0n-1-α4+12+j-α2ζddζ×K2α3α+24,n-αFζ=t-α4-12P2α-α4+12,αFζ.

This completes the proof.

Also, for the symmetry of V 1 + V 2 + V 3, one can write

dt4t=dx2αx+1=du3α-1u,

which yields

ζ=2αx+12αt-α2, ux,t=t3α-14Fζ. (15)

Theorem 3.2. The transformation (15) reduces (1) to the following nonlinear ordinary differential equation of fractional order:

P2α3-α4,αFζ=F+λF+μF2+γF3. (16)

Proof: Similar to the proof of previous theorem.

4. Exact Solutions of TFKPP equation

Symmetry analysis of differential equations gives many information about geometric properties of various differential equations. For example, it is possible to extract vector fields, infinitesimals, conservation laws and reductions of differential equations. Reduction procedure of differential equations allows us to reduce dimension of these equations by one less. In two dimensional partial differential equations (PDEs), reduction procedure gives an ordinary differential equation (ODE). So, solving this ODE concludes exact solution of original PDE. However, in FPDEs with Riemann-Liouville fractional derivatives we get ODEs with the EK derivatives which there is not a systematic method to find their exact solution. Therefore, after reduction of TFKPP equation with the Riemann-Liouville fractional derivative we obtain Eqs. (13) and (16) which it is not possible to find analytical solutions. However, we can obtain exact solution of Eq. (1) with tαu:= tTα(u). In this section, we investigate the exact solutions of TFKPP equation with conformable fractional derivative.

4.1. Simplest equation method and its applications to time fractional differential equations

This approach was proposed in [27,28]. The steps for the approach is stated as follows:

Let the TFDE is given by

Pu, tTαu,ux,uxx,...=0. 0<α1, (17)

Then the modified version of simplest equation method procedure have the following steps:

Step 1: We utilize the following

ux,t=Θξ, ξ=Ax-νtαα, (18)

where A and v are nonzero constants to be determined later.

Consequently we attain with parameters A and v the following

PΘ,-AνΘ',AΘ',A2Θ,...=0. (19)

Step 2: Suppose that Eq. (19) possesses

Θξ=i=0Naizξi, (20)

where a i , i = 0, 1, …, N, are constants to be determined later. The positive value of N in (20), which the pole order for the general solution of Eq. (19), can be determined by substituting Θ(ξ) = ξ -m , (m > 0).

In the present paper, we use the Bernoulli and Riccati equations which their solutions can be expressed by elementary functions. For the Bernoulli equation:i

dzdξ=azξ+bzξk, kN1,

we use the solutions

zξ=a expak-1ξ+ξ01-b expak-1ξ+ξ0k-1,

for the case a < 0; b > 0 and ξ 0 is a constant of integration. For the Riccati equation

dzdξ=a+bzξ2,

which admits the following exact solutions:

zξ=--abbtanh-abξ-ϵlnξ02, ξ0>0, ϵ=±1,

when ab < 0 and

zξ=abbtanabξ+ξ0, ξ0=Const.,

when ab > 0.

Step 3: Plugging (20) into (19) and equating the coefficients of z i to zero, one can obtain an algebraic system in A; v and a i , i = 0; …, N.

4.2. Application to the TFKPP equation

The transformation

ux,t=Θξ, ξ=Ax-νtαα, (21)

changes Eq. (1) with tαu=tTα(u) to:

A2Θ+νAΘ'+λΘ+μΘ2+γΘ3=0. (22)

We suppose that Eq. (22) has solution of the form (20). Balancing the highest order derivative terms with nonlinear terms in Eq. (22), we get N = 1, and hence

Θξ=a0+a1zξ, a10. (23)

Substituting (23) along with (21) into Eq. (22) and then vanishing the coefficients of z i , one can get some algebraic equations about a 0 , a 1 , A and v, which solving them by Maple, concludes:

• Case 1:

a0=0, a1=bμ2|μ|ψ2aμγ,A=-μ2+2λγ±|μ|ψ2aγ,ν=μ2-6λγ|μ|ψ2γ-μ2+2λγ±|μ|ψ,

where ψ=μ2-4λγ. In this case, the exact solutions of Eq. (22) are

Θ(ξ)=b-μ2±μψ expaξ+ξ02γμ(-1+b expaξ+ξ0), a>0, b<0,Θ(ξ)=b-μ2±μψ expaξ+ξ02γμ1+b expaξ+ξ0, a<0, b>0,

and using the substitution in (18) we get the final solutions:

u(x,t)=b-μ2±μψ expηx,t2γμ-1+b expηx,t, a>0, b<0,u(x,t)=b-μ2±μψ expηx,t2γμ1+b expηx,t, a<0, b>0,

where

ηx,t=2xγ-μ2+2γλ±|μ|ψ4γ+-μ2+6γλ±|μ|ψtα+4ξ0γaα4γα.

• Case 2:

a0=ψ-μ2γ, a1=b(3μ2-12λγ-μψ)γa(3ψ-μ),A=2ψ2a-γ, ν=2μ2-γ.

In this case, the exact solutions of Eq. (22) are:

Θ(ξ)=-2μ2+2μψ+6λγ+b6λγ-μ2-μψ expaξ+ξ0γ3ψ-μ(-1+b expaξ+ξ0),a>0,b<0,Θ(ξ)=2μ2+2μψ-6λγ+b6λγ-μ2-μψexpaξ+ξ0γ3ψ-μ1+b expaξ+ξ0,a<0,b>0,

or equivalently

u(x,t)=-2μ2+2μψ+6λγ+b6λγ-μ2-μψ expηx,tγ3ψ-μ-1+b expη(x,t), a>0, b<0,u(x,t)=2μ2+2μψ-6λγ+b6λγ-μ2-μψ expηx,tγ3ψ-μ1+b expηx,t,a<0, b>0,

where

ηx,t=ξ0a+xψ-2γ2γ-μψtα2γα.

• Case 3:

a0=ψ-μ2γ,a1=2bμ2-3λγ-μψγa3ψ-μ,A=2λγ-μ2+μψ2aγ,ν=6λγ-μ2+μψ2γ2λγ-μ2+μψ.

In this case, we can obtain

Θ(ξ)=2-μ2+3λγ+μψγ3ψ-μ-1+b expaξ+ξ0,  a>0, b<0,Θ(ξ)=2μ2-3λγ-μψγ3ψ-μ1+b expaξ+ξ0, a<0, b>0,

and using the substitution in (18) we have

u(x,t)=2-μ2+3λγ+μψγ3ψ-μ-1+b expηx,t,  a>0, b<0,u(x,t)=2μ2-3λγ-μψγ(3ψ-μ)1+b expηx,t,a<0, b>0, (24)

where

ηx,t=ξ0a+x-μ2+μψ+2γλ2γ--μ2+μψ+6γλ4γαtα.

Also, in the use of Riccati equation, substituting (23) along with (21) into Eq. (22) and then vanishing the coefficients of z i , we can obtain some algebraic equations about a 0, a 1, A and ν, that solving them by Computer algebra technique , concludes:

• Case 1:

a0=-μ2γ,a1=±ibψ2γab,A=±2ψ4γab,ν=i2μ2γ.

In this case, the exact solutions of Eq. (22) are:

Θξ=-μ±iψ tanabξ+ξ02γ, ab>0,Θξ=-μ±ψ tanh-abξ-ϵ lnξ022γ,ab<0,

and using the substitution in (18) we get the following final solutions:

ux,t=-12γμ+iψ tan-2γψxα4γξ0α+iμψtα4γα,

when ab > 0 and

ux,t=-12γμ+ψ tanh-i2γψxα±2γϵ lnξ0α-μψtα4γα,

when ab < 0.

• Case 2:

a0=-μ2γ, a1=ibψ2γab,A=±2ψ4γab, ν=-i2μ2γ,

Exact solutions of Eq. (22) extracted from this case are:

Θξ=-μiψ tanabξ+ξ02γ, ab>0,Θξ=-μψ tanh-abξ-ϵ lnξ022γ, ab<0,

or equivalently

ux,t=-12γμ+iψ tan2γψxα±4γξ0α+iμψtα4γα,

when ab > 0 and

ux,t=-12γμ+ψ tanhi2γψxα2γϵ lnξ0α-μψtα4γα,

when ab < 0.

5. Conclusion

In this study, the Lie group analysis method was successfully applied to investigate the reduction and symmetry properties of the TFKPP equation. Moreover, we have arrived to some exact solutions of the conformable TFKPP equation, thanks to the application of simplest equation method. The results of this study undoubtedly offer helpful information about the TFKPP equation.

Acknowledgments

The authors would like to thank the referees for the helpful suggestions.

References

1. K. A. Gepreel, Applied Mathematics Letters, 24 (2011) 1428-1434. [ Links ]

2. G. Hariharan, Journal of Mathematical Chemistry 51 (2013) 992-1000. [ Links ]

3. S. Li-Na, and W. Wei-Guo, Communications in Theoretical Physics, 58 (2012) 182. [ Links ]

4. A. G. Nikitin, and T. A. Barannyk, Central European Journal of Mathematics, 2 (2004) 840-858. [ Links ]

5. D. Baleanu, K. Diethelm, E. Scalas, and J. Trujillo, Fractional Calculus Models and Numerical Methods (Series on Complexity, Nonlinearity and Chaos). World Scientiffic (2012). [ Links ]

6. I. Podlubny, Fractional differential equations: an introduction to fractional derivatives, fractional differential equations, to methods of their solution and some of their applications, Vol. 198. Academic press (1998). [ Links ]

7. M. Caputo, and M. Fabrizio, Progr. Fract. Differ. Appl, 1 (2015) 1-13. [ Links ]

8. A. Atangana, and I. Koca, Chaos, Solitons & Fractals 89 (2016) 447-454. [ Links ]

9. R. Khalil, M. A. Horani, A. Yousef, and M. Sababheh, Journal of Computational and Applied Mathematics, 264 (2014) 65-70. [ Links ]

10. T. Abdeljawad, Journal of Computational and Applied Mathematics, 279 (2015) 57-66. [ Links ]

11. M. Eslami, and H. Rezazadeh, Calcolo, 53 (2016) 475-485. [ Links ]

12. K. Hosseini, P. Mayeli, and R. Ansari, Optik - International Journal for Light and Electron Optics, 130 (2017) 737-742. [ Links ]

13. Y. Çenesiz, D. Baleanu, A. Kurt, and O. Tasbozan, Waves in Random and Complex Media, 27 (2017) 103-116. [ Links ]

14. O. Tasbozan, Y. Çenesiz, andA. Kurt , The European Physical Journal Plus, 131 (2016) 244. [ Links ]

15. D. Zhao, and M. Luo, Calcolo (2017) pp. 1-15. [ Links ]

16. Q. Huang, and R. Zhdanov, Physica A: Statistical Mechanics and its Applications, 409 (2014) 110-118. [ Links ]

17. J. Hu, Y. Ye, S. Shen, and J. Zhang, Appl Math Comput, 233 (2014) 439-444. [ Links ]

18. R. Sahadevan, and T. Bakkyaraj, J Math Anal Appl, 393 (2012) 341-347. [ Links ]

19. M. S. Hashemi, and D. Baleanu, Communications in Theoretical Physics 65 (2016) 11. [ Links ]

20. S. Pashayi, M. S. Hashemi, and S. Shahmorad, Communications in Nonlinear Science and Numerical Simulation, 51 (2017) 66-77. [ Links ]

21. F. Bahrami, R. Najafi, and M. S. Hashemi, Indian Journal of Physics, 91 (2017) 1571-1579. [ Links ]

22. G. W. Wang, and M. S. Hashemi, Pramana, 88 (2017) 7. [ Links ]

23. M. S. Hashemi, F. Bahrami, and R. Najafi, Optik-International Journal for Light and Electron Optics, 138 (2017) 240-249. [ Links ]

24. R. Najafi , F. Bahrami , and M. S Hashemi, Nonlinear Dynamics 87 (2017) 1785-1796. [ Links ]

25. E. Yaşar, Y., Yıdırım, and C. M. Khalique, Results in Physics, 6 (2016) 322-328. [ Links ]

26. M. S. Hashemi, and D. Baleanu, Journal of Optoelectronics and Advanced Materials, 18 (2016). [ Links ]

27. N. Kudryashov, Phys Lett A 342 (2005) 99-106. [ Links ]

28. N. Kudryashov, Chaos Soliton Fract, 24 (2005) 1217-1231. [ Links ]

i In this paper, the case k = 2 has been used to find solutions.

Received: January 17, 2019; Accepted: February 19, 2019

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License