Global convergence of \u3ci\u3ea posteriori\u3c/i\u3e error estimates for a discontinuous Galerkin method for one-dimensional linear hyperbolic problems

In this paper we study the global convergence of the implicit residual-based a posteriori error estimates for a discontinuous Galerkin method applied to one-dimensional linear hyperbolic problems. We apply a new optimal superconvergence result [Y. Yang and C.-W. Shu, SIAM J. Numer. Anal., 50 (2012), pp. 3110-3133] to prove that, for smooth solutions, these error estimates at a fixed time converge to the true spatial errors in the L2 -norm under mesh refinement. The order of convergence is proved to be k + 2, when k-degree piecewise polynomials with k ≥ 1 are used. As a consequence, we prove that the DG method combined with the a posteriori error estimation procedure yields both accurate error estimates and O(hk+2) superconvergent solutions. We perform numerical experiments to demonstrate that the rate of convergence is optimal. We further prove that the global effectivity indices in the L2 -norm converge to unity under mesh refinement. The order of convergence is proved to be 1. 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 k + 3/2 and 1/2, respectively. Our proofs are valid for arbitrary regular meshes using Pk polynomials with k ≥ 1 and for both the periodic boundary condition and the initial-boundary value problem. Several numerical simulations are performed to validate the theory

Analysis of optimal error estimates and superconvergence of the discontinuous Galerkin method for convection-diffusion problems in one space dimension

In this paper, we study the convergence and superconvergence properties of the discontinuous Galerkin (DG) method for a linear convection-diffusion problem in one-dimensional setting. We prove that the DG solution and its derivative exhibit optimal O(hp+1) and O(hp) convergence rates in the L 2 -norm, respectively, when p-degree piecewise polynomials with p ≥ 1 are used. We further prove that the p-degree DG solution and its derivative are O(h2p) superconvergent at the downwind and upwind points, respectively. Numerical experiments demonstrate that the theoretical rates are optimal and that the DG method does not produce any oscillation. We observed optimal rates of convergence and superconvergence even in the presence of boundary layers when Shishkin meshes are used

Superconvergence and \u3ci\u3ea posteriori\u3c/i\u3e error estimates of a local discontinuous Galerkin method for the fourth-order initial-boundary value problems arising in beam theory

In this paper, we investigate the superconvergence properties and a posteriori error estimates of a local discontinuous Galerkin (LDG) method for solving the one-dimensional linear fourth-order initial-boundary value problems arising in study of transverse vibrations of beams. We present a local error analysis to show that the leading terms of the local spatial discretization errors for the k-degree LDG solution and its spatial derivatives are proportional to (k + 1)-degree Radau polynomials. Thus, the k-degree LDG solution and its derivatives are O(hk+2) superconvergent at the roots of (k + 1)-degree Radau polynomials. Computational results indicate that global superconvergence holds for LDG solutions. We discuss how to apply our superconvergence results to construct efficient and asymptotically exact a posteriori error estimates in regions where solutions are smooth. Finally, we present several numerical examples to validate the superconvergence results and the asymptotic exactness of our a posteriori error estimates under mesh refinement. Our results are valid for arbitrary regular meshes and for Pk polynomials with k ≥ 1, and for various types of boundary conditions

Discontinuous Galerkin methods for convection-diffusion equations and applications in petroleum engineering

This dissertation contains research in discontinuous Galerkin (DG) methods applying to convection-diffusion equations. It contains both theoretical analysis and applications. Initially, we develop a conservative local discontinuous Galerkin (LDG) method for the coupled system of compressible miscible displacement problem in two space dimensions. The main difficulty is how to deal with the discontinuity of approximations of velocity, u, in the convection term across the cell interfaces. To overcome the problems, we apply the idea of LDG with IMEX time marching using the diffusion term to control the convection term. Optimal error estimates in Linfinity(0, T; L2) norm for the solution and the auxiliary variables will be derived. Then, high-order bound-preserving (BP) discontinuous Galerkin (DG) methods for the coupled system of compressible miscible displacements on triangular meshes will be developed. There are three main difficulties to make the concentration of each component between 0 and 1. Firstly, the concentration of each component did not satisfy a maximum-principle. Secondly, the first-order numerical flux was difficult to construct. Thirdly, the classical slope limiter could not be applied to the concentration of each component. To conquer these three obstacles, we first construct special techniques to preserve two bounds without using the maximum-principle-preserving technique. The time derivative of the pressure was treated as a source of the concentration equation. Next, we apply the flux limiter to obtain high-order accuracy using the second-order flux as the lower order one instead of using the first-order flux. Finally, L2-projection of the porosity and constructed special limiters that are suitable for multi-component fluid mixtures were used. Lastly, a new LDG method for convection-diffusion equations on overlapping mesh introduced in [J. Du, Y. Yang and E. Chung, Stability analysis and error estimates of local discontinuous Galerkin method for convection-diffusion equations on overlapping meshes, BIT Numerical Mathematics (2019)] showed that the convergence rates cannot be improved if the dual mesh is constructed by using the midpoint of the primitive mesh. They provided several ways to gain optimal convergence rates but the reason for accuracy degeneration is still unclear. We will use Fourier analysis to analyze the scheme for linear parabolic equations with periodic boundary conditions in one space dimension. To investigate the reason for the accuracy degeneration, we explicitly write out the error between the numerical and exact solutions. Moreover, some superconvergence points that may depend on the perturbation constant in the construction of the dual mesh were also found out

AN OPTIMAL A POSTERIORI ERROR ESTIMATES OF THE LOCAL DISCONTINUOUS GALERKIN METHOD FOR THE SECOND-ORDER WAVE EQUATION IN ONE SPACE DIMENSION

In this paper, we provide the optimal convergence rate of a posteriori error estimates for the local discontinuous Galerkin (LDG) method for the second-order wave equation in one space dimension. One of the key ingredients in our analysis is the recent optimal superconvergence result in [W. Cao, D. Li and Z. Zhang, Commun. Comput. Phys. 21 (1) (2017) 211-236]. 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. We use these results 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. These new results enable us to construct residual-based a posteriori error estimates of the spatial errors. We further prove that, for smooth solutions, these a posteriori LDG error estimates converge, at a fixed time, to the true spatial errors in the L 2 -norm at O(h p+2) rate. Finally, we show that the global effectivity indices in the L 2 -norm converge to unity at O(h) rate. The current 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