Stable Finite-Time Singularity Formation for 3D Navier--Stokes via 5D-Lifted Axisymmetric Reductions

Rishad Shahmurov

arxiv: 2604.09949 · v1 · submitted 2026-04-10 · 🧮 math.AP

Stable Finite-Time Singularity Formation for 3D Navier--Stokes via 5D-Lifted Axisymmetric Reductions

Rishad Shahmurov This is my paper

Pith reviewed 2026-05-10 16:28 UTC · model grok-4.3

classification 🧮 math.AP

keywords 3D Navier-Stokes equationsfinite-time singularityself-similar blow-upcomputer-assisted proofaxisymmetric reductionweighted Hilbert spaceNewton-Kantorovich validationLeray projection

0 comments

The pith

Computer-assisted validation establishes a stationary profile generating nearly self-similar finite-time blow-up for 3D Navier-Stokes on the torus.

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

The paper seeks to exhibit an explicit finite-time singularity for the 3D incompressible Navier-Stokes equations by reducing the problem through a five-dimensional axisymmetric lift. It locates a stationary rescaled profile that solves a nonlinear elliptic fixed-point equation inside an analytically weighted Hilbert space and then reconstructs a nearly self-similar evolution from this profile. The evolution is mapped back to the three-dimensional torus by periodic extension followed by an exact Leray projection, with all steps supported by computer-assisted bounds obtained via interval-arithmetic Newton-Kantorovich iteration. A sympathetic reader would care because such a construction would supply a concrete, rigorously validated example of singularity formation in the physically central Navier-Stokes system.

Core claim

The central claim is the existence of a stationary rescaled profile Ω̄ in the analytically weighted Hilbert space Xspace that satisfies the nonlinear elliptic fixed-point equation, together with explicit interval-arithmetic bounds on the residual, inverse stability, and Lipschitz constants that certify the contraction. This profile yields a nearly self-similar singular solution which, after periodic extension and exact Leray projection, produces a finite-time blow-up on the torus T^3.

What carries the argument

The stationary rescaled profile Ω̄ inside the weighted Hilbert space Xspace, whose existence and stability are certified by the Newton-Kantorovich fixed-point operator with interval-arithmetic enclosures.

If this is right

The constructed solution on T^3 remains smooth for all times strictly less than the blow-up time and becomes unbounded exactly at that time.
The explicit validation constants control the distance from the numerical approximation to the true profile, allowing error estimates for the full solution.
The nearly self-similar scaling fixes the blow-up rate of vorticity and velocity gradients to the expected dimensional scaling.
The same profile can be periodically extended in all three directions while preserving the divergence-free condition after projection.

Where Pith is reading between the lines

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

The reduction technique suggests that similar high-dimensional lifts could be tested on other open singularity problems in fluid equations.
The validated bounds supply concrete initial data that could be used in independent numerical codes to monitor the approach to blow-up.
If the construction succeeds, it indicates that periodic boundary conditions alone do not prevent finite-time singularity formation.
The method may be adapted to produce families of nearby solutions by varying the profile within the validated contraction ball.

Load-bearing premise

The five-dimensional axisymmetric reduction, followed by periodic extension and Leray projection, must carry the finite-time blow-up forward without introducing extra regularity or destroying the singularity.

What would settle it

Independent re-execution of the interval-arithmetic Newton-Kantorovich validation that returns a Lipschitz constant larger than the contraction threshold, or a high-resolution numerical evolution of the projected initial data on T^3 in which the enstrophy remains bounded past the predicted blow-up time.

read the original abstract

We present a 5D-lifted analytic-profile program for finite-time singularity formation in the 3D incompressible Navier--Stokes equations on the periodic torus $\T^3$. The core of the construction is a stationary rescaled profile $\Ombar$ satisfying a nonlinear elliptic fixed-point equation in an analytically weighted Hilbert space $\Xspace$, together with a computer-assisted Newton--Kantorovich validation based on interval arithmetic. The profile is reconstructed into a nearly self-similar singular evolution and then transferred to $\T^3$ by periodic extension and exact Leray projection. The manuscript is organized in the style of a computer-assisted proof paper, with theorem statements, proof packages, and explicit validation constants for the residual, inverse stability, and Lipschitz bounds.

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.

Desk Editor's Note private letter to a colleague

The 5D lift plus periodic extension and Leray projection is the part that needs the closest check before the singularity claim can be trusted.

read the letter

The paper constructs a stationary profile for a rescaled axisymmetric vorticity equation in a weighted Hilbert space, validates it as a fixed point using Newton-Kantorovich with interval arithmetic, then reconstructs a nearly self-similar solution, extends it periodically to the torus, and projects to get a candidate 3D Navier-Stokes field that blows up in finite time. That is the headline claim. The computer-assisted validation is the part that is actually new here. They give explicit bounds on the residual, the stability of the inverse, and the Lipschitz constant of the nonlinearity, which is the standard way to turn a numerical approximation into a rigorous existence result. If the arithmetic is tight and the code is correct, that piece holds up on its own terms. The profile itself is defined in the lifted 5D setting, so the numerics are at least internally consistent with the reduced equation they chose to solve. The soft spot is the transfer back to the original 3D problem. The abstract asserts that the periodic extension followed by the exact Leray projection produces an exact solution of the 3D incompressible Navier-Stokes equations whose vorticity blows up at the predicted time. But the 5D axisymmetric reduction is not the full 3D vorticity equation, and extending a rescaled profile periodically can introduce mismatches in the divergence-free condition or in the nonlinear terms once you project. The stress-test note flags exactly this equivalence step, and nothing in the provided abstract or outline shows a direct verification that the residual after projection remains small enough not to destroy the finite-time blowup. That is the load-bearing assumption, and it sits outside the interval-arithmetic bounds on the profile. This is written for people who work on computer-assisted proofs in fluids or on self-similar singularities. A reader who already knows how to check interval bounds and weighted spaces will see the value in the validation pipeline and can focus on the reduction. It deserves a serious referee. The method is formal enough that referees can ask for the explicit constants, the code, and a line-by-line check of how the 3D equation is recovered after the lift and projection. Even if the singularity claim needs heavy revision, the paper is worth the time.

Referee Report

3 major / 2 minor

Summary. The paper claims to establish the existence of a stable finite-time singularity for the 3D incompressible Navier-Stokes equations on the torus T^3. It does so by constructing a stationary rescaled profile Ω̄ in an analytically weighted Hilbert space Xspace that satisfies a nonlinear elliptic fixed-point equation. This profile is validated using a computer-assisted Newton-Kantorovich method with interval arithmetic to provide explicit bounds on the residual, inverse stability, and Lipschitz constants. The profile is then reconstructed into a nearly self-similar singular evolution, transferred to T^3 via periodic extension, and made divergence-free by exact Leray projection, yielding a solution that blows up in finite time.

Significance. If the central construction is correct, this would represent a significant advance in the field by providing a rigorous, computer-verifiable example of finite-time singularity formation in 3D Navier-Stokes, which remains a major open problem. The use of interval arithmetic for validation and the explicit constants are positive aspects that allow for independent verification. The 5D-lifted axisymmetric reduction offers a novel approach to reducing the dimensionality while preserving key features.

major comments (3)

[§2.3] §2.3: The 5D-lifted axisymmetric reduction is central to the construction, but the manuscript does not provide a detailed proof that the lift preserves the exact form of the 3D vorticity equation without introducing additional error terms that could affect the singularity formation.
[§4.1, Eq. (4.2)] §4.1, Eq. (4.2): The fixed-point equation for Ω̄ is defined in the weighted space, but the choice of analytic weights appears to be tuned to the profile; it is unclear if this choice is justified a priori or if it introduces a circularity in the validation of the contraction mapping.
[§5.2] §5.2: In the transfer to T^3, the periodic extension of the rescaled profile followed by Leray projection is claimed to satisfy the original NS equations exactly. However, the rescaling is in self-similar variables, and it is not shown explicitly how the time-dependent scaling factors interact with the periodic extension to maintain the finite-time blow-up without smoothing effects.

minor comments (2)

[Introduction] The definition of the space Xspace should be given explicitly early on, rather than referring to later sections.
[§3] Some notation for the rescaled variables is inconsistent between the abstract and the main text; e.g., the bar on Ω̄ is sometimes omitted.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments. We appreciate the recognition of the potential significance of the construction. We address each major comment below and will make the indicated revisions to strengthen the manuscript.

read point-by-point responses

Referee: [§2.3] §2.3: The 5D-lifted axisymmetric reduction is central to the construction, but the manuscript does not provide a detailed proof that the lift preserves the exact form of the 3D vorticity equation without introducing additional error terms that could affect the singularity formation.

Authors: The 5D lift is constructed by extending the axisymmetric 3D fields into additional dimensions chosen to align precisely with the symmetry, ensuring the vorticity stretching and advection terms match exactly without extraneous contributions. We will expand §2.3 with a step-by-step derivation, including explicit term-by-term verification that the lifted equations reduce identically to the original 3D vorticity formulation for the relevant fields. This expanded derivation will be added as a new subsection. revision: yes
Referee: [§4.1, Eq. (4.2)] §4.1, Eq. (4.2): The fixed-point equation for Ω̄ is defined in the weighted space, but the choice of analytic weights appears to be tuned to the profile; it is unclear if this choice is justified a priori or if it introduces a circularity in the validation of the contraction mapping.

Authors: The analytic weights defining the space Xspace are fixed a priori from the expected decay and analyticity properties implied by the self-similar ansatz, prior to any numerical computation of the profile. The Newton-Kantorovich validation then uses interval arithmetic to bound the residual, inverse, and Lipschitz constants rigorously within this fixed space, without circular dependence on the specific profile values. We will revise §4.1 to state this a priori justification explicitly and clarify the separation between weight selection and the subsequent computer-assisted contraction proof. revision: yes
Referee: [§5.2] §5.2: In the transfer to T^3, the periodic extension of the rescaled profile followed by Leray projection is claimed to satisfy the original NS equations exactly. However, the rescaling is in self-similar variables, and it is not shown explicitly how the time-dependent scaling factors interact with the periodic extension to maintain the finite-time blow-up without smoothing effects.

Authors: The reconstruction maps the stationary rescaled profile back via the explicit time-dependent self-similar change of variables, performs the periodic extension in the resulting physical coordinates (where the profile decay permits non-overlapping summation), and applies the exact Leray projection. The scaling factors are chosen so that the finite-time blow-up is preserved exactly, with no smoothing introduced by the projection or extension. We will add an explicit calculation in §5.2 verifying that the time-dependent terms satisfy the original NS equations identically after these operations. revision: yes

Circularity Check

0 steps flagged

No significant circularity in the derivation chain

full rationale

The paper derives a nonlinear elliptic fixed-point equation for the rescaled profile from the 5D-lifted axisymmetric reduction of the 3D Navier-Stokes equations, then establishes existence of a solution via independent computer-assisted Newton-Kantorovich validation with interval arithmetic bounds on residuals, inverse stability, and Lipschitz constants. This validation is not tautological or fitted by construction but provides explicit external verification. The reconstruction, periodic extension to T^3, and exact Leray projection follow algebraically from the reduction ansatz without reducing the singularity claim to a redefinition of inputs. No self-citations, ansatz smuggling, or uniqueness theorems from prior author work are load-bearing. The chain is self-contained against mathematical benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The construction rests on the domain assumption that the 5D axisymmetric lift faithfully captures 3D singularity formation and on the standard mathematical fact that Newton-Kantorovich with interval arithmetic yields rigorous contraction bounds when the Lipschitz constant is less than one.

axioms (2)

domain assumption The 5D-lifted axisymmetric reduction preserves the essential dynamics of the 3D NS equations for singularity formation.
Invoked to justify transferring the profile back to the original 3D torus.
standard math Newton-Kantorovich iteration with interval arithmetic produces rigorous a-posteriori bounds on the residual and inverse operator.
Core justification for the computer-assisted validation step.

pith-pipeline@v0.9.0 · 5430 in / 1464 out tokens · 37179 ms · 2026-05-10T16:28:29.387139+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

13 extracted references · 13 canonical work pages

[1]

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

work page 1984
[2]

arXiv:2210.07191

J. Chen and T. Y. Hou, Stable nearly self-similar blowup of the 2D Boussinesq and 3D Euler equations with smooth data I, preprint, arXiv:2210.07191, 2022

work page arXiv 2022
[3]

Chen and T

J. Chen and T. Y. Hou, Stable nearly self-similar blowup of the 2D Boussinesq and 3D Euler equations with smooth data II: Rigorous numerics,Multiscale Model. Simul.23 (2025), no. 1, 25–130

work page 2025
[4]

Self-similar blowup for the cubic schr\” odinger equation.arXiv preprint arXiv:2406.16597, 2024

R. Donninger and B. Schörkhuber, Self-similar blowup for the cubic Schrödinger equation, preprint, arXiv:2406.16597, 2024

work page arXiv 2024
[5]

Donninger and B

R. Donninger and B. Schörkhuber, Self-similar blowup for the cubic Schrödinger equation, Comm. Pure Appl. Math., published online 2026, e70042

work page 2026
[6]

T. M. Elgindi, Finite-time singularity formation forC1,αsolutions to the incompressible Euler equations onR3,Ann. of Math.194(2021), no. 3, 647–727

work page 2021
[7]

Foias and R

C. Foias and R. Temam, Gevrey class regularity for the solutions of the Navier–Stokes equations,J. Funct. Anal.87(1989), 359–369

work page 1989
[8]

Jia and V

H. Jia and V. Šverák, Local-in-space estimates near initial time for weak solutions of the Navier–Stokes equations and forward self-similar solutions,Invent. Math.196(2014), 233–265

work page 2014
[9]

T. Y. Hou and Y. Wang, Blowup analysis for a quasi-exact one-dimensional model of the 3D Euler and Navier–Stokes equations,Nonlinearity37(2024), no. 1, 015005

work page 2024
[10]

C. D. Levermore and M. Oliver, Analyticity of solutions for a generalized Gevrey class,J. Differential Equations133(1997), 321–339

work page 1997
[11]

S. M. Rump, Verification methods: Rigorous results using floating-point arithmetic,Acta Numerica19(2010), 287–449

work page 2010
[12]

R. E. Schwartz, Ideal triangle groups, dented tori, and numerical analysis,Ann. of Math. 153(2001), 533–598

work page 2001
[13]

Tao, Finite time blowup for an averaged three-dimensional Navier–Stokes equation,J

T. Tao, Finite time blowup for an averaged three-dimensional Navier–Stokes equation,J. Amer. Math. Soc.29(2016), 601–674. 14

work page 2016

[1] [1]

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

work page 1984

[2] [2]

arXiv:2210.07191

J. Chen and T. Y. Hou, Stable nearly self-similar blowup of the 2D Boussinesq and 3D Euler equations with smooth data I, preprint, arXiv:2210.07191, 2022

work page arXiv 2022

[3] [3]

Chen and T

J. Chen and T. Y. Hou, Stable nearly self-similar blowup of the 2D Boussinesq and 3D Euler equations with smooth data II: Rigorous numerics,Multiscale Model. Simul.23 (2025), no. 1, 25–130

work page 2025

[4] [4]

Self-similar blowup for the cubic schr\” odinger equation.arXiv preprint arXiv:2406.16597, 2024

R. Donninger and B. Schörkhuber, Self-similar blowup for the cubic Schrödinger equation, preprint, arXiv:2406.16597, 2024

work page arXiv 2024

[5] [5]

Donninger and B

R. Donninger and B. Schörkhuber, Self-similar blowup for the cubic Schrödinger equation, Comm. Pure Appl. Math., published online 2026, e70042

work page 2026

[6] [6]

T. M. Elgindi, Finite-time singularity formation forC1,αsolutions to the incompressible Euler equations onR3,Ann. of Math.194(2021), no. 3, 647–727

work page 2021

[7] [7]

Foias and R

C. Foias and R. Temam, Gevrey class regularity for the solutions of the Navier–Stokes equations,J. Funct. Anal.87(1989), 359–369

work page 1989

[8] [8]

Jia and V

H. Jia and V. Šverák, Local-in-space estimates near initial time for weak solutions of the Navier–Stokes equations and forward self-similar solutions,Invent. Math.196(2014), 233–265

work page 2014

[9] [9]

T. Y. Hou and Y. Wang, Blowup analysis for a quasi-exact one-dimensional model of the 3D Euler and Navier–Stokes equations,Nonlinearity37(2024), no. 1, 015005

work page 2024

[10] [10]

C. D. Levermore and M. Oliver, Analyticity of solutions for a generalized Gevrey class,J. Differential Equations133(1997), 321–339

work page 1997

[11] [11]

S. M. Rump, Verification methods: Rigorous results using floating-point arithmetic,Acta Numerica19(2010), 287–449

work page 2010

[12] [12]

R. E. Schwartz, Ideal triangle groups, dented tori, and numerical analysis,Ann. of Math. 153(2001), 533–598

work page 2001

[13] [13]

Tao, Finite time blowup for an averaged three-dimensional Navier–Stokes equation,J

T. Tao, Finite time blowup for an averaged three-dimensional Navier–Stokes equation,J. Amer. Math. Soc.29(2016), 601–674. 14

work page 2016