On the Richtmyer-Meshkov instability evolving from a deterministic multimode planar interface

Abstract

We investigate the shock-induced turbulent mixing between a light and a heavy gas, where a Richtmyer-Meshkov instability (RMI) is initiated by a shock wave with Mach number Ma = 1.5. The prescribed initial conditions define a deterministic multimode interface perturbation between the gases, which can be imposed exactly for different simulation codes and resolutions to allow for quantitative comparison. Wellresolved large-eddy simulations are performed using two different and independently developed numerical methods with the objective of assessing turbulence structures, prediction uncertainties and convergence behaviour. The two numerical methods differ fundamentally with respect to the employed subgrid-scale regularisation, each representing state-of-the-art approaches to RMI. Unlike previous studies, the focus of the present investigation is to quantify the uncertainties introduced by the numerical method, as there is strong evidence that subgrid-scale regularisation and truncation errors may have a significant effect on the linear and nonlinear stages of the RMI evolution. Fourier diagnostics reveal that the larger energy-containing scales converge rapidly with increasing mesh resolution and thus are in excellent agreement for the two numerical methods. Spectra of gradient-dependent quantities, such as enstrophy and scalar dissipation rate, show stronger dependences on the small-scale flow field structures as a consequence of truncation error effects, which for one numerical method are dominantly dissipative and for the other dominantly dispersive. Additionally, the study reveals details of various stages of RMI, as the flow transitions from large-scale nonlinear entrainment to fully developed turbulent mixing. The growth rates of the mixing zone widths as obtained by the two numerical methods are ~t7/12 before re-shock and ~(t - t0)2/7 long after re-shock. The decay rate of turbulence kinetic energy is consistently ~(t-t0)-10/7 at late times, where the molecular mixing fraction approaches an asymptotic limit ≈ 0.85. The anisotropy measure.