Multicomponent Reynolds-averaged Navier-Stokes simulations of reshocked Richtmyer-Meshkov instability-induced mixing
Journal Article
·
· High Energy Density Physics
- Univ. of Michigan, Ann Arbor, MI (United States)
- 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
- Journal Information:
- High Energy Density Physics, Journal Name: High Energy Density Physics Journal Issue: 1 Vol. 9; ISSN 1574-1818
- Publisher:
- ElsevierCopyright Statement
- Country of Publication:
- United States
- Language:
- English
Similar Records
Multi-component Reynolds-averaged Navier–Stokes simulations of Richtmyer–Meshkov instability and mixing induced by reshock at different times
High-Order Weno Simulations of Three-Dimensional Reshocked Richtmyer-Meshkov Instability to Late Times: Dynamics, Dependence on Initial Conditions, and Comparisons to Experimental Data
Richtmyer-Meshkov instability-induced mixing: initial conditions modeling, three-dimensional simulations and comparisons to experiment
Journal Article
·
Sat Nov 16 19:00:00 EST 2013
· Shock Waves
·
OSTI ID:1788353
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 Article
·
Mon Jan 11 23:00:00 EST 2010
· Acta Mathematica Scientia, vol. 30B, no. 2, April 5, 2010, pp. 595-620
·
OSTI ID:978398
Richtmyer-Meshkov instability-induced mixing: initial conditions modeling, three-dimensional simulations and comparisons to experiment
Conference
·
Sun Jan 07 23:00:00 EST 2007
·
OSTI ID:899437