An Exactly Curl-Free Finite-Volume/Finite-Difference Scheme for a Hyperbolic Compressible Isentropic Two-Phase Model

We present a new second order accurate structure-preserving finite volume scheme for the solution of the compressible isentropic two-phase model of Romenski et. al (Romenski J Sci Comput. 42, 68-95 2010), Romenski Q Appl Math 65(2), 259-279 2007) in multiple space dimensions. The governing equations fall into the wider class of symmetric hyperbolic and thermodynamically compatible systems and consist of a set of first-order hyperbolic partial differential equations. In the absence of algebraic source terms, the model is subject to a curl-free constraint for the relative velocity between the two phases. The main objective of this paper is therefore to preserve this structural property exactly also at the discrete level. The new numerical method is based on a staggered grid arrangement where the relative velocity field is stored in the cell vertexes, while all the remaining variables are stored in the cell centers. This allows the definition of discretely compatible gradient and curl operators which ensure that the discrete curl errors of the relative velocity field remain zero up to machine precision. A set of numerical results confirms this property also experimentally.