---
title: Compressible Rayleigh–Taylor Turbulence
url: https://www.emergentmind.com/topics/compressible-rayleigh-taylor-turbulence
type: topic
---

# Compressible Rayleigh–Taylor Turbulence

Compressible Rayleigh–Taylor (RT) turbulence encompasses the nonlinear evolution and turbulent mixing that arises when a heavier fluid overlies a lighter one in the presence of acceleration, accounting for compressibility effects that manifest prominently at finite Atwood numbers, stratification, finite Mach number, and in thermodynamically non-equilibrium regimes. In contrast to the classical Boussinesq RT instability, compressible RT turbulence involves density and pressure variations, baroclinic drivers, and temperature and entropy fluctuations—core features relevant in astrophysical, fusion, and geophysical flows.

## 1. Governing Equations and Physical Regimes

Compressible RT turbulence is governed by the compressible Navier–Stokes–Fourier system, typically expressed as:
- Mass continuity: $\partial_t \rho + \nabla \cdot (\rho \mathbf{u}) = 0$
- Momentum: $\partial_t(\rho \mathbf{u}) + \nabla \cdot (\rho \mathbf{u} \otimes \mathbf{u} + p \mathbf{I} - \boldsymbol{\Pi}) = \rho \mathbf{a}(t)$
- Total energy: $\partial_t E + \nabla \cdot [(E+p)\mathbf{u} - \mathbf{u} \cdot \boldsymbol{\Pi} - \mathbf{q}] = \rho \mathbf{a}(t) \cdot \mathbf{u}$

Here, $p=\rho T$, $E=\frac12\rho |\mathbf{u}|^2 + \rho e$, $\boldsymbol{\Pi}$ is the viscous stress, and $\mathbf{q}$ the heat flux. The acceleration $\mathbf{a}(t)$ can be time-dependent to model variable forcing scenarios. Compressibility is parameterized by the Atwood number $A=(\rho_H-\rho_L)/(\rho_H+\rho_L)$ and Mach number $Ma=U/c_s$.

Key non-dimensional numbers include:
- Reynolds number $Re=UL/\nu$
- Mach number $Ma=U/\sqrt{T}$
- Rayleigh number $Ra=g\beta L^4/(\nu\kappa)$
- Stratification parameter and adiabatic gradient $\beta_{ad}=g/c_p$, leading to an adiabatic length $L_{ad}=c_p\Delta T/g$ [1005.3639, 1009.5819]

Discrete Boltzmann and lattice Boltzmann modeling frameworks, including D2Q16, D2Q37, and higher-order Hermite expansions, recover these macroscopic equations with correct thermodynamic couplings for compressible systems [1005.3639, 2107.03658, 1009.5819, 2504.05128].

## 2. Non-Equilibrium Thermodynamics and Scaling Measures

Compressible RT turbulence is characterized by persistent departures from local thermodynamic equilibrium, especially in regions of strong density and temperature gradients—typically at the bubble and spike tips. Thermodynamic non-equilibrium (TNE) is quantified via kinetic moment deviations:
- Viscous-stress TNE: $\Delta_{2,\alpha\beta}^* = m^\sigma [M_{2,\alpha\beta}^*(f_i^\sigma) - M_{2,\alpha\beta}^*(f_i^{\sigma,eq})]$
- Total TNE: $\Delta^* = |\Delta_2^*| + |\Delta_{3,1}^*| + |\Delta_3^*| + |\Delta_{4,2}^*|$

Domain-averaged quantities such as $\bar D^\sigma$ and fraction of non-equilibrium regions $S_r^\sigma$ provide mesoscopic metrics for turbulence intensity and mixing state [2504.05128, 2107.03658].

Typical temporal patterns for $\bar D^\sigma$ and $|\nabla\rho|$ include an initial dip (due to diffusion), a peak during nonlinear growth (interface stretching, vortex formation), and decay at late times as dissipation dominates. Exponential and linear laws parameterize the dependence of gradients and TNE on acceleration period, amplitude, and phase—allowing predictive control over mixing and non-equilibrium intensity [2504.05128].

## 3. Mixing-Layer Dynamics, Growth Laws, and Compressibility Effects

In classical (incompressible) RT, the mixing-layer width $L_{ml}(t)$ grows quadratically:
$$
\left(\frac{dL}{dt}\right)^2 = 4\alpha A g L,\quad L(t) \sim \alpha A g t^2
$$
with typical $\alpha$ in the $0.03$–$0.05$ range [1005.3639, 1009.5819]. Compressibility introduces key modifications:
- Up–down asymmetry emerges: spikes (denser/cold plumes) grow faster than bubbles as Atwood number increases, with $\alpha_s/\alpha_b$ shifting from $\sim1.05$ at $A\sim0.05$ to $>1.5$ at $A\sim0.4$ [1005.3639, 1009.5819].
- Finite At leads to $O(\Delta\rho/\rho)$ pressure fluctuations, slower homogenization of rms velocities and temperatures, and persistent intermittency in high-order statistics.
- At high stratification, mixing arrests at the finite “adiabatic length” $L_{ad}$, with the mean temperature gradient approaching the adiabatic value and buoyancy vanishing. This is modeled with phenomenological “shutoff” functions $\psi(L/L_{ad})$ modifying the quadratic growth law [1009.5819].

For time-varying acceleration, mixing and TNE measures oscillate in-phase with the forcing, and the timing and magnitude of mixing events can be systematically altered by tuning forcing waveform parameters [2504.05128].

## 4. Numerical and Modeling Methodologies

Several advanced techniques have been developed for compressible RT turbulence:
- **Lattice Boltzmann Methods (LBM):** Accurately recover thermodynamically consistent compressible dynamics, handle large Atwood number and stratification, and capture growth rates, arrest effects, and turbulence statistics [1005.3639, 1009.5819].
- **Kinetic and Discrete Boltzmann Models:** Resolve TNE effects and allow for rigorous quantification of mesoscopic non-equilibrium, including the effects of specific-heat ratio $\gamma$, which amplifies TNE and accelerates vortex formation as $\gamma$ decreases [2107.03658, 2504.05128].
- **Data-Driven Large-Eddy Simulation (LES):** Fourier Neural Operator (FNO) models, trained on filtered DNS data, outperform classical LES subgrid closures (VGM, DSM, ILES) in predicting kinetic energy, mixing-width evolution, concentration and temperature fluctuations, and velocity divergence for 3D compressible RT at $A_t=0.5$, $Sr=1$, $Re=10^4$–$3\times10^4$ [2404.05834]. FNO achieves $O(1\%)$ accuracy on bubble/spike mixing-widths and delivers order-of-magnitude computational speedups over traditional methods.
- **Multiscale Decomposition:** Separation into a compressible bulk (coarse grid) and incompressible–irrotational interface dynamics (“higher-order z-model”) yields speedups of $10^2$–$10^3\times$ relative to fully resolved compressible Euler simulations, with close agreement on classical RT/RM benchmarks [1904.04935].

## 5. Influence of System Parameters and Forcing

The dynamics of compressible RT turbulence are highly sensitive to:
- **Atwood Number ($A$):** Increases amplify asymmetry, pressure fluctuations, and decrease the rate of mixing-layer homogenization.
- **Specific-Heat Ratio ($\gamma$):** Lower $\gamma$ (higher compressibility) increases TNE intensities, proportion of non-equilibrium regions, and advances the onset of vortical breakdown [2107.03658].
- **Time-Dependent Forcing:** By varying acceleration period ($T_0$), amplitude ($A_0$), and phase ($\Phi$), the onset, duration, and intensity of mixing can be delayed or enhanced, following explicit exponential and linear scaling laws for density gradients and TNE measures [2504.05128].

### Example scaling laws from [2504.05128]:

| Parameter | Density Gradient | TNE Strength ($\bar D^A$, $\bar D^B$) |
|-----------|------------------|--------------------------------------|
| $T_0$ @ $t=0.8$ | $404.60 + 167.75\,e^{-1.79 T_0}$ | $0.085+0.313\,e^{-2.032 T_0}$ / $0.088+0.302\,e^{-1.876 T_0}$ |
| $A_0$ (peak) | $458.96-13.36\,e^{-0.01A_0}$ | $0.183-0.031\,e^{-0.775A_0}$ / $0.185-0.026\,e^{-0.902A_0}$ |
| $\Phi$ @ $t=0.2$ | Piecewise linear | Piecewise linear/quadratic |

These data establish quantitative frameworks for regulating RT mixing in engineering and astrophysical applications.

## 6. Applications, Implications, and Open Directions

Compressible RT turbulence is central to inertial confinement fusion (where controlling the growth of RT-induced mixing is critical), core-collapse supernovae and stellar interiors (where stratified, compressible plumes dominate), and aerospace/atmospheric flows. The ability to modulate RT turbulence via time-dependent acceleration or by tuning fluid properties (e.g., $\gamma$) yields potential for flow control strategies in fusion or laboratory experiments [2504.05128, 1005.3639].

Future research challenges include:
- Detailed three-dimensional high-Reynolds-number studies with variable Prandtl and Mach numbers [1009.5819, 2404.05834].
- Extension of data-driven LES and operator-learning frameworks to more complex baroclinic and multi-material interfaces.
- Comprehensive characterizations of intermittency, small-scale cascade properties, and the role of non-equilibrium in transition and late-time mixing.
- Analytical closure and prediction of adiabatically arrested mixing regimes and their dependence on compressibility, stratification, and drive waveform.

Compressibility in the RT context thus fundamentally alters turbulence development, saturation, and energy partition, mandating kinetic, mesoscale, and data-driven approaches for predictive accuracy and computational efficiency.

Source: https://www.emergentmind.com/topics/compressible-rayleigh-taylor-turbulence