Second-order numerical methods for multi-term fractional differential equations: Smooth and non-smooth solutions

Abstract

Starting with the asymptotic expansion of the error equation of the shifted Grünwald--Letnikov formula, we derive a new modified weighted shifted Grünwald--Letnikov (WSGL) formula by introducing appropriate correction terms. We then apply one special case of the modified WSGL formula to solve multi-term fractional ordinary and partial differential equations, and we prove the linear stability and second-order convergence for both smooth and non-smooth solutions. We show theoretically and numerically that numerical solutions up to certain accuracy can be obtained with only a few correction terms. Moreover, the correction terms can be tuned according to the fractional derivative orders without explicitly knowing the analytical solutions. Numerical simulations verify the theoretical results and demonstrate that the new formula leads to better performance compared to other known numerical approximations with similar resolution.

31 pages, 24 figures

References in corpus (3)

• The Galerkin Finite Element Method for A Multi-term Time-Fractional Diffusion equation
• Initial-boundary value problems for multi-term time-fractional diffusion equations with positive constant coefficients
• High-order numerical methods for the Riesz space fractional advection-dispersion equations

Cited by in corpus (7)

• A preconditioning technique for an all-at-once system from Volterra subdiffusion equations with graded time steps
• An unstructured mesh control volume method for two-dimensional space fractional diffusion equations with variable coefficients on convex domains
• A spectrally accurate step-by-step method for the numerical solution of fractional differential equations
• A discrete least squares collocation method for two-dimensional nonlinear time-dependent partial differential equations
• Fast IMEX Time Integration of Nonlinear Stiff Fractional Differential Equations
• Analysis and implementation of collocation methods for fractional differential equations
• Numerical analysis of linear and nonlinear time-fractional subdiffusion equations