Superconvergence of discontinuous Galerkin method for scalar and vector linear advection equations

TL;DR

This paper analyzes the superconvergence properties of the discontinuous Galerkin method for scalar and vector linear advection equations using Fourier analysis, deriving error bounds and confirming results through numerical experiments.

Contribution

It provides a detailed Fourier analysis of superconvergence in DG methods for advection equations, extending results to vector laws and comparing two error computation approaches.

Findings

Superconvergence of order 2k+1 established for scalar advection equations.

Error bounds confirmed by numerical experiments.

Superconvergence extends to vector conservation laws with Lax-Friedrichs flux.

Abstract

In this paper, we use Fourier analysis to study the superconvergence of the semi-discrete discontinuous Galerkin method for scalar linear advection equations in one spatial dimension. The error bounds and asymptotic errors are derived for initial discretization by $L_2$ projection, Gauss-Radau projection, and other projections proposed by Cao et. al. For pedagogical purpose, the errors are computed in two different ways. In the first approach, we compute the difference between the numerical solution and a special interpolation of the exact solution, and show that it consists of an asymptotic error of order $2k+1$ and a transient error of lower order. In the second approach, as by Chalmers and Krivodonova, we compute the error directly by decomposition into physical and nonphysical modes, and obtain agreement with the first approach. We then extend the analysis to vector conservation laws, solved using the Lax-Friedrichs flux. We prove that the superconvergence holds with the same order. The error bounds and asymptotic errors are demonstrated by various numerical experiments for scalar and vector advection equations.