- About this Journal ·
- Abstracting and Indexing ·
- 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 2013 (2013), Article ID 493406, 15 pages
Numerical Algorithms for the Fractional Diffusion-Wave Equation with Reaction Term
Department of Mathematics, Shanghai University, Shanghai 200444, China
Received 24 April 2013; Revised 24 June 2013; Accepted 24 June 2013
Academic Editor: Juan J. Trujillo
Copyright © 2013 Hengfei Ding and Changpin Li. 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.
Two numerical algorithms are derived to compute the fractional diffusion-wave equation with a reaction term. Firstly, using the relations between Caputo and Riemann-Liouville derivatives, we get two equivalent forms of the original equation, where we approximate Riemann-Liouville derivative by a second-order difference scheme. Secondly, for second-order derivative in space dimension, we construct a fourth-order difference scheme in terms of the hyperbolic-trigonometric spline function. The stability analysis of the derived numerical methods is given by means of the fractional Fourier method. Finally, an illustrative example is presented to show that the numerical results are in good agreement with the theoretical analysis.
The realm of fractional differential equations has drawn increasing attention and interest due to their important applications in biology, physics, chemistry, biochemistry, medicine, and engineering [1–8]. Generally speaking, the analytical solutions of most fractional differential equations are not easy, even impossible, to obtain, so seeking numerical solutions of these equations becomes more and more important. Till now, there have been some studies in this respect. For example, Cui  proposed a compact finite difference method for the fractional diffusion equation. Chen et al.  analyzed an implicit difference scheme for the subdiffusion equation and proved the unconditional stability and convergence. Li et al.  constructed some novel methods for the fractional calculus and fractional ordinary differential equation. Liu et al.  considered the numerical method and analytical technique for the modified anomalous subdiffusion equation with a nonlinear source term. Sousa  gave three finite difference schemes for a fractional advection diffusion problem. In , she furthermore presented an improved difference scheme to enhance the spatial accuracy. H. Wang and K. Wang  developed a fast finite difference method for the fractional diffusion equations. Yuste and Acedo  derived an explicit finite difference scheme and gave a new von Neumann-type stability analysis for the fractional diffusion equations. Roop and his coworkers pioneered in considering the fractional differential equation by using the finite element method in [17–19], where the temporal derivative is the classical derivative and the spatial derivative is the fractional derivative. Li et al.  used the finite difference/finite element mixed method to numerically solve the nonlinear subdiffusion and superdiffusion partial differential equation, where both temporal and spatial derivatives are fractional derivatives. Zhang et al.  obtained a numerical method for the symmetric space fractional partial differential equations by using Galerkin finite element method and a backward difference technique. In , Zheng et al. established a fully discrete approximation for the nonlinear spatial fractional Fokker-Planck equation, where the discontinuous Galerkin finite element approach was utilized in time domain and the Galerkin finite element approach was utilized in spatial domain. Yang et al. used the matrix transform method for the Riesz space fractional diffusion and advection-dispersion equations and the time and space fractional Fokker-Planck equations [23, 24], respectively. Recently, the spectral method was used to approximate the fractional calculus and the fractional differential equations [25, 26]. Very recently, Ding and Li  proposed two kinds of finite difference schemes for the reaction-subdiffusion equation.
As far as we know, there are some numerical methods for the subdiffusion equations [9, 10, 12, 16]; however, there seem to be very limited ones for the fractional diffusion-wave (seldom called “superdiffusion”) equation. The possible reason is that there may be no better method to analyze the stability and convergence. And in the studies available, the accuracy of the temporal direction for the subdiffusion and fractional diffusion-wave equations is often less than second order. This motivates us to construct some higher order methods.
In this paper, we study the fractional diffusion-wave equation with a reaction term in the following form:
Roughly speaking, the previous fractional differential equation is obtained from the classical diffusion-wave equation by replacing the second-order time derivative by a fractional derivative of order . It has been proposed to deal with disordered media to comb structures, dielectrics and semiconductors, and viscoelastic problems, for example, in the description of the propagation of stress waves in viscoelastic solids; see [28, 29] and references therein. In particular, in order to carry out the numerical comparisons, we study the following initial boundary values: where , , and are sufficiently smooth functions and and are diffusion and reaction coefficients, respectively. is the Caputo derivative of order , which will be introduced in the following.
The remainder of this paper is organized as follows. In Section 2, we present preliminary knowledge, which is necessary for our study. In Section 3, we develop two numerical methods. Meanwhile, we investigate the solvability and stability by the matrix method and fractional Fourier method, respectively. In Section 4, a numerical experiment is carried out which supports the theoretical analysis. Finally, the paper concludes with a summary in the last section.
2. Preliminary Knowledge
In the present section, we introduce some important definitions and lemmas which are used later on. The following definitions and lemmas can be found in .
Definition 1. The left Riemann-Liouville integral operator of order is defined in the following form: where is the Euler gamma function.
Definition 2. The left Caputo derivative operator of order is defined in the following form:
Definition 3. The left Riemann-Liouville derivative operator of order is defined by
Lemma 4. Suppose , for any , if and is integrable on with respect to , and then
Lemma 5 (see ). The eigenvalues of the tridiagonal matrix with order are given by
Let and , where and are the uniform spatial and temporal mesh sizes, respectively, and , are two positive integers.
We consider a uniform mesh with nodal points on , for each segment and point we define the following hyperbolic-trigonometric spline function which has different form with the mixed spline function in :
where , and are constants and is an arbitrary parameter which can improve the accuracy of the numerical method.
To obtain the expressions of the coefficients in (11), we denote
Through the straightforward calculations, we obtain the following expressions:
in which .
From (11), one has
In view of the continuity of the first-order derivative at the common nodes , that is, , we get the following relation: where
Obviously, from the expressions of , and , we find that .
When , then ; the spline function given by (15) reduces to the ordinary cubic spline function:
If we use and to take the place of ; and in (15), respectively, then one gets
Using Taylor’s series for (18) at , we get the truncation error as follows:
From the previous expression, we find that we can obtain different difference schemes by choosing different parameters , and . In this paper, we take , and , and then . Under this assumption, (18) can be rewritten as which is the same as in , where is the second-order central difference operator with respect to .
3. Numerical Algorithms
3.1. Numerical Algorithm I
The Riemann-Liouville derivative (with homogeneous initial value conditions) can be approximated by the following scheme:
There exist various choices of the generating functions which can lead to different approximation order .
Let be the generating function with coefficients ; that is,
If we take the generating function in the form that is, then we can obtain the following first-order scheme: in which and these coefficients have the following recursive relations:
If we take the generating function as that is, which leads to order , then
In the following, we give another simple but interesting method to obtain the coefficients .
and (26), one gets
Comparing the previous equation with (31), one obtains
which has the following recursive relation: In summary, we have the following lemma.
Lemma 6. The coefficients satisfy
Finally, let be the numerical approximation of ; substituting the expansions (20) and (32) in (22) and removing the higher order term, one obtains Numerical Algorithm as follows: for . The initial and boundary value conditions can be discretized by
where , and .
Next, we study the solvability and the local truncation error.
For convenience, denote
Then we can give the following matrix form of the difference scheme (38): where
Theorem 7. The difference equation (38) is uniquely solvable.
Proof. From Lemma 5, we know that the eigenvalues of the tridiagonal matrix are Due to , , and , hence, That is, the matrix is nonsingular. Thus (41) is uniquely solvable, and so is (38). The proof is finished.
Theorem 8. The local truncation error of the difference scheme (38) is .
In the following, we discuss the stability of Numerical Algorithm by the fractional Fourier method .
Let be the approximate solution of (38) and define
So, we can easily obtain the following round-off error equation:
We suppose that the solution of (48) has the following form: where is any of the real spatial wave numbers supported by the lattice. Substituting (49) into (48) gives If we write and assume that is independent of time, then we obtain the following expression of the amplification factor:
When for some , then the temporal factor of the solution grows to infinity according to (51); the numerical method is unstable. Considering the extreme value , we find that the numerical method is stable when Although depends on , we can estimate it by its limit value In this case, the stability condition becomes
3.2. Numerical Algorithm II
Differentiating (1) with order in the sense of Riemann-Liouville yields Due to it immediately follows that
Next, for the second-order time derivative , we approximate it by the usual three-point central difference formula
Let be the numerical approximation of . Substituting the expansions (20), (32), and (60) in (59) and removing the higher-order terms, one obtains Numerical Algorithm II as follows: where . The initial and boundary value conditions can be discretized by
where , and .
Similar to Numerical Algorithm I, we discuss the solvability and local truncation error of Numerical Algorithm II.
Denote Then the matrix form of Numerical Algorithm II (61) is written as where
Theorem 9. The difference equation (61) is uniquely solvable.
Proof. Obviously, the eigenvalues of the tridiagonal matrix are given by
Therefore, the solution of (61) exists and is uniquely solvable.
In order to get the local truncation error analysis, we need the following lemma.
Lemma 10. Consider
Theorem 11. The local truncation error of difference scheme (61) is .
Next, we study the stability of Numerical Algorithm II (see (61)).
Considering the extreme value , we can obtain the following stability condition from (78):
For large enough , we can estimate by its limit value
It follows that
Because the sine function is bounded by , from (81) one finds that Numerical Algorithm II is to be stable if , that is,
Furthermore, stability condition (82) can be rewritten as
where and are diffusion and reaction coefficients, respectively.
Moreover, from Theorems 8 and 11, we know that Numerical Algorithm I (see (38)) and Numerical Algorithm II (see (61)) are all consistent with the local truncation errors and because as . Therefore, according to Lax’s Equivalence Theorem , Numerical Algorithm I unconditionally converges at the same order , but Numerical Algorithm II converges at the same order under condition (83).
4. Numerical Example
In this section, we give a numerical example to demonstrate the efficiency of the derived numerical algorithms.
Example 1. Consider the following equation:
The analytical solution of this equation is . The initial conditions and the boundary conditions are satisfied with the exact solution given previously.
On one hand, we use Numerical Algorithm I (see (38)) to solve the previous equation. At first, we give the comparisons of the analytical and numerical solutions for different order, , , and in Figures 1 and 2, respectively, and we observe that the numerical solution is in line with the analytical solution and the velocity increases with the increase of order . Meanwhile, the errors are shown in Figure 3 with . Finally, the temporal and spatial convergence orders are listed in Table 1.
Figures 4 and 5 display that the velocity increases with the increase of order for different and . Figure 6 presents the error between the analytical solution and numerical solution of (60). Table 2 lists the temporal convergence orders.
Finally, we checked the stability condition given in (83) in the following way. Firstly, we choose , and , which satisfy stability condition (83), the comparison result of the numerical solution and analytical solution is shown in Figure 7. Secondly, we choose , and , and then these parameters do not satisfy stability condition (83); the comparison result of the numerical solution and analytical solution is shown in Figure 8. From Figures 7 and 8, we declare that stability condition (83) is valid.
In this paper, we construct two difference schemes for solving the fractional diffusion-wave equation with reaction term. It is proved that Numerical Algorithm I (38) is unconditionally stable and Numerical Algorithm II (61) is conditionally stable by using the fractional Fourier method. The local truncation errors of two difference schemes are both , which are the best results till now. Finally, the numerical results given in this paper show the effectiveness of the derived numerical algorithms.
The work was partially supported by the Key Program of Shanghai Municipal Education Commission under Grant no. 12ZZ084 and the grant of “The First-Class Discipline of Universities in Shanghai.”
- E. Barkai, R. Metzler, and J. Klafter, “From continuous time random walks to the fractional Fokker-Planck equation,” Physical Review E, vol. 61, no. 1, pp. 132–138, 2000.
- D. A. Benson, S. W. Wheatcraft, and M. M. Meerschaert, “Application of a fractional advection-dispersion equation,” Water Resources Research, vol. 36, no. 6, pp. 1403–1412, 2000.
- R. Hilfer, Applications of Fractional Calculus in Physics, World Scientific, Singapore, 2000.
- F. Mainardi, “Fractional calculus: some basic problems in continuum and statistical mechanics,” in Fractals and Fractional Calculus in Continuum Mechanics, A. Carpinteri and F. Mainardi, Eds., Springer, New York, NY, USA, 1997.
- K. S. Miller and B. Ross, An Introduction to the Fractional Calculus and Fractional Differential Equations, John Wiley & Sons, New York, NY, USA, 1993.
- A. I. Saichev and G. M. Zaslavsky, “Fractional kinetic equations: solutions and applications,” Chaos, vol. 7, no. 4, pp. 753–764, 1997.
- S. B. Yuste, L. Acedo, and K. Lindenberg, “Reaction front in an reactionsubdiffusion process,” Physical Review E, vol. 69, no. 3, part 2, Article ID 036126, 2004.
- S. B. Yuste and K. Lindenberg, “Subdiffusion-limited reactions,” Physical Review Letters, vol. 87, no. 11, Article ID 118301, 2001.
- M. Cui, “Compact finite difference method for the fractional diffusion equation,” Journal of Computational Physics, vol. 228, no. 20, pp. 7792–7804, 2009.
- C.-M. Chen, F. Liu, I. Turner, and V. Anh, “A Fourier method for the fractional diffusion equation describing sub-diffusion,” Journal of Computational Physics, vol. 227, no. 2, pp. 886–897, 2007.
- C. P. Li, A. Chen, and J. J. Ye, “Numerical approaches to fractional calculus and fractional ordinary differential equation,” Journal of Computational Physics, vol. 230, no. 9, pp. 3352–3368, 2011.
- F. Liu, C. Yang, and K. Burrage, “Numerical method and analytic technique of the modified anomalous subdiffusion equation with a nonlinear source term,” Journal of Computational and Applied Mathematics, vol. 231, no. 1, pp. 160–176, 2009.
- E. Sousa, “Finite difference approximations for a fractional advection diffusion problem,” Journal of Computational Physics, vol. 228, no. 11, pp. 4038–4054, 2009.
- E. Sousa, “Numerical approximations for fractional diffusion equations via splines,” Computers & Mathematics with Applications, vol. 62, no. 3, pp. 938–944, 2011.
- H. Wang and K. Wang, “An alternating-direction finite difference method for two-dimensional fractional diffusion equations,” Journal of Computational Physics, vol. 230, no. 21, pp. 7830–7839, 2011.
- 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.
- V. J. Ervin, N. Heuer, and J. P. Roop, “Numerical approximation of a time dependent, nonlinear, space-fractional diffusion equation,” SIAM Journal on Numerical Analysis, vol. 45, no. 2, pp. 572–591, 2007.
- G. J. Fix and J. P. Roop, “Least squares finite-element solution of a fractional order two-point boundary value problem,” Computers & Mathematics with Applications, vol. 48, no. 7-8, pp. 1017–1033, 2004.
- J. P. Roop, “Computational aspects of FEM approximation of fractional advection dispersion equations on bounded domains in ,” Journal of Computational and Applied Mathematics, vol. 193, no. 1, pp. 243–268, 2006.
- C. P. Li, Z. G. Zhao, and Y. Q. Chen, “Numerical approximation of nonlinear fractional differential equations with subdiffusion and superdiffusion,” Computers & Mathematics with Applications, vol. 62, no. 3, pp. 855–875, 2011.
- H. Zhang, F. Liu, and V. Anh, “Galerkin finite element approximation of symmetric space-fractional partial differential equations,” Applied Mathematics and Computation, vol. 217, no. 6, pp. 2534–2545, 2010.
- Y. Y. Zheng, C. P. Li, and Z. G. Zhao, “A fully discrete discontinuous Galerkin method for nonlinear fractional Fokker-Planck equation,” Mathematical Problems in Engineering, vol. 2010, Article ID 279038, 26 pages, 2010.
- 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.
- Q. Yang, F. Liu, and I. Turner, “Computationally efficient numerical methods for time and space fractional Fokker-Planck equations,” Physica Scripta, vol. 2009, Article ID 014026, 2009.
- C. P. Li, F. H. Zeng, and F. Liu, “Spectral approximations to the fractional integral and derivative,” Fractional Calculus and Applied Analysis, vol. 15, no. 3, pp. 383–406, 2012.
- X. Li and C. Xu, “A space-time spectral method for the time fractional diffusion equation,” SIAM Journal on Numerical Analysis, vol. 47, no. 3, pp. 2108–2131, 2009.
- H. F. Ding and C. P. Li, “Mixed spline function method for reaction-subdiusion equations,” Journal of Computational Physics, vol. 242, pp. 103–123, 2013.
- F. Mainardi, “Fractional relaxation-oscillation and fractional diffusion-wave phenomena,” Chaos, Solitons and Fractals, vol. 7, no. 9, pp. 1461–1477, 1996.
- J. Q. Murillo and S. B. Yuste, “An explicit difference method for solving fractional diffusion and diffusion-wave equations in the Caputo form,” Journal of Computational and Nonlinear Dynamics, vol. 6, no. 2, 6 pages, 2011.
- I. Podlubny, Fractional Differential Equations, Academic Press, San Diego, Calif, USA, 1999.
- 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.
- M. Li, T. Tang, and B. Fornberg, “A compact fourth-order finite difference scheme for the steady incompressible Navier-Stokes equations,” International Journal for Numerical Methods in Fluids, vol. 20, no. 10, pp. 1137–1151, 1995.
- L. Su, W. Wang, and Q. Xu, “Finite difference methods for fractional dispersion equations,” Applied Mathematics and Computation, vol. 216, no. 11, pp. 3329–3334, 2010.