Skip to main content Accessibility help

Lubich Second-Order Methods for Distributed-Order Time-Fractional Differential Equations with Smooth Solutions

  • Rui Du (a1), Zhao-peng Hao (a1) and Zhi-zhong Sun (a1)


This article is devoted to the study of some high-order difference schemes for the distributed-order time-fractional equations in both one and two space dimensions. Based on the composite Simpson formula and Lubich second-order operator, a difference scheme is constructed with convergence in the L1(L)-norm for the one-dimensional case, where τ,h and σ are the respective step sizes in time, space and distributed-order. Unconditional stability and convergence are proven. An ADI difference scheme is also derived for the two-dimensional case, and proven to be unconditionally stable and convergent in the L1(L)-norm, where h1 and h2 are the spatial step sizes. Some numerical examples are also given to demonstrate our theoretical results.


Corresponding author

*Corresponding author. Email (R. Du), (Z. Hao), (Z. Sun)


Hide All
[1]Podlubny,, I.Fractional Differential Equations, Academic Press, New York (1999).
[2]Oldham, K.B. and Spanier, J., The Fractional Calculus, Academic Press, New York (1974).
[3]Kilbas, A., Srivastava, H. and Trujillo, J., Theory and Applications of Fractional Differential Equations, Elsevier, Boston (2006).
[4]Mainardi, F., Fractional Calculus and Waves in Linear Viscoelasticity, Imperial College Press, London (2010).
[5]Huang, F. and Liu, F., The time fractional diffusion equation and the advection-dispersion equation, ANZIAM J. 46, 317330 (2005).
[6]Gorenflo, R., Luchko, Y.and Mainardi, F.Wright functions as scale-invariant solutions of the diffusion-wave equation, J. Comp. Appl. Math. 118, 175191 (2000).
[7]Povstenko, Y., Non-central-symmetric solution to time-fractional diffusion-wave equation in a sphere under Dirichlet boundary condition, Fract. Calc. Appl. Anal. 15, 253266 (2012).
[8]Ray, S.S. and Bera, R.K., Analytical solution of a fractional diffusion equation by Adomian decomposition method, Appl. Math. Comp. 174, 329336 (2006).
[9]Yuste, S. and Acedo, L., An explicit finite difference method and a new von Neumann-type stability analysis for fractional diffusion equations, SIAM J. Num. Anal. 42, 18621874 (2005).
[10]Langlands, T. and Henry, B., The accuracy and stability of an implicit solution method for the fractional diffusion equation, J. Comp. Phys. 205, 719736 (2005).
[11]Chen, C., Liu, F., Turner, I. and Anh, V., A Fourier method for the fractional diffusion equation describing sub-diffusion, J. Comp. Phys. 227, 886897 (2007).
[12]Zhuang, P., Liu, F., Anh, V. and Turner, I., New solution and analytical techniques of the implicit numerical method for the anomalous subdiffusion equation, SIAM J. Num. Anal. 46, 10791095 (2008).
[13]Cao, J. and Xu, C., A high order schema for the numerical solution of the fractional ordinary differential equations, J. Comp. Phys. 238, 154168 (2013).
[14]Li, C. and Ding, H., Higher order finite difference method for the reaction and anomalous-diffusion equation, Appl. Math. Model. 38, 38023821 (2014).
[15]Sun, Z.Z. and Wu, X., A fully discrete difference scheme for a diffusion-wave system, Appl. Num. Math. 56, 193209 (2006).
[16]Wang, Z. and Vong, S., Compact difference schemes for the modified anomalous fractional subdiffusion equation and the fractional diffusion-wave equation, J. Comp. Phys. 277, 115 (2014).
[17]Li, C. and Zeng, F., Finite difference methods for fractional differential equations, Int. J. Bifurcation Chaos 22(04), 1230014 (2012).
[18]Cui, M., Compact alternating direction implicit method for two-dimensional time fractional diffusion equation, J. Comp. Phys. 231, 26212633 (2012).
[19]Cui, M., Convergence analysis of high-order compact alternating direction implicit schemes for the two-dimensional time fractional diffusion equation, Num. Algor. 62, 383409 (2013).
[20]Zhang, Y.N. and Sun, Z.Z., Alternating direction implicit schemes for the two-dimensional fractional subdiffusion equation, J. Comp. Phys. 230, 87138728 (2011).
[21]Zhang, Y.N. and Sun, Z.Z., Error analysis of a compact ADI scheme for the 2D fractional subdiffusion equation, J. Sci. Comp. 59, 104128 (2014).
[22]Dimitrov, Y., Numerical approximations for fractional differential equations, J. Fract. Calc. Appl. 5, suppl. 3S, No. 22, 45 pp. (2014).
[23]Podlubny, I., Skovranek, T., Jara, B.M., Petras, J., Verbitsky, V. and Chen, Y.Q., Matrix approach to discrete fractional calculus III: Non-equidistant grids, variable step length and distributed orders, Phil. Trans. R. Soc. A 371, 20120153 (2012).
[24]Diethelm, K. and Ford, N.J., Numerical analysis for distributed order differential equations, J. Comp. Appl. Math. 225, 96104 (2009).
[25]Liu, F., Meerschaert, M.M., McGough, R.J., Zhuang, P. and Liu, Q., Numerical methods for solving the multi-term time-fractional wave-diffusion equation, Fract. Calc. Appl. Anal. 16, 925 (2013).
[26]Ye, H., Liu, F., Anh, V. and Turner, I., Numerical analysis for the time distributed-order and Riesz space fractional diffusions on bounded domains, IMA J. Appl. Math. 80, 825838 (2015).
[27]Ford, N.J., Morgado, M.L. and Rebelo, M., A numerical method for the distributed order time-fractional diffusion equation, ICFDA'14 Catania, ISBN 978-1-4799-2590-2, pp. 2325 (2014).
[28]Morgado, M.L. and Rebelo, M., Numerical approximation of distributed order reaction-diffusion equations, J. Comp. Appl. Math. 275, 216227 (2014).
[29]Katsikadelis, J.T., Numerical solution of distributed order fractional differential equations, J. Comp. Phys. 259, 1122 (2014).
[30]Lubich, C., Discretizedfractional calculus, SIAM J. Math. Anal. 17, 704719 (1986).
[31]Brunner, H. and Van der Houwen, P.J., The Numerical Solution of Volterra Equations, CWI Monographs, Vol. 3. North-Holland, Amsterdam, (1986).
[32]Zhou, H., Tian, W.Y. and Deng, W., Quasi-compact finite difference schemes for space fractional diffusion equations, J. Sci. Comp. 56, 4566 (2013).
[33]Sun, Z.Z., Numerical Methods of Partial Differential Equations, Science Press, Beijing (2012).
[34]Liao, H.L. and Sun, Z.Z., Maximum norm error bounds of ADI and compact ADI methods for solving parabolic equations, Num. Meth. Partial Differential Eq. 26, 3760 (2010).
[35]Henrici, P., Fast Fourier methods in computational complex analysis, SIAM Rev. 21, 481527 (1979).


Related content

Powered by UNSILO

Lubich Second-Order Methods for Distributed-Order Time-Fractional Differential Equations with Smooth Solutions

  • Rui Du (a1), Zhao-peng Hao (a1) and Zhi-zhong Sun (a1)


Full text views

Total number of HTML views: 0
Total number of PDF views: 0 *
Loading metrics...

Abstract views

Total abstract views: 0 *
Loading metrics...

* Views captured on Cambridge Core between <date>. This data will be updated every 24 hours.

Usage data cannot currently be displayed.