A theoretical one-dimensional model for variable-density Rayleigh-Taylor turbulence

TL;DR

This paper revisits and extends a 1965 theoretical model for variable-density Rayleigh-Taylor turbulence, analyzing solutions of the similarity equations and comparing them with numerical and experimental data.

Contribution

It provides a detailed analysis of the full and simplified similarity equations, demonstrating their ability to capture key RT turbulence features and validating the model against empirical results.

Findings

Full similarity equation captures asymmetric spike and bubble growth.

Comparison shows reasonable agreement with numerical and experimental data.

Global mass correction approximates the full solution effectively.

Abstract

In an early theoretical work published in 1965, Belen'kii & Fradkin proposed a turbulent diffusivity model for Rayleigh--Taylor (RT) mixing. We review its derivation and present alternative arguments leading to the same final similarity equation. The original work then introduced an approximation that led to a simplified ordinary differential equation (ODE), which was used primarily to derive the important scaling result, $h \sim (\ln R)gt^2$. Here, we extend the analysis by examining the solutions to both the full similarity ODE and the simplified ODE in detail. It is shown that the full similarity equation captures many now well-known features of non-Boussinesq RT flows, including asymmetric spike and bubble growth and a systematic shift of velocity statistics toward the light-fluid side. Comparisons of the theoretical model with numerical and experimental studies show reasonable agreement in both spatial profiles and growth trends of mixing layer heights. We further show that a global mass correction applied to the simplified solution closely approximates the full solution, highlighting that, to leading order, RT mixing is governed by the competing dynamics between diffusion of $\ln \bar{\rho}$ and mass conservation.