UPSC 2026 Maths Optional Paper 2 Q8a — Step-by-Step Solution
15 marks · Section B
Laplace equation: Dirichlet/Neumann, separation of variables · PDEs · asked 5× in 14 yrs · Read the full method →
Question
Two-dimensional harmonic equation in plane polar coordinates (r,θ) takes the form:
∂r2∂2V+r1∂r∂V+r21∂θ2∂2V=0
Deduce that it has solutions of the form (Arn+Br−n)e±inθ, where A, B, n are constants. Determine V if it satisfies the two-dimensional harmonic equation in the region 0≤r≤a, 0≤θ≤2π and satisfies the conditions:
(i)V remains finite as r→0,
(ii)V=n∑cncos(nθ) on r=a.
Technique
Separation of variables V=R(r)Θ(θ) turns Laplace’s equation on the disc into an angular oscillator plus a Cauchy–Euler radial equation r2R′′+rR′−n2R=0, whose trial solution R=rm gives m=±n — that is exactly the printed form. Two things then do all the work: single-valuedness in θ forces n to be an integer, and n=0 is a repeated indicial root, so it carries the separate pair {1,logr} that the printed form does not contain. Finiteness at the origin kills r−n and logr; matching Fourier coefficients on r=a fixes the rest.
Solution
Step 1 — Separate the variables
Look for a non-trivial product solution
V(r,θ)=R(r)Θ(θ),R≡0,Θ≡0.
Substituting into Vrr+r1Vr+r21Vθθ=0,
R′′Θ+r1R′Θ+r21RΘ′′=0.
Multiply through by RΘr2 (legitimate wherever RΘ=0) and separate:
Rr2R′′+rR′=−ΘΘ′′.
The left side depends on r alone and the right on θ alone, so each equals a common constant. Call it λ:
Θ′′+λΘ=0,(1)r2R′′+rR′−λR=0.(2)
Step 2 — The radial equation is Cauchy–Euler; deduce r±n
Equation (2) is homogeneous of Cauchy–Euler (equidimensional) type: substituting the trial solution R=rm gives
m(m−1)rm+mrm−λrm=0⟹m2=λ(indicial equation).
Write n=λ, so that λ=n2 and m=±n.
If n=0 the indicial roots m=n and m=−n are distinct, and the two independent radial solutions are rn and r−n:
R(r)=Arn+Br−n.
With λ=n2, equation (1) reads Θ′′+n2Θ=0, whose independent solutions are
Θ(θ)=einθande−inθ.
Multiplying,
V=(Arn+Br−n)e±inθ
which is the form asserted in the question. ■
Replacing n by −n merely interchanges A↔B and the two exponentials, so no generality is lost by taking n≥0.
Step 3 — The exceptional case n=0 (which the printed form omits)
If n=0 the indicial equation m2=0 has a repeated root, and rn,r−n both collapse to the single solution 1. The second solution is obtained from (2) directly: with λ=0,
r2R′′+rR′=0⟺(rR′)′=0⟺rR′=const⟺R=A0+B0logr.
Likewise (1) becomes Θ′′=0, giving Θ=C0+D0θ. So the n=0 family is
V=(A0+B0logr)(C0+D0θ),
not covered by (Arn+Br−n)e±inθ. This matters here, because the boundary data in (ii) contains the term n=0.
Step 4 — Single-valuedness forces n∈Z and kills the θ-linear term
The point (r,θ) and the point (r,θ+2π) are the same physical point of the disc. A potential must therefore be single-valued:
V(r,θ+2π)=V(r,θ)for all r,θ.
For the n=0 family this requires ein(θ+2π)=einθ, i.e.
e2πin=1⟹2πin∈2πiZ⟹n∈Z.
(Note this rules out complex λ automatically: a non-real n would give a factor e−(Imn)θ, which is not 2π-periodic. In particular λ≥0.)
For the n=0 family, single-valuedness forces D0=0, removing the θ term.
Hence, taking n=0,1,2,… and forming real combinations cosnθ,sinnθ of e±inθ, the general single-valued solution on the punctured disc 0<r≤a is
while rn→0 and the constant stays bounded. Since the terms of (3) are mutually orthogonal in θ, each Fourier mode must be separately bounded, and boundedness of V near the origin therefore forces
B0=0,Bn=Bn′=0(n≥1).
So
V(r,θ)=A0+n=1∑∞rn(Ancosnθ+An′sinnθ).(4)
Condition (i) is not decorative — it is what makes the problem well-posed. Without it, log(r/a) and (anr−n−a−nrn)cosnθ are harmonic on 0<r≤a and vanish identically on r=a; they could be added to any solution without disturbing (ii). Finiteness at the origin is exactly the condition that removes them.
Step 6 — Impose condition (ii): match Fourier coefficients on r=a
Take the index set in (ii) to be n=0,1,2,…, the n=0 term being the constant c0 (since cos0=1); assume ∑n∣cn∣<∞ so the boundary series converges uniformly. Setting r=a in (4),
Both sides are uniformly convergent Fourier series on [0,2π], so their coefficients agree. Extracting them by orthogonality — multiply by 1, by cosmθ, by sinmθ and integrate over [0,2π], using
It really is harmonic. For 0≤r≤ρ<a we have cn(r/a)n≤∣cn∣(ρ/a)n, and after k term-by-term differentiations each term is bounded by ∣cn∣nkρ−k(ρ/a)n; since ρ/a<1 and ∑∣cn∣<∞, these majorants are summable. By the Weierstrass M-test the differentiated series converge uniformly on every closed sub-disc, so term-by-term differentiation is valid. Each term rncosnθ satisfies the equation by Step 2, hence so does V.
It is the only such solution.V is harmonic on the open disc and continuous on the closed disc (the boundary series converges uniformly), and by the maximum principle for harmonic functions the difference of two such solutions is harmonic with zero boundary values, hence identically zero. So the answer is unique.
Sanity value at the centre.V(0,θ)=c0, which is precisely the mean of the boundary data 2π1∫02π∑ncncosnθdθ=c0 — the mean-value property of harmonic functions. ■
Remark (closed integral form). Since a2−2arcos(θ−ϕ)+r2a2−r2=1+2∑n≥1(ar)ncosn(θ−ϕ), the same answer can be written as the Poisson integral