There is a qualitative difference between one-dimensional and multidimensional solutions to the Euler equations: new features that arise are vorticity and a nontrivial incompressible (low Mach number) limit. They present challenges to finite volume methods. It seems that an important step in this direction is to first study the new features for the multi-dimensional acoustic equations. There exists an analogue of the low Mach number limit for this system and its vorticity is stationary.It is shown that a scheme that possesses a stationary discrete vorticity (vorticity preserving) also has stationary states that are discretizations of all the analytic stationary states. This property is termed stationarity preserving. Both these features are not generically fulfilled by finite volume schemes; in this paper a condition is derived that determines whether a scheme is stationarity preserving (or, equivalently, vorticity preserving) on a Cartesian grid.Additionally, this paper also uncovers a previously unknown connection to schemes that comply with the low Mach number limit. Truly multi-dimensional schemes are found to arise naturally and it is shown that a multi-dimensional discrete divergence previously discussed in the literature is the only possible stationarity-preserving one (in a certain class).

Finite volume schemes for hyperbolic conservation laws require a numerical intercell flux. In one spatial dimension the numerical flux can be successfully obtained by solving (exactly or approximately) Riemann problems that are introduced at cell interfaces. This is more challenging in multiple spatial dimensions. The active flux scheme is a finite volume scheme that considers continuous reconstructions instead. The intercell flux is obtained using additional degrees of freedom distributed along the cell boundary. For their time evolution an exact evolution operator is employed, which naturally ensures the correct direction of information propagation and provides stability. This paper presents an implementation of active flux for the acoustic equations on two-dimensional Cartesian grids and demonstrates its ability to simulate discontinuous solutions with an explicit time stepping in a stable manner. Additionally, it is shown that the active flux scheme for linear acoustics is low Mach number compliant without the need for any fix.

Based on the Roe solver a new technique that allows to correctly represent low Mach number flows with a discretization of the compressible Euler equations was proposed in [22]. We analyze properties of this scheme and demonstrate that its limit yields a discretization of the continuous limit system. Furthermore we perform a linear stability analysis for the case of explicit time integration and study the performance of the scheme under implicit time integration via the evolution of its condition number. A numerical implementation demonstrates the capabilities of the scheme on the example of the Gresho vortex which can be accurately followed down to Mach numbers of ∼ 10 −10 .

No abstract