← 2026 Paper 2

UPSC 2026 Maths Optional Paper 2 Q6c — Step-by-Step Solution

15 marks · Section B

Two-Dimensional and Axisymmetric Flow · Mechanics & Fluid Dynamics · asked 2× in 14 yrs · Read the full method →

Question

Suppose that there is an infinitely long cylinder of radius 'aa' placed in a uniform stream having velocity −Vi^-V\hat{i}. Then using Milne-Thomson’s circle theorem, what will be its complex velocity potential? If in addition, a circulation round the cylinder of kk is produced, find out the stagnation points. Do they produce a lifting tendency in the vertical direction? Explain how.

Technique

The circle theorem is a machine: give it the potential of the undisturbed stream, and it returns the potential with the cylinder inserted — so the only real work is writing f(z)f(z) for the stream −Vi^-V\hat i correctly (note the minus: this is the mirror image of the textbook +Ui^+U\hat i case, and it flips the sign of the lift at the end). Circulation is then superposed as a line vortex at the centre, which leaves ∣z∣=a|z|=a a streamline because its ψ\psi is a function of rr alone. Stagnation points come from dw/dz=0dw/dz=0, a quadratic in zz whose discriminant changes sign at ∣k∣=4πVa|k|=4\pi Va — the classification by circulation strength is where the marks that most scripts leave on the table are. The lift follows from Blasius (or Kutta–Joukowski), and the physical “how” from Bernoulli applied to the surface speed.

Solution

Throughout: two-dimensional, steady, irrotational flow of an inviscid incompressible fluid of density ρ\rho, in the complex plane z=x+iyz=x+iy; the cylinder’s cross-section is the circle ∣z∣=a|z|=a and the fluid occupies ∣z∣≥a|z|\ge a. The complex potential is w(z)=ϕ+iψw(z)=\phi+i\psi, with the fundamental relation

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

Take V>0V>0, so that the stream runs in the negative xx-direction.

Step 1 — The undisturbed stream.

The stream velocity is q=−Vi^\mathbf q=-V\hat i, i.e. u=−Vu=-V, v=0v=0. Hence dwdz=−V\dfrac{dw}{dz}=-V and, ignoring an irrelevant additive constant,

f(z)=−Vz.(1)f(z)=-Vz . \tag{1}

This has no singularity in the finite plane, so the circle theorem applies with no restriction. (Written out: f=−V(x+iy)f=-V(x+iy) gives ϕ=−Vx\phi=-Vx, ψ=−Vy\psi=-Vy — streamlines y=consty=\text{const} traversed in the sense of decreasing xx. Correct.)

Step 2 — State Milne-Thomson’s circle theorem before using it.

Circle theorem. Let f(z)f(z) be the complex potential of a two-dimensional irrotational flow of incompressible inviscid fluid in which there are no rigid boundaries and all of whose singularities lie in ∣z∣>a|z|>a. Then on introducing the circular cylinder ∣z∣=a|z|=a into the flow, 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.

Two one-line justifications, both worth marks:

(i) 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 therefore

w=f(z)+f(z)‾=2 Re f(z)∈R⟹ψ=Im w=0  on ∣z∣=a,w=f(z)+\overline{f(z)}=2\,\mathrm{Re}\,f(z)\in\mathbb{R}\quad\Longrightarrow\quad \psi=\mathrm{Im}\,w=0\ \text{ on }|z|=a,

i.e. the circle is the streamline ψ=0\psi=0 — exactly the condition for a rigid cylinder.

(ii) The flow far away is unchanged. 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, whose modulus a2/∣z0∣<aa^2/|z_0|<a — inside the cylinder, hence outside the fluid. So no new singularity is introduced into ∣z∣>a|z|>a.

Step 3 — Apply the theorem.

With f(z)=−Vzf(z)=-Vz and VV real, f ⁣(a2zˉ)‾=(−Va2zˉ)‾=−Va2z\overline{f\!\left(\dfrac{a^2}{\bar z}\right)}=\overline{\left(-V\dfrac{a^{2}}{\bar z}\right)}=-V\dfrac{a^{2}}{z}. Therefore

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

Two checks. As ∣z∣→∞|z|\to\infty, dwdz→−V\dfrac{dw}{dz}\to-V, recovering the stream −Vi^-V\hat i. On z=aeiθz=ae^{i\theta}, z+a2z=aeiθ+ae−iθ=2acos⁡θz+\dfrac{a^2}{z}=ae^{i\theta}+ae^{-i\theta}=2a\cos\theta, so w=−2Vacos⁡θw=-2Va\cos\theta is real and ψ=0\psi=0. ✓\checkmark

Step 4 — Superpose the circulation.

The circle theorem produces the acyclic (circulation-free) flow. A circulation kk round the cylinder is obtained by placing a line vortex at the centre z=0z=0 — which is inside the cylinder, hence legitimate.

Convention (state it, because the sign of the final force depends on it): kk is the circulation taken positive in the anticlockwise sense. The corresponding potential is

wv(z)=−ik2πlog⁡z.w_v(z)=-\frac{ik}{2\pi}\log z .

Indeed, with z=reiθz=re^{i\theta}, wv=−ik2πln⁡r+kθ2πw_v=-\dfrac{ik}{2\pi}\ln r+\dfrac{k\theta}{2\pi}, so

ϕv=kθ2π,ψv=−k2πln⁡r,qθ=1r∂ϕv∂θ=k2πr,\phi_v=\frac{k\theta}{2\pi},\qquad \psi_v=-\frac{k}{2\pi}\ln r,\qquad q_\theta=\frac1r\frac{\partial\phi_v}{\partial\theta}=\frac{k}{2\pi r},

and ∮qθ r dθ=k\displaystyle\oint q_\theta\,r\,d\theta=k round any circuit enclosing the origin. Because ψv\psi_v depends on rr alone, it is constant on ∣z∣=a|z|=a, so adding it does not disturb the boundary condition. Hence

w(z)=−V(z+a2z)−ik2πlog⁡z,∣z∣≥a.(3)w(z)=-V\left(z+\frac{a^{2}}{z}\right)-\frac{ik}{2\pi}\log z ,\qquad |z|\ge a. \tag{3}

Check on the circulation. For any circuit CC round the cylinder, ∮Cdwdz dz=Γ+iQ\displaystyle\oint_C\frac{dw}{dz}\,dz=\Gamma+iQ (Γ\Gamma = circulation, QQ = flux). Here dwdz=−V+Va2z2−ik2πz\dfrac{dw}{dz}=-V+\dfrac{Va^{2}}{z^{2}}-\dfrac{ik}{2\pi z} has residue −ik2π-\dfrac{ik}{2\pi} at z=0z=0, so the integral is 2πi(−ik2π)=k2\pi i\left(-\dfrac{ik}{2\pi}\right)=k: circulation kk, zero flux (no source). ✓\checkmark

Step 5 — Stagnation points.

Stagnation points are the zeros of dwdz\dfrac{dw}{dz} lying in the fluid ∣z∣≥a|z|\ge a:

dwdz=−V(1−a2z2)−ik2πz=0.\frac{dw}{dz}=-V\left(1-\frac{a^{2}}{z^{2}}\right)-\frac{ik}{2\pi z}=0 .

Multiplying by −z2-z^{2} (legitimate, since z≠0z\neq0 in the fluid):

Vz2+ik2π z−Va2=0.(4)V z^{2}+\frac{ik}{2\pi}\,z-V a^{2}=0 . \tag{4}

This quadratic has

product of roots=−Va2V=−a2,z=−ik2π±4V2a2−k24π22V=−ik4πV±a2−k216π2V2.(5)\text{product of roots}=\frac{-Va^{2}}{V}=-a^{2},\qquad z=\frac{-\dfrac{ik}{2\pi}\pm\sqrt{4V^{2}a^{2}-\dfrac{k^{2}}{4\pi^{2}}}}{2V} =-\frac{ik}{4\pi V}\pm\sqrt{a^{2}-\frac{k^{2}}{16\pi^{2}V^{2}}} . \tag{5}

Introduce the dimensionless circulation

λ=k4πVa,\lambda=\frac{k}{4\pi V a},

so that (5) becomes the compact

  z±=a(±1−λ2  −  iλ),λ=k4πVa  (6)\boxed{\;z_{\pm}=a\left(\pm\sqrt{1-\lambda^{2}}\;-\;i\lambda\right),\qquad \lambda=\frac{k}{4\pi Va}\;} \tag{6}

Step 6 — Classify by the strength of the circulation. (This is the part answers skip.)

(i) Weak circulation, ∣k∣<4πVa|k|<4\pi Va (i.e. ∣λ∣<1|\lambda|<1). The square root is real and

∣z±∣2=a2[(1−λ2)+λ2]=a2,|z_\pm|^{2}=a^{2}\bigl[(1-\lambda^{2})+\lambda^{2}\bigr]=a^{2},

so both stagnation points lie on the cylinder, at z=aeiθsz=ae^{i\theta_s} with

sin⁡θs=−λ=−k4πVa,cos⁡θs=±1−λ2.(7)\sin\theta_s=-\lambda=-\frac{k}{4\pi Va},\qquad \cos\theta_s=\pm\sqrt{1-\lambda^{2}} . \tag{7}

They are placed symmetrically about the yy-axis and, for k>0k>0, below the horizontal diameter (third and fourth quadrants). For k=0k=0, θs=0,π\theta_s=0,\pi: the points z=±az=\pm a, the front and rear generators of the plain symmetric flow.

(ii) Critical circulation, ∣k∣=4πVa|k|=4\pi Va (∣λ∣=1|\lambda|=1). The two roots coalesce:

z=−ia sgn⁡(k),z=-i a\,\operatorname{sgn}(k),

a single double stagnation point at the lowest point of the cylinder when k>0k>0 (highest when k<0k<0). The critical value is

  kcrit=4πVa  \boxed{\;k_{\text{crit}}=4\pi V a\;}

(iii) Strong circulation, ∣k∣>4πVa|k|>4\pi Va (∣λ∣>1|\lambda|>1). Now a2−k2/16π2V2=iaλ2−1\sqrt{a^2-k^2/16\pi^2V^2}=ia\sqrt{\lambda^{2}-1} is imaginary and

z=ia(−λ±λ2−1),z=ia\left(-\lambda\pm\sqrt{\lambda^{2}-1}\right),

so both roots lie on the yy-axis. Their product is −a2-a^{2}, hence ∣z+∣ ∣z−∣=a2|z_+|\,|z_-|=a^{2}: one is inside the cylinder (it is the inverse point, not part of the fluid and not a physical stagnation point) and one is outside. For k>0k>0 the physical one is

z=−ia(λ+λ2−1),∣z∣=a(λ+λ2−1)>a,z=-ia\left(\lambda+\sqrt{\lambda^{2}-1}\right),\qquad |z|=a\bigl(\lambda+\sqrt{\lambda^{2}-1}\bigr)>a,

a single stagnation point standing off the cylinder, directly below it; the streamline through it closes on itself, trapping a body of fluid that is carried round with the cylinder.

Figure to draw (small, half a text-column, and it earns its space because it makes the classification and the lift argument visible at a glance): draw the circle ∣z∣=a|z|=a with centre OO and the two axes. Along the top and bottom draw three horizontal streamlines with arrowheads pointing in the −x-x direction, labelled "−Vi^-V\hat i"; crowd the lines above the cylinder and space them out below it. Draw a curved arrow anticlockwise round the circle labelled “circulation k>0k>0”. Mark the two stagnation points S1,S2S_1,S_2 on the lower half of the circle, symmetric about the yy-axis, with the angle θs\theta_s from the positive xx-axis to OS2OS_2 and the annotation sin⁡θs=−k/(4πVa)\sin\theta_s=-k/(4\pi Va). Add a thick upward arrow through OO labelled Y=ρVkY=\rho Vk. In a small inset beside it, draw the same circle with a single point below it on the yy-axis, labelled ”∣k∣>4πVa|k|>4\pi Va: stagnation point leaves the cylinder”. Do not add streamlines inside the circle — there is no flow there.

Step 7 — The force on the cylinder: Blasius’ theorem.

Blasius’ theorem. For steady two-dimensional irrotational flow, the force per unit length on a cylinder is given by

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

where CC is any contour enclosing the cylinder and no other singularity.

From (3),

dwdz=−V+Va2z2−ik2πz.\frac{dw}{dz}=-V+\frac{Va^{2}}{z^{2}}-\frac{ik}{2\pi z} .

Only the coefficient of z−1z^{-1} in the square is needed (all other terms integrate to zero). Squaring, the z−1z^{-1} term arises solely from the cross-product of −V-V with −ik2πz-\dfrac{ik}{2\pi z}:

(dwdz)2=⋯+2(−V)(−ik2πz)+⋯=⋯+iVkπ⋅1z+⋯ ,\left(\frac{dw}{dz}\right)^{2}=\cdots+2(-V)\left(-\frac{ik}{2\pi z}\right)+\cdots=\cdots+\frac{iVk}{\pi}\cdot\frac1z+\cdots,

(the terms Va2z2⋅(−V)\frac{Va^2}{z^2}\cdot(-V), Va2z2⋅(−ik2πz)\frac{Va^2}{z^2}\cdot\left(-\frac{ik}{2\pi z}\right) and the two squares contribute only z−2,z−3,z−4z^{-2},z^{-3},z^{-4}). Hence Res⁡z=0(dwdz)2=iVkπ\operatorname{Res}_{z=0}\left(\dfrac{dw}{dz}\right)^{2}=\dfrac{iVk}{\pi} and

X−iY=iρ2⋅2πi⋅iVkπ=iρ2⋅2i2Vk=−iρVk.X-iY=\frac{i\rho}{2}\cdot 2\pi i\cdot\frac{iVk}{\pi}=\frac{i\rho}{2}\cdot 2i^{2}Vk=-i\rho Vk .

Comparing real and imaginary parts:

  X=0,Y=ρVk  (8)\boxed{\;X=0,\qquad Y=\rho V k\;} \tag{8}

X=0X=0 is d’Alembert’s paradox: no drag. The entire force is transverse to the stream — i.e. vertical.

Step 8 — “Explain how”: the physical mechanism, from Bernoulli.

On the cylinder z=aeiθz=ae^{i\theta}, equation (3) gives

ϕ=−2Vacos⁡θ+kθ2π,q=1a∂ϕ∂θ=2Vsin⁡θ+k2πa.(9)\phi=-2Va\cos\theta+\frac{k\theta}{2\pi},\qquad q=\frac1a\frac{\partial\phi}{\partial\theta}=2V\sin\theta+\frac{k}{2\pi a} . \tag{9}

(Setting q=0q=0 in (9) returns sin⁡θ=−k4πVa\sin\theta=-\dfrac{k}{4\pi Va} — an independent confirmation of (7), obtained without touching the quadratic (4).)

Bernoulli’s equation for the steady irrotational flow gives p+12ρq2=p∞+12ρV2p+\tfrac12\rho q^{2}=p_{\infty}+\tfrac12\rho V^{2}, so where the fluid is faster, the pressure is lower. Compare the crown and the keel:

q∣θ=π/2=2V+k2πa,q∣θ=−π/2=−(2V−k2πa),q\big|_{\theta=\pi/2}=2V+\frac{k}{2\pi a},\qquad q\big|_{\theta=-\pi/2}=-\left(2V-\frac{k}{2\pi a}\right), qtop2−qbot2=(2V+k2πa)2−(2V−k2πa)2=4⋅2V⋅k2πa=4Vkπa.q_{\text{top}}^{2}-q_{\text{bot}}^{2}=\left(2V+\tfrac{k}{2\pi a}\right)^{2}-\left(2V-\tfrac{k}{2\pi a}\right)^{2}=4\cdot 2V\cdot\frac{k}{2\pi a}=\frac{4Vk}{\pi a} .

For k>0k>0 this is positive: the flow is faster over the top and slower underneath, hence ptop<pbotp_{\text{top}}<p_{\text{bot}}, hence a net upward force. That is the mechanism.

Integrating it out confirms (8) without Blasius. The pressure force on an element is −p n^ a dθ-p\,\hat n\,a\,d\theta with outward normal n^=(cos⁡θ,sin⁡θ)\hat n=(\cos\theta,\sin\theta), so

Y=−a∫02πpsin⁡θ dθ=ρa2∫02πq2sin⁡θ dθY=-a\int_{0}^{2\pi}p\sin\theta\,d\theta=\frac{\rho a}{2}\int_{0}^{2\pi}q^{2}\sin\theta\,d\theta

(the constant part of pp contributes nothing, since ∫02πsin⁡θ dθ=0\int_0^{2\pi}\sin\theta\,d\theta=0). With (9),

q2sin⁡θ=4V2sin⁡3θ+2Vkπasin⁡2θ+k24π2a2sin⁡θ,q^{2}\sin\theta=4V^{2}\sin^{3}\theta+\frac{2Vk}{\pi a}\sin^{2}\theta+\frac{k^{2}}{4\pi^{2}a^{2}}\sin\theta ,

and ∫02πsin⁡3θ dθ=∫02πsin⁡θ dθ=0\int_0^{2\pi}\sin^{3}\theta\,d\theta=\int_0^{2\pi}\sin\theta\,d\theta=0, ∫02πsin⁡2θ dθ=π\int_0^{2\pi}\sin^{2}\theta\,d\theta=\pi. Hence

Y=ρa2⋅2Vkπa⋅π=ρVk.  ✓Y=\frac{\rho a}{2}\cdot\frac{2Vk}{\pi a}\cdot\pi=\rho Vk . \;\checkmark

The same computation with cos⁡θ\cos\theta in place of sin⁡θ\sin\theta gives X=0X=0, since ∫sin⁡2θcos⁡θ=∫sin⁡θcos⁡θ=∫cos⁡θ=0\int\sin^{2}\theta\cos\theta=\int\sin\theta\cos\theta=\int\cos\theta=0 over a period. ✓\checkmark

Step 9 — Answer the question asked.

Yes — provided k≠0k\neq0, and the tendency is exactly vertical. The reasoning, in the order an examiner wants it:

  1. With k=0k=0 the potential (2) is even in yy: the flow is symmetric about the xx-axis, the stagnation points sit at θ=0,π\theta=0,\pi on the horizontal diameter, the pressure distribution is up–down symmetric, and Y=0Y=0. No lift.
  2. Circulation destroys that symmetry, and the migration of the stagnation points off the horizontal diameter to sin⁡θs=−k/(4πVa)\sin\theta_s=-k/(4\pi Va) is precisely the visible signature of the asymmetry. Their displacement and the lift are two faces of the same broken symmetry — which is the sense in which the stagnation points “produce” the lift.
  3. Quantitatively — the Kutta–Joukowski theorem. In vector form, for circulation Γ=kk^\boldsymbol\Gamma=k\hat k (anticlockwise) in a stream U∞\mathbf U_\infty,
F=ρ U∞×Γ=ρ(−Vi^)×(kk^)=−ρVk (i^×k^)=−ρVk (−j^)=ρVk j^,\mathbf F=\rho\,\mathbf U_{\infty}\times\boldsymbol\Gamma=\rho\bigl(-V\hat i\bigr)\times\bigl(k\hat k\bigr)=-\rho Vk\,(\hat i\times\hat k)=-\rho Vk\,(-\hat j)=\rho V k\,\hat j ,

in agreement with (8). The magnitude is ρV∣k∣\rho V|k| per unit length, the direction is perpendicular to the stream — vertical — upward when k>0k>0 (anticlockwise), downward when k<0k<0. 4. Carry the sign of the stream. For the usual textbook stream +Vi^+V\hat i the same anticlockwise kk gives Y=−ρVkY=-\rho Vk, i.e. downward. Reversing the stream to −Vi^-V\hat i, as this paper does, reverses the lift. Quoting the book formula without this reversal is the single commonest way to lose the last marks here. 5. Validity across the classification. Equation (8) holds for every kk: the lift does not care whether the stagnation points are on the cylinder or standing off below it. What changes at ∣k∣=4πVa|k|=4\pi Va is only their location. ■\qquad\blacksquare

Answer

  Circle theorem: w(z)=−V(z+a2z);with circulation k (anticlockwise): w(z)=−V(z+a2z)−ik2πlog⁡z.Stagnation points: z=−ik4πV±a2−k216π2V2=a(±1−λ2−iλ),λ=k4πVa.∣k∣<4πVa: two points on the cylinder ∣z∣=a,  at sin⁡θs=−k4πVa;∣k∣=4πVa: coincident at z=−iasgn⁡k;∣k∣>4πVa: one point off the cylinder at z=−ia(λ+λ2−1).Force per unit length: X=0  (d’Alembert),Y=ρVk.Yes: a vertical lift of magnitude ρV∣k∣, upward for k>0, downward for k<0.  \boxed{\;\begin{aligned} &\text{Circle theorem: } && w(z)=-V\left(z+\frac{a^{2}}{z}\right);\qquad\text{with circulation } k\ (\text{anticlockwise}):\ w(z)=-V\left(z+\frac{a^{2}}{z}\right)-\frac{ik}{2\pi}\log z.\\[4pt] &\text{Stagnation points: } && z=-\frac{ik}{4\pi V}\pm\sqrt{a^{2}-\frac{k^{2}}{16\pi^{2}V^{2}}}=a\left(\pm\sqrt{1-\lambda^{2}}-i\lambda\right),\qquad \lambda=\frac{k}{4\pi Va}.\\[4pt] & && |k|<4\pi Va:\ \text{two points }\textit{on}\text{ the cylinder } |z|=a,\ \text{ at } \sin\theta_s=-\tfrac{k}{4\pi Va};\\ & && |k|=4\pi Va:\ \text{coincident at } z=-ia\operatorname{sgn}k;\qquad |k|>4\pi Va:\ \text{one point off the cylinder at } z=-ia\bigl(\lambda+\sqrt{\lambda^{2}-1}\bigr).\\[4pt] &\text{Force per unit length: } && X=0\ \ (\text{d'Alembert}),\qquad Y=\rho Vk .\\[2pt] & && \textbf{Yes: a vertical lift of magnitude }\rho V|k|,\ \textbf{upward for }k>0,\ \textbf{downward for }k<0. \end{aligned}\;}
We post more of this — worked solutions, CSAT trap breakdowns, guide chapters — a few times a week on Telegram. Free, no sign-in. Join

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