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.

Euler–Maruyama methods for Caputo tempered fractional stochastic differential equations

Well-posedness and an Euler-Maruyama method for multi-term caputo tempered fractional stochastic differential equations

An Euler–Maruyama method and its fast implementation for multiterm fractional stochastic differential equations

The fast Euler-Maruyama method for solving multiterm Caputo fractional stochastic delay integro-differential equations

An Euler–Maruyama method for Caputo–Hadamard fractional stochastic differential equations on exponential meshes and its fast approximation

Variable-Step Euler–Maruyama Approximations of Regime-Switching Jump Diffusion Processes

A fast Euler-Maruyama method for fractional stochastic differential equations

Semi-implicit Euler–Maruyama method for non-linear time-changed stochastic differential equations

Strong convergence of the truncated Euler–Maruyama method for stochastic functional differential equations

Strong convergence of Euler–Maruyama schemes for doubly perturbed McKean–Vlasov stochastic differential equations

Stochastic fractional integro-differential equations with weakly singular kernels: Well-posedness and Euler–Maruyama approximation

Truncated Euler-Maruyama method for classical and time-changed non-autonomous stochastic differential equations

The semi-implicit Euler-Maruyama method for nonlinear non-autonomous stochastic differential equations driven by a class of Lévy processes

Efficient numerical method for multi-term time-fractional diffusion equations with Caputo-Fabrizio derivatives

Two Mixed Finite Element Methods for Time-Fractional Diffusion Equations

Truncated Euler–Maruyama method for stochastic differential equations driven by fractional Brownian motion with super-linear drift coefficient

Multi-level Monte Carlo methods with the Truncated Euler-Maruyama Scheme for Stochastic Differential Equations

Convergence rates of full-implicit truncated Euler–Maruyama method for stochastic differential equations

Caputo-Hadamard fractional differential equations on time scales: Numerical scheme, asymptotic stability, and chaos

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