Turbulent Richtmyer-Meshkov instability (RMI) is investigated through a series of high resolution three dimensional smulations 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 regime 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 measurement. The asymptotic integral molecular mixing measures Θ = 0.792 ± 0.014, Ξ = 0.800±0.014 and Ψ = 0.782±0.013 which 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 θ. Quantitative data given for both ensemble averages and individual algorithms provide useful benchmark results for future research.

International audienceno abstrac

In this work, we study the large-scale structure of homogeneous flows displaying high density contrasts and small turbulent Mach numbers. Following Batchelor & Proudman [1], we draw an analogy between this configuration and that of a single isolated eddy displaying density non-uniformities. By doing so, we are able to highlight the crucial role played by the solenoidal component of the momentum in the preservation of initial conditions. In particular, we show that the large-scale initial conditions of the spectrum of the solenoidal momentum are preserved provided its infrared exponent is smaller than 4. This condition is reminiscent of the one usually derived in constant density flow, except that it does not apply to the velocity spectrum. The latter is actually shown to be impermanent for infrared exponents larger or equal to 2. The consequences of these properties on the self-similar decay of the flow are discussed. Finally, these predictions are verified by performing large eddy simulations of homogeneous isotropic turbulence.