Stability Framework for the Singularity of the Euler Equations on $\mathbb{R}^3$

Technical Summary: Stability Framework for the Singularity of the Euler Equations on $\mathbb{R}^3$

Problem Statement
The paper addresses the central open problem in fluid dynamics regarding whether smooth initial data for the 3D incompressible Euler equations can lead to finite-time singularity formation (blowup). While finite-time blowup has been established for the Euler equations with boundaries or under non-smooth initial conditions, the existence of a singularity arising from smooth initial data on the unbounded domain $\mathbb{R}^3$ remains unproven. A common strategy to resolve this is to construct an approximate singular profile and prove its nonlinear stability. The authors note that while a high-precision approximate self-similar profile was discovered in a companion numerical study using Physics-Informed Neural Networks (PINNs), a rigorous stability proof for this profile on $\mathbb{R}^3$ had not been established. The specific challenge lies in the lack of a global "outgoing" flow property (which typically aids stability) and the presence of nontrivial fixed points in the meridional flow that can cause perturbations to concentrate.

Methodology
The paper establishes a rigorous framework for proving the nonlinear stability of the approximate self-similar profile discovered in the companion study. The methodology proceeds through the following stages:

1. Approximate Profile Representation: The numerically discovered PINN profile is converted into a piecewise polynomial spline representation. This analytic form allows for exact differentiation and rigorous evaluation of residuals and norms using interval arithmetic (via the Arb library), ensuring that numerical errors are bounded (e.g., $O(10^{-5})$ in $L^2$).
2. Dynamic Rescaling and Linearization: The authors employ a dynamic rescaling formulation where the solution is viewed in a moving frame that expands or contracts to keep the singularity at a fixed scale. The approximate profile becomes a steady state in this rescaled frame. The stability analysis involves linearizing the rescaled Euler equations around this steady state.
3. Scale-Consistent Variables: To handle the differing scaling behaviors of velocity and vorticity, the analysis is formulated using scale-consistent perturbation variables: the vorticity perturbation $\delta\omega$ and the gradient of the velocity perturbation $\nabla \delta u$.
4. Modulation and Normalization: The framework introduces modulation parameters to fix the translation and amplitude of the profile, effectively removing neutral directions (symmetry modes) from the stability analysis. This ensures that the stability question concerns genuine perturbations transverse to these symmetries.
5. Weighted Energy Estimates: The core of the proof relies on constructing a full energy functional $E_k^2 = E_0^2 + \mu_k H_k^2$, combining:
   • Low-order weighted energy ($E_0$): Uses singular weight functions $\Phi_\omega, \Phi_r, \Phi_z$ adapted to the profile to establish linear damping.
   • High-order weighted energy ($H_k$): Uses higher-order weights to control derivatives and pointwise values necessary for closing nonlinear estimates.
6. Computer-Assisted Certification: The proof reduces the infinite-dimensional stability problem to a finite-dimensional optimization problem. The authors derive explicit analytic estimates for linear damping, nonlinear interactions, and PDE residuals. These estimates depend on a large collection of explicit constants (e.g., matrix bounds, interpolation constants, elliptic operator norms). The framework requires these constants to be rigorously certified using interval arithmetic and certified matrix bounds.
7. Formalization: The paper notes a parallel effort (LeanPDE) to formalize the symbolic derivations and proof steps in the Lean theorem prover, connecting them with the certified numerical computations.

Key Contributions

• Stability Framework: The primary contribution is the construction of a detailed, modular framework for proving the nonlinear stability of a candidate self-similar blowup profile for the 3D Euler equations on $\mathbb{R}^3$.
• Reduction to Finite Verification: The authors demonstrate that the stability proof can be reduced to the rigorous certification of a finite set of explicit constants and estimates. This shifts the burden of proof from qualitative analysis to quantitative verification.
• Handling Non-Outgoing Flows: The framework successfully adapts stability techniques to a setting without a global outgoing property, utilizing a weaker local condition to push the flow away from fixed points and employing delicate high-order damping estimates.
• Spline-Based Certification: The use of piecewise polynomial splines to represent the numerical profile allows for the exact evaluation of PDE residuals and derivatives, a necessary step for rigorous interval arithmetic certification.
• Two-Radii Stability Theorem: The paper provides a generalized stability theorem (Theorem 2) that allows for separate bounds on low-order and high-order energies, offering flexibility in the certification process.

Results
The paper does not claim to have completed the final numerical certification of all constants required to close the proof. Instead, it establishes the architecture for such a proof.

• Theoretical Completeness: The authors prove that if the explicit constants (damping margins, nonlinear bounds, residual norms) can be certified to satisfy specific inequalities (e.g., $\Lambda_{stab} > 0$), then the rescaled profile is nonlinearly stable.
• Conditional Stability: Conditional on the rigorous certification of these constants, the framework guarantees that the rescaled profile is stable. Furthermore, via the dynamic rescaling mechanism, this stability implies the existence of an admissible solution in the original physical variables that develops a singularity in finite time.
• Residual Bounds: The paper reports that the spline representation of the profile satisfies the steady-state profile equations with residuals bounded by $O(10^{-5})$ in $L^2$ and $O(10^{-3})$ in $L^\infty$, comparable to the original PINN results.

Significance
The paper claims that the principal remaining obstacles to proving finite-time blowup for the 3D Euler equations with smooth initial data are now computational and quantitative rather than conceptual. By providing a rigorous framework that reduces the problem to a finite collection of verifiable estimates, the authors argue that the path to a complete proof is clear. The work represents a critical step in the "computer-assisted proof" paradigm for PDEs, bridging the gap between numerical discovery (PINNs) and rigorous mathematical proof. The significance lies in demonstrating that the stability of a complex, numerically discovered singular profile can be subjected to a systematic, verifiable analysis, potentially resolving one of the Millennium Prize problems if the remaining constants are successfully certified. The paper emphasizes that the modular structure of the argument allows for targeted refinements (e.g., sharpening specific estimates or adjusting weights) without altering the underlying stability mechanism.