Superconvergence of Discontinuous Galerkin method for linear hyperbolic equations

Cao Waixiang,Zhang Zhimin,Zou Qingsong
DOI: https://doi.org/10.48550/arXiv.1311.6938
2013-11-27
Abstract:In this paper, we study superconvergence properties of the discontinuous Galerkin (DG) method for one-dimensional linear hyperbolic equation when upwind fluxes are used. We prove, for any polynomial degree $k$, the $2k+1$th (or $2k+1/2$th) superconvergence rate of the DG approximation at the downwind points and for the domain average under quasi-uniform meshes and some suitable initial discretization. Moreover, we prove that the derivative approximation of the DG solution is superconvergent with a rate $k+1$ at all interior left Radau points. All theoretical finding are confirmed by numerical experiments.
Numerical Analysis
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the super - convergence property of the discontinuous Galerkin (DG) method when solving one - dimensional linear hyperbolic equations. Specifically, the authors aim to prove the super - convergence order of the DG method at certain special points (such as the downwind point and Radau points) and the regional average when using the upwind flux. ### Main problems: 1. **Super - convergence order**: For any polynomial degree \( k \), the authors hope to prove that the super - convergence order of the DG approximation at the downwind point and the regional average is \( 2k + 1 \) or \( 2k+\frac{1}{2} \). 2. **Derivative super - convergence**: In addition, they also hope to prove that the super - convergence order of the derivative of the DG solution at all internal left Radau points is \( k + 1 \). ### Research background: - **Discontinuous Galerkin method**: The DG method was initially used for neutron transport problems and has attracted wide attention due to its flexibility, high - efficiency parallel implementation ability, ability to handle complex geometric shapes or interfaces, and adaptability to arbitrary h - p self - adaptivity. - **Super - convergence research**: Over the past few decades, a great deal of research has been carried out on the super - convergence properties of the DG method, including ordinary differential equations, multi - dimensional first - order hyperbolic systems, one - dimensional hyperbolic conservation laws, and time - dependent convection - diffusion equations. ### Paper contributions: - **Theoretical proof**: Through rigorous mathematical proofs, the authors have solved the previously unsolved puzzle of the \( 2k + 1 \) super - convergence order. - **Simplified proof**: As a by - product, they provide a simplified method to prove the point - wise \( k + 2 \) - order super - convergence rate at the right Radau point. - **New results**: For the first time, they prove the point - wise \( k + 1 \) - order super - convergence rate of the derivative at the internal left Radau point. ### Method: - **Correction function design**: The authors introduce a correction function \( w \) and use this function to correct the error between the DG solution and the truncated Radau expansion to a high order. - **Special interpolation function**: A special interpolation function \( u_I \) is constructed so that it has a high - order super - closeness with the DG solution. ### Numerical experiments: - The theoretical results are verified by numerical experiments, and the influence of different initial discretization methods on the convergence rate is demonstrated. ### Conclusion: This paper successfully reveals the super - convergence properties of the DG method in one - dimensional linear hyperbolic equations, and provides rigorous mathematical proofs and numerical verifications, laying a solid foundation for future research.