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

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})$, 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.