Vorticity

Definition

The vorticity field ω is defined as the curl of the velocity field u:

ω = ∇ × u

In Cartesian components, with u = (u, v, w):

ω_x = ∂w/∂y − ∂v/∂z
ω_y = ∂u/∂z − ∂w/∂x
ω_z = ∂v/∂x − ∂u/∂y

Vorticity is a vector. Its direction is the axis of local rotation (right-hand rule), and its magnitude is the rate of spin. The crucial physical interpretation is that vorticity equals twice the angular velocity Ω of a fluid parcel rotating as a rigid body:

ω = 2Ω

The factor of two trips up almost everyone the first time. It arises because the antisymmetric part of the velocity-gradient tensor — the rotation rate — contributes to both off-diagonal terms in the curl. A solid body spinning at Ω = 1 rad/s carries vorticity 2 s⁻¹ everywhere inside it.

Rotation is not the same as going in circles

The single most common misconception is that fluid moving in curved paths must have vorticity, and fluid moving in straight lines must not. Both halves are wrong. Vorticity is about local rotation of a parcel, not the curvature of its trajectory. Imagine freezing a tiny cross — two perpendicular material lines — painted on the fluid:

Flow	Velocity profile	Streamlines	Does the cross rotate?	Vorticity
Rigid (forced) vortex	v = Ωr (rises with r)	circles	yes, at Ω	ω = 2Ω (uniform)
Free (potential) vortex	v = Γ/(2πr) (falls with r)	circles	no — only orbits	ω = 0 (except center)
Uniform shear flow	u = γ̇·y (straight lines)	parallel lines	yes	ω_z = −γ̇ ≠ 0
Uniform translation	u = const	parallel lines	no	ω = 0

The free vortex is the showpiece. Streamlines are perfect circles, yet the flow is irrotational everywhere except the singular core. A paddle wheel placed off-center has its inner edge in faster fluid and its outer edge in slower fluid; the two torques cancel exactly, so the wheel orbits the center without ever turning. Conversely, simple shear flow has dead-straight streamlines but nonzero vorticity, because the top of the cross is dragged forward faster than the bottom.

The vorticity transport equation

Taking the curl of the incompressible Navier–Stokes equations eliminates the pressure term entirely (the curl of a gradient is zero) and yields the vorticity transport equation:

Dω/Dt = (ω·∇)u + ν ∇²ω

where D/Dt = ∂/∂t + (u·∇) is the material derivative and ν is the kinematic viscosity. The two terms on the right carry all the physics:

• (ω·∇)u — vortex stretching and tilting. When a vortex tube is stretched along its axis, conservation of angular momentum forces it to spin faster, amplifying vorticity. This term is absent in two-dimensional flow, which is precisely why 2D and 3D turbulence behave so differently — 3D turbulence cascades energy to small scales through vortex stretching.
• ν ∇²ω — viscous diffusion. Vorticity spreads out and decays, much like heat. This is the only term that can create or destroy vorticity in the bulk; in an inviscid fluid it vanishes.

For a compressible or stratified fluid, a third term appears — the baroclinic source (∇ρ × ∇p)/ρ² — which generates vorticity whenever surfaces of constant density and constant pressure are misaligned. This drives sea breezes, thunderstorm outflows, and the Rayleigh–Taylor instability.

Circulation and Stokes' theorem

Vorticity is local; its global partner is circulation Γ, the line integral of velocity around a closed loop:

Γ = ∮_C u · dl = ∬_S (∇ × u) · dA = ∬_S ω · dA

Stokes' theorem makes vorticity literally "circulation per unit area." Two theorems built on this dominate inviscid vortex dynamics:

• Kelvin's circulation theorem. In an inviscid, barotropic flow with conservative forces, the circulation around any material loop is constant in time: DΓ/Dt = 0. Vorticity cannot appear from nothing.
• Helmholtz's theorems. Vortex lines move with the fluid, the strength of a vortex tube is constant along its length, and a vortex tube cannot end inside the fluid — it must close on itself (a smoke ring) or terminate at a boundary.

Numerical examples

System	Vorticity / circulation
Stirred coffee, rigid rotation at 2 rev/s	Ω ≈ 12.6 rad/s → ω ≈ 25 s⁻¹
Atmospheric synoptic flow (mid-latitudes)	ω ~ 10⁻⁵ – 10⁻⁴ s⁻¹
Planetary (Coriolis) vorticity at 45° latitude	f = 2Ω sin φ ≈ 1.0×10⁻⁴ s⁻¹
Tornado core (~50 m/s at 30 m radius)	ω ~ 3 s⁻¹ (10⁵× ambient)
Boeing 747 wingtip vortex circulation	Γ ≈ 500 m²/s; lingers minutes
Bathtub vortex (lab, controlled)	Coriolis effect real but ~10⁻³× the residual swirl

Vorticity and lift

The reason aerospace engineers care so much about vorticity is the Kutta–Joukowski theorem: the lift per unit span on any 2D body equals the fluid density times the freestream speed times the circulation bound to the body,

L' = ρ V Γ

A lifting wing carries bound vorticity. Because Kelvin's theorem demands the total circulation stay zero, the wing sheds an equal and opposite starting vortex the instant it begins to move, and the bound vorticity trails off the tips as a pair of counter-rotating wingtip vortices. These trailing vortices are responsible for induced drag and for the wake turbulence that forces aircraft to maintain spacing.

JavaScript — computing vorticity on a grid

// 2D vorticity ω_z = ∂v/∂x − ∂u/∂y on a regular grid via central differences.
// u, v are Nx-by-Ny arrays (u[i][j]); dx, dy are grid spacings.
function vorticity2D(u, v, dx, dy) {
  const Nx = u.length, Ny = u[0].length;
  const w = Array.from({ length: Nx }, () => new Float64Array(Ny));
  for (let i = 1; i < Nx - 1; i++) {
    for (let j = 1; j < Ny - 1; j++) {
      const dvdx = (v[i + 1][j] - v[i - 1][j]) / (2 * dx);
      const dudy = (u[i][j + 1] - u[i][j - 1]) / (2 * dy);
      w[i][j] = dvdx - dudy;
    }
  }
  return w;
}

// Free (potential) vortex: v_theta = Γ / (2π r). Sample its velocity field...
const Gamma = 1.0, dx = 0.05, dy = 0.05, N = 81, c = (N - 1) / 2;
const U = [], V = [];
for (let i = 0; i < N; i++) {
  U.push(new Float64Array(N)); V.push(new Float64Array(N));
  for (let j = 0; j < N; j++) {
    const x = (i - c) * dx, y = (j - c) * dy;
    const r2 = x * x + y * y;
    if (r2 < 1e-6) continue;                 // skip the singular core
    const vth = Gamma / (2 * Math.PI * Math.sqrt(r2));
    U[i][j] = -vth * (y / Math.sqrt(r2));    // tangential -> Cartesian
    V[i][j] =  vth * (x / Math.sqrt(r2));
  }
}

const W = vorticity2D(U, V, dx, dy);
// Away from the center the curl is ~0 even though streamlines are circles:
console.log('off-center vorticity:', W[20][40].toExponential(2)); // ~0 (machine epsilon)

// Rigid rotation v_theta = Ω r gives uniform vorticity ω = 2Ω instead:
const Omega = 1.5;
for (let i = 0; i < N; i++) for (let j = 0; j < N; j++) {
  const x = (i - c) * dx, y = (j - c) * dy;
  U[i][j] = -Omega * y;  V[i][j] = Omega * x;
}
const Wrigid = vorticity2D(U, V, dx, dy);
console.log('rigid-rotation vorticity:', Wrigid[40][40].toFixed(3), '≈ 2Ω =', (2 * Omega).toFixed(3));

Where vorticity shows up

• Aerodynamics. Bound vorticity sets lift; trailing wingtip vortices set induced drag and wake-turbulence spacing rules (heavy jets need up to 8 km separation).
• Meteorology. Absolute vorticity (relative + planetary) is conserved on isobaric surfaces; potential-vorticity maps are the workhorse of weather forecasting. Tornado and hurricane cores are intense vorticity concentrations.
• Oceanography. Eddies, the Gulf Stream meanders, and Rossby waves are all governed by potential-vorticity conservation on the rotating Earth.
• Turbulence. Vortex stretching cascades energy from large eddies to small dissipative scales; vorticity is the natural variable for the Kolmogorov picture.
• Sports and propulsion. The Magnus effect on a spinning ball, the leading-edge vortex on a swift's wing, and fish-fin thrust all manipulate shed vorticity.
• Computation. Vortex methods and vorticity–streamfunction formulations solve incompressible flow without ever computing pressure explicitly.

Common mistakes

• Equating circular streamlines with vorticity. The free vortex circulates but is irrotational. Curvature of a path and rotation of a parcel are independent.
• Forgetting the factor of two. Vorticity is 2Ω, not Ω. Mixing the two halves your spin rate.
• Assuming straight flow is vorticity-free. Any shear (a transverse velocity gradient) carries vorticity even with perfectly straight streamlines.
• Ignoring vortex stretching in 3D. The (ω·∇)u term has no 2D analog; without it you cannot explain why 3D turbulence drains energy to small scales or why tornado funnels intensify.
• Believing vorticity is conserved like mass. Only circulation around a material loop is conserved (and only in inviscid, barotropic flow). Vorticity magnitude changes by stretching; viscosity diffuses and destroys it.
• Overstating the bathtub Coriolis myth. Earth's rotation does impart vorticity, but in a basin it is roughly a thousand times weaker than residual swirl from filling — it does not reliably set the drain direction.