Analysis of Godunov type schemes applied to the compressible Euler system at low Mach number

Abstract

We propose a theoretical framework to clearly explain the inaccuracy of Godunov type schemes applied to the compressible Euler system at low Mach number on a Cartesian mesh. In particular, we clearly explain why this inaccuracy problem concerns the 2D or 3D geometry and does not concern the 1D geometry. The theoretical arguments are based on the Hodge decomposition, on the fact that an appropriate well-prepared subspace is invariant for the linear wave equation and on the notion of first-order modified equation. This theoretical approach allows to propose a simple modification that can be applied to any colocated scheme of Godunov type or not in order to define a large class of colocated schemes accurate at low Mach number on any mesh. It also allows to justify colocated schemes that are accurate at low Mach number as, for example, the Roe-Turkel and the AUSM+-up schemes, and to find a link with a colocated incompressible scheme stabilized with a Brezzi-Pitkäranta type stabilization. Numerical results justify the theoretical arguments proposed in this paper. © 2009 Elsevier Inc. All rights reserved.