Error analysis of a FVEM for fractional order evolution equations with nonsmooth initial data

In this paper, a finite volume element (FVE) method is considered for spatial approximations of time fractional diffusion equations involving a Riemann-Liouville fractional derivative of order α ∈ (0, 1) in time. Improving upon earlier results [Karaa et al., IMA J. Numer. Anal. 37 (2017) 945-964], error estimates in L²(Ω)- and H¹(Ω)-norms for the semidiscrete problem with smooth and mildly smooth initial data, i.e., v ∈ H²(Ω) ∩ H¹₀ (Ω) and v ∈ H¹₀ (Ω) are established. For nonsmooth data, that is, v ∈ L²(Ω), the optimal L²(Ω)-error estimate is shown to hold only under an additional assumption on the triangulation, which is known to be satisfied for symmetric triangulations. Super-convergence result is also proved and as a consequence, a quasi-optimal error estimate is established in the L^∞(Ω)-norm. Further, two fully discrete schemes using convolution quadrature in time generated by the backward Euler and the second-order backward difference methods are analyzed, and error estimates are derived for both smooth and nonsmooth initial data. Based on a comparison of the standard Galerkin finite element solution with the FVE solution and exploiting tools for Laplace transforms with semigroup type properties of the FVE solution operator, our analysis is then extended in a unified manner to several time fractional order evolution problems. Finally, several numerical experiments are conducted to confirm our theoretical findings.