Take Home Midterm 2013

Problem 1:

A 1D advection-diffusion equation is of the form

∂u ∂u ∂2u
+v =k 2
∂t ∂x ∂x
(a) Solve this equation on the interval x ∈ [−10, 10] with initial and boundary conditions:

u(−10, t) = u(10, t) = 0
2
u(x, 0) = e−x /2
− e−50

(b) Plot the solution for t = 0, 1, 2, 3, 4, 5 for the 3 following cases:

Problem 2:

Your instructor is currently writing a paper (Oglethorpe, Garaud and Wood, in prep) in which she has to
study the following problem. On the interval z ∈ (0, H), we have
2 2
∂4v ∂2v
 
∂v 2Ω(t)R 2Ω(t)
+ κ − κ =0 (1)
∂t Nλ ∂z 4 N ∂z 2
 −1/2
t
where Ω(t) = Ω0 t0 , and where N , λ, R, κ, Ω0 and t0 are just constants. The boundary conditions
on v(z, t) are:

∂v ∂3v
(0, t) = 0 and (0, t) = 0,
∂z ∂z 3
∂v ∂3v
(H, t) = 0 and v(H, t) = K 3 (H, t),
∂z ∂z
where K is a positive constant. The initial conditions are arbitrary, so that

v(z, t0 ) = v0 (z).

1
The paper contains the following statements. For each point marked (**N **) where N = 1..10, prove the
corresponding statement and/or perform the required calculation.

”To solve Equation (1), we note that it is separable in z and t, and write
X
v(z, t) = Vn (t)Zn (z),
n

where the vertical modes Zn (z) satisfy

d4 Zn λ2 d2 Zn  µ 2
n
L(Zn ) ≡ − = Zn (**1**) (2)
dz 4 R2 dz 2 R2
We see that the boundary conditions on the eigenfunctions are homogeneous (**2**). It can easily
be shown that the operator L on the left-hand-side of (2) is symmetric (**3**) with these boundary
conditions, which then implies that the vertical eigenmodes are orthogonal, with
Z H Z H
Zn (z)Zm (z)dz = δnm Zn2 (z)dz. (**4**)
0 0
RH
It can also be shown by considering the integral 0 Zn L(Zn )dz, suitably integrating it by parts, and
applying the boundary conditions, that the eigenvalues associated with the operator L and our boundary
conditions must be strictly positive (**5**), hence our choice of writing them as (µn /R2 )2 in Equation (2).

Since (2) is an equation with constant coefficients, we seek solutions of the form eσn z and find four
solutions for σn : ±σ1,n and ±iσ2,n , where
"r #1/2 "r #1/2
1 λ4 λ2 1 λ4 λ2
σ1,n = µ2n + + , σ2,n = µ2n + − . (**6**)
R 4 2 R 4 2

Applying the boundary conditions we then have

σ2,n sin(σ2,n H)
Zn (z) = cosh (σ1,n (z)) + cos (σ2,n (z)) , (**7**)
σ1,n sinh(σ1,n H)

while the µn coefficients satisfy: (**8** : write the equation that the µn coefficients satisfy).

Projecting (1) onto the Zn functions gives (**9**: prove the equation below and write what is τn (t)
in terms of all the other parameters introduced )

dVn Vn
+ = 0,
dt τn (t)

Solving this equation, using the solution to reconstruct v(z, t), and applying initial conditions finally yields
(**10**: write the complete formal solution). ”

Problem 3:

What are the eigenfrequencies of a rectangular drum of size H × L, assuming that the wave-speed of
the drum is c? Justify all the steps of your calculation very carefully.