Superconvergence of a discontinuous Galerkin method for fractional diffusion and wave equations

Abstract

We consider an initial-boundary value problem for ${\partial }_{t}u-{\partial }_t^{-\alpha }{\nabla }^{2}u=f(t)$, that is, for a fractional diffusion ($-1<\alpha <0$) or wave ($0<\alpha <1$) equation. A numerical solution is found by applying a piecewise-linear, discontinuous Galerkin method in time combined with a piecewise-linear, conforming finite element method in space. The time mesh is graded appropriately near $t=0$, but the spatial mesh is quasiuniform. Previously, we proved that the error, measured in the spatial $L_2$-norm, is of order $k^{2+{\alpha }_{-}}+h^2\ell (k)$, uniformly in $t$, where $k$ is the maximum time step, $h$ is the maximum diameter of the spatial finite elements, ${\alpha }_{-}=\min (\alpha ,0)\le 0$ and $\ell (k)=\max (1,|\log k|)$. Here, we generalize a known result for the classical heat equation (i.e., the case $\alpha =0$) by showing that at each time level $t_n$ the solution is superconvergent with respect to $k$: the error is of order $(k^{3+2{\alpha }_{-}}+h^2)\ell (k)$. Moreover, a simple postprocessing step employing Lagrange interpolation yields a superconvergent approximation for any $t$. Numerical experiments indicate that our theoretical error bound is pessimistic if $\alpha <0$. Ignoring logarithmic factors, we observe that the error in the DG solution at $t=t_n$, and after postprocessing at all $t$, is of order $k^{3+{\alpha }_{-}}+h^2$.

24 pages, 2 figures

References in corpus (5)

• On an explicit finite difference method for fractional diffusion equations
• From Diffusion to Anomalous Diffusion: A Century after Einstein's Brownian Motion
• Fast and oblivious convolution quadrature
• A spectral order method for inverting sectorial Laplace transforms
• Fast summation by interval clustering for an evolution equation with memory

Cited by in corpus (10)

• The Galerkin Finite Element Method for A Multi-term Time-Fractional Diffusion equation
• Second-order numerical methods for multi-term fractional differential equations: Smooth and non-smooth solutions
• A Two-Grid Finite Element Approximation for A Nonlinear Time-Fractional Cable Equation
• Sharp $H^1$-norm error estimates of two time-stepping schemes for reaction-subdiffusion problems
• Time-stepping error bounds for fractional diffusion problems with non-smooth initial data
• A high-order nodal discontinuous Galerkin method for nonlinear fractional Schrödinger type equations
• A hybridized discontinuous Galerkin method for 2D fractional convection-diffusion equations
• A semidiscrete finite element approximation of a time-fractional Fokker-Planck equation with nonsmooth initial data
• Implementation of high-order, discontinuous Galerkin time stepping for fractional diffusion problems
• Uniform stability for a spatially-discrete, subdiffusive Fokker-Planck equation