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.

Analysis of the L1 scheme for fractional wave equations with nonsmooth data

An analysis of the L1 Scheme for the subdiffusion equation with nonsmooth data

The Uniform L1 Long-Time Behavior of Time Discretization for Time-Fractional Partial Differential Equations with Nonsmooth Data

Pointwise-in-time error analysis of the corrected L1 scheme for a time-fractional sine-Gordon equation

A New Analytical Technique of the L-type Difference Schemes for Time Fractional Mixed Sub-Diffusion and Diffusion-Wave Equations.

The Stability and Convergence Analysis of Finite Difference Methods for the Fractional Neutron Diffusion Equation

Analysis of a temporal discretization for a semilinear fractional diffusion equation

A Sharp $Α$-Robust $L1$ Scheme on Graded Meshes for Two-Dimensional Time Tempered Fractional Fokker-Planck Equation

Two Finite Difference Schemes for Multi-Dimensional Fractional Wave Equations with Weakly Singular Solutions

A Second-Order L 2-1 Σ Crank-Nicolson Difference Method for Two-Dimensional Time-Fractional Wave Equations with Variable Coefficients.

A Fractional Adams–Simpson-type Method for Nonlinear Fractional Ordinary Differential Equations with Non-Smooth Data

Energy stability and convergence of variable-step L1 scheme for the time fractional Swift-Hohenberg model

An Optimal Error Estimate for the Two-Dimensional Nonlinear Time Fractional Advection-Diffusion Equation with Smooth and Non-Smooth Solutions

<i>L2</i>-1<sub></sub> Finite Element Method for Time-Fractional Diffusion Problems with Discontinuous Coefficients

Unconditional analysis of the linearized second-order time-stepping scheme combined with a mixed element method for a nonlinear time fractional fourth-order wave equation

Mittag-Leffler Stability of Numerical Solutions to Linear Homogeneous Time Fractional Parabolic Equations

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

Numerical Blow-Up Analysis of the Explicit L1-scheme for Fractional Ordinary Differential Equations

Finite Element Approximations of a Class of Nonlinear Stochastic Wave Equations with Multiplicative Noise

Stability analysis of a second-order difference scheme for the time-fractional mixed sub-diffusion and diffusion-wave equation

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