Approaches to the Navier-Stokes Problem

The most fundamental approach uses energy. A moving fluid has kinetic energy, and viscosity dissipates it — like friction slowing things down. The total energy can only decrease over time (assuming no external forcing).

This was Leray's key insight in 1934: use the energy bound to prove that some kind of solution must exist. His method constructs approximate solutions (with artificial smoothing), proves they all satisfy the energy bound, and then takes a limit.

The limitation: the energy bound is too coarse to guarantee smoothness. It tells you the fluid has finite total energy, but not that the velocity stays finite everywhere.

Paper links: Leray, Sur le mouvement d'un liquide visqueux emplissant l'espace (1934); Hopf, Über die Anfangswertaufgabe für die hydrodynamischen Grundgleichungen (1950).

The Leray-Hopf construction proceeds via Galerkin approximation or mollification. Key steps:

1. Approximate: Solve the mollified system $\partial_t u_\varepsilon + (J_\varepsilon u_\varepsilon \cdot \nabla) u_\varepsilon = \nu \Delta u_\varepsilon - \nabla p_\varepsilon$ on finite-dimensional subspaces.
2. Energy bound: The a priori estimate $\|u_\varepsilon(t)\|_{L^2}^2 + 2\nu \int_0^t \|\nabla u_\varepsilon\|_{L^2}^2 \leq \|u_0\|_{L^2}^2$ holds uniformly in $\varepsilon$.
3. Compactness: Extract a weakly convergent subsequence $u_\varepsilon \rightharpoonup u$ in $L^2_t \dot{H}^1_x$ using the Aubin-Lions lemma.
4. Pass to limit: The nonlinear term converges by strong $L^2_{\text{loc}}$ convergence of $u_\varepsilon$.

The resulting weak solution satisfies the energy inequality (not equality — energy may be lost at irregular times). The gap between the energy class $L^\infty_t L^2_x \cap L^2_t \dot{H}^1_x$ and smoothness is precisely the regularity problem.

Paper links: Leray, Sur le mouvement d'un liquide visqueux emplissant l'espace (1934); Hopf, Über die Anfangswertaufgabe für die hydrodynamischen Grundgleichungen (1950).

The Caffarelli-Kohn-Nirenberg approach (1982) doesn't try to prove full smoothness. Instead, it asks: how bad can the singularities be?

The answer: remarkably tame. Their $\varepsilon$-regularity theorem says that if the fluid's energy is sufficiently small in a small space-time region, then the solution is smooth there. Since the total energy is finite, there simply aren't enough "budget" for many singular points.

This is like proving that a wall might have cracks, but the total length of all cracks combined is zero — they can only be isolated points.

Paper links: Caffarelli-Kohn-Nirenberg, Partial regularity of suitable weak solutions of the Navier-Stokes equations (1982); Albritton-Barker-Prange, Epsilon regularity for the Navier-Stokes equations via weak-strong uniqueness.

The CKN $\varepsilon$-regularity theorem: there exists $\varepsilon_{\text{CKN}} > 0$ such that if $(u,p)$ is a suitable weak solution and

$$\frac{1}{r} \int_{Q_r(z_0)} |\nabla u|^2 \, dx \, dt < \varepsilon_{\text{CKN}}$$

then $u$ is regular (Hölder continuous) at $z_0 = (x_0, t_0)$. Here $Q_r(z_0) = B_r(x_0) \times (t_0 - r^2, t_0)$ is a parabolic cylinder.

The proof combines the local energy inequality with a Campanato-type iteration: if the scale-invariant energy is small, a boot-strap argument shows $u$ is bounded, then Hölder, then smooth by classical Schauder theory.

The dimensional estimate $\mathcal{P}^1(\Sigma) = 0$ follows by a Vitali covering: if $\Sigma$ had positive $\mathcal{P}^1$ measure, infinitely many disjoint parabolic cylinders would each carry $\varepsilon_{\text{CKN}}$ energy, contradicting finite total energy.

Paper links: Caffarelli-Kohn-Nirenberg, Partial regularity of suitable weak solutions of the Navier-Stokes equations (1982); Albritton-Barker-Prange, Epsilon regularity for the Navier-Stokes equations via weak-strong uniqueness.

A classic continuation criterion, first proved by Beale, Kato, and Majda for the 3D Euler equations and later adapted in several forms to Navier-Stokes, says that blowup can only happen if vorticity control is lost.

Vorticity measures how much the fluid is spinning locally. The message of the BKM-type criteria is: if you can keep the maximum spin under control in the right norm, then the solution continues smoothly. Other dangerous quantities are then forced to stay controlled as well.

This reduced the problem to a single family of quantities. Unfortunately, controlling them has proved just as hard as the original problem.

Paper links: Beale-Kato-Majda, Remarks on the breakdown of smooth solutions for the 3-D Euler equations (1984); Kozono-Taniuchi, Bilinear estimates in BMO and the Navier-Stokes equations (2000).

The original Beale-Kato-Majda theorem (1984) is for the 3D Euler equations. For Navier-Stokes, analogous continuation criteria imply that a smooth solution $u$ on $[0, T^*)$ extends beyond $T^*$ whenever

$$\int_0^{T^*} \|\omega(\cdot, t)\|_{L^\infty} \, dt < \infty,$$

where $\omega = \nabla \times u$ is the vorticity. Refinements include:

• Kozono-Taniuchi (2000): $\|\omega\|_{L^\infty}$ can be replaced by $\|\omega\|_{\mathrm{BMO}}$
• Besov-space variants: critical or borderline Besov control can also serve as a continuation criterion
• Direction-restricted criteria: Da Veiga (1995) showed that certain scale-invariant bounds on $\nabla u$ already suffice

These criteria connect to the vortex-stretching picture: any finite-time singularity must force the vorticity to accumulate too quickly for the time integral above to stay finite.

Paper links: Beale-Kato-Majda, Remarks on the breakdown of smooth solutions for the 3-D Euler equations (1984); Kozono-Taniuchi, Bilinear estimates in BMO and the Navier-Stokes equations (2000); Chemin-Planchon, Self-improving bounds for the Navier-Stokes equations (2012).

A modern approach works with special mathematical spaces that sit right at the boundary of what the scaling symmetry allows. These are called critical spaces.

The idea: if you can show that a solution stays within certain critical-space bounds, smoothness follows automatically. Several teams have established this, creating a menu of "regularity criteria" — conditions that, if verified, guarantee smoothness.

The challenge remains going from what we can prove (subcritical bounds from energy) to what we need (critical bounds). This gap is narrow but has resisted all attempts to bridge it.

Paper links: Kenig-Koch, An alternative approach to regularity for the Navier-Stokes equations in critical spaces; Gallagher-Koch-Planchon, A profile decomposition approach to the $L^\infty_t(L^3_x)$ Navier-Stokes regularity criterion.

Major programs in critical-space regularity:

• Koch-Tataru (2001): Global well-posedness for small data in $\text{BMO}^{-1}$, the largest critical space where the bilinear estimate $\|\mathbb{P}\nabla \cdot (u \otimes v)\|_{\text{BMO}^{-1}} \lesssim \|u\|_{\text{BMO}^{-1}} \|v\|_{\text{BMO}^{-1}}$ holds. This is essentially optimal for perturbative methods.
• Gallagher-Koch-Planchon (2013): Profile decomposition for Navier-Stokes in $\dot{H}^{1/2}$. Any sequence of solutions with bounded critical norm has a subsequence decomposing into asymptotically decoupled profiles.

The fundamental obstruction: no known coercive functional is both controlled by the evolution and critical with respect to the Navier-Stokes scaling.

Paper links: Kenig-Koch, An alternative approach to regularity for the Navier-Stokes equations in critical spaces; Gallagher-Koch-Planchon, A profile decomposition approach to the $L^\infty_t(L^3_x)$ Navier-Stokes regularity criterion.

Modern PDE theory uses tools from harmonic analysis — the mathematics of breaking functions into waves at different frequencies (like a musical chord into individual notes).

By decomposing the fluid velocity into components at different spatial scales and tracking how energy moves between scales, mathematicians can make the vague intuition of "energy cascade" precise. These techniques, called Littlewood-Paley decomposition, have produced the sharpest known results about regularity criteria and blowup rates.

Paper links: Cannone-Meyer, Littlewood-Paley decomposition and Navier-Stokes equations (1995); Gallagher-Koch-Planchon, A profile decomposition approach to the $L^\infty_t(L^3_x)$ Navier-Stokes regularity criterion.

Littlewood-Paley theory decomposes $u = \sum_j \Delta_j u$ where $\Delta_j$ localizes to frequencies $|\xi| \sim 2^j$. For Navier-Stokes:

The paraproduct decomposition of the nonlinearity $(u \cdot \nabla)u$ splits into low-high, high-low, and high-high frequency interactions:

$$(u \cdot \nabla)u = T_u \nabla u + T_{\nabla u} u + R(u, \nabla u)$$

where $T$ is the paraproduct and $R$ the remainder. Each piece has different regularity properties in Besov spaces $\dot{B}^s_{p,q}$.

Key results using this machinery:

• Chemin-Lerner spaces: $\widetilde{L}^\rho_T \dot{B}^s_{p,q}$ provide the natural framework for critical well-posedness: the Navier-Stokes bilinear form maps $\widetilde{L}^\infty_T \dot{B}^{-1+3/p}_{p,q} \times \widetilde{L}^1_T \dot{B}^{1+3/p}_{p,q} \to \widetilde{L}^1_T \dot{B}^{-1+3/p}_{p,q}$.
• Cannone-Meyer: Littlewood-Paley methods give a clean wavelet/Besov formulation of the small-data theory.

Paper links: Cannone-Meyer, Littlewood-Paley decomposition and Navier-Stokes equations (1995); Gallagher-Koch-Planchon, A profile decomposition approach to the $L^\infty_t(L^3_x)$ Navier-Stokes regularity criterion.

A less traditional but increasingly powerful approach uses the geometry of the flow. Instead of tracking numbers (norms, energies), these methods study the shape of the solution — how vortex tubes bend, how regions of intense rotation are arranged in space.

The insight is that blowup isn't just about something getting big — it's about the fluid organizing itself in a very specific geometric configuration. If you can show that configuration is impossible (because it leads to a contradiction with, say, energy conservation or incompressibility), you've ruled out blowup.

This is one of the styles of argument explored by the proof notes on this site.

Paper links: Constantin, Geometric statistics in turbulence (1994); Albritton-Barker-Prange, Localized smoothing and concentration for the Navier-Stokes equations in the half space.

Geometric-topological approaches to regularity exploit structural constraints that are invisible to purely analytical methods:

• Vortex line geometry: Constantin (1994) showed that if the vorticity direction field $\hat{\omega} = \omega/|\omega|$ is Lipschitz in regions of high vorticity, the solution is regular. Blowup requires the vorticity direction to develop a singularity simultaneously with the magnitude.
• Incompatibility arguments: if a blowup configuration is geometrically constrained (e.g., via packing bounds on the number of independent concentration regions that fit within energy and dissipation budgets), one can derive a contradiction without directly estimating critical norms.
• Case partition: by classifying each spatial region as belonging to one of finitely many scenarios (e.g., locally regular, Type-I-like, Type-II-like, densely packed) and showing each scenario either gives regularity or transfers the problem to a bounded counting argument, the blowup can be ruled out combinatorially.

The proof page on this site is intended to explore arguments of this flavor; it is not presented here as a completed formal proof of the Millennium problem.

Paper links: Constantin, Geometric statistics in turbulence (1994); Albritton-Barker-Prange, Localized smoothing and concentration for the Navier-Stokes equations in the half space.

Each of these approaches has made progress, but none has solved the problem completely. See our proof notes for in-progress notes on a geometric approach that aims to combine several of these techniques.

Or go back to the basics: what exactly is the Millennium Problem?

For in-progress notes on a geometric-topological strategy that aims to combine CKN ε-regularity, vorticity criteria, and packing arguments, see Proof Notes.