Finite element methods based on two families of second-order numerical formulas for the fractional Cable model with smooth solutions

Abstract

We apply two families of novel fractional $\theta$-methods, the FBT-$\theta$ and FBN-$\theta$ methods developed by the authors in previous work, to the fractional Cable model, in which the time direction is approximated by the fractional $\theta$-methods, and the space direction is approximated by the finite element method. Some positivity properties of the coefficients for both of these methods are derived, which are crucial for the proof of the stability estimates. We analyse the stability of the scheme and derive an optimal convergence result with $O({\tau }^2+h^{r+1})$ for smooth solutions, where $\tau$ is the time mesh size and $h$ is the spatial mesh size. Some numerical experiments with smooth and nonsmooth solutions are conducted to confirm our theoretical analysis. To overcome the singularity at initial value, the starting part is added to restore the second-order convergence rate in time.

arXiv comment24 pages, 9 figures

References in corpus (5)

• A Two-Grid Finite Element Approximation for A Nonlinear Time-Fractional Cable Equation
• Fast algorithm based on TT-M FE system for space fractional Allen-Cahn equations with smooth and non-smooth solutions
• High-order numerical methods for the Riesz space fractional advection-dispersion equations
• A new second-order midpoint approximation formula for Riemann-Liouville derivative: algorithm and its application
• Two families of novel second-order fractional numerical formulas and their applications to fractional differential equations