A positivity-preserving Eulerian two-phase approach with thermal relaxation for compressible flows with a liquid and gases

TL;DR

This paper introduces a positivity-preserving Eulerian two-phase numerical method for compressible flows involving liquids and gases, ensuring physical consistency and robustness in simulations with phase transitions.

Contribution

It develops a fractional algorithm combining hyperbolic model steps with thermal relaxation, enabling high-order schemes to preserve physical properties in complex flow simulations.

Findings

The algorithm maintains positivity of densities and sound speed.

It guarantees conservative updates suitable for phase transition problems.

Numerical tests demonstrate high accuracy and robustness.

Abstract

A positivity-preserving fractional algorithm is presented for solving the four-equation homogeneous relaxation model (HRM) with an arbitrary number of ideal gases and a liquid governed by the stiffened gas equation of state. The fractional algorithm consists of a time step of the hyperbolic five-equation model by Allaire et al. and an algebraic numerical thermal relaxation step at an infinite relaxation rate. Interpolation and flux limiters are proposed for the use of high-order Cartesian finite difference or finite volume schemes in a general form such that the positivity of the partial densities and squared sound speed, as well as the boundedness of the volume fractions and mass fractions, are preserved with the algorithm. A conservative solution update for the four-equation HRM is also guaranteed by the algorithm which is advantageous for certain applications such as those with phase transition. The accuracy and robustness of the algorithm with a high-order explicit finite difference weighted compact nonlinear scheme (WCNS) using the incremental-stencil weighted essentially non-oscillatory (WENO) interpolation, are demonstrated with various numerical tests.