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.

Preservation of Stability and Oscillation of Euler-Maclaurin Method for Differential Equation with Piecewise Constant Arguments of Alternately Advanced and Retarded Type

Stochastic Stability of Duffing Oscillator with Fractional Derivative Damping under Combined Harmonic and White Noise Parametric Excitations

Stochastic stability of the harmonically and randomly excited Duffing oscillator with damping modeled by a fractional derivative

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

Adapted Falkner-type Methods Solving Oscillatory Second-Order Differential Equations.

Mittag-Leffler stability of complete monotonicity-preserving schemes for time-dependent coefficients sub-diffusion equations

Polynomial stability of exact solution and a numerical method for stochastic differential equations with time-dependent delay

Oscillation of damped second order quasilinear wave equations with mixed arguments

Stability and Numerical Analysis of a Coupled System of Piecewise Atangana–Baleanu Fractional Differential Equations with Delays

Oscillation of modified Euler type half-linear differential equations via averaging technique

Convergence and stability of the balanced Euler method for stochastic pantograph differential equations with Markovian switching

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

Existence and k-Mittag Leffler Ulam stabilities of implicit pantograph differential equation via -Hilfer fractional derivative

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

Numerical Stability for Differential Equations with Memory

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

Functional Differential Equations with Distributed Deviating Arguments: Oscillatory Features of Solutions

Oscillation and Asymptotic Behavior of a Special Delay Third Order Nonlinear Neutral Functional Differential Equation

Stability for variable-coefficient wave equation with local frictional damping and memory type dynamic boundary conditions

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

Oscillation for Systems of a Class of Second Order Partial Eifferential Equations with Deviating Arguments