Pith. sign in

REVIEW 4 major objections 5 minor 42 references

Nonlinear bubble behaviours of compressible Rayleigh-Taylor instability with isothermal stratification in cylindrical geometry

T0 review · 4 major / 5 minor · reviewed 2026-08-09 · deepseek-v4-flash

Pith's one-line read The paper proposes an additive three-part velocity model for bubble growth in compressible cylindrical Rayleigh-Taylor instability and verifies it against direct numerical simulation across Atwood and Mach numbers.

desk verdict Careful DNS gives a credible phase diagram for bubble acceleration in cylindrical RTI; the additive model is diagnostic, not predictive, but the paper deserves refereeing. read the letter →

arxiv 2502.00625 v1 pith:UUQ2XTVX submitted 2025-02-02 physics.flu-dyn

classification physics.flu-dyn
keywords Rayleigh-Taylorinstabilitycylindricalgeometrycompressibleflowisothermalstratificationbubbleaccelerationvorticityaccumulationcompressibilitydirectnumericalsimulation
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

The paper claims that the nonlinear velocity of a single-mode bubble in compressible, isothermally stratified Rayleigh-Taylor instability in cylindrical geometry is the sum of three contributions: a potential-flow velocity, a vorticity-accumulation velocity, and a compressibility velocity. It shows by direct numerical simulation that in convergent geometry the late-time bubble acceleration is dominated by vorticity accumulation at low Atwood and Mach numbers and by flow compressibility at high Atwood and Mach numbers, while in divergent geometry only flow compressibility produces robust acceleration. The improved model, which reads the simulated vorticity and dilatation fields, reproduces bubble-velocity histories from the linear through the highly nonlinear regime for Atwood numbers 0.1-0.9 and Mach numbers 0.1-0.9, for perturbation modes 8 and 16. If correct, the result gives a mechanistic bookkeeping for late-time cylindrical Rayleigh-Taylor instability and shows that planar-based models can misjudge convergent geometry.

What carries the argument

The carrying objects are the three velocity contributions in Eq. (3.8) and the diagnostic fields behind them: the potential-flow amplitude ODE of the incompressible cylindrical theory with a local Atwood number, the bubble-averaged vorticity $\bar{\omega}$ with a density-ratio correction and efficiency factor, and the mean dilatation $\bar{\theta}$ in the annulus between the bubble tip and the undisturbed front. The compressibility velocity is derived from a Green's-formula identity relating the line integral of the dilatational velocity to the area integral of $\nabla\cdot\mathbf{u}$, and the factor $(1+r_\infty/r_b)/2$ encodes the interfacial curvature, recovering the planar formula when $r_b\to r_\infty$. The model is not predictive from initial conditions alone; it reads the vorticity, the dilatation, and the front position from the simulation.

What would settle it

Measure the angular dependence of the radial velocity along the heavy-fluid front immediately ahead of the bubble tip in DNS at $A_T=0.9$ and $Ma=0.9$; if $u_r$ varies substantially with angle, the uniform-piston assumption behind $v_d$ is violated and the model's compressibility term would need revision.

Watch

Extended reading notes

Core claim

The central claim is that the bubble velocity in two-dimensional single-mode cylindrical Rayleigh-Taylor instability with isothermal stratification can be modelled as $v_b = v_p + v_{\mathrm{vort}} + v_d$ (Eq. 3.8), where $v_p$ solves the incompressible cylindrical potential-flow ODE with the local Atwood number evaluated at the bubble tip, $v_{\mathrm{vort}} = \eta\, r_\rho\, \bar{\omega} / (2n/r_b)$ accounts for accumulated vorticity inside the bubble, and $v_d = -\bar{\theta}\, (r_\infty - r_b)\,(1 + r_\infty/r_b)/2$ accounts for the compression of the heavy fluid ahead of the bubble through the mean dilatation. The paper verifies this additive decomposition against DNS for convergent and divergent accelerations across $A_T=0.1$-$0.9$ and $Ma=0.1$-$0.9$. In convergent cases the highly nonlinear acceleration is dominated by vorticity accumulation at low $A_T$ and low $Ma$ and by flow compressibility at high $A_T$ and high $Ma$; in divergent cases robust acceleration is purely due to flow compressibility. Density stratification acts to suppress acceleration at low $A_T$ and high $Ma$.

Load-bearing premise

The model's compressibility term assumes that a rising bubble of light fluid acts like a piston that uniformly compresses the heavy fluid in front of it, so the radial velocity at the bubble tip is uniform along the angular direction and equals $v_d$; if this piston picture fails, the compressibility contribution and the conclusion that flow compressibility dominates at high Atwood and Mach numbers lose their basis.

Editorial extensions

If this is right

  • Convergent and divergent bubbles follow different acceleration mechanisms: vorticity accumulation drives robust acceleration in convergent cases at low Atwood and low Mach numbers, while it does not produce robust acceleration in divergent cases.
  • At high Atwood and high Mach numbers, flow compressibility becomes the dominant destabilizing effect and causes transient acceleration in divergent cases that is absent from the incompressible potential-flow description.
  • The phase diagram in the Atwood-Mach plane for cylindrical Rayleigh-Taylor instability differs from the planar stratified case, so planar-based late-time models should not be applied directly to cylindrical or spherical configurations.
  • The model reproduces DNS bubble-velocity histories for both $n=8$ and $n=16$ perturbation modes, supporting the additive decomposition as a general characterization for single-mode cylindrical compressible Rayleigh-Taylor instability.
  • Because the model requires simulated vorticity and dilatation fields, it is a diagnostic decomposition rather than a self-contained predictive model.

Reading between the lines

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

  • Beyond the paper: the curvature factor $(1+r_\infty/r_b)/2$ implies that the compressibility contribution should strengthen as the convergent bubble approaches the centre and weaken as a divergent bubble moves outward; this predicted scaling could be tested directly from the DNS data by plotting $v_d$ against $1/r_b$.
  • Beyond the paper: the vorticity efficiency factor $\eta=0.4$ differs from the planar value 0.45, suggesting that vortex attenuation at the bubble tip depends on geometry; a systematic sweep over mode number, Reynolds number, and initial amplitude could calibrate this factor and reveal whether it is universal.
  • Beyond the paper: if the same additive decomposition carries over to spherical geometry with an appropriately modified curvature factor, it would offer a route to modelling supernova and inertial-confinement-fusion mixing layers without fully three-dimensional simulation.
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

4 major / 5 minor

Summary. The paper reports direct numerical simulations of two-dimensional single-mode compressible Rayleigh–Taylor instability in cylindrical geometry with isothermal stratification, covering Atwood numbers A_T = 0.1–0.9 and Mach numbers Ma = 0.1–0.9. It documents the nonlinear and highly nonlinear bubble dynamics for convergent and divergent accelerations, classifies the late-time response into robust acceleration, transient acceleration, and non-acceleration phases, and proposes a composite model, Eq. (3.8), that adds a vorticity-accumulation term (3.5) and a compressibility term (3.7) to the incompressible potential-flow ODE of Zhao et al. (2020b), with a density-ratio correction at the bubble tip. The model is compared with DNS for mode numbers B = 8 and 16 and reported to reproduce the bubble-velocity evolution from linear to highly nonlinear regimes.

Significance. The DNS dataset and the documented differences between convergent and divergent cylindrical RTI are a useful contribution, and the paper is careful in several respects: a grid-convergence check is performed at the most demanding parameters (A_T = 0.9, Ma = 0.9), the incompressible limit is validated against the nonlinear theory of Zhao et al. (2020b), and the model is exercised for two mode numbers. The proposed additive decomposition, if it can be supported by an independent validation, would provide a valuable physical interpretation of late-time bubble acceleration in convergent geometries. However, the model as presented is not predictive—the authors themselves state this in Section 4—and the validation is partly circular because the added terms are computed from the same DNS.

major comments (4)
  1. [Sec. 3.2, Eq. (3.8), Figs. 6–7] The central validation claim is weakened by circularity between the model terms and the DNS. The vorticity velocity v_vort in Eq. (3.5) is computed from the DNS-measured mean vorticity and the density ratio at the bubble tip, and the compressibility velocity v_d in Eq. (3.7) is computed from the DNS-measured mean dilatation. Therefore, comparing v_p + v_vort + v_d with the DNS bubble velocity in Figs. 6 and 7 mainly checks the internal consistency of the definitions, not the predictive power of the model. The statement in Section 4 that the model 'requires the simulation and thus is not predictive' is an honest acknowledgment, but it is in tension with the abstract's claim that the model is 'verified by numerical results.' The paper should either provide an independent test (e.g., calibrating eta on a subset of cases and validating on the remaining ones) or explicitly frame the model as a diagnostic decomposition.
  2. [Sec. 3.2, Eq. (3.7)] The derivation of v_d from Green's formula assumes that the dilatational radial velocity at the bubble-tip radius is angularly uniform and equal to v_d (the 'piston' picture). In cylindrical geometry with mode B = 8 the interface is curved, and u_r^d at r_b can vary with angle phi; the boundary integral in Eq. (3.6) equals -2 pi r_b v_d only if the angular mean of u_r^d happens to be v_d. No angular-uniformity check is reported, and no sensitivity study on the choice of r_infinity (the computational boundary rather than the actual hydrostatic front) is reported. Since v_d is central to the compressibility-dominated phase classification in Fig. 5, the model agreement shown in Figs. 6 and 7 is conditional on this unverified assumption. I ask the authors to measure the angular distribution of u_r^d at r_b and to test the sensitivity of v_d to r_infinity.
  3. [Eq. (3.5), Sec. 4] The efficiency factor eta = 0.4 in Eq. (3.5) is introduced as an empirical constant 'to account for the attenuation of vortices in cylindrical geometry.' No derivation, prior calibration, or sensitivity analysis is provided. Since Figs. 6 and 7 use the same eta, the B = 16 cases provide a consistency check but not an independent test of this parameter. Please discuss the physical basis for eta = 0.4 and quantify how the model agreement degrades as eta is varied.
  4. [Fig. 5, Sec. 3.1] The phase diagram in Fig. 5 classifies the dominant mechanism in the highly nonlinear stage based on the time-averaged magnitudes of v_vort and v_d, both of which are defined from DNS fields. The finding that 'acceleration in divergent cases is purely induced by flow compressibility' is therefore a statement about the relative sizes of the model terms, not a direct measurement of the physical mechanism. The paper should clarify this distinction, and ideally corroborate the classification with an independent flow diagnostic (e.g., the correlation between local dilatation and bubble acceleration, or the work done by pressure dilatation).
minor comments (5)
  1. [Sec. 2.2, Eq. (2.6)] The expression for the initial density and pressure is difficult to parse after typesetting; please rewrite it with proper parentheses, for example by explicitly writing the exponent as sgn(a*) Ma^2 (1 +/- A_T) (r - eta(phi)) B/(2 pi).
  2. [Sec. 3.2, above Eq. (3.5)] The definition of the averaged vorticity uses 'the volume V inside the bubble'; since the simulations are two-dimensional, please clarify whether this is a surface area and how the bubble interior is delimited.
  3. [Fig. 3] The text refers to 'figure 3(a,b)' and 'figure 3(c,d)' when describing vortex positions; according to the caption, the convergent cases at Ma = 0.1 are panels (a) and (c), and the divergent cases are panels (b) and (d). Please correct the panel references.
  4. [Sec. 2.2, Fig. 1] The grid-convergence study is reported only for the most demanding parameters (A_T = 0.9, Ma = 0.9). A second convergence test at, for example, A_T = 0.1 and Ma = 0.9 would help establish that the resolution is adequate across the parameter range.
  5. [Throughout] The word 'initial' is misspelled as 'intial' in Section 2.2 and Section 3.1; please proofread the manuscript for typographical errors.

Circularity Check

2 steps flagged · score 8.0 of 10

Model (3.8) reproduces DNS using DNS-measured vorticity, density ratio, and dilatation, with a chosen efficiency factor; the compressibility term imports a planar piston ansatz from self-cited work.

  1. fitted input called prediction [Sec. 3.2, Eqs. (3.5), (3.7), (3.8); Sec. 4 conclusion]
    "The efficiency factor η = 0.4 is used here to account for the attenuation of vortices in cylindrical geometry. ... By adding (3.5) and (3.7) to v_p, the bubble velocity can be completed as the following model: v_b = v_p + v_vort + v_d. (3.8) ... It is noted that current model inheriting the illuminating idea of the Betti-Sanz model requires the simulation and thus is not predictive."

    The model's three right-hand-side terms are all constructed from data taken from the same DNS whose bubble velocity is then said to be reproduced. v_vort uses the DNS-measured mean vorticity and bubble-tip density ratio, with an efficiency factor η = 0.4 chosen for cylindrical geometry; v_d uses the DNS-measured mean dilatation and DNS-determined radii r_b and r_∞; and v_p is obtained by solving the ODE with a local Atwood number A'_T = (1 - r_F)/(1 + r_F) measured from the DNS density field. The agreement in Figs. 6-7 is therefore a diagnostic reconstruction rather than an independent prediction. The paper's own concluding statement that the model 'requires the simulation and thus is not predictive' makes this fitted-input nature explicit.

  2. ansatz smuggled in via citation [Sec. 3.2, Eqs. (3.6)-(3.7)]
    "Given that a rising bubble of light fluid acts like a piston uniformly compressing its front heavy fluid (Luo et al. 2020; Fu et al. 2022, 2023), the compressing velocity u_r^d at the radial position r_b is almost uniform along the periodic φ-direction and can be approximated as v_d (Fu et al. 2023)."

    Equation (3.7) is obtained by replacing the boundary integral ∮(-u_r^d r dφ) with -2π r_b v_d, which is valid only if u_r^d is independent of φ at r_b and equals v_d. In the cylindrical setup the interface is curved, r = r0 + η0 cos(Bφ), so u_r^d at r_b is not demonstrated to be angularly uniform. The uniform-piston assumption is imported, via citation, from the same authors' planar compressible RTI papers, and no angular-profile check or sensitivity to the choice of r_∞ is reported. Thus the compressibility contribution and the subsequent 'flow-compressibility-dominated' phase classification rest on a self-cited planar ansatz rather than on a verified cylindrical derivation.

full rationale

The base nonlinear ODE (3.1)-(3.4) is independently validated in Fig. 1(a) against DNS at Ma = 0.1, so invoking Zhao et al. (2020b) is not by itself circular. The circularity appears in the construction and verification of the 'improved model' (3.8). The vorticity and compressibility contributions are evaluated from DNS-measured fields and a chosen efficiency factor, and the model is then 'verified' against the very DNS that supplied those fields. The authors' admission in the conclusion that the model 'requires the simulation and thus is not predictive' confirms that the central claim is a diagnostic description, not an out-of-sample prediction. Additionally, the derivation of v_d depends on an angular-uniform piston ansatz that is inherited from planar self-cited work and is not checked in cylindrical geometry. These issues together make the central quantitative claim substantially circular, so the score is 8.

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

The model rests on the cylindrical potential-flow ODE plus four ad-hoc extensions: local-Atwood substitution, linear superposition, piston-compression for the dilatation term, and a hand-tuned vortex efficiency eta=0.4. Crucially, the vorticity and compressibility inputs are measured from the same DNS fields used for validation, so the ledger is not closed independently.

free parameters (1)
  • efficiency factor eta = 0.4
    Chosen in Eq. (3.5) to account for attenuation of vortices in cylindrical geometry; no derivation or sensitivity study is provided, and model agreement depends on it.
assumptions (5)
  • domain assumption The Zhao et al. (2020b) potential-flow nonlinear ODE, Eqs. (3.1)-(3.4), is a valid base description for cylindrical RTI bubble evolution.
    The improved model builds directly on this ODE, which was derived under uniform-density, incompressible, potential-flow assumptions.
  • ad hoc to paper Bubble velocity is a linear superposition v_b = v_p + v_vort + v_d, Eq. (3.8).
    The three contributions are added without a derivation showing that the mechanisms are independent; no cross-coupling or nonlinear interaction terms are included.
  • domain assumption A rising bubble acts as a piston uniformly compressing its front heavy fluid, so the radial velocity at r_b is uniform and equals v_d, Eq. (3.7).
    Used to relate the contour integral in Green's formula to v_d; if this piston picture is inaccurate, v_d is not a reliable measure of compressibility-induced velocity.
  • ad hoc to paper The efficiency factor eta=0.4 in Eq. (3.5) correctly represents vortex attenuation in cylindrical geometry.
    Hand-chosen, not derived from flow parameters or validated independently; model agreement is sensitive to this value.
  • ad hoc to paper Density variation at the bubble tip can be represented by replacing A_T with a local Atwood number A'_T in Eqs. (3.1)-(3.4).
    The original ODE has a constant background density; using a time-varying local density ratio to inject stratification is a plausible but non-derived extension.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Nonlinear bubble behaviours of compressible Rayleigh-Taylor instability with isothermal stratification in cylindrical geometry." pith.science (2026). https://pith.science/paper/UUQ2XTVX

@misc{pith2026250200625,
  author       = {Pith},
  title        = {Pith review of: Nonlinear bubble behaviours of compressible Rayleigh-Taylor instability with isothermal stratification in cylindrical geometry},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/UUQ2XTVX}},
  note         = {Machine review of arXiv:2502.00625}
}
abstract

Nonlinear evolutions of two-dimensional single-mode compressible Rayleigh--Taylor instability (RTI) with isothermal stratification are investigated in cylindrical geometry via direct numerical simulation for different Atwood numbers ($A_T=0.1-0.9$) and Mach numbers ($Ma=0.1-0.9$). It is found that the nonlinear bubble growth involves the effects of density stratification, vorticity accumulation and flow compressibility and shows considerable differences between convergent (acceleration acting radially inward) and divergent (acceleration acting radially outward) cases. Specifically, the density stratification leads to non-acceleration at low $A_T$ and high $Ma$. The accelerations in convergent cases are dominated by vorticity accumulation at low $A_T$ and low $Ma$ and by flow compressibility at high $A_T$ and high $Ma$ whereas the accelerations in divergent cases are purely induced by flow compressibility at high $A_T$ and high $Ma$. Based on the nonlinear theory of incompressible cylindrical RTI with uniform-density background~(Zhao et al., J. Fluid Mech., vol. 900, 2020, A24), an improved model is proposed by taking the density variation, vorticity accumulation and flow compressibility into consideration. This model is verified by numerical results and well reproduces the bubble evolution for different $A_T$ and $Ma$ from linear to highly nonlinear regimes.

Figures

Figures reproduced from arXiv: 2502.00625 by the authors.

Figure 1
Figure 1. Time-varying positions of bubble tips for convergent ( > 1) and divergent ( < 1) cases at () = 0.9 and = 0.1 with symbols representing the simulation data and at () = 0.9 and = 0.9 with three grid resolutions: 2562 (dashed lines), 5122 (dot-dashed lines) and 10242 (long-dashed lines). The solid lines denote the nonlinear theory (Zhao et al. 2020b). dynamics in cylindrical RTI evolution. To obtain fine flow field cha… view at source ↗
Figure 2
Figure 2. Temporal bubble velocity at different of () = 0.1 and of () = 0.9. The solid and long-dashed lines in each panel represent the bubble velocities of convergent and divergent cases for = 8, respectively. The horizontal black dot-dashed lines in each panel denote the saturation velocity obtained from potential flow theory (Goncharov 2002). of vortex generation around the sinking spike drift further away from the bubble… view at source ↗
Figure 3
Figure 3. Vorticity contours for the cases into highly nonli [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figures from the paper (4 more)
Figure 4
Figure 4. Figure 4: Density profiles along the radial lines across bubb [PITH_FULL_IMAGE:figures/full_fig_p009_4.png]
Figure 5
Figure 5. Figure 5: Bubble acceleration phase diagram of − of () convergent and of () divergent cases for = 8 in highly nonlinear stage. Triangle, circle and square represent the phases of non-acceleration, transient acceleration and robust acceleration, respectively. Red and blue symbols…
Figure 6
Figure 6. Figure 6: Bubble velocities obtained from DNS results (symb [PITH_FULL_IMAGE:figures/full_fig_p012_6.png]
Figure 7
Figure 7. Figure 7: Bubble velocities obtained from DNS results (symb [PITH_FULL_IMAGE:figures/full_fig_p013_7.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

42 extracted references · 40 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]

    & Menon, S

    Balakrishnan, K. & Menon, S. 2010 On turbulent chemical explosions into dilute aluminum particle clouds . Combust. Theory Model. 14 (4), 583--617

  4. [4]

    Bell, G. I. 1951 Taylor instability on cylinders and spheres in the small amplitude approximation . Report No. LA-1321, LANL

  5. [5]

    Bell, J. B. , Day, M. S. , Rendleman, C. A. , Woosley, S. E. & Zingale, M. 2004 Direct numerical simulations of type Ia supernovae flames. II. The Rayleigh--Taylor instability . Astrophys. J. 608 (2), 883

  6. [6]

    & Hurricane, O

    Betti, R. & Hurricane, O. A. 2016 Inertial-confinement fusion with lasers . Nat. Phys. 12 (5), 435--448

  7. [7]

    & Sanz, J

    Betti, R. & Sanz, J. 2006 Bubble acceleration in the ablative Rayleigh--Taylor instability . Phys. Rev. Lett. 97 (20), 205002

  8. [8]

    , Aluie, H

    Bian, X. , Aluie, H. , Zhao, D. , Zhang, H. & Livescu, D. 2020 Revisiting the late-time growth of single-mode Rayleigh--Taylor instability and the role of vorticity . Physica D 403 , 132250

Show all 42 references
  1. [9]

    , Lebedev, V

    Chertkov, M. , Lebedev, V. & Vladimirova, N. 2009 Reactive Rayleigh--Taylor turbulence . J. Fluid Mech. 633 , 1--16

  2. [10]

    , Zhao, Z

    Fu, C.-Q. , Zhao, Z. , Wang, P. , Liu, N.-S. , Wan, Z.-H. & Lu, X.-Y. 2023 Bubble re-acceleration behaviours in compressible Rayleigh--Taylor instability with isothermal stratification . J. Fluid Mech. 954 , A16

  3. [11]

    , Zhao, Z

    Fu, C.-Q. , Zhao, Z. , Xu, X. , Wang, P. , Liu, N.-S. , Wan, Z.-H. & Lu, X.-Y. 2022 Nonlinear saturation of bubble evolution in a two-dimensional single-mode stratified compressible Rayleigh--Taylor instability . Phys. Rev. Fluid 7 (2), 023902

  4. [12]

    2017 Compressible Rayleigh--Taylor turbulent mixing layer between Newtonian miscible fluids

    Gauthier, S. 2017 Compressible Rayleigh--Taylor turbulent mixing layer between Newtonian miscible fluids . J. Fluid Mech. 830 , 211--256

  5. [13]

    Ge, J. , Li, H. , Zhang, X. & Tian, B. 2022 Evaluating the stretching/compression effect of Richtmyer--Meshkov instability in convergent geometries . J. Fluid Mech. 946 , A18

  6. [14]

    Goncharov, V. N. 2002 Analytical model of nonlinear, single-mode, classical Rayleigh--Taylor instability at arbitrary Atwood numbers . Phys. Rev. Lett. 88 (13), 134502

  7. [15]

    Hester, J. J. 2008 The Crab Nebula: an astrophysical chimera . Annu. Rev. Astron. Astrophys. 46 (1), 127--155

  8. [16]

    Hicks, E. P. 2014 A shear instability mechanism for the pulsations of Rayleigh--Taylor unstable model flames . J. Fluid Mech. 748 , 618--640

  9. [17]

    Houseman, G. A. & Molnar, P. 1997 Gravitational (Rayleigh--Taylor) instability of a layer with non-linear viscosity and convective thinning of continental lithosphere . Geophys. J. Int. 128 (1), 125--150

  10. [18]

    , Miyagoshi, T

    Isobe, H. , Miyagoshi, T. , Shibata, K. & Yokoyama, T. 2005 Filamentary structure on the Sun from the magnetic Rayleigh--Taylor instability . Nature 434 (7032), 478--481

  11. [19]

    & Shvarts, D

    Kishony, R. & Shvarts, D. 2001 Ignition condition and gain prediction for perturbed inertial confinement fusion targets . Phys. Plasmas 8 (11), 4925--4936

  12. [20]

    1955 On the instability of superposed fluids in a gravitational field

    Layzer, D. 1955 On the instability of superposed fluids in a gravitational field . Astrophys. J. 122 , 1

  13. [21]

    , Landen, O

    Lindl, J. , Landen, O. , Edwards, J. , Moses, E. & team, N.I.C. 2014 Review of the national ignition campaign 2009-2012 . Phys. Plasmas 21 (2), 020501

  14. [22]

    2004 Compressibility effects on the Rayleigh--Taylor instability growth between immiscible fluids

    Livescu, D. 2004 Compressibility effects on the Rayleigh--Taylor instability growth between immiscible fluids . Phys. Fluids 16 (1), 118--127

  15. [23]

    2013 Numerical simulations of two-fluid turbulent mixing at large density ratios and applications to the Rayleigh--Taylor instability

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

  16. [24]

    & Wang, J

    Luo, T. & Wang, J. 2021 Effects of Atwood number and stratification parameter on compressible multi-mode Rayleigh--Taylor instability . Phys. Fluids 33 (11), 115111

  17. [25]

    & Wang, J

    Luo, T. & Wang, J. 2022 Mixing and energy transfer in compressible Rayleigh--Taylor turbulence for initial isothermal stratification . Phys. Rev. Fluid 7 (10), 104608

  18. [26]

    , Wang, J

    Luo, T. , Wang, J. , Xie, C. , Wan, M. & Chen, S. 2020 Effects of compressibility and Atwood number on the single-mode Rayleigh--Taylor instability . Phys. Fluids 32 (1), 012110

  19. [27]

    1883 Investigation of the character of the equilibrium of an incompressible heavy fluid of variable density

    Rayleigh, L. 1883 Investigation of the character of the equilibrium of an incompressible heavy fluid of variable density . Proc. Lond. Math. Soc. s1--14 (1), 170--177

  20. [28]

    Reckinger, S. J. , Livescu, D. & Vasilyev, O. V. 2016 Comprehensive numerical methodology for direct numerical simulations of compressible Rayleigh--Taylor instability . J. Comput. Phys. 313 , 181--208

  21. [29]

    & Nishihara, K

    Sakagami, H. & Nishihara, K. 1990 Rayleigh--Taylor instability on the pusher--fuel contact surface of stagnating targets . Phys. Fluids B 2 (11), 2715--2730

  22. [30]

    Sykes, J. P. , Gallagher, T. P. & Rankin, B. A. 2021 Effects of Rayleigh--Taylor instabilities on turbulent premixed flames in a curved rectangular duct . Proc. Combust. Inst. 38 (4), 6059--6066

  23. [31]

    Taylor, G. I. 1950 The instability of liquid surfaces when accelerated in a direction perpendicular to their planes. I. Proc. R. Soc. Lond. A 201 (1065), 192--196

  24. [32]

    , Wu, J.-F

    Wang, L.-F. , Wu, J.-F. , Ye, W.-H. , Zhang, W.-Y. & He, X.-T. 2013 Weakly nonlinear incompressible Rayleigh--Taylor instability growth at cylindrically convergent interfaces . Phys. Plasmas 20 (4), 042708

  25. [33]

    Weir, S. T. , Chandler, E. A. & Goodwin, B. T. 1998 Rayleigh--Taylor instability experiments examining feedthrough growth in an incompressible, convergent geometry . Phys. Rev. Lett. 80 (17), 3763

  26. [34]

    , Liu, H

    Wu, J. , Liu, H. & Xiao, Z. 2021 Refined modelling of the single-mode cylindrical Richtmyer--Meshkov instability . J. Fluid Mech. 908 , A9

  27. [35]

    & Livescu, D

    Yu, H. & Livescu, D. 2008 Rayleigh--Taylor instability in cylindrical geometry with compressible fluids . Phys. Fluids 20 (10), 104103

  28. [36]

    , Zhao, Z

    Yuan, M. , Zhao, Z. , Liu, L. , Wang, P. , Liu, N.-S. & Lu, X.-Y. 2023 Instability evolution of a shock-accelerated thin heavy fluid layer in cylindrical geometry . J. Fluid Mech. 969 , A6

  29. [37]

    , Betti, R

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

  30. [38]

    , Liu, N.-S

    Zhao, Z. , Liu, N.-S. & Lu, X.-Y. 2020 a\/ Kinetic energy and enstrophy transfer in compressible Rayleigh--Taylor turbulence . J. Fluid Mech. 904 , A37

  31. [39]

    , Wang, P

    Zhao, Z. , Wang, P. , Liu, N.-S. & Lu, X.-Y. 2020 b\/ Analytical model of nonlinear evolution of single-mode Rayleigh--Taylor instability in cylindrical geometry . J. Fluid Mech. 900 , A24

  32. [40]

    , Wang, P

    Zhao, Z. , Wang, P. , Liu, N.-S. & Lu, X.-Y. 2021 Scaling law of mixing layer in cylindrical Rayleigh--Taylor turbulence . Phys. Rev. E 104 (5), 055104

  33. [41]

    2017 a\/ Rayleigh--Taylor and Richtmyer--Meshkov instability induced flow, turbulence, and mixing

    Zhou, Y. 2017 a\/ Rayleigh--Taylor and Richtmyer--Meshkov instability induced flow, turbulence, and mixing. I. Phys. Rep. 720--722 , 1--136

  34. [42]

    2017 b\/ Rayleigh--Taylor and Richtmyer--Meshkov instability induced flow, turbulence, and mixing

    Zhou, Y. 2017 b\/ Rayleigh--Taylor and Richtmyer--Meshkov instability induced flow, turbulence, and mixing. II. Phys. Rep. 723--725 , 1--160

Pith tools

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