Why the Navier-Stokes Problem Is Hard

Many of the equations people first meet in physics are linear: double the input and the response doubles. Navier-Stokes is not like that.

Navier-Stokes? Nonlinear. The fluid's velocity affects its own rate of change, which means the fluid pushes itself. Imagine trying to predict where a crowd will go when every single person's movement depends on what everyone around them is doing, and what those people are doing depends on everyone around them, spiraling outward forever. That's the situation you're staring at.

The culprit is the self-interaction term $(u \cdot \nabla)u$. It creates feedback loops where small disturbances amplify into large ones, and it's why fluid turbulence is so wildly complex (see subproblems for more).

The convective nonlinearity $(u \cdot \nabla)u$ is the fundamental obstacle. In the vorticity formulation $\omega = \nabla \times u$, the equation becomes

$$\partial_t \omega + (u \cdot \nabla)\omega = (\omega \cdot \nabla)u + \nu \Delta \omega$$

The vortex stretching term $(\omega \cdot \nabla)u$ has no sign. It can amplify vorticity without bound. In 2D this term vanishes (since $\omega$ is a scalar perpendicular to the flow), which is why 2D global regularity is known (Ladyzhenskaya, 1969). In 3D, vortex stretching is the primary candidate mechanism for finite-time blowup.

Here's the crux: the nonlinearity is quadratic in $u$. The $H^1$ energy estimate gives $\|\nabla u\|_{L^2}$, but controlling $(u \cdot \nabla)u$ in $L^2$ typically needs stronger information such as $u \in L^\infty$ together with $\nabla u \in L^2$, or equivalent critical-scale control; the energy class alone does not provide this.

The Navier-Stokes equations have a scaling symmetry. Zoom in on a solution, make everything smaller and faster by the right amounts, and you get another perfectly valid solution. That symmetry is mathematically natural, but analytically dangerous.

Why? The only quantity we can reliably control is the total energy of the fluid, and it sits at completely the wrong scale, telling us about the big picture but saying absolutely nothing about what's happening at the microscopic scales where a blowup would actually form.

Think of monitoring a city's total electricity usage to detect one spark. Useful? Sure. Fine-grained enough? Not even close. That gap is the whole problem.

Under the natural scaling $u_\lambda(x,t) = \lambda u(\lambda x, \lambda^2 t)$, the critical Sobolev space is $\dot{H}^{1/2}(\mathbb{R}^3)$ (equivalently $L^3$). A norm is scale-critical if it is invariant under this transformation. If a controlled norm becomes smaller under zoom-in rescaling, the equation is supercritical relative to that control: the bound is below the scale needed to see concentration.

The energy inequality gives control of $u \in L^\infty_t L^2_x \cap L^2_t \dot{H}^1_x$. At the energy level, the 3D problem is supercritical:

$$\|u_\lambda\|_{L^2} = \lambda^{-1/2} \|u\|_{L^2}, \quad \|u_\lambda\|_{L^2_t \dot{H}^1_x} = \lambda^{-1/2} \|u\|_{L^2_t \dot{H}^1_x}$$

As $\lambda \to \infty$ (zooming in), these bounds weaken relative to the invariant quantities $\|u\|_{L^3}$ and $\|u\|_{\dot{H}^{1/2}}$. The nonlinear and viscous terms scale together; the problem is that the known a priori estimate sits below the critical scale. Getting from the energy class to a critical norm is the central difficulty.

Watch a river. Fluid motion goes chaotic. Turbulent. Big eddies shatter into smaller ones, which shatter into even smaller ones, cascading all the way down to microscopic scales where viscosity finally smooths things out.

Kolmogorov described this energy cascade in 1941, and the Navier-Stokes equations capture it beautifully. But here's what keeps people up at night: what if energy concentrates into smaller and smaller regions faster than viscosity can dissipate it? That's a blowup.

Can it actually happen? Or does viscosity always win? That's the open question, full stop. For the physical bridge from Reynolds number to this small-scales picture, see Reynolds Number, Turbulence, and Why Small Scales Matter.

Kolmogorov's K41 theory predicts an energy spectrum $E(k) \sim \varepsilon^{2/3} k^{-5/3}$ in the inertial range $k_f \ll k \ll k_\eta$, where $k_\eta \sim (\varepsilon/\nu^3)^{1/4}$ is the Kolmogorov dissipation wavenumber. Constant energy flux across scales. That's the clean picture.

The regularity question asks something darker: can this cascade degenerate? Can the dissipation scale $k_\eta^{-1}$ shrink to zero in finite time, with $\|\nabla u\|_{L^2} \to \infty$ while total energy stays finite?

Onsager's dissipation anomaly conjecture (1949) says yes, in a sense: in the vanishing viscosity limit $\nu \to 0$, energy dissipation persists, and weak solutions of Euler can dissipate energy. Confirmed for Hölder exponents below $1/3$ (Isett, 2018; Buckmaster et al., 2018). But what this means for Navier-Stokes regularity remains genuinely unclear. For regime-level intuition, see Reynolds Number, Turbulence, and Why Small Scales Matter.

Pressure in the Navier-Stokes equations is strange. It's not an independent variable at all; the velocity completely determines it through a single constraint: the fluid is incompressible, so it can't be squeezed.

This makes pressure nonlocal. In the incompressible model, changing the velocity in one region affects the pressure field globally, because the pressure is determined by the whole velocity field at that time.

For analysis, that's devastating. You can't study what happens at a single point without accounting for the entire fluid at once. Local reasoning? Forget it.

The incompressibility constraint $\nabla \cdot u = 0$ determines the pressure through the Poisson equation

$$-\Delta p = \nabla \cdot ((u \cdot \nabla)u) = \partial_i \partial_j (u_i u_j)$$

so $p = (-\Delta)^{-1} \partial_i \partial_j (u_i u_j)$, involving Riesz transforms (singular integral operators). The pressure is a nonlocal function of velocity, and this nonlocality is the core obstacle to pointwise or local-in-space estimates.

Standard maximum principle arguments fail here: even though the viscous term $\nu \Delta u$ is dissipative, the pressure gradient $-\nabla p$ can concentrate energy from distant regions. The Caffarelli-Kohn-Nirenberg theory handles this via local energy inequalities on parabolic cylinders, but extracting pointwise regularity from these remains the hard step.

For the 2D incompressible Navier-Stokes equations in the standard settings, global smooth solutions are known; this was established in classical work including Ladyzhenskaya's (1969). See why 2D is easier.

Three dimensions? Everything falls apart, and the reason comes down to one mechanism: vortex stretching. In 2D, vortices can spin and merge but they can't stretch. In 3D, fluid can grab vortex tubes and pull them thinner and thinner and thinner, potentially concentrating every last bit of energy into an infinitely thin filament.

Can this stretching run away to infinity in finite time, or does viscosity always step in? That's the million-dollar question. Literally.

The dichotomy is sharp.

2D: Vorticity $\omega$ is a scalar satisfying $\partial_t \omega + u \cdot \nabla \omega = \nu \Delta \omega$. The maximum principle gives $\|\omega(t)\|_{L^\infty} \leq \|\omega_0\|_{L^\infty}$, BKM implies global regularity, and the vortex stretching term $(\omega \cdot \nabla)u$ is identically zero.

3D: Vorticity $\omega \in \mathbb{R}^3$ satisfies $\partial_t \omega + (u \cdot \nabla)\omega = (\omega \cdot \nabla)u + \nu \Delta \omega$, where the stretching term $(\omega \cdot \nabla)u$ can amplify $|\omega|$ superlinearly (formally $\sim |\omega|^2$ via Biot-Savart). No maximum principle available.

The enstrophy $\|\omega\|_{L^2}^2$ satisfies

$$\frac{d}{dt}\|\omega\|_{L^2}^2 \leq C\|\omega\|_{L^2}^2 \|\nabla u\|_{L^\infty} - 2\nu \|\nabla \omega\|_{L^2}^2$$

Controlling $\|\nabla u\|_{L^\infty}$ requires $\omega \in L^\infty$, which requires controlling $\|\nabla u\|_{L^\infty}$. Circular. No existing technique has broken this loop.

This article is part of The Problem.

These obstacles have led mathematicians to decompose the problem into subproblems and develop specialized approaches for each.

For context on why the viscous term helps but isn't enough, see Euler vs. Navier-Stokes.

Exact solutions exist for symmetric flows, and engineers solve the equations numerically every day. Why neither settles the problem is the subject of Solving the Navier-Stokes Equations.

This article is part of The Problem.

The regularity problem breaks into tractable pieces; see Subproblems. The analytical tools built to attack these obstacles are covered in Approaches. Why does the viscous term help but still not close the gap? Euler vs. Navier-Stokes. For why closed-form and numerical solvability do not imply well-posedness, see Solving the Navier-Stokes Equations.