Mathematical Methods (Second Year) MT 2009: Problem Set 5: Partial Differential Equations
Mathematical Methods (Second Year) MT 2009: Problem Set 5: Partial Differential Equations
Mathematical Methods (Second Year) MT 2009: Problem Set 5: Partial Differential Equations
2. “Method of Characteristics”
A function u(r, t) satisfies the equation
µ ¶
1 ∂ 2 ∂u 1 ∂2u
r = 2 2
r2 ∂r ∂r c ∂t
where c is a constant. By introducing the new function v(r, t) = ru(r, t), and writing ξ =
r + ct, η = r − ct, reduce this equation to
∂2v
= 0.
∂ξ∂η
Hence show that the general solution u(r, t) has the form u(r, t) = 1r [f (r +ct)+g(r −ct)], where f
and g are arbitrary (twice-differentiable) functions. What does the solution represent physically
if t is the time and r is the 3-D polar variable?
∂2y ∂2y
= .
c2 ∂t2 ∂x2
Find the solution satisfying each of the following initial (t = 0) conditions:-
(i) y(x, 0) = L sin(πx/a), ∂y
∂t (x, 0) = 0;
(ii) y(x, 0) = 0, ∂y
∂t (x, 0) = V sin(2πx/a);
1
5. 2-D Laplace equation in polars
Laplace’s equation in two dimensions may be written, using plane polar coordinates r, θ, as
µ ¶
∂ ∂V (r, θ) ∂ 2 V (r, θ)
r r + = 0.
∂r ∂r ∂θ2
Find all separable solutions of this equation which have the form V (r, θ) = R(r)S(θ), which are
single valued for all r, θ. What property of the equation makes any linear combination of such
solutions also a solution?
A continuous potential V (r, θ) satisfies Laplace’s equation everywhere except on the concentric
circles r = a, r = b where b > a.
(i) Given that V (r = a, θ) = V0 (1 + cos θ), and that V is finite as r → 0, find V in the region
r ≤ a.
(ii) Given, separately, that V (r = b, θ) = 2V0 sin2 θ, and V is finite as r → ∞, find V for r ≥ b.
where ωm,n is the frequency of the mode and L is the length of a side of the square. Express
ωm,n in terms of m and n and give any restrictions which apply to m and n.
q
Show that the second lowest frequency is a factor of 52 larger than the frequency of the lowest
mode, that two modes have this frequency and that by combining them it is possible to have
a node along either diagonal of the square. Deduce the fundamental mode for a triangular
membrane obtained by constraining the membrane along a diagonal.
ω2
∇2 Ez + Ez = 0
c2
for oscillatory solutions. When this is separated in cylindrical polar coordinates the solution has
the form R(ρ)Φ(φ)Z(z) where
µ ¶
d2 Z 2 d2 Φ d dR ω2
2
= −p Z, 2
= −m2 Φ, ρ ρ + [( 2 − p2 )ρ2 − m2 ]R = 0.
dz dφ dρ dρ c
Write ω 2 /c2 − p2 = γ 2 and introduce x = γρ. Show that the R equation is equivalent to
µ ¶
d dJm (x)
x x + (x2 − m2 )Jm (x) = 0.
dx dx
2
ω2 2
αmr n2 π 2
= +
c2 a2 h2
where αmr is the rth root of Jm , i.e. Jm (αmr ) = 0.
d2 P 2 1 d dT m2
= −m P, − sin θ + T = λT.
dφ2 sin θ dθ dθ sin2 θ
(i) Why is m an integer or zero? (ii) Verify that for m = 0, T = cos θ is a solution for a certain
value of λ, and find that value. (iii) Show that the equation for R is
d 2 dR
(r ) = λR.
dr dr
Find the solutions for R corresponding to the particular value of λ in (ii). (iv) What is the
most general solution for V having this θ, φ dependence? (v) An uncharged conducting sphere
of radius a is placed at the origin in an intially uniform electrostatic field (0,0,E). Show that it
behaves as an electric dipole. [Hint: V satisfies the boundary condition V (r = a, θ) = 0 for all
θ on r = a, since it’s a conductor. The other condition required to fix V uniquely is the one at
r → ∞. At very large distances from the sphere, the field must be essentially the same as in the
absence of the field, giving V (r → ∞) = −Ez = −Er cos θ. Use these two boundary conditions
to determine the arbitary constants introduced in (iv)].
9. The temperature T is a one-dimensional bar whose sides are perfectly insulated obeys the heat
flow equation
∂T ∂2T
=κ 2
∂t ∂x
where κ is a positive constant. The bar extends from x = 0 to x = L and is perfectly insulated
at x = L. At t < 0 the temperature is 00 C throughout the bar and at t = 0 the uninsulated
end is placed in contact with a heat bath at 1000 C. Show that the temperature of the bar at
subsequent times is given by
∞ µ ¶ ( µ ¶2 )
T X 4 (2n + 1)πx (2n + 1)π
=1− sin exp −κ t .
100 (2n + 1)π 2L 2L
n=0
Estimate the time taken for the average temperature of the bar to reach 900 C if L = 1m and
κ = 0.001m2 s−1 .
∂ 2 ∂V 1 ∂ ∂V
(r )+ (sin θ ) = 0.
∂r ∂r sin θ ∂θ ∂θ
(i) Show that there is a solution of the form
∞
X αn
V (r, θ) = Pn (cos θ)( + βn rn )
rn+1
n=0
3
where Pn (x) satisfies
µ ¶
d 2 dPn (x)
(1 − x ) + n(n + 1)Pn (x) = 0.
dx dx
(ii) By first calculating the potential along the z-axis, find the potential at arbitrary θ and r > a
for the following charge distributions:
a) A dipole lying along the z-axis and consisting of a charge +1 at z = +a and a charge -1 at
z = −a.
b) A quadrupole lying along the z-axis and consisting of a charge -1 at z = +a, a charge +2 at
z = 0, and a charge -1 at z = −a.
You may assume Pn (1) = 1.