Fluid Mechanics on Nam Le

Navier–Stokes Existence and Smoothness

Fri, 29 May 2026 00:00:00 +0000

The motion of a viscous incompressible fluid is described by the Navier–Stokes equations, first written down by Claude-Louis Navier in 1822 and given their modern form by George Gabriel Stokes. Whether smooth solutions to these equations can always be continued for all time (or whether they can spontaneously develop a singularity at some finite time) is one of the deepest open problems in mathematics, and one of the seven Clay Millennium Prize Problems, carrying a 1,000,000$ prize for a solution.

Problem (Clay Millennium Prize, Fefferman 2000)

Let $u_0 : \mathbb{R}^3 \to \mathbb{R}^3$ be a smooth divergence-free vector field. Does there exist a smooth solution $u(x,t)$, $p(x,t)$ to the 3D incompressible Navier–Stokes equations $$\partial_t u + (u \cdot \nabla)u - \nu\Delta u + \nabla p = 0, \qquad \nabla \cdot u = 0, \qquad u(\cdot,0) = u_0$$ defined for all $t > 0$ and satisfying $\int_{\mathbb{R}^3}|u(x,t)|^2,dx < C$ for all $t \geq 0$? A solution or a counterexample (a smooth $u_0$ for which no such smooth solution exists) both qualify for the prize.

The Equations and Their Scaling #

Compared to the Euler equations (which describe inviscid flow), the Navier–Stokes equations add the viscous term $\nu\Delta u$, where $\nu > 0$ is the kinematic viscosity. This term dissipates energy and regularises the flow locally. The central tension is that the nonlinear term $(u\cdot\nabla)u$ can concentrate energy at small spatial scales faster than viscosity can diffuse it away.

Scaling symmetry. The Navier–Stokes equations are invariant under the rescaling $$u(x,t) \mapsto \lambda u(\lambda x,, \lambda^2 t), \qquad p(x,t) \mapsto \lambda^2 p(\lambda x,, \lambda^2 t).$$ A norm is critical (or scale-invariant) if it is preserved by this rescaling. The critical norm in $L^p(\mathbb{R}^3)$ is $L^3$, since $|\lambda u(\lambda\cdot)| _{L^3} = |u| _{L^3}$. The energy norm $|u| _{L^2}$ is subcritical: it scales as $\lambda^{1/2}|u| _{L^2}$, which shrinks under the rescaling $\lambda \to \infty$ (i.e., zoom into small scales). This mismatch is the core of the difficulty: global energy control does not prevent concentration at arbitrarily small scales.

2D global regularity. In two dimensions the scaling is different: the enstrophy $|\nabla u|_{L^2}^2$ is scale-invariant and is controlled by the energy. Global regularity in 2D follows from this enstrophy estimate, a fact known since the 1960s. In 3D no analogous critical quantity is controlled globally, and the problem is open.

The Hierarchy of Known Results #

Leray–Hopf Weak Solutions (1934) #

Theorem (Leray 1934, Hopf 1951)

For any $u_0 \in L^2(\mathbb{R}^3)$ divergence-free, there exists a global weak solution $u \in L^\infty(0,\infty;, L^2) \cap L^2(0,\infty;, H^1)$ satisfying the energy inequality $$|u(t)| _{L^2}^2 + 2\nu\int _0^t |\nabla u| _{L^2}^2, ds \leq |u_0| _{L^2}^2.$$

Leray’s construction, via a compactness argument on regularised equations, produces a solution that is globally defined but potentially not smooth, and the term “weak” refers to the fact that the equations are satisfied only in an integral (distributional) sense, not pointwise. The energy inequality is the only bound available globally. Whether Leray–Hopf solutions are unique, or whether they are the same as smooth solutions when the initial data is smooth, is unknown.

Partial Regularity: The CKN Theorem #

The best known result limiting the size of potential singularities is the following.

Theorem (Caffarelli–Kohn–Nirenberg, 1982)

For any suitable weak solution to the 3D Navier–Stokes equations, the set of space-time singular points has parabolic Hausdorff dimension at most 1. In particular, at any given time the spatial singular set has Hausdorff dimension at most $\dfrac{1}{2}$.

A “suitable weak solution” is a weak solution satisfying a local energy inequality. The CKN theorem proves that singularities, if they exist, cannot fill a curve or surface: they can occupy at most a set of dimension one in space-time. This is the most quantitative partial regularity result available and was simplified by Lin (1998). Scheffer (1977) had earlier shown singular times have Hausdorff dimension at most $\dfrac{1}{2}$.

Conditional Regularity: Ladyzhenskaya–Prodi–Serrin #

Theorem (Ladyzhenskaya 1967, Prodi 1959, Serrin 1962)

If a weak solution additionally satisfies $u \in L^r(0,T;, L^s(\mathbb{R}^3))$ with $\dfrac{2}{r} + \dfrac{3}{s} = 1$ and $3 < s \leq \infty$, then $u$ is smooth on $(0,T]$.

The condition $\dfrac{2}{r} + \dfrac{3}{s} = 1$ is precisely the scale-invariant line in the $(r,s)$ plane: membership in any of these spaces implies regularity. The family ranges from $(r,s)=(\infty, 3)$ (critical $L^3$ control in space, uniform in time) to $(r,s)=(2,\infty)$ (square-integrable $L^\infty$ control in time). These are conditional results: they do not prove that a weak solution lies in such a space, only that if it does, it must be smooth.

The Critical Endpoint: Escauriaza–Seregin–Šverák #

Theorem (Escauriaza–Seregin–Šverák, 2003)

If $u$ is a Leray–Hopf weak solution with $\sup _{t \in [0,T^*)} |u(\cdot,t)| _{L^3(\mathbb{R}^3)} < \infty$, then $u$ can be extended as a smooth solution past $T^*$.

The endpoint case $s=3$ of the LPS family is the critical one: $L^3(\mathbb{R}^3)$ is exactly the scale-invariant norm for Navier–Stokes. The ESS proof is substantially harder than the subcritical cases; it uses a compactness argument to reduce to a smooth, backwards self-similar solution and then invokes a backwards uniqueness theorem for parabolic equations to rule it out.

Tao’s Quantitative Criterion #

Theorem (Tao, 2019)

If a smooth finite-energy solution first becomes singular at time $T^*$, then $$\limsup_{t \uparrow T^*} \dfrac{|u(\cdot,t)| _{L^3(\mathbb{R}^3)}}{\bigl(\log\log\log\tfrac{1}{T^*-t}\bigr)^c} = \infty$$ for some absolute constant $c>0$. In particular, the critical $L^3$ norm must blow up at least as fast as a triple-logarithm in $(T^*-t)^{-1}$.

Tao’s result is the first supercritical regularity criterion for Navier–Stokes: it gives quantitative information about the blowup rate that goes (by a triple logarithm) beyond what scaling alone can detect. The proof quantifies the compactness arguments in the ESS proof, replacing each use of a compactness method by an explicit Carleman inequality, and propagates lower bounds for the vorticity across dyadic annuli. The triple-exponential dependence in Tao’s bound has since been localised and sharpened by Barker–Prange (2021) and others.

The Supercriticality Problem #

The fundamental analytical obstruction is that Navier–Stokes is supercritical with respect to the only globally controlled norm ($L^2$): the energy.

Define the critical regularity index as the Sobolev exponent $s$ such that $\dot{H}^s(\mathbb{R}^3)$ is scale-invariant. For Navier–Stokes, $s = 1/2$. The energy controls $\dot{H}^0 = L^2$ (subcritical), and regularity theory requires control at $\dot{H}^1$ (critical viscous norm) or $L^3$ (critical Lebesgue norm). There is a regularity gap between what is globally available ($L^2$) and what is needed ($L^3$ or $\dot{H}^1$). Every known approach to closing this gap runs into the same obstruction: the nonlinearity can create structure at arbitrarily small scales that the subcritical $L^2$ bound cannot see.

Tao (2016) made this gap precise by constructing an averaged Navier–Stokes system, where the bilinear nonlinearity $(u\cdot\nabla)u$ is replaced by a carefully designed convex average of related nonlinearities, for which finite-time blowup can be rigorously proved. This construction does not produce a counterexample to the true Navier–Stokes equations, but it demonstrates that the specific algebraic structure of the nonlinearity is load-bearing: any proof of global regularity must use something specific about $(u\cdot\nabla)u$ that is not shared by its averages.

Research Directions #

1. Improving the Quantitative Blowup Rate #

Tao’s triple-logarithmic rate is the sharpest known lower bound on blowup of the critical $L^3$ norm. Scaling considerations suggest that the true rate, if blowup occurs, should be much faster; conjecturally $|u|_{L^3} \sim (T^*-t)^{-\delta}$ for some $\delta > 0$, analogous to Type I blowup in nonlinear heat equations. The gap between the triple-logarithmic lower bound and the conjectured power-law rate represents the frontier of quantitative regularity theory. Closing even part of this gap, for instance establishing a single-logarithmic or power-of-log lower bound, would require new ideas beyond Carleman estimates.

2. Type I vs. Type II Blowup #

A blowup is called Type I if the scale-invariant norm $|u(\cdot,t)|_{L^3}$ grows no faster than $O((T^-t)^{-1/2})$ near $T^$. It is Type II otherwise. For the Navier–Stokes equations, ruling out Type I blowup would be a significant advance: all self-similar singularities (where $u(x,t) = (T^*-t)^{-1/2}U(x/(T^*-t)^{1/2})$) are of Type I, and several results (including work of Ružička and Seregin) already rule them out under mild additional assumptions. Whether all Type I blowup can be excluded, leaving only the less structured Type II, is open.

3. Uniqueness of Weak Solutions #

Leray–Hopf weak solutions exist globally, but they may not be unique. This is a separate, equally deep question: even if all smooth solutions extend globally, one must also ask whether weak solutions coincide with smooth ones when started from smooth data. Recent work of Buckmaster and Vicol (2019) showed that weak solutions below the Ladyzhenskaya–Prodi–Serrin threshold are indeed non-unique, using convex integration techniques developed for the Euler equations (De Lellis–Székelyhidi). Whether Leray–Hopf solutions with the energy inequality are unique is still open and is perhaps the central problem in the weak solution theory.

4. Self-Similar and Discretely Self-Similar Solutions #

Self-similar solutions of the form $u(x,t) = (T^*-t)^{-1/2} U(x/(T^*-t)^{1/2})$ satisfy a nonlinear elliptic system for the profile $U$. Several non-existence theorems show that backward self-similar solutions with certain integrability must be trivial (Nečas–Ružička–Šverák, 1996). The case of discretely self-similar solutions, where $u(x,t) = \lambda u(\lambda x, \lambda^2 t)$ for a fixed $\lambda \neq 1$, is less understood and was recently revisited. Whether the set of self-similar profiles that could appear as blowup limits is empty is not known.

5. Computer-Assisted Proofs via Rigorous Numerics #

The Chen–Hou approach to Euler singularities (2025) used a computer-assisted proof framework: construct a numerical approximate profile, then verify its stability rigorously using interval arithmetic. For Navier–Stokes the presence of viscosity complicates such an approach (the profile is dissipated rather than transported), but the same framework (dynamical rescaling plus nonlinear stability verification) might in principle detect or rule out singularities in specific axi-symmetric geometries. Applying and adapting the Hou group’s methods to the viscous problem is an active direction.

6. The Zero-Viscosity Limit and Euler–Navier–Stokes Connection #

As $\nu \to 0$, Navier–Stokes formally converges to Euler. The precise relationship is subtle: in the presence of boundaries (Prandtl layers) or after a potential Euler singularity, the zero-viscosity limit can fail to hold in strong norms. If Euler develops a finite-time singularity at time $T^*_E$ from smooth data (as Chen–Hou suggest for bounded domains), then for small $\nu$ the Navier–Stokes solution must either also develop a near-singularity or be regularised by viscosity before $T^*_E$. Whether viscosity is always sufficient to regularise an Euler singularity, or whether a Navier–Stokes singularity can arise from a nearby Euler one, is entirely open.

References #

Fefferman, C. L. (2000). Existence and smoothness of the Navier–Stokes equation. Clay Mathematics Institute Millennium Prize Problems. https://www.claymath.org/wp-content/uploads/2022/06/navierstokes.pdf
Leray, J. (1934). Sur le mouvement d’un liquide visqueux emplissant l’espace. Acta Mathematica, 63, 193–248.
Hopf, E. (1951). Über die Anfangswertaufgabe für die hydrodynamischen Grundgleichungen. Mathematische Nachrichten, 4(1–6), 213–231.
Caffarelli, L., Kohn, R., & Nirenberg, L. (1982). Partial regularity of suitable weak solutions of the Navier–Stokes equations. Communications on Pure and Applied Mathematics, 35(6), 771–831.
Ladyzhenskaya, O. A. (1967). On uniqueness and smoothness of generalized solutions to the Navier–Stokes equations. Zapiski Nauchnykh Seminarov LOMI, 5, 169–185.
Escauriaza, L., Seregin, G. A., & Šverák, V. (2003). $L_{3,\infty}$-solutions of the Navier–Stokes equations and backward uniqueness. Russian Mathematical Surveys, 58(2), 211–250.
Tao, T. (2019). Quantitative bounds for critically bounded solutions to the Navier–Stokes equations. arXiv:1908.04958. Published in Nine Mathematical Challenges, AMS, 2021, pp. 149–193.
Tao, T. (2016). Finite time blowup for an averaged three-dimensional Navier–Stokes equation. Journal of the American Mathematical Society, 29(3), 601–674.
Buckmaster, T. & Vicol, V. (2019). Nonuniqueness of weak solutions to the Navier–Stokes equation. Annals of Mathematics, 189(1), 101–144.
Barker, T. & Prange, C. (2021). Localized quantitative estimates and potential blow-up rates for the Navier–Stokes equations. Communications in Mathematical Physics, 385, 717–792.

Navier–Stokes Regularity: The Uniqueness of Weak Solutions

Fri, 29 May 2026 00:00:00 +0000

The companion post on Navier–Stokes existence and smoothness asked whether smooth solutions can break down in finite time. This post asks the opposite question: when a solution is only weakly defined, satisfying the equations in an integral sense rather than pointwise, is it uniquely determined by its initial data? The answer, developed over the last two decades through a dramatic series of results, is a resounding no in many regimes. The frontier is now whether the physically natural class of Leray–Hopf weak solutions retains uniqueness.

Question (Weak Uniqueness)

Are Leray–Hopf weak solutions of the 3D incompressible Navier–Stokes equations $$\partial_t u + (u\cdot\nabla)u - \nu\Delta u + \nabla p = 0, \qquad \nabla\cdot u = 0$$ uniquely determined by their initial data $u_0 \in L^2(\mathbb{R}^3)$?

The question is one of the most urgent open problems in the PDE theory of fluid dynamics. It is logically independent of the blowup question: Leray–Hopf solutions exist globally for all time regardless of whether smooth solutions break down. What is not known is whether two Leray–Hopf solutions started from the same data must coincide.

Nash’s h-Principle: The Conceptual Ancestor #

The story begins not in fluid mechanics but in differential geometry. In 1954, John Nash proved that any Riemannian manifold admits a $C^1$ isometric embedding into Euclidean space, a result that contradicted the expectation, based on the rigid behaviour of $C^2$ embeddings (Cauchy), that the metric should impose strong constraints. The key insight is that $C^1$ embeddings are flexible: one can deform them by adding high-frequency oscillations that are invisible at the large scale but locally produce any prescribed metric tensor.

Gromov formulated this phenomenon as the h-principle: for certain underdetermined differential relations, the topological (homotopy-theoretic) obstructions are the only ones, and any formal solution can be deformed into an actual solution. The h-principle is a flexibility result: it says geometry is surprisingly unconstrained below a critical regularity threshold.

De Lellis and Székelyhidi recognised in the mid-2000s that the incompressible Euler equations are formally analogous to Nash’s embedding problem. The Euler system is underdetermined (more unknowns than equations), and one can attempt to construct wild solutions by adding high-frequency oscillations. The crucial observation is that the nonlinearity $u\otimes u$ in the Reynolds stress tensor plays the role of the metric tensor in Nash’s problem.

Wild Euler Solutions #

The first step was to show that the Euler equations possess infinitely many weak solutions for given initial data.

Theorem (De Lellis–Székelyhidi, 2009–2013)

For any divergence-free $u _0 \in L^2(\mathbb{T}^3)$ and any prescribed energy profile $e(t) \in C^\infty([0,T])$ with $e(t) > |u _0| _{L^2}^2$ for all $t > 0$, there exist infinitely many weak solutions $u \in C_t^0 L_x^2$ of the 3D Euler equations with $u(\cdot,0) = u _0$ and $|u(\cdot,t)| _{L^2}^2 = e(t)$.

In particular, the Euler equations admit weak solutions that spontaneously gain or lose kinetic energy for no reason: wild solutions. The construction proceeds by convex integration: one builds the solution iteratively, at each stage adding a high-frequency perturbation (a Beltrami wave) that corrects the error in the momentum equation while staying nearly invisible in the velocity field.

Earlier, Scheffer (1993) and Shnirelman (1997) had shown the existence of weak Euler solutions with compact support in space-time: the fluid is at rest, then spontaneously moves, then returns to rest; but their constructions were indirect. De Lellis and Székelyhidi’s convex integration scheme gave the first systematic and quantitative approach.

Onsager’s Conjecture #

The De Lellis–Székelyhidi results raise an immediate question: at what regularity does the fluid behaviour transition from flexible (wild, non-unique) to rigid (energy-conserving, unique)? This is precisely what Lars Onsager conjectured in 1949.

Onsager's Conjecture (1949)

For the 3D incompressible Euler equations, the threshold regularity for energy conservation is the Hölder exponent $1/3$:

If $u \in C^{0,\alpha}$ with $\alpha > 1/3$, then every weak solution conserves kinetic energy.
For every $\alpha < 1/3$, there exist weak solutions in $C^{0,\alpha}$ that dissipate energy.

The positive direction (conservation above $1/3$) was proved by Constantin–E–Titi (1994). The negative direction (dissipation possible below $1/3$) required much more work and was fully resolved only recently.

Theorem (Isett, 2018)

For every $\alpha < 1/3$ there exist weak solutions $u \in C^{0,\alpha}(\mathbb{T}^3\times[0,T])$ of the 3D Euler equations that fail to conserve kinetic energy.

Isett’s proof, published in the Annals of Mathematics in 2018, was the culmination of a decade of refinements of the De Lellis–Székelyhidi scheme. The key difficulty at regularity exactly $1/3$ is that the high-frequency perturbations must be sized to cancel the Reynolds stress error while staying in $C^{1/3-}$; this requires a delicate interplay of oscillation and concentration (intermittency). De Lellis, Székelyhidi, Buckmaster, and Vicol also obtained solutions attaining any prescribed energy profile in $C^{1/3-}$. Onsager’s conjecture is now a theorem.

Viscous Non-Uniqueness: Buckmaster–Vicol #

Adapting the convex integration scheme from Euler to Navier–Stokes requires overcoming the viscous term $\nu\Delta u$, which smooths out high-frequency oscillations. The intermittent Beltrami waves used by Isett concentrate energy at sparse spatial sets, reducing their interaction with the Laplacian. Buckmaster and Vicol exploited this idea to bring convex integration into the viscous setting.

Theorem (Buckmaster–Vicol, 2019)

There exist infinitely many weak solutions $u \in C_t^0 L_x^2(\mathbb{T}^3)$ of the 3D Navier–Stokes equations, belonging to the same regularity class as Leray–Hopf solutions, that do not satisfy the global energy inequality. In particular, weak solutions of 3D Navier–Stokes are not unique in the class $C_t^0 L_x^2$.

The Buckmaster–Vicol solutions, published in the Annals of Mathematics 189 (2019), 101–144, are weak in both the PDE sense and the energy sense: they satisfy the equations distributionally and have finite kinetic energy, but they can gain energy spontaneously, violating the natural dissipation law $\partial _t|u| _{L^2}^2 \leq -2\nu|\nabla u| _{L^2}^2$.

This non-uniqueness is striking but also limited: the Buckmaster–Vicol solutions are not Leray–Hopf solutions, because Leray–Hopf solutions are required to satisfy the energy inequality $|u(t)| _{L^2}^2 \leq |u _0| _{L^2}^2$. Whether this single additional constraint, that energy does not increase, suffices to restore uniqueness is the open question.

Crossing the Energy Barrier: Albritton–Brué–Colombo #

The energy inequality distinguishing Leray–Hopf solutions from Buckmaster–Vicol wild solutions seemed for a long time to be a genuine barrier to non-uniqueness. The following result crossed this barrier, but required introducing an external force.

Theorem (Albritton–Brué–Colombo, 2022)

There exists a body force $f \in L^1(0,T;, L^2(\mathbb{R}^3))$ and two distinct Leray–Hopf weak solutions of the forced 3D Navier–Stokes equations $\partial_t u + (u\cdot\nabla)u - \nu\Delta u + \nabla p = f$ with the same initial data $u_0 \equiv 0$ and the same force $f$.

Published in the Annals of Mathematics 196 (2022), 415–455, the proof uses a completely different mechanism from convex integration. The key ingredient is an unstable background solution: using Vishik’s construction of spectrally unstable steady states of the 2D Euler equations, Albritton–Brué–Colombo lift a 2D unstable vortex ring to an axisymmetric 3D solution and embed it into the Navier–Stokes flow via a self-similar change of variables. The force $f$ is chosen precisely to make this background exactly solve the forced equations; the instability then allows two different solutions to branch from the same initial data.

The force is singular; it belongs to $L^1_t L^2_x$ but is not smooth, and is concentrated near the initial time $t=0$. Whether the same non-uniqueness can be achieved with a smooth or zero force is the remaining open problem.

The Unforced Case: Current Frontier #

Non-uniqueness of Leray–Hopf solutions for the unforced Navier–Stokes equations remains open. The route to the unforced case requires finding a self-similar background profile that solves the unforced equations exactly and has an unstable eigenvalue, a far more demanding task than the forced case, where the profile can be any divergence-free function.

Open Problem (Jia–Šverák Programme)

Do there exist two distinct Leray–Hopf solutions of the 3D Navier–Stokes equations with the same initial data and no external force?

Jia and Šverák (2013–2014) showed that non-uniqueness would follow from a spectral assumption: if there exists a forward self-similar Navier–Stokes solution whose linearised operator has an eigenvalue with positive real part, then Leray–Hopf solutions are non-unique. Guillod and Šverák (2017) provided compelling numerical evidence that such an unstable self-similar profile exists.

In September 2025, Giri and Kwon posted a preprint (arXiv:2509.25116) claiming a computer-assisted proof of the existence of an unstable self-similar profile for the unforced equations, which, via the Jia–Šverák mechanism, would establish non-uniqueness of Leray–Hopf solutions. The proof uses rigorous interval arithmetic to verify the existence of an unstable eigenvalue. As of this writing the preprint is under review by the community.

The Regularity Threshold #

The accumulated results suggest the following picture of the flexibility-rigidity dichotomy for the Euler and Navier–Stokes equations.

Regularity class	Euler	Navier–Stokes
$C^{0,\alpha}$, $\alpha < 1/3$	non-unique, dissipative (Isett 2018)	n/a
$C^{0,\alpha}$, $\alpha > 1/3$	energy-conserving (Constantin–E–Titi 1994)	n/a
$L^2$ (global energy inequality)	non-unique	open (unforced); non-unique forced (ABC 2022)
$L^\infty_t L^3_x$ (LPS regularity)	n/a	unique and smooth (ESS 2003)

The Leray–Hopf class sits precisely at the boundary where uniqueness is expected to break down but has not yet been proved to do so in the unforced case.

Research Directions #

1. Resolving the Jia–Šverák Spectral Condition #

The most direct path to unforced Leray–Hopf non-uniqueness is to rigorously confirm or refute the spectral condition of Jia–Šverák: find (or prove the nonexistence of) a forward self-similar Navier–Stokes profile with an unstable linearised eigenvalue. The 2025 Giri–Kwon computer-assisted preprint claims this is now done. If confirmed, the consequence is striking: Leray’s 1934 existence theorem cannot be supplemented by uniqueness, and the Navier–Stokes Cauchy problem is ill-posed in the Leray–Hopf class.

2. Selection Principles and Physical Solutions #

If Leray–Hopf solutions are indeed non-unique, a fundamental question becomes which solution is the physically correct one, the one observed in experiments and computed in simulations. Several selection criteria have been proposed: the vanishing viscosity limit of the Navier–Stokes solution as $\nu\to 0$ from above, entropy conditions analogous to those for hyperbolic conservation laws, and renormalisation group or statistical ensemble approaches motivated by turbulence theory. None of these has been rigorously validated as a selection criterion that distinguishes a unique Leray–Hopf solution from the others.

3. Sharp Regularity Thresholds for Navier–Stokes #

For Euler, Onsager’s conjecture identifies $C^{1/3}$ as the sharp regularity threshold for energy conservation. What is the analogous threshold for Navier–Stokes? The Buckmaster–Vicol solutions are in $C_t^0 L_x^2$ (very rough), while the Ladyzhenskaya–Prodi–Serrin class gives uniqueness. The precise exponent at which uniqueness breaks down, if it does, is not known. Determining the sharp Sobolev or Hölder regularity threshold for Navier–Stokes uniqueness, analogous to Onsager’s $1/3$, is a central open problem.

4. Uniqueness for Axisymmetric Initial Data #

A natural restricted problem is whether Leray–Hopf solutions with axisymmetric, swirl-free initial data are unique. Such data imposes a strong geometric constraint that eliminates most of the degrees of freedom available to convex integration. Partial results are known (e.g., global regularity for axisymmetric data without swirl is not proved but no counterexamples exist), but uniqueness in this class has not been established. If the Giri–Kwon instability is confirmed, understanding whether the instability mechanism survives axisymmetric perturbations is an immediate question.

5. Stochastic Regularisation #

There is a well-studied phenomenon, regularisation by noise, in which adding a stochastic forcing term to an ill-posed deterministic PDE restores well-posedness. For the Navier–Stokes equations, Hofmanová–Zhu–Zhu (2023) showed non-uniqueness persists even under multiplicative noise for certain body forces, by adapting the Albritton–Brué–Colombo construction. Whether a generic stochastic perturbation can restore uniqueness of Leray–Hopf solutions, and what the appropriate notion of “generic” should be, is a rich open direction combining convex integration with stochastic analysis.

References #

Nash, J. (1954). $C^1$ isometric imbeddings. Annals of Mathematics, 60(3), 383–396.
De Lellis, C. & Székelyhidi, L. (2009). The Euler equations as a differential inclusion. Annals of Mathematics, 170(3), 1417–1436.
De Lellis, C. & Székelyhidi, L. (2013). Dissipative continuous Euler flows. Inventiones Mathematicae, 193(2), 377–407.
Constantin, P., E, W., & Titi, E. S. (1994). Onsager’s conjecture on the energy conservation for solutions of Euler’s equation. Communications in Mathematical Physics, 165(1), 207–209.
Isett, P. (2018). A proof of Onsager’s conjecture. Annals of Mathematics, 188(3), 871–963.
Buckmaster, T. & Vicol, V. (2019). Nonuniqueness of weak solutions to the Navier–Stokes equation. Annals of Mathematics, 189(1), 101–144.
Buckmaster, T. & Vicol, V. (2019). Convex integration and phenomenologies in turbulence. EMS Surveys in Mathematical Sciences, 6(1–2), 1–88.
Albritton, D., Brué, E., & Colombo, M. (2022). Non-uniqueness of Leray solutions of the forced Navier–Stokes equations. Annals of Mathematics, 196(1), 415–455.
Jia, H. & Šverák, V. (2014). Local-in-space estimates near initial time for weak solutions of the Navier–Stokes equations and forward self-similar solutions. Inventiones Mathematicae, 196(1), 233–265.
Giri, V. & Kwon, H. (2025). Nonuniqueness of Leray–Hopf solutions to the unforced incompressible 3D Navier–Stokes equation. arXiv:2509.25116.

The Regularity Problem for the 3D Euler Equations

Fri, 29 May 2026 00:00:00 +0000

Leonhard Euler wrote down the equations governing the motion of an ideal incompressible fluid in 1757. Whether smooth solutions to these equations can develop a singularity in finite time, a point at which derivatives of the velocity blow up, has been an open problem ever since, and remains one of the central questions in mathematical fluid dynamics.

Problem (Euler Regularity)

Let $u_0 : \mathbb{R}^3 \to \mathbb{R}^3$ be a smooth, divergence-free initial velocity field with sufficient decay at infinity. Does the unique local smooth solution $u(x,t)$ to the 3D incompressible Euler equations $$\partial_t u + (u \cdot \nabla)u + \nabla p = 0, \qquad \nabla \cdot u = 0, \qquad u(\cdot,0)=u_0$$ remain smooth for all time $t > 0$?

The problem is rated L4 on UnsolvedMath, reflecting its depth, and is closely related to the Clay Millennium Prize Problem on the Navier–Stokes equations. The two questions are linked through the zero-viscosity limit, but neither implies the other.

The Equations and What Regularity Means #

The Euler equations express conservation of momentum (first equation) and incompressibility (second equation) for an inviscid fluid. The unknowns are the velocity field $u(x,t) \in \mathbb{R}^3$ and pressure $p(x,t) \in \mathbb{R}$; the pressure is determined implicitly by incompressibility via an elliptic equation.

Vorticity. The central quantity for singularity analysis is the vorticity $\omega = \nabla \times u$, which satisfies the vorticity equation $$\partial_t \omega + (u \cdot \nabla)\omega = (\omega \cdot \nabla)u.$$ The right-hand side, the vortex stretching term, is the essential source of difficulty. It creates a quadratic feedback: large $\omega$ produces large $(\omega \cdot \nabla)u$, which can further amplify $\omega$.

Local well-posedness. For $u_0 \in H^s(\mathbb{R}^3)$ with $s > 5/2$, there exists a unique smooth solution on a time interval $[0, T^*)$ for some $T^* > 0$ depending on $|u _0| _{H^s}$ (Kato, 1972). The question is whether $T^*$ can be taken equal to $+\infty$.

Why 2D is easy, 3D is not. In two dimensions the vortex stretching term $(\omega \cdot \nabla)u$ vanishes identically by antisymmetry. The scalar vorticity $\omega = \partial_1 u_2 - \partial_2 u_1$ is then simply transported along fluid particle paths without amplification, and $|\omega|_{L^\infty}$ is conserved. Global regularity in 2D follows immediately. In 3D no such conservation holds, and the problem is genuinely open.

The Beale–Kato–Majda Criterion #

The first major structural result reduces the regularity problem to a single quantity.

Theorem (Beale–Kato–Majda, 1984)

A smooth solution $u$ of the 3D Euler equations loses regularity at time $T^*$ if and only if $$\int _0^{T^*} |\omega(\cdot,t)| _{L^\infty(\mathbb{R}^3)}, dt = +\infty.$$ In particular, if the vorticity remains bounded in $L^\infty$ on $[0,T]$ for every finite $T$, the solution remains smooth globally.

The BKM criterion redirects the problem: one must show that the vorticity magnitude $|\omega|_{L^\infty}$ cannot accumulate to infinity in finite time. Since $\omega$ satisfies a transport-stretching equation, this requires understanding the geometric structure of the vorticity field under its own evolution.

Geometric Conditions and Depletion of Stretching #

The vortex stretching term $(\omega \cdot \nabla)u$ can be decomposed as $$(\omega \cdot \nabla)u = |\omega|^2 (\hat\omega \cdot \nabla)\hat u,$$ where $\hat\omega = \omega/|\omega|$ is the unit vorticity direction. The key observation is that stretching is governed not only by the magnitude of $\omega$ but also by the geometry of the vorticity field.

Theorem (Constantin–Fefferman–Majda, 1996)

If the unit vorticity direction $\hat\omega = \omega/|\omega|$ is uniformly Lipschitz in a neighbourhood of the set ${|\omega| > \lambda}$ for all $t \in [0, T]$ and some $\lambda > 0$, then the solution remains smooth on $[0,T]$.

This result says that blowup, if it occurs, must be accompanied by violent geometric irregularity of vortex lines, not just large vorticity magnitude, but also loss of Lipschitz regularity of the vorticity direction. It has motivated a line of research on the geometric structure of vortex tubes near potential singularities.

Blowup for Less Regular Data #

Recent years have seen dramatic progress on singularity formation for initial data that is smooth except at isolated points.

Theorem (Elgindi, 2021)

There exist axisymmetric, swirl-free initial velocity fields $u_0 \in C^{1,\alpha}(\mathbb{R}^3)$ for sufficiently small $\alpha > 0$ such that the corresponding solution to the 3D Euler equations develops a finite-time singularity.

Elgindi’s proof, published in the Annals of Mathematics 194 (2021), 647–727, constructs a self-similar blowup profile and establishes its nonlinear stability using a dynamical rescaling formulation. The initial data is not smooth: it belongs to $C^{1,\alpha}$ but not to $C^2$. The singularity forms at the axis of symmetry $r=0$.

This was a breakthrough, but it left open the smooth case. Elgindi himself noted the next target: constructing blowup from initial data that is non-smooth only at a single point, or eventually from fully smooth data.

Extending Elgindi’s construction. Chen and Hou (2022) proved the same type of $C^{1,\alpha}$ blowup for the 3D axisymmetric Euler equations with boundary (inside a periodic cylinder), realising the Hou–Luo blowup scenario numerically proposed in 2014. Subsequent work by Córdoba, Martínez-Zoroa, and Zheng (2025, Annals of PDE) showed that the singularity can be formed from initial data in $C^\infty(\mathbb{R}^3 \setminus {0}) \cap C^{1,\alpha}$, with non-smoothness at a single point, a further step toward the smooth case.

The 2025 Breakthrough: Smooth Blowup with Boundary #

The most significant recent development is the following result, which provides a rigorous proof of finite-time singularity from smooth initial data.

Theorem (Chen–Hou, PNAS 2025)

There exists a family of smooth, finite-energy initial data for the 3D axisymmetric Euler equations in a smooth bounded domain (periodic cylinder) such that the corresponding solutions develop a finite-time singularity. The blowup is nearly self-similar and occurs at the intersection of the boundary $r=1$ and the symmetry plane $z=0$.

The paper, contributed by Thomas Hou and published in PNAS in June 2025 (reviewed by Caflisch, Gómez-Serrano, Sverak, and Tao), provides a computer-assisted proof. The strategy is to:

construct a numerical approximate self-similar blowup profile via the dynamical rescaling formulation,
prove rigorously that the true solution remains close to this profile using energy estimates with carefully verified error bounds (computed with interval arithmetic), and
conclude nonlinear stability of the blowup via a bootstrap argument.

This resolves the problem affirmatively in the setting of smooth data and a smooth bounded domain. The boundary plays a crucial role: it creates an antisymmetric flow pattern driving azimuthal vorticity toward a critical ring, generating intense vortex stretching at a hyperbolic saddle point on the wall.

The remaining open case. The problem in $\mathbb{R}^3$ (or on the periodic torus $\mathbb{T}^3$) without boundary remains open. It is not known whether smooth initial data in free space can produce a singularity, or whether the absence of a boundary provides a genuine stabilising mechanism.

Research Directions #

1. Removing the Boundary #

The most pressing open question is whether the Chen–Hou construction can be extended to $\mathbb{R}^3$ or $\mathbb{T}^3$. The boundary in the 2025 result acts as a geometric catalyst: it enforces a no-flow condition that concentrates vorticity at a specific ring on the wall. Without a boundary, the antisymmetric flow structure that drives the singularity must be sustained entirely by the initial data and the nonlinear dynamics. Whether a comparable mechanism can persist in free space, without the reflective constraint of the wall, is the central open question.

2. Self-Similar Blowup in Full 3D #

All current singularity results are for axisymmetric flows, which reduce the problem from 3 spatial dimensions to 2 (the $rz$-plane). In full 3D, the angular variable $\theta$ is active, and perturbations in the azimuthal direction can either stabilise or destabilise the singularity. Elgindi, Ghoul, and Masmoudi (2021) proved stability of the $C^{1,\alpha}$ blowup under axisymmetric perturbations. Whether the singularity survives fully 3D (non-axisymmetric) perturbations, a question Elgindi posed as open, is crucial: a blowup that is destroyed by any non-symmetric perturbation has limited physical relevance.

3. Quantitative Vortex Stretching and the Role of Geometry #

The BKM criterion and the Constantin–Fefferman–Majda theorem both express the same idea from opposite directions: blowup is controlled by the magnitude and geometry of the vorticity. Current research asks whether a quantitative version can be made sharp. Specifically: if the vorticity direction $\hat\omega$ becomes Hölder-continuous but not Lipschitz, does blowup necessarily follow? Or is there a finer scale invariant quantity, perhaps involving the Hessian of the velocity or the curvature of vortex lines, that governs the problem?

4. Weak Solutions and Non-Uniqueness #

Separate from the question of whether smooth solutions blow up is the question of what happens after a potential singularity. De Lellis and Székelyhidi (2009–2013) proved that the Euler equations have infinitely many weak $L^\infty$ solutions for generic initial data, via convex integration. Isett (2018) proved that weak solutions can dissipate energy, confirming Onsager’s 1949 conjecture. These results show that the solution concept must be carefully chosen. After a smooth blowup, the system likely enters a regime of non-unique weak solutions, and identifying the physically relevant selection criterion, entropy conditions, vanishing viscosity, $h$-principle, is a major open problem.

5. Vanishing Viscosity and the Navier–Stokes Connection #

The Navier–Stokes equations add a viscous term $\nu \Delta u$ to the right-hand side. For any $\nu > 0$, global regularity of Navier–Stokes in 3D is itself open (the Clay Millennium Problem). For the zero-viscosity limit $\nu \to 0$, the central question is whether Navier–Stokes solutions converge to Euler solutions uniformly in time, a question tied to boundary layer behaviour (the Prandtl conjecture) and to the regularity of the Euler solution. If Euler develops a singularity at time $T^*$, the behaviour of Navier–Stokes solutions near $T^*$ as $\nu \to 0$ is completely unknown.

References #

Euler, L. (1757). Principes généraux du mouvement des fluides. Mémoires de l’Académie des Sciences de Berlin, 11, 274–315.
Beale, J. T., Kato, T., & Majda, A. (1984). Remarks on the breakdown of smooth solutions for the 3-D Euler equations. Communications in Mathematical Physics, 94(1), 61–66.
Constantin, P., Fefferman, C., & Majda, A. J. (1996). Geometric constraints on potentially singular solutions for the 3-D Euler equations. Communications in Partial Differential Equations, 21(3–4), 559–571.
Elgindi, T. M. (2021). Finite-time singularity formation for $C^{1,\alpha}$ solutions to the incompressible Euler equations on $\mathbb{R}^3$. Annals of Mathematics, 194(3), 647–727.
Elgindi, T. M., Ghoul, T.-E., & Masmoudi, N. (2021). On the stability of self-similar blow-up for $C^{1,\alpha}$ solutions to the incompressible Euler equations. Cambridge Journal of Mathematics, 9(4), 1035–1075.
Chen, J. & Hou, T. Y. (2023). Finite time blowup of 2D Boussinesq and 3D Euler equations with $C^{1,\alpha}$ velocity and boundary. Communications in Mathematical Physics, 383, 4827–4890.
Chen, J. & Hou, T. Y. (2025). Singularity formation in 3D Euler equations with smooth initial data and boundary. Proceedings of the National Academy of Sciences, 122(27). https://doi.org/10.1073/pnas.2500940122
Córdoba, D., Martínez-Zoroa, L., & Zheng, F. (2025). Finite time singularities to the 3D incompressible Euler equations for solutions in $C^\infty(\mathbb{R}^3\setminus{0})\cap C^{1,\alpha}\cap L^2$. Annals of PDE. https://doi.org/10.1007/s40818-025-00214-2
Isett, P. (2018). A proof of Onsager’s conjecture. Annals of Mathematics, 188(3), 871–963.
Majda, A. J. & Bertozzi, A. L. (2002). Vorticity and Incompressible Flow. Cambridge University Press.