4.1 Irrotational flow

Irrotational flow describes fluid motion where individual fluid particles translate and deform but do not rotate about their own axes. Because the vorticity is everywhere zero, the entire velocity field can be derived from a single scalar function (the velocity potential), which collapses the vector problem into a much simpler scalar one governed by Laplace's equation.

This section covers the mathematical foundations of irrotational flow, the elementary building-block solutions, how to combine them via superposition, and the key theorems (Kutta-Joukowski, Kelvin, Bernoulli) that connect these ideas to lift and pressure.

Definition of irrotational flow

A flow is irrotational when every fluid element has zero net angular velocity about its own center. Formally, the curl of the velocity field vanishes everywhere:

$\nabla \times \vec{V} = 0$

Fluid particles in an irrotational flow can still translate along streamlines and undergo linear deformation (stretching or compression). What they cannot do is spin. This distinction matters: a fluid element moving along a curved streamline is not necessarily rotating in the vorticity sense.

Mathematical representation

Velocity potential function

Because the curl of any gradient is identically zero, the condition $\nabla \times \vec{V} = 0$ is automatically satisfied whenever the velocity field is written as the gradient of a scalar:

$\vec{V} = \nabla \phi$

The scalar $\phi$ is called the velocity potential. Its existence is both necessary and sufficient for irrotationality. Once you know $\phi$, you recover every component of velocity by differentiation, which is far easier than solving a coupled vector system.

Laplace's equation for irrotational flow

Combining the velocity potential with the incompressibility condition $\nabla \cdot \vec{V} = 0$ gives:

$\nabla^2 \phi = 0$

This is Laplace's equation, a linear, second-order PDE. Linearity is the key feature: it means you can add solutions together and still have a valid solution (the superposition principle). Solving Laplace's equation with the appropriate boundary conditions (typically no flow through solid surfaces and a known freestream at infinity) determines $\phi$ and therefore the entire velocity field.

Properties of irrotational flow

Zero vorticity

The vorticity vector is defined as:

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

In irrotational flow, $\vec{\omega} = 0$ at every point. This is the defining property and the one you should check first when asked whether a given velocity field is irrotational.

Path independence of velocity potential

Because $\vec{V} = \nabla \phi$, the line integral of velocity between two points depends only on the endpoints:

$\int_A^B \vec{V} \cdot d\vec{l} = \phi(B) - \phi(A)$

The path you take doesn't matter. This is analogous to how a conservative force has a path-independent work integral, and it guarantees that $\phi$ is a well-defined, single-valued function throughout the domain (with a caveat for multiply-connected regions containing vortices).

Circulation in irrotational flow

Circulation is the line integral of velocity around a closed curve:

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

For any closed curve that can be continuously shrunk to a point without leaving the irrotational region, $\Gamma = 0$. This follows directly from Stokes' theorem and zero vorticity. However, if the curve encloses a singularity (like a point vortex), the circulation can be nonzero. This distinction is central to how potential flow generates lift.

Irrotational vs rotational flow

Most real flows have both irrotational and rotational regions. The irrotational assumption works well away from walls and wakes, where viscous effects are negligible and vorticity hasn't been generated or diffused into the flow.

Potential flow theory

Applicability to irrotational flow

Potential flow theory is the mathematical framework built on three assumptions: the flow is inviscid, incompressible, and irrotational. Under these assumptions, the problem reduces to solving Laplace's equation for $\phi$, from which you can extract velocity fields, pressure distributions (via Bernoulli), and aerodynamic forces.

Limitations of potential flow theory

• It ignores viscosity, so it cannot predict boundary-layer behavior, skin-friction drag, or flow separation.
• It assumes irrotationality everywhere, which breaks down in wakes, recirculation zones, and anywhere vorticity is significant.
• It predicts zero drag on a closed body in steady flow (d'Alembert's paradox), which contradicts reality.
• It cannot predict when or where separation occurs.

Despite these limitations, potential flow gives accurate pressure distributions over the forward portions of streamlined bodies and remains the starting point for panel methods and thin-airfoil theory.

Elementary flows in irrotational flow

Each elementary flow is a solution to Laplace's equation. The formulas below are given in their most common coordinate forms.

Uniform flow

The simplest irrotational flow: constant velocity $U_\infty$ in the $x$-direction.

$\phi = U_\infty x$

Every streamline is a straight horizontal line. Uniform flow serves as the "freestream" onto which other elementary flows are superimposed.

Source/sink flow

A source (positive $Q$) emits fluid radially outward from a point; a sink (negative $Q$) draws fluid radially inward.

• 2-D: $\phi = \frac{Q}{2\pi} \ln r$
• 3-D: $\phi = -\frac{Q}{4\pi r}$

Here $Q$ is the volumetric flow rate (source strength) and $r$ is the distance from the source/sink. The velocity is purely radial and decays as $1/r$ in 2-D or $1/r^2$ in 3-D.

Doublet flow

Place a source and a sink of equal strength an infinitesimal distance apart while increasing their strength so the product remains finite. The result is a doublet with strength $\mu$.

• 2-D: $\phi = -\frac{\mu \cos\theta}{2\pi r}$
• 3-D: $\phi = -\frac{\mu \cos\theta}{4\pi r^2}$

Doublets are the key ingredient for modeling flow around closed bodies (cylinders, spheres).

Vortex flow

A point vortex with circulation $\Gamma$ produces purely tangential velocity that decays as $1/r$:

$\phi = \frac{\Gamma}{2\pi}\theta$

Note that $\phi$ is multi-valued (it increases by $\Gamma$ each time you go around the vortex). The flow is irrotational everywhere except at the vortex center itself, where the velocity is singular. Adding a vortex to a flow around a body introduces circulation and, through the Kutta-Joukowski theorem, lift.

Superposition principle for irrotational flows

Because Laplace's equation is linear, you can add any number of elementary solutions:

$\phi_{total} = \phi_1 + \phi_2 + \phi_3 + \cdots$

The velocity field of the combined flow is simply:

$\vec{V}_{total} = \nabla\phi_{total} = \nabla\phi_1 + \nabla\phi_2 + \cdots$

This is what makes potential flow so powerful. Complex flow patterns around arbitrary shapes can be built up from simple, well-understood pieces.

Irrotational flow around simple geometries

Flow past a cylinder

Superimpose a uniform flow and a 2-D doublet (axis aligned with the freestream). The velocity potential is:

$\phi = U_\infty \left(r + \frac{a^2}{r}\right) \cos\theta$

where $a$ is the cylinder radius. The surface $r = a$ becomes a streamline (no flow through it), so it acts as a solid cylinder. The resulting streamline pattern is symmetric fore-and-aft, which means zero net drag (d'Alembert's paradox). Adding a point vortex of strength $\Gamma$ at the center breaks the vertical symmetry and produces lift.

Flow past a sphere

The 3-D analog uses a uniform flow plus a 3-D doublet:

$\phi = U_\infty \left(r + \frac{a^3}{2r^2}\right) \cos\theta$

where $a$ is the sphere radius. The flow structure is qualitatively similar to the cylinder case, with symmetric streamlines dividing at the front stagnation point and reconnecting at the rear.

Kutta-Joukowski theorem

Lift generation in irrotational flow

The Kutta-Joukowski theorem states that the lift per unit span on a 2-D body in a uniform stream is:

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

where $\rho_\infty$ is the freestream density, $U_\infty$ is the freestream speed, and $\Gamma$ is the circulation around the body. This result is exact within potential flow and applies regardless of body shape.

Circulation and lift relationship

• Positive (counterclockwise) circulation produces upward lift; negative (clockwise) circulation produces downward lift.
• Lift scales linearly with $\Gamma$, $U_\infty$, and $\rho_\infty$.
• For a real airfoil, the Kutta condition (smooth flow departure at the sharp trailing edge) selects the unique value of $\Gamma$ that the physical flow adopts. Without this condition, the potential-flow solution is non-unique.

Kelvin's circulation theorem

Conservation of circulation in irrotational flow

Kelvin's theorem states that the circulation around a material (fluid-following) closed curve remains constant in time for an inviscid, barotropic flow with only conservative body forces:

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

If the flow starts irrotational ($\Gamma = 0$ for every material loop), it stays irrotational. This is why irrotational flow is such a robust assumption for inviscid analysis far from boundaries.

Implications for lift generation

Kelvin's theorem says you can't create net circulation out of nothing in an inviscid fluid. So where does the circulation around a wing come from? When the wing starts moving, a starting vortex of equal and opposite circulation is shed from the trailing edge. The total circulation of the system (wing + starting vortex) remains zero, satisfying Kelvin's theorem, while the bound circulation around the wing provides lift.

Bernoulli's equation in irrotational flow

Pressure-velocity relationship

For steady, incompressible, irrotational flow, Bernoulli's equation holds between any two points in the field (not just along the same streamline):

$\frac{p}{\rho} + \frac{1}{2}V^2 + gz = \text{constant throughout the flow}$

This is stronger than the streamline-restricted form that applies to rotational flows. The constant is the same everywhere because the flow is irrotational.

Applications of Bernoulli's equation

• Pressure distribution on bodies: Once you know the velocity from $\phi$, Bernoulli directly gives the surface pressure. Integrating that pressure yields lift and (in potential flow, zero) pressure drag.
• Stagnation points: Where $V = 0$, the pressure reaches its maximum value, $p_0 = p + \frac{1}{2}\rho V^2$.
• Speed-pressure tradeoff: Faster flow means lower pressure. This explains suction peaks on the upper surface of airfoils and the pressure drop in constricted channels.

Streamlines and equipotential lines

Orthogonality of streamlines and equipotential lines

In 2-D irrotational flow, the stream function $\psi$ and the velocity potential $\phi$ satisfy the Cauchy-Riemann equations. As a result:

• Lines of constant $\phi$ (equipotential lines) and lines of constant $\psi$ (streamlines) intersect at right angles everywhere.
• Together they form an orthogonal flow net.

A common point of confusion: equipotential lines are lines of constant $\phi$, not lines of constant velocity magnitude. The spacing between equipotential lines (and between streamlines) indicates the local speed: closer spacing means higher velocity.

Visualization of irrotational flow patterns

Drawing the orthogonal network of streamlines and equipotential lines gives a complete picture of the flow. Regions where both sets of lines are densely packed correspond to high-speed flow (and, by Bernoulli, low pressure). Regions where they spread apart correspond to low-speed, high-pressure zones. Singularities such as sources, sinks, and vortices appear as points where the regular grid pattern breaks down.

Conformal mapping techniques

Transformation of irrotational flows

Conformal mapping exploits the fact that the 2-D potential flow problem can be written in terms of a complex potential $w(z) = \phi + i\psi$, where $z = x + iy$. An analytic function $\zeta = f(z)$ maps the physical plane to a transformed plane while preserving:

• The orthogonality of streamlines and equipotential lines
• Local angles between any two curves

The strategy is to map a complicated geometry (like an airfoil) onto a simple one (like a circle), solve the easy problem, and then map back.

Examples of conformal mapping applications

• Joukowski transformation: Maps a circle (with an off-center shift) to an airfoil-shaped profile. This is the classic method for obtaining analytic lift predictions on cambered, thick airfoils.
• Kármán-Trefftz transformation: A generalization of the Joukowski map that allows trailing-edge angles other than zero, producing more realistic airfoil shapes.
• Schwarz-Christoffel transformation: Maps the upper half-plane onto the interior of a polygon, useful for analyzing flows around sharp corners, steps, and channel expansions.

These techniques have largely been superseded by numerical panel methods for design work, but they remain valuable for building physical intuition and for validating computational results.