A compressible Reynolds-averaged mixing model considering turbulent composition and heat fluxes

In typical nature and engineering scenarios, such as supernova explosion and inertial confinement fusion, mixing flows induced by hydrodynamic interfacial instabilities are essentially compressible. Despite their significance, accurate predictive tools for these compressible flows remain scarce. For engineering applications, the Reynolds-averaged Navier–Stokes (RANS) simulation stands out as the most practical approach due to its outstanding computational efficiency. However, existing RANS studies focus primarily on cases where the compressible effect plays an insignificant role in mixing development, with quite limited attention given to scenarios with significant compressibility influence. Moreover, most of the existing RANS mixing models demonstrate significantly inaccurate predictions for the latter. This study develops a novel compressible RANS mixing model by incorporating physical compressibility corrections into the $K$ – $L$ – $\gamma$ mixing transition model recently proposed by Xie et al. ( J. Fluid Mech. 1002, 2025, A31). Specifically, taking the density-stratified Rayleigh–Taylor mixing flows as representative compressible cases, we first analyse the limitations of the existing model for compressible flows, based on high-fidelity data and local instability criteria. Subsequently, the equation of state for a perfect gas is employed to derive comprehensive compressibility corrections. The crucial turbulent composition and heat fluxes are integrated into the closure of the key turbulent mass flux term of the turbulent kinetic energy equation. These corrections enable the model to accurately depict compressible mixing flows. Systematic validations confirm the efficacy of the proposed modelling scheme. This study offers a promising strategy for modelling compressible mixing flows, paving the way for more accurate predictions in complex scenarios.