Molecular dynamics simulation of cylindrical Richtmyer-Meshkov instability

Abstract

The microscopic-scale Richtmyer-Meshkov (RM) instability of a single-mode Cu-He interface subjected to a cylindrically converging shock is studied through the classical molecular dynamics simulation. An unperturbed interface is first considered to examine the flow features in the convergent geometry, and notable distortions at the circular inhomogeneity are observed due to the atomic fluctuation. Detailed processes of the shock propagation and interface deformation for the single-mode interface impacted by a converging shock are clearly captured. Different from the macroscopic-scale situation, the intense molecular thermal motions in the present microscale flow introduce massive small wavelength perturbations at the single-mode interface, which later significantly impede the formation of the roll-up structure. Influences of the initial conditions including the initial amplitude, wave number and density ratio on the instability growth are carefully analyzed. It is found that the late-stage instability development for interfaces with a large perturbation does not depend on its initial amplitude any more. Surprisingly, as the wave number increases from 8 to 12, the growth rate after the reshock drops gradually. The distinct behaviors induced by the amplitude and wave number increments indicate that the present microscopic RM instability cannot be simply characterized by the amplitude over wavelength ratio ($\eta$). The pressure history at the convergence center shows that the first pressure peak caused by the shock focusing is insensitive to $\eta$, while the second one depends heavily on it.