Fractional Integral Transform and ApplicationView this Special Issue
Research Article | Open Access
Yuxin Zhang, "Numerical Treatment of the Modified Time Fractional Fokker-Planck Equation", Abstract and Applied Analysis, vol. 2014, Article ID 282190, 10 pages, 2014. https://doi.org/10.1155/2014/282190
Numerical Treatment of the Modified Time Fractional Fokker-Planck Equation
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.