A hybrid scheme to compute contact discontinuities in one-dimensional Euler systems

The present paper is devoted to the computation of single phase or two phase flows using the single-fluid approach. Governing equations rely on Euler equations which may be supplemented by conservation laws for mass species. Emphasis is given on numerical modelling with help of Godunov scheme or an approximate form of Godunov scheme called VFRoe-ncv based on velocity and pressure variables. Three distinct classes of closure laws to express the internal energy in terms of pressure, density and additional variables are exhibited. It is shown first that a standard conservative formulation of above mentioned schemes enables to predict “perfectly” unsteady contact discontinuities on coarse meshes, when the equation of state (EOS) belongs to the first class. On the basis of previous work issuing from literature, an almost conservative though modified version of the scheme is proposed to deal with EOS in the second or third class. Numerical evidence shows that the accuracy of approximations of discontinuous solutions of standard Riemann problems is strengthened on coarse meshes, but that convergence towards the right shock solution may be lost in some cases involving complex EOS in the third class. Hence, a blend scheme is eventually proposed to benefit from both properties (“perfect” representation of contact discontinuities on coarse meshes, and correct convergence on finer meshes). Computational results based on an approximate Godunov scheme are provided and discussed.

[1]  J. Smoller Shock Waves and Reaction-Diffusion Equations , 1983 .

[2]  P. Roe Approximate Riemann Solvers, Parameter Vectors, and Difference Schemes , 1997 .

[3]  S. Osher,et al.  A Non-oscillatory Eulerian Approach to Interfaces in Multimaterial Flows (the Ghost Fluid Method) , 1999 .

[4]  Jean-Marc Hérard,et al.  Some recent finite volume schemes to compute Euler equations using real gas EOS , 2002 .

[5]  Rémi Abgrall,et al.  A Simple Method for Compressible Multifluid Flows , 1999, SIAM J. Sci. Comput..

[6]  Arun In,et al.  Numerical Evaluation of an Energy Relaxation Method for Inviscid Real Fluids , 1999, SIAM J. Sci. Comput..

[7]  Smadar Karni,et al.  Hybrid Multifluid Algorithms , 1996, SIAM J. Sci. Comput..

[8]  Keh-Ming Shyue,et al.  Regular Article: A Fluid-Mixture Type Algorithm for Compressible Multicomponent Flow with van der Waals Equation of State , 1999 .

[9]  Rémi Abgrall,et al.  Computations of compressible multifluids , 2001 .

[10]  S. Clerc,et al.  Accurate computation of contact discontinuities in flows with general equations of state , 1999 .

[11]  Jean-Marc Hérard,et al.  A sequel to a rough Godunov scheme: application to real gases , 2000 .

[12]  Smadar Karni,et al.  Multicomponent Flow Calculations by a Consistent Primitive Algorithm , 1994 .

[13]  Raymond F. Bishop,et al.  Thermo-fluid Dynamic Theory of Two-Phase Flow , 1975 .

[14]  R. LeVeque Numerical methods for conservation laws , 1990 .

[15]  S. Clerc,et al.  Numerical Simulation of the Homogeneous Equilibrium Model for Two-Phase Flows , 2000 .

[16]  Michel Barret,et al.  Computation of Flashing Flows In Variable Cross-Section Ducts , 2000 .

[17]  T. Hou,et al.  Why nonconservative schemes converge to wrong solutions: error analysis , 1994 .

[18]  M. Ishii Thermo-fluid dynamic theory of two-phase flow , 1975 .

[19]  Steven F. Son,et al.  Two-phase modeling of DDT: Structure of the velocity-relaxation zone , 1997 .

[20]  Grégoire Allaire,et al.  A five-equation model for the numerical simulation of interfaces in two-phase flows , 2000 .

[21]  Michel Barret,et al.  Schemes to Compute Unsteady Flashing Flows , 2002 .

[22]  P. Raviart,et al.  Numerical Approximation of Hyperbolic Systems of Conservation Laws , 1996, Applied Mathematical Sciences.

[23]  Grégoire Allaire,et al.  A five-equation model for the simulation of interfaces between compressible fluids , 2002 .

[24]  B. Perthame,et al.  Relaxation of Energy and Approximate Riemann Solvers for General Pressure Laws in Fluid Dynamics , 1998 .

[25]  Frédéric Lagoutière,et al.  Modelisation mathematique et resolution numerique de problemes de fluides compressibles a plusieurs constituants , 2000 .

[26]  J. Kestin,et al.  Physical aspects of the relaxation model in two-phase flow , 1990, Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences.

[27]  R. Abgrall How to Prevent Pressure Oscillations in Multicomponent Flow Calculations , 1996 .

[28]  Richard Saurel,et al.  Mathematical and numerical modeling of two-phase compressible flows with micro-inertia , 2002 .

[29]  Stanley Osher,et al.  Level Set Methods , 2003 .

[30]  R. LeVeque Approximate Riemann Solvers , 1992 .

[31]  E. Toro Riemann Solvers and Numerical Methods for Fluid Dynamics , 1997 .

[32]  T. Gallouët,et al.  ON THE USE OF SYMMETRIZING VARIABLES TO DEAL WITH VACUUM , 2001 .

[33]  I. Bohachevsky,et al.  Finite difference method for numerical computation of discontinuous solutions of the equations of fluid dynamics , 1959 .

[34]  R. Abgrall,et al.  A Multiphase Godunov Method for Compressible Multifluid and Multiphase Flows , 1999 .

[35]  Rémi Abgrall,et al.  Ghost-Fluids for the Poor: A Single Fluid Algorithm for Multifluids , 2001 .