A slope limiting procedure in discontinuous Galerkin finite element method for gasdynamics applications
Shuangzhang Tu 1, Shahrous Allibadi 11 223 James P. Brawley Dr. SW, Clark Atlanta University, Department of Engineering, Atlanta, GA 30314, USA
Received by the editors January 1, 2004 and, in revised form, April 22, 2004
In this paper we demonstrate the performance of a slope limiting procedure combined with a discontinuous Galerkin (DG) finite element solver for 2D compressible Euler equations. The slope limiter can be categorized into van Albada type and is differentiable. This slope limiter is modified from a similar limiter used in finite volume solvers to suit the needs of the DC solver. The gradient in an element is limited using the weighted average of the face gradients. The face gradients are obtained from the area-weighted average of the gradient on both sides of the faces. The slope limiting process is very suitable for meshes discretized by triangle elements. The HLLC (Harten, Lax and van Leer) or the local Lax-Friedrich (LLF) flux functions is used to compute the interface fluxes in the DG formulation. The second order TVD Runge-Kutta scheme is employed for the time integration. The numerical examples including transonic, supersonic and hypersonic flows show that the current slope limiting process together with the DG solver is able to remove overshoots and undershoots around high gradient regions while preserving the high accuracy of the DC method. The convergence histories of all examples demonstrate that the limiting process does not stall convergence to steady state as many other slope limiters do.
AMS subject classifications: 35L65, 76M10, 76N15
Key words: discontinuous Galerkin (DC); slope limiting and 2D inviscid compressible flows
Email: firstname.lastname@example.org (S. Tu), email@example.com (S. Allibadi)