Research Article  Open Access
Numerical Solutions of Fractional FokkerPlanck Equations Using Iterative Laplace Transform Method
Abstract
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 FokkerPlanck 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 spacetime fractional FokkerPlanck equations. The method provides a promising tool for solving spacetime 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 [1–4]. 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 [5–8], fractional differential transform method [9], 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, DaftardarGejji 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 [20] and fractional boundary value problem [21]. Most recently, Jafari et al. firstly applied Laplace transform in the iterative method and proposed a new direct method called iterative Laplace transform method [22] 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 spacetime 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 [23].
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 spacetime fractional FokkerPlanck equations. The fractional derivatives described here are in the Caputo sense.
FokkerPlanck equation has been applied in various natural science fields such as quantum optics, solidstate physics, chemical physics, theoretical biology, and circuit theory. It is firstly proposed by Fokker and Planck to characterize the Brownian motion of particles [24]. The general form of FokkerPlanck equation is as follows: with the initial condition where are called diffusion coefficient and drift coefficient, respectively. If depend on and the time , then (1) becomes the following generalized nonlinear form [25]:
The spacetime fractional FokkerPlanck is as follows [26]: which is the generalized fractional form of (3). Here are the Caputo fractional derivative with respect to and defined in Section 2. When , (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 FokkerPlanck equations. In Section 5, some conclusions and discussions are given.
2. Preliminaries
Definition 1. The Caputo fractional derivative [27, 28] of function is defined as where denotes the gamma function.
Definition 2. The Laplace transform of is defined as [27, 28]
Definition 3. Laplace transform of is given as [27, 28] where is the order derivative of at .
Further information about fractional derivative and its properties can be found in [27–29].
3. The Iterative Laplace Transform Method
To illustrate the basic idea of the iterative Laplace transform method [22], we consider the general spacetime fractional partial differential equation with initial conditions of the form with initial value conditions where is a linear or nonlinear operator of , and is the unknown function that will be determined later.
Taking Laplace transfer of both sides of (8) results in Equivalently,
Operating with Laplace inverse (denoted by throughout the present paper) on both sides of (11) gives which can be rewritten as the form where .
The iterative Laplace transform method represents the solution as an infinite series: where the terms are to be recursively computed. The linear or nonlinear operator can be decomposed as follows: Substituting (14) and (15) into (13) yields
We set
Therefore the term numerical solution of (8)(9) is given by
4. Numerical Solutions of Fractional FokkerPlanck Equations
The iterative Laplace transform method, described in Section 3, will be applied to solve three special cases of spacetime fractional FokkerPlanck equations with initial conditions.
Example 1. Consider the FokkerPlanck equation in the case that [30]
subject to
Taking Laplace transform on both sides of (19) gives
Operating with Laplace inverse on both sides of (22) results in
Substituting (14) and (15) into (23) and applying (17), we obtain the components of the solution as follows:
Therefore, the solution of (19)(20) in a closed form can be obtained as follows:
If we take , then (25) can be reduced to
which is exactly the same as that obtained by homotopy perturbation method in [30].
It should be pointed out that the iterative Laplace transform method is the generalization algorithm of iterative method proposed by DaftardarGejji and Jafari [18]. 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 spacetime fractional FokkerPlanck equation with initial value conditions [26]:
subject to
Taking Laplace transform on both sides of (27) gives
Operating with Laplace inverse on both sides of (30), we obtain the following Laplace equation:
Following the algorithm given in (17), the first three components of the solution are as follows:
The solution in series form is then given by
Setting in (33), we get the solution of the problem by
and in a closed form by
which is in full agreement with the results by homotopy perturbation method in [26].
Example 3. Consider the following spacetime fractional nonlinear initial value problem that describes FokkerPlanck equation [12]:
subject to
Taking Laplace transform on both sides of (35) gives
Operating with Laplace inverse on both sides of (38), we obtain the following Laplace equation:
Following the algorithm given in (17), the first few components of the solution are as follows:
The solution in series form is then given by
If we take , the first few components of the solution are as follows:
For this special case, the exact solution of (36) and (37) is therefore given by
which is exactly the result obtained by homotopy perturbation transformation method in [12].
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 . 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 FokkerPlanck 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.

5. Conclusions
With the aid of the symbolic computation system Mathematica, the iterative Laplace transform method is first successfully applied to solve fractional FokkerPlanck 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.
References
 J. Sabatier, O. P. Agrawal, and J. A. Tenreiro Machado, Advances in Fractional Calculus: Theoretical Developments and Applications in Physics and Engineering, Springer, 2007. View at: Publisher Site  MathSciNet
 D. Baleanu, K. Diethelm, E. Scalas, and J. J. Trujillo, Fractional Calculus, vol. 3 of Series on Complexity, Nonlinearity and Chaos, World Scientific, Singapore, 2012. View at: Publisher Site  MathSciNet
 Y. Liu and B. Xin, “Numerical solutions of a fractional predatorprey system,” Advances in Difference Equations, vol. 2011, Article ID 190475, 11 pages, 2011. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 J. Ma and Y. Liu, “Exact solutions for a generalized nonlinear fractional FokkerPlanck equation,” Nonlinear Analysis: Real World Applications, vol. 11, no. 1, pp. 515–521, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. Rehman and R. Ali Khan, “The Legendre wavelet method for solving fractional differential equations,” Communications in Nonlinear Science and Numerical Simulation, vol. 16, no. 11, pp. 4163–4173, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Y. Li, “Solving a nonlinear fractional differential equation using Chebyshev wavelets,” Communications in Nonlinear Science and Numerical Simulation, vol. 15, no. 9, pp. 2284–2292, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J. L. Wu, “A wavelet operational method for solving fractional partial differential equations numerically,” Applied Mathematics and Computation, vol. 214, no. 1, pp. 31–40, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Ü. Lepik, “Solving fractional integral equations by the Haar wavelet method,” Applied Mathematics and Computation, vol. 214, no. 2, pp. 468–478, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 V. S. Ertürk and S. Momani, “Solving systems of fractional differential equations using differential transform method,” Journal of Computational and Applied Mathematics, vol. 215, no. 1, pp. 142–151, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Y. Li and N. Sun, “Numerical solution of fractional differential equations using the generalized block pulse operational matrix,” Computers & Mathematics with Applications, vol. 62, no. 3, pp. 1046–1054, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 A. Saadatmandi and M. Dehghan, “A new operational matrix for solving fractionalorder differential equations,” Computers & Mathematics with Applications, vol. 59, no. 3, pp. 1326–1336, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Y. Liu, “Approximate solutions of fractional nonlinear equations using homotopy perturbation transformation method,” Abstract and Applied Analysis, vol. 2012, Article ID 752869, 14 pages, 2012. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Y. Liu, “Variational homotopy perturbation method for solving fractional initial boundary value problems,” Abstract and Applied Analysis, vol. 2012, Article ID 727031, 10 pages, 2012. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 N. H. Sweilam, M. M. Khader, and R. F. AlBar, “Numerical studies for a multiorder fractional differential equation,” Physics Letters A, vol. 371, no. 12, pp. 26–33, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. Das, “Analytical solution of a fractional diffusion equation by variational iteration method,” Computers & Mathematics with Applications, vol. 57, no. 3, pp. 483–487, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 H. Jafari, C. M. Khalique, and M. Nazari, “Application of the Laplace decomposition method for solving linear and nonlinear fractional diffusionwave equations,” Applied Mathematics Letters, vol. 24, no. 11, pp. 1799–1805, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. Y. Ongun, “The Laplace adomian decomposition method for solving a model for HIV infection of CD4^{+} cells,” Mathematical and Computer Modelling, vol. 53, no. 56, pp. 597–603, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 V. DaftardarGejji and H. Jafari, “An iterative method for solving nonlinear functional equations,” Journal of Mathematical Analysis and Applications, vol. 316, no. 2, pp. 753–763, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 H. Jafari, Iterative Methods for solving system of fractional differential equations [Ph.D. thesis], Pune University, 2006. View at: Zentralblatt MATH
 S. Bhalekar and V. DaftardarGejji, “Solving evolution equations using a new iterative method,” Numerical Methods for Partial Differential Equations, vol. 26, no. 4, pp. 906–916, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 V. DaftardarGejji and S. Bhalekar, “Solving fractional boundary value problems with Dirichlet boundary conditions using a new iterative method,” Computers & Mathematics with Applications, vol. 59, no. 5, pp. 1801–1809, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 H. Jafari, M. Nazari, D. Baleanu, and C. M. Khalique, “A new approach for solving a system of fractional partial differential equations,” Computers & Mathematics with Applications, vol. 66, no. 5, pp. 838–843, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 H. Jafari and S. Seifi, “Solving a system of nonlinear fractional partial differential equations using homotopy analysis method,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 5, pp. 1962–1969, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 H. Risken, The FokkerPlanck Equation, Methods of Solution and Applications, Springer, Berlin, Germany, 1989. View at: Publisher Site  MathSciNet
 T. D. Frank, “Stochastic feedback, nonlinear families of Markov processes, and nonlinear FokkerPlanck equations,” Physica A, vol. 331, no. 34, pp. 391–408, 2004. View at: Publisher Site  Google Scholar  MathSciNet
 A. Yiidirim, “Analytical approach to FokkerPlanck equation with space and time fractional derivatives by means of the homotopy perturbation method,” Journal of King Saud University, vol. 22, no. 4, pp. 257–264, 2010. View at: Publisher Site  Google Scholar
 K. S. Miller and B. Ross, An Introduction to the Fractional Calculus and Fractional Differential Equations, John Wiley & Sons, New York, NY, USA, 1993. View at: MathSciNet
 I. Podlubny, Fractional Differential Equations, vol. 198, Academic Press, New York, NY, USA, 1999. View at: MathSciNet
 R. Hilfer, Applications of Fractional Calculus in Physics, World Scientific, Singapore, 2000. View at: Publisher Site  MathSciNet
 J. Biazar, K. Hosseini, and P. Gholamin, “Homotopy perturbation method FokkerPlanck equation,” International Mathematical Forum, vol. 3, no. 19, pp. 945–954, 2008. View at: Google Scholar  Zentralblatt MATH  MathSciNet
Copyright
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.