Late-time growth rate, mixing, and anisotropy in the multimode narrowband Richtmyer–Meshkov instability: The θ-group collaboration
- Univ. of Sydney, NSW (Australia). School of Aerospace, Mechanical and Mechatronic Engineering
- Alternative Energies and Atomic Energy Commission (CEA), Arpajon (France)
- Univ. of North Carolina, Charlotte, NC (United States). Mechanical Engineering and Engineering Science
- Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
- AWE, Aldermaston (United Kingdom)
- Russian Academy of Sciences (RAS), Moscow (Russian Federation). Keldysh Inst. of Applied Mathematics; Russian Academy of Sciences (RAS), Moscow (Russian Federation). P.N. Lebedev Physical Inst.
- Russian Academy of Sciences (RAS), Moscow (Russian Federation). Keldysh Inst. of Applied Mathematics
- Russian Academy of Sciences (RAS), Moscow (Russian Federation). P.N. Lebedev Physical Inst.
- Univ. of Strathclyde, Glasgow (Scotland). Dept. of Mechanical and Aerospace Engineering
Turbulent Richtmyer–Meshkov instability (RMI) is investigated through a series of high resolution three-dimensional simulations of two initial conditions with eight independent codes. The simulations are initialised with a narrowband perturbation such that instability growth is due to non-linear coupling/backscatter from the energetic modes, thus generating the lowest expected growth rate from a pure RMI. By independently assessing the results from each algorithm and computing ensemble averages of multiple algorithms, the results allow a quantification of key flow properties as well as the uncertainty due to differing numerical approaches. A new analytical model predicting the initial layer growth for a multimode narrowband perturbation is presented, along with two models for the linear and non-linear regimes combined. Overall, the growth rate exponent is determined as θ=0.292±0.009, in good agreement with prior studies; however, the exponent is decaying slowly in time. Also, θ is shown to be relatively insensitive to the choice of mixing layer width measurements. The asymptotic integral molecular mixing measures Θ=0.792±0.014, Ξ=0.800±0.014, and Ψ=0.782±0.013 are lower than some experimental measurements but within the range of prior numerical studies. The flow field is shown to be persistently anisotropic for all algorithms, at the latest time having between 49% and 66% higher kinetic energy in the shock parallel direction compared to perpendicular and does not show any return to isotropy. The plane averaged volume fraction profiles at different time instants collapse reasonably well when scaled by the integral width, implying that the layer can be described by a single length scale and thus a single θ. In conclusion, quantitative data given for both ensemble averages and individual algorithms provide useful benchmark results for future research.
- Research Organization:
- Lawrence Livermore National Laboratory (LLNL), Livermore, CA (United States)
- Sponsoring Organization:
- USDOE National Nuclear Security Administration (NNSA)
- Grant/Contract Number:
- AC52-07NA27344; AC52-07NA2734; AC52-06NA2-5396
- OSTI ID:
- 1468910
- Alternate ID(s):
- OSTI ID: 1402116
- Report Number(s):
- LLNL-JRNL-733858; 885603
- Journal Information:
- Physics of Fluids, Vol. 29, Issue 10; ISSN 1070-6631
- Publisher:
- American Institute of Physics (AIP)Copyright Statement
- Country of Publication:
- United States
- Language:
- English
Similar Records
Simulations and Models for the Richtmyer-Meshkov Instability with Broadband Perturbations
Simulations and models for the Richtmyer– Meshkov instability with broadband perturbations