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.

Long Time Numerical Behaviors of Fractional Pantograph Equations.

Numerical Approach for Solving the Fractional Pantograph Delay Differential Equations

A New Numerical Approach for Solving the Fractional Nonlinear Multi-pantograph Delay Differential Equations

Results on Implicit Fractional Pantograph Equations with Mittag-Leffler Kernel and Nonlocal Condition

Computational and Numerical Analysis of the Caputo-Type Fractional Nonlinear Dynamical Systems via Novel Transform

A Solution to Parabolic System with the Fractional Laplacian

Finite difference methods for nonlinear fractional differential equation with ψ -Caputo derivative

Efficient spectral collocation method for nonlinear systems of fractional pantograph delay differential equations

A Fast Linearized Finite Difference Method for the Nonlinear Multi-Term Time-Fractional Wave Equation

A Corrected L1 Method for a Time-Fractional Subdiffusion Equation

Existence and Hyers–Ulam Stability for Random Impulsive Stochastic Pantograph Equations with the Caputo Fractional Derivative

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

An Efficient Numerical Scheme for a Time-Fractional Burgers Equation with Caputo-Prabhakar Derivative

Spectral collocation methods for fractional multipantograph delay differential equations*

A New High-Order Fractional Parallel Iterative Scheme for Solving Nonlinear Equations

Numerical Simulation and Solutions for the Fractional Chen System via Newly Proposed Methods

Separation method of semi-fixed variables together with dynamical system method for solving nonlinear time-fractional PDEs with higher-order terms

An Efficient Numerical Method for Solving a Class of Nonlinear Fractional Differential Equations and Error Estimates

Theoretical and numerical aspects for the longtime behavior of nonlinear delay time Caputo fractional reaction–diffusion equations

A novel numerical method to solve fractional ordinary differential equations with proportional Caputo derivatives

Numerical Analysis of a High-Order Scheme for Nonlinear Fractional Differential Equations with Uniform Accuracy