Pith. sign in

REVIEW 3 major objections 5 minor 13 references

Laminar gaps mirror turbulent puffs in pipe flow

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

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

Pith's one-line read Turbulent pipe flow contains stable laminar gaps, called anti-puffs, that disappear near a critical Reynolds number.

desk verdict New DNS shows laminar gaps in pipe flow are coherent structures with a preferred width and exponential lifetimes; the self-tuning stability claim is plausible but not directly demonstrated, and the Re~2900 onset is a conjecture. read the letter →

arxiv 2608.04823 v1 pith:SMIRAHJU submitted 2026-08-05 physics.flu-dyn cond-mat.stat-mechnlin.CD

Shai Kapon , Tobias Grafke , Anna Frishman This is my paper · ORCID
classification physics.flu-dyncond-mat.stat-mechnlin.CD
keywords pipeflowlaminar-turbulentintermittencyanti-puffsturbulentpuffsslugfrontspeedsself-tuningcriticalReynoldsnumberdirectnumericalsimulation
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 argues that the laminar holes observed inside turbulent pipe flow for Reynolds numbers 2400 to 2550 are not leftover fragments of an expanding slug, but a distinct dynamical state, the anti-puff. An anti-puff persists because its two laminar-turbulent fronts match speeds: the downstream strong front produces a blunted mean flow profile that slows turbulence production enough to keep pace with the upstream weak front. The result is a traveling laminar pocket of roughly constant width embedded in turbulence, whose lifetime is exponentially distributed and shortens as the Reynolds number rises. The authors propose that anti-puffs lose stability near Re≈2900, making the onset of homogeneous turbulence a sharp transition rather than a gradual crossover. A sympathetic reader cares because this fills the least understood segment of the route to fully developed turbulence in pipe flow.

What carries the argument

The central object is the anti-puff: a laminar gap bounded downstream by a strong front with a blunted mean flow profile and upstream by a weak front adjacent to streamwise-homogeneous turbulence. The carrying mechanism is front-speed matching: writing $c_+(\mathrm{Re},u_{\mathrm{front}})$ for the strong-front speed and $c_-(\mathrm{Re},u_{\mathrm{turbulence}})$ for the weak-front speed, a too-small gap gives a flatter downstream profile, slower turbulence production, and expansion, while a too-large gap gives a less blunted profile and contraction; the stable width is where $c_+ - c_- \approx 0$. The exponential lifetime statistics identify closure as escape from the chaotic saddle surrounding the anti-puff state.

What would settle it

Track laminar gaps in direct numerical simulation or experiment at Re between 2900 and 3000 for times long compared with the measured lifetimes of roughly 40–90D/Ū: observation of near-constant-width gaps with c_+ - c_- ≈ 0 would refute the sharp disappearance, as would direct measurement inside a real gap showing c_+ ≥ c_- at these Reynolds numbers.

Watch

Extended reading notes

Core claim

Direct numerical simulations at Re = 2400, 2450, 2500 and 2550 show that the centreline velocity at the strong front of a laminar gap is bimodally distributed: one peak sits at the full Hagen–Poiseuille value 2Ū and belongs to slugs, while a second, Reynolds-dependent peak at lower values belongs to a new type of strong front, the anti-puff front. Histograms of front-speed difference versus gap width show small gaps clustered at $c_+ - c_- \approx 0$ with a negative slope, meaning size fluctuations are corrected by a self-tuning feedback loop. Mean profiles show the strong-front centreline velocity becoming more blunted as Re increases, matching the slowing of relaminarization at the weak front. Survival probabilities of anti-puffs are exponential, with mean lifetimes decreasing from roughly 90D/Ū at Re = 2400 to roughly 40D/Ū at Re = 2550, indicating memoryless closure by escape from a chaotic saddle. The paper also argues that above Re≈2900 no speed matching is possible, so laminar gaps always contract and statistically homogeneous turbulence sets in.

Load-bearing premise

The predicted sharp loss of stability at Re≈2900 rests on assuming that the speed of the weak front of an anti-puff can be read off from the measured speed of the downstream front of a slug; if that transfer between geometries fails, no simulation in the paper directly constrains the critical Reynolds number.

Editorial extensions

If this is right

  • Anti-puffs, not slugs, are the building blocks of intermittency for Re≳2300, so the turbulence fraction is set by their nucleation and decay rates.
  • Gap closing is a Poissonian process, so the mean lifetime is the only parameter needed to describe anti-puff death.
  • The anti-puff strong front is a genuinely new front type whose blunting increases with Reynolds number, linking mean flow profile to front speed.
  • The transition to homogeneous turbulence near Re≈2900 is sharp: above it, laminar gaps contract on average and no stable anti-puff exists.
  • At Re≈2300 the distinction between puffs, slugs and anti-puffs dissolves, with a jammed transitional range 2250<Re<2350 that blurs the boundary.

Reading between the lines

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

  • The paper leaves implicit that if gap closing is escape from a chaotic saddle, anti-puff lifetime should scale with the distance to the gap edge state; this could be tested by computing the edge state at several Re and comparing with the measured τ.
  • A testable extension is that anti-puff stability depends only on the centreline velocity at the strong front, so an experiment measuring centreline velocity and front speeds along the pipe could confirm the self-tuning mechanism without full three-dimensional fields.
  • The same front-speed-matching logic may transfer to plane Couette flow, where laminar gaps with exponential lifetimes are seen; comparing the c_+ and c_- curves across geometries would tell whether anti-puffs are a universal feature of subcritical wall-bounded transition.
  • If the sharp onset at Re≈2900 is right, laminar gaps should be entirely absent in long-time simulations or experiments at Re=3000, including after strong finite-amplitude perturbations.
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 / 5 minor

Summary. The manuscript presents direct numerical simulations of pipe flow at Re=2400, 2450, 2500, and 2550 and argues that laminar gaps embedded in otherwise turbulent flow form a distinct dynamical state, termed 'anti-puffs.' These gaps are proposed to be the mirror image of puffs: a laminar pocket of roughly constant width, stabilized by a self-tuning matching between the speed of a weak upstream front and a strong downstream front whose mean flow profile is blunted relative to the laminar profile. The authors characterize the front velocity distributions, the width-dependent front-speed difference, the mean spatial profiles, and the lifetime statistics, reporting exponentially distributed lifetimes with a mean that decreases with Reynolds number. They further suggest that anti-puffs lose stability above Re≈2900, marking the onset of statistically homogeneous turbulence.

Significance. If the anti-puff picture is correct, it would provide a unified description of the intermittent laminar-turbulent phase in pipe flow and would identify a new dynamical state with a self-regulation mechanism reminiscent of puffs. The paper's strengths include the use of open-source DNS, a transparent front-detection procedure with a stated threshold, explicit reporting of the front velocity and width histograms, and quantitative lifetime fits. The exponential lifetime distributions and the Reynolds-number trends are valuable and reproducible ingredients. However, the central dynamical claim—self-tuned stability of anti-puffs—is not directly demonstrated, and the sharp-onset prediction at Re≈2900 rests on an assumption that fronts of slugs and gaps are interchangeable. These issues are load-bearing for the paper's main conclusions.

major comments (3)
  1. [Sec. 2, Fig. 5] The negative slope of the conditional mean c(w)=c+−c− and its zero crossing are presented as evidence for the self-tuning feedback, but for a stationary stochastic width process the stationary width distribution and the conditional drift are linked by the Fokker-Planck equation (for constant noise, d(w)∝d ln ρ/dw). A peak in the width histogram therefore implies a negative-slope zero crossing of d(w) regardless of whether gaps are dynamically stable. Because gaps are selected for T>30D/U and front velocities are averaged over the same five-snapshot tracks used to define widths, the negative slope in Fig. 5 may be an artifact of conditioning rather than evidence of a restoring mechanism. The paper itself concedes in Sec. 2 that the self-tuning mechanism is not directly demonstrated, yet Sec. 5 concludes that anti-puffs 'remain dynamically stable due to a self-regulated matching of speeds.' This is an overstatement; the authors should either soften the conclusion or provide a direct test, e.g., showing that width perturbations relax to the preferred width in a Lagrangian frame, or comparing the measured c(w) with that expected from a stationary process with the observed width distribution.
  2. [Sec. 4] The prediction that anti-puffs disappear at Re≈2900 rests on the sentence 'We can deduce the speed of the upstream front of an anti-puff from measurements of the speed of the downstream front of a slug (as the two are equivalent).' This equivalence is an assumption that is not tested or referenced with a direct comparison in the present DNS. If the two fronts are not equivalent (e.g., because the mean flow profile or the turbulence intensity upstream of a gap differs from that downstream of a slug), no simulation in this paper constrains the existence or absence of anti-puffs above Re=2550. I ask the authors to either validate this transfer using their own data at 2400–2550 (e.g., comparing measured upstream gap front speeds with the downstream slug front speeds at the same Re) or to reframe the Re≈2900 statement as a conjecture rather than a supported scenario.
  3. [Sec. 2/3] The operational definition of an anti-puff as a laminar gap smaller than 20D for at least 80% of its lifetime, with lifetimes restricted to T>30D/U and with the first and final three snapshots removed, introduces conditioning that may bias the statistics. In particular, the peak in the width histogram at c≈0 in Fig. 5 could simply reflect the fact that only long-lived gaps are included, while transient gaps (which would show larger |c|) are excluded. The authors report ~100–400 gaps per Reynolds number but do not report how the mean lifetimes or the exponential fits change when the size cutoff or the minimum-lifetime threshold is varied. Please provide robustness checks and, if possible, statistics over all detected gaps (including shorter-lived ones) to substantiate the claim that 'most laminar gaps' are anti-puffs.
minor comments (5)
  1. [Sec. 4] The symbol c is used for the turbulence advection speed in Sec. 4 after being defined as c=c+−c− in Sec. 2; this notation clash should be fixed to avoid confusion.
  2. [Fig. 6(d)] The survival probability fits start at T0=50D/U, but the text does not state how many gaps contribute to each fit or how the fitted lifetimes and the goodness of fit vary with T0. Reporting this would strengthen the exponential-lifetime claim.
  3. [Fig. 6 caption] The shaded region is described only as 'variance'; please specify whether this is the standard error, standard deviation, or another quantile measure.
  4. [Sec. 2] The sentence 'Increasing axial resolution does not change the results' is stated without a supporting figure or quantitative check; a brief resolution-convergence statement or supplementary figure would be useful.
  5. [Sec. 2] A quantitative threshold-sensitivity test for q_th (e.g., varying q_th over 3×10^-2 to 10^-1) would support the claim that the front statistics are robust to the detection threshold.

Circularity Check

0 steps flagged · score 2.0 of 10

No significant circularity: the anti-puff concept is carried over from the authors' own Frishman and Grafke (2022) paper, but the present DNS evidence, exponential lifetime fits, and the Re≈2900 argument rest on independent simulation data and external measurements.

full rationale

The derivation chain is not circular. The paper introduces 'anti-puffs' via a self-citation to Frishman and Grafke (2022), but that citation is used to name and motivate the hypothesis, not to establish the present results; the existence and stability of the structures are then tested with new DNS at Re=2400–2550, with the front-velocity and centreline-velocity statistics shown in Figures 4 and 5. The exponential lifetime distribution is fitted to survival data after conditioning on T0=50 D/U; this is a statistical characterization, not a quantity predicted from the fit, and the fitted lifetime tau is not used to define an anti-puff. The Re≈2900 onset scenario is based on external measurements by Song et al. (2017) of the downstream slug-front speed, not on parameters fitted here; its validity depends on the explicitly stated assumption that an anti-puff's upstream front is equivalent to a slug's downstream front, so a failure of that assumption would be a correctness risk, not a circular reduction. The paper itself flags, in Section 2, that the observations 'do not explain, nor directly demonstrate, the self-tuning mechanism which keeps u_front≈u_anti-puff on average at the downstream front'; this is an acknowledged limitation in evidence, not a step that equates output with input. The negative-slope evidence for self-tuning could plausibly be affected by conditioning on long-lived gaps, but that is a statistical-validity concern rather than circularity by construction. No equation is shown to reduce to its own definition, and no fitted parameter is renamed as a prediction. The only self-citation is non-load-bearing, so the circularity score is 2.

Assumptions & free parameters 3 free parameters · 4 assumptions · 2 invented entities

The central claims rest on the standard Navier-Stokes model for pipe flow, a set of front-speed modeling assumptions adopted from prior work, and an operational classification of anti-puffs. No physically new free constants are fitted, but several analysis thresholds and cutoffs shape the reported statistics.

free parameters (3)
  • Turbulence threshold q_th = 5 x 10^-2
    Used to detect laminar-turbulent fronts; authors state varying it yields qualitatively identical results, but the value is chosen by hand.
  • Anti-puff size cutoff = 20D for at least 80% of lifetime
    Operational classification of small gaps as anti-puffs; no physical derivation of this cutoff is given.
  • Exponential fit start T0 = 50 D/Ubar
    Starting point for the survival probability fit; affects the reported mean lifetimes, as the authors note threshold sensitivity.
assumptions (4)
  • domain assumption Front speeds c+(Re,u_front) and c-(Re,u_turbulence) are well-defined functions of Re and centreline velocity alone.
    Invoked in Sec. 2 to define the speed-matching condition that stabilizes both puffs and anti-puffs.
  • domain assumption Turbulence in the bulk is advected at a constant speed c in the lab reference frame.
    Used in Sec. 4 to argue c+ < c < c- for Re>2900, based on Song et al. (2017).
  • domain assumption The upstream front of an anti-puff is equivalent to the downstream front of a slug, so its speed can be taken from slug measurements.
    Central to the Re≈2900 instability argument in Sec. 4; not independently measured in this paper.
  • domain assumption Exponentially distributed gap lifetimes imply escape from a chaotic saddle.
    Analogous to puff decay arguments (Hof et al., 2006); used to interpret the lifetime statistics in Sec. 3.
invented entities (2)
  • Anti-puff independent evidence
    purpose: Named dynamical state describing persistent laminar gaps in turbulent pipe flow, with self-tuned front speeds.
    Directly observed in DNS as gaps with a preferred width and bimodal strong-front velocity; the entity is a classification, not a new force or field.
  • Gap edge state
    purpose: Hypothesized critical laminar nucleus on the basin boundary between turbulent and anti-puff states, mediating nucleation and decay.
    Proposed in Sec. 3 as a qualitative explanation; no direct observation or separate falsifiable handle is provided.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Laminar gaps mirror turbulent puffs in pipe flow." pith.science (2026). https://pith.science/paper/SMIRAHJU

@misc{pith2026260804823,
  author       = {Pith},
  title        = {Pith review of: Laminar gaps mirror turbulent puffs in pipe flow},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/SMIRAHJU}},
  note         = {Machine review of arXiv:2608.04823}
}
abstract

Pipe flow at intermediate Reynolds numbers, between the laminar and fully turbulent regimes, takes the form of several spatially and temporally intermittent phases in which turbulent and laminar states coexist. In the lower range, $Re\in (1750,2300)$, turbulence appears in the form of localized traveling structures called "puffs", which form long-lived chaotic dynamical states, whose stochastic decays and splits control the steady state intermittency. At the other end, $Re\in (2300,3000)$, puffs are replaced by an extended turbulent state, with laminar pockets intermittently forming and disappearing within it. Using direct numerical simulations of pipe flow at $Re = 2400, 2450, 2500, 2550$, we provide evidence that these laminar gaps form a distinct dynamical state analogous to puffs: a traveling laminar pocket in a turbulent surrounding, stabilized by a shear-dependent self-tuning mechanism. We analyse the mean spatial profile of these gaps and show that their lifetimes are exponentially distributed, suggesting that gap closing corresponds to an escape from a chaotic saddle. Finally, we suggest these laminar gaps become unstable and disappear at a finite Reynolds number, $Re\sim 2900$, which can be interpreted as the onset point of spatially and temporally homogeneous turbulence.

Figures

Figures reproduced from arXiv: 2608.04823 by the authors.

Figure 1
Figure 1. Space–time visualization of the flow for four different Reynolds numbers using a colour map of the [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Schematic illustration of (a) a puff/slug and (b) an anti-puff. The upper panel shows the streamwise [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Comparison between an anti-puff (a) and a laminar complement of a slug (b) at [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figures from the paper (3 more)
Figure 4
Figure 4. Figure 4: Distributions of the centreline velocity at laminar–turbulent fronts, [PITH_FULL_IMAGE:figures/full_fig_p005_4.png]
Figure 5
Figure 5. Figure 5: Two-dimensional histograms of the front-velocity difference [PITH_FULL_IMAGE:figures/full_fig_p006_5.png]
Figure 6
Figure 6. Figure 6: Mean structure and lifetime statistics of anti-puffs for [PITH_FULL_IMAGE:figures/full_fig_p007_6.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

13 extracted references · 12 canonical work pages

  1. [13]

    doi: 10.1017/jfm.2023.1015

    ISSN 0022- 1120, 1469-7645. doi: 10.1017/jfm.2023.1015. URLhttps://www.cambridge.org/core/product/ identifier/S0022112023010157/type/journal_article. I Wygnanski, Mo Sokolov, and D Friedman. On transition in a pipe. part

    work page arXiv 2023

  2. [1883]

    doi: 10.1098/ rstl.1883.0029

    ISSN 0261-0523. doi: 10.1098/ rstl.1883.0029. URLhttps://doi.org/10.1098/rstl.1883.0029. J Rotta. An experimental contribution to the transition from laminar to turbulent flow in a pipe.IXe Congr. intern. Mecan. Appl., 3,

    work page arXiv

  3. [1986]

    doi: 10.1017/S0022112086002379

    ISSN 0022-1120, 1469-7645. doi: 10.1017/S0022112086002379. URLhttps://www.cambridge.org/core/product/identifier/ S0022112086002379/type/journal_article. Dwight Barkley. Simplifying the complexity of pipe flow.Physical Review E—Statistical, Nonlinear, and Soft Matter Physics, 84(1):016309,

    work page doi:10.1017/s0022112086002379

  4. [2008]

    doi: 10.1017/S0022112008003315

    ISSN 0022-1120, 1469-7645. doi: 10.1017/S0022112008003315. URLhttps://www.cambridge.org/core/product/ identifier/S0022112008003315/type/journal_article. Yves Pomeau. Front motion, metastability and subcritical bifurcations in hydrodynamics.Physica D: Nonlinear Phenomena, 23(1-3):3–11,

    work page doi:10.1017/s0022112008003315

  5. [2010]

    URLhttps://www.science.org/doi/10.1126/science.1186091

    doi: 10.1126/science.1186091. URLhttps://www.science.org/doi/10.1126/science.1186091. Pavan V. Kashyap, Yohann Duguet, and Olivier Dauchot. Linear Instability of Turbulent Channel Flow. Physical Review Letters, 129(24):244501, December

    work page doi:10.1126/science.1186091

  6. [2011]

    doi: 10.1063/1.3580263

    ISSN 1070-6631. doi: 10.1063/1.3580263. URLhttps: //doi.org/10.1063/1.3580263. Laurette S. Tuckerman, Tobias Kreilos, Hecke Schrobsdorff, Tobias M. Schneider, and John F. Gibson. Turbulent-laminar patterns in plane Poiseuille flow.Physics of Fluids, 26(11):114103, November

    work page doi:10.1063/1.3580263

  7. [2014]

    doi: 10.1063/1.4900874

    ISSN 1070-6631, 1089-7666. doi: 10.1063/1.4900874. URLhttps://pubs.aip.org/pof/article/ 26/11/114103/315350/Turbulent-laminar-patterns-in-plane-Poiseuille. Ashley P Willis. The Openpipeflow Navier–Stokes solver.SoftwareX, 6:124–127,

    work page doi:10.1063/1.4900874

  8. [2015]

    Duguet, P

    Y. Duguet, P. Schlatter, and D. S. Henningson. Formation of turbulent patterns near the onset of transi- tion in plane Couette flow.Journal of Fluid Mechanics, 650:119–129, May 2010a. ISSN 0022-1120, 1469-7645. doi: 10.1017/S0022112010000297. URLhttps://www.cambridge.org/core/product/ identifier/S0022112010000297/type/journal_article. Y. Duguet, A. P. Wil...

    work page doi:10.1017/s0022112010000297

Show all 13 references
  1. [2017]

    URLhttp://arxiv.org/abs/1705.03838

    doi: https: //doi.org/10.1016/j.softx.2017.05.003. URLhttp://arxiv.org/abs/1705.03838. Haoyang Wu and Baofang Song. Scaling and mechanism of the propagation speed of the upstream turbulent front in pipe flow.Journal of Fluid Mechanics, 977:R4, December

    2017 arXiv

  2. [2018]

    doi: 10.1017/jfm.2017.923

    ISSN 0022-1120, 1469-7645. doi: 10.1017/jfm.2017.923. URLhttps://www.cambridge.org/core/journals/journal-of-fluid-mechanics/ article/abs/critical-point-of-the-transition-to-turbulence-in-pipe-flow/ 2FD81552A7908EE1E695EBF2005DC9E5. Mina Nishi, B¨ ulent¨Unsal, Franz Durst, and ...

    2017 doi

  3. [2019]

    doi: 10.1103/PhysRevFluids.4.113903

    ISSN 2469-990X. doi: 10.1103/PhysRevFluids.4.113903. URL https://link.aps.org/doi/10.1103/PhysRevFluids.4.113903. Joseph D Skufca, James A Yorke, and Bruno Eckhardt. Edge of chaos in a parallel shear flow.Physical review letters, 96(17):174101,

    doi

  4. [2022]

    URL https://link.aps.org/doi/10.1103/PhysRevLett.129.244501

    doi: 10.1103/PhysRevLett.129.244501. URL https://link.aps.org/doi/10.1103/PhysRevLett.129.244501. Gr´egoire Lemoult, Vasudevan Mukund, Hong-Yan Shih, Gaute Linga, Joachim Mathiesen, Nigel Goldenfeld, and Bj ¨orn Hof. Directed percolation and puff jamming near the transition to...

    doi

  5. [2023]

    doi: 10.1017/jfm.2023.289

    ISSN 0022- 1120, 1469-7645. doi: 10.1017/jfm.2023.289. URLhttps://www.cambridge.org/core/product/ identifier/S0022112023002896/type/journal_article. Bj¨orn Hof, Jerry Westerweel, Tobias M Schneider, and Bruno Eckhardt. Finite lifetime of turbulence in shear flows.Nature, 443(7...

    2023 doi

Pith tools

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