Navier–Stokes energy and weak solutions
Derive viscous dissipation, construct the logic of Leray–Hopf compactness, and distinguish global energy-class existence from strong regularity.
On this page
The question for today
What you’ll understand
- Derive the unforced smooth Navier–Stokes energy equality with the viscosity coefficient.
- State a standard global Leray–Hopf energy class on the three-dimensional torus.
- Explain why weak convergence needs compactness to pass the nonlinear product and why an energy inequality does not imply smoothness.
Viscosity adds a rigorous energy budget to incompressible motion. That budget is strong enough to build global weak solutions through approximation and compactness, even when high derivatives are not controlled. It is not the same as a global smooth-solution theorem. Understanding the gap requires tracking the topology of convergence, the nonlinear product, and the exact norm supplied by dissipation. A simple shear solution makes the coefficients tangible without pretending every nonlinear fluid evolves by scalar heat flow.
Viscosity dissipates gradients
For unforced Navier–Stokes with positive viscosity nu, add nu times the Laplacian of velocity to the Euler right-hand side. On a periodic torus, dot the smooth equation with velocity and integrate. Incompressible transport and pressure give no net contribution, while integration by parts of the Laplacian gives minus nu times the squared gradient norm. With kinetic energy defined as half the squared L² norm, the coefficient is nu, not two nu; the latter appears only when the equality is written for the squared norm itself.
The integrated identity bounds velocity uniformly in L² and its gradient in L² over every finite time interval. In two dimensions these bounds interact with stronger estimates to support global smooth evolution from smooth data. In three dimensions the energy budget alone does not bound the high regularity used to continue a classical solution. The distinction is analytic, not a failure of the formal cancellation: the cancellation is exact, but its controlled quantity is too weak to supply all desired derivative bounds.
Compactness produces an energy-class solution
On the three-dimensional torus, take divergence-free L² initial data. A standard Leray–Hopf existence theorem gives a global distributional solution that is weakly continuous in L², belongs to L-infinity in time with values in L², and to L² in time with values in H¹ on every finite interval. It attains the initial data strongly in L² and satisfies the energy inequality from the initial time, together with the standard almost-every-starting-time form. The pressure can be recovered distributionally; divergence-free test fields remove it from the basic weak formulation.
One construction projects onto finitely many divergence-free Fourier modes. The finite-dimensional equations retain the energy cancellation, giving bounds uniform in the projection size. A time-derivative bound in a suitable negative space, together with spatial compactness on the torus, yields strong local-in-time L² convergence along a subsequence. Mere weak convergence of u would not identify u tensor u. Strong convergence of the needed factors lets the nonlinear term pass to the limit, while lower semicontinuity turns the approximate energy identities into an inequality.
Keep the solution classes distinct
A weak existence theorem supplies neither uniqueness nor regularity merely by producing a limit. A weak–strong uniqueness theorem says that, while a sufficiently regular strong solution exists, an energy-class Leray–Hopf solution with the same data agrees with it. This compares classes on an interval where the strong hypotheses hold; it does not extend that interval by definition. Additional spacetime bounds such as Prodi–Serrin control can bridge the regularity gap, and their exponents will be checked in the next lesson.
The setting also matters for research claims. External forcing contributes a work term to the energy balance, and a bounded-domain problem adds boundary conditions that affect estimates. Whole-space and periodic formulations have related but distinct compactness and pressure issues. The equations in this lesson are explicitly unforced and periodic. A dated announcement about another regularity or breakdown formulation should not be substituted for their energy-class theorem, nor should the energy-class theorem be mistaken for an assessment of that announcement. The logical content here is global weak existence under the stated assumptions and the limits of the resulting bounds.
Work it through
Check the viscous budget on an exact shear
On [0,2pi]³, set nu=1/2 and u(t,x,y,z)=(2 exp(−t/2) sin y,0,0). Verify the unforced equation with constant pressure and calculate its energy and dissipation rate.
The divergence and nonlinear transport vanish because the sole nonzero component is independent of x. Its Laplacian is minus itself, so u_t=−u/2=nu Delta u, verifying the full unforced equation.
Integrating its squared magnitude gives 16pi³ exp(−t), hence kinetic energy E(t)=8pi³ exp(−t). The gradient squared integral has the same value 16pi³ exp(−t), since cosine squared and sine squared integrate equally.
Multiplying that gradient integral by nu gives 8pi³ exp(−t), exactly minus E prime. At time 2 log two, the velocity amplitude is one and the energy is 2pi³, one quarter of its initial value.
The exact shear obeys the smooth energy equality. It reduces to heat flow because its nonlinear transport vanishes, a special property not shared by general solutions.
The living notebook
Explore the model
Interactive lab · heat
Match an exact viscous shear mode
Pause and reset, then set both controls to n = 1 and κ = 0.5. At time one, predict normalized shear amplitude exp(−0.5); use two Steps. For a fresh comparison, pause and reset, then set n = 2 and κ = 0.5, predicting exp(−2) at the same time. For velocity with one component proportional to sin(ny), independent of its own flow direction, these are the exact amplitude factors of unforced periodic Navier–Stokes shear solutions. The instrument displays a scalar heat mode, not general nonlinear transport or pressure coupling.
- Current amplitude
- 1
- Spatial L² norm
- 1.772
- Decay rate · κn²
- 3.6 /s
The exact periodic heat-equation mode is u(x,t) = exp(−κn²t) sin(nx), on [0,2π]. Playback advances physical time. The initial sine is shown separately, and the unnormalized L² norm is √π exp(−κn²t). This analytic instrument does not run a finite-difference solver or test discretization stability. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.
Look forHeat decay exactly describes this special shear family, while the general energy inequality allows nonlinear transfer between scales.
Pause & reason
Check your understanding
Your turn
Try it yourself
An unforced Leray–Hopf solution has initial kinetic energy ten, viscosity 1/4, and kinetic energy four at a later permitted time T. Bound the time integral of its squared gradient norm up to T. Can this bound determine the supremum norm of the gradient at time T?
Take a moment to make a prediction or work it out on paper. Open a hint when you need one.
Give me a hint
Use the energy inequality with the same half-norm convention, and distinguish time-integrated L² control from a pointwise-in-time L∞ quantity.
Compare with the solution
The inequality gives four plus one quarter of the gradient integral at most ten. Thus the integral is at most twenty-four. Equality would follow for a smooth solution satisfying the energy identity, but cannot be assumed from the weak inequality. This integrated spatial L² bound does not determine a gradient supremum at T, or even a uniform high Sobolev bound throughout the interval. Short, spatially concentrated gradient growth is not excluded by that single finite budget.
Keep the idea
Take it with you
- The smooth viscous energy identity gives velocity L² control and time-integrated gradient L² control.
- Compactness must identify nonlinear products when constructing a weak limit.
- Leray–Hopf existence, weak–strong uniqueness, and global classical regularity are distinct statements.
Read further
Original sources and deeper explanations to follow at your own pace.
Optional progress, saved in this browser when available. Revisit any lesson at any time.
Lesson not marked complete.