Numerical Solution of Multiterm Fractional Differential Equations Using the Matrix Mittag–Leffler Functions

Multiterm fractional differential equations (MTFDEs) nowadays represent a widely used tool to model many important processes, particularly for multirate systems. Their numerical solution is then a compelling subject that deserves great attention, not least because of the difficulties to apply general purpose methods for fractional differential equations (FDEs) to this case. In this paper, we first transform the MTFDEs into equivalent systems of FDEs, as done by Diethelm and Ford; in this way, the solution can be expressed in terms of Mittag–Leffler (ML) functions evaluated at matrix arguments. We then propose to compute it by resorting to the matrix approach proposed by Garrappa and Popolizio. Several numerical tests are presented that clearly show that this matrix approach is very accurate and fast, also in comparison with other numerical methods.

[1]  Igor Moret,et al.  On the time-fractional Schrödinger equation: Theoretical analysis and numerical solution by matrix Mittag-Leffler functions , 2017, Comput. Math. Appl..

[2]  Igor Moret,et al.  Solving the time-fractional Schrödinger equation by Krylov projection methods , 2015, J. Comput. Phys..

[3]  K. Diethelm,et al.  Fractional Calculus: Models and Numerical Methods , 2012 .

[4]  R. Gorenflo,et al.  AN OPERATIONAL METHOD FOR SOLVING FRACTIONAL DIFFERENTIAL EQUATIONS WITH THE CAPUTO DERIVATIVES , 1999 .

[5]  R. Gorenflo,et al.  Fractional Calculus: Integral and Differential Equations of Fractional Order , 2008, 0805.3823.

[6]  Yury F. Luchko.  Initial-boundary-value problems for the generalized multi-term time-fractional diffusion equation , 2011 .

[7]  F. Mainardi.  Fractional Relaxation-Oscillation and Fractional Diffusion-Wave Phenomena , 1996 .

[8]  Neville J. Ford,et al.  Systems-based decomposition schemes for the approximate solution of multi-term fractional differential equations , 2009 .

[9]  N. Higham.  Functions Of Matrices , 2008 .

[10]  R. Bagley,et al.  On the Appearance of the Fractional Derivative in the Behavior of Real Materials , 1984 .

[11]  N. Ford,et al.  The numerical solution of linear multi-term fractional differential equations: systems of equations , 2002 .

[12]  T. Kaczorek,et al.  Fractional Differential Equations , 2015 .

[13]  Roberto Garrappa,et al.  Trapezoidal methods for fractional differential equations: Theoretical and computational aspects , 2015, Math. Comput. Simul..

[14]  Francesco Mainardi,et al.  Fractional Calculus: Some Basic Problems in Continuum and Statistical Mechanics , 2012, 1201.0863.

[15]  Igor Moret.  A Note on Krylov Methods for Fractional Evolution Problems , 2013 .

[16]  N. Ford,et al.  Pitfalls in fast numerical solvers for fractional differential equations , 2006 .

[17]  I. Podlubny.  Fractional differential equations , 1998 .

[18]  K. Diethelm,et al.  The Fracpece Subroutine for the Numerical Solution of Differential Equations of Fractional Order , 2002 .

[19]  R. F. Cameron,et al.  Product integration methods for second-kind Abel integral equations , 1984 .

[20]  Roberto Garrappa,et al.  Explicit methods for fractional differential equations and their stability properties , 2009 .

[21]  Roberto Garrappa,et al.  Numerical Solution of Fractional Differential Equations: A Survey and a Software Tutorial , 2018 .

[22]  Arak M. Mathai,et al.  Mittag-Leffler Functions and Their Applications , 2009, J. Appl. Math..

[23]  Shahrokh Esmaeili.  THE NUMERICAL SOLUTION OF THE BAGLEY-TORVIK EQUATION BY EXPONENTIAL INTEGRATORS , 2017 .

[24]  Roberto Garrappa,et al.  Stability-Preserving High-Order Methods for Multiterm fractional differential equations , 2012, Int. J. Bifurc. Chaos.

[25]  F. Mainardi,et al.  Models of dielectric relaxation based on completely monotone functions , 2016, 1611.04028.

[26]  Roberto Garrappa,et al.  Numerical Evaluation of Two and Three Parameter Mittag-Leffler Functions , 2015, SIAM J. Numer. Anal..

[27]  Roberto Garrappa,et al.  Evaluation of generalized Mittag–Leffler functions on the real line , 2013, Adv. Comput. Math..

[28]  Roberto Garrappa,et al.  On accurate product integration rules for linear fractional differential equations , 2011, J. Comput. Appl. Math..

[29]  C. Lubich,et al.  Runge-Kutta theory for Volterra and Abel integral equations of the second kind , 1983 .

[30]  C. Lubich.  Discretized fractional calculus , 1986 .

[31]  Mohammed Al-Refai,et al.  Maximum principle for the multi-term time-fractional diffusion equations with the Riemann-Liouville fractional derivatives , 2015, Appl. Math. Comput..

[32]  C. Lubich,et al.  Fractional linear multistep methods for Abel-Volterra integral equations of the second kind , 1985 .

[33]  N. Ford,et al.  Numerical Solution of the Bagley-Torvik Equation , 2002, BIT Numerical Mathematics.

[34]  Roberto Garrappa,et al.  On some explicit Adams multistep methods for fractional differential equations , 2009 .

[35]  Alan D. Freed,et al.  Detailed Error Analysis for a Fractional Adams Method , 2004, Numerical Algorithms.

[36]  N. Ford,et al.  Analysis of Fractional Differential Equations , 2002 .

[37]  Igor Moret,et al.  On the Convergence of Krylov Subspace Methods for Matrix Mittag-Leffler Functions , 2011, SIAM J. Numer. Anal..

[38]  Roberto Garrappa,et al.  Effect of perturbation in the numerical solution of fractional differential equations , 2017 .

[39]  Roberto Garrappa,et al.  Computing the Matrix Mittag-Leffler Function with Applications to Fractional Calculus , 2018, Journal of Scientific Computing.

[40]  Kai Diethelm,et al.  Multi-order fractional differential equations and their numerical solution , 2004, Appl. Math. Comput..

[41]  M. T. Cicero.  FRACTIONAL CALCULUS AND WAVES IN LINEAR VISCOELASTICITY , 2012 .