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

Journal Article · Sun Mar 01 00:00:00 EST 2009 · Journal of Computational Physics

DOI:https://doi.org/10.1016/j.jcp.2008.11.002· OSTI ID:949777

Saurel, Richard; Petitpas, Fabien; Berry, Ray A

Numerical approximation of the five-equation two-phase flow of Kapila et al. [A.K. Kapila, R. Menikoff, J.B. Bdzil, S.F. Son, D.S. Stewart, Two-phase modeling of deflagration-to-detonation transition in granular materials: reduced equations, Physics of Fluids 13(10) (2001) 3002–3024] is examined. This model has shown excellent capabilities for the numerical resolution of interfaces separating compressible fluids as well as wave propagation in compressible mixtures [A. Murrone, H. Guillard, A five equation reduced model for compressible two phase flow problems, Journal of Computational Physics 202(2) (2005) 664–698; R. Abgrall, V. Perrier, Asymptotic expansion of a multiscale numerical scheme for compressible multiphase flows, SIAM Journal of Multiscale and Modeling and Simulation (5) (2006) 84–115; F. Petitpas, E. Franquet, R. Saurel, O. Le Metayer, A relaxation-projection method for compressible flows. Part II. The artificial heat exchange for multiphase shocks, Journal of Computational Physics 225(2) (2007) 2214–2248]. However, its numerical approximation poses some serious difficulties. Among them, the non-monotonic behavior of the sound speed causes inaccuracies in wave’s transmission across interfaces. Moreover, volume fraction variation across acoustic waves results in difficulties for the Riemann problem resolution, and in particular for the derivation of approximate solvers. Volume fraction positivity in the presence of shocks or strong expansion waves is another issue resulting in lack of robustness. To circumvent these difficulties, the pressure equilibrium assumption is relaxed and a pressure non-equilibrium model is developed. It results in a single velocity, non-conservative hyperbolic model with two energy equations involving relaxation terms. It fulfills the equation of state and energy conservation on both sides of interfaces and guarantees correct transmission of shocks across them. This formulation considerably simplifies numerical resolution. Following a strategy developed previously for another flow model [R. Saurel, R. Abgrall, A multiphase Godunov method for multifluid and multiphase flows, Journal of Computational Physics 150 (1999) 425–467], the hyperbolic part is first solved without relaxation terms with a simple, fast and robust algorithm, valid for unstructured meshes. Second, stiff relaxation terms are solved with a Newton method that also guarantees positivity and robustness. The algorithm and model are compared to exact solutions of the Euler equations as well as solutions of the five-equation model under extreme flow conditions, for interface computation and cavitating flows involving dynamics appearance of interfaces. In order to deal with correct dynamic of shock waves propagating through multiphase mixtures, the artificial heat exchange method of Petitpas et al. [F. Petitpas, E. Franquet, R. Saurel, O. Le Metayer, A relaxation-projection method for compressible flows. Part II. The artificial heat exchange for multiphase shocks, Journal of Computational Physics 225(2) (2007) 2214–2248] is adapted to the present formulation.

• Research Organization: Idaho National Lab. (INL), Idaho Falls, ID (United States)
• Sponsoring Organization: USDOE
• DOE Contract Number: DE-AC07-99ID-13727
• OSTI ID: 949777
• Report Number(s): INL/JOU-07-12966; JCTPAH; TRN: US0901787
• Journal Information: Journal of Computational Physics, Vol. 228, Issue 5; ISSN 0021-9991
• Country of Publication: United States
• Language: English

Related Subjects

99 GENERAL AND MISCELLANEOUS
45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEFENSE
22 GENERAL STUDIES OF NUCLEAR REACTORS
33 ADVANCED PROPULSION SYSTEMS
ACOUSTICS
ALGORITHMS
APPROXIMATIONS
COMPRESSIBLE FLOW
ENERGY CONSERVATION
EXACT SOLUTIONS
FLOW MODELS
GRANULAR MATERIALS
MIXTURES
MULTIPHASE FLOW
NEWTON METHOD
PHYSICS
RELAXATION
RESOLUTION
SHOCK WAVES
TWO-PHASE FLOW
VELOCITY
WAVE PROPAGATION
cavitation
Godunov
hyperbolic system
multifluid
multiphase
multiphysics
real gases