New Preprint ‘Entropy-Stable and Well-Balanced Discontinuous Galerkin Methods for the Compressible Euler Equations in Vector-Invariant Form’ on arXiv

Marco Artiano, Kieran Ricardo, Oswald Knoth, Peter Spichtinger, and I have published a new preprint Entropy-Stable and Well-Balanced Discontinuous Galerkin Methods for the Compressible Euler Equations in Vector-Invariant Form on arXiv.

We develop structure-preserving methods for the compressible Euler equations with gravity in vector-invariant form and potential temperature as a prognostic variable within the flux-differencing discontinuous Galerkin spectral element method (DGSEM) framework. By discretizing the nonconservative terms as symmetric and antisymmetric products, we derive two-point numerical fluxes that conserve both the thermodynamic entropy and the total energy. Moreover, we design an entropy-stable numerical flux that is well-balanced for both isothermal and isentropic background states. All properties are shown to carry over to the high-order DGSEM on general curvilinear meshes. Several numerical examples confirm the theoretical findings and show the robustness and accuracy of the scheme for use in modern dynamical cores for atmospheric flows.

The reproducibility repository is available on GitHub and Zenodo.