1.6 Vorticity and circulation

Definition of vorticity

Vorticity quantifies the local rotation of a fluid element at any point in the flow. It connects directly to circulation, lift generation, and the behavior of boundary layers and wakes, making it one of the most important quantities in aerodynamics.

Mathematical representation

Vorticity is defined as the curl of the velocity field:

$\boldsymbol{\omega} = \nabla \times \mathbf{u}$

In Cartesian coordinates, this produces a vector with three components: $\omega_x$, $\omega_y$, and $\omega_z$. The magnitude of $\boldsymbol{\omega}$ tells you the strength of local rotation, and its direction gives the axis about which the fluid element spins.

For two-dimensional flows, vorticity reduces to a single scalar component ($\omega_z$ or simply $\omega$), since rotation can only occur about the axis perpendicular to the flow plane.

Physical interpretation

Think of vorticity as twice the angular velocity of a fluid element spinning about its own center. A small fluid parcel can translate through space and rotate simultaneously; vorticity captures only the rotational part, independent of the translational motion.

• Vorticity is a local property: it can vary from point to point within the flow.
• High-vorticity regions typically correspond to strong shear flows, boundary layers, and turbulent eddies.
• A flow with zero vorticity everywhere is called irrotational (more on this below).

Circulation in fluid dynamics

Circulation gives you a single number that captures the net rotational strength of the flow around a closed path. Where vorticity describes local rotation at a point, circulation describes the cumulative rotational effect over a finite region.

Definition and properties

Circulation, denoted $\Gamma$, is the line integral of velocity around a closed contour $C$:

$\Gamma = \oint_C \mathbf{u} \cdot d\mathbf{l}$

• $\Gamma$ is a scalar with SI units of $\text{m}^2/\text{s}$.
• Its value depends on what vorticity the contour encloses, not on the specific shape of the contour itself.
• In inviscid, barotropic flows, circulation is conserved (Kelvin's circulation theorem, discussed next).

Relationship to vorticity

Circulation and vorticity are linked through Stokes' theorem. The circulation around a closed contour equals the surface integral of vorticity over any surface $S$ bounded by that contour:

$\Gamma = \iint_S \boldsymbol{\omega} \cdot d\mathbf{S}$

In 2D, this simplifies to:

$\Gamma = \iint_A \omega \, dA$

This is the bridge between the local picture (vorticity at each point) and the global picture (total rotation around a region). If there's no vorticity inside a contour, the circulation around it is zero.

Kelvin's circulation theorem

Kelvin's circulation theorem states that in an inviscid, barotropic flow with only conservative body forces, the circulation around a closed material contour (one that moves with the fluid) does not change over time:

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

Here $\frac{D}{Dt}$ is the material derivative, tracking a specific set of fluid particles.

Conservation of circulation

The theorem requires three conditions:

1. Inviscid flow (no viscosity)
2. Barotropic fluid (density is a function of pressure only)
3. Conservative body forces (e.g., gravity, but no Coriolis or electromagnetic forces)

When these hold, the contour can stretch, deform, and twist as it moves with the fluid, but the circulation around it stays fixed.

Implications for vorticity

• Vorticity cannot be created or destroyed within the interior of an inviscid, barotropic fluid. It can only enter or leave through boundaries (solid walls, free surfaces).
• Vortex lines move with the fluid. A material line that starts as a vortex line remains a vortex line for all time.
• Vortex lines cannot terminate inside the fluid; they must close on themselves or end at boundaries.

These results explain why, in real flows, viscosity at solid surfaces is the primary source of new vorticity.

Vortex lines and tubes

Vortex lines and tubes are geometric tools for visualizing the structure of vorticity in a flow, much like streamlines visualize velocity.

Definitions and characteristics

• A vortex line is a curve that is everywhere tangent to the local vorticity vector at a given instant. Vortex lines cannot cross each other, since $\boldsymbol{\omega}$ has a unique direction at each point.
• A vortex tube is a bundle of vortex lines forming a tubular surface. The strength of a vortex tube equals the circulation $\Gamma$ computed around any cross-section of the tube.
• Physical examples of vortex tubes include wingtip vortices, tornado cores, and bathtub drain vortices.

Helmholtz's vortex theorems

Helmholtz formulated three theorems for inviscid, barotropic flows:

1. Constant strength: The circulation (strength) of a vortex tube is the same at every cross-section along its length.
2. No free ends: A vortex line cannot terminate inside the fluid. It must form a closed loop or extend to a boundary.
3. Material identity: Fluid particles that lie on a vortex line at one instant remain on a vortex line at all later times, and the tube's strength stays constant as it moves with the fluid.

These theorems are the geometric counterparts of Kelvin's circulation theorem and reinforce why vorticity is such a persistent feature of fluid flows.

Generation of vorticity

In real (viscous) fluids, vorticity is not conserved. It can be created, amplified, and destroyed. Understanding these mechanisms is critical for predicting boundary layer behavior, separation, and wake structure.

Role of viscosity

Viscosity is the primary agent of vorticity generation and dissipation:

• At solid boundaries, the no-slip condition forces the fluid velocity to match the wall velocity, creating steep velocity gradients and therefore vorticity.
• Away from walls, viscous diffusion spreads vorticity outward and smooths sharp gradients, analogous to how thermal diffusion spreads heat.
• The balance between generation at walls and diffusion into the flow determines the overall vorticity distribution.

Boundary layer vorticity

Boundary layers are thin regions near solid surfaces where viscous effects dominate.

• The steep velocity gradient normal to the surface (from zero at the wall to the freestream value) produces intense vorticity.
• This vorticity can be shed into the outer flow when the boundary layer separates from the surface.
• Boundary layer separation on airfoils and bluff bodies is a major source of vorticity entering the wake.

Vorticity in wakes and shear layers

• Wakes form downstream of objects and contain velocity deficits and elevated turbulence, both driven by the vorticity shed from the body.
• Shear layers are regions of high velocity gradient, found at wake edges or where two streams of different speed meet. Vorticity is concentrated here.
• Shear layers are prone to Kelvin-Helmholtz instability, where small perturbations grow into coherent, rolling vortical structures.

Vorticity equation

The vorticity equation governs how vorticity evolves in time and space. It is derived by taking the curl of the Navier-Stokes equations.

Derivation from Navier-Stokes equations

For an incompressible, Newtonian fluid:

$\frac{\partial \boldsymbol{\omega}}{\partial t} + (\mathbf{u} \cdot \nabla)\boldsymbol{\omega} = (\boldsymbol{\omega} \cdot \nabla)\mathbf{u} + \nu \nabla^2 \boldsymbol{\omega}$

Each term has a distinct physical meaning:

Transport and diffusion of vorticity

• The advection term $(\mathbf{u} \cdot \nabla)\boldsymbol{\omega}$ moves vorticity from place to place without changing its magnitude or orientation. It simply carries vorticity along with the flow.
• The viscous diffusion term $\nu \nabla^2 \boldsymbol{\omega}$ smooths out vorticity gradients over time. Higher kinematic viscosity $\nu$ means faster diffusion.

Vortex stretching and tilting

The term $(\boldsymbol{\omega} \cdot \nabla)\mathbf{u}$ is unique to three-dimensional flows (it vanishes identically in 2D) and has two effects:

• Vortex stretching: When the velocity gradient elongates a vortex tube along its axis, the tube narrows and spins faster (conservation of angular momentum). This amplifies vorticity.
• Vortex tilting: Velocity gradients perpendicular to a vortex line can reorient it, transferring vorticity from one component to another.

These mechanisms are the primary way turbulence generates small-scale, intense vortical structures from larger-scale motions.

Potential flow vs. rotational flow

These are two fundamentally different categories of fluid motion, distinguished by whether vorticity is present.

Irrotational vs. rotational flow

Irrotational (potential) flow:

• Vorticity is zero everywhere: $\boldsymbol{\omega} = \nabla \times \mathbf{u} = 0$
• Fluid elements translate and deform but do not rotate.
• The velocity field can be written as the gradient of a scalar potential: $\mathbf{u} = \nabla \phi$
• Useful for modeling flow far from solid surfaces where viscous effects are negligible.

Rotational flow:

• Vorticity is nonzero in at least part of the domain.
• Fluid elements undergo local rotation.
• The velocity field cannot be expressed purely as $\nabla \phi$.
• Found in boundary layers, wakes, and any region where viscous effects have introduced vorticity.

Velocity potential and stream function

In irrotational flow, two scalar functions simplify the analysis:

• Velocity potential $\phi$: velocity components are $u = \frac{\partial \phi}{\partial x}$, $v = \frac{\partial \phi}{\partial y}$, $w = \frac{\partial \phi}{\partial z}$. For incompressible flow, $\phi$ satisfies Laplace's equation $\nabla^2 \phi = 0$.
• Stream function $\psi$ (2D only): $u = \frac{\partial \psi}{\partial y}$, $v = -\frac{\partial \psi}{\partial x}$. Curves of constant $\psi$ are streamlines, the paths fluid particles follow in steady flow.

Both $\phi$ and $\psi$ satisfy Laplace's equation in 2D irrotational, incompressible flow, which is why potential flow problems are mathematically tractable.

Kutta-Joukowski theorem

The Kutta-Joukowski theorem is the central theoretical result connecting circulation to aerodynamic lift.

Lift generation by circulation

For a 2D airfoil in steady, inviscid, incompressible flow, the lift per unit span is:

$L' = \rho_\infty U_\infty \Gamma$

where $\rho_\infty$ is the freestream density, $U_\infty$ is the freestream velocity, and $\Gamma$ is the circulation around the airfoil.

How does circulation arise if the freestream is irrotational? Through two mechanisms:

1. Viscosity in the boundary layer generates vorticity at the airfoil surface.
2. The Kutta condition requires smooth flow departure at the sharp trailing edge, which fixes a unique value of $\Gamma$.

When the airfoil first starts moving, a starting vortex of equal and opposite circulation is shed from the trailing edge, so the total circulation in the flow remains zero (consistent with Kelvin's theorem applied to a large contour enclosing everything).

Applications in airfoil design

• Airfoil camber (curvature of the mean line) directly increases circulation and therefore lift at zero angle of attack.
• Increasing angle of attack also increases circulation, up to the point of stall.
• Designers tailor the camber and thickness distribution to achieve target lift while minimizing drag for specific flight regimes (subsonic, transonic, supersonic).
• The theorem also explains unsteady phenomena like dynamic stall, where rapid changes in angle of attack temporarily alter the circulation and shed additional vortices.

Vortex-induced vibrations

Vortex-induced vibrations (VIV) occur when periodic vortex shedding from a bluff body excites structural oscillations. VIV is a major design concern for bridges, offshore risers, chimneys, and heat exchanger tubes.

Mechanism of vortex shedding

1. Flow approaches a bluff body (e.g., a circular cylinder).
2. Boundary layers develop on the upper and lower surfaces.
3. The boundary layers separate near the rear of the body.
4. Separated shear layers roll up into discrete vortices that detach alternately from each side.
5. The result is a von Kármán vortex street: a staggered pattern of counter-rotating vortices in the wake.

Each vortex shedding event produces a fluctuating pressure force on the body, with a lateral (lift) component oscillating at the shedding frequency and a streamwise (drag) component oscillating at twice that frequency.

Strouhal number and lock-in phenomena

The Strouhal number relates shedding frequency to flow conditions:

$St = \frac{fD}{U}$

where $f$ is the shedding frequency, $D$ is the body's characteristic dimension (e.g., cylinder diameter), and $U$ is the freestream velocity. For a circular cylinder over a wide range of Reynolds numbers (roughly $300 < Re < 3 \times 10^5$), $St \approx 0.2$.

Lock-in occurs when the shedding frequency approaches a natural frequency of the structure:

• The shedding frequency synchronizes with the structural vibration frequency.
• Vibration amplitudes grow dramatically, sometimes reaching one diameter peak-to-peak for cylinders.
• The synchronized state persists over a range of flow velocities, not just a single speed.
• This can cause fatigue failure if not accounted for in design (e.g., through helical strakes or dampers).

Vorticity in turbulent flows

Turbulent flows contain chaotic, multi-scale motions, and vorticity is at the heart of their dynamics. Understanding vorticity in turbulence connects the small-scale dissipative structures to the large-scale energy-containing motions.

Vorticity dynamics in turbulence

• Vorticity in turbulence is concentrated in intense, slender structures often called vortex tubes or worms, with diameters on the order of the Kolmogorov length scale.
• These tubes are continuously stretched, tilted, and twisted by the surrounding velocity field, amplifying vorticity and generating ever-smaller scales.
• Complex interactions occur: vortex tubes can reconnect (break and rejoin with neighbors), merge, or break down into smaller structures.
• The vorticity field is highly intermittent: patches of very strong vorticity coexist with large quiescent regions.

Enstrophy and energy cascade

Enstrophy measures the total squared vorticity in the flow:

$\mathcal{E} = \frac{1}{2} \int \omega^2 \, dV$

Enstrophy is directly proportional to the rate of viscous dissipation of kinetic energy in incompressible flow, making it a key diagnostic quantity.

The energy cascade describes how turbulent kinetic energy flows through scales:

1. Energy enters at large scales (large eddies driven by mean flow instabilities or external forcing).
2. Large eddies break down into smaller eddies through vortex stretching and instability.
3. This continues through the inertial subrange, where energy transfers without significant dissipation.
4. At the Kolmogorov microscale, viscosity converts kinetic energy into heat.

Vortex stretching is the physical mechanism driving this cascade: it transfers energy from large scales to small scales by intensifying vorticity in progressively thinner tubes.