AAA Abstract and Applied Analysis 1687-0409 1085-3375 Hindawi Publishing Corporation 465160 10.1155/2013/465160 465160 Research Article Numerical Solutions of Fractional Fokker-Planck Equations Using Iterative Laplace Transform Method Yan Limei http://orcid.org/0000-0002-8921-6662 Baleanu Dumitru School of Mathematical Sciences Dezhou University Dezhou 253023 China dzu.edu.cn 2013 4 12 2013 2013 03 08 2013 17 10 2013 02 11 2013 2013 Copyright © 2013 Limei Yan. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

A relatively new iterative Laplace transform method, which combines two methods; the iterative method and the Laplace transform method, is applied to obtain the numerical solutions of fractional Fokker-Planck equations. The method gives numerical solutions in the form of convergent series with easily computable components, requiring no linearization or small perturbation. The numerical results show that the approach is easy to implement and straightforward when applied to space-time fractional Fokker-Planck equations. The method provides a promising tool for solving space-time fractional partial differential equations.

1. Introduction

Fractional calculus has attracted much attention for its potential applications in various scientific fields such as fluid mechanics, biology, viscoelasticity, engineering, and other areas of science . So it becomes important to find some efficient methods for solving fractional differential equations. A great deal of effort has been spent on constructing of the numerical solutions and many effective methods have been developed such as fractional wavelet method , fractional differential transform method , fractional operational matrix method [10, 11], fractional improved homotopy perturbation method [12, 13], fractional variational iteration method [14, 15], and fractional Laplace Adomian decomposition method [16, 17].

In 2006, Daftardar-Gejji and Jafari proposed a new iterative method to seek numerical solutions of nonlinear functional equations [18, 19]. By now, the iterative method has been used to solve many nonlinear differential equations of integer and fractional order  and fractional boundary value problem . Most recently, Jafari et al. firstly applied Laplace transform in the iterative method and proposed a new direct method called iterative Laplace transform method  to search for numerical solutions of a system of fractional partial differential equations. The method is based on Laplace transform, iterative method, Caputo fractional derivative, and symbolic computation. By using this method, Jafari and Seifi successfully obtained the numerical solutions of two systems of space-time fractional differential equations. It has been shown that, with this method, one can discover some solutions found by the existing methods such as homotopy perturbation method, Laplace Adomian decomposition method, and variational iterative method .

It is well known that the choice of an appropriate ansatz is of great importance when a method is applied to search for numerical solutions of nonlinear partial differential equations. In the present paper, we will use the iterative Laplace transform method to solve space-time fractional Fokker-Planck equations. The fractional derivatives described here are in the Caputo sense.

Fokker-Planck equation has been applied in various natural science fields such as quantum optics, solid-state physics, chemical physics, theoretical biology, and circuit theory. It is firstly proposed by Fokker and Planck to characterize the Brownian motion of particles . The general form of Fokker-Planck equation is as follows: (1)ut=[-xA(x)+2x2B(x)]u(x,t) with the initial condition (2)u(x,0)=f(x), where A(x),B(x)>0 are called diffusion coefficient and drift coefficient, respectively. If A(x),B(x)>0 depend on u(x,t) and the time t, then (1) becomes the following generalized nonlinear form : (3)ut=[-xA(x,t,u)+2x2B(x,t,u)]u(x,t).

The space-time fractional Fokker-Planck is as follows : (4)Dtαu=[-DxβA(x,t,u)+Dx2βB(x,t,u)]u(x,t) which is the generalized fractional form of (3). Here Dtα(·),Dxβ(·),Dx2β(·) are the Caputo fractional derivative with respect to t and x defined in Section 2. When α=β=1, (4) reduces to (3).

The rest of this paper is organized as follows. In Section 2, we review some basic definitions of Caputo fractional derivative and Laplace transform. In Section 3, we describe the iterative Laplace transform method for solving fractional partial differential equations. In Section 4, we give three applications of the method to Fokker-Planck equations. In Section 5, some conclusions and discussions are given.

2. Preliminaries Definition 1.

The Caputo fractional derivative [27, 28] of function u(x,t) is defined as (5)Dtαu(x,t)=1Γ(m-α)0t(t-η)m-α-1u(m)(x,η)dη,hhhhhhhhhhhhhhhhhhhm-1<αm,mN, where Γ(·) denotes the gamma function.

Definition 2.

The Laplace transform of f(t) is defined as [27, 28] (6)F(s)=[f(t)]=0e-stf(t)dt.

Definition 3.

Laplace transform of Dtαu(x,t) is given as [27, 28] (7)[Dtαu(x,t)]=sα[u(x,t)]-k=0m-1u(k)(x,0)sα-1-k,hhhhhhhhhhhhhhhhhhhhhhhhhim-1<αm, where u(k)(x,0) is the k-order derivative of u(x,t) at t=0.

Further information about fractional derivative and its properties can be found in .

3. The Iterative Laplace Transform Method

To illustrate the basic idea of the iterative Laplace transform method , we consider the general space-time fractional partial differential equation with initial conditions of the form (8)Dtαu=𝒜(u,Dxβu,Dx2βu,),m-1<αm,n-1<βn,m,nN, with initial value conditions (9)u(k)(x,0)=hk(x),k=0,1,,m-1, where 𝒜(u,Dxβu,Dx2βu,) is a linear or nonlinear operator of u,Dxβu,Dx2βu,, and u=u(x,t) is the unknown function that will be determined later.

Taking Laplace transfer of both sides of (8) results in (10)sα[u(x,t)]-k=0m-1sα-1-ku(k)(x,0)=[𝒜(u,Dxβu,Dx2βu,)]. Equivalently, (11)[u(x,t)]=k=0m-1s-1-ku(k)(x,0)+s-α[𝒜(u,Dxβu,Dx2βu,)].

Operating with Laplace inverse (denoted by -1 throughout the present paper) on both sides of (11) gives (12)u(x,t)=-1[k=0m-1s-1-ku(k)(x,0)]+-1[s-α[𝒜(u,Dxβu,Dx2βu,)]], which can be rewritten as the form (13)u(x,t)=-1[k=0m-1s-1-ku(k)(x,0)]+(u,Dxβu,Dx2βu,), where (u,Dxβu,Dx2βu,)=-1[s-α[𝒜(u,Dxβu,Dx2βu,)]].

The iterative Laplace transform method represents the solution as an infinite series: (14)u(x,t)=n=0un, where the terms un are to be recursively computed. The linear or nonlinear operator (u,Dxβu,Dx2βu,) can be decomposed as follows: (15)(n=0un,Dxβn=0un,Dx2βn=0un,)=(u0,Dxβu0,Dx2βu0,)+j=1(k=0juk,Dxβk=0juk,Dx2βk=0juk,)-j=1(k=0j-1uk,Dxβk=0j-1uk,Dx2βk=0j-1uk,). Substituting (14) and (15) into (13) yields (16)n=0un=-1[k=0m-1s-1-ku(k)(x,0)]+(u0,Dxβu0,Dx2βu0,)+j=1[(k=0juk,Dxβk=0juk,Dx2βk=0juk,)-(k=0j-1uk,Dxβk=0j-1uk,Dx2βk=0j-1uk,)].

We set (17)u0=-1[k=0m-1s-1-ku(k)(x,0)],u1=(u0,Dxβu0,Dx2βu0,),um+1=(k=0muk,Dxβk=0muk,Dx2βk=0muk,)-(k=0m-1uk,Dxβk=0m-1uk,Dx2βk=0m-1uk,),m1.

Therefore the m-term numerical solution of (8)-(9) is given by (18)u(x,t)u0(x,t)+u1(x,t)++um(x,t),m=1,2,.

4. Numerical Solutions of Fractional Fokker-Planck Equations

The iterative Laplace transform method, described in Section 3, will be applied to solve three special cases of space-time fractional Fokker-Planck equations with initial conditions.

Example 1.

Consider the Fokker-Planck equation in the case that  (19)Dtαu=ux+2ux2,0<α1, subject to (20)u(x,0)=x.

Taking Laplace transform on both sides of (19) gives (21)sα[u(x,t)]-sα-1u(x,0)=[ux+2ux2],(22)[u(x,t)]=xs+1sα[ux+2ux2].

Operating with Laplace inverse on both sides of (22) results in (23)u(x,t)=x+-1[1sα[ux+2ux2]].

Substituting (14) and (15) into (23) and applying (17), we obtain the components of the solution as follows: (24)u0(x,t)=u(x,0)=x,u1(x,t)=-1[1sα[u0x+2u0x2]]=tαΓ(1+α),u2(x,t)=0,u3(x,t)==un(x,t)==0.

Therefore, the solution of (19)-(20) in a closed form can be obtained as follows: (25)u(x,t)=x+tαΓ(1+α).

If we take α=1, then (25) can be reduced to (26)u(x,t)=x+t, which is exactly the same as that obtained by homotopy perturbation method in .

It should be pointed out that the iterative Laplace transform method is the generalization algorithm of iterative method proposed by Daftardar-Gejji and Jafari . When these two methods are used to solve differential equations with integer order derivatives, especially for linear cases, they are not different from each other.

Example 2.

Consider the following space-time fractional Fokker-Planck equation with initial value conditions : (27)Dtαu=-Dxβ(ux6)+Dx2β(ux212),t>0,x>0,0<α,β1, subject to (28)u(x,0)=x2.

Taking Laplace transform on both sides of (27) gives (29)sα[u(x,t)]-sα-1u(x,0)=[-Dxβ(ux6)+Dx2β(ux212)],(30)[u(x,t)]=u(x,0)s+1sα[-Dxβ(ux6)+Dx2β(ux212)].

Operating with Laplace inverse on both sides of (30), we obtain the following Laplace equation: (31)u(x,t)=u(x,0)+-1[1sα[-Dxβ(ux6)+Dx2β(ux212)]]. Following the algorithm given in (17), the first three components of the solution are as follows: (32)u0(x,t)=u(x,0)=x2,u1(x,t)=-1[1sα[-Dxβ(u0x6)+Dx2β(u0x212)]]=(-1Γ(4-β)x3-β+2Γ(5-2β)x4-2β)·tαΓ(1+α),u2(x,t)=[4-β6Γ(5-2β)x4-2β-5-3β3Γ(6-3β)x5-3β-(5-β)(4-β)12Γ(6-3β)x5-3β+(6-2β)(5-2β)6Γ(7-4β)×x6-4β4-β6Γ(5-2β)]·t2αΓ(1+2α).

The solution in series form is then given by (33)u(x,t)=x2+(-1Γ(4-β)x3-β+2Γ(5-2β)x4-2β)·tαΓ(1+α)+[4-β6Γ(5-2β)x4-2β-5-3β3Γ(6-3β)x5-3β-(5-β)(4-β)12Γ(6-3β)x5-3β+(6-2β)(5-2β)6Γ(7-4β)x6-4β]·t2αΓ(1+2α)+.

Setting α=β=1 in (33), we get the solution of the problem by (34)u(x,t)=x2(1+t2+(t/2)22!+) and in a closed form by (35)u(x,t)=x2et/2 which is in full agreement with the results by homotopy perturbation method in .

Example 3.

Consider the following space-time fractional nonlinear initial value problem that describes Fokker-Planck equation : (36)Dtαu=-Dxβ(4u2x-xu3)+Dx2βu2,0<α,β1, subject to (37)u(x,0)=x2.

Taking Laplace transform on both sides of (35) gives (38)sα[u(x,t)]-sα-1u(x,0)=[-Dxβ(4u2x-xu3)+Dx2βu2],(39)[u(x,t)]=x2s+1sα[-Dxβ(4u2x-xu3)+Dx2βu2].

Operating with Laplace inverse on both sides of (38), we obtain the following Laplace equation: (40)u(x,t)=x2+-1[1sα[-Dxβ(4u2x-xu3)+Dx2βu2]].

Following the algorithm given in (17), the first few components of the solution are as follows: (41)u0(x,t)=x2,u1(x,t)=-1[1sα[-Dxβ(4u02x-xu03)+Dx2βu02]]=(-22Γ(4-β)x3-β+24Γ(5-2β)x4-2β)·tαΓ(1+α),u2(x,t)=-1[(4(u0+u1)2x-x(u0+u1)3)1sα[-Dxβ(4(u0+u1)2x-x(u0+u1)3)+Dx2β(u0+u1)2(4(u0+u1)2x-x(u0+u1)3)]]--1[1sα[-Dxβ(4u02x-xu03)+Dx2βu02]]=-Γ(1+2α)Γ(8-4β)Γ2(1+α)Γ(1+3α)Γ(8-5β)Γ(5-2β)·(2304Γ(5-2β)+1056Γ(4-β))t3αx7-5β-1Γ(1+2α)Γ(6-3β)(184Γ(6-2β)Γ(5-2β)+44Γ(6-β)Γ(4-β))·t2αx5-3β+Γ(1+2α)Γ2(1+α)Γ(1+3α)Γ(7-4β)Γ(4-β)×(484Γ(7-2β)Γ(4-β)+4224Γ(7-3β)Γ(5-2β))·t3αx6-4β+506Γ(5-β)3Γ(1+2α)Γ(5-2β)Γ(4-β)·t2αx4-2β-1936Γ(1+2α)Γ(6-2β)Γ2(1+α)Γ(1+3α)Γ(6-3β)Γ2(4-β)·t3αx5-3β+576Γ(1+2α)Γ(9-4β)Γ2(1+α)Γ(1+3α)Γ(8-5β)Γ2(5-2β)·t3αx8-6β+48Γ(7-2β)Γ(1+2α)Γ(7-4β)Γ(5-2β)·t2αx6-4β+(22Γ(4-β)x3-β-24Γ(5-2β)x4-2β)·tαΓ(1+α), The solution in series form is then given by (42)u(x,t)=u0(x,t)+u1(x,t)+u2(x,t)+. If we take α=β=1, the first few components of the solution are as follows: (43)u0(x,t)=x2,u1(x,t)=x2t,u2(x,t)=x2t22,

For this special case, the exact solution of (36) and (37) is therefore given by (44)u(x,t)=x2et which is exactly the result obtained by homotopy perturbation transformation method in .

Table 1 shows the numerical solutions for (36) and (37) by using iterative Laplace transform method, homotopy perturbation transform method, Adomian decomposition method, and the exact solution as α=β=1. It should be pointed out that only three terms of these methods are used to evaluate the numerical solutions in Table 1. It is obvious that the iterative Laplace transform method used in the present paper has the same convergence as the convergence of homotopy perturbation transform method and Adomian decomposition method for solving this fractional nonlinear Fokker-Planck equation. Therefore, iterative Laplace transform method is an effective method for solving fractional partial differential equations just as homotopy perturbation transform method and Adomian decomposition method.

Several approximate values and exact solutions for (36) and (37) when α=β=1.

x t Solutio n ILTM Solutio n HPM Solutio n ADM Exact solution
0.25 0.1 0.0690729 0.0690729 0.0690729 0.0690732
0.2 0.0763333 0.0763333 0.0763333 0.0763377
0.4 0.0931667 0.0931667 0.0931667 0.093239
0.8 0.137833 0.137833 0.137833 0.139096

0.5 0.1 0.276292 0.276292 0.276292 0.276293
0.2 0.305333 0.305333 0.305333 0.305351
0.4 0.372667 0.372667 0.372667 0.372956
0.8 0.551333 0.551333 0.551333 0.556385

0.75 0.1 0.621656 0.621656 0.621656 0.621659
0.2 0.687 0.687 0.687 0.687039
0.4 8385 0.8385 0.8385 0.839151
0.8 1.2405 1.2405 1.2405 1.25187
5. Conclusions

With the aid of the symbolic computation system Mathematica, the iterative Laplace transform method is first successfully applied to solve fractional Fokker-Planck equations. The results obtained by the iterative Laplace transform method are the same as those obtained by homotopy perturbation transform method and Adomian decomposition method. The method finds the solutions without unnecessary linearization, small perturbation and other restrictive assumptions. Therefore, the method considerably reduces the computational work to a great extent. It is worth mentioning that the method can also be applied to solve other nonlinear fractional differential equations with initial value conditions.

Acknowledgments

The author would like to express her sincere thanks to the referees for their valuable suggestions which lead to an improved version. This research is supported by Excellent Young Scientist Foundation of Shandong Province under Grant no. BS2013HZ026.

Sabatier J. Agrawal O. P. Tenreiro Machado J. A. Advances in Fractional Calculus: Theoretical Developments and Applications in Physics and Engineering 2007 Springer 10.1007/978-1-4020-6042-7 MR2432163 Baleanu D. Diethelm K. Scalas E. Trujillo J. J. Fractional Calculus 2012 3 Singapore World Scientific Series on Complexity, Nonlinearity and Chaos 10.1142/9789814355216 MR2894576 Liu Y. Xin B. Numerical solutions of a fractional predator-prey system Advances in Difference Equations 2011 2011 11 190475 MR2780665 ZBL1217.35205 Ma J. Liu Y. Exact solutions for a generalized nonlinear fractional Fokker-Planck equation Nonlinear Analysis: Real World Applications 2010 11 1 515 521 10.1016/j.nonrwa.2009.01.006 MR2570570 ZBL1181.35293 Rehman M. Ali Khan R. The Legendre wavelet method for solving fractional differential equations Communications in Nonlinear Science and Numerical Simulation 2011 16 11 4163 4173 10.1016/j.cnsns.2011.01.014 MR2806727 ZBL1222.65063 Li Y. Solving a nonlinear fractional differential equation using Chebyshev wavelets Communications in Nonlinear Science and Numerical Simulation 2010 15 9 2284 2292 10.1016/j.cnsns.2009.09.020 MR2602712 ZBL1222.65087 Wu J. L. A wavelet operational method for solving fractional partial differential equations numerically Applied Mathematics and Computation 2009 214 1 31 40 10.1016/j.amc.2009.03.066 MR2537437 ZBL1169.65127 Lepik Ü. Solving fractional integral equations by the Haar wavelet method Applied Mathematics and Computation 2009 214 2 468 478 10.1016/j.amc.2009.04.015 MR2541683 ZBL1170.65106 Ertürk V. S. Momani S. Solving systems of fractional differential equations using differential transform method Journal of Computational and Applied Mathematics 2008 215 1 142 151 10.1016/j.cam.2007.03.029 MR2400624 ZBL1141.65088 Li Y. Sun N. Numerical solution of fractional differential equations using the generalized block pulse operational matrix Computers & Mathematics with Applications 2011 62 3 1046 1054 10.1016/j.camwa.2011.03.032 MR2824692 ZBL1228.65135 Saadatmandi A. Dehghan M. A new operational matrix for solving fractional-order differential equations Computers & Mathematics with Applications 2010 59 3 1326 1336 10.1016/j.camwa.2009.07.006 MR2579495 ZBL1189.65151 Liu Y. Approximate solutions of fractional nonlinear equations using homotopy perturbation transformation method Abstract and Applied Analysis 2012 2012 14 752869 10.1155/2012/752869 MR2910714 ZBL1242.65136 Liu Y. Variational homotopy perturbation method for solving fractional initial boundary value problems Abstract and Applied Analysis 2012 2012 10 727031 MR2926899 ZBL1246.65191 Sweilam N. H. Khader M. M. Al-Bar R. F. Numerical studies for a multi-order fractional differential equation Physics Letters A 2007 371 1-2 26 33 10.1016/j.physleta.2007.06.016 MR2419274 ZBL1209.65116 Das S. Analytical solution of a fractional diffusion equation by variational iteration method Computers & Mathematics with Applications 2009 57 3 483 487 10.1016/j.camwa.2008.09.045 MR2488620 ZBL1165.35398 Jafari H. Khalique C. M. Nazari M. Application of the Laplace decomposition method for solving linear and nonlinear fractional diffusion-wave equations Applied Mathematics Letters 2011 24 11 1799 1805 10.1016/j.aml.2011.04.037 MR2812215 ZBL1231.65179 Ongun M. Y. The Laplace adomian decomposition method for solving a model for HIV infection of CD4+ cells Mathematical and Computer Modelling 2011 53 5-6 597 603 10.1016/j.mcm.2010.09.009 MR2769431 ZBL1217.65164 Daftardar-Gejji V. Jafari H. An iterative method for solving nonlinear functional equations Journal of Mathematical Analysis and Applications 2006 316 2 753 763 10.1016/j.jmaa.2005.05.009 MR2207344 ZBL1087.65055 Jafari H. Iterative Methods for solving system of fractional differential equations [Ph.D. thesis] 2006 Pune University ZBL1087.65055 Bhalekar S. Daftardar-Gejji V. Solving evolution equations using a new iterative method Numerical Methods for Partial Differential Equations 2010 26 4 906 916 10.1002/num.20463 MR2642326 ZBL1194.65117 Daftardar-Gejji V. Bhalekar S. Solving fractional boundary value problems with Dirichlet boundary conditions using a new iterative method Computers & Mathematics with Applications 2010 59 5 1801 1809 10.1016/j.camwa.2009.08.018 MR2595954 ZBL1189.35357 Jafari H. Nazari M. Baleanu D. Khalique C. M. A new approach for solving a system of fractional partial differential equations Computers & Mathematics with Applications 2013 66 5 838 843 10.1016/j.camwa.2012.11.014 MR3089391 Jafari H. Seifi S. Solving a system of nonlinear fractional partial differential equations using homotopy analysis method Communications in Nonlinear Science and Numerical Simulation 2009 14 5 1962 1969 10.1016/j.cnsns.2008.06.019 MR2474456 ZBL1221.35439 Risken H. The Fokker-Planck Equation, Methods of Solution and Applications 1989 Berlin, Germany Springer 10.1007/978-3-642-61544-3 MR987631 Frank T. D. Stochastic feedback, nonlinear families of Markov processes, and nonlinear Fokker-Planck equations Physica A 2004 331 3-4 391 408 10.1016/j.physa.2003.09.056 MR2044918 Yiidirim A. Analytical approach to Fokker-Planck equation with space- and time fractional derivatives by means of the homotopy perturbation method Journal of King Saud University 2010 22 4 257 264 10.1016/j.jksus.2010.05.008 Miller K. S. Ross B. An Introduction to the Fractional Calculus and Fractional Differential Equations 1993 New York, NY, USA John Wiley & Sons MR1219954 Podlubny I. Fractional Differential Equations 1999 198 New York, NY, USA Academic Press MR1658022 Hilfer R. Applications of Fractional Calculus in Physics 2000 Singapore World Scientific 10.1142/9789812817747 MR1890104 Biazar J. Hosseini K. Gholamin P. Homotopy perturbation method Fokker-Planck equation International Mathematical Forum 2008 3 19 945 954 MR2403410 ZBL1165.82020