Numerical simulation of a compressible two-layer model: a first attempt with an implicit-explicit splitting scheme

摘要

This work is devoted to the numerical simulation of the compressible two-layer model developed in [16]. The latter is an hyperbolic two-fluid two-pressure model dedicated to gas-liquid flows in pipes, especially stratified air-water flows. Using explicit schemes, one obtains a CFL condition based on the celerity of (fast) acoustic waves which typically brings large numerical diffusivity for the (slow) material waves and small time steps. In order to overcome these drawbacks, the proposed scheme involves an operator splitting and an implicit-explicit time discretization. Thus, the full system is split into two hyperbolic subsystems. The first one deals with the transport equation on the liquid height using an explicit scheme and upwind fluxes. The second one deals with the averaged mass and momentum conservation equations of both phases using an implicit scheme which handles the propagation of acoustic waves. At last, the positivity of heights and densities is ensured under a CFL condition which involves material velocities. Numerical experiments are performed using acoustic as well as material time steps. Adding the Rusanov scheme for comparison, the best accuracy is obtained with the proposed scheme used with acoustic time steps. When used with material time steps, efficiency on the slow waves and stability are obtained regarding analytical solutions of the convective part.