We propose an entropy stable and positivity preserving discontinuous Galerkin (DG) scheme for the Euler equations with gravity, which is also well-balanced for hydrostatic equilibrium states. To achieve these properties, we utilize the nodal DG framework and carefully design the source term discretization using entropy conservative fluxes. Furthermore, we demonstrate that the proposed methodology is compatible with a positivity preserving scaling limiter, ensuring positivity of density and pressure under an appropriate CFL condition. To the best of our knowledge, this is the first DG scheme to simultaneously achieve these three properties with theoretical justification. Numerical examples further demonstrate its robustness and efficiency.
Liu et al. (Sat,) studied this question.