Research Article  Open Access
Finite Difference and SincCollocation Approximations to a Class of Fractional DiffusionWave Equations
Abstract
We propose an efficient numerical method for a class of fractional diffusionwave equations with the Caputo fractional derivative of order . This approach is based on the finite difference in time and the global sinc collocation in space. By utilizing the collocation technique and some properties of the sinc functions, the problem is reduced to the solution of a system of linear algebraic equations at each time step. Stability and convergence of the proposed method are rigorously analyzed. The numerical solution is of order accuracy in time and exponential rate of convergence in space. Numerical experiments demonstrate the validity of the obtained method and support the obtained theoretical results.
1. Introduction
Fractional diffusionwave equations have attracted considerable attention as generalization of the classical diffusion/wave equation by replacing the integerorder time derivative with a fractional derivative of order () [1]. It can be derived from the physical case of a general random timedependent velocity function equipped with an algebraic timecorrelation function (but with a Gaussian space correlation) [2]. Owing to the anomalous superdiffusion of particles, many electromagnetic, acoustic, mechanical, and biological responses can be modelled accurately by fractional diffusionwave equations, for example, the propagation of stress waves in viscoelastic solids [3, 4]. There exist some analytical methods to find exact solutions of fractional diffusionwave equations [5–9]. However, it is usually difficult or even impossible to achieve the exact solutions in many cases, and ones have to resort to numerical methods.
Compared with the considerable literature on numerical solutions of fractional diffusion equations [10–15], only a few works have been carried out in the numerical methods of fractional diffusionwave equations. The homotopy analysis method [16] and the Adomian decomposition method [17] were used to construct analytical approximate solutions of fractional diffusionwave equations, respectively. Sun and Wu [18] presented a full discrete difference scheme for the initialboundaryvalue problem of a diffusionwave equation by introducing two new variables and obtained the stability and convergence properties by using energy method. In [19], under the weak smoothness conditions, two finite difference schemes with firstorder accuracy in temporal direction and secondorder accuracy in spatial direction were proposed to solve the same problem. Murillo and Yuste [20] developed an explicit difference method where the L2 discretization formula is used for solving fractional diffusionwave equations. Due to the nonlocal nature of fractional derivatives, all previous solutions have to be saved to compute the solution at the current time level. It makes the storage very expensive when loworder methods are employed for spatial discretization.
In this paper, we propose a numerical method mixing finite difference with sinc collocation to solve a class of fractional diffusionwave equations. The sinc method is widely used to solve integral, integrodifferential, ordinary differential, and partial differential equations [21–27]. As far as we know, there are very limited papers on solving fractional differential equation by the sinc method [28, 29]. The present method has great advantage in storage requirement because the sinccollocation method needs fewer grid points to produce highly accurate solution when it is compared with a loworder method.
The remainder of this paper is organized as follows. In Section 2, we introduce some necessary definitions and relevant results for developing this method. Section 3 is devoted to the finite difference and sinccollocation discretization for the fractional diffusionwave equations. As a result, a system of linear algebraic equations is formed, and the solutions of the considered problems are obtained. Section 4 is concerned with the stability and convergence analysis. In Section 5, some numerical experiments are given to demonstrate the effectiveness of the proposed method and the obtained theoretical results. A brief conclusion ends this paper in the final section.
2. Notations and Some Preliminary Results
In this section, we introduce some basic definitions and relevant results of the fractional calculus [30, 31] and sinc functions.
Definition 1 (see [32]). Let . The operator defined on by for is called the RiemannLiouville fractional integral operator of order . For , one sets , that is, the identity operator.
Definition 2 (see [32]). Let and . The Caputo fractional differential operator for is defined as
The sinc functions and its properties are discussed thoroughly in [21, 23]. For any , the translated sinc functions with equidistant space nodes are given as where the sinc functions are defined on the whole real line by If is defined on , then for any the series is called the Whittaker cardinal expansion of whenever this series converges. can be approximated by truncating (5). To construct our needed approximations on the interval , we choose which maps a finite interval onto . is a conformal map which maps the eyeshaped domain in the complex plane onto the infinite strip of the complex plane
The basis functions on are taken to be the composite translated sinc functions Thus one may define the inverse image of the equidistant space node as The class of functions such that the known exponential convergence rate exists for the sinc interpolation is denoted by and defined in the following.
Definition 3 (see [21, 23]). Let be the class of functions which are analytic in and satisfy where and on the boundary of (denoted ).
The following theorem gives the error which results from differentiating the truncated cardinal series.
Theorem 4 (see [21, 22]). Let and let for ; here is a weight function and is a constant depending only on , , and . Assume that there exists a constant such that here is a positive constant. Then taking , one has for all ; here is a constant depending only on , , , , , and .
The above expressions show that the sinc interpolation on converges exponentially. We also require the following derivatives of the composite translated sinc functions evaluated at the nodes:
3. The Finite Difference and SincCollocation Discretization
We consider the following fractional diffusionwave equations with a nonhomogeneous field [33]: with the initial conditions and the boundary conditions where and are space and time variables, and are constants, and denotes the field variable. Here the timefractional derivative is defined as the Caputo fractional derivative, and . Many authors refer to the fractional equation (19) as the fractional diffusionwave equation when , which is expected to interpolate the diffusion equation and the wave equation [7, 18].
3.1. Temporal Discretization by a Finite Difference Scheme
First, we derive a finite difference scheme for temporal discretization of this equation. Let , , where is the time stepsize. In order to discretize the timefractional derivative by using a finite difference approximation [18, 34], we introduce the following lemmas.
Lemma 5 (see [18]). Suppose . Then where and
It is direct to check that
Let and define the temporal grid functions By the Taylor expansion, it follows that where , are constants. Based on Lemma 5, one has Hence, where are constants. Let . In view of (27), one has Then substituting (31) into (28) and observing , one obtains where
It follows from (32) that one can construct easily the following semidiscrete finite difference scheme for (19): which is equivalent to where According to (33), this scheme is of order accuracy. A rigorous analysis of the convergence rate will be provided later. The above scheme can be rewritten into Specially, for the case , the scheme simply reads So (37) and (38) together with the initial condition and the boundary conditions form a complete semidiscrete problem.
3.2. Space Discretization by the SincCollocation Method
Next we consider space discretization to (37) by the sinccollocation method. We select the collocation points by (10). The space discretization proceeds by approximating the solution based on the composite translated sinc functions (9) The unknown parameters will be determined by the collocation method. It is noted that the approximation in (41) satisfies the boundary conditions in (40) since , are zero when tends to and .
Now substituting (41) into (37) and collocating in points , we obtain where . Based on (17) and (18), we let Hence, (42) is reduced immediately by (16) and (43) to where . To obtain a matrix representation of the above equation, we let where the matrix is the identity matrix and . Therefore, at each time step, we get the following system of linear equations with unknown parameters , and this system can be expressed in a matrix form where For the initial condition (39), we have Consequently (46) can be solved easily for the unknown coefficients . Hence the approximation solution given in (41) can be obtained.
4. Stability and Convergence Analysis of the Derived Method
To analyse the stability and convergence of the derived method, the inner product is defined by with the corresponding norm , which will be used thereafter. We first give the following lemmas.
Lemma 6 (see [18]). For any and , here , it can be verified that where is defined in (23).
Lemma 7. Suppose that , , satisfy Then one has the estimate
Proof. Based on (50), we have Applying Lemma 6 to the term in left side of (54), we obtain On the other side, based on the fact that , a straightforward calculation of the right terms in (54) gives Substituting (55) and (56) into (54) yields inequality (53).
Suppose that , , satisfy the conditions of Theorem 4. Substituting the approximation solution given in (41) into (34) leads to Subtracting (57) from (34), we obtain Based on Theorem 4, we have where are constants.
Theorem 8. Let , , be the approximation solution given in (41) and let satisfy the conditions of Theorem 4. Then one has where is constant.
Proof. Considering (57) and noting that and , , we can easily obtain this result from Lemma 7 and the Poincaré inequality.
Theorem 9. Let the problem (19)–(21) have the exact solution and let be the approximation solution given in (41), where . If satisfy the conditions of Theorem 4, then for , one has where .
Proof. Subtracting (57) from (32), we obtain the error equation where . Noting that and , , and applying Lemma 7, we have Based on (33) and (59), we obtain Applying the Poincaré inequality yields the needed result.
Remark 10. In Theorem 9, the error estimate formula contains (the second term in the righthand side), where the error contribution from the spatial approximation is affected by the inverse of the time step. It is worthwhile noting that similar results are also found for the classical diffusion/wave equation. However for large , can be much smaller than ; therefore, this affection generally would not reduce the global accuracy.
5. Numerical Examples
To validate the effectiveness of the proposed method for the problem (19)–(21), we consider the example given in [18]. Consider the following: The exact solution of the above problem is [33] where which is oneparameter MittagLeffler function.
For solving the above problem (with and , resp.) by using the method described in Section 3, we choose and , and this leads to . We will report the efficiency and accuracy of the given method based on the errors and errors. Figure 1 gives the 3D diagrams of the numerical solutions on the whole computational domain with , . In Figure 2, we plot the curves of the numerical solutions and exact solutions for several fixed time instants. A good agreement of the numerical solution with the exact one is achieved. Furthermore, Figure 3 shows the absolute error obtained by the present method with , .
(a)
(b)
(a)
(b)
(a)
(b)
To check the convergence behavior of the numerical solution with respect to the time step and the parameter about the number of collocation points, we represent the errors in two discrete norms: and . All the numerical results reported in Figures 4 and 5, and Tables 1 and 2 are evaluated at .


(a)
(b)
(a)
(b)
Firstly, the computational investigation is concerned with the temporal convergence rate. We choose the parameter and it is a value large enough such that the errors coming from the spatial approximation are negligible [35]. In Tables 1 () and 2 (), we list the temporal errors when decreases from 1/10 to 1/320 and the convergence order which is very close to the expectation order . We also plot the errors in the and norms as a function of the time stepsize for several in Figure 4, where a logarithmic scale has been used for both axis and error axis. From Figure 4, it is clear that, for and 1.7, the slopes of the error curves in these loglog plots approach 1.7 and 1.3, respectively. So the proposed method yields a temporal approximation order which is close to as forecasted by the theoretical estimate.
Now we check the spatial accuracy with respect to the parameter about the number of collocation points. For the similar reason mentioned above, we fix the time step sufficiently small to avoid contamination of the temporal error. In Figure 5, we present the spatial errors as a function of with , where a logarithmic scale is now used for the spatial errors axis. It is clearly observed that the spatial errors appear in an exponential decay. In this semilog representations, we observe that the error variations are essentially linear versus [36].
6. Conclusion
In this paper, we have developed and analyzed an efficient numerical method for a class of the initialboundaryvalue problems of fractional diffusionwave equations. Based on a finite difference scheme in time and a global sinccollocation method in space, the problem is reduced to the solution of a system of linear algebraic equations at each time step. The proposed method leads to order accuracy in time and exponential convergence in space. The theoretical results are perfectly confirmed by the numerical experiments.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The present research is supported by the National Natural Science Foundation of China (Grant nos. 11271311 and 11371016), the Chinese Program for Changjiang Scholars and Innovative Research Team in University (PCSIRT) (Grant no. IRT1179), and the Hunan Province Innovation Foundation for Postgraduate (Grant no. CX2013B252). The authors are grateful to the anonymous reviewers for their careful comments and valuable suggestions on this paper.
References
 F. Mainardi, “Fractional diffusive waves in viscoelastic solids,” in Nonlinear Waves in Solids, J. L. Wegner and F. R. Norwood, Eds., pp. 93–97, ASME/AMR, Fairfield, NJ, USA, 1995. View at: Google Scholar
 R. Balescu, “VLangevin equations, continuous time random walks and fractional diffusion,” Chaos, Solitons & Fractals, vol. 34, no. 1, pp. 62–80, 2007. View at: Publisher Site  Google Scholar
 F. Mainardi, “Fractional relaxationoscillation and fractional diffusionwave phenomena,” Chaos, Solitons & Fractals, vol. 7, no. 9, pp. 1461–1477, 1996. View at: Publisher Site  Google Scholar
 A. A. Kilbas, H. M. Srivastava, and J. J. Trujillo, Theory and Applications of Fractional Differential Equations, Elsevier, Amsterdam, The Netherland, 2006.
 F. Mainardi, “The fundamental solutions for the fractional diffusionwave equation,” Applied Mathematics Letters, vol. 9, no. 6, pp. 23–28, 1996. View at: Publisher Site  Google Scholar  MathSciNet
 S. Momani, Z. Odibat, and V. S. Erturk, “Generalized differential transform method for solving a space and timefractional diffusionwave equation,” Physics Letters A, vol. 370, no. 56, pp. 379–387, 2007. View at: Publisher Site  Google Scholar  MathSciNet
 O. P. Agrawal, “Solution for a fractional diffusionwave equation defined in a bounded domain,” Nonlinear Dynamics, vol. 29, no. 1–4, pp. 145–155, 2002. View at: Publisher Site  Google Scholar  MathSciNet
 Y. Luchko, “Fractional wave equation and damped waves,” Journal of Mathematical Physics, vol. 54, no. 3, Article ID 031505, 16 pages, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 Y. Povstenko, “Axisymmetric solutions to fractional diffusionwave equation in a cylinder under Robin boundary condition,” The European Physical Journal Special Topics, vol. 222, no. 8, pp. 1767–1777, 2013. View at: Publisher Site  Google Scholar
 I. Podlubny, A. Chechkin, T. Skovranek, Y. Chen, and B. M. Vinagre Jara, “Matrix approach to discrete fractional calculus II: partial fractional differential equations,” Journal of Computational Physics, vol. 228, no. 8, pp. 3137–3153, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 S. B. Yuste, “Weighted average finite difference methods for fractional diffusion equations,” Journal of Computational Physics, vol. 216, no. 1, pp. 264–274, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 P. Zhuang, F. Liu, V. Anh, and I. Turner, “New solution and analytical techniques of the implicit numerical method for the anomalous subdiffusion equation,” SIAM Journal on Numerical Analysis, vol. 46, no. 2, pp. 1079–1095, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. M. Khader, “On the numerical solutions for the fractional diffusion equation,” Communications in Nonlinear Science and Numerical Simulation, vol. 16, no. 6, pp. 2535–2542, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 C. Piret and E. Hanert, “A radial basis functions method for fractional diffusion equations,” Journal of Computational Physics, vol. 238, pp. 71–81, 2013. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 A. Mohebbi, M. Abbaszadeh, and M. Dehghan, “A highorder and unconditionally stable scheme for the modified anomalous fractional subdiffusion equation with a nonlinear source term,” Journal of Computational Physics, vol. 240, pp. 36–48, 2013. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 H. Jafari, A. Golbabai, S. Seifi, and K. Sayevand, “Homotopy analysis method for solving multiterm linear and nonlinear diffusionwave equations of fractional order,” Computers and Mathematics with Applications, vol. 59, no. 3, pp. 1337–1344, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 Z. M. Odibat and S. Momani, “Approximate solutions for boundary value problems of timefractional wave equation,” Applied Mathematics and Computation, vol. 181, no. 1, pp. 767–774, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Z. Sun and X. Wu, “A fully discrete difference scheme for a diffusionwave system,” Applied Numerical Mathematics, vol. 56, no. 2, pp. 193–209, 2006. View at: Publisher Site  Google Scholar  MathSciNet
 J. Huang, Y. Tang, L. Vázquez, and J. Yang, “Two finite difference schemes for time fractional diffusionwave equation,” Numerical Algorithms, vol. 64, no. 4, pp. 707–720, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 J. Q. Murillo and S. B. Yuste, “An explicit difference method for solving fractional diffusion and diffusionwave equations in the Caputo form,” Journal of Computational and Nonlinear Dynamics, vol. 6, no. 2, Article ID 021014, 6 pages, 2011. View at: Publisher Site  Google Scholar
 F. Stenger, Numerical Methods Based on Sinc and Analytic Functions, Springer, New York, NY, USA, 1993.
 F. Stenger, Handbook of Sinc Numerical Methods, CRC Press, New York, NY, USA, 2011.
 J. Lund and K. L. Bowers, Sinc Methods f or Quadrature and Differential Equations, SIAM, Philadelphia, Pa, USA, 1992.
 K. Parand, M. Dehghan, and A. Pirkhedri, “The sinccollocation method for solving the ThomasFermi equation,” Journal of Computational and Applied Mathematics, vol. 237, no. 1, pp. 244–252, 2013. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 K. AlKhaled, “Numerical study of Fisher's reactiondiffusion equation by the Sinc collocation method,” Journal of Computational and Applied Mathematics, vol. 137, no. 2, pp. 245–255, 2001. View at: Publisher Site  Google Scholar  MathSciNet
 A. Secer, “Numerical solution and simulation of secondorder parabolic PDEs with SincGalerkin method using Maple,” Abstract and Applied Analysis, vol. 2013, Article ID 686483, 10 pages, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 M. Dehghan and F. EmamiNaeini, “The Sinccollocation and SincGalerkin methods for solving the twodimensional Schrödinger equation with nonhomogeneous boundary conditions,” Applied Mathematical Modelling, vol. 37, no. 22, pp. 9379–9397, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 A. Saadatmandi, M. Dehghan, and M. Azizi, “The SincLegendre collocation method for a class of fractional convectiondiffusion equations with variable coefficients,” Communications in Nonlinear Science and Numerical Simulation, vol. 17, no. 11, pp. 4125–4136, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 G. Baumann and F. Stenger, “Fractional calculus and Sinc methods,” Fractional Calculus and Applied Analysis, vol. 14, no. 4, pp. 568–622, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 I. Podlubny, Fractional Differential Equations, Academic Press, New York, NY, USA, 1999. View at: MathSciNet
 D. Baleanu, K. Diethelm, E. S, and J. J. Trujillo, Fractional Calculus: Models and Numerical Methods, World Scientific, Singapore, 2012.
 M. Weilbeer, Efficient Numerical Methods for Fractional Differential Equations and their Analytical Background, Papierflieger, ClausthalZellerfeld, Germany, 2006.
 O. P. Agrawal, “Response of a diffusionwave system subjected to deterministic and stochastic fields,” Zeitschrift für Angewandte Mathematik und Mechanik, vol. 83, no. 4, pp. 265–274, 2003. View at: Publisher Site  Google Scholar  MathSciNet
 K. Diethelm, The Analysis of Fractional Differential Equations, Springer, Berlin, Germany, 2010. View at: MathSciNet
 Y. Lin and C. Xu, “Finite difference/spectral approximations for the timefractional diffusion equation,” Journal of Computational Physics, vol. 225, no. 2, pp. 1533–1552, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 X. Li and C. Xu, “A spacetime spectral method for the time fractional diffusion equation,” SIAM Journal on Numerical Analysis, vol. 47, no. 3, pp. 2108–2131, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
Copyright
Copyright © 2014 Zhi Mao et al. 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.