Numerical Investigation of Atwood Number Effects on Shock-Driven Single-Mode Stratified Heavy Fluid Layers

This work presents a numerical investigation of Richtmyer–Meshkov instability (RMI) in shock-driven single-mode stratified heavy fluid layers, with emphasis on the influence of the Atwood number. High-order modal discontinuous Galerkin simulations are carried out for Atwood numbers ranging from A=0.30 to 0.72, allowing a systematic study of interface evolution, vorticity dynamics, and mixing. The analysis considers diagnostic quantities such as interface trajectories, normalized interface length and amplitude, vorticity extrema, circulation, enstrophy, and kinetic energy. The results demonstrate that the Atwood number plays a central role in instability development. At low A, interface deformation remains smooth and coherent, with weaker vorticity deposition and delayed nonlinear roll-up. As A increases, baroclinic torque intensifies, leading to rapid perturbation growth, stronger vortex roll-ups, and earlier onset of secondary instabilities such as Kelvin–Helmholtz vortices. Enstrophy, circulation, and interface measures show systematic amplification with increasing density contrast, while the total kinetic energy exhibits relatively weak sensitivity to A. Overall, the study highlights how the Atwood number governs the transition from linear to nonlinear dynamics, controlling both large-scale interface morphology and the formation of small-scale vortical structures. These findings provide physical insight into shock–interface interactions and contribute to predictive modeling of instability-driven mixing in multicomponent flows.