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

TL;DR

This paper introduces novel fractional $ heta$-methods combined with finite element techniques to efficiently solve the fractional Cable model, achieving optimal stability and convergence for smooth solutions, validated by numerical experiments.

Contribution

The authors develop and analyze two new fractional $ heta$-methods integrated with finite element discretization for the fractional Cable model, ensuring stability and optimal convergence.

Findings

Proved positivity properties of the methods' coefficients.

Established stability and optimal convergence rate of $O( au^2+h^{r+1})$.

Numerical experiments confirm theoretical results and address initial singularity issues.

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.