Numerical resolution of a two-component compressible fluid model with interfaces

Abstract

We study an algorithm for moving interfaces in a two-component compressible fluid model. We propose to use the limited downwind scheme developed in Despres and Lagoutiere (1999, 2002) to avoid artificial numerical spreading of interfaces. The numerical treatment of the mixture is shown to be free of spurious pressure and velocity oscillations near the contact discontinuity. The algorithm is conservative for partial masses, total impulse and total energy. Various numerical simulations show the interest of this approach, for interfaces in dimensions 1, 2 and 3. Simplicity of the coding is an important feature of the algorithm.