Richtmyer-Meshkov Instability (RMI) in a spherical geometry is studied via direct numerical simulation using a high-order three-dimensional in-house solver. Specifically, a six-order compact difference scheme coupled with localized artificial diffusivity method is adopted in order to capture discontinuities with high accuracy. A pure converging shock propagation in a sphere is simulated and the result agrees well with Guderley's theory. For RMI in a spherical geometry, the development of mixing width and its growth rate at different stages are examined and the underlying mechanism is also briefly analyzed. Particularly addressed is the effect of Mach number on the growth rate of perturbations and turbulent mixing process.
Bibliographical noteFunding Information:
We would like to thank Xisheng Luo for many valuable discussions on this work. Numerical simulations were carried out on the Tianhe-2 supercomputing facility at National Supercomputer Center in Guangzhou, China. We acknowledge the financial supports provided by National Natural Science Foundation of China (Grants No. U1630138 and No. 91852112). This work was also supported by the Challenge Program (Grant No. JCKY2016212A501).
© 2019 Global Science Press.
- Direct numerical simulation
- Mach number
- Richtmyer-Meshkov instability
- Spherical geometry