Fast IMEX Time Integration of Nonlinear Stiff Fractional Differential Equations

Abstract

Efficient long-time integration of nonlinear fractional differential equations is significantly challenging due to the integro-differential nature of the fractional operators. In addition, the inherent non-smoothness introduced by the inverse power-law kernels deteriorates the accuracy and efficiency of many existing numerical methods. We develop two efficient first- and second-order implicit-explicit (IMEX) methods for accurate time-integration of stiff/nonlinear fractional differential equations with fractional order $\alpha \in (0,1]$ and prove their convergence and linear stability properties. The developed methods are based on a linear multi-step fractional Adams-Moulton method (FAMM), followed by the extrapolation of the nonlinear force terms. In order to handle the singularities nearby the initial time, we employ Lubich-like corrections to the resulting fractional operators. The obtained linear stability regions of the developed IMEX methods are larger than existing IMEX methods in the literature. Furthermore, the size of the stability regions increase with the decrease of fractional order values, which is suitable for stiff problems. We also rewrite the resulting IMEX methods in the language of nonlinear Toeplitz systems, where we employ a fast inversion scheme to achieve a computational complexity of $\mathcal{O}(N\log N)$, where $N$ denotes the number of time-steps. Our computational results demonstrate that the developed schemes can achieve global first- and second-order accuracy for highly-oscillatory stiff/nonlinear problems with singularities.

References in corpus (9)

• Trapezoidal methods for fractional differential equations: theoretical and computational aspects
• Second-order numerical methods for multi-term fractional differential equations: Smooth and non-smooth solutions
• A Finite Difference Method with Non-uniform Timesteps for Fractional Diffusion Equations
• A Petrov-Galerkin Spectral Element Method for Fractional Elliptic Problems
• A Unified Spectral Method for FPDEs with Two-sided Derivatives; A Fast Solver
• A Unified Spectral Method for FPDEs with Two-sided Derivatives; Stability, and Error Analysis
• Petrov-Galerkin Method for Fully Distributed-Order Fractional Partial Differential Equations
• Vibration Analysis of Geometrically Nonlinear and Fractional Viscoelastic Cantilever Beams
• An Automated Singularity-Capturing Scheme for Fractional Differential Equations

Cited by in corpus (3)

• Equivalence between a time-fractional and an integer-order gradient flow: The memory effect reflected in the energy
• Solving time-fractional differential equation via rational approximation
• A Fractional Subgrid-scale Model for Turbulent Flows: Theoretical Formulation and a Priori Study