Simple and efficient relaxation methods for interfaces separating compressible fluids, cavitating flows and shocks in multiphase mixtures

A hyperbolic two-phase flow model involving five partial differential equations is constructed for liquid-gas interface modelling. The model is able to deal with interfaces of simple contact where normal velocity and pressure are continuous as well as transition fronts where heat and mass transfer occur, involving pressure and velocity jumps. These fronts correspond to extra waves in the system. The model involves two temperatures and entropies but a single pressure and a single velocity. The closure is achieved by two equations of state that reproduce the phase diagram when equilibrium is reached. Relaxation toward equilibrium is achieved by temperature and chemical potential relaxation terms whose kinetics is considered infinitely fast at specific locations only, typically at evaporation fronts. Thus, metastable states are involved for locations far from these fronts. Computational results are compared to the experimental ones. Computed and measured front speeds are of the same order of magnitude and the same tendency of increasing front speed with initial temperature is reported. Moreover, the limit case of evaporation fronts propagating in highly metastable liquids with the Chapman-Jouguet speed is recovered as an expansion wave of the present model in the limit of stiff thermal and chemical relaxation.

/pdf/modelling-phase-transition-in-metastable-liquids-application-3q040zs8m5.pdf

Modelling phase transition in metastable liquids: Application to cavitating and flashing flows.

https://www.math.u-bordeaux.fr/~rabgrall/Articles/dipha.pdf

Discrete equations for physical and numerical compressible multiphase mixtures

A robust and accurate approach to computing compressible multiphase flow: Stratified flow model and AUSM+-up scheme

Modeling phase transition for compressible two-phase flows applied to metastable liquids

The aim of this paper is the derivation of a multiphase model of compressible fluids. Each fluid has a different average translational velocity, density, pressure, internal energy as well as the energies related to rotation and vibration. The main difficulty is the description of these various translational, rotational and vibrational motions in the context of a one-dimensional model. The second difficulty is the determination of closure relations for such a system: the 'drag' force between inviscid fluids, pressure relaxation rate, vibration and rotation creation rates, etc. The rotation creation rate is particularly important for turbulent flows with shock waves. In order to derive the one-dimensional multiphase model, two different approaches are used. The first one is based on the Hamilton principle. We use the second approach, in which the pure fluid equations are discretized at the microscopic level and then averaged. In this context, the flow is considered to be the annular flow of two turbulent fluids. We also derive the continuous limit of this model which provides explicit formulae for the closure laws

A multiphase model with internal degrees of freedom: application to shock–bubble interaction

This paper is devoted to the modeling and numerical resolution of a non-dissipative compressible turbulent plasma flow model involving three temperatures (turbulence, ions and electrons). The first step is to derive such a model. To do this, an analysis of the Reynolds averaged Euler equations (the k-model) is carried out. It is shown that thermodynamic requirements enable the derivation of an equation of state for turbulent variables. This equation of state is of the same type as those of an ideal gas. In this context, the various thermodynamic variables of turbulence can be obtained (energy, pressure, temperature etc.). This hyperbolic conservative model has exactly the same structure as the two temperatures plasma model of Zel’dovich. Thanks to the clear structure of these two models, the turbulent plasma model is derived and involves three temperatures. The second step is to derive an accurate numerical scheme for its solution. A linearized Riemann solver and a positive HLLC type solver are derived and embedded into a conventional Godunov scheme. It is shown that this method requires important corrections to preserve contact discontinuities and temperatures monotonicity. The corrections are based upon a non-conservative formulation of the turbulence and electrons energy equations, while total energy conservation is preserved. The modified method behaves correctly with contacts and shocks.

Thermodynamic analysis and numerical resolution of a turbulent - fully ionized plasma flow model

This paper is devoted to the numerical approximation of a hyperbolic non-equilibrium multiphase flow model with different velocities on moving meshes. Such goal poses several difficulties. The presence of different flow velocities in conjunction with cell velocities poses difficulties for upwinding fluxes. Another issue is related to the presence of non-conservative terms. To solve these difficulties, the discrete equations method (J. Comput. Phys. 2003; 186(2):361–396; J. Fluid. Mech. 2003; 495:283–321; J. Comput. Phys. 2004; 196:490–538; J. Comput. Phys. 2005; 205:567–610) is employed and generalized to the context of moving cells. The complementary conservation laws, available for the mixture, are used to determine the velocities of the cells boundaries. With these extensions, an accurate and robust multiphase flow method on moving meshes is obtained and validated over several test problems with exact or experimental solutions. Copyright © 2007 John Wiley & Sons, Ltd.

François Renaud

Papers

A multiphase model with internal degrees of freedom: application to shock–bubble interaction

Thermodynamic analysis and numerical resolution of a turbulent - fully ionized plasma flow model

A numerical method for one‐dimensional compressible multiphase flows on moving meshes