KD2 calculates its
values for thermal conductivity, resistivity, and diffusivity by
monitoring the dissipation of heat from a line heat source given a known
voltage. The equation for
radial heat conduction in a homogeneous, isotropic medium is given by (1) where T is
temperature (”ĘC), t is time (s), k
is thermal diffusivity (m2 s-1), and r
is radial distance (m). When a long, electrically heated
probe is introduced into a medium, the rise in temperature from an
initial temperature, T0, at some distance, r,
from the probe is (2) where q
is the heat produced per unit length per unit time (W m-1)
and Ei is the exponential integral function (3) with a = r2/4kt
and g
is Euler's constant (0.5772”¦). When t is large, the
higher order terms can be ignored, so combining Eqs. (2) and (3) yields (4) where
lh
is the thermal conductivity of the medium (W m-1C-1).
It is apparent from the relationship between thermal conductivity and DT
= T-T0, shown in Eq. (4), that DT
and ln(t) are linearly related with a slope m = (q/4plh).
Linearly regressing DT
on ln(t) yields a slope that, after rearranging, gives the
thermal conductivity as (5) where q
is known from the power supplied to the heater. The diffusivity can also
be obtained from Eq. (4). The intersection of the regression line with
the t axis (DT
= 0) gives (6). From the calculated t0(from
the intercept of DT
vs. ln(t)) and finite r, Eq. (6) gives the
diffusivity. Because the higher
order terms of Eq. (3) have been neglected, Eq. (4) is not exact.
However, if the slope and intercept are computed only for DT
and ln(t) values, where t is large enough to
ignore the higher order terms, Eq. (5) and (6) give correct values for lh
and k.
To verify these relationships, realistic values of lh
and k
were supplied to Eq. (2), varying both lh
and volumetric heat capacity (rcp),
and the resulting slope and intercept tabulated for t
ranging from 1 to 30 s. Plots of slope vs. theoretical lh
and ln(intercept) vs. ln(theoretical k) show
an exact linear relationship (Fig. 1 and 2, respectively) with low
cross-covariance. Figure 1: calculated lh
vs. theoretical lh ”” ”” Figure 2: calculated k
vs. theoretical k The experimental
analysis differs from the theoretical shown in Eq. (2) in that the
heater and sensor have their own conductivity and diffusivity, which, in
general, differ from those of the medium being measured. We have shown
experimentally that the relationships in Eq. (5) and (6) still allow
calculations of lh
and k,
but empirical factors must be introduced to correct for heater thermal
properties. Assumptions: The
thermal conductivity measurement assumes several things: the long heat
source can be treated as a infinitely long heat source, the medium is
both homogeneous and isotropic, and a uniform initial temperature, T0.
Although these assumptions are not true in the strict sense, they are
adequate for accurate thermal properties measurements. Further Readings: Bristow, K.L.,
White, R.D., Kluitenberg, G.J., 1994 Comparison of Single and Dual
Probes for Measuring Soil Thermal Properties with Transient Heating.
Australian Journal of Soil Research 32, 447-464. Bruijn, P.J, van
Haneghem, I.A., Schenk, J. 1983 An Improved Nonsteady-State Probe
Method for Measurements in Granular Materials. Part 1: Theory. High
Temperatures - High Pressures 15, 359-366 Shiozawa, S.,
Campbell, G.S., 1990. Soil Thermal Conductivity. Remote Sensing Rev.
5, 301-310. van Haneghem, I.A.,
Schenk, J., Boshoven, H.P.A., 1983. An Improved Nonsteady-State
Probe method for Measurements in Granular Materials. Part II:
Experimental Results. High Temperatures - High Pressures 15, 367-374.



