4.5 Kelvin's circulation theorem

Kelvin's circulation theorem overview

Kelvin's circulation theorem tells you something powerful: in an inviscid, barotropic flow, the circulation around a closed material curve doesn't change over time. This result constrains how vorticity can be created and transported, and it underpins much of classical aerodynamics and geophysical fluid dynamics.

The theorem connects three ideas you already know (velocity fields, vorticity, and material derivatives) into a single conservation statement. From here, you can explain why starting vortices form behind airfoils, why large-scale ocean eddies persist, and what conditions actually can generate new vorticity.

Circulation in fluid dynamics

Definition of circulation

Circulation, $\Gamma$, is a scalar quantity that measures the net rotational motion of fluid around a closed curve. Think of it as asking: "If you walked along a closed loop in the flow, how much does the velocity field push you along that path?"

It's a macroscopic measure of rotation, as opposed to vorticity, which captures rotation at a point.

Circulation as a line integral

Circulation is defined as the line integral of velocity around a closed curve $C$:

$\Gamma = \oint_C \vec{v} \cdot d\vec{l}$

where $\vec{v}$ is the velocity field and $d\vec{l}$ is the infinitesimal line element along $C$.

• Positive $\Gamma$ corresponds to counterclockwise rotation (by the right-hand rule convention).
• Negative $\Gamma$ corresponds to clockwise rotation.
• The value of $\Gamma$ depends on which fluid particles the curve encloses. Two different curves enclosing the same fluid will give the same circulation only if they bound the same vorticity distribution (via Stokes' theorem).

Circulation and vorticity

Circulation and vorticity are linked through Stokes' theorem. Vorticity is the local measure of rotation, defined as the curl of velocity:

$\vec{\omega} = \nabla \times \vec{v}$

Stokes' theorem lets you convert the line integral into a surface integral:

$\Gamma = \oint_C \vec{v} \cdot d\vec{l} = \iint_S \vec{\omega} \cdot d\vec{S}$

where $S$ is any surface bounded by $C$. So circulation equals the total vorticity flux through the enclosed area. If there's no vorticity inside the curve, the circulation is zero.

Mathematical formulation

Kelvin's circulation theorem equation

The theorem states that for a closed material curve $C(t)$ (one that moves with the fluid), the circulation doesn't change in time:

$\frac{D\Gamma}{Dt} = 0$

Here $\frac{D}{Dt}$ is the material derivative, meaning you're tracking the same fluid particles as they move. This holds under three conditions:

1. The flow is inviscid (no viscous forces).
2. The flow is barotropic (pressure is a function of density alone, $p = p(\rho)$).
3. Body forces are conservative (derivable from a potential, like gravity).

If all three conditions are met, whatever circulation the material curve starts with, it keeps forever.

Proof sketch

The derivation starts from Euler's equation for inviscid flow and computes $\frac{D}{Dt}\oint_C \vec{v} \cdot d\vec{l}$. You need to differentiate both $\vec{v}$ and $d\vec{l}$ (since the curve itself is deforming). The pressure gradient term vanishes when $p = p(\rho)$ because $\nabla p / \rho$ becomes the gradient of a scalar (the pressure function $\int dp/\rho$), and the integral of any exact gradient around a closed loop is zero. The body force term vanishes for the same reason if the force is conservative. The remaining terms from differentiating $d\vec{l}$ cancel exactly, leaving $\frac{D\Gamma}{Dt} = 0$.

Inviscid vs. viscous flows

• Inviscid flows: Circulation is conserved on material curves. Vorticity can only enter the flow through boundaries or through violations of the barotropic assumption.
• Viscous flows: Viscous diffusion (the $\nu \nabla^2 \vec{v}$ term in Navier-Stokes) adds a non-zero contribution to $\frac{D\Gamma}{Dt}$. Vorticity can be generated, diffused, and dissipated within the fluid interior.

The inviscid assumption is reasonable at high Reynolds numbers away from boundaries, which is why Kelvin's theorem is so useful in external aerodynamics and large-scale geophysical flows.

Barotropic vs. baroclinic flows

• Barotropic: $p = p(\rho)$, so surfaces of constant pressure and constant density are parallel. No new circulation is generated.
• Baroclinic: Pressure depends on both density and another variable (e.g., temperature), so $\nabla p$ and $\nabla \rho$ are misaligned. This misalignment produces a torque (the baroclinic torque, $\frac{\nabla \rho \times \nabla p}{\rho^2}$) that generates vorticity and changes circulation.

A classic example: sea breezes. Differential heating creates misaligned pressure and density gradients near a coastline, generating circulation that Kelvin's theorem (in its basic form) can't account for.

Physical interpretation

Conservation of circulation

The core physical message is straightforward: in an ideal barotropic fluid, rotation doesn't appear from nothing. If a material loop starts with zero circulation, it stays at zero. If it starts with some finite $\Gamma$, that value is locked in as the loop stretches, deforms, and advects with the flow.

This is a strong constraint. It means that any circulation you observe in an inviscid, barotropic flow must trace back to initial conditions or to boundary effects.

Relationship to vorticity

Through Stokes' theorem, conservation of circulation implies conservation of the vorticity flux through any material surface. In an inviscid barotropic flow:

• Vortex lines move with the fluid (they are "frozen in").
• Vorticity cannot be created or destroyed in the fluid interior.
• Vortex tubes maintain their strength (circulation) as they stretch and deform.

This "frozen-in" property of vorticity is one of the most important consequences of Kelvin's theorem. It's closely related to Helmholtz's vortex theorems, which describe how vortex lines and tubes behave in inviscid flow.

Role in fluid motion

Circulation and vorticity govern much of what you see in real flows:

• Vortices and eddies form and persist because of vorticity conservation.
• Lift on airfoils is directly proportional to circulation (via the Kutta-Joukowski theorem).
• Large-scale atmospheric and oceanic currents are shaped by vorticity dynamics and conservation principles that descend from Kelvin's theorem.

Applications of Kelvin's theorem

Vortex dynamics

Kelvin's theorem is the starting point for understanding how vortices behave:

• Aircraft wake vortices: The bound circulation on a finite wing sheds into trailing vortices. Kelvin's theorem requires that the total circulation in any large loop around the wing system remains constant, so the trailing vortices must carry circulation equal and opposite to changes in the bound circulation.
• Vortex rings: A smoke ring maintains its circulation as it propagates, consistent with Kelvin's theorem.
• Vortex shedding: Behind bluff bodies, vortices shed alternately from each side (von Kármán vortex street). The theorem constrains the total circulation budget.

Aerodynamics and lift generation

This is where Kelvin's theorem has its most famous consequence. When an airfoil starts from rest:

1. Initially, circulation around any material curve is zero (fluid at rest).
2. As the airfoil accelerates, the Kutta condition demands a specific circulation $\Gamma$ around the airfoil to keep the flow smooth at the trailing edge.
3. Kelvin's theorem requires total circulation to remain zero, so a starting vortex of strength $-\Gamma$ is shed from the trailing edge.
4. The bound circulation $\Gamma$ around the airfoil produces lift via the Kutta-Joukowski theorem: $L = \rho U \Gamma$ (per unit span).

The starting vortex is real and observable. It's a direct, physical consequence of Kelvin's theorem.

Geophysical fluid dynamics

On planetary scales, Kelvin's theorem (and its extensions) explain large-scale circulation patterns:

• Ocean eddies persist for months because the surrounding flow is approximately inviscid at large scales.
• Potential vorticity conservation (Ertel's theorem) extends Kelvin's ideas to rotating, stratified fluids and governs the behavior of atmospheric jet streams and Rossby waves.
• Ekman pumping and other boundary-layer effects are the mechanisms by which vorticity enters the large-scale flow, consistent with the theorem's requirement that vorticity generation needs viscosity or baroclinicity.

Limitations and extensions

Assumptions and limitations

Kelvin's theorem requires:

• Inviscid flow: No viscous stresses. Violated near walls and in boundary layers.
• Barotropic fluid: $p = p(\rho)$ only. Violated whenever temperature or composition gradients create misaligned $\nabla p$ and $\nabla \rho$.
• Conservative body forces: Forces like gravity (derivable from a potential) satisfy this. Coriolis force in a rotating frame requires careful treatment but doesn't violate the theorem when handled properly.

In real fluids, all three assumptions break down to some degree, especially near solid boundaries.

Viscous effects and dissipation

Viscosity modifies the circulation equation to:

$\frac{D\Gamma}{Dt} = \oint_C \nu \nabla^2 \vec{v} \cdot d\vec{l}$

This term is generally non-zero, meaning circulation on a material curve can increase or decrease. Physically, viscosity diffuses vorticity (spreading it out) and ultimately dissipates it into heat. Near walls, viscosity is the primary mechanism for injecting new vorticity into the flow.

Generalized theorems

Two important extensions build on Kelvin's result:

• Bjerknes circulation theorem: Accounts for baroclinic effects by adding a term proportional to the number of $(\rho, p)$ solenoids (regions where density and pressure contours intersect) enclosed by the curve. This is essential in meteorology for understanding cyclone formation.
• Ertel's potential vorticity theorem: Combines Kelvin's theorem with the thermodynamic equation in a rotating, stratified fluid. The conserved quantity is potential vorticity, $PV = \frac{\vec{\omega}_a \cdot \nabla \theta}{\rho}$, where $\vec{\omega}_a$ is absolute vorticity and $\theta$ is potential temperature. This is one of the most powerful tools in geophysical fluid dynamics.

Historical context

Lord Kelvin's contributions

William Thomson (Lord Kelvin) published the circulation theorem in 1869. It was part of a broader effort to understand vortex motion in ideal fluids, motivated partly by his (ultimately incorrect) vortex atom hypothesis, which proposed that atoms were knotted vortex tubes in the ether. Despite the wrong physical motivation, the mathematical results were lasting.

Development of circulation concepts

The theorem builds on earlier work:

• Leonhard Euler (1750s) established the equations of inviscid fluid motion.
• Hermann von Helmholtz (1858) proved his vortex theorems, showing that vortex lines move with the fluid and that vortex tubes conserve their strength in inviscid flow. Kelvin's theorem can be seen as a more general statement that encompasses Helmholtz's results.

Modern relevance

Kelvin's circulation theorem remains a foundational tool. It's used in:

• Computational vortex methods (tracking discrete vortex elements)
• Turbulence theory (understanding the cascade and conservation of enstrophy in 2D turbulence)
• Flow control and aerodynamic design (manipulating circulation for lift and drag optimization)

Numerical simulations now allow researchers to study how closely real viscous flows approximate the inviscid predictions of Kelvin's theorem, and where the deviations matter most.