pith. sign in

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

Euler Singularities II: Interior Quadrupole Blow-Up for Smooth Axisymmetric Euler with Swirl in texorpdfstring{mathbb R³}

Rishad Shahmurov This is my paper

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

classification 🧮 math.AP
keywords axisymmetric Euler equationsfinite-time singularityquadrupole blow-upswirlinterior mechanismBiot-Savart strainbootstrap estimatescomparison system
0
0 comments X

The pith

Smooth axisymmetric Euler flows with swirl develop finite-time singularities through an interior quadrupole mechanism.

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

The paper constructs smooth, divergence-free initial data for the axisymmetric incompressible Euler equations with swirl in three-dimensional space that lead to finite-time blow-up of the velocity gradient. It localizes the construction away from the axis and introduces quadrupole profiles for the active vorticity and swirl that interact through the Biot-Savart law to generate positive hyperbolic strain. A full four-quadrant quadrupole score is tracked alongside a coercive subscore, and a collection of bootstrap estimates closes to produce a comparison system of differential inequalities. The system forces the score to blow up in finite time, which in turn implies blow-up of the L^∞ norm of the velocity gradient.

Core claim

For suitable smooth initial data the axisymmetric Euler equations with swirl admit solutions in which the vorticity profile G and swirl profile Γ maintain approximate quadrupole forms G ≈ a(t)xy and Γ ≈ Γ_*(t) + (1/2)b(t)xy². These forms produce positive interior strain via the Biot-Savart integral while the Euler source term regenerates the same quadrupole shape. The resulting master estimates reduce to the comparison system Q' ≥ cC, C' ≥ cQC, C ≥ κQ², so that the tracked quadrupole score Q(t) becomes infinite at a finite comparison time and the strain lower bound forces ||∇u(t)||_L^∞ to blow up.

What carries the argument

The interior quadrupole mechanism consisting of the active profiles G(x,y,t) ≈ a(t)xy and Γ(x,y,t) ≈ Γ_*(t) + (1/2)b(t)xy² together with the full four-quadrant quadrupole score and its coercive subscore, which together close the bootstrap estimates and drive the comparison system.

If this is right

  • The L^∞ norm of the velocity gradient blows up in finite time.
  • The singularity forms in the interior away from the axis of symmetry.
  • Smooth decaying initial data exist that enter and remain in the required bootstrap regime.
  • The full set of named estimates on kernel signs, source compatibility, and profile defects persist up to the blow-up time.

Where Pith is reading between the lines

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

  • The same localized quadrupole construction could be tested numerically by evolving the explicit initial data to observe whether the predicted score growth occurs.
  • The comparison system offers a template that might be adapted to other symmetric or reduced fluid models where strain amplification is present.
  • Because the mechanism is interior and self-contained, it raises the possibility of constructing singularities in settings without boundaries or with different symmetry reductions.

Load-bearing premise

The initial data must enter and remain inside the quadrupole bootstrap regime with every listed geometric and analytic constraint holding throughout the evolution.

What would settle it

A numerical solution of the constructed initial data in which the quadrupole score Q(t) stays bounded while the strain remains bounded below by a positive multiple of Q(t), or in which the source term ceases to regenerate the quadrupole shape.

read the original abstract

We present a self-contained interior quadrupole mechanism for finite-time singularity formation in the axisymmetric three-dimensional incompressible Euler equations with swirl in the whole space. The construction is localized away from the axis. In local variables \[ x=r-r_*(t),\qquad y=z, \] centered at a tracked radial point, the active vorticity and swirl profiles are \[ G(x,y,t)\approx a(t)xy, \qquad \Gamma(x,y,t)\approx \Gamma_*(t)+\frac12 b(t)xy^2, \qquad \Gamma_*(t)>0. \] The first profile produces a positive interior Biot--Savart hyperbolic strain; the second profile makes the Euler source term in the equation for \(G=\omega^\theta/r\) regenerate the same quadrupole shape. The active quantity is the full four-quadrant quadrupole score, while a narrow diagonal sector is used only as a coercive subscore. We give the notation and the 5D recovery formula connecting the 3D axisymmetric variables to the lifted elliptic problem, construct explicit smooth decaying divergence-free data, verify their initial entry into the quadrupole bootstrap, prove the master propagation estimates, and derive the comparison system \[ Q'(t)\ge cC(t),\qquad C'(t)\ge cQ(t)C(t),\qquad C(t)\ge \kappa Q(t)^2. \] Consequently the tracked quadrupole score blows up in finite comparison time, and the strain lower bound gives blow-up of \(\norm{\nabla u(t)}_{L^\infty}\). All geometric and analytic constraints used by the construction are stated as named estimates: the interior quadrupole kernel sign expansion, source compatibility, swirl-jet amplification, full-score/coercive-subscore comparison, angular-profile defect persistence, radial-center tracking, neutral-jet hierarchy, and two-sided Dini bounds. This is Part II of a two-paper Euler series; Part I treats boundary blow-up in a periodic cylinder.

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

0 major / 2 minor

Summary. The manuscript constructs explicit smooth, decaying, divergence-free initial data for the axisymmetric 3D incompressible Euler equations with swirl in R^3 that enters an interior quadrupole bootstrap regime localized away from the axis. In local variables (x = r - r_*(t), y = z), it prescribes quadrupole profiles G ≈ a(t)xy and Γ ≈ Γ_*(t) + (1/2)b(t)xy², establishes named estimates (interior kernel sign expansion, source compatibility, swirl-jet amplification, full-score/coercive-subscore comparison, angular defect persistence, radial-center tracking, neutral-jet hierarchy, two-sided Dini bounds), proves master propagation inequalities, and derives the comparison system Q' ≥ cC, C' ≥ cQC, C ≥ κQ². Solutions of this system blow up in finite time, yielding ||∇u||_∞ blow-up via the strain lower bound. The argument is self-contained on its own terms and uses a 5D recovery formula to connect to the lifted elliptic problem.

Significance. If the derivations hold, the result supplies an explicit, verifiable example of finite-time interior singularity formation for smooth axisymmetric Euler flows with swirl, complementing the boundary mechanism of Part I. The explicit initial-data construction, parameter-free comparison system derived from Biot-Savart expansions, and complete list of bootstrap constraints as named estimates are genuine strengths that render the argument falsifiable and potentially machine-checkable in principle. This advances the program of constructing singularities in 3D Euler without relying on ad-hoc fitting or hidden circularity.

minor comments (2)
  1. The 5D recovery formula is invoked repeatedly but its precise statement and derivation from the 3D axisymmetric variables appear only in outline; expanding the formula with explicit coordinate expressions would aid verification of the lifted elliptic problem.
  2. Notation for the full quadrupole score versus the coercive subscore is introduced in the abstract and used throughout, yet a single consolidated definition table or equation block would reduce cross-referencing.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive and detailed summary of our manuscript, which accurately reflects the construction of the interior quadrupole blow-up mechanism, the named estimates, the comparison system, and the connection to the 5D recovery formula. The recommendation for minor revision is noted, and we will prepare a revised version accordingly. No specific major comments were listed in the report.

Circularity Check

0 steps flagged

No significant circularity; derivation self-contained from explicit data and derived estimates

full rationale

The paper constructs explicit smooth decaying divergence-free initial data, verifies its entry into the quadrupole bootstrap regime, proves master propagation estimates from the Euler equations and Biot-Savart kernel expansions, and derives the comparison system Q' ≥ cC, C' ≥ cQC, C ≥ κQ² directly from those estimates. The finite-time blow-up of the tracked quadrupole score then follows by standard ODE comparison, yielding the strain lower bound on ||∇u||_∞. All geometric and analytic constraints are stated as named estimates proved within the paper itself. No step reduces by construction to a fitted input, self-definition, or load-bearing self-citation; the argument is independent of external fitted parameters or unverified prior results.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The construction relies on standard properties of the Biot-Savart law in axisymmetric coordinates and on the existence of smooth decaying divergence-free vector fields; no new physical entities are postulated, but several comparison constants and profile coefficients are introduced via the named estimates.

free parameters (1)
  • comparison constants c and κ
    Positive constants appearing in the differential inequalities for the quadrupole score Q and coercive subscore C; derived from kernel sign expansions rather than fitted to data.
axioms (2)
  • domain assumption The Biot-Savart law produces a positive interior hyperbolic strain from the quadrupole vorticity profile G ≈ a(t)xy.
    Invoked in the description of the active quantity and the strain lower bound.
  • domain assumption The swirl profile Γ ≈ Γ_* + (1/2)b(t)xy² regenerates the quadrupole shape via the Euler source term.
    Central to the self-reinforcing mechanism stated in the abstract.

pith-pipeline@v0.9.0 · 5679 in / 1644 out tokens · 50980 ms · 2026-05-08T17:31:20.896353+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

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

    Constantin, C

    P. Constantin, C. Fefferman, and A. Majda,Geometric constraints on potentially singular solutions for the 3-D Euler equations, Comm. Partial Differential Equations 21 (1996), 559–571

    work page 1996

  5. [5]

    A. J. Majda and A. L. Bertozzi,Vorticity and Incompressible Flow, Cambridge University Press, 2002

    work page 2002

  6. [6]

    Luo and T

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

    work page 2014

  7. [7]

    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

  8. [8]

    T. Y. Hou and G. Luo,Potentially singular solutions of the 3D axisymmetric Euler equations, SIAM Review 61 (2019), 661–708

    work page 2019

  9. [9]

    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. USA 122 (2025), e2500940122

    work page 2025

  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]

    T. M. Elgindi, T.-E. Ghoul, and N. Masmoudi,Stable self-similar blow-up for a family of nonlocal transport equations, Anal. PDE 14 (2021), 891–908

    work page 2021

  12. [12]

    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

  13. [13]

    S. C. Preston and A. Sarria,Lagrangian aspects of the axisymmetric Euler equation, J. Math. Fluid Mech. 17 (2015), 85–101. Cellular Products Research and Development, Roswell, Georgia 30075, USA

    work page 2015