Bridging the local and the global: a physically constrained buoyancy--drag model for unified prediction of Rayleigh--Taylor and Richtmyer--Meshkov mixing widths across density ratios
Abstract
Accurate prediction of the macroscopic width of Rayleigh--Taylor (RT) and Richtmyer--Meshkov (RM) turbulent mixing layers is central to inertial confinement fusion and supernova dynamics. However, bubble--spike asymmetry, density-ratio dependence and unsteady forcing pose a persistent closure challenge: existing low-order buoyancy--drag models struggle to describe different mixing problems accurately with one model and coefficient set. We combine local front dynamics with global mass conservation in separate buoyancy--drag equations for the bubble and spike fronts. Rather than imposing shared or fixed empirical coefficients, the model retains separate inertia, buoyancy and drag coefficients on the two sides and allows them to vary independently with density ratio. Given the bubble-side state scalings, RT/RM similarity relations, a mean-composition profile and endpoint asymptotics jointly constrain all six effective coefficients without case-by-case fitting. A profile-shape parameter c labels distinct internal composition states and is selected a priori from RT spike scaling measured in linear-electric-motor experiments. The model then cross-predicts the RM spike exponent without recalibration to RM spike data and, by construction, recovers low-Atwood-number bubble--spike symmetry and the high-density-ratio free-fall RT-spike and ballistic RM-spike limits. Tests against constant- and variable-acceleration RT mixing, post-impulse RM evolution and Nova laser deceleration show that one closure describes mixing-width evolution across density ratios and acceleration histories without case-specific retuning, while reducing excessive spike growth at high density ratio. This physically interpretable, asymptotically consistent framework enables cross-problem prediction of wide-density-ratio RT and post-impulse RM mixing.
Related
- Early-time dynamics of the Richtmyer-Meshkov instability at Kelvin-Voigt viscoelastic interfaces
- The Richtmyer-Meshkov Instability of Thermal, Isotope, and Species Interfaces in a five-moment multi-fluid plasma
- Physics-Constrained Neural Surrogate for Domain Growth Prediction in Systems with Conserved Kinetics
- Dynamic control of carbon-ion radiation dose distribution amid respiratory motion challenges: A review of prediction research using geant4 and deep learning physics-constrained models.
- Differential diffusion effects on buoyancy-driven instabilities of acid–base fronts: the case of a color indicator
BibTeX
@misc{zhang2026bridging,
title = {{Bridging the local and the global: a physically constrained buoyancy--drag model for unified prediction of Rayleigh--Taylor and Richtmyer--Meshkov mixing widths across density ratios}},
author = {You-Sheng Zhang and Ya-Feng Li and Meng-Juan Xiao and Yu-Hui Wang},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.02128},
eprint = {2609.02128},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.02128},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Bridging the local and the global: a physically constrained buoyancy--drag model for unified prediction of Rayleigh--Taylor and Richtmyer--Meshkov mixing widths across density ratios
AU - You-Sheng Zhang
AU - Ya-Feng Li
AU - Meng-Juan Xiao
AU - Yu-Hui Wang
PY - 2026
DA - 2026/09/02
PB - arXiv
DO - 10.48550/arxiv.2609.02128
UR - https://arxiv.org/abs/2609.02128
AB - Accurate prediction of the macroscopic width of Rayleigh--Taylor (RT) and Richtmyer--Meshkov (RM) turbulent mixing layers is central to inertial confinement fusion and supernova dynamics. However, bubble--spike asymmetry, density-ratio dependence and unsteady forcing pose a persistent closure challenge: existing low-order buoyancy--drag models struggle to describe different mixing problems accurately with one model and coefficient set. We combine local front dynamics with global mass conservation in separate buoyancy--drag equations for the bubble and spike fronts. Rather than imposing shared or fixed empirical coefficients, the model retains separate inertia, buoyancy and drag coefficients on the two sides and allows them to vary independently with density ratio. Given the bubble-side state scalings, RT/RM similarity relations, a mean-composition profile and endpoint asymptotics jointly constrain all six effective coefficients without case-by-case fitting. A profile-shape parameter c labels distinct internal composition states and is selected a priori from RT spike scaling measured in linear-electric-motor experiments. The model then cross-predicts the RM spike exponent without recalibration to RM spike data and, by construction, recovers low-Atwood-number bubble--spike symmetry and the high-density-ratio free-fall RT-spike and ballistic RM-spike limits. Tests against constant- and variable-acceleration RT mixing, post-impulse RM evolution and Nova laser deceleration show that one closure describes mixing-width evolution across density ratios and acceleration histories without case-specific retuning, while reducing excessive spike growth at high density ratio. This physically interpretable, asymptotically consistent framework enables cross-problem prediction of wide-density-ratio RT and post-impulse RM mixing.
ER -