WOLFRAM|DEMONSTRATIONS PROJECT

Transient Heat Conduction with Temperature-Dependent Thermal Conductivity

​
number of collocation points
10
δ parameter
0.5
Consider transient heat conduction in a finite slab with temperature-dependent thermal conductivity. The phenomenon is governed by the following dimensionless equation:
∂θ
∂τ
=f(θ,δ)
2
∂
θ
∂
2
x
+f'(θ,δ)
2
∂θ
∂x
,
where
f(θ,δ)
is the dimensionless thermal diffusivity and
δ
is an empirical parameter to be chosen by the user. The variables
τ
,
θ
, and
x
are the dimensionless time, temperature, and position. The initial condition and boundary conditions are:
θ(x,τ=0)=0
,
θ(x=1,τ)=1
, and
∂θ
∂x
x=0
=0
.
This Demonstration plots the dimensionless temperature versus the dimensionless position in the slab for
τ=0.01
,
0.05
,
0.2
, and
0.5
(the orange, green, red, and cyan curves, respectively). The colored dots correspond to the solution using orthogonal collocation. You can specify the number of collocation points for the calculation as well as the functional dependency for the thermal diffusivity through the parameter
δ
. The colored curves correspond to the dimensionless temperature obtained using the Mathematica built-in function NDSolve. The solid curves correspond to
f(θ,δ)=exp(δθ)
; the dashed curves correspond to
f(θ,δ)=1+δθ
, and the dotted curves are obtained using
f(θ,δ)=1
. The calculations show that for small values of
δ
the linearized form for the thermal diffusivity is adequate for short time periods, but less so for large time periods.