Direct numerical simulation of the multimode narrowband Richtmyer–Meshkov instability. (15th November 2019)

• Record Type: Journal Article
• Authors: Groom, M.; Thornber, B.
• Is Part Of: Computers & fluids. Volume 194(2019)
• Publication Date: 2019-11-15
• Subjects: Shock wave -- Turbulent mixing -- Compressible -- Turbulence -- Multispecies; Fluid dynamics -- Data processing -- Periodicals; 532.050285
• DOI: 10.1016/j.compfluid.2019.104309
• Languages: English
• ISSNs: 0045-7930

Abstract

Highlights: Direct numerical simulations performed of a planar mixing layer induced by Richtmyer–Meshkov instability. Grid convergence study demonstrates that all integral and spectral quantities of interest are fully resolved. Comparisons are made with previous implicit large eddy simulations of the same initial condition in the high Reynolds number limit. Viscous suppression of secondary instabilities reduces transfer of turbulent kinetic energy to transverse components. Reduced energy transfer gives rise to increased levels of anisotropy and greater rate of decay at late time. Abstract: Early to intermediate time behaviour of the planar Richtmyer–Meshkov instability (RMI) is investigated through direct numerical simulation (DNS) of the evolution of a deterministic interfacial perturbation initiated by a Ma = 1.84 shock. The model problem is the well studied initial condition from the recent θ -group collaboration [16]. A grid convergence study demonstrates that the Kolmogorov microscales are resolved by the finest grid for the entire duration of the simulation, and that both integral and spectral quantities of interest are converged. Comparisons are made with implicit large eddy simulation (ILES) results from the θ -group collaboration, generated using the same numerical algorithm. The total amount of turbulent kinetic energy (TKE) is shown to be decreased in the DNS compared to the ILES, particularly in the transverse directions, giving rise to a greater level of anisotropy in the flow (70% vs. 40% more TKE in the shock parallel direction at the latest time considered). This decrease in transfer of TKE to the transverse components is shown to be due to the viscous suppression of secondary instabilities. Overall the agreement in the large scales between the DNS and ILES is very good and hence the mixing width and growth rate exponent θ are very similar. There are substantial differences in the small scale behaviour however, with a 38% difference observed in the minimum values obtained for the mixing fractions Θ and Ξ. Differences in the late time decay of TKE are also observed, with decay rates calculated to be τ − 1.41 and τ − 1.25 for the DNS and ILES respectively.