Abstract

We consider the numerical approximation of a time-fractional cable equation involving two Riemann–Liouville fractional derivatives. We investigate a semidiscrete scheme based on the lumped mass Galerkin finite element method (FEM), using piecewise linear functions. We establish optimal error estimates for smooth and middly smooth initial data, i.e., v∈Hq(Ω)∩H01(Ω), q=1,2. For nonsmooth initial data, i.e., v∈L2(Ω), the optimal L2(Ω)-norm error estimate requires an additional assumption on mesh, which is known to be satisfied for symmetric meshes. A quasi-optimal L∞(Ω)-norm error estimate is also obtained. Further, we analyze two fully discrete schemes using convolution quadrature in time based on the backward Euler and the second-order backward difference methods, and derive error estimates for smooth and nonsmooth data. Finally, we present several numerical examples to confirm our theoretical results.

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call

Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.