skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Multicomponent Reynolds-averaged Navier-Stokes simulations of reshocked Richtmyer-Meshkov instability-induced mixing

Journal Article · · High Energy Density Physics
 [1];  [2]
  1. Univ. of Michigan, Ann Arbor, MI (United States)
  2. Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)

A third-order weighted essentially nonoscillatory (WENO) finite-difference implementation of a two-equation K–ε multicomponent Reynolds-averaged Navier–Stokes (RANS) model is used to simulate reshocked Richtmyer–Meshkov turbulent mixing of air and sulfur hexafluoride at incident shock Mach numbers Mas = 1.24, 1.50, 1.98 with Atwood number At = 0.67 and Mas = 1.45 with At = –0.67. The predicted mixing layer width evolutions are compared with experimental measurements of the width before and after reshock [M. Vetter, B. Sturtevant, Shock Waves 4 (1995) 247; F. Poggi, M.H. Thorembey, G. Rodriguez, Phys. Fluids 10 (1998) 2698] and with the analytical self-similar power-law solution of the simplified model equations before reshock. A new procedure is introduced for the specification of the initial turbulent kinetic energy and its dissipation rate, in which these quantities are related by the linear instability growth rate. The predicted mixing layer widths before reshock are shown to be sensitive to changes in the initial turbulent kinetic energy and its dissipation rate, while the widths after reshock are sensitive to changes in the model coefficients Cε0 and σρ appearing in the buoyancy (shock) production terms in the turbulent kinetic energy and dissipation rate equations. Here, a set of model coefficients and initial conditions is shown to predict mixing layer widths in generally good agreement with the pre-reshock experimental data, and very good agreement with the post-reshock data for all cases. Budgets of the turbulent kinetic energy equation just before and after reshock for the Mas = 1.24 case are used to identify the principal physical mechanisms generating turbulence in reshocked Richtmyer–Meshkov instability: buoyancy production (pressure work) and shear production. Numerical convergence of the mixing layer widths under spatial grid refinement is also demonstrated for each of the Mach numbers considered.

Research Organization:
Lawrence Livermore National Laboratory (LLNL), Livermore, CA (United States)
Sponsoring Organization:
USDOE National Nuclear Security Administration (NNSA)
Grant/Contract Number:
AC52-07NA27344; FC52-08NA28616
OSTI ID:
1788339
Report Number(s):
LLNL-JRNL-560761; 622132; TRN: US2210585
Journal Information:
High Energy Density Physics, Vol. 9, Issue 1; ISSN 1574-1818
Publisher:
ElsevierCopyright Statement
Country of Publication:
United States
Language:
English

References (19)

High-resolution simulations and modeling of reshocked single-mode Richtmyer-Meshkov instability: Comparison to experimental data and to amplitude growth model predictions journal February 2007
Two-dimensional Navier–Stokes simulations of gaseous mixtures induced by Richtmyer–Meshkov instability journal July 2000
A k ‐ε model for turbulent mixing in shock‐tube flows induced by Rayleigh–Taylor instability journal September 1990
Development and validation of a turbulent-mix model for variable-density and compressible flows journal October 2010
High Order Weighted Essentially Nonoscillatory Schemes for Convection Dominated Problems journal February 2009
Physics of reshock and mixing in single-mode Richtmyer-Meshkov instability journal August 2007
PLIF flow visualization and measurements of the Richtmyer–Meshkov instability of an air/SF 6 interface journal August 2002
A second-order turbulence model for gaseous mixtures induced by Richtmyer—Meshkov instability journal January 2005
High-order WENO simulations of three-dimensional reshocked Richtmyer–Meshkov instability to late times: dynamics, dependence on initial conditions, and comparisons to experimental data journal March 2010
The k-L turbulence model for describing buoyancy-driven fluid instabilities journal September 2006
Velocity measurements in turbulent gaseous mixtures induced by Richtmyer–Meshkov instability journal November 1998
Direct numerical simulation of canonical shock/turbulence interaction journal December 2009
Weighted Essentially Non-oscillatory Schemes journal November 1994
Effects of WENO flux reconstruction order and spatial resolution on reshocked two-dimensional Richtmyer–Meshkov instability journal February 2007
K-L turbulence model for the self-similar growth of the Rayleigh-Taylor and Richtmyer-Meshkov instabilities journal August 2006
Large-eddy simulation and multiscale modelling of a Richtmyer–Meshkov instability with reshock journal June 2006
Simulations of Richtmyer–Meshkov instabilities in planar shock-tube experiments journal March 2011
Experiments on the Richtmyer-Meshkov instability of an air/SF6 interface journal March 1995
Instability of Taylor-Sedov blast waves propagating through a uniform gas journal May 1991

Cited By (6)

Two-length-scale turbulence model for self-similar buoyancy-, shock-, and shear-driven mixing journal January 2018
Late-time growth rate, mixing, and anisotropy in the multimode narrowband Richtmyer–Meshkov instability: The θ -group collaboration journal October 2017
Modeling of Rayleigh-Taylor mixing using single-fluid models journal January 2019
Turbulent transport and mixing in the multimode narrowband Richtmyer-Meshkov instability journal September 2019
Turbulent mixing and transition criteria of flows induced by hydrodynamic instabilities journal August 2019
The transition to turbulence in shock-driven mixing: effects of Mach number and initial conditions journal May 2019