Two-Dimensional and Axisymmetric Flow

At a Glance

Why This Chapter Matters

The connection between 2D irrotational flow and complex analysis is one of the most elegant results in applied mathematics: the velocity potential and stream function satisfy the Cauchy-Riemann equations, making the complex potential w=ϕ+iψw = \phi + i\psi an analytic function of z=x+iyz = x + iy. UPSC 2015 tested this with a derivation of the complex potential for a standard flow — a Section B question requiring both theoretical justification and explicit computation.

Minimum Theory

2D Irrotational Incompressible Flow

Consider a 2D flow in the xyxy-plane with velocity (u,v)(u, v):

Velocity Potential and Stream Function

Velocity potential ϕ\phi (exists for irrotational flow):

u=∂ϕ∂x,v=∂ϕ∂y.u = \frac{\partial \phi}{\partial x}, \qquad v = \frac{\partial \phi}{\partial y}.

Irrotationality is automatic: ∂v/∂x−∂u/∂y=ϕxy−ϕyx=0\partial v/\partial x - \partial u/\partial y = \phi_{xy} - \phi_{yx} = 0. Incompressibility gives ∇2ϕ=0\nabla^2\phi = 0 (Laplace’s equation for ϕ\phi).

Stream function ψ\psi (exists for incompressible flow):

u=∂ψ∂y,v=−∂ψ∂x.u = \frac{\partial \psi}{\partial y}, \qquad v = -\frac{\partial \psi}{\partial x}.

Continuity is automatic. Irrotationality gives ∇2ψ=0\nabla^2\psi = 0 (Laplace’s equation for ψ\psi).

The Complex Potential — Key Result

Define the complex potential:

w=ϕ(x,y)+i ψ(x,y).w = \phi(x,y) + i\,\psi(x,y).

Claim. ww is an analytic function of z=x+iyz = x + iy.

Proof. Check the Cauchy-Riemann equations for w=ϕ+iψw = \phi + i\psi:

∂ϕ∂x=u=∂ψ∂y✓\frac{\partial \phi}{\partial x} = u = \frac{\partial \psi}{\partial y} \quad \checkmark

∂ϕ∂y=v=−(−∂ψ∂x)  ⟹  ∂ϕ∂y=−(−v)\frac{\partial \phi}{\partial y} = v = -\left(-\frac{\partial \psi}{\partial x}\right) \implies \frac{\partial \phi}{\partial y} = -\left(-v\right)

More carefully: ϕx=u=ψy\phi_x = u = \psi_y and ϕy=v=−ψx\phi_y = v = -\psi_x. Thus:

ϕx=ψyandϕy=−ψx.\phi_x = \psi_y \quad \text{and} \quad \phi_y = -\psi_x.

These are exactly the Cauchy-Riemann equations for w=ϕ+iψw = \phi + i\psi. Hence ww is analytic. ■\blacksquare

Complex Velocity

The derivative of the complex potential gives the complex velocity:

dwdz=∂ϕ∂x+i∂ψ∂x=u+i(−v)=u−iv.\frac{dw}{dz} = \frac{\partial \phi}{\partial x} + i\frac{\partial \psi}{\partial x} = u + i(-v) = u - iv.

Thus if you know w(z)w(z), the velocity components are recovered by:

u−iv=dwdz.u - iv = \frac{dw}{dz}.

The speed is ∣dw/dz∣=u2+v2|dw/dz| = \sqrt{u^2 + v^2}.

Standard Complex Potentials

Flow typeComplex potential w(z)w(z)Complex velocity dw/dzdw/dz
Uniform flow at speed UU (along xx-axis)UzUzUU
Uniform flow at angle α\alphaUe−iαzUe^{-i\alpha}zUe−iαUe^{-i\alpha}
Line source/sink of strength mm at originm2πlog⁡z\dfrac{m}{2\pi}\log zm2πz\dfrac{m}{2\pi z}
Line vortex of strength κ\kappa at origin−iκ2πlog⁡z\dfrac{-i\kappa}{2\pi}\log z−iκ2πz\dfrac{-i\kappa}{2\pi z}
Doublet of strength μ\mu at originμ2πz\dfrac{\mu}{2\pi z}−μ2πz2\dfrac{-\mu}{2\pi z^2}
Flow past cylinder of radius aa, speed UUU ⁣(z+a2z)U\!\left(z + \dfrac{a^2}{z}\right)U ⁣(1−a2z2)U\!\left(1 - \dfrac{a^2}{z^2}\right)

Warning about the last row. It is a consequence of the circle theorem below, not an independent fact. If a question says “using Milne-Thomson’s circle theorem, find the complex potential”, quoting this row answers a different question and forfeits the marks attached to the theorem. Derive it.

Milne-Thomson’s Circle Theorem

Statement. Let f(z)f(z) be the complex potential of a 2D irrotational flow of incompressible inviscid fluid with no rigid boundaries, all of whose singularities lie in ∣z∣>a|z|>a. On introducing the circular cylinder ∣z∣=a|z|=a, the complex potential becomes

w(z)=f(z)+f ⁣(a2zˉ)‾,w(z)=f(z)+\overline{f\!\left(\frac{a^2}{\bar z}\right)},

and ∣z∣=a|z|=a is a streamline of the new flow.

Why it works — two one-line justifications, both examinable.

  1. The boundary condition holds. On ∣z∣=a|z|=a we have zzˉ=a2z\bar z=a^2, so a2/zˉ=za^2/\bar z=z and w=f(z)+f(z)‾=2 Re f(z)w=f(z)+\overline{f(z)}=2\,\mathrm{Re}\,f(z), which is real. Hence ψ=Im w=0\psi=\mathrm{Im}\,w=0 on the circle: it is a streamline, which is exactly the condition for a rigid cylinder.
  2. The far field is undisturbed. If ff is singular at z0z_0 with ∣z0∣>a|z_0|>a, the added term is singular only at the inverse point a2/zˉ0a^2/\bar z_0, of modulus a2/∣z0∣<aa^2/|z_0|<a — inside the cylinder, hence outside the fluid. No new singularity enters the flow region.

Adding circulation. The circle theorem produces the acyclic (circulation-free) flow. A circulation κ\kappa round the cylinder is added as a line vortex at the centre — legitimate because the centre is inside the cylinder, and harmless to the boundary condition because the vortex’s ψ=−κ2πln⁡r\psi=-\tfrac{\kappa}{2\pi}\ln r depends on rr alone and so is constant on ∣z∣=a|z|=a:

w(z)=f(z)+f ⁣(a2zˉ)‾−iκ2πlog⁡z,κ>0 anticlockwise.w(z)=f(z)+\overline{f\!\left(\frac{a^2}{\bar z}\right)}-\frac{i\kappa}{2\pi}\log z,\qquad \kappa>0 \text{ anticlockwise}.

Blasius’ theorem (force on the cylinder). For steady 2D irrotational flow, the force per unit length is

X−iY=iρ2∮C(dwdz)2dz,X-iY=\frac{i\rho}{2}\oint_C\left(\frac{dw}{dz}\right)^{2}dz,

CC any contour enclosing the cylinder and no other singularity. Only the residue at the origin survives. Equivalently, Kutta–Joukowski: F=ρ U∞×Γ\mathbf{F}=\rho\,\mathbf{U}_\infty\times\boldsymbol\Gamma, of magnitude ρU∣κ∣\rho U|\kappa| perpendicular to the stream — the drag is always zero (d’Alembert’s paradox).

Axisymmetric Flow (Stokes Stream Function)

For flow with symmetry about the zz-axis (cylindrical coordinates r,θ,zr, \theta, z with no θ\theta-dependence), a Stokes stream function Ψ(r,z)\Psi(r,z) is defined so that:

Unlike the 2D case, Ψ\Psi does not generally satisfy Laplace’s equation; instead it satisfies the Stokes stream function equation.

Question Archetypes

ArchetypeRecognition
find-complex-potential”Find the complex potential for the given 2D irrotational flow”
circle-theorem-cylinder”Using Milne-Thomson’s circle theorem…”; a cylinder is inserted into a stream, usually with circulation and a lift question
derive-cr-structure”Show that ϕ\phi and ψ\psi satisfy the C-R equations / ww is analytic”
combine-flows”Superpose source, vortex, uniform stream; find the combined ww“
axisymmetric-properties”State properties of the Stokes stream function; find flux”

find-complex-potential (1 question; 2015)

Recognition Cues

Solution Template

  1. Identify the flow components (uniform stream, source, vortex, doublet, etc.) from the problem description.
  2. Write the complex potential for each component using the standard table.
  3. Superpose: w=w1+w2+…w = w_1 + w_2 + \ldots (complex potentials add linearly).
  4. Extract ϕ=Re(w)\phi = \text{Re}(w) and ψ=Im(w)\psi = \text{Im}(w) if asked.
  5. Compute dw/dz=u−ivdw/dz = u - iv to find velocity components.
  6. If asked to prove ww is analytic: verify C-R equations ϕx=ψy\phi_x = \psi_y, ϕy=−ψx\phi_y = -\psi_x using the flow relations u=ϕx=ψyu = \phi_x = \psi_y, v=ϕy=−ψxv = \phi_y = -\psi_x.

Worked Example

2015 Paper 2, 2015-P2-QMF (20 marks)

Find the complex potential for the 2D irrotational flow consisting of a uniform stream of speed UU parallel to the xx-axis superposed with a line source of strength mm at the origin. Find the velocity at the stagnation point and the equation of the streamline through the stagnation point.

Step 1: Write the complex potential.

Uniform stream: w1=Uzw_1 = Uz.

Line source of strength mm at origin: w2=m2πlog⁡zw_2 = \dfrac{m}{2\pi}\log z.

Superposed complex potential:

w(z)=Uz+m2πlog⁡z.w(z) = Uz + \frac{m}{2\pi}\log z.

Step 2: Find the complex velocity.

dwdz=U+m2πz=u−iv.\frac{dw}{dz} = U + \frac{m}{2\pi z} = u - iv.

Step 3: Find the stagnation point.

At a stagnation point, u=v=0u = v = 0, i.e., dw/dz=0dw/dz = 0:

U+m2πz=0  ⟹  z=−m2πU.U + \frac{m}{2\pi z} = 0 \implies z = -\frac{m}{2\pi U}.

This is the point (−m2πU, 0)\left(-\dfrac{m}{2\pi U},\, 0\right) on the negative xx-axis. The velocity there is zero (stagnation).

Step 4: Stream function and stagnation streamline.

Write z=reiθz = re^{i\theta} (polar form). Then log⁡z=ln⁡r+iθ\log z = \ln r + i\theta, so:

w=U(x+iy)+m2π(ln⁡r+iθ),w = U(x + iy) + \frac{m}{2\pi}(\ln r + i\theta),

ψ=Im(w)=Uy+mθ2π=Ursin⁡θ+mθ2π.\psi = \text{Im}(w) = Uy + \frac{m\theta}{2\pi} = Ur\sin\theta + \frac{m\theta}{2\pi}.

At the stagnation point z=−m/(2πU)z = -m/(2\pi U): here r=m/(2πU)r = m/(2\pi U), θ=π\theta = \pi, so:

ψstag=U⋅m2πU⋅sin⁡π+m⋅π2π=0+m2=m2.\psi_{\text{stag}} = U \cdot \frac{m}{2\pi U} \cdot \sin\pi + \frac{m\cdot\pi}{2\pi} = 0 + \frac{m}{2} = \frac{m}{2}.

The stagnation streamline is ψ=m/2\psi = m/2:

Ursin⁡θ+mθ2π=m2.\boxed{Ur\sin\theta + \frac{m\theta}{2\pi} = \frac{m}{2}.}

This streamline divides the flow: above and below it the stream passes around the “body” formed by the stagnation dividing streamline.

Step 5: Verify ww is analytic.

w(z)=Uz+m2πlog⁡zw(z) = Uz + \dfrac{m}{2\pi}\log z is a sum of analytic functions (linear function, which is entire, plus log⁡z\log z, which is analytic on C∖(−∞,0]\mathbb{C}\setminus(-\infty,0]). Hence ww is analytic on C∖{0}∩C∖(−∞,0]\mathbb{C}\setminus\{0\} \cap \mathbb{C}\setminus(-\infty,0], i.e., on C\mathbb{C} cut along the negative real axis. The corresponding ϕ\phi and ψ\psi automatically satisfy the C-R equations:

ϕx=ψy=u,ϕy=−ψx=v.■\phi_x = \psi_y = u, \quad \phi_y = -\psi_x = v. \quad \blacksquare

w(z)=Uz+m2πlog⁡z\boxed{w(z) = Uz + \frac{m}{2\pi}\log z}

Common Traps


circle-theorem-cylinder (1 question(s); 2026)

The separating feature from find-complex-potential is that a rigid boundary is introduced into an existing flow. Superposition from the standard table cannot do that; the circle theorem is a machine that does it for you. These questions almost always continue into circulation, stagnation points and lift, and that continuation is where most of the marks sit.

Recognition Cues

Solution Template

  1. Write f(z)f(z) for the undisturbed stream and check it has no finite singularity, so the theorem applies unrestricted. For q=−Vi^\mathbf q=-V\hat i: dw/dz=u−iv=−Vdw/dz=u-iv=-V, so f(z)=−Vzf(z)=-Vz.
  2. State the circle theorem before using it, with the two justifications (boundary condition; inverse-point singularity). Then apply it.
  3. Superpose the circulation as a central vortex, and declare the sign convention for κ\kappa — the sign of the final force depends on it and the paper will not tell you.
  4. Stagnation points: solve dw/dz=0dw/dz=0. Multiplying up gives a quadratic in zz, not a linear equation. Introduce the dimensionless λ=κ/(4πVa)\lambda=\kappa/(4\pi Va) to make the roots readable.
  5. Classify by circulation strength — this is the part scripts skip. The discriminant changes sign at ∣κ∣=4πVa|\kappa|=4\pi Va: two points on the cylinder below it, a coincident double point at it, and beyond it two points on the vertical axis of which one lies inside the cylinder and is not a stagnation point of the flow at all (the root product is −a2-a^2, so ∣z+∣∣z−∣=a2|z_+||z_-|=a^2). Discard it.
  6. Lift by Blasius (or Kutta–Joukowski), then explain the mechanism from Bernoulli: compare the surface speed at crown and keel.

Worked Example

2026 Paper 2, 2026-P2-Q6c (15 marks)

An infinitely long cylinder of radius aa is placed in a uniform stream of velocity −Vi^-V\hat i. Using Milne-Thomson’s circle theorem, find its complex velocity potential. If a circulation κ\kappa round the cylinder is produced, find the stagnation points. Do they produce a lifting tendency in the vertical direction? Explain how.

Source: analysis/solutions/2026-P2-Q6c.md

Step 1 — Undisturbed stream. u=−Vu=-V, v=0v=0, so dw/dz=−Vdw/dz=-V and f(z)=−Vzf(z)=-Vz (constant discarded). No finite singularity, so the theorem applies with no restriction.

Step 2 — Circle theorem. Since VV is real, f(a2/zˉ)‾=−Va2/zˉ‾=−Va2/z\overline{f(a^2/\bar z)}=\overline{-V a^2/\bar z}=-Va^2/z, giving

w(z)=−V ⁣(z+a2z),∣z∣≥a.\boxed{w(z)=-V\!\left(z+\frac{a^2}{z}\right)},\qquad |z|\ge a.

Checks: dw/dz→−Vdw/dz\to-V as ∣z∣→∞|z|\to\infty ✓; on z=aeiθz=ae^{i\theta}, z+a2/z=2acos⁡θz+a^2/z=2a\cos\theta is real so ψ=0\psi=0 ✓.

Step 3 — Circulation. Taking κ>0\kappa>0 anticlockwise,

w(z)=−V ⁣(z+a2z)−iκ2πlog⁡z.w(z)=-V\!\left(z+\frac{a^2}{z}\right)-\frac{i\kappa}{2\pi}\log z.

Step 4 — Stagnation points. dwdz=−V ⁣(1−a2z2)−iκ2πz=0\dfrac{dw}{dz}=-V\!\left(1-\dfrac{a^2}{z^2}\right)-\dfrac{i\kappa}{2\pi z}=0. Multiplying by −z2-z^2:

Vz2+iκ2πz−Va2=0  ⟹  z±=a ⁣(±1−λ2−iλ),λ=κ4πVa.Vz^2+\frac{i\kappa}{2\pi}z-Va^2=0\;\Longrightarrow\;z_\pm=a\!\left(\pm\sqrt{1-\lambda^2}-i\lambda\right),\qquad \lambda=\frac{\kappa}{4\pi Va}.

Step 5 — The three cases.

Step 6 — Lift, by Blasius. With dwdz=−V+Va2z2−iκ2πz\dfrac{dw}{dz}=-V+\dfrac{Va^2}{z^2}-\dfrac{i\kappa}{2\pi z}, the only z−1z^{-1} term in the square comes from the cross-product 2(−V) ⁣(−iκ2πz)=iVκπz2(-V)\!\left(-\tfrac{i\kappa}{2\pi z}\right)=\tfrac{iV\kappa}{\pi z}, so Res⁡0(dw/dz)2=iVκπ\operatorname{Res}_0(dw/dz)^2=\tfrac{iV\kappa}{\pi} and

X−iY=iρ2⋅2πi⋅iVκπ=−iρVκ⟹X=0,Y=ρVκ.X-iY=\frac{i\rho}{2}\cdot2\pi i\cdot\frac{iV\kappa}{\pi}=-i\rho V\kappa\quad\Longrightarrow\quad\boxed{X=0,\qquad Y=\rho V\kappa.}

Step 7 — “Explain how”. On the cylinder, ϕ=−2Vacos⁡θ+κθ2π\phi=-2Va\cos\theta+\tfrac{\kappa\theta}{2\pi}, so the surface speed is q=1a∂θϕ=2Vsin⁡θ+κ2πaq=\tfrac1a\partial_\theta\phi=2V\sin\theta+\tfrac{\kappa}{2\pi a}. (Setting q=0q=0 returns sin⁡θs=−λ\sin\theta_s=-\lambda — an independent confirmation of Step 5.) Then

qtop2−qbot2=(2V+κ2πa)2−(2V−κ2πa)2=4Vκπa>0(κ>0),q^2_{\text{top}}-q^2_{\text{bot}}=\left(2V+\tfrac{\kappa}{2\pi a}\right)^2-\left(2V-\tfrac{\kappa}{2\pi a}\right)^2=\frac{4V\kappa}{\pi a}>0\quad(\kappa>0),

so by Bernoulli (p+12ρq2p+\tfrac12\rho q^2 constant) the pressure is lower over the top and there is a net upward force. With κ=0\kappa=0 the potential is even in yy, the stagnation points sit at θ=0,π\theta=0,\pi and the lift vanishes: circulation is what breaks the up–down symmetry, and the migration of the stagnation points off the horizontal diameter is the visible signature of the same broken symmetry that produces the lift. Cross-check by Kutta–Joukowski: ρ(−Vi^)×(κk^)=+ρVκ j^\rho(-V\hat i)\times(\kappa\hat k)=+\rho V\kappa\,\hat j ✓.

Common Traps

Marks-Aware Writing

This is a 20-mark Section B derivation question. UPSC expects:

  1. Explicit statement of each component’s complex potential — do not just write the answer.
  2. Superposition step clearly shown.
  3. dw/dzdw/dz computed and interpreted as u−ivu - iv.
  4. Stagnation point found by setting dw/dz=0dw/dz = 0, solving for zz.
  5. Stream function extracted as ψ=Im(w)\psi = \text{Im}(w) — shown in polar coordinates.
  6. Stagnation streamline: evaluate ψ\psi at the stagnation point, then write the streamline equation.
  7. Analyticity of ww: mention C-R, even if briefly.

Marks are distributed across all these steps. A bare final answer without working receives minimal credit.

Practice Set

Only one historical question on this atom (shown above).

Ready to drill what you just read?

Daily Practice turns these patterns into one adaptive set a day — practised daily until they're automatic, free for everyone.

See Daily Practice →

This chapter is part of the Maths Coverage Map — 14 years, mapped. Get the take-away PDF free.