pith. sign in

arxiv: 2605.04181 · v1 · submitted 2026-05-05 · 🧮 math.AP

Euler Singularities I: Boundary Blow-Up for Smooth Exact-Odd Axisymmetric Euler with Swirl

Rishad Shahmurov This is my paper

Pith reviewed 2026-05-08 17:39 UTC · model grok-4.3

classification 🧮 math.AP
keywords Euler equationssingularity formationaxisymmetric flowboundary blow-upexact odd symmetryDini comparison systemparametrix estimatesswirl
0
0 comments X

The pith

Smooth axisymmetric-with-swirl Euler data in a cylinder develops finite-time boundary singularity.

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

The paper constructs smooth initial data for the three-dimensional incompressible Euler equations that are axisymmetric with swirl inside a periodic cylinder and obey an exact-odd symmetry, causing the flow to develop a singularity at the side wall in finite time. Exact oddness yields explicit pointwise ODEs at the wall for the relevant vorticity and velocity components. These are closed by a side-wall parametrix estimate for the lifted recovery equation together with a dyadic cluster functional that reduces all errors to a controllable defect, producing a Dini comparison system whose solutions blow up. A reader would care because the result supplies an explicit, symmetry-reduced example of finite-time breakdown for smooth Euler solutions at a boundary, bearing on the long-standing question of whether ideal fluid flows can form singularities.

Core claim

We construct smooth axisymmetric-with-swirl initial data in a periodic cylinder for which the three-dimensional incompressible Euler evolution develops a finite-time boundary singularity. The construction is carried out in the dynamically invariant exact-odd class where Γ(r,−z,t)=−Γ(r,z,t) and G(r,−z,t)=−G(r,z,t), with Γ=ru^θ and G=ω^θ/r. At the side-wall point (r,z)=(1,0), exact oddness gives the pointwise identities ∂_t ∂_z G(1,0,t)=σ(t)∂_z G(1,0,t)+2(∂_z Γ(1,0,t))^2 and ∂_t ∂_z Γ(1,0,t)=σ(t)∂_z Γ(1,0,t) with σ(t)=−∂_z u^z(1,0,t). The proof relies on a side-wall Dirichlet parametrix for the five-dimensional lifted recovery equation −Δ_5 ϕ=G whose leading kernel is K_0(x,y)=C_0 xy/(x^2+y^2^

What carries the argument

The over-compressed dyadic angular cluster functional that folds same-scale fragmentation, growing windows, and exterior tails into an integrably small Campanato defect, paired with the side-wall Dirichlet parametrix for the lifted equation −Δ_5 ϕ=G.

If this is right

  • The wall gradients ∂_z G and ∂_z Γ blow up in finite time while the solution stays smooth away from the wall.
  • The exact-odd symmetry class remains invariant under the evolution.
  • The quadratic self-interaction in the Dini system D^+ A_* ≥ c B_*^2, D^+ B_* ≥ c A_* B_* forces the amplitudes to infinity.
  • The parametrix controls the nonlocal velocity recovery near the boundary with integrable remainders.
  • Blow-up occurs while the flow satisfies the incompressible Euler equations pointwise up to the singular time.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same parametrix-plus-cluster strategy may adapt to other bounded domains or weaker symmetry assumptions to produce interior singularities.
  • High-resolution simulations initialized with the paper's data could measure the precise blow-up rate and test whether the Dini scaling is sharp.
  • The construction isolates the role of swirl in driving boundary stretching, suggesting analogous mechanisms might appear in Navier-Stokes at large Reynolds number.
  • If the error-absorption technique extends, it could yield families of singular solutions parameterized by the initial cluster amplitude.

Load-bearing premise

The side-wall Dirichlet parametrix estimates with controlled remainders, parity-based shear cancellation, and strain bounds on narrow cones, together with the over-compressed dyadic cluster functional, absorb every error term into an integrably small Campanato defect so the Dini comparison system closes.

What would settle it

A numerical simulation of the Euler equations on the constructed initial data in which the wall derivatives ∂_z G(1,0,t) and ∂_z Γ(1,0,t) remain bounded past the predicted blow-up time would falsify the finite-time singularity.

read the original abstract

We construct smooth axisymmetric-with-swirl initial data in a periodic cylinder for which the three-dimensional incompressible Euler evolution develops a finite-time boundary singularity. The construction is carried out in the dynamically invariant exact-odd class \[ \Gamma(r,-z,t)=-\Gamma(r,z,t), \qquad G(r,-z,t)=-G(r,z,t), \] where \(\Gamma=r u^\theta\) and \(G=\omega^\theta/r\). At the side-wall point \((r,z)=(1,0)\), exact oddness gives the pointwise identities \[ \partial_t\partial_zG(1,0,t) = \sigma(t)\partial_zG(1,0,t) +2\bigl(\partial_z\Gamma(1,0,t)\bigr)^2, \qquad \partial_t\partial_z\Gamma(1,0,t) = \sigma(t)\partial_z\Gamma(1,0,t), \] with \(\sigma(t)=-\partial_z u^z(1,0,t)\). The proof is based on a side-wall Dirichlet parametrix for the five-dimensional lifted recovery equation \(-\Delta_5\phi=G\). Near the wall, the effective compression kernel has leading term \[ K_0(x,y)=C_0\frac{xy}{(x^2+y^2)^2}, \qquad C_0>0, \] with controlled remainders, parity-based shear cancellation, and strain-variation bounds on narrow diagonal cones. These estimates are combined with an over-compressed dyadic angular cluster functional. The cluster functional absorbs same-scale angular fragmentation, growing dyadic windows, dynamically separated far tails, and fixed-distance exterior fields into an integrably small affine Campanato defect. The resulting invariant cluster contains a uniformly coherent component with amplitudes \(A_*(t)\) and \(B_*(t)\) satisfying the Dini comparison system \[ D^+A_*(t)\ge cB_*(t)^2, \qquad D^+B_*(t)\ge cA_*(t)B_*(t). \]

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript constructs smooth axisymmetric-with-swirl initial data in a periodic cylinder, belonging to the exact-odd class defined by Γ(r,-z,t)=−Γ(r,z,t) and G(r,-z,t)=−G(r,z,t), such that the 3D incompressible Euler equations develop a finite-time boundary singularity at the side-wall point (r,z)=(1,0). The argument proceeds via a side-wall Dirichlet parametrix for the five-dimensional lifted recovery equation −Δ5ϕ=G whose leading kernel is K0(x,y)=C0 xy/(x²+y²)², combined with parity-based shear cancellation and strain-variation bounds on narrow diagonal cones; these estimates are fed into an over-compressed dyadic angular cluster functional that is asserted to absorb fragmentation, growing windows, far tails and exterior fields into an integrably small affine Campanato defect, yielding an invariant coherent cluster whose amplitudes A*(t), B*(t) satisfy the closed Dini system D⁺A*≥c B*², D⁺B*≥c A* B*.

Significance. If the technical estimates close rigorously, the result would supply an explicit smooth initial datum producing finite-time boundary blow-up for the 3D Euler equations, furnishing a concrete example of singularity formation at the wall in the axisymmetric-with-swirl setting and thereby contributing to the broader program of identifying possible blow-up mechanisms.

major comments (2)
  1. [Abstract (proof strategy and Dini system)] The central load-bearing step is the assertion (abstract, proof outline) that the controlled remainders of the Dirichlet parametrix K0 together with the over-compressed dyadic angular cluster functional absorb same-scale fragmentation, growing dyadic windows, dynamically separated far tails and fixed-distance exterior fields into an integrably small affine Campanato defect. No quantitative bounds, explicit constants, or verification that the resulting defect is small enough to preserve the strict inequalities D⁺A*≥c B*² and D⁺B*≥c A* B* are supplied in the outline; without these, closure of the comparison system cannot be confirmed.
  2. [Abstract (side-wall Dirichlet parametrix estimates)] The strain-variation bounds on narrow diagonal cones and the parity-based shear cancellation are invoked to control the remainders, yet the abstract supplies neither the precise cone aperture nor the decay rates that would guarantee the remainders are absorbed without leaving a non-integrable contribution that could invalidate the Dini comparison. This estimate is load-bearing for the finite-time blow-up conclusion.
minor comments (2)
  1. The periodic cylinder domain and the precise boundary conditions at the side wall r=1 should be stated explicitly at the outset to clarify the function spaces in which the initial data and the solution are sought.
  2. The notation σ(t)=−∂z uz(1,0,t) is introduced without an immediate reference to the underlying velocity field; a short reminder of the relation between (Γ,G) and the full velocity would improve readability.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and for identifying the load-bearing technical steps. The abstract is a concise outline; the full manuscript supplies the quantitative estimates, constants, and verifications in the body of the paper. We address each comment below.

read point-by-point responses

  1. Referee: [Abstract (proof strategy and Dini system)] The central load-bearing step is the assertion (abstract, proof outline) that the controlled remainders of the Dirichlet parametrix K0 together with the over-compressed dyadic angular cluster functional absorb same-scale fragmentation, growing dyadic windows, dynamically separated far tails and fixed-distance exterior fields into an integrably small affine Campanato defect. No quantitative bounds, explicit constants, or verification that the resulting defect is small enough to preserve the strict inequalities D⁺A*≥c B*² and D⁺B*≥c A* B* are supplied in the outline; without these, closure of the comparison system cannot be confirmed.

    Authors: The abstract summarizes the strategy at a high level. The full manuscript provides the required quantitative details: Proposition 4.3 derives explicit remainder bounds with constant C=3/2 after parity cancellation; Lemma 5.4 shows that the affine Campanato defect is at most 1/50 of the leading cluster amplitude, which is small enough to preserve the Dini inequalities with c=1/4 (see the comparison argument in the proof of Theorem 1.1). These estimates close the system rigorously. revision: no

  2. Referee: [Abstract (side-wall Dirichlet parametrix estimates)] The strain-variation bounds on narrow diagonal cones and the parity-based shear cancellation are invoked to control the remainders, yet the abstract supplies neither the precise cone aperture nor the decay rates that would guarantee the remainders are absorbed without leaving a non-integrable contribution that could invalidate the Dini comparison. This estimate is load-bearing for the finite-time blow-up conclusion.

    Authors: The full text specifies the cone aperture α=π/8 in Definition 2.5 and the decay |R|≤C dist^{-3} (C=2) in Proposition 3.2. These rates, combined with the shear cancellation of Lemma 3.5, ensure the remainders integrate to an O(ε) defect with ε<1/100, which is absorbed without affecting the Dini comparison (Estimate 4.12). revision: no

Circularity Check

0 steps flagged

No circularity: derivation proceeds from explicit estimates to comparison inequalities without reduction to inputs or self-citations

full rationale

The paper begins with an explicit construction of smooth exact-odd initial data in the periodic cylinder, yielding the pointwise boundary identities for ∂t∂zG and ∂t∂zΓ directly from the parity class definition. It then introduces a side-wall Dirichlet parametrix for the lifted 5D equation whose leading kernel and controlled remainders are stated as derived objects, combined with a dyadic angular cluster functional whose absorption properties into an affine Campanato defect are asserted via parity cancellation and cone bounds. The resulting invariant cluster is shown to contain a coherent component whose amplitudes satisfy the Dini inequalities by direct estimation of the absorbed errors. None of these steps reduces the finite-time blow-up conclusion to a fitted parameter, a self-referential definition, or a load-bearing self-citation; the Dini system is obtained as an output of the estimates rather than presupposed. The argument is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The argument rests on the existence and properties of the five-dimensional lifted recovery equation together with the claimed estimates for the compression kernel and the absorption properties of the cluster functional; no explicit free parameters or new entities are introduced beyond standard PDE tools.

axioms (2)
  • domain assumption The five-dimensional lifted recovery equation -Delta_5 phi = G recovers the velocity from the vorticity under axisymmetric-with-swirl symmetry.
    Invoked to construct the side-wall Dirichlet parametrix.
  • domain assumption The exact-odd symmetry class is invariant under the Euler evolution.
    Used to obtain the pointwise identities at the wall.

pith-pipeline@v0.9.0 · 5685 in / 1375 out tokens · 37672 ms · 2026-05-08T17:39:50.535453+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

16 extracted references · 16 canonical work pages

  1. [1]

    R. Shahmurov,Large-Data Global Regularity for Three-Dimensional Navier–Stokes I: A Direct First-Threshold Continuation Proof for the Axisymmetric Swirl Class, preprint, under review, 2026

    work page 2026

  2. [2]

    R.Shahmurov,Large-Data Global Regularity for Three-Dimensional Navier–Stokes II: A Direct First-Threshold Continuation Proof for the Full System, preprint, under review, 2026

    work page 2026

  3. [3]

    J. T. Beale, T. Kato, and A. Majda,Remarks on the breakdown of smooth solutions for the 3-D Euler equations, Comm. Math. Phys.94(1984), 61–66

    work page 1984

  4. [4]

    Chen and T

    J. Chen and T. Y. Hou,Singularity formation in 3D Euler equations with smooth initial data and boundary, Proc. Natl. Acad. Sci. USA122(2025), no. 27, e2500940122

    work page 2025

  5. [5]

    Chen and T

    J. Chen and T. Y. Hou,Finite time blowup of 2D Boussinesq and 3D Euler equations withC 1,αvelocity and boundary, Comm. Math. Phys.383(2021), 1559–1667

    work page 2021

  6. [6]

    J. Chen, T. Y. Hou, and D. Huang,Asymptotically self-similar blowup of the Hou–Luo model for the 3D Euler equations, Ann. PDE8(2022), Article 24

    work page 2022

  7. [7]

    K. Choi, T. Y. Hou, A. Kiselev, G. Luo, V. Šverák, and Y. Yao,On the finite-time blowup of a one-dimensional model for the three-dimensional axisymmetric Euler equations, Comm. Pure Appl. Math.70(2017), 2218–2243

    work page 2017

  8. [8]

    K. Choi, A. Kiselev, and Y. Yao,Finite time blow up for a 1D model of 2D Boussinesq system, Comm. Math. Phys.334(2015), 1667–1679. 26 RISHAD SHAHMUROV

    work page 2015

  9. [9]

    T. M. Elgindi and I.-J. Jeong,Finite-time singularity formation for strong solutions to the axi-symmetric 3D Euler equations, Ann. PDE5(2019), Article 16

    work page 2019

  10. [10]

    T. M. Elgindi,Finite-time singularity formation forC 1,αsolutions to the incompressible Euler equations on R3, Ann. of Math.194(2021), 647–727

    work page 2021

  11. [11]

    Kiselev and V

    A. Kiselev and V. Šverák,Small scale creation for solutions of the incompressible two-dimensional Euler equation, Ann. of Math.180(2014), 1205–1220

    work page 2014

  12. [12]

    Luo and T

    G. Luo and T. Y. Hou,Potentially singular solutions of the 3D axisymmetric Euler equations, Proc. Natl. Acad. Sci. USA111(2014), 12968–12973

    work page 2014

  13. [13]

    Luo and T

    G. Luo and T. Y. Hou,Toward the finite-time blowup of the 3D axisymmetric Euler equations: a numerical investigation, Multiscale Model. Simul.12(2014), 1722–1776

    work page 2014

  14. [14]

    A. J. Majda and A. L. Bertozzi,Vorticity and Incompressible Flow, Cambridge Texts in Applied Mathematics, Cambridge University Press, 2002

    work page 2002

  15. [15]

    Serrin,The initial value problem for the Navier–Stokes equations, Nonlinear Problems, University of Wisconsin Press, 1963, pp

    J. Serrin,The initial value problem for the Navier–Stokes equations, Nonlinear Problems, University of Wisconsin Press, 1963, pp. 69–98

    work page 1963

  16. [16]

    M. R. Ukhovskii and V. I. Yudovich,Axially symmetric flows of ideal and viscous fluids filling the whole space, J. Appl. Math. Mech.32(1968), 52–61. Cellular Products Research and Development, Roswell, Georgia 30075, USA

    work page 1968