Solutions To Assignment 4 PDF
Solutions To Assignment 4 PDF
Following the process on Page. 85, we get the solution formula is the following
∞
X nπ 2 nπx
u(x, t) = An e−( l
) kt
sin .
n=1
l
where r is a constant, 0 < r < 2πc/l. Write down the series expansion of the solution.
Answer: Let u(x, t) = T (t)X(x), we get
T 00 + rT 0 X 00
= = −λ.
c2 T X
Hence
nπ 2 nπx
λn = ( ) , X(x) = sin , n = 1, 2, ....
l l
Since 0 < r < 2πc/l, we get
p p
Tn (t) = An cos( −∆n t/2) + Bn sin( −∆n t/2), n = 1, 2, ...,
1
MATH 4220 (2012-2013) Partial differential equations CUHK
where ∆n = r2 − ( 2nπc
l
)2 relative to the equation
nπc 2
λ2 + rλ + ( ) = 0.
l
Therefore ∞
X p p nπx
u(x, t) = [An cos( −∆n t/2) + Bn sin( −∆n t/2)] sin .
n=1
l
5. Do the same for 2πc/l < r < 4πc/l.
Answer: Let u(x, t) = T (t)X(x). As above, we get
T 00 + rT 0 X 00
= = −λ.
c2 T X
Hence
nπ 2 nπx
λn = (
) , X(x) = sin , n = 1, 2, ....
l l
For n = 1, since 2πc/l < r < 4πc/l,
+ −
T1 (t) = A1 eλ1 t + B1 eλ1 t ,
where q
−r ± r2 − ( 2πc
l
)2
λ±
1 =
2
be the roots of the equation
πc 2
λ2 + rλ + ( ) = 0.
l
For n ≥ 2, as before,
p p
Tn (t) = An cos( −∆n t/2) + Bn sin( −∆n t/2), n = 2, 3, ...,
where ∆n = r2 − ( 2nπc
l
)2 relative to the equation
nπc 2
λ2 + rλ + ( ) = 0.
l
Therefore
∞
λ+
1 t λ−
1 t
πx X p p nπx
u(x, t) = [A1 e + B1 e ] sin + [An cos( −∆n t/2) + Bn sin( −∆n t/2)] sin .
l n=2
l
6. Separate variables for the equation tut = uxx + 2u with the boundary conditions u(0, t) =
u(π, t) = 0. Show that there are an infinite number of solutions that satisfy the initial condition
u(x, 0) = 0. So uniqueness is false for this equation!
Answer: Let u(x, t) = T (t)X(x), we get
tT 0 − 2T X 00
= = −λ.
T X
2
MATH 4220 (2012-2013) Partial differential equations CUHK
As before we have
λn = n2 , Xn (x) = sin nx, n = 1, 2, 3, ....
By the initial condition, we have
tT 0 − 2T = −λT, T (0) = 0.
Therefore
u(x, t) = ct sin x, for any constant c,
are solutions (You can check them directly). So uniqueness is false for this equation!
Suggested Answer to Exercise 4.2
1. Solve the diffusion problem ut = kuxx in 0 < x < l, with the mixed boundary conditions
u(0, t) = ux (l, t) = 0.
Answer: Let u(x, t) = T (t)X(x), we get
T0 X 00
= = −λ.
kT X
By the boundary condition, we have the following eigenvalue problem
Since the eigenvalues be [(n + 12 )π]2 /l2 and the eigenfunctions be sin[(n + 12 )πx/l] for n =
0, 1, 2, ..., then we have
∞
X
−[
1 )π
(n+ 2
]2 kt (n + 12 )πx
u(x, t) = An e l sin .
n=0
l
2. Consider the equation utt = c2 uxx for 0 < x < l, with the boundary conditions ux (0, t) =
0, u(l, t) = 0 (Neumann at the left, Dirichlet at the right).
(a) Show that the eigenfunctions are cos[(n + 21 )πx/l].
(b) Write the series expansion for a solution u(x, t).
Answer: (a) This can be proved as before (Please see it in the textbook). Here we give
another proof.
Since X 0 (0) = 0, by the even expansion X(−x) = X(x) for −l ≤ x < 0, then X satisfies
Hence
1 1
λn = [(n + )π]2 /l2 , Xn (x) = cos[(n + )πx/l], n = 0, 1, 2, ....
2 2
(b) As before, using (a) we have
∞ h
X (n + 12 )πct (n + 12 )πct i (n + 12 )πx
u(x, t) = An cos + Bn sin cos .
n=0
l l l
3
MATH 4220 (2012-2013) Partial differential equations CUHK
3. Consider diffusion inside an enclosed circular tube. Let its length (circumference) be
2l. Let x denote the arc length parameter where −l ≤ x ≤ l. Then the concentration of the
diffusing substance satisfies
ut = kuxx for − l ≤ x ≤ l
u(−l, t) = u(l, t) and ux (−l, t) = ux (l, t).
Answer: (a) As before, let λ be any complex number and γ be either one of the two square
roots of −λ; the other one is −γ. Then
Assume that a 6= 0.
Answer: Firstly, let’s look for the positive eigenvalues λ = β 2 > 0. As usual, the general
solution of the ODE is
X(x) = C cos βx + D sin βx.
By the boundary conditions,
4
MATH 4220 (2012-2013) Partial differential equations CUHK
Secondly, let’s look for the zero eigenvalue, i.e., X(x) = Ax+B, by the boundary conditions,
al + 1 = 0. Hence λ = 0 is an eigenvalue if and only if al + 1 = 0.
Thirdly, let’s look for the negative eigenvalues λ = −γ 2 < 0. As usual, the solution of the
ODE is
X(x) = C cosh γx + D sinh γx.
By the boundary conditions,
−X 00 = λX
X 0 (0) − a0 X(0) = 0, X 0 (l) + al X(l) = 0.
A − a0 B = 0, and A + al (Al + B) = 0,
which is equivalent to
a0 + al = −a0 al l.
Therefore λ = 0 is an eigenvalue if and only if a0 + al = −a0 al l.
(b) By (a), we have X(x) = B(a0 x + 1), here B is constant.
3. Derive the eigenvalue equation (16) for the negative eigenvalues λ = −γ 2 and the formula
(17) for the eigenfunctions.
Answer: If λ = −γ 2 < 0, we have
Thus
X 0 (x) = Cγ sinh γx + Dγ cosh γx.
Hence by the boundary conditions,
which implies
(a0 + al )γ
tanh γl = − .
γ 2 + a0 al
And the corresponding eigenfunction
a0
X(x) = C cosh γx + C sinh γx,
γ
where C is constant.
5
MATH 4220 (2012-2013) Partial differential equations CUHK
6
MATH 4220 (2012-2013) Partial differential equations CUHK
a + bβ = 0, a cos β + b sin β = 0.
a + b + aβ − bβ = 0, aeβ + be−β = 0.
7
MATH 4220 (2012-2013) Partial differential equations CUHK
1 l −2 2
Z
E= (c ut + u2x ) dx.
2 0
(Compare this definition with Section 2.2.)
(b) Do the same for the Neumann BCs.
(c) For the Robin BCs, show that
1 l −2 2
Z
1 1
ER = (c ut + u2x ) dx + al [u(l, t)]2 + a0 [u(0, t)]2
2 0 2 2
is conserved. Thus, while the total energy ER is still a constant, some of the internal energy is
“lost” to the boundary is a0 and al are positive and “gained” from the boundary if a0 and al
are negative.
Answer: (a) By the wave equation,
Z l
dE 1
= [ 2 ut utt + ux uxt ] dx
dt 0 c
Z l
= [ut uxx + ux uxt ] dx
0
= (ut ux )|l0 = ut (l, t)ux (l, t) − ut (0, t)ux (0, t).
(d) Use part (c) to find half of the eigenvalues explicitly and half of them graphically.
(e) Assuming that all the eigenvalues are nonnegative, make a list of all the eigenfunctions.
(f) Solve the problem ut = kuxx for 0 < x < l, with the BCs given above, and with
u(x, 0) = φ(x).
8
MATH 4220 (2012-2013) Partial differential equations CUHK
(g) Show that, as t → ∞, lim u(x, t) = A + Bx, assuming that you can take limits term by
term.
Answer: (a) Let λ = 0, we get −vxx = 0, which implies v(x) = Ax+B. Since v(x) = Ax+B
satisfy the boundary conditions for any numbers A and B, so λ = 0 is a double eigenvalue.
(b) Let λ = β 2 > 0 and suppose β > 0, −vxx = λv, then v(x) = C cos βx + D sin βx. Hence
by the boundary conditions, we have
C cos βl + D sin βl − C
Dβ = −Cβ sin βl + Dβ cos βl = .
l
Hence λ = β 2 , where β is a root of the following equation
where κ(x) = κ21 for x < a, κ(x) = κ22 for x > a, ρ(x) = ρ21 for x < a, and ρ(x) = ρ22 for x > a.
All these constants are positive and 0 < a < l.
Answer: Let λ = β 2 , then
βρ1 x βρ1 x
X(x) = A cos + B sin , 0 < x < a;
k1 k1
βρ2 x βρ2 x
X(x) = C cos + D sin , a < x < l.
k2 k2
Hence by the boundary conditions,
βρ2 l βρ2 l
A = 0, C cos + D sin = 0;
k2 k2
9
MATH 4220 (2012-2013) Partial differential equations CUHK
Then by using the boundary condition and X(x) and with its first derivative should be contin-
uous at x = a. Such A and B can not exist except A = B = 0.
Let λ = −β 2 , then we have
βρ1 x βρ1 x
X(x) = A cosh + B sinh , 0 < x < a,
k1 k1
βρ2 x βρ2 x
X(x) = C cosh + D sinh , a < x < l.
k2 k2
By using the boundary condition and X(x) and with its first derivative should be continuous
at x = a. We have
βρ2 l βρ2 l
A = 0, C cosh + D sinh ,
k2 k2
βρ1 a βρ2 a βρ2 a
B sinh = C cosh + D sinh ,
k1 k2 k2
βρ1 βρ1 a βρ2 βρ2 a βρ2 βρ2 a
B cosh =C sinh +D cosh .
k1 k1 k2 k2 k2 k2
After computations, we find β is a root of the following equation
ρ1 βρ1 a ρ2 βρ2 (l − a)
coth + coth = 0.
k1 k1 k2 k2
However, since the left handside is always positive unless β = 0. Therefore, it has no negative
eigenvalues.
16. Find the positive eigenvalues and the corresponding eigenfunctions of the fourth-order
operator +d4 /dx4 with the four boundary conditions
Answer: Let λ = β 4 > 0, where β > 0 and X(x) = A cosh(βx) + B sinh(βx) + C cos(βx) +
D sin(βx).
Hence by the boundary conditions,
nπ nπ nπx
β= , λ = ( )4 , Xn (x) = sin( ), n = 1, 2, . . . .
l l l
10
MATH 4220 (2012-2013) Partial differential equations CUHK
17. Solve the fourth-order eigenvalue problem X 0000 = λX in 0 < x < l, with the four
boundary conditions
X(0) = X 0 (0) = X(l) = X 0 (l) = 0,
where λ > 0. (Hint: First solve the fourth-order ODE.)
Answer: Let λ = β 4 > 0 where β > 0, and X(x) = A cosh βx + B sinh βx + C cos(βx) +
D sin(βx). Hence by the boundary conditions,
A + C = 0, B + D = 0,
A cosh βl + B sinh βl + C cos(βl) + D sin(βl) = 0,
and A sinh βl + B cosh βl − C sin(βl) + D cos(βl) = 0.
4
Hence λ = β , where β is a root of the equation
cosh βl cos βl = 1,
and the corresponding eigenfunction is
X(x) = (sinh βl − sin βl)(cosh βx − cos βx) − (cosh βl − cos βl)(sinh βx − sin βx).
Suggested Answer to Extra Problems Problem 10. Solve the following
problem by using separation of variables.
utt + a2 uxxxx = 0, 0 < x < l, t > 0,
u(x, 0) = φ(x), ut (x, 0) = ψ(x),
u(0, t) = uxx (0, t) = u(l, t) = uxx (l, t) = 0.
00 0000
Answer: Let u = X(x)T (t), then −aT 2 T = XX = λ.
Let λ = 0, X(x) = ax3 + bx2 + cx + d, using the boundary condition, we obtain a = b =
c = d = 0. Hence zero can not be an eigenvalue.
Let λ = β 4 , X(x) = A cosh βx + B sinh βx + C cos βx + D sin βx.
Following a similar process as in (17). If we want to obtain a nontrivial solution, we must
have l = kπ
β
. The corresponding eigenfunction is
kπ
X(x) = D sin .
l
Function T (t) suffices the following equation
T (t) = C cos aβ 2 t + D sin aβ 2 t.
Using integration by parts, we can exactly prove this eigenvalue problem does not admit
negative eigenvalue. Therefore the solution to this problem is the following
∞
X nπx
u(x, t) = (Cn cos aβn2 t + Dn sin aβn2 t) sin .
n=1
l
Using the initial condition
∞ ∞
X nπx X nπx
φ(x) = Cn sin , ψ(x) = aβn2 Dn sin .
n=1
l n=1
l
We can decide the Cn , Dn , n = 1, 2, 3, · · · .
11