WOLFRAM|DEMONSTRATIONS PROJECT

Thermal Diffusivity of a Sphere

​
dimensionless time
0.35
Cp
f
ρ
f
V
f
Cp
s
ρ
s
V
s
2.5
This Demonstration describes a method for determining the thermal diffusivity of a sphere. The sphere is suddenly immersed in a cooling fluid and the temperature of the fluid is observed as a function of time.
Consider a homogeneous solid sphere of radius
R
, initially at uniform temperature
T
1
that is immersed at time
t=0
in a volume
V
f
of fluid in a well-stirred, insulated tank at temperature
T
0
.
It is convenient to define the following dimensionless variables:
Θ
s
(ξ,τ)=
T
s
-
T
0
T
1
-
T
0
= dimensionless solid temperature,
Θ
f
(τ)=
T
f
-
T
0
T
1
-
T
0
= dimensionless fluid temperature,
ξ=
r
R
= dimensionless radial coordinate,
τ=
α
s
t
2
R
= dimensionless time,
where
r
is the radial coordinate,
R
is the radius of the sphere,
α
s
is the thermal diffusivity
k
s
ρ
s
Cp
s
,
k
s
is the thermal conductivity,
ρ
s
and
ρ
f
are the solid and fluid densities, and
Cp
s
, and
Cp
f
are the specific heat capacities of the solid and liquid respectively.
The heat equation, in terms of these dimensionless variables, can be written [1]:
∂
Θ
s
∂τ
=
1
2
ξ
∂
∂ξ
2
ξ
∂
Θ
s
∂ξ
,
with the boundary conditions:
Θ
s
(ξ,0)=1
,
Θ
s
(1,τ)=
Θ
f
(τ)
.
Note that
Θ
s
(0,τ)
is finite, with

Θ
f
τ
=
-3
β
∂
Θ
s
(1,τ)
∂ξ
,
and
Θ
f
(0)=0
, where
β=
ρ
f
C
p
f
V
f
ρ
s
C
p
s
V
s
is the ratio of the heat capacities of the fluid to the solid and
V
s
is the volume of the solid.
The following solution is obtained in [2]:​
Θ(τ)=
1
1+β
-6β
∞
∑
n=1
-
2
b
n
τ
e
9(1+β)+
2
β
2
b
n
,
where the
b
n
are the nonzero roots of
tan
b
n
=
3
b
n
3+β
2
b
n
.
The only place where the thermal diffusivity of the solid
a
s
appears is in the dimensionless time
τ=
α
s
t
2
R
, so that the temperature rise of the fluid can be used to determine the thermal diffusivity of the solid. Interestingly, the temperature history of the fluid is obtained without requiring the temperature profiles for the solid.