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.

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

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

Stability of Time-Marching MPS–MFS for Wave Equations

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

A Symplectic Algorithm For The Stability Analysis Of Nonlinear Parametric Excited Systems

POLYNOMIAL STABILITY FOR WAVE EQUATIONS WITH WENTZELL BOUNDARY CONDITIONS

An Analysis of the Numerical Stability of the Immersed Boundary Method

Stability Analysis for Discretized Steady Convective-Diffusion Equation

Von Neumann Stability Analysis for Multi-level Multi-step Methods

The continuation and stability analysis methods for quasi-periodic solutions of nonlinear systems

Stability and Convergence of Mixed Discontinuous Finite Element Methods for Second-Order Differential Problems

Linearized Stability and Instability of Nonconstant Periodic Solutions of Lagrangian Equations

Stability and Numerical Simulation of a Nonlinear Hadamard Fractional Coupling Laplacian System with Symmetric Periodic Boundary Conditions

Stability analysis of an extended quadrature method of moments for kinetic equations

Mittag-Leffler Stabilization of an Unstable Time Fractional Hyperbolic PDE

Characteristics Analysis of Two-dimensional Parabolized Stability Equations

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

A numerical method with properties of consistency in the energy domain for a class of dissipative nonlinear wave equations with applications to a Dirichlet boundary-value problem

Characteristic Stabilized Nonconforming Finite Element Metho d for the Unsteady Incompressible Navier-Stokes Equations

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

Equilibrium Stability Analysis of Hyperbolic Shallow Water Moment Equations