Optimal a Posteriori Error Estimates of the Local Discontinuous Galerkin Method for Convection-Diffusion Problems in One Space Dimension

Abstract

In this paper, we derive optimal order a posteriori error estimates for the local discontinuous Galerkin (LDG) method for linear convection-diffusion problems in one space dimension. One of the key ingredients in our analysis is the recent optimal superconvergence result in [Y. Yang and C.-W. Shu, J. Comp. Math., 33 (2015), pp. 323-340]. We first prove that the LDG solution and its spatial derivative, respectively, converge in the L 2-norm to (p + 1)-degree right and left Radau interpolating polynomials under mesh refinement. The order of convergence is proved to be p + 2, when piecewise polynomials of degree at most p are used. These results are used to show that the leading error terms on each element for the solution and its derivative are proportional to (p + 1)-degree right and left Radau polynomials. We further prove that, for smooth solutions, the a posteriori LDG error estimates, which were constructed by the author in an earlier paper, converge, at a fixed time, to the true spatial errors in the L 2-norm at O(h p+2) rate. Finally, we prove that the global effectivity indices in the L 2-norm converge to unity at O(h) rate. These results improve upon our previously published work in which the order of convergence for the a posteriori error estimates and the global effectivity index are proved to be p + 3/2 and 1/2, respectively. Our proofs are valid for arbitrary regular meshes using P p polynomials with p ≥ 1. Several numerical experiments are performed to validate the theoretical results.

A Bound-Preserving Numerical Scheme for Space–Time Fractional Advection Equations

“…However, both methods have a rate of time convergence of 1. To obtain a higher accuracy in the numerical method, the discontinuous Galerkin method [41][42][43] can be adopted for space and/or time derivatives, and this can lead to interesting outcomes. It could be that a compact finite difference method [44,45] can be adopted for time derivatives in our future work, which could thus lead to a better convergency.…”

Numerical solution of the Whitham-Broer-Kaup shallow water equation by quartic B-spline collocation method

“…It offers a robust and efficient numerical approach that is able to achieve high levels of accuracy. In our future research, we plan to explore the continuous and discontinuous Galerkin methods in both spatial and temporal dimensions [43][44][45][46][47][48][49][50][51] for solving the WBK problem, with the goal of achieving higher levels of accuracy in our solutions. Additionally, we intend to broaden the application of the B-spline collocation method to various challenging real-world problems including, timefractional coupled WBK problems, see [52,53], integro-differential beam problems [54], nonlinear Bratu problems [55], and delayed-differential equations [56].…”

Superconvergence Analysis of Discontinuous Galerkin Methods for Systems of Second-Order Boundary Value Problems

“…He proved that the p-degree DG solution achieves an O(h p+2 ) superconvergence rate at the roots of these particular polynomials. Baccouch and Temimi [20] further extended the DG error analysis to secondorder BVPs and showed that when using p-degree piecewise polynomials, the UWDG solution and its derivative exhibit a superconvergence rate of O(h 2p ) at the upwind and downwind endpoints. Moreover, Baccouch and Temimi [21] proposed a novel DG scheme for solving the wave equation based on the method of lines.…”