Pith. sign in

REVIEW 2 major objections 4 minor 50 references

A compressible Reynolds-averaged mixing model considering turbulent entropy and heat flux

T0 review · 2 major / 4 minor · reviewed 2026-08-06 · deepseek-v4-flash

Pith's one-line read Extending an incompressible K-L-gamma mixing-transition RANS model to stratified compressible Rayleigh-Taylor flows, this paper replaces the baroclinic density-gradient source with entropy-gradient and counter-gradient heat-flux terms and…

desk verdict A coherent compressible RANS closure with real derivation value, but the validation is in-sample because C4=C3 is tuned on the same HiFi data used as ground truth. read the letter →

arxiv 2506.16296 v1 pith:JMT5WEMU submitted 2025-06-19 physics.flu-dyn

classification physics.flu-dyn MSC 76F2576F5576N1576E17 PACS 47.20.Ma47.27.-i47.40.-x
keywords compressibleRayleigh-TaylorinstabilityRANSturbulencemodelingK-L-gammamixingtransitionmodelturbulentmassfluxclosureheatcounter-gradientdiffusiondensitystratificationinterfacial
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper extends the incompressible $K$-$L$-$\gamma$ mixing-transition RANS model to density-stratified compressible Rayleigh-Taylor flows, which matter for supernova and inertial-confinement-fusion applications. It identifies the failure of the baseline model in the closure of the turbulent mass flux: the baroclinic density-gradient product that feeds turbulent kinetic energy production switches off when stratification reverses the density gradient, even though the flow remains unstable. Using the perfect-gas equation of state and the Gibbs relation, the authors rewrite the turbulent mass flux as an entropy-gradient term plus a counter-gradient turbulent heat-flux term and fold both into the local Atwood number closure (2.38). The corrected closure keeps all original model coefficients, reduces to the incompressible form as the sound speed becomes large, and matches high-fidelity mixing widths and density and mass-fraction profiles for stratification parameters $S_r = 0.5$ to $3$. If the claim holds, compressible mixing predictions in engineering-scale settings gain a practical RANS route without a new empirical-coefficient calibration campaign.

What carries the argument

The central object is the modified local Atwood number closure (2.38): $A_{ssi} = \frac{\tilde{L}}{\bar{\rho}+\tilde{L}} |\partial\bar{\rho}/\partial x_i| [ C_A\,\partial\bar{\rho}/\partial x_i + \frac{C_A(\bar{\Gamma}-1)}{\bar{\Gamma}} \frac{1}{\bar{c}^2} (\bar{\rho}\bar{c}_p\,\partial\bar{T}/\partial x_i - \frac{2\bar{\Gamma}-1}{\bar{\Gamma}-1}\,\partial\bar{p}/\partial x_i) ]$, with $C_3=C_4=C_A(\bar{\Gamma}-1)/\bar{\Gamma}$. In one algebraic change, the raw density-gradient factor is replaced by the compressible instability combination plus a counter-gradient heat-flux term. This routes the compressibility corrections into the buoyancy production term $S_{Kf}$ without adding model coefficients, and it reduces exactly to the original closure as $\bar{c}\to\infty$.

What would settle it

Measure the turbulent heat flux $T'u'_i$ in the high-fidelity $S_r=3$ case at $t/\tau=4$ and test whether the counter-gradient closure (2.33) with $C_4=C_A(\bar{\Gamma}-1)/\bar{\Gamma}$ reproduces its sign and magnitude in the mixing zone; if the counter-gradient term is mis-signed, the corrected model's agreement with mixing widths would have to come from error cancellation rather than the claimed mechanism. A complementary test is to run the same closure on a compressible Richtmyer-Meshkov case with unchanged coefficients, an extension the paper leaves for future work.

Watch

Extended reading notes

Core claim

On its own terms, the paper claims that the failure of incompressible RANS models on density-stratified compressible Rayleigh-Taylor flows has a specific cause: the buoyancy source in the TKE equation is closed through the baroclinic product $\partial\bar{p}/\partial x_i\,\partial\bar{\rho}/\partial x_i$, and in stratified flows the density gradient reverses sign too early, so the production term switches off while the flow is still unstable. Using the perfect-gas equation of state and the Gibbs relation, the authors derive the turbulent mass flux velocity $u''_i$ as a sum of an entropy-gradient contribution, which mirrors the compressible instability criterion, and a turbulent heat-flux contribution that exhibits counter-gradient transport. Folding these into the local Atwood number produces closure (2.38), with $C_3=C_4=C_A(\bar{\Gamma}-1)/\bar{\Gamma}$ so that no new model coefficients enter. With this closure, the $K$-$L$-$\gamma$ mixing-transition model matches high-fidelity bubble and spike mixing widths and mean density and mass-fraction profiles for stratification parameters $S_r = 0.5, 1, 2, 3$, including the strongly stratified cases where the original baroclinic source vanishes.

Load-bearing premise

The load-bearing premise is the fitted equality $C_4=C_3$: the heat-flux coefficient was calibrated to the same high-fidelity simulations that later serve as validation, so the model's predictive accuracy outside that calibration set remains an assumption.

Editorial extensions

If this is right

  • For stratification parameters $S_r=0.5$, $1$, $2$, and $3$, the corrected model reproduces high-fidelity bubble and spike mixing widths and mean density and mass-fraction profiles, while the baseline model stagnates at late times.
  • The original incompressible coefficients are retained unchanged, so existing calibrated parameter sets transfer directly to the compressible regime.
  • The closure supplies turbulent kinetic energy production in regions where the baroclinic product vanishes, which is the specific failure mode identified in the baseline model.
  • In the incompressible limit $\bar{c}\to\infty$, the modified local Atwood number reduces exactly to the original form, preserving the baseline model's incompressible behavior.
  • The derivation identifies turbulent entropy flux and turbulent heat flux as the two channels through which compressibility enters the dominant production term of the TKE equation.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Because $C_4=C_3$ was inferred by optimal matching to high-fidelity data rather than measured independently, a sharper test would extract $C_4$ directly from turbulent heat-flux statistics in the same or different flows.
  • The entropy-gradient/heat-flux mechanism is likely not specific to Rayleigh-Taylor flow; applying closure (2.38) unchanged to compressible Richtmyer-Meshkov and Kelvin-Helmholtz mixing would test that generalization, which the paper explicitly leaves open.
  • The derivation neglects fluctuations of the adiabatic index, so flows with strongly dissimilar molecular structures or extreme temperatures, such as some inertial-confinement-fusion mixtures, may need a $\Gamma'$ term; this is an untested edge of the parameter-free claim.
  • The counter-gradient heat flux correlates with organized large-scale structures at high $S_r$, which suggests an observable link between the closure's success and flow organization that conditional statistics of high-fidelity data could check.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

2 major / 4 minor

Summary. The paper extends the incompressible K-L-gamma mixing transition model of Xie et al. (2025) to density-stratified compressible Rayleigh-Taylor flows. Using the perfect-gas equation of state and the Gibbs relation, the authors derive a closure for the turbulent mass flux that replaces the density-gradient term in the local Atwood number with an entropy-gradient term and adds a counter-gradient turbulent heat-flux term. The final closure in Eq. (2.38) preserves the original coefficient C_A by setting C3 = C_A(Gamma-1)/Gamma and C4 = C3. The model is validated against the Luo-Wang HiFi simulations for Sr = 0.5, 1, 2, and 3, showing good agreement for mixing widths and mean species/density profiles.

Significance. If the predictive claim survives, this is a useful contribution to RANS modeling of compressible interfacial mixing: the derivation is systematic, the incompressible limit of Eq. (2.35) correctly constrains C3, and the final form (2.38) retains the baseline coefficient structure on paper. The emphasis on counter-gradient turbulent heat flux in stratified compressible RT flows is a physically interesting feature that standard gradient-diffusion closures miss. However, the central evidence is weakened by the in-sample calibration of C4 and by the absence of any out-of-sample test, so the general predictive claim for compressible mixing flows is not yet established.

major comments (2)
  1. [Sections 2.3 and 3, Eqs. (2.33)-(2.38) and Figs. 8-10] The closure coefficient C4 is not fixed by the incompressible-limit constraint; the paper states in Section 2.3 that 'optimal predictions' are obtained when C4=C3, 'based on calibrations with HiFi data'. This makes C4 effectively a fitted parameter, despite the claim in Sections 2.3 and 4 that the model 'maintains the original model coefficient C_A without introducing additional parameters'. Because the same HiFi data are then used as ground truth in all four validation cases in Section 3, the agreement shown in Figs. 8-10 is an in-sample measure of the calibration, not an independent test of the closure. A decisive out-of-sample test, such as a different Atwood number, a non-isothermal initial stratification, or an RM or KH case, is needed before the predictive claim can be sustained.
  2. [Section 3, Figs. 8-10] All four validated cases share the same flow configuration: isothermal density-stratified RT with Atwood number 0.5, Gamma=5/3, and Sr = 0.5, 1, 2, and 3. Since C4=C3 is calibrated against these very cases, the close agreement in mixing widths, mass-fraction profiles, and density profiles is expected and does not demonstrate generalization. The abstract and title claim a model for 'compressible mixing flows' generally, but no RM, KH, or different-Atwood case is shown. The final paragraph of Section 4 acknowledges this limitation, but the absence of out-of-sample validation is a load-bearing gap for the central claim of the paper.
minor comments (4)
  1. [Section 2.3, after Eq. (2.38)] The sentence 'In the incompressible limit c to 0' should read 'c to infinity'.
  2. [Section 2.2] The phrase 'It it attributed' is a typo and should read 'It is attributed'.
  3. [Section 2.3] The term 'contour-gradient diffusion' appears to be a typo for 'counter-gradient diffusion'.
  4. [Eq. (2.31)] The speed of sound c-bar is used before being defined; the definition c-bar^2 = Gamma-bar p-bar / rho-bar should be stated near its first use.

Circularity Check

1 steps flagged · score 6.0 of 10

The C4=C3 relation is fixed by calibration on the same HiFi cases used for validation, so the Section 3 agreement is an in-sample fit check rather than an independent prediction.

  1. fitted input called prediction [Section 2.3, after Eq. (2.37); validated in Section 3, Figs. 8-10]
    "Additionally, it is found that the present model can achieve optimal predictions when C4=C3, based on calibrations with HiFi data. Consequently, C4 = C3 = CA(Γ−1)/Γ is adopted in this study."

    C4 is introduced as a coefficient 'to be determined' in Eq. (2.33), and it is not derived from the instability criterion or the thermodynamics chain; it is set equal to C3 because that choice gave optimal agreement with the HiFi data. The validation in Section 3 then compares the model against those same HiFi cases (Sr = 0.5, 1, 2, 3 from Luo & Wang 2022), with the flow configuration 'remain consistent with the baseline model' as stated in Section 2.2. The agreement in Figs. 8-10 therefore reflects how well the fitted relation C4=C3 reproduces the calibration data, not an out-of-sample test of the closure.

full rationale

Most of the algebraic derivation is self-contained and not circular. The entropy-gradient term in Eq. (2.31) is derived from the perfect-gas EOS, the Gibbs relation, and the compressible instability criterion of Eq. (2.18) with an external reference (Gamalii et al. 1980); the constraint C3 = CA(Γ−1)/Γ follows from requiring the incompressible limit to restore the baseline closure. These steps do not assume the target result. The circularity is localized to the heat-flux coefficient: C4=C3 is explicitly calibrated against HiFi data, and the validation in Figures 8-10 uses the same HiFi dataset as ground truth. Thus the central predictive claim for the four Sr cases is partially an in-sample demonstration. The self-citation to the baseline K-L-gamma model (Xie et al. 2025) is not load-bearing in a circular sense, because the baseline coefficients are tabulated in Table 1 and the baseline is used as a starting point rather than as evidence for the new correction. The paper also honestly defers compressible RM and KH validation to future work, which reinforces that generalizability beyond the calibrated RT cases is not yet established. Overall, the derivation is not vacuous, but the fitted C4=C3 relation and the in-sample validation justify a partial-circularity score of 6.

Assumptions & free parameters 1 free parameters · 5 assumptions · 0 invented entities

The central claim rests on the perfect-gas EOS, the Gibbs relation linearization, the borrowed Cloutman heat-flux closure, the neglect of dynamic compressibility, and the calibrated relation C4=C3. The only genuinely new free choice is C4=C3, calibrated against the same HiFi data used for validation.

free parameters (1)
  • C4 (heat-flux model coefficient, set equal to C3) = C4 = C3 = C_A(Gamma-1)/Gamma, approximately 2.835 for Gamma=5/3
    The equality C4=C3 is not derived; it is chosen after calibrating against HiFi mixing-width data. All four validation cases are in the calibration set, so this is a free parameter of the validation claim.
assumptions (5)
  • domain assumption Perfect-gas equation of state p=(Gamma-1)rho e and mixture EOS with isothermal, partial-pressure assumptions (Livescu 2013).
    Used in Section 2.3 to derive the density fluctuation relation (2.21)-(2.22). Standard for these flows but restricts the result to perfect gases with negligible Gamma-prime fluctuations.
  • standard math Gibbs relation T ds = dh - dp/rho applied to small perturbations between adjacent equilibrium states.
    Used to derive Eq. (2.26); assumes fluctuations small enough for the relation to hold linearly.
  • ad hoc to paper Nonlinear fluctuation correlations and Gamma-prime fluctuations are neglected in the derivation.
    Required for Eqs. (2.27)-(2.31); the mixing layer fluctuations may be finite-amplitude, so the linearization is an unverified assumption.
  • domain assumption Cloutman (2003) closure (2.33) for turbulent heat flux with counter-gradient diffusion.
    Borrowed closure; no independent verification in this paper.
  • domain assumption Static compressibility dominates and dynamic compressibility is neglected.
    Stated in Section 1 following prior literature; restricts scope to density-stratified RT.

how reviews work

0 comments
Cite this review

Pith. "Pith review of A compressible Reynolds-averaged mixing model considering turbulent entropy and heat flux." pith.science (2026). https://pith.science/paper/JMT5WEMU

@misc{pith2026250616296,
  author       = {Pith},
  title        = {Pith review of: A compressible Reynolds-averaged mixing model considering turbulent entropy and heat flux},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/JMT5WEMU}},
  note         = {Machine review of arXiv:2506.16296}
}
read the original abstract

In typical nature and engineering scenarios, such as supernova explosion and inertial confinement fusion, mixing flows induced by hydrodynamics 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, the majority of RANS mixing studies reported have concentrated on incompressible scenarios, with quite limited attention given to compressible cases. Moreover, most of the existing RANS mixing models demonstrate significantly inaccurate predictions for compressible mixing flow. This study develops a novel compressible RANS mixing model by incorporating physical compressibility corrections into the incompressible K-L-y mixing transition model recently proposed by Xie et al. (J. Fluid Mech., 1002, A31, 2025). Specifically, taking the density-stratified Rayleigh-Taylor mixing flows as representative compressible cases, we firstly analyze 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 and the thermodynamic Gibbs relation are employed to derive comprehensive compressibility corrections. The crucial turbulent entropy and heat flux 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 modeling scheme. This study offers a promising strategy for modeling compressible mixing flows, paving the way for more accurate predictions in complex scenarios.

Figures

Figures reproduced from arXiv: 2506.16296 by the authors.

Figure 1
Figure 1. Schematic diagrams of the initial density distributions for (a) the incompressible [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Initial distributions for density and pressure fields. The symbols and lines [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗
Figure 3
Figure 3. Temporal evolution of the mixing width predicted by (a) the [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗
Figures from the paper (7 more)
Figure 4
Figure 4. Figure 4: Budgets for terms 𝐼 ∼ 𝑉 of the TKE equation (2.16). Subfigures (a)∼(c) and (d)∼(f) correspond to moments of 𝑡/𝜏 =1, 2.5 and 4 for cases of 𝑆𝑟 = 0.5 and 3 respectively. simulations, budgets of the terms 𝐼 ∼ 𝑉 at 𝑡/𝜏 =1, 2.5, and 4 are given in figure 4, where the cases …
Figure 5
Figure 5. Figure 5: Identification of unstable regions based on the local instability criteria for [PITH_FULL_IMAGE:figures/full_fig_p012_5.png]
Figure 6
Figure 6. Figure 6: Schematic diagrams of two adjacent equilibrium states. The subfigure (a) gives [PITH_FULL_IMAGE:figures/full_fig_p013_6.png]
Figure 7
Figure 7. Figure 7: Contours of instantaneous temperature in the central [PITH_FULL_IMAGE:figures/full_fig_p015_7.png]
Figure 8
Figure 8. Figure 8: Temporal evolutions of the mixing widths. Subfigures (a) to (d) correspond to [PITH_FULL_IMAGE:figures/full_fig_p017_8.png]
Figure 9
Figure 9. Figure 9: Spatial profiles of the mean species at three different moments [PITH_FULL_IMAGE:figures/full_fig_p018_9.png]
Figure 10
Figure 10. Figure 10: Spatial profiles of the mean density at three different moments [PITH_FULL_IMAGE:figures/full_fig_p019_10.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

50 extracted references · 48 canonical work pages

  1. [1]

    , " * write output.state after.block = add.period write newline

    ENTRY address author booktitle chapter edition editor howpublished institution journal key month note number organization pages publisher school series title type volume year eprint label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts #0 'before.all := #1 'mid.sentence ...

  2. [2]

    write newline

    " write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 global.max substring 't := if while FUNCTION word.in bbl.in capitalize " " * FUNCT...

  3. [3]

    1961 Hydrodynamic and hydromagnetic stability\/

    Chandrasekhar, S. 1961 Hydrodynamic and hydromagnetic stability\/ . Oxford- C larendon P ress and N ew Y ork- O xford U niv. P ress

  4. [4]

    , Antonia, R

    Chassaing, P. , Antonia, R. , Anselmet, F. , Joly, L. & Sarkar, S. 2002 Variable Density Fluid Turbulence\/ . Springer Netherlands: Imprint: Springer

  5. [5]

    Cloutman, L. D. 1987 A new estimate of the mixing length and convective overshooting in massive stars . Astrophys. J. 313 , 699--710

  6. [6]

    Cloutman, L. D. 2003 Compressibility corrections to closure approximations for turbulent flow simulations . Tech. Rep.\/ . Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)

  7. [7]

    Craxton, R. S. , Anderson, K. S. , Boehly, T. R. , Goncharov, V. N. , Harding, D. R. , Knauer, J. P. , McCrory, R. L. , McKenty, P. W. , Meyerhofer, D. D. , Myatt, J. F. , Schmitt, A. J. , Sethian, J. D. , Short, R. W. , Skupsky, S. , Theobald, W. , Kruer, W. L. , Tanaka, K. , Betti, R. , Collins, T. J. B. , Delettrez, J. A. , Hu, S. X. , Marozas, J. A. ,...

  8. [8]

    , Rollin, B

    Denissen, N. , Rollin, B. , Reisner, J. & Andrews, M. 2014 The tilted rocket rig: A R ayleigh- T aylor test case for RANS models . J. Fluids Eng. 136 , 091301

Show all 50 references
  1. [9]

    & Tipton, R

    Dimonte, G. & Tipton, R. 2006 K-L turbulence model for the self-similar growth of the R ayleigh- T aylor and R ichtmyer- M eshkov instabilities . Phys. Fluids 18 (8), 85101

  2. [10]

    Fu, Chengquan , Zhao, Zhiye , Xu, Xin , Wang, Pei , Liu, Nansheng , Wan, Zhenhua & Lu, Xiyun 2022 Nonlinear saturation of bubble evolution in a two-dimensional single-mode stratified compressible R ayleigh- T aylor instability . Phys. Rev. Fluids 7 , 023902

  3. [11]

    Gamalii, E. G. , Rozanov, V. B. , Samarskii, A. A. , Tishkin, V. F. , Tyurina, N. N. & Favorskii, A. P. 1980 Hydrodynamic stability of compression of spherical laser targets . Sov. Phys. JETP 52 , 230--237

  4. [12]

    2013 Compressibility effects in R ayleigh- T aylor flows: influence of the stratification

    Gauthier, S. 2013 Compressibility effects in R ayleigh- T aylor flows: influence of the stratification . Phys. Scr. 2013 (T155), 014012

  5. [13]

    Gauthier, Serge 2017 Compressible R ayleigh- T aylor turbulent mixing layer between newtonian miscible fluids . J. Fluid Mech. 830 , 211–256

  6. [14]

    & Le Creurer, B

    Gauthier, S. & Le Creurer, B. 2010 Compressibility effects in R ayleigh- T aylor instability-induced flows . Philos. T. R. Soc. A 368 (1916), 1681--1704

  7. [15]

    & Menon, S

    G\' e nin, F. & Menon, S. 2010 Simulation of turbulent mixing behind a strut injector in supersonic flow . AIAA J. 48 (3), 526--539

  8. [16]

    & Glimm, J

    George, E. & Glimm, J. 2005 Self-similarity of R ayleigh- T aylor mixing rates . Phys. Fluids 17 (5), 054101

  9. [17]

    1868 On discontinuous movements of fluids

    Helmholtz, V. 1868 On discontinuous movements of fluids . Lond. Edinb. Dubl. Phil. Mag. J. Sci. 36 (244), 337--346

  10. [18]

    , Liu, X

    Jin, H. , Liu, X. F. , Lu, T. , Cheng, B. , Glimm, J. & Sharp, D. H. 2005 Rayleigh- T aylor mixing rates for compressible flow . Phys. Fluids 17 (2), 024104

  11. [19]

    1871 Hydrokinetic solutions and observations

    Kelvin, L. 1871 Hydrokinetic solutions and observations . Lond. Edinb. Dubl. Phil. Mag. J. Sci. 42 (281), 362--377

  12. [20]

    Kokkinakis, I. W. , Drikakis, D. & Youngs, D. L. 2019 Modeling of R ayleigh- T aylor mixing using single-fluid models . Phys. Rev. E 99 , 013104

  13. [21]

    Kokkinakis, I. W. , Drikakis, D. & Youngs, D. L. 2020 Two-equation and multi-fluid turbulence models for R ichtmyer- M eshkov mixing . Phys. Fluids 32 (7), 074102

  14. [22]

    Kokkinakis, I. W. , Drikakis, D. , Youngs, D. L. & Williams, R. J. R. 2015 Two-equation and multi fluid turbulence models for R ayleigh- T aylor mixing . Int. J. Heat Fluid Fl. 56 , 233--250

  15. [23]

    Li, Haifeng , He, Zhiwei , Zhang, Yousheng & Tian, Baolin 2019 On the role of rarefaction/compression waves in R ichtmyer- M eshkov instability with reshock . Phys. Fluids 31 (5), 054102

  16. [24]

    2004 Compressibility effects on the R ayleigh- T aylor instability growth between immiscible fluids

    Livescu, D. 2004 Compressibility effects on the R ayleigh- T aylor instability growth between immiscible fluids . Phys. Fluids 16 (1), 118--127

  17. [25]

    2013 Numerical simulations of two-fluid turbulent mixing at large density ratios and applications to the R ayleigh- T aylor instability

    Livescu, D. 2013 Numerical simulations of two-fluid turbulent mixing at large density ratios and applications to the R ayleigh- T aylor instability . Philos. Trans. R. Soc. A 371 (2003), 20120185

  18. [26]

    Luo, Tengfei , Li, Zhijie , Yuan, Zelong , Peng, Wenhui , Liu, Tianyuan , Wang, Liangzhu (Leon) & Wang, Jianchun 2024 Fourier neural operator for large eddy simulation of compressible R ayleigh- T aylor turbulence . Phys. Fluids 36 (7), 075165

  19. [27]

    Luo, Tengfei & Wang, Jianchun 2022 Mixing and energy transfer in compressible R ayleigh- T aylor turbulence for initial isothermal stratification . Phys. Rev. Fluids 7 , 104608

  20. [28]

    Meshkov, E. E. 1969 Instability of the interface of two gases accelerated by a shock wave . Fluid Dyn. 4 (5), 101--104

  21. [29]

    & Schilling, O

    Mor \'a n-L \'o pez, J. & Schilling, O. 2013 Multicomponent R eynolds-averaged N avier- S tokes simulations of reshocked R ichtmyer- M eshkov instability-induced mixing . High Energy Density Phys. 9 , 112--121

  22. [30]

    & Greenough, J

    Morgan, B. & Greenough, J. 2015 Large-eddy and unsteady RANS simulations of a shock-accelerated heavy gas cylinder . Shock Waves 26

  23. [31]

    Qi, Han , He, Zhi-wei , Xu, Ai-guo & Zhang, You-sheng 2024 The vortex structure and enstrophy of the mixing transition induced by R ayleigh- T aylor instability . Phys. Fluids 36 (11), 114107

  24. [32]

    1882 Investigation of the C haracter of the E quilibrium of an I ncompressible H eavy F luid of V ariable D ensity

    Rayleigh, L. 1882 Investigation of the C haracter of the E quilibrium of an I ncompressible H eavy F luid of V ariable D ensity . Proc. Lond. Math. Soc. 201 (1), 170--177

  25. [33]

    Richtmyer, R. D. 1960 Taylor instability in shock acceleration of compressible fluids . Commun. Pure Appl. Math. 13 (2), 297--319

  26. [34]

    1950 The I nstability of L iquid S urfaces when A ccelerated in a D irection P erpendicular to their P lanes

    Taylor, G. 1950 The I nstability of L iquid S urfaces when A ccelerated in a D irection P erpendicular to their P lanes. I . Proc. R. Soc. Lond., Ser. A 201 (1065), 192

  27. [35]

    , Drikakis, D

    Thornber, B. , Drikakis, D. , Williams, R. J. R. & Youngs, D. 2008 a\/ On entropy generation and dissipation of kinetic energy in high-resolution shock-capturing schemes . J. Comput. Phys. 227 (10), 4853--4872

  28. [36]

    , Mosedale, A

    Thornber, B. , Mosedale, A. , Drikakis, D. , Youngs, D. & Williams, R. J. R. 2008 b\/ An improved reconstruction method for compressible flows with low mach number features . J. Comput. Phys. 227 (10), 4873--4894

  29. [37]

    Toro, E. F. , Spruce, M. & Speares, W. 1994 Restoration of the contact surface in the HLL - R iemann solver . Shock Waves 4 , 25--34

  30. [38]

    Xiao, M. , Hu, Z. , Dai, Z. & Zhang, Y. 2022 Experimentally consistent large-eddy simulation of re-shocked R ichtmyer- M eshkov turbulent mixing . Phys. Fluids 34 (12), 125125

  31. [39]

    Xiao, Mengjuan , Qi, Han & Zhang, Yousheng 2025 Local transition indicator and modelling of turbulent mixing based on the mixing state . J. Fluid Mech. 1002 , A4

  32. [40]

    , Zhang, Y

    Xiao, M. , Zhang, Y. & Tian, B. 2020 Modeling of turbulent mixing with an improved K-L model . Phys. Fluids 32 (9), 092104

  33. [41]

    Xie, Hansong , Qi, Han , Xiao, Mengjuan , Zhang, Yousheng & Zhao, Yaomin 2025 An intermittency based R eynolds-averaged transition model for mixing flows induced by interfacial instabilities . J. Fluid Mech. 1002 , A31

  34. [42]

    , Zhao, Y

    Xie, H. , Zhao, Y. & Zhang, Y. 2023 Data-driven nonlinear K-L turbulent mixing model via gene expression programming method . Acta Mech. Sin. 39 (2), 322315

  35. [43]

    AIP Advances 11 (7), 075213

    Xie, Han-song , Xiao, Meng-juan & Zhang, You-sheng 2021 a\/ Predicting different turbulent mixing problems with the same k- model and model coefficients . AIP Advances 11 (7), 075213

  36. [44]

    , Xiao, M.-J

    Xie, H.-S. , Xiao, M.-J. & Zhang, Y.-S. 2021 b\/ Unified prediction of turbulent mixing induced by interfacial instabilities via B esnard- H arlow- R auenzahn-2 model . Phys. Fluids 33 (10), 105123

  37. [45]

    Xue, Chuang & Ye, Wenhua 2010 Destabilizing effect of compressibility on R ayleigh- T aylor instability for fluids with fixed density profile . Phys. Plasmas 17 (4), 042705

  38. [46]

    Zhang, Y. , He, Z. , Xie, H. , Xiao, M. & Tian, B. 2020 Methodology for determining coefficients of turbulent mixing model . J. Fluid Mech. 905 , A26

  39. [47]

    , Betti, R

    Zhao, D. , Betti, R. & Aluie, H. 2022 Scale interactions and anisotropy in R ayleigh- T aylor turbulence . J. Fluid Mech. 930 , A29

  40. [48]

    Zhao, Zhiye , Liu, Nansheng & Lu, Xiyun 2020 Kinetic energy and enstrophy transfer in compressible R ayleigh- T aylor turbulence . J. Fluid Mech. 904 , A37

  41. [49]

    2017 a\/ Rayleigh- T aylor and R ichtmyer- M eshkov instability induced flow, turbulence, and mixing

    Zhou, Y. 2017 a\/ Rayleigh- T aylor and R ichtmyer- M eshkov instability induced flow, turbulence, and mixing. I . Phys. Rep. 720-722 , 1--136

  42. [50]

    2017 b\/ Rayleigh- T aylor and R ichtmyer- M eshkov instability induced flow, turbulence, and mixing

    Zhou, Y. 2017 b\/ Rayleigh- T aylor and R ichtmyer- M eshkov instability induced flow, turbulence, and mixing. II . Phys. Rep. 723-725 , 1--160

Pith tools

Reviewed August 6, 2026 · model on record in the stance chip above.