- About this Journal ·
- Abstracting and Indexing ·
- Advance Access ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents
Abstract and Applied Analysis
Volume 2014 (2014), Article ID 282190, 10 pages
Numerical Treatment of the Modified Time Fractional Fokker-Planck Equation
School of Mathematics and Statistics, Tianshui Normal University, Tianshui 741001, China
Received 26 December 2013; Revised 8 February 2014; Accepted 16 February 2014; Published 27 March 2014
Academic Editor: Adem Kilicman
Copyright © 2014 Yuxin Zhang. 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 numerical method for the modified time fractional Fokker-Planck equation is proposed. Stability and convergence of the method are rigorously discussed by means of the Fourier method. We prove that the difference scheme is unconditionally stable, and convergence order is , where and are the temporal and spatial step sizes, respectively. Finally, numerical results are given to confirm the theoretical analysis.
Fractional differential equations have attracted considerable attention due to their frequent appearance in various applications in fluid mechanics, biology, physics, and engineering [1, 2]. Usually, fractional differential equations do not have analytic solutions and can only be solved by some semianalytical and numeric techniques. Recently, several semianalytical methods, such as variational iteration method, homotopy perturbation method, Adomian decomposition method, homotopy analysis method, and collocation method, have been used to solve fractional differential equations [3–7]. Meanwhile, some effective numerical techniques are developed; see [8–17].
In the present paper, we are motivated to study the following modified time fractional Fokker-Planck equation : where is a fractional diffusion coefficient, and is a fractional friction coefficient. denotes the temporal Riemann-Liouville derivative operator defined as  where is the gamma function.
The outline of the paper is as follows. In Section 2, an effective numerical method for solving the modified time fractional Fokker-Planck equation is proposed. The solvability, stability, and convergence of the numerical method are discussed in Sections 3 and 4, respectively. In Section 5, we give some numerical results demonstrating the convergence orders of the numerical method. Also a conclusion is given in Section 6.
2. The Construction of Numerical Method
Let and , where and are the uniform spatial and temporal mesh sizes, respectively, and , are two positive integers.
Lemma 1. Suppose that then a fourth-order difference scheme for the above equation is given as follows: where and are two difference operators and are defined by in which is an unit operator and and are average central and second central difference operators with respect to and are defined by
Proof. In view of Taylor expansion, we can obtain
Noting (4), we easily obtain This completes the proof.
Using the relation of the Riemann-Liouville fractional derivative and Grünwald-Letnikov fractional derivative , we can approximate the Riemann-Liouville fractional derivative by  where .
For first-order derivative , we apply the following backward difference scheme:
The discrete form of above system is where , , , and .
The initial and boundary conditions can be discretized as
Obviously, the local truncation error of difference scheme (14) is .
3. The Solvability of the Difference Scheme
Firstly, let us denote
Then we can give the compact form of difference scheme (14) as follows: where
Theorem 2. Difference equation (17) is uniquely solvable.
Proof. It is well known that the eigenvalues of the matrix are
Note that and ; if then we easily know that .
At the moment, we obtain ; that is to say, the matrix is invertible. Hence, difference equation (17) has a unique solution.
4. Stability and Convergence Analysis
Lemma 3. The coefficients satisfy  as follows:(1), , , ;(2) .
Let be the approximate solution of (14) and define respectively.
So, we can easily obtain the following roundoff error equation: Now, we define the grid functions then can be expanded in a Fourier series: where
We introduce the following norm: and according to the Parseval equality we obtain Through the above analysis, we can suppose that the solution of (25) has the following form: where .
Substituting the above expression into (25) one gets
Lemma 4. The following relation holds: where
Proof. Because of
Furthermore, we can rewrite the above inequality as that is, This completes the proof of Lemma 4.
Lemma 5. Supposing that is the solution of (33), then we have
Lemma 6. Difference scheme (14) is unconditionally stable.
Proof. According to Lemma 5, we obtain
which means that difference scheme (14) is unconditionally stable.
Next, we give the convergence analysis. Suppose and denote
From (14), we obtain
Similar to the stability analysis method, we define the grid functions then and can be expanded to the following Fourier series, respectively: where The same as before, we also have Based on the above analysis, we can assume that and have the following forms: respectively. Substituting the above two expressions into (48) yields
Lemma 7. Let be the solution of (55); then there exists a positive constant , so that
Proof. From , we have
In view of the convergence of the series of (53), there is a positive constant , such that
For , from (55), we have
Noticing (58), then
Now, we suppose that Then when , we obtain This completes the proof.
Theorem 8. Difference scheme (14) is convergent, and the convergence order is .
Remark 9. From above discussion, we know that difference scheme is an implicit scheme and it is unconditionally stable and convergent. If we take in (1), then we can obtain an explicit scheme and it is conditionally stable and convergent.
5. Numerical Example
In this section, a numerical example is presented to confirm our theoretical results.
Example 10. Consider the following equation:
with the initial and boundary conditions
where . The analytical solution of this equation is .
The maximum error, temporal, and spatial convergence orders by difference scheme (14) for various are listed in Table 1. From the obtained results, we can draw the following conclusions: the experimental convergence orders are approximately 1 and 4 in temporal and spatial directions, respectively. Figures 1 and 2 show the comparison of the numerical solution with the analytical solution at and for different temporal and spatial mesh sizes.
By Table 1 and Figures 1 and 2, it can be seen that the numerical solution is in excellent agreement with the analytical solution. These results confirm our theoretical analysis.
In this paper, a computationally effective numerical method is proposed for simulating the modified time fractional Fokker-Planck equation. It has proven the unconditional stability and solvability of proposed scheme. Also we showed that the method is convergent with order . The numerical results demonstrate the effectiveness of the proposed scheme.
Conflict of Interests
The author declares that there is no conflict of interests regarding the publication of this paper.
The work was partially supported by Tianshui Normal University “QingLan” Talent Engineering Funds.
- K. S. Miller and B. Ross, An Introduction To the Fractional Calculus and Fractional Differential Equations, John Wiley, New York, NY, USA, 1993.
- I. Podlubny, Fractional Differential Equations, Academic Press, 1999.
- A. Barari, M. Omidvar, A. R. Ghotbi, and D. D. Ganji, “Application of homotopy perturbation method and variational iteration method to nonlinear oscillator differential equations,” Acta Applicandae Mathematicae, vol. 104, no. 2, pp. 161–171, 2008.
- S. Das, “Analytical solution of a fractional diffusion equation by variational iteration method,” Computers and Mathematics with Applications, vol. 57, no. 3, pp. 483–487, 2009.
- D. D. Ganji and A. Sadighi, “Application of homotopy-perturbation and variational iteration methods to nonlinear heat transfer and porous media equations,” Journal of Computational and Applied Mathematics, vol. 207, no. 1, pp. 24–34, 2007.
- E. Hesameddini and F. Fotros, “Solution for time-fractional coupled Klein-Gordon Schrodinger equation using decomposition method,” International Mathematics Olympiad, vol. 7, pp. 1047–1056, 2012.
- G.-C. Wu and E. W. M. Lee, “Fractional variational iteration method and its application,” Physics Letters A, vol. 374, no. 25, pp. 2506–2509, 2010.
- C.-M. Chen, F. Liu, and K. Burrage, “Finite difference methods and a fourier analysis for the fractional reaction-subdiffusion equation,” Applied Mathematics and Computation, vol. 198, no. 2, pp. 754–769, 2008.
- C. Çelik and M. Duman, “Crank-Nicolson method for the fractional diffusion equation with the Riesz fractional derivative,” Journal of Computational Physics, vol. 231, no. 4, pp. 1743–1750, 2012.
- X. Hu and L. Zhang, “Implicit compact difference schemes for the fractional cable equation,” Applied Mathematical Modelling, vol. 36, pp. 4027–4043, 2012.
- M. M. Meerschaert and C. Tadjeran, “Finite difference approximations for two-sided space-fractional partial differential equations,” Applied Numerical Mathematics, vol. 56, no. 1, pp. 80–90, 2006.
- D. K. Salkuyeh, “On the finite difference approximation to the convection-diffusion equation,” Applied Mathematics and Computation, vol. 179, no. 1, pp. 79–86, 2006.
- E. Sousa, “Numerical approximations for fractional diffusion equations via splines,” Computers and Mathematics with Applications, vol. 62, no. 3, pp. 938–944, 2011.
- E. Sousa, “Finite difference approximations for a fractional advection diffusion problem,” Journal of Computational Physics, vol. 228, no. 11, pp. 4038–4054, 2009.
- C. Tadjeran and M. M. Meerschaert, “A second-order accurate numerical method for the two-dimensional fractional diffusion equation,” Journal of Computational Physics, vol. 220, no. 2, pp. 813–823, 2007.
- S. B. Yuste and L. Acedo, “An explicit finite difference method and a new von Neumann-type stability analysis for fractional diffusion equations,” SIAM Journal on Numerical Analysis, vol. 42, no. 5, pp. 1862–1874, 2005.
- Q. Yang, F. Liu, and I. Turner, “Numerical methods for fractional partial differential equations with Riesz space fractional derivatives,” Applied Mathematical Modelling, vol. 34, no. 1, pp. 200–218, 2010.
- B. I. Henry, T. A. M. Langlands, and P. Straka, “Fractional Fokker-Planck equations for subdiffusion with space- and time-dependent forces,” Physical Review Letters, vol. 105, no. 17, Article ID 170602, 2010.