Abstract:The asymptotic stable region and long-time decay rate of solutions to linear homogeneous Caputo time fractional ordinary differential equations (F-ODEs) are known to be completely determined by the eigenvalues of the coefficient matrix. Very different from the exponential decay of solutions to classical ODEs, solutions of F-ODEs decay only polynomially, leading to the so-called Mittag-Leffler stability, which was already extended to semi-linear F-ODEs with small perturbations. This work is mainly devoted to the qualitative analysis of the long-time behavior of numerical solutions. By applying the singularity analysis of generating functions developed by Flajolet and Odlyzko (SIAM J. Disc. Math. 3 (1990), 216-240), we are able to prove that both $\mathcal{L}$1 scheme and strong $A$-stable fractional linear multistep methods (F-LMMs) can preserve the numerical Mittag-Leffler stability for linear homogeneous F-ODEs exactly as in the continuous case. Through an improved estimate of the discrete fractional resolvent operator, we show that strong $A$-stable F-LMMs are also Mittag-Leffler stable for semi-linear F-ODEs under small perturbations. For the numerical schemes based on $\alpha$-difference approximation to Caputo derivative, we establish the Mittag-Leffler stability for semi-linear problems by making use of properties of the Poisson transformation and the decay rate of the continuous fractional resolvent operator. Numerical experiments are presented for several typical time fractional evolutional equations, including time fractional sub-diffusion equations, fractional linear system and semi-linear F-ODEs. All the numerical results exhibit the typical long-time polynomial decay rate, which is fully consistent with our theoretical predictions.

On the Sheffer-type polynomials related to the Mittag-Leffler functions: applications to fractional evolution equations

THE FUNDAMENTAL SOLUTIONS FOR MULTI-TERM MODIFIED POWER LAW WAVE EQUATIONS IN A FINITE DOMAIN Vol. 1(1) Jan. 2013, No. 7, Pp. 55-66 H. JIANG, F. LIU, M. M. MEERSCHAERT, R. J. MCGOUGH, Q. LIU

Extended Mittag-Leffler function and associated fractional calculus operators

Theory on linear L-fractional differential equations and a new Mittag-Leffler-type function

Mittag-leffler-type function of arbitrary order and their application in the fractional kinetic equation

On some new properties of fractional derivatives with Mittag-Leffler kernel

Fractional Helmholtz and fractional wave equations with Riesz-Feller and generalized Riemann-Liouville fractional derivatives

Application of Fractional Derivative Operators to Anomalous Diffusion and Propagation Problems

Theory on Linear L-Fractional Differential Equations and a New Mittag–Leffler-Type Function

On the Equivalence between Differential and Integral Forms of Caputo-Type Fractional Problems on Hölder Spaces

Applications of the Mittag–Leffler function in solvability and stability of a class of fractional Langevin equations with two fractional orders

On the fractional abstract Schrodinger-type evolution equations on the Hilbert space and its applications to the fractional dispersive equations

Analysis of fractional Cauchy problems with some probabilistic applications

Mittag--Leffler stability of numerical solutions to time fractional ODEs

Quasi-Linear (Stochastic) Partial Differential Equations with Time-Fractional Derivatives

Completely monotone multinomial mittag-leffler type functions and diffusion equations with multiple time-derivatives

Spectral Collocation Method for Fractional Differential/Integral Equations with Generalized Fractional Operator

Laguerre-Type Bernoulli and Euler Numbers and Related Fractional Polynomials

Solutions of fractional reaction-diffusion equations in terms of Mittag-Leffler functions

Hermite, Higher order Hermite, Laguerre type polynomials and Burgers like equations

Fractional Calculus Operator Emerging from the 2D Biorthogonal Hermite Konhauser Polynomials