07.04.2013 Views

Fourier Series and Partial Differential Equations Lecture Notes

Fourier Series and Partial Differential Equations Lecture Notes

Fourier Series and Partial Differential Equations Lecture Notes

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Chapter 3. The heat equation 31<br />

3.6 Non-zero steady state<br />

It may be that the temperatures of the ends x = 0 <strong>and</strong> x = L are prescribed <strong>and</strong> constant<br />

but not equal to zero.<br />

Example 3.3 Solve the IBVP<br />

subject to the initial condition<br />

<strong>and</strong> the boundary conditions<br />

∂T<br />

∂t = κ∂2 T<br />

∂x2, 0 < x < L, t > 0, (3.58)<br />

T(x,0) = 0, 0 ≤ x ≤ L, (3.59)<br />

T(0,t) = T0 <strong>and</strong> T(L,t) = T1 for t > 0. (3.60)<br />

We cannot use separation of variables <strong>and</strong> <strong>Fourier</strong> series right at the outset. However, we<br />

conjecture that, as t → ∞, T(x,t) → U(x), where<br />

i.e.<br />

κ d2 U<br />

dx 2 = 0, U(0) = T0 <strong>and</strong> U(L) = T1, (3.61)<br />

U(x) = T0<br />

<br />

1− x<br />

L<br />

<br />

+T1<br />

<br />

x<br />

<br />

. (3.62)<br />

L<br />

If we now put S(x,t) := T(x,t)−U(x), we find that S is a solution of the IBVP<br />

with<br />

∂S<br />

∂t = κ∂2 S<br />

∂x2, 0 < x < L, t > 0, (3.63)<br />

S(0,t) = 0 <strong>and</strong> S(L,t) = 0 for t > 0, (3.64)<br />

<strong>and</strong><br />

<br />

S(x,0) = −T0 1− x<br />

<br />

x<br />

<br />

−T1 . (3.65)<br />

L L<br />

In view of the form of the boundary conditions, this IBVP can be solved by our previous<br />

methods. The solution is<br />

<strong>and</strong> so<br />

T(x,t) = T0<br />

<br />

S(x,t) = 2<br />

π<br />

1− x<br />

L<br />

<br />

+T1<br />

∞<br />

n=1<br />

1<br />

n [−T0 +(−1) n T1]sin<br />

<br />

x<br />

<br />

+<br />

L<br />

2<br />

π<br />

∞<br />

n=1<br />

1<br />

n [−T0 +(−1) n T1]sin<br />

<br />

nπx<br />

<br />

e<br />

L<br />

−n2π2κt/L2 , (3.66)<br />

<br />

nπx<br />

<br />

e<br />

L<br />

−n2π2κt/L2 . (3.67)

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!