Hostname: page-component-78c5997874-t5tsf Total loading time: 0 Render date: 2024-11-18T16:36:48.832Z Has data issue: false hasContentIssue false

Amplitude growth rate of a Richtmyer–Meshkov unstable two-dimensional interface to intermediate times

Published online by Cambridge University Press:  17 February 2003

NORMAN J. ZABUSKY
Affiliation:
Department of Mechanical and Aerospace Engineering and CAIP Center, Rutgers University, 98 Brett Road, Piscataway, NJ 08854, USA
ALEXEI D. KOTELNIKOV
Affiliation:
Department of Mechanical and Aerospace Engineering and CAIP Center, Rutgers University, 98 Brett Road, Piscataway, NJ 08854, USA
YURIY GULAK
Affiliation:
Department of Mechanical and Aerospace Engineering and CAIP Center, Rutgers University, 98 Brett Road, Piscataway, NJ 08854, USA
GAOZHU PENG
Affiliation:
Department of Mechanical and Aerospace Engineering and CAIP Center, Rutgers University, 98 Brett Road, Piscataway, NJ 08854, USA

Abstract

The Richtmyer–Meshkov instability in an incompressible and compressible stratified two-dimensional ideal flow is studied analytically and numerically. For the incompressible problem, we initialize a single small-amplitude sinusoidal perturbation of wavelength λ, we compute a series expansion for the amplitude a in powers of t up to t(11) with the MuPAD computer algebra environment. This involves harmonics up to eleven. The simulations are performed with two codes: incompressible, a vortex-in- cell numerical technique which tracks a single discontinuous density interface; and compressible, PPM for a shock-accelerated case with a finite interfacial transition layer (ITL). We identify properties of the interface at time t = tM at which it first becomes ‘multivalued’. Here, we find the normalized width of the ‘spike’ is related to the Atwood number by (wm/λ)−0.5 = −0.33A. A high-order Pad approximation is applied to the analytical series during early time and gives excellent results for the interface growth rate a˙. However, at intermediate times, t > tM, the agreement between numerical results and different-order Padé approximants depends on the Atwood number. During this phase, our numerical solutions give a˙∝O(t−1) for small A and a˙∝O(t−0.4) for A = 0.9. Experimental data of Prasad et al. (2000) for SF6 (post shock Atwood number = 0.74) shows an exponent between −0.68 and −0.72 and we obtain −0.683 for the compressible simulation. For this case, we illustrate the important growth of vortex-accelerated (secondary) circulation deposition of both signs of vorticity and the complex nature of the roll-up region.

Type
Research Article
Copyright
© 2003 Cambridge University Press

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)
Supplementary material: PDF

ZABUSKY et al. supplementary material

Appendix

Download ZABUSKY et al. supplementary material(PDF)
PDF 698.2 KB