A relaxation scheme for a hyperbolic multiphase flow model. Part I: barotropic eos

TL;DR

This paper develops a relaxation finite volume scheme for multiphase hyperbolic flow models with barotropic equations of state, demonstrating improved accuracy, stability, and efficiency over traditional schemes like Rusanov's.

Contribution

It extends a relaxation scheme from two-phase to N-phase flow models, ensuring positivity, stability, and better computational performance.

Findings

Relaxation scheme is more accurate than Rusanov's for the same refinement.

The scheme maintains positivity of phase fractions and densities.

It remains stable even with vanishing phases.

Abstract

This article is the first of two in which we develop a relaxation finite volume scheme for the convective part of the multiphase flow models introduced in the series of papers [10, 9, 4]. In the present article we focus on barotropic flows where in each phase the pressure is a given function of the density. The case of general equations of state will be the purpose of the second article. We show how it is possible to extend the relaxation scheme designed in [8] for the barotropic Baer-Nunziato two-phase flow model to the multiphase flow model with N-where N is arbitrarily large-phases. The obtained scheme inherits the main properties of the relaxation scheme designed for the Baer-Nunziato two phase flow model. The approximated phase fractions and phase densities are proven to remain positive and a discrete energy inequality is also proven under a classical CFL condition. For the same level of refinement, the relaxation scheme is shown to be much more accurate than Rusanov's scheme, and for a given level of approximation error, the relaxation scheme is shown to perform much better in terms of computational cost than Rusanov's scheme. Moreover, contrary to Rusanov's scheme which develops strong oscillations when approximating vanishing phase solutions, the numerical results show that the relaxation scheme remains stable in such regimes.