WOLFRAM|DEMONSTRATIONS PROJECT

Heat Transfer along a Rod

​
conductivity parameter k
0.1
number of terms n
1
This Demonstration shows the solution to the heat equation for a one-dimensional rod. The rod is initially submerged in a bath at 100 degrees and is perfectly insulated except at the ends, which are then held at 0 degrees. This is a a Sturm–Liouville boundary value problem for the one-dimensional heat equation
∂u(x,t)
∂t
-k
2
∂
u(x,t)
∂
2
x
=0
,
with boundary conditions
u(0,t)=0
,
u(L,t)=0
, and
u(x,0)=f(x)
, where
t
is time,
x
is distance along the rod,
L
is the length of the rod, and
f(x)=100
.
The solution is of the form
u(x,t)=
∞
∑
i=1
A
n
sin
nπ
L
-
2
k
nπ
L
t
e
,
where
k
is the conductivity parameter (a product of the density, thermal conductivity, and specific heat of the rod) and
A=
2
L
L
∫
0
f(x)sin
nπ
L
xdx
.
If you increase the number of terms
n
, the solution improves as long as the time is small. As
t∞
(the final state), the entire rod approaches a temperature of 0 degrees. You can see the effect of the thermal properties by varying the conductivity parameter
k
.