Finite volume element method for two-dimensional fractional subdiffusion problems

TL;DR

This paper develops and analyzes a finite volume element method for two-dimensional fractional subdiffusion equations, providing optimal error estimates, superconvergence results, and a fully discrete scheme with proven convergence rates.

Contribution

It introduces a semi-discrete finite volume method for 2D fractional subdiffusion problems with rigorous error analysis and a fully discrete scheme combining FV and discontinuous Galerkin methods.

Findings

Optimal error estimates in $L^ abla(L^2)$ norm.

Superconvergence results with quasi-optimal $L^ abla(L^ abla)$ convergence.

Numerical experiments confirm theoretical convergence rates.

Abstract

In this paper, a semi-discrete spatial finite volume (FV) method is proposed and analyzed for approximating solutions of anomalous subdiffusion equations involving a temporal fractional derivative of order $\alpha \in (0,1)$ in a two-dimensional convex polygonal domain. Optimal error estimates in $L^\infty(L^2)$- norm is shown to hold. Superconvergence result is proved and as a consequence, it is established that quasi-optimal order of convergence in $L^{\infty}(L^{\infty})$ holds. We also consider a fully discrete scheme that employs FV method in space, and a piecewise linear discontinuous Galerkin method to discretize in temporal direction. It is, further, shown that convergence rate is of order $O(h^2+k^{1+\alpha}),$ where $h$ denotes the space discretizing parameter and $k$ represents the temporal discretizing parameter. Numerical experiments indicate optimal convergence rates in both time and space, and also illustrate that the imposed regularity assumptions are pessimistic.