On the order reduction of entropy stable DGSEM for the compressible Euler equations

DOI10.1007/978-3-030-39647-3_2MaRDI QIDQ2056464zbMATH OpenOpenAlexFDO

Authors Florian J. Hindenlang, Gregor J. Gassner

Publication date 2 December 2021

Copyright license Creative Commons Attribution 4.0 International

Full work available at URL https://arxiv.org/abs/1901.05812

zbMATH Keywords

compressible Euler equationsentropy stabilityorder reductiondiscontinuous Galerkin spectral element methodsummation-by-parts

Mathematics Subject Classification ID

Euler equations (35Q31) Multistep, Runge-Kutta and extrapolation methods for ordinary differential equations (65L06) Stability and convergence of numerical methods for initial value and initial-boundary value problems involving PDEs (65M12) Spectral, collocation and related methods for boundary value problems involving PDEs (65N35) Spectral methods applied to problems in fluid mechanics (76M22) Existence, uniqueness, and regularity theory for compressible fluids and gas dynamics (76N10) Spectral, collocation and related methods for initial value and initial-boundary value problems involving PDEs (65M70)

Abstract: Is the experimental order of convergence lower when using the entropy stable DGSEM-LGL variant? Recently, a debate on the question of the convergence behavior of the entropy stable nodal collocation discontinuous Galerkin spectral element method (DGSEM) with Legendre-Gauss-Lobatto nodes has emerged. Whereas it is well documented that the entropy conservative variant with no additional interface dissipation shows an odd-even behavior when testing its experimental convergence order, the results in the literature are less clear regarding the entropy stable version of the DGSEM-LGL, where explicit Riemann solver type dissipation is added at the element interfaces. We contribute to the ongoing discussion and present numerical experiments for the compressible Euler equations, where we investigate the effect of the choice of the numerical surface flux function. In our experiments, it turns out that the choice of the numerical surface flux has an impact on the convergence order. Penalty type numerical fluxes with high dissipation in all waves, such as the LLF and the HLL flux, appear to affect the convergence order negatively for odd polynomial degrees $N$, in contrast to the entropy conserving variant, where even polynomial degrees $N$ are negatively affected. This behavior is more pronounced in low Mach number settings. In contrast, for numerical surface fluxes with less dissipative behavior in the contact wave such as e.g. Roe's flux, the HLLC flux and the entropy conservative flux augmented with 5-wave matrix dissipation, optimal convergence rate of $N+1$ independent of the Mach number is observed.

Recommendations

• An analysis of discontinuous Galerkin methods for the compressible Euler equations: entropy and \(L_2\) stability
• On the entropy conserving/stable implicit DG discretization of the Euler equations in entropy variables
• Convergence of discontinuous Galerkin schemes for the Euler equations via dissipative weak solutions
• On discretely entropy conservative and entropy stable discontinuous Galerkin methods
• Entropy stable discontinuous Galerkin schemes on moving meshes for hyperbolic conservation laws

Cited in

(3)

• Entropy-stable discontinuous Galerkin approximation with summation-by-parts property for the incompressible Navier-Stokes/Cahn-Hilliard system
• Skew-symmetric entropy stable modal discontinuous Galerkin formulations
• Convergence of discontinuous Galerkin schemes for the Euler equations via dissipative weak solutions

This page was built for publication: On the order reduction of entropy stable DGSEM for the compressible Euler equations

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2056464)