Pith. sign in

REVIEW 5 major objections 5 minor 15 references

Variational determination of minimal absorbing zones in incompressible shear flows

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

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

Pith's one-line read For plane Couette and Poiseuille flow, the paper computes the minimal absorbing zone: a state-space sphere that every trajectory must eventually enter and never leave.

desk verdict Novel variational method for minimal absorbing zones in shear flows, with an honest negative result on mean-flow prediction, but the rigorous-radius claim is undercut by the one-dimensional restriction. read the letter →

arxiv 2507.18359 v1 pith:ZBSC2CBO submitted 2025-07-24 physics.flu-dyn

Péter Tamás Nagy This is my paper · ORCID
classification physics.flu-dyn MSC 76D0576E30
keywords absorbingzoneReynolds-OrridentityshearflowstabilityvariationaloptimisationplaneCouettePoiseuilleturbulentmeanglobal
topics Navier-Stokes Regularity
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 sets out to show that incompressible shear flows possess a computable absorbing zone: a hypersphere in state space, centred on an arbitrary 'shift flow', such that every trajectory starting outside the sphere eventually enters it and then stays inside. Because the zone attracts all trajectories, it necessarily contains every attractor, unstable periodic orbit, and edge state of the flow. Existence rests on the Reynolds–Orr identity, which removes the nonlinear advection terms from the kinetic-energy equation and makes energy a Lyapunov-like functional everywhere outside the zone. The paper turns this existence statement into a numerical method, using gradient-based optimisation to shrink the zone to minimal radius for plane Couette and plane Poiseuille flows. The hoped-for bonus, that the zone's centre approximates the turbulent mean flow, is tested against direct numerical simulation data and found quantitatively wanting, but the minimal zone itself survives as a rigorous attractor-containing set and a tool for global-stability arguments.

What carries the argument

The load-bearing identity is the Reynolds–Orr kinetic-energy balance: after decomposing the velocity into a stationary shift flow and fluctuations, the nonlinear self-interaction of the fluctuations makes no net contribution to $de/dt$, so outside a sufficiently large sphere the energy derivative is dominated by negative viscous dissipation. Existence of the zone is governed by an infinite-perturbation eigenvalue problem whose largest eigenvalue must be negative; the zone boundary is then set by the largest energy level at which the constrained Euler–Lagrange equations admit a solution with $de/dt=0$, the critical finite perturbation. Minimality is achieved by gradient-based optimisation, where the sensitivity of the zone radius to the shift flow is computed by solving an adjoint linear system; in one-dimensional wall-normal profiles the eigenvalue problem is simplified by a Fourier ansatz in the two homogeneous directions.

What would settle it

For one of the reported optimal shift flows, say Poiseuille flow at $Re=2000$, solve the infinite-perturbation eigenvalue problem on a grid of nonzero streamwise and spanwise wavenumbers. If any oblique wavenumber pair yields a growth rate larger than the $k_1=k_3=0$ value used in the paper, the minimal zone is not absorbing. A complementary dynamical test is to initialise a resolution-converged simulation just outside the claimed zone and check that the kinetic energy never increases above its initial value before entering; one crossing would contradict the zone's definition.

Watch

Extended reading notes

Core claim

The central claim is that, for any shift flow whose largest infinite-perturbation growth rate is negative, there is a critical energy level at which the maximal kinetic-energy growth rate changes sign; the energy sphere at that level is an absorbing zone. The boundary is found by solving a constrained Euler–Lagrange problem with a prescribed energy level, and the true boundary is the maximal energy at which the energy derivative vanishes. Varying the shift flow under the constraint that the infinite-perturbation eigenvalue remain negative, the paper minimises the zone radius and obtains explicit minimal absorbing zones for Couette and Poiseuille flow as functions of Reynolds number. On the paper's own evidence, the optimal shift-flow profiles are not quantitatively accurate turbulent mean profiles: they have sharper near-wall gradients and, in viscous units, sit well above the law of the wall. What the computation does deliver is a finite, explicit, provable set that contains all attractors, plus an absorbing zone around the laminar state whenever the laminar state lies inside the minimal zone.

Load-bearing premise

The load-bearing assumption is that for every optimized shift flow the fastest-growing infinite perturbation has no variation in the flow direction and the critical finite perturbation is uniform in the spanwise direction; no proof is given for arbitrary optimized shift flows, and if an oblique perturbation grew faster, the computed zone radius would be too small.

Editorial extensions

If this is right

  • Below the classical energy-stability limits ($Re_E=20.6625$ for Couette and $Re_E=49.6035$ for Poiseuille), the minimal absorbing zone collapses to the laminar point, recovering the classical result that no other attractor can exist.
  • Above those limits the zone has finite, computable radius comparable to the laminar energy, so all turbulence-related attractors are guaranteed to live inside an explicitly known finite-energy set.
  • Any larger set containing the minimal zone is also absorbing; in particular one can build an absorbing zone around the laminar base flow, and if that zone sits inside the laminar region of attraction, global stability of the base flow is proven.
  • At high Reynolds number the optimal shift-flow profiles become nearly Reynolds-number-independent and resemble dissipation-maximising variational profiles, but they do not match measured turbulent mean profiles; the paper attributes the mismatch to the use of plain kinetic energy and to non-uniform exploration of state space by turbulent trajectories.
  • Directly minimising the laminar-centred absorbing zone rather than the shift-flow zone gives about 5 percent smaller guaranteed bounds for Poiseuille flow, while the zone around the shift flow grows by about 10 percent.

Reading between the lines

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

  • If the minimal zone radius is read as a lower bound for energy-driven transition, any perturbation with energy below that radius that still becomes turbulent must rely on non-normal transient growth; comparing the zone radius with known minimal seed energies would calibrate how much of transition is linear transient growth rather than direct energy growth.
  • The zone's failure to track the turbulent mean may be partly a metric artefact: replacing kinetic energy with a weighted norm that still cancels nonlinear terms would yield a different, possibly tighter zone, and the paper's own discussion invites this test.
  • Applying the same variational construction to pipe flow or boundary layers would show whether the two patterns seen here, near-wall steepening and a plateau below the law of the wall, are universal properties of energy-based absorbing zones or artefacts of plane-channel geometry.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

5 major / 5 minor

Summary. The paper introduces a variational method for computing 'minimal absorbing zones' in the state space of incompressible shear flows. A shift flow U is taken as the centre of a ball measured by perturbation kinetic energy; using the Reynolds–Orr identity, the author shows that if the largest eigenvalue of the infinite-amplitude variational problem is negative, the perturbation energy decays outside a sufficiently large ball. The finite-amplitude Euler–Lagrange equations are then solved to find the critical energy e_AZ at which the maximum growth rate vanishes. A gradient-based optimisation over one-dimensional shift flows is performed for plane Couette and Poiseuille flows, and the resulting optimised profiles are compared with DNS mean profiles and the Spalding law; the comparison is unfavourable, leading the author to reject the random-wandering hypothesis for the turbulent mean. The paper argues that an absorbing zone around the laminar base flow can be constructed and used for global-stability proofs.

Significance. If the computed zones are indeed absorbing, the method provides a rigorous, fully nonlinear containment set for all attractors in wall-bounded shear flows, complementing edge-state and unstable-periodic-orbit analyses. A clear strength is that no parameters are fitted to turbulent-mean data: the shift flow is optimised against the zone-size objective, and the DNS/Spalding comparison is a post-hoc benchmark. The explicit test and rejection of the turbulent-mean hypothesis is scientifically honest and useful. However, the numerical value of e_AZ and the base-flow-zone construction rest on several unproved technical assumptions that must be resolved before the central claims can be accepted.

major comments (5)
  1. [§3, Eq. (3.5)] The reduction of the finite-amplitude problem to spanwise-uniform perturbations (k1=k3=0) is not rigorously justified. The claim that a non-oscillatory right-hand side forces all oscillatory Fourier components to vanish is only valid when μ_L is not an eigenvalue of the operator in (2.8) for any oblique wavenumber pair. At resonance, homogeneous oblique solutions can be added to the particular solution while preserving the Euler–Lagrange equations, and such branches may yield a larger critical energy. The statement that the maximum infinite-amplitude growth rate is 'consistently found' at k1=0 does not settle the finite-energy zero crossing. The author should either prove that the critical μ_L lies outside the spectrum for all k1,k3 or perform a full wavenumber search over the finite-amplitude branches.
  2. [§4, Eq. (4.1)] The radius of the base-flow-centred absorbing zone is not e_AZ,opt + e_d. In the L2 norm associated with kinetic energy, if a ball of radius sqrt(2e_AZ,opt) around U_SF is contained in a ball of radius R around U_BF, the triangle inequality gives R ≥ sqrt(2e_AZ,opt) + sqrt(2e_d), so the required kinetic-energy radius is (sqrt(e_AZ,opt)+sqrt(e_d))^2, not e_AZ,opt+e_d. The printed formula understates the size of the base-flow zone, and the subsequent conclusion that the minimal zone contains the laminar state is not established by the inequality e_d ≤ e_AZ,BF (which is tautological). This affects the quantitative global-stability discussion and should be corrected.
  3. [§2.2, Figure 2] Selecting the largest-energy stationary solution with dE/dt=0 does not by itself identify the absorbing boundary, because the Euler–Lagrange equations are only necessary conditions for a constrained extremum. The envelope theorem implies that for the global-maximizer branch the derivative dF/dβ equals μ_L; the reported critical value μ_L=+0.024 would imply that the maximum growth rate increases with energy at the crossing, so that F>0 just outside the claimed zone, contradicting the absorbing-zone property. This suggests either that the selected branch is not the global maximizer or that a sign error is present (Eq. (2.10) has −μ_L while Eq. (3.5) has +μ_L). The author should clarify the sign convention and verify directly that max dE/dt < 0 for all β > e_AZ over a fine scan of μ_L and wavenumbers.
  4. [§2.1–§2.2] The proof that trajectories 'eventually enter' the absorbing zone requires more than strict negativity of dE/dt outside the ball. If the maximum growth rate approaches zero as the energy approaches the boundary from above, trajectories may approach the zone asymptotically without entering it in finite time. To establish an absorbing set in the standard sense, one needs a uniform negative bound on dE/dt for energies above some slightly enlarged radius, or an explicit differential inequality guaranteeing finite-time entry. The current argument only gives monotone decrease of energy outside the zone, which is insufficient for the 'eventually enters' claim as stated.
  5. [§2.3] The optimisation is performed with a gradient-based local method (fmincon) initialised from a continuation path; no evidence is given that the computed shift flow is the global minimiser of e_AZ over all admissible shift flows. Since the objective is non-convex, the reported zone should be described as a locally minimal absorbing zone unless a global search or convexity argument is supplied. The abstract's phrase 'minimal radius' is stronger than what is demonstrated.
minor comments (5)
  1. [§3, Eq. (3.5)] There is a sign inconsistency between Eq. (2.10), which contains −μ_L u_i on the left-hand side, and Eq. (3.5), which contains +μ_L u_1. The derivation in the text should be checked and the sign convention stated consistently throughout.
  2. [§4, p. 11] The sentence 'it can be seen that e_d ≤ e_AZ,BF ... therefore, the minimal absorbing zone includes the laminar state' is logically insufficient: containing the laminar state requires e_d ≤ e_AZ,opt, not e_d ≤ e_AZ,BF. The comparison in Figure 7 should be made against e_AZ,opt.
  3. [Figure 2] The axis labels and the sign of μ_L in Figure 2 need clarification, especially because the caption reports a positive critical μ_L while the text and Eq. (3.5) suggest opposite conventions.
  4. [§1 and §4] The reference 'Olivier Dauchot & Paul Manneville 1997' should be formatted consistently with the journal style, and the definition of 'eventually enters' should be made precise.
  5. [§4, Figure 3] The abrupt change in behaviour near Re_E and the statement about sharp velocity profiles at larger Reynolds numbers could be documented more quantitatively (e.g., by reporting the wall-normal grid resolution tests in a table).

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the absorbing-zone radius is obtained by variational optimization and benchmarked externally; the random-wandering hypothesis is tested and rejected.

full rationale

The paper's central derivation is self-contained. The absorbing zone around a shift flow is defined by the Reynolds-Orr energy balance, and the radius e_AZ is obtained by solving the constrained variational problem (2.9)-(2.10) for the critical finite perturbation, with no parameter fitted to the turbulent-mean data that appear later. The shift flow is optimized directly against the zone-size objective (2.12), so the comparison with DNS profiles and the Spalding law is an external benchmark, not an input to the derivation. The auxiliary hypothesis that the zone centroid approximates the turbulent mean flow is an additional interpretive assumption, and the paper explicitly tests it and rejects it as an oversimplification. The only self-citation, Nagy (2025), is invoked for an optional global-stability extension and is not load-bearing for the existence or size of the minimal absorbing zone; moreover, it concerns a finite-dimensional Lyapunov construction rather than the absorbing-zone calculation itself. The restriction to spanwise-uniform critical perturbations is a potential correctness gap, not a circularity, because the radius is not fitted to the quantity it is later claimed to bound. No step in the derivation reduces by construction to its own inputs, so the circularity score is 0.

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

No physical entities are postulated. The shift flow and absorbing zone are mathematical constructions in state space, not new particles, forces, or conserved quantities. There are no fitted constants in the theory; the Lagrange multiplier mu_L is a continuation parameter, not a data-fitted parameter.

assumptions (5)
  • domain assumption Global existence and infinite-time trajectories of the Navier-Stokes solutions in the considered channel flows are assumed.
    The absorbing-zone argument requires every trajectory to evolve for all time so that 'eventually enters' is meaningful. This is standard in fluid-dynamics practice but unproven for 3D Navier-Stokes.
  • domain assumption The shift flow U_i is time-independent, satisfies the boundary conditions, and is divergence-free.
    Section 2 introduces the shift flow as a state-space point, not necessarily a solution. The derivation of eq. (2.5) requires the boundary and divergence properties used in the Gauss divergence theorem steps.
  • ad hoc to paper For the one-dimensional reduction, the critical perturbations are assumed to have k1=k3=0 and the maximum energy growth rate is assumed to occur at k1=0.
    Section 3 states that this is 'consistently found' and cites an analogy to Joseph and Carmi, but no proof is given for arbitrary optimized shift flows. This is load-bearing because the computed zone radii depend on the critical perturbation being spanwise-uniform.
  • ad hoc to paper Gradient-based optimization converges to the global minimum of e_AZ[U].
    Section 2.3 uses fmincon with interior-point method, which only guarantees local minima. The paper calls the result 'minimal' without a global-optimality proof, while noting that multiple absorbing zones may exist.
  • ad hoc to paper The turbulent attractor explores the state space randomly and uniformly, so the centroid of the minimal absorbing zone approximates the turbulent mean flow.
    This is the primary hypothesis of the paper, stated in the abstract and Section 1, and later explicitly identified as likely oversimplified. It is not used to derive the zone but is the premise for the mean-flow prediction.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Variational determination of minimal absorbing zones in incompressible shear flows." pith.science (2026). https://pith.science/paper/ZBSC2CBO

@misc{pith2026250718359,
  author       = {Pith},
  title        = {Pith review of: Variational determination of minimal absorbing zones in incompressible shear flows},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/ZBSC2CBO}},
  note         = {Machine review of arXiv:2507.18359}
}
read the original abstract

The dynamical analysis of shear flows remains challenging, as turbulence generation and evolution are not fully understood. Here, a lesser-explored feature of incompressible shear flows-the absorbing zone-is investigated. This region in the infinite dimensional state space is shown to act as an attractor for all trajectories: any solution initialised outside eventually enters and remains inside. Consequently, the absorbing zone must contain all possible attractors, both chaotic and non-chaotic. Existence is established through the Reynolds-Orr identity, which indicates that nonlinear terms do not directly influence the temporal evolution of kinetic energy. The zone is constructed around a so-called shift flow, and multiple such regions may exist, even when the laminar state is linearly unstable. Gradient based optimisation is employed to identify the absorbing zone of minimal radius. Assuming a chaotic trajectory explores state space randomly, it is hypothesised that the centroid of this minimal zone approximates the turbulent mean flow. This central state is computed for plane Poiseuille and Couette flows and compared with established turbulent mean profiles. Although the proposed hypothesis is found to be an oversimplification-yielding profiles that qualitatively resemble but do not quantitatively reproduce the turbulent means-the methodology provides novel insights into shear flow dynamics. Furthermore, it offers a promising foundation for refining estimates of global stability thresholds. With continued development, this framework may facilitate the direct computation of turbulent mean states without reliance on empirical turbulence models or time-dependent numerical simulations.

Figures

Figures reproduced from arXiv: 2507.18359 by the authors.

Figure 1
Figure 1. The sketch of the minimal absorbing zone around the optimal shift flow [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Time derivative of the kinetic energy (a,b) and energy level (c,d) of the solutions [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. Kinetic energy of the critical finite perturbation, normalised by the kinetic [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗
Figures from the paper (5 more)
Figure 4
Figure 4. Figure 4: Optimal shift-flow profiles (a,c) and their corresponding critical finite [PITH_FULL_IMAGE:figures/full_fig_p010_4.png]
Figure 5
Figure 5. Figure 5: Velocity profiles non-dimensionalised by the friction velocity versus wall [PITH_FULL_IMAGE:figures/full_fig_p011_5.png]
Figure 6
Figure 6. Figure 6: The sketch of the absorbing zones around the optimal shift flow and the laminar [PITH_FULL_IMAGE:figures/full_fig_p011_6.png]
Figure 7
Figure 7. Figure 7: The kinetic energy of critical finite perturbations of the optimal shift flow [PITH_FULL_IMAGE:figures/full_fig_p012_7.png]
Figure 8
Figure 8. Figure 8: Poiseuille flow: comparison of the original optimisation (dashed, pale curves) [PITH_FULL_IMAGE:figures/full_fig_p013_8.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

15 extracted references · 15 canonical work pages

  1. [1]

    2023 Identifying invariant solutions of wall-bounded three- dimensional shear flows using robust adjoint-based variational techniques

    Ashtari, Omid & Schneider, Tobias M. 2023 Identifying invariant solutions of wall-bounded three- dimensional shear flows using robust adjoint-based variational techniques. Journal of Fluid Mechanics 977, A7

    work page 2023

  2. [2]

    Avila, Marc, Barkley, Dwight & Hof, Bj¨orn 2023 Transition to turbulence in pipe flow.Annual Review of Fluid Mechanics 55 (Volume 55, 2023), 575–602

    work page 2023

  3. [3]

    Cavalieri, Andr´e V. G. & Nogueira, Petr ˆonio A. S. 2022 Reduced-order Galerkin models of plane Couette flow. Phys. Rev. Fluids 7, L102601

    work page 2022

  4. [4]

    Doering, Charles R.2009 The 3D Navier-Stokes Problem.Annual Review of Fluid Mechanics41 (Volume 41, 2009), 109–128

    work page 2009

  5. [5]

    2013 Minimal transition thresholds in plane Couette flow

    Duguet, Yohann, Monokrousos, Antonios, Brandt, Luca & Henningson, Dan S. 2013 Minimal transition thresholds in plane Couette flow. Physics of Fluids 25 (8)

    work page 2013

  6. [6]

    Annual Review of Fluid Mechanics 39 (Volume 39, 2007), 447–468

    Eckhardt, Bruno, Schneider, Tobias M., Hof, Bjorn & Westerweel, Jerry2007 Turbulence transition in pipe flow. Annual Review of Fluid Mechanics 39 (Volume 39, 2007), 447–468. Jim´enez, Javier2018 Coherent structures in wall-bounded turbulence.Journal of Fluid Mechanics842, P1

    work page 2007

  7. [7]

    Joseph, D. D. & Carmi, S 1969 Stability of Poiseuille flow in pipes, annuli, and channels. Quarterly of Applied Mathematics 26 (4), 575–599

    work page 1969

  8. [8]

    Annual Review of Fluid Mechanics 44 (Volume 44, 2012), 203–225

    Kawahara, Genta, Uhlmann, Markus & van Veen, Lennaert2012 The significance of simple invariant solutions in turbulent flows. Annual Review of Fluid Mechanics 44 (Volume 44, 2012), 203–225

    work page 2012

Show all 15 references
  1. [9]

    Kerswell, R. R. 1999 Variational principle for the Navier-Stokes equations.Phys. Rev. E 59, 5482–5494

    1999

  2. [10]

    Kerswell, R R, Pringle, C C T & Willis, A P 2014 An optimization approach for analysing nonlinear stability with transition to turbulence in fluids as an exemplar.Reports on Progress in Physics77 (8), 085901

    2014

  3. [11]

    Journal of Fluid Mechanics 1014, A25

    Nagy, P´eter Tam´as 2025 Constructing quadratic Lyapunov functions for conditionally stable fluid dynamics systems. Journal of Fluid Mechanics 1014, A25. Olivier Dauchot & Paul Manneville 1997 Local 𝑣𝑒𝑟𝑠𝑢𝑠 global concepts in hydrodynamic stability theory. J. Phys. II France 7 ...

    2025

  4. [12]

    Plasting, S. C. & Kerswell, R. R. 2003 Improved upper bound on the energy dissipation rate in plane Couette flow: the full solution to Busse’s problem and the Constantin–Doering–Hopf problem with one-dimensional background field. Journal of Fluid Mechanics 477, 363–379

    2003

  5. [13]

    Plasting, S. C. & Kerswell, R. R. 2004 A friction factor bound for transitional pipe flow. Physics of Fluids 17 (1), 011706, arXiv: https://pubs.aip.org/aip/pof/article- pdf/doi/10.1063/1.1828103/14691402/011706 1 online.pdf

    2004 doi

  6. [14]

    T., Willis, Ashley P

    Pringle, Chris C. T., Willis, Ashley P. & Kerswell, Rich R. 2012 Minimal seeds for shear flow turbulence: using nonlinear transient growth to touch the edge of chaos.Journal of Fluid Mechanics 702, 415–443

    2012

  7. [15]

    2024 A canonical Hamiltonian formulation of the Navier–Stokes problem

    Berlinghieri, Joel C. 2024 A canonical Hamiltonian formulation of the Navier–Stokes problem. Journal of Fluid Mechanics 984, A27

    2024

Pith tools

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