Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Mathematical Methods (Second Year) MT 2009: Problem Set 5: Partial Differential Equations

Download as pdf or txt
Download as pdf or txt
You are on page 1of 4

Mathematical Methods (Second Year) MT 2009

Problem Set 5: Partial Differential Equations

1. ‘Linear’ and ‘Homogeneous’


What is meant by the terms (a) linear and (b) homogeneous as applied to differential equations?
Give a physical example of a second order linear non-homogeneous ordinary differential equation,
and of a second order linear homogeneous partial differential equation. In both examples, discuss
whether the linearity is believed to be exact, or is the result of an approximation (and if so, say
what it is). In the case of a second order partial differential equation, describe the properties of
the solutions that follow when it is (a) linear (b) linear and homogeneous.

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?

3. Modes of a 1-D string


The transverse displacement y(x, t) (assumed small) of a string stretched between the points
x = 0 and x = a satisfies the equation

∂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);

(iii) y(x, 0) = L sin(πx/a), ∂y


∂t (x, 0) = V sin(2πx/a);
πx ∂y
(iv) y(x, 0) = 2L sin( 3πx
2a ) cos( 2a ), ∂t (x, 0) = 0.
Sketch the displacement of the string in part (iv) at time t = a/2c.

4. Damped thermal wave


The temperature T in a one-dimensional bar whose sides are perfectly insulated obeys the heat
flow equation
∂T (x, t) ∂ 2 T (x, t)
=K
∂t ∂x2
where K is a constant. The bar extends from x = 0 to x = ∞. The temperature at the end
x = 0 oscillates in time according to T (x = 0, t) = T0 cos ωt. By looking for solutions that are
separated in x and t (T = X(x)F (t)) find the solution for all x ≥ 0, and t, which matches the
ω
boundary condition at x = 0 . Sketch T versus x for ωt = π/2, given that 2K = a12 .

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.

6. Modes of a square membrane


The equation governing the motion of a square vibrating membrane has the form
µ ¶
∂2u ∂2u ∂2u
= c2 + 2
∂t2 ∂x2 ∂y
with u = 0 on the boundary.
Show that the normal modes are given by

u(x, y, t) = Am,n sin(mπx/L) sin(nπy/L) cos(ωm,n t)

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.

7. Cylindrical symmetry → Bessel functions


A resonant electromagnetic cavity is in the shape of a right circular cylinder of radius a and
height h. Taking the z-axis along the axis of the cylinder, the electric field in this direction
satisfies

ω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

The boundary conditions on z = 0, h require Z ∼ cos(nπz/h) and that R(ρ = a) = 0. Show


that the frequencies of allowed modes are given by

2
ω2 2
αmr n2 π 2
= +
c2 a2 h2
where αmr is the rth root of Jm , i.e. Jm (αmr ) = 0.

8. Spherical symmetry → Spherical Harmonics


When Laplace’s equation ∇2 V = 0 is separated in spherical polar coordinates, the solutions
have the form V = R(r)T (θ)P (φ) where

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 .

10. Use of other eigenfunction equations


Laplace’s equation in spherical polar coordinates (r, θ, φ) for a potential V (r, θ) which is inde-
pendent of φ is

∂ 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.

You might also like