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.

NUMERICAL STABILITY ANALYSIS FOR EQUATION $ \mathit{x'}\left( \mathit{t} \right) = \mathit{ax}\left( \mathit{t} \right) + \mathit{bx}\left( {3\left[ {\left( {\mathit{t }}+ 1 \right)/3} \right]} \Right)$

Numerical Simulation for the System of Ordinary Differential Equations

Stability analysis of a mixed differential equation with piecewise continuous arguments

Stability of Runge–Kutta Methods in the Numerical Solution of Linear Impulsive Differential Equations

Explicit Unconditionally Stable Methods for the Heat Equation Via Potential Theory

Exact multiplicity and stability of solutions of second-order Neumann boundary value problem

Numerical Stability for Differential Equations with Memory

Stability analysis of multi-term fractional-differential equations with three fractional derivatives

Stochastic C-stability and B-consistency of explicit and implicit Euler-type schemes

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

Stability Analysis for Discretized Steady Convective-Diffusion Equation

Quantitative Stability and Numerical Analysis of Markovian Quadratic BSDEs with Reflection

Existence and Stability Analysis of Solution for Mathieu Fractional Differential Equations with Applications on Some Physical Phenomena

The Existence and Ulam Stability Analysis of a Multi-Term Implicit Fractional Differential Equation with Boundary Conditions

Euler Methods with Different Step-Sizes for Different Equations for Stochastic Differential Equations

Stability of Θ-Schemes in the Numerical Solution of a Partial Differential Equation with Piecewise Continuous Arguments

Stability Analysis for Nonlinear Neutral Stochastic Functional Differential Equations

Stability Equivalence Between the Neutral Delayed Stochastic Differential Equations and the Euler–Maruyama Numerical Scheme

Stability and qualitative dynamics of a class of three‐dimensional difference equation system of exponential type

Stability and Robustness of Time-discretization Schemes for the Allen-Cahn Equation via Bifurcation and Perturbation Analysis

Stability analysis of the numerical Method of characteristics applied to a class of energy-preserving systems. Part I: Periodic boundary conditions