A comparative study of the turbulent Rayleigh–Taylor instability using high-resolution three-dimensional numerical simulations: The Alpha-Group collaboration

• Published in: Physics of fluids (1994) Vol. 16; no. 5; pp. 1668 - 1693
• Main Authors: Dimonte, Guy, Youngs, D. L., Dimits, A., Weber, S., Marinak, M., Wunsch, S., Garasi, C., Robinson, A., Andrews, M. J., Ramaprabhu, P., Calder, A. C., Fryxell, B., Biello, J., Dursi, L., MacNeice, P., Olson, K., Ricker, P., Rosner, R., Timmes, F., Tufo, H., Young, Y.-N., Zingale, M.
• Format: Journal Article
• Language: English
• Published: 01.05.2004
• ISSN: 1070-6631; 1089-7666
• DOI: 10.1063/1.1688328

The turbulent Rayleigh–Taylor instability is investigated in the limit of strong mode-coupling using a variety of high-resolution, multimode, three dimensional numerical simulations (NS). The perturbations are initialized with only short wavelength modes so that the self-similar evolution (i.e., bubble diameter D b ∝ amplitude  h b ) occurs solely by the nonlinear coupling (merger) of saturated modes. After an initial transient, it is found that h b ∼α b Agt 2 , where A=Atwood number, g=acceleration, and t=time. The NS yield D b ∼h b /3 in agreement with experiment but the simulation value α b ∼0.025±0.003 is smaller than the experimental value α b ∼0.057±0.008. By analyzing the dominant bubbles, it is found that the small value of α b can be attributed to a density dilution due to fine-scale mixing in our NS without interface reconstruction (IR) or an equivalent entrainment in our NS with IR. This may be characteristic of the mode coupling limit studied here and the associated α b may represent a lower bound that is insensitive to the initial amplitude. Larger values of α b can be obtained in the presence of additional long wavelength perturbations and this may be more characteristic of experiments. Here, the simulation data are also analyzed in terms of bubble dynamics, energy balance and the density fluctuation spectra.