A simple HLLC-type Riemann solver for compressible non-equilibrium two-phase flows

Damien Furfaro (1) , Richard Saurel (2, 3)

1 AMU - Aix Marseille Université
2 IUSTI - Institut universitaire des systèmes thermiques industriels
3 M2P2 - Laboratoire de Mécanique, Modélisation et Procédés Propres

A simple, robust and accurate HLLC-type Riemann solver for two-phase 7-equation type models is built. It involves 4 waves per phase, i.e. the three conventional right- and left-facing and contact waves, augmented by an extra “interfacial” wave. Inspired by the Discrete Equations Method (Abgrall and Saurel, 2003), this wave speed (uIuI) is assumed function only of the piecewise constant initial data. Therefore it is computed easily from these initial states. The same is done for the interfacial pressure PIPI. Interfacial variables uIuI and PIPI are thus local constants in the Riemann problem. Thanks to this property there is no difficulty to express the non-conservative system of partial differential equations in local conservative form. With the conventional HLLC wave speed estimates and the extra interfacial speed uIuI, the four-waves Riemann problem for each phase is solved following the same strategy as in Toro et al. (1994) for the Euler equations. As uIuI and PIPI are functions only of the Riemann problem initial data, the two-phase Riemann problem consists in two independent Riemann problems with 4 waves only. Moreover, it is shown that these solvers are entropy producing. The method is easy to code and very robust. Its accuracy is validated against exact solutions as well as experimental data.

Mots clés: Multiphase; Hyperbolic; Non-conservative; HLLC solver

Domaines: Mécanique des fluides [physics.class-ph]

Damien Furfaro, Richard Saurel. A simple HLLC-type Riemann solver for compressible non-equilibrium two-phase flows. Computers and Fluids, 2015, 111, pp.159-178. ⟨10.1016/j.compfluid.2015.01.016⟩. ⟨hal-01278892⟩