We present a numerical approximation of the solutions of the Euler equations with a gravitational source term. On the basis of a Suliciu type relaxation model with two relaxation speeds, we construct an approximate Riemann solver, which is used in a first order Godunov-type finite volume scheme. This scheme can preserve both stationary solutions and the low Mach limit to the corresponding incompressible equations. In addition, we prove that our scheme preserves the positivity of density and internal energy, that it is entropy satisfying and also guarantees not to give rise to numerical checkerboard modes in the incompressible limit. Later we give an extension to second order that preserves these properties. Finally, the theoretical properties are investigated in numerical experiments.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.