Modelling evaporation fronts with reactive Riemann solvers
- IUSTI, UMR CNRS 6595, Universite Aix-Marseille I, Technopole de Chateau Gombert, 5 Rue E. Fermi, 13453 Marseille Cedex 13 (France) and Projet SMASH, INRIA, 2004 route des Lucioles, 06902 Sophia Antipolis (France)
- IUSTI, UMR CNRS 6595, Universite Aix-Marseille I, Technopole de Chateau Gombert, 5 Rue E. Fermi, 13453 Marseille Cedex 13 (France) and Projet SMASH, INRIA, 2004 route des Lucioles, 06902 Sophia Antipolis (France) and Institut Universitaire de France, Technopole de Chateau Gombert, 5 Rue E. Fermi, 13453 Marseille Cedex 13 (France)
This work deals with the modelling of permeable fronts and the building of a numerical method allowing the multi-dimensional propagation of such fronts. A particular attention is given to evaporation waves that appear in cavitating systems. These ones are considered as discontinuities through which a non-equilibrium liquid turns to a liquid-vapor mixture at thermodynamic equilibrium. Such transformation occurs at finite rate. In order to determine this kinetics, the evaporation front is assumed to propagate at the maximum admissible speed corresponding to the Chapman-Jouguet deflagration point [J.R., Simoes-Moreira, J.E., Shepherd, Evaporation waves in superheated dodecane, J. Fluid Mech. 382 (1999) 63-86]. Using this particular kinetic relation, Rankine-Hugoniot relations are closed at such fronts. Then it is possible to solve the associated reactive Riemann problem. However, another difficulty is present to solve the multi-dimensional propagation of permeable fronts. This kind of front is subsonic and a conventional averaging scheme (such as Godunov scheme) is inappropriate. To overcome this difficulty, the reactive Riemann problem solution is embedded into the discrete equations method (DEM) [R., Abgrall, R., Saurel, Discrete equations for physical and numerical compressible multiphase mixtures, J. Comp. Phys. 186 (2003) 361-396; R., Saurel, S., Gavrilyuk, F., Renaud, A multiphase model with internal degrees of freedom: application to Shock-Bubble Interaction, J. Fluid. Mech., 495 (2003) 283-321]. This numerical method necessitates deep extensions that are detailed herein. Numerical results are shown and validated over experimental data. Some examples show that the same method may be applied to the propagation of detonation fronts.
- OSTI ID:
- 20687234
- Journal Information:
- Journal of Computational Physics, Vol. 205, Issue 2; Other Information: DOI: 10.1016/j.jcp.2004.11.021; PII: S0021-9991(04)00481-4; Copyright (c) 2004 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved; Country of input: International Atomic Energy Agency (IAEA); ISSN 0021-9991
- Country of Publication:
- United States
- Language:
- English
Similar Records
The Discrete Equation Method (DEM) for Fully Compressible, Two-Phase Flows in Ducts of Spatially Varying Cross-Section
Simple and efficient relaxation methods for interfaces separating compressible fluids, cavitating flows and shocks in multiphase mixtures