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

JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS

卷:346

Article

Demay, Charles^1,2,3  Bourdarias, Christian^2,3  de Meux, Benoit de laage¹  Gerbi, Stephane^2,3  Herard, Jean-Marc^1,4,5 

[1] EDF R&D, 6 Quai Watier, F-78401 Chatou, France

[2] CNRS, UMR 5127, Lab Math, F-73376 Le Bourget Du Lac, France

[3] Univ Savoie Mt Blanc, F-73376 Le Bourget Du Lac, France

[4] Univ Aix Marseille, CNRS, UMR 7373, Inst Math Marseille, F-13453 Marseille, France

[5] Cent Marseille, F-13453 Marseille, France

关键词: Implicit-explicit scheme;    Operator splitting;    Two-layer model;    Gas-liquid flows;   

DOI  :  10.1016/j.cam.2018.06.027

来源: Elsevier

【 摘 要 】

This paper is devoted to the numerical simulation of the compressible two-layer model developed in Demay and Herard (2017). The latter is a 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 sub-systems. 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. Focusing on material waves of the convective system, the efficiency of the latter is improved when using material steps. However, considering the whole system with relaxation source terms, an efficient approximation of slow dynamics, typically a gravity driven flow, is still challenging. (C) 2018 Elsevier B.V. All rights reserved.

【 授权许可】

Free