Heat Conduction in a Rod
Heat Conduction in a Rod
Consider the problem of unsteady-state heat conduction in a rod, as governed by the heat equation:
∂T
∂t
2
∂
∂
2
x
The initial and boundary conditions are:
t=0
T(x,0)=2x
0≤x≤1/2
T(x,0)=2(1-x)
1/2≤x≤1
x=0
T(0,t)=0
x=1
T(1,t)=0
where is the temperature, is time, and is the position.
T
t
x
This problem has an analytical solution in the form of a Fourier series after separation of variables:
T(x,t)=sinsin(nπx)
8
2
π
∞
∑
n=1
1
2
n
nπ
2
-t
2
n
2
π
e
This Demonstration plots the solution . The numerical solution obtained using Chebyshev orthogonal collocation is given by the red dots. The analytical solution is given by the blue curve. Excellent agreement between the two solutions is observed. You can vary the value of as well as the number of Chebyshev collocation points, .
T(x,t)
t
N+1