Towards a fully well-balanced and entropy-stable scheme for the Euler equations with gravity: preserving isentropic steady solutions

TL;DR

This paper develops a numerical scheme for the Euler equations with gravity that is both well-balanced for equilibrium solutions and entropy-stable, ensuring accurate and physically admissible simulations.

Contribution

The paper introduces a novel scheme that preserves all steady and moving equilibrium solutions of the Euler equations with gravity while maintaining entropy stability and positivity.

Findings

The scheme exactly preserves all equilibrium solutions.

It satisfies discrete entropy inequalities.

Numerical experiments confirm the scheme's effectiveness.

Abstract

The present work concerns the derivation of a numerical scheme to approximate weak solutions of the Euler equations with a gravitational source term. The designed scheme is proved to be fully well-balanced since it is able to exactly preserve all moving equilibrium solutions, as well as the corresponding steady solutions at rest obtained when the velocity vanishes. Moreover, the proposed scheme is entropy-preserving since it satisfies all fully discrete entropy inequalities. In addition, in order to satisfy the required admissibility of the approximate solutions, the positivity of both approximate density and pressure is established. Several numerical experiments attest the relevance of the developed numerical method. An extension to two-dimensional problems is given, applying the one-dimensional framework direction by direction on Cartesian grids.