Pith. sign in

REVIEW 3 major objections 3 minor 23 references

Nonlinear asymptotic bubble growth in single-mode spherical Rayleigh-Taylor instability

T0 review · 3 major / 3 minor · reviewed 2026-08-03 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read This paper claims that a Rayleigh–Taylor bubble on a spherical interface does not saturate to a terminal velocity; instead, in the nonlinear regime it enters a uniform-acceleration stage governed by |η0'|/sqrt(r0+η0) → constant and η0'' → c

desk verdict Worth a serious look: novel spherical-geometry acceleration law with honest DNS, but the printed potential ansatz in (6)-(7) is not harmonic and the central derivation needs fixing before the result can be trusted. read the letter →

arxiv 2607.29260 v1 pith:WDN2SM2S submitted 2026-07-31 physics.flu-dyn physics.plasm-ph

De-Hua Zhang , Shi-Heng Wang , Ke-Jian Qian , Zhu-Jun Li , Rui Yan , Hang Ding This is my paper · ORCID
classification physics.flu-dynphysics.plasm-ph
keywords Rayleigh-Taylorinstabilitysphericalgeometrybubblegrowthasymptoticaccelerationpotential-flowmodelconverginggravitydivergingsingle-modeperturbation
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

Rayleigh–Taylor bubbles on spherical interfaces behave differently from their flat-interface counterparts: instead of settling to a constant terminal velocity, they enter a stage of uniform acceleration. The paper claims that in the nonlinear regime the bubble-tip satisfies |η0'|/sqrt(r0+η0) → constant and η0'' → constant, with the constant fixed by mode degree, Atwood number, and gravity. This matters for long-time predictions in spherical inertial-confinement-fusion configurations and other converging/diverging flows, where assuming a flat-interface terminal speed would misstate late-time bubble growth.

What carries the argument

The central object is a two-unknown local model of the bubble tip: the interface near the pole is written as η(θ,t)=η0(t)+η2(t)θ^2, and each fluid's velocity potential is represented by the single Laplace-equation eigenmode selected for the spherical mode l, with the light-fluid potential allowed to be singular at the center but used only near the tip. Imposing the kinematic and dynamic boundary conditions yields a closed pair of ordinary differential equations for η0 and η2. Their large-time integration produces the asymptotic balance (r0+η0)A = (r0+η0)η0''B + (η0')^2 C, from which the constant-acceleration result follows. All of the geometry and density-ratio dependence is carried in the t

What would settle it

Run a single-mode spherical RTI simulation or experiment with, say, l=10 and Atwood number 0.5 and measure the bubble-tip acceleration η0'' long after the linear stage. The paper's claim requires it to approach A/(B+2C) computed from their coefficients; observing a terminal velocity (η0'' → 0 with η0' constant) or an acceleration that drifts with radius would falsify the central claim.

Watch

Extended reading notes

Core claim

Spherical geometry introduces an expanding (converging-gravity) or shrinking (diverging-gravity) intrinsic length scale, so the natural asymptotic quantity is the ratio of bubble-tip velocity to the square root of the instantaneous radius, |η0'|/sqrt(r0+η0), and this ratio—not the velocity itself—approaches a constant. Equivalently, the tip acceleration tends to a positive constant A/(B+2C), meaning outward bubbles keep accelerating and inward bubbles keep decelerating without saturation. The coefficients A, B, C are explicit algebraic functions of the mode number l, Atwood number, and gravity for the two configurations. The same local-expansion model reproduces exponential linear growth and

Load-bearing premise

The load-bearing premise is that the bubble-tip interface remains adequately described by the truncated local form η0(t)+η2(t)θ^2 with single-mode Legendre potentials; if higher harmonics become dynamically important in the asymptotic stage, the closed two-equation model and the constant-acceleration conclusion would not hold.

Editorial extensions

If this is right

  • Long-time bubble-radius or mixing-width estimates for spherical RTI should use uniform acceleration, not a terminal velocity; flat-interface models would understate late-time growth in converging gravity and overstate it in diverging gravity.
  • For the same effective wavenumber, the sphere's asymptotic acceleration is larger than the cylinder's, and the ratio approaches 3 for high modes independently of Atwood number, so spherical convergence is a stronger late-time driver than cylindrical convergence.
  • The model unifies linear and nonlinear stages, recovering the exponential growth of spherical RTI in the linear phase and the new acceleration law after several e-folding times.
  • In the high-mode limit the accelerating solution reduces to the classical constant bubble-velocity law for a tube, so the new result is consistent with the flat/tube picture at large l.

Reading between the lines

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

  • If the 3:1 sphere-to-cylinder asymptotic ratio persists at higher modes, late-time bubble growth in implosion geometries could be significantly faster than cylindrical prototypes suggest; extending the DNS to l=12–16 and larger Atwood numbers would test that directly.
  • The model omits surface tension and viscosity; these may cut off the acceleration when the bubble tip approaches the center in the diverging-gravity case, since the singular light-fluid potential is only a local approximation. A regular inner-region model could reveal where the acceleration law breaks.
  • A practical consequence left implicit: for double-cone ignition targets, estimates of bubble-spike impact timing near the axis should be revised relative to planar estimates, because the bubble keeps accelerating instead of coasting.
  • The constant-acceleration law suggests a simple late-time similarity: bubble-tip displacement grows like t^2 with a coefficient set by A/(B+2C); this could be used as a reduced-order input for mix models.
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

3 major / 3 minor

Summary. The paper develops a potential-flow model for the nonlinear evolution of a single-mode Rayleigh-Taylor bubble in spherical geometry, covering both converging-gravity (CG) and diverging-gravity (DG) configurations. The interface near the bubble tip is represented by η(θ,t)=η0(t)+η2(t)θ², and velocity potentials are written in single-mode Legendre form with three unknown coefficients. This leads to a closed two-variable ODE system for η0 and η2: Eqs. (8)-(9) for CG and (11)-(12) for DG. An asymptotic analysis of this system gives |η̇0|/√(r0+η0) → √(2A/(B+2C)) and η̈0 → A/(B+2C), i.e., a constant-acceleration/deceleration regime rather than a planar-style terminal velocity. The model is compared with axisymmetric Euler DNS for AT=0.3,0.6 and l=4,8, with reported favorable agreement. The paper also claims that spherical geometry enhances bubble growth relative to cylindrical geometry with the same effective wavenumber in CG, and mitigates it in DG, with the acceleration ratio approaching 3 for large l.

Significance. If the derivation is correct, the result is a meaningful extension of earlier cylindrical-geometry analyses: it provides closed-form asymptotic coefficients depending only on l, AT, and g, with no fitted parameters, and it makes falsifiable predictions about the late-time scaling of the bubble velocity and acceleration. The recovery of Layzer's planar terminal-velocity limit in the large-l limit is a nice consistency check. The paper also credits and builds on prior work by Goncharov, Zhao et al., and Layzer. However, as printed, the analytical derivation contains internal errors that undermine the central claim until they are repaired: the velocity potentials in Eqs. (6)-(7) are not harmonic, and the ODEs in Eqs. (9)/(12) are dimensionally inconsistent with the coefficient definitions in Eqs. (10)/(13).

major comments (3)
  1. [§II.B, Eqs. (6)-(7)] The printed potentials do not satisfy ∇²φ=0. For a mode f(r)P_l(cosθ), Laplace's equation requires f∝r^l or r^{-(l+1)}. In the CG case, φ_H with S_H=l+1 is r^{-(l+1)}P_l (harmonic), but φ_L with S_L=-l is r^{-l}P_l, whose Laplacian is -2l r^{-l-2}P_l ≠ 0 for l>0. In the DG case, S_H=l gives the non-harmonic r^{-l}P_l, while S_L=-(l+1) is harmonic. The isotropic b2 η0/r term cannot cancel this angularly dependent failure. Since Eqs. (8)-(13) are derived by substituting these potentials into (3)-(5), the ODE system and the asymptotic laws (20)-(21) are not consequences of the written ansatz. The likely correction is a sign change in the inner-region exponent (S_L=l for CG, S_H=-l for DG); the favorable DNS comparison suggests the implemented model may already use the corrected exponents, but the manuscript as printed is internally inconsistent.
  2. [§II.B, Eqs. (9)-(13)] There is a separate dimensional inconsistency. In Eq. (9), H1 as defined in (10a) has dimensions of length² (it contains (r0+η0)², η2², and mixed products), and [8η2-l(r0+η0)] has dimensions of length, so H1/8 [8η2-l(r0+η0)] η̈0 has L⁴/T². The final term AT g η2 has L²/T². The H2 term in (9) has L⁶/T². The same problem appears in Eq. (12) with (13). After substituting the asymptotic η2 relation (15), the terms in (9) scale as (r0+η0)³ and (r0+η0)⁴, whereas the claimed reduced equation (16) is linear in R=r0+η0 for the η̈ term and independent of R for the η̇² term. The printed equations can only be valid if factors of 1/(r0+η0)² and 1/(r0+η0)⁴ are missing from the H1 and H2 terms. This also affects the claimed recovery of the linear growth rate (1). The authors should correct the equations and re-verify the derivation.
  3. [§II.B, before Eqs. (8)-(12)] The model rests on the near-tip truncation η=η0+η2θ² together with single-mode Legendre potentials. This is a standard Layzer-type ansatz, but the paper provides no estimate of the neglected θ⁴ terms, which could become dynamically important in the asymptotic stage. The asymptotic relation (15) and the constant-acceleration law (20)-(21) depend directly on this truncation. I recommend adding a consistency check—for example, retaining a θ⁴ coefficient and verifying that it remains small compared to η2, or comparing the predicted local interface shape with the DNS interface in the late-time regime—so that the validity of the asymptotic claim does not rest solely on the visual agreement in Figs. 4-5.
minor comments (3)
  1. [§III, Figs. 4-5] The quantitative support for the central claim is weak: only four parameter combinations are shown, and the agreement is assessed visually. Please report a quantitative measure of agreement (e.g., relative L2 error of η0(t) and of |η̇0|/√(r0+η0) relative to the predicted asymptotic value), and provide the grid-convergence results mentioned in the text.
  2. [§II.C, Eq. (19)-(20)] The sign α in (20) is introduced without explanation. It would be helpful to state explicitly that α=-1 in the DG case follows from η̇0<0 for an inward-moving bubble, and that the square-root branch is chosen accordingly.
  3. [References] The Chandrasekhar reference appears as an unnumbered entry after the conclusion; the reference list otherwise starts with Goncharov. Please format consistently.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the asymptotic acceleration law is derived analytically from the local potential-flow ansatz and is checked against DNS, not fitted or assumed.

full rationale

Walking the derivation chain: Section II.B begins with Laplace's equation and the boundary conditions (3)-(5), then introduces an explicit local expansion for the interface and assumed forms for the velocity potentials (6)-(7). The coupled ODEs (8)-(12) follow by substitution, and the coefficients A, B, C in (17)-(18) are computed analytically from the model inputs (l, A_T, g, r0). The asymptotic relation (20) and constant-acceleration result (21) are obtained by integrating the resulting ODE (16) and taking the stated limits; no parameter is fitted to DNS or to the target asymptotic value. In Section III, DNS is used as an external validation with independently chosen initial conditions and Froude number, not as a source of model coefficients. The self-citations (e.g., Qian et al. 2026) are contextual and not load-bearing for the derivation; the cited linear growth rates and cylindrical comparison come from distinct external works. The borrowings from Layzer/Goncharov and Zhao et al. are explicit methodological precedents rather than hidden inputs. The possible typo in the printed exponents of (6)-(7) would be a mathematical correctness issue, not circularity: the claimed result is not identical to an input by construction. Overall, the central claim is self-contained analytically and externally checked.

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

The model is parameter-free apart from physical inputs r0, l, AT, g; the burden lies in the local-flow representation and truncation, not in fitted constants. No new physical entities are postulated.

assumptions (5)
  • domain assumption Both fluids are inviscid, incompressible, irrotational, with no surface tension, so velocity potentials satisfy Laplace's equation.
    Invoked in Section II.B before equation (2); excludes compressibility, viscosity and surface-tension effects present in real ICF/droplet flows.
  • ad hoc to paper The interface shape near the bubble tip is truncated as eta(theta,t)=eta0(t)+eta2(t)theta^2, and velocity potentials take the single-mode eigenfunction forms (6)-(7).
    Used to close the system in equations (8)-(12); no bound is given on neglected theta^4 or higher-harmonic terms, so this is the main modeling burden.
  • domain assumption The light-fluid potential phi_L in (7) may be singular at r=0 because only its local form near the bubble tip is used.
    Explicitly noted by the authors in Section II.B; the global unphysical divergence is ignored.
  • domain assumption Gravity is uniform and radial with potential Psi = +- g r.
    Used in the Bernoulli boundary condition (5); constant g is a simplification for spherical RTI.
  • domain assumption In the late-time limit, [r0+eta0(0)]/[r0+eta0(t)] tends to 0 for CG and infinity for DG.
    Stated in Section II.C before equation (15); in the DG case the bubble approaches r=0, where the local expansion may eventually break down.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Nonlinear asymptotic bubble growth in single-mode spherical Rayleigh-Taylor instability." pith.science (2026). https://pith.science/paper/WDN2SM2S

@misc{pith2026260729260,
  author       = {Pith},
  title        = {Pith review of: Nonlinear asymptotic bubble growth in single-mode spherical Rayleigh-Taylor instability},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/WDN2SM2S}},
  note         = {Machine review of arXiv:2607.29260}
}
read the original abstract

We present an analytical model for the nonlinear growth of a single-mode Rayleigh-Taylor instability (RTI) bubble in spherical geometry. The model captures the bubble growth along the polar axis, spanning the linear to nonlinear regimes, for arbitrary Atwood numbers and under both converging- and diverging-gravity configurations. The model predicts that the bubble acceleration approaches an asymptotic value in the nonlinear stage. The spherical geometry is found to enhance the RTI bubble growth relative to planar and cylindrical configurations with the same effective perturbation wavenumber in the converging-gravity cases, whereas it mitigates the bubble growth in the diverging-gravity cases. The model predictions show favorable agreement with direct numerical simulations.

Figures

Figures reproduced from arXiv: 2607.29260 by the authors.

Figure 1
Figure 1. FIG. 1. Schematic of a single-mode RTI in spherical geometry. (a) Initial single-mode perturbation [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. Ratio of the asymptotic bubble accelerations with [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3. Density contours from the simulations with [PITH_FULL_IMAGE:figures/full_fig_p010_3.png] view at source ↗
Figures from the paper (2 more)
Figure 4
Figure 4. Figure 4: FIG. 4. Radial trajectory of the bubble tip in the CG ( [PITH_FULL_IMAGE:figures/full_fig_p011_4.png]
Figure 5
Figure 5. Figure 5: FIG. 5. Evolution of [PITH_FULL_IMAGE:figures/full_fig_p012_5.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

23 extracted references

  1. [1]

    I , author=

    The Instability of Liquid Surfaces when Accelerated in a Direction Perpendicular to their Planes. I , author=. Proceedings of the Royal Society of London , volume=

  2. [2]

    1961 , publisher=

    Hydrodynamic and hydromagnetic stability , author=. 1961 , publisher=

    1961

  3. [3]

    Plesset, M. S. , title =. Journal of Applied Physics , volume =. 1954 , month =

    1954

  4. [4]

    and Hurricane, Omar A

    Zhou, Ye and Sadler, James D. and Hurricane, Omar A. Instabilities and Mixing in Inertial Confinement Fusion. Annual Review of Fluid Mechanics. 2025

    2025

  5. [5]

    2005 , volume =

    Drake, R P , title =. 2005 , volume =

    2005

  6. [6]

    Inogamov, N. A. , journal=

  7. [7]

    2003 , publisher=

    Sohn, Sung-Ik , journal=. 2003 , publisher=

    2003

  8. [8]

    2024 , publisher=

    Fluid mechanics , author=. 2024 , publisher=

    2024

Show all 23 references
  1. [9]

    , title = "

    Layzer, D. , title = ". The Astrophysical Journal , year = 1955, volume =

    1955

  2. [10]

    Goncharov, V. N. , journal =

  3. [11]

    Journal of Fluid Mechanics , volume =

    Drop fragmentation by laser-pulse impact , author =. Journal of Fluid Mechanics , volume =. 2020 , publisher=

    2020

  4. [12]

    Qian, Ke-Jian and Li, Zhu-Jun and Tao, Tao and Zhang, De-Hua and Yan, Rui and Ding, Hang , journal =

  5. [13]

    and Wang, W

    Zhang, J. and Wang, W. M. and Yang, X. H. and Wu, D. and Ma, Y. Y. and Jiao, J. L. and Zhang, Z. and Wu, F. Y. and Yuan, X. H. and Li, Y. T. and Zhu, J. Q. , title =. Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences , volume =

  6. [14]

    2020 , pages=

    Journal of Fluid Mechanics , author=. 2020 , pages=

    2020

  7. [15]

    , journal =

    Mikaelian, Karnig O. , journal =

  8. [16]

    Physics of Plasmas , volume =

    Li, Xilai and Wu, Yilin and Chen, Zhengnuo and Yang, Mengqi and Zhang, Jie , title =. Physics of Plasmas , volume =

  9. [17]

    and Wang, L

    Zhang, J. and Wang, L. F. and Ye, W. H. and Wu, J. F. and Guo, H. Y. and Zhang, W. Y. and He, X. T. , title =. Physics of Plasmas , volume =

  10. [18]

    and Wang, L

    Zhang, J. and Wang, L. F. and Ye, W. H. and Wu, J. F. and Guo, H. Y. and Ding, Y. K. and Zhang, W. Y. and He, X. T. , title =. Physics of Plasmas , volume =

  11. [19]

    and Wang, L

    Zhang, J. and Wang, L. F. and Wu, J. F. and Ye, W. H. and Zou, S. Y. and Ding, Y. K. and Zhang, W. Y. and He, X. T. , title =. Physics of Plasmas , volume =

  12. [20]

    2018 , volume =

    Guo, Hong-Yu and Wang, Li-Feng and Ye, Wen-Hua and Wu, Jun-Feng and Li, Ying-Jun and Zhang, Wei-Yan , title =. 2018 , volume =

    2018

  13. [21]

    Goncharov, V. N. and Li, D. , journal =

  14. [22]

    1996 , author =

    Journal of Computational Physics , volume =. 1996 , author =

    1996

  15. [23]

    and Spruce, Mark and Speares, William , journal=

    Toro, Eleuterio F. and Spruce, Mark and Speares, William , journal=. 1994 , publisher=

    1994

Pith tools

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