research36 min studyLesson 428 in Mathematics

Incompressible Euler evolution and vorticity

Derive pressure and vorticity constraints, compare two-dimensional transport with three-dimensional stretching, and separate energy control from derivative control.

On this page

What you’ll understand

  • Derive the pressure Poisson equation from incompressibility.
  • Explain why scalar two-dimensional vorticity is transported while three-dimensional vorticity can stretch.
  • State the Sobolev threshold of a standard classical local theory without treating its endpoint as automatic.

An ideal incompressible fluid transports velocity while pressure continuously enforces the divergence constraint. The pressure is determined nonlocally, so a fluid parcel does not evolve as an isolated ordinary differential equation. Taking curl replaces pressure by a vorticity equation that exposes a major dimensional difference: in two dimensions the scalar rotation is transported, while in three dimensions the flow can stretch that rotation. Energy conservation remains compatible with increasingly sharp gradients, making these two kinds of control fundamentally different.

Pressure enforces a nonlocal constraint

For an unforced fluid of constant density, normalized to one, write velocity evolution as u_t plus u dot gradient u plus gradient pressure equals zero, with divergence of u zero. Work on a periodic torus to avoid boundary fluxes, or on whole space with sufficient decay for the integrations used. Taking divergence gives minus the Laplacian of pressure equal to the sum of partial_i u_j times partial_j u_i. Thus pressure is recovered by an elliptic equation, up to its irrelevant spatial constant in the periodic setting.

For divergence-free initial velocity in H to the s with s strictly greater than one plus d over two, the standard classical Sobolev local theory gives existence, uniqueness, and persistence in that regularity, together with appropriate time differentiability in one lower Sobolev order. The strict threshold makes velocity Lipschitz through Sobolev embedding, controlling transport. This statement covers dimensions two and three in the usual periodic or whole-space settings. It does not assert the same conclusion at the borderline exponent, where endpoint phenomena need separate analysis.

∂tu+(u⋅∇)u+∇p=0,∇⋅u=0,−Δp=∑i,j(∂iuj)(∂jui)\partial_tu+(u\cdot\nabla)u+\nabla p=0,\qquad \nabla\cdot u=0,\qquad -\Delta p=\sum_{i,j}(\partial_i u_j)(\partial_j u_i)
The pressure equation follows by differentiating the divergence constraint. Periodic pressure is determined up to a spatial constant when the compatibility condition is satisfied.↔ Scroll wide equations horizontally.

Curl reveals transport and stretching

In three dimensions define omega as curl u. Taking curl of the smooth Euler equation yields material derivative of omega equal to omega dot gradient u. The right-hand side stretches and tilts vorticity. Its sign is not fixed: a direction aligned with extending strain can grow, while compression in that direction can weaken it. The antisymmetric rotation part of the velocity gradient does not change the squared vorticity magnitude through this quadratic pairing; the symmetric strain supplies that change.

For a genuinely two-dimensional flow, vorticity is a scalar perpendicular to the plane and the stretching term vanishes. The scalar is transported by the incompressible flow, preserving its L to the p norms for smooth periodic solutions. This mechanism underlies the stronger global theory for smooth two-dimensional data. A planar projection of a three-dimensional flow does not remove stretching, and a numerical picture of nearly planar motion does not justify replacing the three-dimensional vorticity equation by a scalar transport law.

Dtω=(ω⋅∇)u(d=3),Dtω=0(d=2),Dt=∂t+u⋅∇D_t\omega=(\omega\cdot\nabla)u\quad(d=3),\qquad D_t\omega=0\quad(d=2),\qquad D_t=\partial_t+u\cdot\nabla
The three-dimensional equation uses vector vorticity; the two-dimensional equation uses its scalar out-of-plane component. Both are for smooth unforced incompressible Euler flow.↔ Scroll wide equations horizontally.

Conserved kinetic energy does not bound every derivative

Dotting the smooth Euler equation with velocity and integrating gives constant kinetic energy: transport is a divergence and pressure does no net work under the stated periodic or decaying conditions. This L² velocity quantity does not directly control the supremum of vorticity or the high Sobolev norm required by local theory. A family can keep its L² size while moving variation to shorter scales. In three dimensions, stretching offers a mechanism for such derivative amplification without contradicting the energy identity.

A continuation criterion identifies what would prevent or be necessary for breakdown in a specified regularity class. It must not be confused with a proof that the full equation realizes a proposed stretching scenario. Smoothness, finite Hölder regularity, external forcing, the whole-space versus periodic domain, and the presence of a physical boundary distinguish different singularity questions. Current research announcements should be read with those hypotheses intact. The calculations here establish the classical equations and local-control mechanisms; they do not classify every announced three-dimensional construction or turn a local strain example into a global singular solution.

A stationary periodic shear transports without changing

On the two-dimensional torus [0,2pi]², take u=(2 sin y,0). Verify incompressibility and stationary Euler evolution with constant pressure, calculate vorticity and kinetic energy, and follow one parcel.

  1. The divergence is zero because the first component has no x dependence and the second is zero. The nonlinear directional derivative is 2 sin y times partial_x u, also zero, so the time-independent field solves Euler with constant pressure.

  2. The scalar vorticity is partial_x u₂ minus partial_y u₁, equal to −2 cos y. Its supremum magnitude is two, and it is constant along each parcel because y does not change.

  3. The energy is one half times the integral of 4sin² y over the torus, giving 4pi². A parcel satisfies y(t)=y₀ and x(t)=x₀+2t sin y₀ modulo 2pi, illustrating shear transport without vortex stretching.

The result

This is one explicit globally smooth solution, not an argument that every three-dimensional Euler flow has the same behavior.

Explore the model

Interactive lab

Interactive lab · fields

Separate rotation from compressibility

Pause and reset, then set both controls to q = 0 and Ω = 1. The displayed local field is (−y,x); predict divergence zero and scalar curl two. For a fresh comparison, pause and reset, then set q = 0 and Ω = 2; predict divergence zero again and curl four. This is a local solid-rotation vector-field illustration, not the periodic shear above, a finite-energy whole-space Euler solution, or a pressure solver. In particular, larger rotation alone is not a demonstration of three-dimensional vortex stretching.

Separate rotation from compressibilityThe softened illustration is q·r/(r² + 0.2) + Ω(−y, x). Moving markers indicate local direction, not integrated particle paths or a Maxwell solution. Current readings: Fx at (1, 0): 0.833; Fy at (1, 0): 0.5; Field magnitude there: 0.972. All parameters are adjustable with the labeled controls below.−3.2−1.601.63.2−1.5−0.7500.751.5xyq = 1
Local field directionSource strength
Fx at (1, 0)
0.833
Fy at (1, 0)
0.5
Field magnitude there
0.972
t = 0 s

The softened illustration is q·r/(r² + 0.2) + Ω(−y, x). Moving markers indicate local direction, not integrated particle paths or a Maxwell solution. 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 forVorticity measures rotation, divergence measures volume change, and stretching is an additional three-dimensional mechanism.

Check your understanding

Why does smooth two-dimensional Euler vorticity obey a simpler equation than three-dimensional vorticity?

Try it yourself

At a point of a smooth three-dimensional incompressible flow, let the velocity gradient A, with entries A_ij=partial_j u_i, be [[−1,−1,0],[1,−1,0],[0,0,2]]. Compute divergence, vorticity, its instantaneous stretching vector, and the material derivative of half its squared magnitude.

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

First take the trace and antisymmetric curl combinations. Then use D_t omega=A omega, not the trace alone.

Compare with the solution

The trace is −1−1+2=0. Curl is (0,0,A₂₁−A₁₂)=(0,0,2). Multiplying gives A omega=(0,0,4). Therefore D_t(|omega|²/2)=omega dot A omega=8, a positive instantaneous amplification. This local algebra is consistent with incompressibility and the curl of the given gradient. It is not a construction of a globally stationary linear flow or a finite-time singular Euler solution; the full gradient and pressure must evolve according to the PDE.

Take it with you

  • Pressure is an elliptic nonlocal constraint, not an independent parcel force chosen arbitrarily.
  • Two-dimensional scalar vorticity is transported, while three-dimensional strain can stretch vector vorticity.
  • Energy, classical regularity, and current singularity results require different norms and exact hypotheses.

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.