A positive and entropy-satisfying finite volume scheme for the Baer-Nunziato model

JOURNAL OF COMPUTATIONAL PHYSICS

卷:330

Article

Coquel, Frederic¹  Herard, Jean-Marc²  Saleh, Khaled³ 

[1] Ecole Polytech, CMAP, UMR 7641, CNRS, Route Saclay, F-91128 Palaiseau, France

[2] EDF R&D, Dept MFEE, 6 Quai Watier, F-78401 Chatou, France

[3] Univ Lyon 1, CNRS UMR 5208, Inst Camille Jordan, 43 Bd 11 Novembre 1918, F-69622 Villeurbanne, France

关键词: Compressible multi-phase flows;    Hyperbolic PDEs;    Energy-entropy duality;    Entropy-satisfying methods;    Relaxation techniques;    Riemann problem;    Riemann solvers;    Finite volumes;   

DOI  :  10.1016/j.jcp.2016.11.017

来源: Elsevier

【 摘 要 】

We present a relaxation scheme for approximating the entropy dissipating weak solutions of the Baer-Nunziato two-phase flow model. This relaxation scheme is straightforwardly obtained as an extension of the relaxation scheme designed in [16] for the isentropic Baer-Nunziato model and consequently inherits its main properties. To our knowledge, this is the only existing scheme for which the approximated phase fractions, phase densities and phase internal energies are proven to remain positive without any restrictive condition other than a classical fully computable CFL condition. For ideal gas and stiffened gas equations of state, real values of the phasic speeds of sound are also proven to be maintained by the numerical scheme. It is also the only scheme for which a discrete entropy inequality is proven, under a CFL condition derived from the natural sub characteristic condition associated with the relaxation approximation. This last property, which ensures the non-linear stability of the numerical method, is satisfied for any admissible equation of state. We provide a numerical study for the convergence of the approximate solutions towards some exact Riemann solutions. The numerical simulations show that the relaxation scheme compares well with two of the most popular existing schemes available for the Baer-Nunziato model, namely Schwendeman-Wahle-Kapila's Godunov-type scheme [39] and Tokareva-Toro's HLLC scheme [44]. The relaxation scheme also shows a higher precision and a lower computational cost (for comparable accuracy) than a standard numerical scheme used in the nuclear industry, namely Rusanov's scheme. Finally, we assess the good behavior of the scheme when approximating vanishing phase solutions. (C) 2016 Elsevier Inc. All rights reserved.

【 授权许可】

Free