Thermal Conductivity of Soils
Thermal Conductivity of Soils
Thermal Conductivity of Soils
_ X Draft
..... Translation
July 1977 637
- iL
I C THERMAL CONDUCTIVITY OF SOILS:
o O. Johansen
SEP 12 1977
q=
•CORPS OF ENGINEERS, U.S. ARMY
"COLD REGIONS RESEARCH AND ENGINEERING LABORATORY
HANOVER, NEW HAMPSHIRE
""i
Unclassified
S EC
U R ITY C L ASSI FI
C AT ION O F T1HIS P A G E("' On D ataEntered)
JJ Johansen
P EIF'ORMING ORGANIZATION NAME AND ADDRESS 10. PROGRAM £LEMNT1. PROJECT, TASK
AREA, A ORK UNIT NuMaiR$
U.S. Army Cold Regions Research and .-
Engineering Laboratory
Hanover, New Hampshire
I I. CONTROLL.ING OPICE NAME AND ADDRESS .1 . ReowVOWT0Av
/iJul1.r'0
14. MONITORING AOENCY N AME 6 ADORCS(It diflftert Iom Ce.illing 0ics) S. SECURITY CL'X . (of t piw t)
17, 0IUTRIdUTION STATEMENT (.tho bettract ontord In Sleok 20, It diffren trm Repopte)
/• / ' 'j .,
S.:. tI. KEY WOROS (ContInue an ?over*@ side it n.t laff aid Identify by btlok number)
20. AGSTRACT CrI•,n•u. on rover" Oldi It n•w.cae•. ind identify by block iu•.be)
The aim of this investigation has been to create a mathematical model for
calculating thermal conductivity of soils with ordinary soil parametirs as
~J~A input data. One part of this work has been devoted to literature studies on
heat-transfer mechanisms in moist materials. These studies have made it pos-
sible to give bounds for the different domains where the various mechanisms
have an appreciable influence on the total heat transfer.
AUTHOR: 0. Johansen
IAK
Translated by Rusetta Stone, Nashua, N.H. for U.S. Army Cold Regions
Research and Engineering Laboratory, 1977, 2 9 1p.
d "'
NOTICE
The contents of this publication have been translated as presented in the
original text. No attempt has been made to verify the accuracy of any
statement contained herein. This translation is published with a minimum
of copy editing and graphics preparation in order to expedite the dissemination
of information. Requests for additional copies of this document should be
addressed to the Defense Documentation Center, Cameron Station, Alexandria,
Virginia 22314.
',.
TRANSLATOR'S
NOTES: To the translation of "THERMAL CONDUCTIVITY OF SOILS",
by Oistein Johansen.
3. Expressions .of the type "Hazens formel" can mean either "Hazen's
formula" or "Hazens' formula". The translation may contain some
errors in this respect.
Eq 27:
L•vv•._•pv•._--e '-4I
(~1b
~r
Eq 28: -A" ---
Heat flow 4
4-.4
THERMAL CONDUCTIVITY OF SOILS *
OF
Siviling. 0istein Johansen
Institutt for kjoleteknikk
703 4-Trondheim-NTH
1975
We
I~i
0istein Johansen
FOREWORD
S~iii
the G:'Qund "or the final form of the report. Work on presenta-
tions and papers has taken considerable time but also provided
has also been of great importance for the results presented here.
In particular, I wish to thank Mr. Per Erling Frivik, M. Sc.,
who has directed the work of the group, as well as other
co-workers who have performed much of the practical work in
connection with the development of test equipment and who have
given valuable help in performing the experiments.
The report was written while I was employed at the National Road
Laboratory. I wish to thaak the management and colleagues here
for the consideration they have shown despite the fact that this
work has taken more time than originally estimated.
iv
INDEX
P age
Foreword iii
Index v
Symbols viii
1. Latin letters viii
2. Numbers independent of units x
3. Greek letters xi
4. Mathematical symbols xiii
5. Often used indexes (subscripts) xiii
6. Conversion factors xiv
Bibliography xv
INTRODUCTION 1
CHAPTER I: HEAT TRANSPORT IN MOIST SOIL 3
1. Thermal Radiation 4
"2. Combined Moisture and Heat Transport 11
A. Water vapor diffusion 11
B. Capillary water transport 14
C. Contribution of vapor diffusion
to thermal conductivity 17
3, Convection 20
A. Heat Transport 20
B. Liquid transport in poro.us materials- 23
C. Forced convection 26
D. Free convection 31
4. Heat Transport in Freezing Soils 38
5. The Effects of Various Heat Transport
Mechanisms 42
REFERENCES 47
CHAPTER II: THEORIES FOR HEAT CONDUCTION IN COMPOSITE
MATERIALS 50
1. Exact Theory 50
"2. Evaluation of Limits 54
A. Limits according to Hashin-Shtrikman 54
B. Limits based on point correlations 57
C. Other evaluations of limits 62
D. conclusion 65
3. Mathematical Models 67
A. Introduction 67
B. "Exact" calculations 68
C. Approximate methods based on elec-
trical field theory 70
D. Approximate calculations based on
elements connected in series or shunt 76
E. The geometric mean model 80
F. Conclusion 83
REFERENCES 85
V
Page
CHAPTER III: SOIL PARAMETERS 88
1. Thermal Conductivity of Soil Generating
Materials 88
A. Soil generating minerals 88
B. Quartz content in loose deposits
C. Thermal conductivity for the compo-
nents 10
2. Volumetric Composition and Texture 104
A. Volume ratios in soil materials 104
B. The importance of texture 107
C. Unfrozen water 113
3. Soil Material Classification 118
A. Geotechnical analysis 118
B. Supplemental analysis 123
REFERENCES 125
CHAPTER IV: EXPERIMENTAL METHODS FOR DETERMINING
CONDUCTIVITY IN SOIL MATERIALS 127
1. Thermal Conductivity Measurements
in Materials Containing Moisture 127
2. Stationary Methods 131
A. Kersten's cylindrical test unit 131
B. Planar test unit 136
3. Thermal Conduction Probes 143
A. Influence of finite size and varying 143
thermal properties
B. The effects of moisture transport 149
C. Experimental arrangement 152
4. Other Transient Methods 157
A. Linear temperature rise 157
B. Constant thermal flux 161
C. Periodic variation of temperature 165
D. Conclusion 173
REFERENCES 174
: vi ]4
Page
3. Moist Mineral Oils 205
A. Introduction 205
B. Study by Kersten 207
C. Investigation of Norwegian soil
matefials 214
D. Conc us on 217
4. Soil Materials Containing Humus 218
A. Introduction 218
B. Measurements on peat and bog 219
REFERENCES 222
CHAPTER VI: METHODS FOR COMPUTING THERMAL CONDUCTIVITY OF
SOIL MATERIALS 224
1. Existing Methods 224
A. Analytical methods 224
B. Empirical methods 230
C. Conclusion 246
2. New Methods for Determining Thermal
Conductivity in Soil Materials 248
A. Introduction 248
B. Mathematical model for mineral soils 248
3. Effects of Variations in Soil Parameters 254
SA. Effects on conductivity calculations 254
B. Conclusions 266
4. Effects of Conductivity Variations ... 270
A. Introduction 270
B. Effects of conductivity, variations 271
C. Conclusion 275
REFERENCES 276
SUMMARY 278
POST-SCRIPT 284
APPENDIX 286
k), REFERENCES 291
vii
SYMBOLS
i. Latin letters
A attenuation ratio
a X/cp
X thermal diffusivity m2 /s
A Xe/(Cp)f thermal diffusivity, included in
M Rayleigh's number m2 /s
a,bc,d constants in Kersten's equations -
viii I'
A
gi' i = a,b,c form factor. Used in Fricke's
equations.
h thickness
i potential gradient m/m
K hydraulic conductivity m/s
k permeability m2
1 length m
n porosity
N/r 2
p saturation pressure
q quartz content
q heat flux W/m2
ql water flux kg/m2s
r radius m
Sr degree of saturation
T temperature K, °C
t time s
ix
V flow velocity m/s
V5 , V1 0 water content which gives 5 (10)
mm penetration (sinking)of cone in
consistency measurements
(relative to dry weight)
w relative water content by dry
weight
wv relative water volume
ww relative water content by wet
(total) weight
w25 1 wl0 0 water content which causes flowing
for 25 (100) blows in Lassagrande's
flow limit test-set (ratio by
weight)
wL (-W 2 5 ) flow limit
wu unfrozen water ratio by weight
(water/dry material)
I
3. Greek letters
v chinematic viscosity m2 /s
frost depth m
mathematical symbol for product,
p
a ai - a1 a 2 .... ap
i-I
0 density
mathematical symbol for sum, e.g.
V
0 relative humidity
dimension of cubical unit cell
Oi(r) function which assumes values 0 or
1 if the (position) vector r falls
outside or inside component "i",
respectively
S"' • - •--• • --- i .... "" ' , ",• • -- • '-• • = • -mxii• ,'
I,
-- 4
moisture potential (suction) N/m2 ,m
angular velocity rad/s
4. Mathematical symbols
ill ,f liquid,
liquid (cp) thermal capacitance by volume for a
e.g. f.
xiii
6. Conversion factors
xii
xiv S
BIBLIOGRAPHY
xv
J.H. Blackwell: A. Transient Heat Flow Method for Determining
the Thermal Constants. J..Appl. Phys. 25, (2),
1954, pp. 137-144.
J.H. Blackwell: The Axial-Flow Error in the Thermal Conductivity
Probe, Can. J. Phys. 34, 1956, pp. 412-417.
C.A. Bower and J.0. Goertzen: Surface area of soils and clays by
an Equilibrium Ethylene Glycol Method. Soil
Science, 87, 1959, pp. 289-292.
E. Brendeng and P.E. Frivik: New Development in Design of Equipment
for Measur•.a Thermal Conductivity and Heat Flow.
Institutt for kjoleteknikk, Trondheim, NTH, 1973.
W.F. Brown: Solid Mixture Permittivities. The Journal of
Chemical Physics, 23, (8), 1955, pp. 1514-1517.
C
J.W. Cary: water Flux in Moist Soil: Thermal Versus Suction
Gradients. Soil Science, 100, (3), lq55,
pp. 168-175.
H.S. Carslaw and J.C. Jaeger: Conduction of Heat in Solids. Sec.
ed. Oxford Press, 1959, p. 345.
J.C. Champoussin: Determination Simultance des caracteristiques
termocinetiques des solides. Int. J. Heat.
Mass, 15, 1972, pp. 1407-1418.
S.C. Cheng and R.I. Vachon: The Prediction of the Thernal Conductivity
of Two and Three Phase Solid Heterogeneous Mixtures.
Int. J. Heat Mass. Transfer. 12, 1969, pp. 249-264.
S.C. Cheng and R.I. Vachon: A Technique for Predicting the Thermal
Conductivity of Suspensions, Emolusions and Porous
Materials. Int. J. Heat Mass. Transfer, 13, 1970,
pp. 537-546.
xvi
C. Codegone et al: Resonant Thermal Waves in Insulating Slabs.
1IR. 4
Corn.x..ssion 2. (Trondheiin) 1966. (Saertrykk)
S.R. Coriell and J.L. Jackson: Bounds on Transport Coefficients of
Two-Phase Materials. J. Appl. Phys. 39, (10),
1968, pp. 4733-4736.
P.B. Corson: Correlation functions for predicting properties
of heterogeneous materials. I. Experimental
nmeasurements of spatial correlation functions in
multiphase solids, J. Appl. Phys. 45, (7),
1974, pp. 3159-3164.
P.B. Corson: Correlation functions ... 1. Empirical construction
of spatial correlation functions for two-phase
solids. J. Appl. Phys. 45, (7), 1974, pp. 3165-3170.
P.B. Corson: Correlation functions ... IV. Effective thermal
conductivity of two-phase solids. J. Appl. Phys.
45, (7), 1974, pp. 3180-3182.
D. Croney and J.D. Coleman: Pore Pressure and Suction in Soil. Fra
"Pore Pressure and Suction in Soils". Butterworths.
London 1961, pp. 31-37.
J.A. Currie: Gaseous Diffusion in Porous Media. Part 2 Dry
sc
Uanylar telfs,
pp.
Bri jurnal of Applied
3 ?h
xvii
R.G. Deissler and J.S. Boegli: An Investigation of Effective Thermal
Conductivities of Powders in Various Gases.
Transactions of the ASME, October 1058, pp. 1417-
1425.
P.B. Desphande and J.R. Cooper: Thermal Conductivity of Two-Phase
Systems. J. Heat Transfer. Trans. of the ASME,
May 1972, pp. 249-264.
E
J.W. Elder: Steady Free Convection in a Porous MediutL Heated
from Below. J. Fluid Mech. 1967, 27, pp. 29-48.
M.A. z0lsayed: Bounds for effective thermal, electrical, and
magnetic properties of heterogeneous materials
using higher order statistical information. J.
Math. Physics, 15, (ii), 1974, pp. 2001-2015.
I.M. Eltantawy and P.W. Arnold% Reappraisal of Ethylene Glycol
Mono-.Ethyl Ether (RGME) Method for Surface area
Estimation of Clays. Journal-Soil Science, 24,
1973, pp. 232-238.
S. ERgun: Fl~asigkeitsstr6mung durch Fallk5rpers~ulen
(Festhetten). Aeferert i Chemie Ing.-Techn,, 2ý,
1953, p,. 262.
xviii
G
•i, I, J
*II.
xix
-----------------------------------
1
0. Johansen: Unders~kelse av varmetransporten i kult og pukklag.
Det store ekuamensarbeid. Instituitt for kjoletek-
nikk, NTH, Trondheim, Host 1968.
0. JOhansen: Innflytelse av fri konveksjon i grovkornede
vegbyggings-materialer. Institutt for kjoietek-
nikk, NTH, Trondheim, 1970.
W.D. Joshua and E. De Jong: Soil Moisture Movement under Temperature
Gradients. Can. J. Soil. Scie. 53, 1973, pp. 49-57.
xx
R.J. Kunze et al: Factors Important in the Calculation of
Hydraulic Conductivity. Soil Sci. Soc. Amer.
Proc., 32, 1968, pp. 760-765.
1
E.R. Lapwood: Convection of A Fluid in a Porous Medium.
Proc. Cambr. Soc. 44, 1948, pp. 508-521.
P.F. Low et al: Some Thermodynamic Relationship for Soill at or
Below the Freezing Point: Freezing Po nt
Depression and Heat Capacity. Water Resources
Research, 4, 1968, pp. 379-394.
A.V. Luikov: Analytical Heat Diffusion Theory. Academic Press,
New York, 1968, p. 312.
M
M.W. Makowski and K. Mochlinski: An Evaluation of Two Rapid Methods
of Assessing the Thermal Relativity of Soil.
Proceedings. Institution of Electrical Engineers,
103, A (11), 1956, pp. 443-470.
S. Masamune and J.M. Smith: Thermal Conductivity of Beds of Spherical
Particles, Ind. and Eng. Chem. Fundamentals, 2, (2).
"Mc Gaw: Thermal Conductivity of Compacted Sand/Ice Mixtures.
,' Highway Res. Rec. 215, 1968, pp. 35-47.
R. McGaw: Heat Conduction in Saturated Granular Materials,
Highway Res. Board. Spec. Rep. No. 103, 1969,
Pp. 114-131.
R.E. Meredith and C.W. Tobias: Resistance to Potential Flow through
a Cubical Array of Spheres. J. Appl. Phys. 31,
(7), pp. 1270-1273.
M.N. Miller: Bounds for Effective Electrical, Thermal and
Magnetic Properties of Heterogeneous Materials,
J. Math. Physics, 10, (11), 1969, pp. 1988-2004.
A. Missenard: Conductivitd Thermique. Editions Eyrolles, Paris
1965, p. 356.
N
Y. Nakano and J. Brown: Mathematical Modeling and Validation of the
Thermal Regimes in Tundra Soils. Arctic and Alpine
Research 4, (1), 1972, p. 19-38.
xxi
w4
O
K. Ofuchi and D. Kunii: Heat-Transfer Characteristics of Packed
Beds with Stagnant Fluids, Int. J. Heat Mass
Transfer, 8, 1965, pp. 749-757.
R
E.H. Ratcliffe: The Thermal Conductivities of Ocean Sediments.
J. Geophys. Res. 65 (5), 1960, pp. 1535-1541.
Lord Rayleigh (1892). Referert av Meredith (1969) (18).
J.A. Reynolds and J.M. Hough: Formulae for Dielectric Constant of
Mixtures. Proc. Phys. Soc. (London) B70, 1957,
pp. 769-775. (Referert H. Fricke (1924).
R.L. Rollins et al: Movement of Water in Soil due to a Temperature
Gradient. Highway Research Board Proc. 33, 1954,
pp. 492-508.
C.W. Rose: Water Transport in Soil with a Daily Temperature
Wave. I. Theory and Experiment. Aust. J. Soil
Res., G. 1968, pp. 31-44.
S
E. Saare and C.G. Wenner: V~rmeledningstal hos olida jordarter.
Statens n&mnd f6r byggnadsforskning. Handlingar
no. 31, 1957.
xxii
14 J.H. Sass et al: Thermal Conductivity of Rocks from Measurements
on Fragments. J. Geophys. Res. 76. (4), pp.
3391-3401.
Xx~ii
(
W.O. Smith: The Thermal Conductivity of Dry Soil. Soil
Science, 53, 1942, pp. 435-459.
H. Sveian: Norske jordarters kvartsinnhold. Hovedoppgave i
Ingeni~rgeologi, Geologisk Institutt, NTH.
H~st 1972.
H. Sveian: Under~kelse av kvartsinnhold i losmasser ved hJelp
av DTA. Utf~rt for Institutt for kj~leteknikk,
Trondheim, 1973.
1. Szanto and J. Auirre Puente: Etude des caract~ristiques des
miyieux poreux fins humides.
K. Terzaghi and R.B. Peck: Soil Mechanics. John Wiley & Sons,
New York, 1948, p. 32.
J. V. Thue: Red 7sert fundamentdybde, - termiske problemer. Del I.
".nstitutt for husbyggingsteknikk. Trondheim, NTH,
p. 82.
A.R. Tics et al: The prediction of Unfrozen Water Contents in
Frozen Soils from Liquid Limit Determinations.
Symposium on Frost Action on Roads. 1. OECD.
Paris 1973, pp. 329-344.
G.T. Tsao: Thermal Conductivity of Two-Phase Materials. Ind.
Eng. Chem. 53, (5), 1961, pp. 395-397.
U,V,W
E.M.F. van der Held and F.G. van Drunen: A Method of Measuring
the Thermal Conductivity. Physica 15, (10), 1949,
pp. 865-881.
M. van Rooyen.and H.F. Winterkorn: Structural and Textural Influences
on Thermal Conductivity of Soils. Highway Research
Board. Proceedings. Washington, D.C. 1959, pp.
576-621.
L. van Zee and C.L. Babcock: A Method for the Measurement of Thermal
Diffusivity of Molten Glass. J. Am. Ceram. Soc.,
34 (8) 1951.
Vegnormaler. Statens vegvesen. Vegbygging. Overbygning med
bituminest dekke. 1970. Kapittel VI. Avan...3.
xxiv
K. Wakao and K. Kato: Effective Thermal Conductivity of Packed
Beds. Journal of Chemical Engineering of Japan.
2, (1), 1969, pp. 24-33.
N. Wakao et al: View Factor between two Hemispheres in Contact
and Radiation Heat-Transfer Coefficient in Packed
Beds. Int. J. Heat Mass Transfer, 12, 1969, pp.
118-120.
1x
, i
x'Y, z
xxvi
INTRODUCTION
Methods which have been mostly used previously to treat this pro-
blem have been based on empirical correlation between a limited
V set of mineral parameters and results from systematic heat transfer
measurements with selected minerals. This has given results with
relatively wide tolerance limits. These methods mpsy have been
satisfactory for approximate determination of frost penetration.
' However, while using modern computer programs for finding tempera-
ture distribution in the ground, researchers have felt a need for
improved accuracy of thermal parameters.
a -
4 undertaken, one has to clarify the most fundamental relations
between mineral parameters and thermal properties, in order to
arrive at the most conclusive test program possible.
2
4/
geotechnical analysis. A mathematical model developed along these
lines will consequently give increased accurary in proportion to
the information available on the minerals involved. NThen using
mathematical models it is important to know the relation between
accuracy and level of information. Part of this report is there-
fore devoted to a study of this topic. The effects of uncertain-
ties in thermal conductivity values on thermal computations is also
discussed,
CHAPTER I
3
S.
. . . . . . . . ... . ... . . . . . . . . .. . . .
•5
" M. 11 L ,r .• . r
i rur
%a ý- jw If
3, 4 t nlg
1. THERMAL RADIATION
The share of thermal radiation in the effective heat conductivity
increases with increased pore area and increasing temperature.
This can be illustrated by the simple case of radiation across a
gap between two parallel.i surfaces at temperaturei T1 and T2 .
If the emissivity of both surfaces is c, the transferred heat flucv
can be expressed as
q12 " a (T 4 _ T24) (W/m2 ) 1
r - 4o T3 (W/m2 K) 4
Wakao and Kato used this heat transfer coefficient for numerical
analysis of regular lattices of identical spheres (4,5,6). Two
packing configurations were studied, cubic lattice with n equi-
valent porosity n = 0.476 and orthorhombic lattice with an equivalent
porosity n a 0.395. Figure 2 shows the two packing configurations
1)
and the assumed heat flow is indicated by arrows.
5
, •: - •-• .- • . ..- ... ,.. . , • • """"";w" - • . •---•• :-•::: ""• •:•" . 4•
A'JBISK PAKY4NG n.Qovs OROPI5 noons
CAOiC ,ATT1CE
II
The results of their computations were presented in the form of
curves with a dimensionless N4usselt number for radiation N4ur as
a parameters
Nur " OrI"JIWO11V
dp/Xsof Iso CA.",i "m ".PkW
i;Ya 6
where ais defined by Eq 5, d~ is the particle diameter and s'
is the thermal c.onductivity of each particle. Figure 3 shows the
result obtained for orthorhombic configuration (n =0.395).
r 1.00i
-
'4
Wakao and Kato (5) have used the results shown in Figure 3 to
derive an empirical equation of the form
Xe - ÷ 1.3
+O Or ' dp (W/mK) 7
where Ae is the effective thermal conductivity of the material
when no radiation takes place. This equation was derived on the
basis of both packing configurations and is said to be approxi-
mately valid also outside the considered range of porosity
(n - 0.395 to 0.476).
Eq 7 and the results shown for Nur - 0 in Figure 3 were used for
calculating the following examples. Since the study of dry soil
is of interest, thermal properties for gas and solid material are
assumed to be those of air and granite, respectively. See Table 1.
7
S__,
, • | • .. . .*- ...... ,i "= . . ; ' ' • "• , " _ • T ••
It should be noted that the conductivity for air increases with
temperature, while the opposite occurs for granite. The same is
true for most minerals, except for some alkaline materials such
as gabbro, for which the conductivity is nearly constant or
20 -
0 o
RIO3I~
N 11nI M1 7 ri• p.l , M, L
A@ C',8
qv .T
Da 7p (kg/m 2 s)
S1)
SF, (!/ 2) 9
3
0 is the density of water (kg/m )
R , gas constant of water vapor (461 J/kgK)
T , absolute temperature (K)
9
It is evident that the relative humidity falls off significantly
only for comparativ.?ly low water content, even for finer grain
soils.
aT
However, Philip and de Vries (15) have pointed out that water va-
por and gas diffusion in moist materials are different, since
water and water vapor interact.
10
2. COMBINED MOISTURE AND HEAT TRANSPORT
A. Water vapor diffusion
100N/Mr
JLE
-- LOAM 4 sO
I-
1 i \
0 ,
0 2i
20 20 */ 40
VANNINNHOLD. VYKT
WATIERCONtN!, W?~iOmT
I'P
1±J
08
0.4[
POROSITY
D~ a 00fftitfln'
Q8 n, if~ ~f ;Q~d'od?'0,4¶nar
'01
02
Based on the previous discussion one can define two diffusion co-
efficients for vapor diffusion due to moisture and temperature
variations, D and D, respectively
qV/0Dv D v
* - ye - DT 7 (m3 /m 2 s) 11
"V V
where both DO W/o) and DT (m2 /sK) contain correction factors
V V
which account for the previously mentioned increase in diffusion
,1 rate, relative to the case of gas diffusion in comparable materials.
13
B. Capillary water transport.
Se (M/m) 13
De u K (m2 /s)
A. Medium a&nd
0 10 20 Lo io I. t•LP| silt
.,:, YAN. INNHO4LO,
V61..*/ S• it |&l,111i1
'3 2 0
3C 41: so
VAN'h44'O.O, ',4 *1.
- "BP! I6 I, R DII'III
PUlgrl i en~l'l
to l over
l or "ierUt o
i nL ' Ii....
apff
7 tV 16
According to Pick's law, the amount of diffused vapor is given by
*i.Da 17"
q a "a.~ (kg/m2.) 17
17
= rDa "p"t
RT ýT
a Xv 20
T sa.nw~'UM11".nn
l . Med,'••
1" {01'i 'I i .. x v.1
'm •) igjIInl..nN"
ing:,soynen ni'L
4V ý-Url fD)
W/m k 10 o, vwf•i 'A). Appae
0 10 20 30
A
01'
o 20 1.0 60 4
?"WPWRA!UR
18
mited, the dependence on temperature is wea, while this depen-
dence is strong for low to moderate degrees of saturation, where
the thermal diffusion is significant.
Cx ,T) •ee)
7 +rV7, 21
a2T
a•• 22.
19
c is specific heat (J/kgK)
0 is density (kg/m 3 )
3. CONVECTION
A. Heat transport
20
addition, convection between particles and liquid ma:, contribute,
as well as heat transport due to "turbulent" mixing of liquids
with different temperatures within the "pore strata". The lat-
ter effect is often termed "turbulent diftussion" or dispersion
and can be included as a velocity dependent constribution to the
thermal conductivity of the liquid (25).
, aT + (co) v VT tm 23
21
heat flow and planar horizontal liquid flow in porous materials,
similar to horizontal ground water flow in a soil layer. There
is still reason to assume that the previously mentioned increase
in the effective vertical conductivity also should occur in such
cases, provided that the flow rate is sufficiently high.
A
Goo
* 601
/- /ASCW'G RINGER'
* -• SANOSTEIN
t 0002STIROMId•NOSMASTIGMIT
000 .00rollk_0.00 2.0 IWO
h , t,.ow (21),
BESTLAVAILABLE COPY'
22
The contribution x. is given by
xt PP( 25
Xf
Pe m- d (co)±f v/xf
v - P 26
23
*@ 0 is the density of the liquid
P is the pressuire in the liquid
g is the gravitation (acceleration) constant
k 3.
5S2(1-n)e (m ) 27
S = 6/dr (ma/m 3
) a8
k ....0 (m 2)
(i0()-n)
29
d - j. m-J •i 30
24 "
As an example, Hazen's formula for water permeability of loosely
packed sand consisting of particles of equal sizeI) reads (28)
2
K = 116 d (0.7 + 0.03t) (cm/s) 31
v =
* 32
25
velocities, the first term wvill be negligible in
comparison to
the second. The equation will then approach the form suggested
by Burke and Plummer (30), where the flow velocity is proportio-
nal to the square root of both the pressure gradient and the par-
ticle size.
d-
J C. Forced Convection
In the case of flowing ground water, the liquid flow will usually
be nearly perpendicular to the direction of the heat flow. The
266I
contribution to heat flow'by the liquli can thor. be ascribed
principally to the previously mentioned dispersion. If one
assumes that the case of vertical heat flow in the presence of
(horizontal) planar liquid flow is analogous to the case of
radial heat flow in a column with axial flow of liquid, the con-
tribution due to dispersion -an be derived from the empirical
relations of Yagi and Kunil , established for the latter case,
see Eqs 24 and 25.
For sandy materials one finds that the effect of ground water
flow is slight except in the case of very coarse sands, where
the dispersion effects may increase the conductivity by up to
20 percent. For still coarser materials, such as gravel or rock,
this effect may be of an entirely different magnitude. However,
ground water flow seldom occurs in coarse gravel, due to its ex-
cellent drainage properties.
m I
.I o b Pi
;I 00,0 saila
io.
1*l,"
27
ioo
too istI 0
•t MANN M•/t
U tI /
I.*tll
%5,
1?, Ie.ativ andel hy dim.rsimonun til 4n verikal.
111n tMLfeznsofm I rJo ar•e .JqnnoIfstrpAaMt \v Yuhn.
The two fans located at one end of the rock bed gave an air velo-
city of about 2 m/s across the fill. The air temperature was
28
10C, while the temperature at the bottom of the fill ranged
from 2 to 4*C. Figure 17 shows the isotherm surface for 9*C,
obtained through temperature measurements in the fill after
stabilization. Temperature profiles along the center line cut
0of the rock fill are shown in Figure 18. Both these represen-
tations show that there may occurred a strong influx of room
tkis of tn v 6aei-rv.
rock
.S
1~.L L~t5ar%tVQ0! (IOAi oool t u5. ~
2.
..
yarrl4og.
i.r.tf
Neate.
Z ah
pa n@.
e ';.9
,
j -sd.
/k.mj1ror.
;.-r.."orntod peso.,
air in the upper layers of the fill and that this influx helps
to maintain a temeprature close to that of the air in the room
deep down in the rock fill. The average effective conductivity
of the fill was determined to be 1.46 V/mK, as compared to 0.45
SW/mK for the same conditions but with no forced air current.
This illustrates that even small wind velocities result in large
changes in temperature distribution and heat transport in coarse
rock fill. For high wind velocities, the demonstrated air in-
flux can practically eliminate the thermal resistance, particu-
larly if the air can penetrate completely into the fill. In
the wintertime, the outer surfaces of such fill may be
29 " f
covered1) by ice or snow, but during the autipnnonths and in cold
periods with no snow cover this mechanism can result in extreme-
ly heavy loss of heat from the ground, The importance of cove.
ring top and side surfaces of coarse fills wvth some-what finer
gravel or crushed rock is obvious.
I3
30
D. Free Convection
am (c0)ff
31
Lapwood (33) has calculated critical Rayleigh's numbers for
layers of liquid-filled porous materials for various boundary
conditions:
32
20 to 80 mm. The rock was rather loosely packed and had a dry
density of about 1,500 kg/m3. The layer was 0.48 m thick and
covered an area of 1.8 "2.2 m. The arrangement in the refrige-
rated room is shown in Figure 16, see the previous section,
During the free convection experiments, a cover plate was in-
stalled at a certain height above the layer to protect it from air
movements caused by the air cooler. As stated before, heat trans-
port through the layer was measured with 9 heat flow gauges even-
ly distributed across the area, while temperature distributions
were recorded in 7 horizontal cuts, spaced 8 cm apart and with 22
sensors in each. The latter allowed generation of temperature
profiles and isotherms within the rock layer.
33
Table II: Experiments with free convection in crushed rock (32).
6 m2
Permeability K a 0.97 *10"
i3
34
S,--
i• '.-•
, ••k -,...
• • m. .. '--?" '• r.j• •,-••.•, .,•.,•,• . . . ..... •-.:... A.
values for the Rayleigh number correspond to temperature differences
of 7.8 and 11.60C for exposed and covered upper boundary, respec-
tively, with an average temperature of O°C. For thicker layers,
the critical temperature difference will be less, since Tk is in-
versely proportional to the thickness. The effect of this in the
thermal resistance of rock fill is outlined in Figure 20. The ma-
terial is assumed to be dry and to have the same permeability as
the rock used in the experiments.
It is clear that doubling the depth, e.g. from 0.5 m to 1.0 m, may
in the worst case result in the same thermal resistance, provided
that the temperature difference is above the critical limit in
both cases.
A-1-
71 /"
L I -
/ - '"
t 0 ..... 4 04 r r i ...
,,g * u O SO ld 10
* td t
• l L<
e, IIiO'rI..ad
fio
d~lI~iije * 1 ,J@ i
I-JP .
l•35
P (15J
.
35
*
I-I
Thisill-%p
wil cas a lowrin ofte teprt ure' at te bondar be
eatihe
prgrses u-airs temperaturewilnrmallyntedt inhies creaserdue
Jt uayo
thenyearl var Al temperature difference abov
ition, cded,
tex -
XO c ~ 0 a
W
....
J/k$oC ..kg/.
.........
1.25.. 3.
m../s 0-... ...
Dry 1/273 0.45 1000 . .6 o'4 3.6 10o6
Sa~turated 2 1 1.50 4187 1000 0.36 10-6 1.792 10"6
with water
*) tmm 20aC
37
.. hTEAT TRANSPORT IN FREEZING SOILS
S(t) = 2F '37
t
F t (T- To)dt
0
where To is the freezing point
t is the time at which freezing occurs
38
IL. more fine grain materials, the assumption that all of the pore
water freezes at the normal freezing point for water will no longer
be valid. As an example, Figure 21 shows the portion of un-frozen
water in some silt and clay materials at temperatures below the
(normal) freezing point. The finer the grain, the more water will,
as a rule, remain unfrozen at a certain temperature. This sub-
ject will be discussed in more detail in Chapter II. Only the
effect of this condition on heat transport will be discussed here.
3i
20d
~3
YI10 2. Eksemnler at raiatt vann~rnhQ'.J i Yanr%1V1V
e•
O
4t f *
1-e g " -i
120, ,-
10,4
40/
41
as well as (liquid) mass balance versus time, through measurements
of both heat flow and water transport velocity.
The water flow towards the freezing zone referred to above will
clearly tend to reduce frost penetration and will to some degree
counteract the effects of un-frozen water content.
For fine xrain materials, such as silt and clay) there is a region
at the lower saturation levels where moisture content gradients
will cause vapor diffusion due to changes in relative humidity in
the pore air as function of water content. From Figure 5 in this
chapter one can conclude that large variations in relative humidity
42
T piau tIoc$.t ~wt~jo onm
0' hea
of~ in '0i
44
~.. ... i'
permeability must be about 0.5 -10-6 m2 for convection to occur
at a temperature difference of 100C. Permeabilities of this high
magnitude will only be present in very coarse gravel materials or
crushed rock. By comparison, a permeability of 10-9 m2 is requir-
ed for free convection to take place at the same temperature dif-
ference in a material saturated with water and for an average tem-
perature of 20C. Such permeability values can be expected in
medium sand with equal particle size. However, up-ward heat flow
will primarily occur during the cold (frost) season when convection
in water can not take place.
Forced convection due to ground water flow may affect heat trans-
port in materials having permeabilities in the range just mention-
ed (see Figure 15). For axial liquid flow in particle filled
* columns having radial host conduction (flow) it has been demonstra-
ted that the flow rate (velocity) affects radial heat transport,
which may be explained by "lateral mixing" or dispersion. Corre-
sponding studies involving horizontal liquid flow in fills (layers)
having vertical heat conduction (flow) are not known (to the author).
S..However, such conditions will normally not occur in conventional
road structures due to flow within the coarser materials in the
* structure itself.
Ii
45
S~_ ____ ____ __-A
effective heat conductivity by 5 percent in dry gravel with a uni-
form particle size of 10 mm (see Figure 4). This contribution
increases in proportion to particle size. Reduced pore cross-section
due to a certain water content and an increase in conductivity due
to increased saturation levels will tend to reduce this contribu-
tion in moist materials. The boundary of the region where ther-
mal radiation is important was derived under certain assumptions
regarding the effe:t of saturation level on the effective pore
cross-section and the increase in conductivity due to increased
degree of saturation. The boundary corresponds to a contribution
to the conductivity, due to heat radiation, of 5 percent.
46
S..
. . •. ....i.....•j ... . "' nm
Iq lua ~la l~m , ,,,• • lilii wh• .. ..... 4
REFERENCES - CHAPTER I
47
"----. --
---------------------------
,'-,-r---
16. W. 0. Smith: Thermal Conductivities in Moist Soils. Soil Sci.
Soc. Am. Proc, 4, 1939, pp. 32-40.
17. 0. Krischer: Die Leitflhigkeit des Erdbodens. Beiheft zum
Gesundheits-Ingenieur, Reihe I Heft 33, MUnchen (Munich) 1934.
18. C. G. Gurr et al: MOvement of Water in Soil due to a Temperature
Gradient. Soil Science 74, (5), 1952, pp. 335-345.
19. H. L. Rollins et al: Movement of Water in Soil due to a Tempe-
rature Gradient, Highway Research Board Proc. 33 1954, pp. 492-508.
20. J. W. Cary: Water Flux in Moist Soil: Thermal Versus function
Gradients. Soil Science. -CO, (3) 1965, pp. 168-175.
21. W. D. Joshua, E. De Jongi Soil Moisture Movement under Tempera-
ture Gradients, Can. J. Soil.-Scie., 53, 1973, pp. 49-57.
22. R. J. Kunze et al: Factors Important in the Calculation of
Hydraulic Conductivity, Soil Sci. Soc. Amer. Proc., 32, 1968,
pp. 760-765.
23. D. A. de Vries: Simultaneous Transfer of Heat and Moisture in
Porous Media. Transactions, American Geophysical Union, 39 (5)t
1958.
24. D. A. do Vries: Hat warmtgeleidingsvermogen van grond, Med.
Landbouwhogeschool, Wageningen, 52, 1952.
25. P. Adivarahan at al: Heat Transfer in Porous Rocks, Society
of Petroleum Engineers Journal 1962, pp. 290-296.
30. Burke & Plummer: Referext i Chemie Ing.-Techn., 25, 1953, p. 262.
31. 0,. Johansen: Unders0kolse av varmetransporten i kult og pukklag.
Det store ekuamensarteid, Inst. for kjoleteknikk, NTH, Trondheim,
1970.
32. 0. Johansen: Innfl teise av fri konveksjon i grovkornede
vegbyggings-materialer. Institutt for Kj~leteknikk, NTH,
Trondheim, 1970.
48
-- ... ~-
33. .R:. Lapwood: Convection of a Fluid in a Porous Medium.
c3 Cambr. Soc. 44, 1948, pp. 508-521.
34. J. W. Elder: Steady Free Convection in a Porous Medium
Heated from Below. J. Fluid Mech., 1967, 27, pp. 29-48.
35. E. Palm et al: On steady convection in a porous medium.
Preprint series, Matematisk Institutt, Universitetet i
Oslo, 1971.
36. Sv. Skaven haug: The Design of Frost Foundations. Norwegian
Geotechnical Institute, Publ. 90, Oslo 1971.
37. A. F. Knutson: Upubliserte resultater, 1973.
38. E. Angen: Analyse av vanninnhold i vegkonstruksjoner, Intern
rapport, SINTEF, NTH, Tronheim, 1973.
39. A. K. Fleming: Applications of a Computer Program to Freezing
Processes. Proceedings of the XII International Congress of
Refrigeration, Washington D.C. 1971.
49
CHAPTIM II
Many of the methods (to be) mentioned are developed for determining
transport coefficients which are analogous to thermal conductivity,
e.g. parameters such as electrical conductivity, dielectric constant,
magnetic permeability, permeabilty associated with liquid flow in
porous materials and diffusion coefficients. To simplify matters,
all these analogue transport coefficients will in the (following)
text be denoted by "conductivity", without specific mention of
the nature of the transport coefficient. That will be evident from
the list of references.
1. EXACT THEORY
5o
Ideally (in principle), this problem can be solved by inteqra-
ting the heat flow equation for steady state (stationary) condi-
tions over a region of the material, large enough to give a re-
presentative picture of the average microgeomnetrical conditions(i).
The volume average for the function oi (rj will correspond to the
relative volume (volume share) of the component "il" in the material,
provided that the integration is carried over a sufficiently large
volume in a macroscopically homogenous and isotropic material.
That relative volume also gives the probability that a point selec-
ted at random will fall in component "1'02).
51
In the same manner, a volume average of *the type <Pi ('r)1
S(r),.....
0 (' *)> expresses the probability that all the
points r1 through rp are located within the same component. Such
a volume average is called a p-point correlation. As an examp2e,
a two-point correlation will define the probability that both end-
-points of a line segment fall within the same component. In a
macroscopically homogenous and isotropic material, this probability
will be a function of the length of the line segment. Three-point
correlations will give the probability for the three corners of a
triangle to fall within the same component. This probability is
a function of the angle between two sides in the triangle and their
lengths.
1 * 7.A A (30o 1
XB
2 XI (conductivity of the components)
52
tivity of each component and their relative volumes. This may
appear to be a trivial result, but has nonetheless been over-
looked in many studies within this field: "The flow of papers
seeking an1) universal formula (for the conductivity of compo-
site materials) is unending"t. Quote from P.J. Beran (1971)(4).
" _( -n) .
1) Sorry about this one, but that is exactly what the original
says. (Translator's note)
2) Probably some-what overstated (Translator's note).
53
However, in a more recent report, Brown has pointed out the possi-
bility of determining limits for the conductivity in a certain
system, based on knowledge of the p-point correlations up to a
given order (5). The more information one has of the statistical
properties of the microstructure in the material, the narrower
these limits are.
2. EVALUATION OF LIMITS
SXH X A 3
F Awhere
A E n
ixI (arithmetic mean)
E ni/Ai)"
E.( ( harmonic mean)
The proof for this (relation) is usually based on (the work by)
Wiener (1912)(6), so the limits could be called Wiener limits.
These limits are general and will apply even if the material is
not homogenous or isotropic.
54
Using calculus of variations, Z. Hashin and S. Shtrikman (1962)(7)
developed limits (limiting values) for macroscopically homogenous
and isotropic materials, which represent a further refinement (im-
provement) of limits for the conductivity of a composite material.
2 X1 - ý2~ S
m 7
jME56N
B. Limits based on point correlations
x . ~.08
,0 ,I•"• i '..
.. •I•. LoI, EI0 ,, .,IId3'iB ,it,.'-,o
'0 0 2 04, "S ý"S'•1`0 '0 -. ' 10 0 0
'" .. •...
.... ~I t•
" Z ... -*' '
57dotit%,l
N
7OP
An article from 1974 by P.B. Corson (9) gives the results of ex-
perimental determination of three-point correlations in selected
alloys (lead-aluminum, iron-lead). From these experimental re-
sults, Corson derives empirical formulas which express the three-
point correlations for the investigated materials as a function of
the length of the three sides in the triangle and the relative
volume for the components.
58
same average geometry, while in the other types the component
geometries are different. The most important restriction for the
cell materials is that the conductivity of one cell must be sta-
tistically independent of the conductivities of all other cells.
In addition, the material must be homogenous and isotropic in a
macroscopic sense.
r -.
"* .( ~ -'
r'4.** ' .5'. -: : i r):I
where X /X I X2 '1
59
Porous materials (G1 known)
One finds that the largest improvement is obtained for small con-
centrations of particles or pores, where both Miller limits asymp-
totically approach one of the Hashin-Slitrikman limits.
f/
60
, * -I | -t T " '
CAA.
0 00 . 00 000 0 0 0
B:SI AVA'LABL cO
1TY
., I Q•
C)• : '. Os 0 00 10 1.0 0• i o 0
k. 1
61
thus contribute to better accuracy for estimates of conductivity
in certain systems.
.on
Si,,. • • ,, S ,. .- €
The general expressions for these limits are given by the follow-
ing integrals:
62
As illustrated by Figure 3, these expressions are obtained by
assuming that each elementary cell is divided into an infinite
number of filaments parallel tojor plates perpendicular to the
direction of heat flow, respectively. The relative volume (share)
of the lower conductivity component in each of these filaments or
plates can be regarded asrespectively, one and two-dimensional
porosity.
-i /
"4.
A.1
I/
63
- 1,•-. • 1 . • .. : ., ...... -...
:.. .....-i .... .-..-- I,..-*w .- -1,..v .- q• . i---. i -m•
. 4.i,
•
In principle it may be possible to improve on the limits (further),
provided that one can account for the statistical distributions
of one and two-dimensional porosities. However, no such efforts
have been reported in the literature.
For the case where the (net) conductivity is assumed known for
/ 125, the Prager limita fall outside the Hashin-Shtrikman
limits over part of the range. It is also evident, that the Prager
limits diverge rapidly for values of x 2/AI higher than that where
the conductivity is assumed to be known.
64
porosity). The limits (bounds) for the conductivity are here
given as the lowest upper limit and the highest lower limit for
the two cases, respectively. As shown by the Figure (5c), the
resulting limits (bounds) constitue a very substantial improve-
ment relative to the Hashin-Shtrikman limits. For conductivity
ratios below 125, the coaductivity is defined within a narrow
band with a doviation from the average which is less than t1o per-
cent.
D. Conclusion
65
(.tpdrt'~dM~a1 uo!ud givenl' It ',/X 2 , M zi 0,40
ko,,
5 4,070 4C 6
7, . 0 0
!o~s.twoah
'e~e.vrd ed ýuei,! a! a nsa!prer
In
such casess onl _ exermetalivsiainswleuti
by mans
f th Praer lmits
Bet/a/lbeCP
66
3. MATHEMATICAL MODELS
A. Introduction
67
B. "Exact" calculations
68
17
,,• , 'I// I
'-.. 4, 1'-,4-L'L.
o.lot
69
,
69.
",j
These results show that the expected variations in conductivity
of granular materials with dispersed particles is considerably
less than the theoretical limits (derived) for spherical particles.
The numerically determined values all fall within the empirical
limits.
Based on these studies one can state that the conductivity of ma-
terials containing spherical particles can not be expressed as a
function of porosity and conductivity ratio only. In addition to
these parameters, the geometrical (packing) configuration is also
significant. However, the latter parameter has less importance
(influence) than indicated by theoretical limitsderived for sphe-
rical particles. However, these computed values are valid for a
rather narrow range of porosities. This condition, as well as the
relatively wide incertainties that may be expected, mean that an
empirical determination of conductivity will be much superior to
this (analytical) method in most cases.
.X<aT> 14
70
<VT> = n <vT>I + (1-n) <"T >a 16
cq> = X <c7T> 17
<q> 2 = xp <vT 18
71
14 Inserted in Equation 19, this results in the following expression
for the conductivity:
(X + 2)n 12
(-n)
(2. 2X 1 )n + 34XI2 ( -n)
11
a,b,c I + (A 2 /X 1)gj
discshaped particles: C a 11 gb go a 0
72
1 += 2 " 1 )" (1-n) < T, 2 / ev T> 23
<VT>2 / <VT> 24
n + (1-n)
+-*-.-A 0 25
Not even this assumption will be valid for the case of concentra-
ted particles, since the temperature field around each particle
(then) must depend on neighbouring particles.
n_ 26
nY2 '
Also in this case it is in effect assumed that the field inside
the few added particles does not depend on adjacent particles.
73
a given article geometry it is evident that a general expression
for conductivity in a material containing particles cannot be
formulated. This was also shown, by means of a study of limits,
in the previous main section.
74
•... ... ...* .. l +r m •m,• i • • A.
which gives the conductivity as a function of conductivity ratio
for constant porosity. The result of sucha fitting to a numeri-
cally computed functional dependence of the conductivity for shperes
in an orthorombic (packing) configuration is shown in Figure 9.
This discussion shows that equations of this type only have theo-
retical relevance for rather low particle concentrations, where
the theoretical limits for a given particle shape coincide. In
other cases, they may give acceptable results when the conductivity
ratio is low and the theoretical limits for the particle shape at
hand define a relatively narrow range of conductivities. For cases
other than these, the equations are only relevant as semi-empirical
(mathematical) models which can be fitted to experimental conduc-
tivity curves within limited ranges of conductivity ratios and po-
rosities, as well as for a certain choice of geometrical (form)
factors.
7,+
11 IIV
5,kLt olidl4
q Ij 6t Ol
"1!4.11 I,
PIL
J * 0 ! t) , u0 ol II 0 0• p L O ., '17
715 "'dl
D. Approximate calculations based on elements
connected in series or shunt
... ~~
', . .' ~I
S(1
11
-0) / X+
1 ( - ) /27
76
/ /x (0-0) / + X2 ) 28
These equations represent lower and upper bounds for the conduc-
tivityrespectively, for this particular (type of) cell. These
bounds are identical to those given by Coriell and Jackson, dis-
cussed in the previous section. In this case, the difference be-
tween limits (bounds) is relatively small, but for other configu-
rations (e.g. that of Figure 10 D), the difference may be signifi-
cant.
Kv t R4, ,
0J J
77
A variable geometry parameter) such as 0 in Figure 10p was used
by " ii eteal. (33,34,35). These reports give the geometry para-
meter in graphical form, as function of conductivity ratio, for
two porosities. The authors recommend linear interpolation for
other values of porosity. The cell geometry is selected on the
basis of approximate calculations of conductivity in a (packing)
configuration of npheres, using the principle of parallel heat flow
lines. By considering the number of contacting points in other
types of configurations, an effort was made (by the authors) to ex-
tend these results to other values of porosity.
7 8a
pirical correlations derived from experimental data for the mate-
rial at hand.
79
The agreement with experimental data gives this model a certain
credibi±)1y as a toll for computing conductivity in granular and
dispersed materials. However, the uncertainty becomes large for
relatively small to moderate porosities, as shown in the corre-
sponding figure. In addition, the model does not allow fitting
to moderate porosities, as shown in the corresponding figure.
In addition, the model does not allow fitting to experimental func-
tion and is thus worthless as a semi-empirical relation.
80
individual elements in the lattice was made random. For each such
configuration, the effective conductivity was calculated by nume-
rical solution of the heat flow equation in the three-dimensional
lattice. The configura'ions were changed a sufficient number of
times to give the most probable (expected) conductivity with rea-
sonable accuracy
The geometric mean equation has also been suggested by other authors
as a model for conductivity calculations for two component materials.
81 i
•••••••••••••••••••••••
For example, W. Woodside and J.H. Messmer (1961)(30) found that
this model gave good agreement with conductivity measurement in
water-saturated sand at various porosities. Also, J.H. Sass et.al.
(1971) recommended this equation for use in connection with water
saturated crushed rock (41).
F
&
Vachon, de Vries and Bruggeman) have been shown by many to give
good agreement with experimental conductivities,
vity ratios (23, 31, 42).
for low conducti-
82
i. 4641
-W
v.v• • •
: .
." tuni;.'s
'. equa iu.
0 .,. 3*S. -
Fla.
. 2. S . 1-r..Ing melllmP maocellcr r'or ,-. *r.';
er ¢nir.. av 'drinlrevv'fl
F. Conclusion
For low conductivity ratios (<10), where the different models give
almost identical values, it seems logical to prefer the simple
geometric mean model as a basis for conductivity estimateso
"li4
3'
i8
REFERENCES -. CHAPTER II
85
14. J. L. Jackson, S. R. Coriell: Transport Coefficients of
Composite Materials. J. Appl. Phys. 39, (5), 1968, pp. 2349-2354.
15. S. R. Coriall and J. L. Jackson: Bounds on Transport Coeffi-
cients of Two-Phase Materials. J. Appl. Phys. 39, (10),
1968, pp. 4733-4736.
16. S. Prager: Improved Variationa'l Bounds on Some Bulk Properties
of a Two-Phase Random Medium. J. Chem. Phys. 50, (10),
1969, pp. 4305-4312.
17. Lord Rayleigh (1892). Referent av Meredith (1969) (18).
18. R. Z. Meredith and Array
C.W. of
Tobias: Resistance
through a Cubca- Spheres. J. Appl.to Phys.
Potential Flow
31, (7),
pp. 1270-1273.
19. D.A. Do Vries: Hot warmtegeleidingsvermogen van grond. Med.
Landbouwhogeschool. Wageningen, 52, 1952.
20. G.K.H. Walther: Zgr makehemAtischen Behandlung der Stationaren
Waireleitung in stuckwelso homogenen Median. aernforschungs-
anlage Juilich, 657-MA, 1970.
21. R.S. Deissler, J.S. Boegli: An Investigation of Effective
Thermal Conductivities of Powders in Various Gases. Transac-
tions of the ASME, October 1958, pp. 1417-1425.
22. N. Wakao, K. Kate: Effective Thermal Conductivities of Packed
Beds. Journal of Chemical Engineering of Japan, 2, (1),
1969, pp. 24-33.
86
29. W.O. SMith: The Thermal Conductivity of Dry Soil. Soil
Science, 53, 1942, pp. 435-459.
30. W. Woodside, J.H. Messmer: Thermal Conductivity of Porous
Media. I. Unconsolidated Sand., J. Appl. Phys,, 32, (9),
1961, pp. 1688-1699.
31. R. Mc Caw: Heat Conduction in SAturated Granular Materials.
Highway Res. Board., Spec. Rep. no. 103, 1969, pp. 114-131.
32. G.S.G. Beveridge, D.P. Haughey: Axial Heat Transfer in
Packed Beds. Stagnant Beds Between 20 and 7500 C. Int. J.
Heat Mass, 14, 1971, pp. 1093-1113.
33. S. Yagi, D. Kunii: Studies on Effective Thermal Conductivities
in Packed Beds. A.I. Ch. E. Journal 3, (3), 1957, pp. 373-381.
34. D. Kunii, J.M. Smith: Heat Transfer Characteristics of Porous
Rocks. A.I. Ch. E. Journal, 6, (1), 1969, pp. 71-78.
35. K. Ofuchi, D. Kuniii Heat-Transfer Characteristics of
Packed Beds with Stagnant Fluids. Int. J. Heat Mass Transfer,
8, 1965, pp. 749-757.
36. S. Masamune, J.M. Smith: Thermal Conductivity of Beds of
Spherical Particles. Ind. and Engl. Chem. Fundamentals, 2,
1963, pp. 136-143.
37. G.T. Tsao: Thermal Conductivity of Two-Phase Materials. Ind.
Eng. Chem. 53, (5), 1961, pp. 395-397.
38. P.B. Desphande, J.R. Cooper: Thermal Conductivity of Two-Phase
Systems. J. Heat Transfer. Trans. of the ASME, May 1972, pp.
249-264.
39. S.C. Cheng, R.I, VAchon: The Prediction of the Thermal Conduc-
tivity of Two and Three Phase Solid Heterogeneous Mixtures,
Int. J. Heat Mass Transfer, 12, 1969, pp. 249-264.
40. J.E. Warren, H.S. Price: Flow in Heterogeneous Porous Media.
Society of Petroleum Engineers Journal, Sept. 1961, pp. 153-169.
41. J. H. Sass et al: Thermal Conductivity of Rocks from Measure-
ments on Frarents. J. Geophysical Research, 76, 014), 1971,
pp. 3391-3401.
42. K, Gotoht Thermal Conductivity of Two-Phase Heterogeneous
Substances. Int. J. Heat Mass Transfer, 14, 1971, pp. 645-646.
87
CHAPTER III
Soil parameters
8'
sits (sand, coarse silt). In addition, finely crushed feldspar
is particularly prone to chemical withering, which gives rise to
new types of minerals. These minerals will occur in the form of
rather small flakes, rich in Aluminum, while the content of the
more easily soluble elements Na and Ca is reduced. Hydromica
is such a residual product while different types of clay (kaolinite,
montmorrillonite) represent an even more complete change during
which all Na or Ca has been lost (1).
quarts 12 35
Feldspar 59.5 16
Mica 3.8 15
Miscellaneou 24.7 34
I WA
W".4 %A
89
In Norway, most older loose deposits were removed during the ice
age. Present soils were primarily deposited during the last mel-
ting period, between 8,500 and 12,0OO years ago, or in later times.
Thus, the previously mentioned (chemical) changes in the most fine
grain materials have had a short time in which to occur and Nor-
wegian clays in general contain small amounts of the typical clay
minerals, except for hydromica.
90
The most common minerals occurring in the loose deposits will
thus, as indicatedp be quartz, feldspar, mica1 amphibolite and
pyroxene. The firstthree of these will normally dominate, while
in limited areas amphibolite and pyroxene may occur in signifi-
cant amounts, together with feldspar.
For a small number of selected samples, the quartz content was de-
termined for several fractions (by particle oize) of the samples.
Based on tnese experiments, one can plot quartz content versus
91
0 0 0 0 C) 00O 0
0 0
o - tt a - M. C.~--t p u-O0 0 0 Q u
y4~
u to
Id 0
(Oi
u0
* 0 0)
hotot\N~O 00000000o 000010 9
-0 0
8 14
4
$4 0
tho I U
0 04 0 -t 0 4-
.0 5 k .
00 0 t
(00
00
(0
6 * 0(0 cc H .w 1
4J (0 0 0
N N I, co0
92~ 4
grain size for these materials. These curves are shown in Figure
3. Most of these samples were sand or gravel materials with low
content of silt or finer fractions. The results from these samples
are drawn as bold solid lines, while results from samples of moraine
and lake-shore sediments having a more even particle size distri-
bution are shown as broken lines. The thin solid lines define the
range of quartz content variations in more fine grained materials
found by Selmer and Olsen (3).
SIT ,6!! SAND
I!ý I
*OEN IVtOADR
40/ . ....
.
i ./ . .- ** /"**•, -- '- '" *- .•-•""- -
// -0',..., . . -•r:
_ . . . .. __ I.
Jr.,
IL,*-d'L ;s, :'L n... I : , &f" * +,uVp,..j
L l? 1 . 'j• (S'a, ' ;j|I'j/• (a.'
All told, this study shows that the quartz content in Norwegian
loose deposits can vary between wide limits. These limits de-
pend partially on the particle size distribution within the ma-
terials and partially on variations in mineralogical composition
of the original (rock) material, The first of these "mechanisms"
93
will cause fine grain materials with relatively high content of
materials finer than average silt (d< 0.02 mm) to, in general,
have lower quartz contents than coarser materials where silt and
finer fractions are scarce (d>0.06 mm). In regions where the
bed-rock contains little quartz, the second mechanism will come
into play and cause low quartz content for all fractions.
*, / ,,-/
, 'o rV. / .. ,,
i , /
,• ; // /",
.a
20
i 0 'h li .
41A111
I.
O-
1
1,/
'.N
.2 *I
r.IAt, t
f"A'. 'I.E ' . 1 1 I. '
ii /,Vri
!IVm. f . bV1' W- 1 '.
rra ,.u i,:,c ji :
l.'e .
l 4YJ,:ýQI 0:, J., 14t,ýJ fph aP ~j'P V
94
C. Thermal conductivity for the components
N OR
lGTA In..i.
ii
95
Figure 6 shows how the conductivity of some other soil materials
vary with temperature. The contribution from the temperature do-
pendent vapor diffusion is included (aee chapter I). Little data
exist about the conductivity of organic materials which may be
present in the soil (plant residue, humus). However, Watzinger
(10) gave a value of 0.45 W/m°C for the solid organic components
in peat or bog* The temperature dependence is not known.
i,. ~~~W/mK<.
wimxg
tO
0,0 2
T 1,0.00421 0C.
TINP,'AT'.R4
AV
K
:Table III shows the distribution of results for the most important
* rock-generating minerals. As can be seen, the variation is rela-
tively large for the same minerals However, both feldspar and mica
show relatively small variations around a common mean of about
2.0 W/mOC. For pyroxene-amphibolite the variations are more sig-
nificant, with means of about 4.4 and 3.5 W/moC, respectively.
quartz has a notably higher conductivity of 7.7 WV/mC. The latter
result is also in agreement with measurements by Birch and Clark
on pure quartz, provided that the following rule (formula) is used
for computing the conductivity of randomly oriented crystals
"97
tions are lowest for the most common minerals. Thus, the errors
in conductivity estimates for a rock material, based on these
averages, should normally be relatively small.
.98
9I.
t4
C)
~~~ 40''
'0v.0 0 0 00 0
O0 N 0 4N N N 0000000
. I.l
5.4*4),
Nu~l,4U 0 0 00 00 00 0 00 0
A 0 X's N UI a% N UIA I ON E% I N
99
'pý-N N 4 N N - ý N N U -~
4- -
(1
U 04
0 UN \UND Lý
N 00
q4 ~' 0 .
0o a*
oc
00
.44J
4- 4
100
W~4
41
104
1.101
Among other things, it is evident that rocks containing quartz
have small or no quantities of pyroxene, amphibolite or olivine
These can thus be approximately represented as a two component
system consisting of quartz and feldspar - mica with conductivi-
ties 7.7 W/mK and 2.0 W/mK (feldspar, etc.), respectively. The
conductivity of such two component system can be determined ap-
proximately from the geometric average equation
• • - ~,7.7 qq 2.01-
-q
7 7 32 0
Figure 7 shows results obtained by Birch and Clark for rock ma-
terials containing quartz, plotted versus the relative quartz con-
tent. The values obtained from Eq 2are drawn in for comparison.
The data indicate that knowledge of the quartz content will make
it possible to estimate the conductivity of quartz-bearing rocks
with errors not exceeding 20 percent.
j10
C''h
2C el
102
A corresponding mathematical model could be used for eruptive
rocks not containing quartz. One component would be feldspar-
mica (conductivity 2.0 W/mK) while the other would be pyroxene,
amphibolite, etc., (with a conductivity of 4.0 W/mK). However,
the large variation in conductivity within the latter group will
result in relatively large uncertainties.
For loose deposits where the quartz content is un-known, one may
use the results obtained by Sveian to estimate quartz content
based on the grain size distribution in the material. In the pre-
vious section it was suggested that the following values may be
used
d <0.002 mm 0 percent
0.002,d 0.02 mm 15 percent
d 0.02 mm 45 percent
PFl • 40 2,,
"0 -0
0 103
As shown in Figure 4 of the previous chapter, this mathematical
model normally gives uncertainites for the quartz content which
are less than 30 percent. This means a some-what higher uncer-
tainty (about 40 percent) in the estimated particle conductivity,
when calculated as a weighted geometric mean o-f conductivities for
quartz and feldspar/mica. To this one must add the uncertainty
due to any occurance of minerals such as pyroxene, amphibolite,
etc. Overall, this method will therefore give an uncertainty in
the particle conductivity of up towards 50 percent of the esti-
mated value, ise. twice that which seems possible to obtain when
the quartz content is known.
104
the sand fraction and about 2800 kg/m 3 for the finest fractions
where the content of mica and heavier minerals usually is larger
(see Table V). For organic materials (plant residues, etc.), Y8
will vary between 1460 and 1: Ig/m 3 (13). In soils containing
both organic materials and minerals, the volume ratio between the
two types of material is given by the average specific weight
of solid components. Thus, the relative volume (share) of organic
miterials is given by
Y Ys 7
0v • Yd Y " Y0
LOS
S•.-L-
w Yd Ys 5
ASr Yd A6s
r Ys '*Yd Y's
where S and LY are the relative errors in Sr and Ys, respective-
ly. Figure 9 presents this relation in graphical form.
.,*' /
.a.
For the most tightly packed materials, the relative error in de-
termining the degree of saturation can thus be 3 to 4 times greater
than the relative uncertainty in specific weight of the mineral
particles. If the latter uncertainty is assumed to be about _5
percent (corresponding to Y. 2700 + 150) the error in Sr can thus
amount to 15 - 20 percent in tightly packed gravel materials. The
corresponding error in determining the relative particle volume is
+5 percent i.e. the same as the specific weight uncertainty.
106
Higher accuracy in determining volume ratios in soil materials
requires knowledge about the specific weight of the mineral par-
ticles. This is normally obtained by means of a pyknometer (14).
,25 I
> 1.
1,5 I1,9
107
From curves of this type one can determine such parameters as the
size distribution coefficient Cu a d 6 0 /do 1) or the sorting number
S d7 5 /d 2 5 , which together with the median Md a d50 can be used
to obtain a numerical approximation to the particle size distribu-
tion or to represent this distribution as one point in a diagram
(17). This can also be achieved by separation into three fractions
and plotting the result in a triangular diagram as indicated in
Figure 8.
1) The subscripts for d indicate the percentage of particles which
are smaller than the diameter in question. Thus, 60 percent
of the particles have a diameterid 6 0 , etc. (Translator's note)
108
$AND
40,1 - -. q-.--j , . 4 4 ; . i. j
. 1.
I s
. .. ,......
...... ,," ..... . ', ,•
**, A -- u II4,*~ .,w
109
In connection with thermal properties of soil materials it may
seem as if the particle size distribution curve plays a secondary
role as "parameter" for the more central factor of specific sur-
face area. For a certain particle shape (e.g. spherical), the
ratio between surface area and volume (specific surface area) is
inversely proportional to the particle size.
1...*. .-.. . ........
....
"S f/d 7
Discs 1 (2 d/t
Cylinders sa'
110
while a logarithmic distribution gives
f1/
-ev di
UI 1
5%
111
"ini n li i" - .... .
The logarithmic distribution will clearly always give lower
equivalent diameters than the linear distribution. For wide
size variations one also finds that the equivalent diameter is
close to the minimum diameter dI,
At any rate, these results indicate that the specific surface area
of a soil material is largely detemined by the content of finer
fractions, with the clay content being the most important parameter.
Within the past few years, a method for determining specific sur-
face area has been developed based on adsorption of a monomolecu-
lar layer of ethylene glycol (19) or ethylene glycol and mono-
ethyl ether (20). These methods have been used by Anderson and
Tics (1972) in connection with a study of un-frozen water in soil
materials (21). The results from this study will be discussed
in the next section.
112 I
C. Un-frozen water
The water in afine grain soil material with a water content be-
low full saturation will freeze at a temperature below the normal
freezing point for water. The decrease in freezing point is often
related to the moisture potential in the soil material at the ac-
tual degree of saturation. For example, Williams (19?3)(22)
gives the following relation
vT
6T * Eo12
&T is the decrease in freezing point (°C)
To is absolute temperature at the normal freezing point (K)
L is the latent heat of water (J/kg)
v is the specific volume of water (m3 /kg)
2
Sis the moisture potential (N/m )
AT a 0.0815 * (bar) 13
Low, Anderson and Hoekstra (1968) point out that this equation is
valid only for soil materials with relatively high degrees of
saturation (23). Using thermodynamic theory they formulated a
more general relation between freezing point decrease, moisture
potential and latent heat. This complicated expression for AT
was also used to generate an extensive table from which either
of the parameters mentioned can be found when the other two are
known.
113.3
mine the content of un-frozen water as function of temperature when
performing freezing experiments on samples with different water
content (degree of saturation). For relatively low water contents,
this method can be expected to give low accuracy. For this reason
calorimetric methods are mostly used to-day.
wu ( AT )8 (percent by weight) 14
In a 0.552 in S + 0.262 15
in (- 8) a 0.264 in S * 0.371 16
I1
114
weighing 60 g to sink 10 mm (sometimes 5 mm) into the sample
(v 10 , v 5 ). Selmer Olsen (1954) show that the water content de-
fined in this manner agrees very well with previous definitions
of the flow limit (26).
w100 - 3.72 18
u2 R 0.338
where wul and wu2 are the un-frozen water contents at -1 and -2'C,
respectively.
By plotting wu and wU2 in a log - log diagram for each of the two
temperatures one can determine the contours for un-frozen water
content by drawing straight lines between these values, due to the
exponential nature of the equation. This results in a graphical
method for determining the un-frozen water content at various
temperatures, as shown in Figure 13.
11 5
L_= . . .. .. . ........
i 0c I I II . ..
800
•- ,, .. 2. .-•. • _•
42------------ _-
In many cases only the flow limit w2 5 (V1 o) :ls determined. However,
Selmer Olsen's consistence limit determinations for a number of
Norweg~an clay: show an avorage ratio of about 1.2 between ViO ad
i
Consequently, if
only one of the consistence limits is given, the
other can be found from the relations
v 10 1.2 v 5 20
or w2 5 1•2 w1 0 0 21
This forms the basis for the thin un-frozen water contours in
Figure 13. If both consistence limits are given, the two scales
may be used as described previously.
80
I ,
1(.0
W ,W,
>"20 j"
- VO . VS
0 ?0 40 60 1(
V S W'10' %
11"7
Until determinations of un-frozen water content in Norwegian soil
materials become available, the approach taken by Anderson and
Tice may be applied also to Norwegian conditions. Figure 14 also
indicates that the normal range of flow limits for Norwegian clays
is well covered by the materials studied by Anderson and Tics.
A. Geotechnical analysis
118
normally uses a sampling tube (diame.ter 54 mm) to extract "lndis-
turbed" samples. The samples are taken for esch metre (a metre
at a time) to a depth which depends in soil conditionsý. From 4;hese
samples one determines volumetric conditions (density, water con-
tent), consistence limits, screrning curvel) and various other
properties of importance to stdbility and compreasibilit'.
, I119
T _
... , .. . ...
• .... ;' '
0! *r d . f
01,J
.- • .s
- aj
U. . ,. i ai.•
" *a*•
:,;. .- . . .
V-.. 4 ,
- V!2 . U ;12 0 ner'- o
• ra, :•a r
¶ a *
:' a. J.• .
Proiline 5k i •sc av 3p r.'r. '., av e' n v e .,oV •' r bY • •• Upe-e1 . a n r iv surfac
ze . ~5 -'-
Of .z :O,2, D'J"p:,~.j'I .
,Jott i '
1-W i c
]--. -'rr I.,,9'P ri tv'Ik" 1 c./]..
"k'i.o.
L .l. P (LIi LradI
W OW"
r' ''~~ I , L
Dv Lil~ll,
ýa*'Ll t m "', ,yi.•y 'J P
4,n
I -),
,T,/'r', r,)1tVL,'k soGm rjf J i•U+mti
ay vuxnninmnnoý. yu-1+r la,~Tn.,
9
,,,.. .... ... T ." j .........
• -....... ,
.4 ,P • •_ _ • _ . ,
1211
I~1.3. 3.1. KI'I,iV '.UJ Io~ (Or,ifotJ I ,¶ItA PI..r J~t I'r d;5l 1..' ivrl+.+
ia|,JUr.,
121.
- .- ,- -. I
Requirements in wear resistance for the materials are checked
through fall tests (crush tests) and passing through a rod screen.
A petrographic determination of the rock fraction is also included
(30). However, this evaluation is primarily of a qualitative na-
ture. Artificially crushed rock materials are subjected to a more
quantative evaluation, including the percentage of different mine-
rals.
Type of analysis
-'~~
94v.r-4~
10)r0
1) Given by proctor
122
-.- I--,-;y*-~.'i
B. Supplemental analysis
1VARfl4_SD'4: NGSEVNE ,
{ WE4IWIA.,%. '• 'C,,t
0• SO,'.
T.'S ,VA NN
".'JF !PA .
It'CLE COI.NVO AC','N, ,EY ,'AF
I.'r
'775
T'S!vSR
'VTs!'I4N": .7 E
. :S S ..... i Q GANIS•
I !r~~r;
,"IS
LYTr:,3,1ENSr7,L~eu
c~
123
Together, T$. and Figure 18 give an overview of how compre-
hensive the • normally performed during road construction
are for the factors which determine thermal conductivity. One
finds that the tests listed for fine grain sub-soils are almost
* equivalent to level II, except for knowledge of quartz content
in the material. However, this parameter can be estimated from
the particle size distribution (screening curve) with sufficient
accuracy for this case. For coarser (sandy) sub-soils, informa-
tion about dry density is lacking. This makes results of common
analysis incomplete for all the levels shown in Figure 18. For
materials used in the upper levels, informption on water content
will be lacking. One must therefore make an estimate based on
experience from similar road bed designs and take into account
expected variations. The petrographic analysis, which is limited
to rock fractions and also is rather superficial, gives little
information about quartz content in the material at hand. For
"* that case, the sample analysis may thus at best give sufficient
information for the lowest level of conductivity determination.
124
REFERENCES - CHAPTER III
125
........................................... I
16. R. Selmer-Olsen: Ingeni~rgeologi. Del I. Generell geologi.
Tapir forlag. Tronheim 1971, p. 8.
126
CHAPTER IV
EXPERIMENTAL METHODS FOR DETERMINING THERMAL
CONDUCTIVITY IN SOIL MATERIALS
* q- VT
q is the thermal flux (l/ma)
127
aT L V2 T2
. 2 3
For variable conductivity:
*T v( 7 T) 4
Scp
For the case where no water transport occurs during freezing, the
thermal conduction equation for the partially frozen zone in a soil
material can be written:
128
4
( CO + u YL) BT
-v(X T)
ST
The latent heat when freezing can thus be included in the apparent
thermal capacity/volume, but this parameter will vary extremely
fast with temperature near the freezing point. An example is shown
in Figure 1. The curves apply for a silt-clay material having a
dry density of 14090 kg/m3 and a flow limit of 30 percent by weight
(see Figure 13 of chapter III).
I.
" .t , j j
TEPERA%~ I
129
expected. The approach will then be to determine an average ther-
mal conductivity, based on heat flux and temperature difference.
This average conductivity can then be related to the average tem-
perature in the samples. Non-linearities in the conductivity due
to varying amounts of un-frozen water can be expected to give cer-
tain errors in this relation. If the temperature coefficient for
the conductivity can be assumed constant within the given tempe-
rature interval, the errors should be insignificant. However, this
will be true only if relatively small temperature differences are
used across the samples.
130
One must clearly impose certain requirements on methods used for
measuring thermal conductivity within the region of critical sa-
turation. The goal must be to maintain changes in the moisture
distribution within certain acceptable limits. This could be
accomplished by lowering the temperature gradient or by reducing
the time during which samples are subjected to this gradient, or
by combining these approaches.
2. STATIONARY METIHODS
131
&VRE SKJERM • '
, -. . .,.
SENI. ALT i
vARFE,=LEMENT "•.•, l
NEU.'..O !, ,..,
SPRI~T --.
SPR!Ti
I nim*!.n
Xr 1 j•.. P.i
PLUGO FMI
-* Pt~~rrnon &p~r~
anip or r av o:yi.r var*mn"-
qOZ
1500w'" 50i Diameters 150 mm, 50 mm
SKJEM SEKSJON
•VRE
.I Cro,"ecton. upper screen
JORDPROVE
SENTRALT VARMEELEMENT
MALESEKSJON
Soil sample
:Center heating
S.. element
Measurement zone
. ...
NEDKJ0LT SPRIT 2Cooled alcohol
NEDRE SKJZRM SEKSJOX Crou..sectiong lower screen
SPRIT INTAK Alcohol inlet
BORD TBoard (table)
PLUGO FOR TOMNING :Drain valve
132
latter correspond to the ring-shaped screens used in planar arrange-
ments )4 One tries to minimize the temperature differences across
the gaps between sections by adjusting the (heating) power in the
screen sections. These adjustments are made manually.
Kersten does not report typical stabilization times, but his des-
cription of the procedure used at the end of each experiment gives
an indication of the conditions.
Temperatures at the inner and outer surfaces of the sample were re-
corded each time, as well as current and voltage for the center hea-
ting element. If the thermal conductivity calculated from these
values was found to be constant within error limits of one percent
during three to five hours, the final conductivity was computed as
the mean of results obtained for the final hours.
The long time required for each experiment is perhaps the biggest
drawback with the stationary methods. This is also the reason why
non-stationary (transient) methods (particularly thermal conduction
probes) have come into use, despite the considerably higher accu-
racy obtained with stationary methods. However, the long duration
133
of each experiment may also be a questionable factor due to the
previously mentioned moisture transport.
In r 2 /r 1 - in r 2 /r 1 - 1
/1
Sr
2 /rI is the ratio between outer and inner radii
134
Under the assumptions made, the conductivity obtained for even mois-
ture distribution is given by
Xm = 1 ý2 +7
0.5. -
0.2 -
1X x,
2.5,
2.0-• //
IX Y
z-, 1.5 ,-
"-~.2- 05- I /Y
S0.5 .
0 0 A- I I
I
0 0.5 %.0 I.5 1.0 2.5
F 6RSTE M;\L'NG. Wire
m
It can be seen, that the final (repeated) test with few exceptions
gave a lower value for the conductivity than the first measurement
at 20 0 C. However, the deviation is in the majority of cases less
than 5 percent.
Experiments with moist soil in a planar (plate) test unit show con-
_ , 60
"A SA !-'LCW
. ' ,6 f
"•-'-• 1 ? 2. S ~'0 o2 • •0 'O ?
T'9 ! n':!,r'
One may also note the large thermal flux at the warm surface of the
sample during the early phase of the experiment. This indicates
that large temperature gradients occur at the warm sample surface
following'a step-wise increase in heater temperature at the start
of the experiment. Correspondingly large gradients will also
occur initially when experiments are performed in Kersten's circu-
lar (cylindrical) test unit. These large gradients will obviously
accelerate the moisture transport.
Best Available Cpy
One could possibly reduce these gradients by using a "softer" ini-
tial temperature transient (step). However, this would lengthen
the stabilization period and allow the moisture transport more time
to re-distribute the moisture content.
I:.C = T )
138
program named "Frost in the Ground", directed by NTNF and the
Highway Authoroty with participation of the "Research Group for
Thermal Analysis of Frost in the Ground" at the Institute for Cold
Technology.
Figure 6 also shows that the heating plate is divided into a mea-
surement area and a screen area. The power applied to the screen
is automatically controlled by means of thermal sensors placed
across the gap between measurement and screen areas so that the
temperature difference across that gap is maintained at zero.
139
A. Cross-section of the heating plate
unit.
B, Top view of heating plate unit.
a. Central heating plate with hea-
ting foil 0.2 X0.2 m2
- -
.. .b, Ring-shaped screen with heating
SL ',foil
*, .
,- - \- \.. . . . .. r i i
. *' *i
~14
C. Assembly sketch of plate unit.
h. Heating plate
±. Cold plate with pipe. for circulating alcohol and heating foil
j. Gradient screen. Heating cable at top.
k. Upper cold plate.
1. Mixing tank with magnetic valve and pump.
The power applied to this element is also controlled by thermal sen-
sors so that its temperature is maintained equal to that of the
central measurement area. The stainless steel sample holder con-
nects warm and cold plates thermally such that one obtains approxi-
mately the same temperature profile in the connecting member as in
the sample material (gradient screen). This prevents heat leakage
from the sides of the samples.
- -.
141.
1
of temperature stability that can be obtained with this test equip.
ment.
F let .. -. .. ..... . P t a
'F'!
- -r
7, Ot 20.7.
-Z 4
-o .- 0 7'.
' ,2
2 40
hV
1.42
In addition to the plate units described (3 units), equipment was
developed and built for transient measurements by means of thermal
conduction probes. This will be described in the following sec-
tion.
143
T - To q rEi r' )I~
4 at 9
- El (-x) * - y - in x + x - x2 + .... 10
For 4 at/r 2 >1 one can thus write the equation describing the
temperature field in the approximate form
S~q
T- To -Y + lnt + In 4 a/r 2 ) 11
145
Kierkus et.al. (1973) (19) have performed an analysis for a thin
filament of finite length, with some-what different boundary con-
ditions. The study was performed in conjunction with conductivity
measurements in liquids, when the probe consists of a thin metal
heating wire. Since the conductors required for supplying the
heating current have significantly larger cross-sections (than the
heater wire itself), the temperatures at the ends of the probe are
assumed constant during the entire experiment and equal to the ini-
tial temperatures* Using modern computer methods for numerical in-
tegration, the temperature rise in such heating wires of finite
extent was computed for different length/radius ratios and diffe-
rences in thermal properties between heating wire and the surroun-
ding medium. The results obtained (from these computations) give
a good illustration to the dependence of different parameters.
Some of these results are shown in Figure 9, where temperature rise
and time are given as normalized1) numbers.
These results show that an infinite long thinZ) cylinder has the
same response as an infinitely thin filament when the cylinder has
the same heat capacity as the surrounding medium. If the infinitely
long cylinder and the medium have different heat capacities, the
temperature rise will initially be steeper than for a filament but
after a certain time the response will flatten out and coincide
with that of a filament.
146
- Lit
Lit
z_ LS- LI .r,00 /
;" | X~),II, C
;',2
fl, '0i
02 "
CII,,0
"00
1
/
0.7sF
2000
0.2Si
10
L V0 a 0
e tsi t0e
0n 1en0.,Ie- 'o,Ü r". tt
'~c /oj.
VariLd.th.,dusr'n•.
eer~.~
i oi
.ýn C~rrr1o'.I
give .
:ion vartob.:c:;pac:'
co,/ ..?7I
n
The finite extent (length) of the probe will, after a certain time,
cause another bend in the temperature rise curve towards values be-
low those for the infinitely long probe. This bend will occur ear-
lier for lower values of C/C 1 . Thus, for different heat capacities
and finite length one may at best expect that the probe will have
the same response as a thin filament over a limited time interval,
148
Comparative investigations initiated by the Institute for Cold
Technology indicate that probes having the previously given di-
mensions (d = 2 mm, L = 20 cm) give up to 20 percent higher con-
ductivity values than steady state measurements performed in a
planar test unit. This applies both for dry and moist soil samples.
However, experiments performed with longer probes (d a 1 mm,
L = 100 cm) gave much better agreement with results obtained in the
planar test unit. Experiments using "probes" consisting of a hea-
ting wire (alone) also gave results in good agreement with planar
measurements in dry and saturated soils. Even if these studie6
are still going on, they still serve to support the analytical re-
sults discussed here, including the new requirements for length/radius
ratios developed here (20).
8T 14
ar 21r X
149
~
. *"A -,
If one wants to use a power which causes the temperature to rise
by IOC over a 10:1 time interval (e.g. from t = 1 min to t 2 = 10
min) one finds from Eq 13 that the ratio q/(2Hx ) should be very
close to 1. During the measurement period, the temperature gra-
dient near the probe then becomes:
3T 1
- = __ 15
ar r
The gradients which occur are significantly larger than those oc-
curring during stationary experiments in a planar test unit. It is
doubtful whether the much shorter measurement times (using the probe
method) can compensate for this difference.
I, 150
is no similar effect and the effect is also relatively small for
the. highest degree of saturation.
,- ". 20 C it
b I
16 cm rvtaC/.
C 0~c
0 I ! ,t ~ II IO 2
20
20 3oo00o
T
Tl0 PTTER STAR . SEKUNDER
151
C. Experimental arrangement
152
15 kg. Preparation of six samples normally requires one man-day.
The next day, following the thermal conductivity measurements,
tests are made of the moisture content in the containers.
1I5
:53
r14- ITr
A- r7rI I \I i <
IIr
IC 2?7ým i,,t
Fit J,
4
ýo~I ~ t Ag eorvaih pd a'niriin
i-r for '~tN~.; I. i
154
Legend, u Ii.
A -Fan
B - Cooling olement
C - Oven
D - Pt element (Thermal sensor)
E - Probe
F - Material
i•XPERIMDT :o 5
PROBE No 4
INITIAL PROBE T*1,PERATURZ : 4.29"C
155
I.. . . ....... ............ ......................
.......................................
*444*4
fl. ~ ~ .. 444
....... . *4...14
.. ......... **4**4*4
F44181. list
.1. 141 1
k I.
!0 t ,dt I i
hall ..*!t3 a t, v is t
156
4. OTT-I" TRAUISIE•T METHODS
The various method which have been utilized are based on either
a programmed temperature variation on the cylindrical surfaces or
constant thermal flux. The latter case is most aprlicable for a
hollow cylinder when hoat is applied from its inside. As will be
evident in the following, thermal diffusivity can be derived from
data on temperature only, while computation of thermal conductivity
and heat capacity by volume always required knowledCe about both
thermal flu:: and temperature distribution.
157
If the temperature rise on the surface is given by
T (R, t) - To = Bt 16
T (r,t) - To - -8 (R 2 - r 2) 17
4a
a-. BR2
44T
TSR
- )R - 19
3r
2a
158 5.
"-A
2C - 'C ,._
~~T1
. ' ... 2TC
Soo 100
Og0.5? 2CO 0
T!20, IS-K.i NT.F.
~'I
IIa'JrI'OIU•
[, 'i.al,
e.: I.rILt!. C p ip(rrI.:. -V
I- !'1 ( t .U1 P I P.&I
I.4~
in
".Wr~j r I.- it, )I if r' it C
For example, when considering a fine grain soil material within the
previously discussed critical range of saturation levels, the ther-
mal diffusivity can be estimated to 2-10-7 m2 /s.
FO M at/R2 - 0.2 20
159
Fo is Fourier's number
6C/sek
e. r I4000
sek
< 5/ 3000
Cr
2000 a. zT Io
cc L 4U: O
2 ,000
J -C 0 1 2 3 4 5 6 7 8!0.2 m
-u,,,'e u
wit r: o .j;'.5j. j;:. a = 2 I(-/
' ,
!oLIfl .. ,..; ., 3 ?,C 2J. IL'*J'S.: (
One finds that the smallest radius which satisfies both conditions
is 2 cm. This leads to a quasi-stationary period of 400 seconds.
The rate of temperature rise which satisfies the requirements is
2 • 10-30C/sec, corresponding to about 7 0 C/hour.
This effect means that this method can be used for soil samples
taken with the 5-.4 mm diameter sampler normally used for geotech-
nical studies. Some trimming around the periphery -mill probably
be advantageous, but this is common practice in most normal studies
of samples of this kind.
When the thermal flux is knovrm one can determine both thermal capa-
city and conductivity for the sample material, as mell as diffusi-
vity. Aguirre-Puente give the following expressions for c0 and
co= q 21
O - q
22
411A T
161
concentric cylinders. The heating elements are attached to the
inside of the hollow cylinder. The outer cylindrical surface faces
a heat screen (guard) which is controlled by a signal f~rom a tempe-
rature sensor at the outside of the sample. In this manner, adiaba-
tic conditions (no exchange of heat) are ensured on the surface of
the samples. A sketch illustrating this method is shown in Figure 15.
II ,
, , SENTRAL VA,•M(,UOE
7/ -'SCN'Ut. HEAT SOUICE
ra * CONSTANT
AAT
/ , A IA-A11C GUARD
j 162
The method utilizes the principle of constant heat flux. However,
also in this case a quasi-stationary state will be reached after
some time where the temperature rises linearly (with time) and a
constant temperature difference AT exists between inner and outer
surfaces of the sample. In this phase, the temperatures are given
by 2
q R2 2Fo - 1 r2
T (r,t) - To = R 1 (1- 2
R - R 4 R2
22 33 2
R2 r R2
- 2 + 1 R 3
(in 2 2 i + )j
2
R R2 R1 R2 R1 4
The thermal parameters for the sample can be computed from the rate
of temperature rise a and the temperature difference 4T:
XR 2 R2 R2 2R2 I R
R2
2 (R 2 - k2 ) AT
1 2
2 2 2
a at (R2 R -2R 1I in
4 AT a2 25
) cP 8 q R2.. .
2 2
26
163
The temperature rise within the samples is determined by taking the
derivative of Eq 23 and setting r RI.
aT
-qbX
Dr R 227
where q is the thermal flux per unit area of the inner surface.
If the heat flux is measured per unit length of the cylinder (like
for the thermal conduction probe) one obtains
R 28
ar 2n R X
1
• 00
20 1.
I , a
164
As mentioned previously, this method has the advantage that the
thermal flux can be measured directly as current and voltage for
the heating elements. An axial heat screen similar to that used
in Kersten's cylindrical test unit should contribute to increased
accuracy in these measurements and ensure radial heat flow. The
fact that temperature sensors can be mounted directly on the walls
of the sample holder also represents an advantaged in comparison
to the previously described method where one of the temperature
sensors had to be centered within the sample itself. However, the
geometry is somewhat less favourable than the small (compact) cyl-
inders in terms of sample preparation. In addition, undisturbed
samples can not be used.
165
IkiNt -LA
_T. Av.dy1
qjtl-550
n?r(v1rA r
IIl'Pi i'a
. 6C ,ma
It m.*1 r? r0 01
Fhaseshift tan bei X 29
bor X
X -R wa 31
bei and ber are Kelvin functions, which are tabulated in mathe-
matical handbooks (e.g. (31), page 430).
-• .- O '4INGA J' *
- = - *4*~ ~or
Wt- V1NX1:E1.AS!!GHC!T raid./%
3 f IP•~ P~(Y'fA'NE
iGIl'F,'AA: DIF•'US WIVI
P1R
.If.RAOIUS
v
1AM•I.t$L
RADIUS
I 2 3 . 5 a '7 0
167
Such a diagram can be used for determining × from measured phase-
shift or amplitude attenuation. Thermal diffus.vity can then be
calculated from Eq 32, if period (frequency) and sample radius are
known. In cases where the attenuation is strcn:,, calculations
based on phaseshift will give the best results, as evident from
Figure 18.
2000
000, L ,,t
.1WWI
c ' I
The sample radius can evidently not be much larger than about 2.5 cm
if one wants the period to be shorter than 10 minutes.
A prototype unit of this kind was built at the Institute for Cold
Technology in the spring of 1972, as part of the study reported here.
Figure 20 shows a sketch of the set-up, as well as a modified ver-
1.68
-- - - -- -- - --
C,@*'0 04 4y
00*1
of *0kIo1-0 4 "f
II0,'
.. ..... ...
porlatilfl autoati
regiuin a tm~lqltre
ved tslakeilkuomatiSng
sgg~ectd mesiiod foi' auto'-
matio camYpalratdre Oontrol
b'ao~d om feed-bacOk.
'
169
sion with improved temperature control system. The sample holder
was a steel tube with inner diameter 60 mm, wall thickness 2 mm and
a length of 150 mm. A heating foil of type ESVA1) was wrapped a-
round the tube. To ensure constant average temperature in the sam-
ples, a 1/4 inch copper pipe was also wrapped around the outside,
separated from the sample holder by a layer of pipe insulation. In
addition, the sample holder was insulated from the environment by
means of pipe Insulation on the outside of the cooling pipe.
~qgt nt
of the power level is required (via the auto-transformer) in order
to obtain a stable mean temperature in the samples. This balance
is reached when the power in the heating foil equals the heat loss
due to the cooling pipe. A better solution has been proposed, by
which the signal from the function generator is used to control the
surface temperature via a feed-back loop containing a proportional
integrating regulator, as well as the previously mentioned arrange-
ment with a thyristor control and auto-transformer. This will permit
automatic setting of the mean temperature by adjusting the signal
source (function generator).
170
then correspond to the equipment developed by Aguirre-Puente. As
mentioned previously, such a test unit will yield information about
thermal diffusivity, thermal conductivity and thermal capacity in
the samples. This has also been achieved by Champoussin (1972) in
a planar test unit using periodic temperature variations (27). The
cylindrical configuration using sample holders made from a material
with known thermal properties can also be expected to yield this
(same) information.
171
is suitable for measuring ti:erzmal Qiffusiv_.ty of fino grain mate-
rials within the freezing zone, despite the non-linearities which
occur. This is based in the possibility of using extremely small
temperature variations for the tests. For example, it is stated
that amplitudes as low as 0.20C can be used without loss of accu-
racy. (Thermistors were used as temperature sensors).
a (M(T) aT ) u co(T) aT 33
ax ax rt
OX aT +X ~2r' cp MT aT
The fact that the thermal conductivity and capacitance are tempera-
ture dependent will cause deviations from the linear equation. How-
ever, if the variations due to temperature are sufficiently small
one can possibly neglect also these effects. In this context, varia-
tions in the apparent heat capacitance will be most important.
Since this parameter may vary by more than a factor of 10 for a tem-
perature variation of only 0.5'C near the freezing point (0 to -1 0 C)
it would seem that this condition poses unrealistic requirements on
(sufficiently) small temperature variations.
172
A more complete mathematical analysis of the effects of these non-
linearities may give a better understanding (clarification) of
these matters.
D. Conclusion
173
- ~ 4 _j
4 REFERENCES - CHAPTER IV
174
16. D.A. de Vries, A.J. Peck: On the Cylindrical Probe Method of
Measuring Thermal Conductivity. Aust. J. Physics, 11, 1958,
pp. 255-271.
17. J.H. Blackwell: A Transient Heat Flow Method for Determining
the Thermal Constants. J. Appl. Phys. 25, (2), 1954, pp. 137-144.
18. J.H. Blackwell: The Axial-Flow Error in the Thermal Conductivity
Probe. Can. J. Phys. 34, 1956, pp. 412-417.
19. W.T. Kierkus et al: Radial-Axial Transient Heat Conduction in
a Region Bounded Internally by a Circular Cylinder of Finite
Length.
Dept. of Mechanical Engineering University of Calgary, Calgary,
Alberta (Preprint).
20. P.E. Privik: Institutt for kjoleteknikk, NTH 1974. Upubliserte
resultater.
21. W. Woodside: Probe for Thermal Conductivity Measurement of
Dry and Moist Materials. Heating, Piping & Air Conditioning,
Sept. 1958, pp. 163-170.
22. P.E. Privik et al: Sondeapparat. Rapport fra Gruppen for
termisk analyse av Frost i Jord. Institutt for kj~leteknikk,
Trondheim, NTH, 1975.
23. 0. Krischer, H. Esdorn: Einfaches Kurzzeitverfarhren zur
gleichzeitigen Bestimmung der Warmeleitzahl, der Warmekapazitat
und der Warmeeindringzahl fester Stoffe, VDI-Forschungsheft.
450, 1955, pp. 28-39.
175
30. L. van Zee, C.L. Babcock:A Method for the Measurement of
Thermal Diffusivity of Molten Glass. J. Am. Ceram. Soc. 34
(8) 1951.
31. M. Ambramowitz, I* Stegum: Ed. Handbook of Mathematical
Functions. Nat. Bur. of Standards. Appl. Math. Ser.,
New York 1969, p. 374.
32. P. Hoekstra et al: Measuring the Thermal Properties of
Cylindrical Specimens by use of Sinusoidal Temperature Waves.
Cold Regions Research and Engineering Laboratory, Technical
Report 244, Oct. 1973.
J~J
176
-A
CHAPT2R V
i. -. ,
mental results pertaining to Norwegian soil materials is now being
improved on by a series of experiments initiated at the Institute
for Cold Technology, NTH. This investigation is part of the acti-
vities within the committe for frost on the ground. Some of the
results will be discussed in the following section, which is devo-
ted to measurements on dry soil materials.
A. Introduction
•,' 50 r
'I0
20
0
0 0.2 0,4 0.5 0.• •1.0
PART IIXKELV0LUW ANOEL
VOb.diMgF*ACTIO OiP PADICIII•|
178l.
:,~~':•• O , ' tl 0 •••:r~• :L.U"L.
For a relative particle volume 1-n = 0.60, the conductivity at the
upper limit is at least 7 times that at the lower limit. For lar-
ger conductivity ratios, the distance (spread) between limits will
be even larger (wider). As mentioned earlier, these limits repre-
sent extreme differences in microstructure. Even if such differen-
ces are not likely to occur in soil materials, the wide spread be-
tween these limits still indicates a high sensitivity to differen-
ces in the microstructure for dry soil materials.
'A,
'001
GC
VWiC!V'TY "ANO
IVIO N i
a-lk biv ll U.14 I twa ~ v r
PIz O, j# , (•,,"%I
•k a,
l . IA. , to , I ":
'! . , '
179
i.
creases by 50 percent. By comparison, when the conductivity ratio
is 10, the same relative increase in conductivity requires leas than
a three-fold increase in particle conductivity,
04
I S
0.2
•.+.0 * ' 4*
O e Slf*lM 4 *
0.1
0.2 0.3 0. 03.
. of 0.1 10
18,
180
j
As can be seen, these results show a variation of about t20 per-
cent around a mean curve which can be drawn as a straight line in
the log/log diagram. Variations among the results can be ascrib-
ed to differences in mineralogic composition as well as different
textures (microstructure).
Figure 8 shows results for four natv-ral soil materialn, all having
quartz contents in the range 4i to 53 percent. For those four
materials one can expect the particle conductivities to be about
the same, about 4.0 W/mK. For each material, the results follow
a straight line in the log/loi. diagram rather closely. Variations
from material to iaterial are .i.rnitficantly larger than the varia-
181
i i i , i i.i
'1. 1V,r I 7
I d 10
ow nvr .I~
0
I ":ht7
?r720£roo 0'4o
TABLE 11
192 02ý
Translation of terms:
JORDART Soil material
KVARTSINNH. 0, Quartz content %
SPES. VEI,:T kg/m3 Spec. weight kg/m3
TORR R. VEKT Dry density
POROSITET Porosity
LEDN. EVNE \V/mK Conductivity W/mI,
Tyholt leire TyholtI) clay
(pulverisert) (pulverized)
Folio silt Folio )silt
pulveriser t pulverized
Klumpet Lumps
Vormsun silt Vormsund ) silt
Tyholt fincand Tyholt )1fine sand
Jessheim JessheimI
finsand fine sand
Skanir Skan~rI)
finsand fine sand
Boskorp BoskorpI )
finsand fine sand
Tiller sand Tiller1 )sand
Tyholt TyholtI)
fin grus fine gravel
0Otfold ostfoldI)
TrUS gravel
Troms grus Troms I)gravel
_ __ _ _ _ _ '
183I
TABLE I. B. CRUSHED ROCK MATERIALS.
MATERIALE Material
KVARTSINNH. % quartz content %
SPEC. VEKT kg/m 3 Spec. weight %
TORR R. VEKT Dry density
POROSITET Porosity
LEDN. EVNE W/an Conductivity W/mK
Kunst kvarts Crushed quartz
Knust glimmerskifer Crushed mica shale
1
Crushed greenstone )
Xnust gronnstein
Oppland pukk Oppland macadamZ)
0atfold pukik ,0'stfold macadam2 )
ai07
,4 ii.r ok I.for50 17'0 *,,,,
179C 0,.I,'
,. l j ,!0 ,, •" 0.•9
A,0-Ito n I 0 e 3
PIa
. 6. • O .s60 )2•
Io ,31K33)
KN 312680 1"920 1( 1'21
1 0 ,.,.1.
190C 3,:i1
184
r140
01` W0411 aa to to h a 1I 1,4" 1 1"
03
I- Z
451410 Al.4)81
Fin
. f ikt~exury-ir Cor U-4s f Inkornsoivho J 2:,. Vr~a
UP1f n 1mL~
Ulf .1 8 1 ~i011 L.*
iona 1." J Ll
94hro j jpýA4."W'.
ff-M
)11 lvls ia e1l01,1 ___p~~ j&0 __1J ' dJ4
V1, O -f-l'-
UM10 4/jou to 30 I V 1 . . 1 1,
i -' -l T ** T - 8-
CY.
ACA
sO00,110 Al 4101
195
ý014 IN
*di.. too ~
%.1... 3%)o 1A
L; q j
go 0
.. I
ne. Ovate ,Agpl4 jisel1e k4104,P*4
*300 ~r~rr7!
~ 4
14I
toI.
`7t~
0-. IOU
0.001, e w '0.149 .4 C4 1).. 441)
! L0AN
.0IJ80l KI'.PW
, ,ýN T fu
3
7 .¶iktakurvur for- remn knustfv bur~art or tri~ nannip ng. uAr-
* 1 ke.06 aIO
;.J J8Ift Af4djo jcpjb1 g ,,I IIIe
186
4 ~0 5'
0
0.2
0.0
S 0.1" 0. 3 0.*l
.&!.enint-sevnen
Iz. av t~rr* j0:IILrter OK knust# btfrg:.r'tel,
p~eI.I.CI Sam~ Vi'nKSjonf ay pr-.1111tuIrl. (.'&J.14cIt C v t t a
of dr s fI 4
Puti 0,i1-n~h14 'O1# ...
ro : nrit4 1: JnDto I~ r
*4aturJige joJar4t4er, &Knustbo0FI!o
NlaturaiJ uoý..C in'
1-0~iid
S% 14VARTS
1 *.3
0.2~I0.
K ~ ~ ~ ~~1.87-'-....--
tions around these straight lines for each material.
188
As shown, the experiments give no basis for expecting variations
within each of these groups due to differences in microstructure.
D. Empirical relations
n - 1.0
X X'2 n - 0.0
*n + (1-n) F 2X
n + (1-n) F
1* + g2 + 93 1.0
189
iA
Values of the factors g1' g. and g 3 formally refer to definite
geometric relations for ellipsoid-shaped particles. However,
as mentioned in chapter II, the equation is based on approxima-
tions which make it impossible to derive these geometry related
parameters from the particle shapes in an actual material. Thus,
g1' g. and g 3 should in practice be regarded as parameters avail-
able for fitting to empirical data.
The latter set of g-values are thus more representative for the
conductivity of natural soils. If the conductivity ratio is set to
125, one obtains the following semiempirical expression for the
conductivity of dry, natural soil materials*
n + (1-n) 6.65
/0.024 -n + (1-n) .053 20%
190
A0 -
XA -200
0.35
y + 6_.7 20% S
I/
20.04
.
U.'!I
/ / /'K,
0.1. r*
l*-O2; / ," / 1. I'
0 '
F.
IQI ? 3•.3, ui n iknirlij t e .. UPI ! ± rLnin l.uVru .i ,l
One finds that this equation gives a relatively good fit to the
,• experimental conductivity data for crushed dry materials, if a
conductivity ratio of 125 is selected. However, due to the
complicated structure of the equation it is more convenient to
use a purely empirical relation for the mean curve. The follow-
ing exponential function is found to give a good approximation to
the conductivity in dry crushed rocks, within the porosity range
considered
192
However, not consistent (with all
this equation is parameter val-
ues).. Thus, it should be regarded as an approximation to Br-ugge-
man's equation (Eq 5) which meets the physical boundary condi-
tions mentioned.
W/m10 K --
0.4 'l
120
0.1
*'P•
13., lI' • .mai•(l11fl,,-,',;, t'u, I ury11 r'. vriP l~arpil~ 1 Ii.'li
,tU•
I'. Ill. IIt ,'C'l' J
•a ".I l~t 'll. iii J llfl~'ltl. I M jllllli ~ ni
E. Stab.lilzed materilals
van Rooyern and Wi-nterkorn (1959) (3) and later Farouki. (1966) (11)
have published (results from) conductivity measurements on dried
soil materials. For some materials, these results show much
higher conductivities than even the "highest values obtained at the
S.........
.. . '.. . ..... .. . ,, 0. ..--0 -- '- S-- 0 -. O .- :i',s-
1 ,~.; -0, -- , -•-: • • .•
Institute for Cold Technology. Those measurements were performed
with thermal conduction probes having significantly smaller diame-
ters (about 1 mm) than those used in the previously described ex-
periments (about 2 mm). The experiments were performed at various
degrees of saturation and for different dry densities. Only the
measurements performed in dry materials will be discussed here.
For such materials, the small probe diameters used should not affect
the results.
193
van Rooyen found exceptionally high (conductivity) values for
materials with some-what special properties, where a certain
amount of clay was mixed into otherwise sandy materials rich in
quartz. These materials were used as fill around high voltage ca-
bles and high thermal conductivity was specified in order to pre-
vent overheating even when high powers were transmitted. The
materials were called "thermal sand". Figure 14 shows results
obtained for these materials, along with (results from) van Rooyen's
measurements on one type of quartz sand and crushed quartz, as
well as measurements on crushed quartz performed at the Institute
for Cold Technology.
As can be seen, van Rooyen's results for quartz sand and crushed
quartz are close to the results for crushed quartz describet ear-
lier, while his results for the two special materials show sig'nifi-
cantly higher conductivities. Of these materials, "thermal sand I"
had the largest clay content (6.0 percent particles < 5 um), while
"thermal sand II" contained 3.2 percent clay. Particle size distri-
"butions for these two materials are shown in Figure 15.
•, 0.... .. . , "-
2.0
1'o2
0.4
"0. a
i,•]
,l, L~utninj:=•w~nvn ity "Vhvrnt). ,,nt" sl¢,,,un•.ikiiet rMtj
of
;,,,.aj .I.t .'t-'100
l 4rd '1''aLSV.V ,"'Iifrt td a I
(I 4Q r I '11111'1& O j¼ 4'I I- it
q~
I~i'Li
194 J
I m. ', P I . ' 1 . T
, ,o 106 ". 1f
0 A. , 11,4 NI 1
' .....- I
,, . . .-
r~'.a' LOLi~i
"tiliriD CU III
195
In all these cases, the conductivity was determined for materials
which were dried after kaolinite was added. Farouki also reports
ono test performed on a dry material into which dry kaolinite
was mixed. This mixture showed a significantly lower conductivity
than the corresponding materials which were mixed under wet con-
ditions and subsequently dried. Measured values for a dry density
of 2100 kg/m 3 were 1.25 and 2.5 W/mK, respectively.
1!
iiw
4
ItO 16- Z10 aoa v t, irustru Wurmn i toft%*t Rin met i*-
8"11011r~ ~i r
UV fi'd L , SIt~ Ia) ~
00 *i)P1J tP'"W. L., r.
196
random distribution of clay particles. The occurrence of an opti-
mum clay content can possibly be explained by assuming that in-
crease of the clay content above a certain limit will result in
the sand particles becoming separated by "cushions" of clay par-
ticles. This hypothesis is supported by "packing" experiments per-
formed by van Rooyen on samples with different clay contents. When
the latter increased beyond 8 percent, the densities decreased con-
tinously for the entire sample (mixture) and thus also for the en-
tire "sand skeleton'"2). Up to that limit, the density of the entire
mixture increased, while dry density based on the sand alone was
nearly constant. In this range one can thus assume that the sand
skeleton is independent of the added clay, since the clay particles
are deposited in the pores. When the clay content exceeds 8 per-
cent, the clay disturbs the sand skeleton and separates the cand
particles from each other.
F. Conclusion
197
to the minerals with high quartz content for which the measurements
were made.
300 -
, JOK • 7 .---
10
0
0 0: 4 g O6 l 01
198
A"
2. SATURATED MINERAL SOIL MATERIALS
A. Introduction
-
; 11 0............
Ila .li v ih tl- irlJ•, F..tn"v',•l r'r fill•• va rro tot. | 1ovi t-
I '+ ' ;• . 1W &,11.
+ L nL I
lt~lJ+ .'lu .l
Vio '.1 l'r I , luril
199
Evidently, the limits are particularly narrow when the material
is complotely saturated by iCe, but also for water saturated ma-
terials the limits form a basis for relatively accurate estimates
of the conductivity, without regard for differences in microstruc-
ture.
In chapter II, results were also presented which indicate that the
geometric mean of component conductivittes, calculated on the basis
of the relative volume for each component, will give a good appro-
ximation to the conductivity in such materials. As shown in the
figure values calculated from this model fall within the Hashin-
Shtrikman .limits over the entire range considered. The validity
of this model will be discussed in the following, based on results
from experimental investigations of conductivities in saturatoA soil
materials.
200
1/,
0 .. .. .. ... .. .
4t
Lr*ýlC..-I *s.Ift0 :1 11
of )at~e I I
r'? om
aCt..b~ 81; k",~O
NIext to Figure 19, the soil materials concerned are listed accor-
ding to decreasing quartz content, with the symbols used in the
figure. By comparing the plotted data ,.nd the table one finds that,
with few exceptions, the conductivities are grouped according to
increasing quartz content,
201
A, A
,
, //"/ le//
0 Z/
) O• 40 4. 20 2• 20
MALT LIONIPdO!SV'4E
S')l•~iý)5 '•r,!i v • :1
genv,•u •). ! .",)I Vtl~, '.
lr. b•,s1 ! n,) " Or!
o
'As can be seen, the calculated values fall largely inside error
limits of between 15 and 20 percert of measured values. These
error limits contain both the measurement uncertainties and the
uncertainty of the mathematLcal model used.
sand, one can assume that all the water is frozen. For the four
soil materials studied by Kersten, which meet this requirement, one
obtains a difference of less than 10 percent between calculated
values and those extrapolated to full saturation by Kersten.
C. Ocean sediments
202
Earlier investigations of thermal conductivity in ocean sediments
have resulted in rather stable values which mainly vary with water
content in the samples (Ratcliffe 1960 (16)). Similar values were
found during this investigation in the eastern and central parts
of the region covered, while measurements from the western part
showed significatly wider variations. A few samples gave thermal
conductivities more than twice as Iigh as for normal ocean sedi-
ments. The authors refer to geological surveys in the W/es- Atlan-
tic which revealed thick deposits of quartz sand and silt with high
quartz content also in particles smaller than 2 ,m.
0177
o i/ VAN
20 40 so 40 00
VAN4IN IHOLO AV VArT VV(
~ ~~~
.1,1~~~ IIv. . e 11 ,.
Yd Yw (1-Ww) 9
Y8 YW (1wW) + Ys ww
203 I
where ww is relative water content by wet weight
yd is dry density
D. Conclusion
204
However, this mathematical model should also be compared against
experimental data from a representative selection of Norwegian
soil materials. To this end, measurements have been initiated
by the Institute for Cold Technology, which will cover a wide
variety of soil materials. This research program includes measure-
ments for varying degrees of saturation (as well as dry materials) as
will be described in the next section.
A. Introduction
2 0 2 0 %-4 '
J-1
.. iis i -h r ua do
iin ra ir iso v
02 0
VU•
[,l', U I.••4eldltl•Pl'-)l[
lYII/tl!l••'ll'. {•I.|'
"1.1|ll"l 205Oi
As mentioned previously, the wide separation between the limits
is due to high sensitivity of the conductivity to variations in
microstructure. Since this effect can not be incorporated into
analytical models for conductivity calculation, due to the compli-
cated geometry within natural soil materials, one must here resort
to purely empirical methods for describing the conductivity depen-
dence on the degree of saturation.
X2 10
206
B. Study by Kersten
of silt and clay fractions. This would indicate that this kind
of representation also may be used for studying the effects of
different textures on the conductivity of moist woils.
207
06
;,! .
03 //
-
"+'.
010. of 00 O ./a*
0I t1
0 ---
-'6 .II 0 f; it,1,1.i..a.....a.If
0 .US ' 0, -" 0 3•
0
WOIWI I. IOIA
*I'•lj .rtlL,,.., r ' i20
.'.,.43'/l•I : ,'sU .'i P.t I
If4+ m ('P at-Ik 1, 1 .ýfft~~
W0,, 4 *1C,+
, )• *+l•
208l,,• ,l l{~1 + i ;A~+1., . .f 1
For the first group - fine grain materials - the spread is rela-
tively small around a mean curve given by
Ke % log Sr + 1.0 11
Sr > 0.1
Sr > 0.05 12
• 06'*
0-
0 0.4 06 08 '0
METNIN050AA0 I
209
As evident from Figure 25, all the results fall near a common
mean curve (line), regardless of the previously used separation
according to texture. The mean can be represented by the simple
relation
Ke - Sr 13
with a variation somewhat less than +tKe a 0.1. This simple re-
lation follows directly from the nearly linear relation between
conductivity and degree of saturation shown by Karsten's data
from measurements on frozen soilso
210
G;0"VK0ftNET
J * FINKORNET
#0
09%
6, ZiIlnrt
1Z. lktfi Av tw"pUOLtull.Iti A WIinini!,swinui Il
*21
!
flNKORNZr
EP ~
1 2 0 20 % 406*01
M!CTN164GS(RAO
Dt('Rq C)" 5AI.'ATION
Also in this case there seems to be different trends for the two
groups of materials. However, the decreasing conductivity of air
with decreasing temperature seems to be the dominant factor at low
and moderate degrees of saturation. For higher degrees of satu-
ration, the increasing conductivities of mineral particles and
ice with decreasing temperatures become more and more-important.
In the fine grain materials, the latter effect is amplified by
the fact that the relative content of u~n-frozen water decreases
both with decreasing temperature and increasing degree of satura-
ti~on, For the lowest degrees of saturation, all the water may be
frozen at both temperature levels at hand, For higher degrees of
saturation, a significantly larger portion of the water will be
frozen at the lower temperature (-300 C) tihan at the higher (-.46C),
which contributes to an increase in conductivity with falling
temperatures, Also in this case one finds that the increase in
conductivity seldom exceeds 10 percent of the value at -4"C,
212
.0,in~
Ket'u.
%(elit
e.07 wl st
I,,*JQ e I~* 'A l
" "M4
"IL 00 100.0 00 'JA 2 47 77.2 20.61 . 3i.0 -
Niti me
11nI0uUgq
berOl 10.2 17.0 21 W to -) 7.71 60 1 13.0 - 17.1 174
4
14 .Y
W rikilhall'O no 0.4 S3.61 2) 4 IRi T.111 51.3 11.4 16A)
0 'l*I.~Vq 1,11hlkifl tria 1.0 21.0 tK4. 1301 7.)n1 i.s 3i,1 19.5 71 '1
i t 1 1
jf IllOt~i1
1
IluVl IM 0.0 7.0 Vill) 11 5 2 10 04 1. 1
,'. ',i...
'lI.I1,10yn 0.0 91.2 43 U 2117 '.71 1141 - 3.? 2.2 21h.0
213
C. Investigation of Norwegian soil materials.
214
Table III. Summary of quartz content in the mineral library
at the Institute for Cold Technology
CLAY 7 10 - a8 18.8
SILT 6 21 - 53 33.5
SANJD 13 10 - 100 44.9
GRAVEL 13 2 - 57 38.8
MACADAM 10 2 - 53 12.7
TOTALS 49 2 - 100
As shown by Figure 28, the selected coil samples cover marine se-
diments and glacial flow materials rather well, while moraine clays
and moraine gravels (not separated into layers) are poorly repre-
sented. In terms of quartz content, the ran!;e of interest seoms
to be well covered for all types oi soil, porhaps with tho excep-
tion of quartz-poor loose deposits from regions characterized by
special petrographical conditions. lowever, the macadams include
a relatively large selection of materials with low quartz content.
215
'0
1s
20- I
180
- ..- 2 2
0.2 2 20 0 2 -Tim 2 20
M E0 AN Nid=d
PIG. ~ eri~iv a.~ jori- tf.ri~lt~r f rit~ r 1~b
v&, r.: u ctLr k. .,:tkvi~kk,;,TH{. ~ 'rI2 ru:
i YM1 ti]. .atrialbibll±otzk~t,.
216
4
D. Conclusion
217
4. SOIL MATERIALS CO'ITAINI.NG 1iUIUS
A. Introduction
On the other hand, certain reliable results are available for pure-
ly organic materials such as pure peat and bog. These will be dis-
cussed in the following.
218
.........
. .....
the conductivity of hardwood (18)IJ, Zincc tee conductivitjy for
water is 0.57 W/mK (at about t46C), the previously mentioned
Hashin-Shtrikman limits will be much narrower for water-saturated
peat and the geometric mean equation can be used without major
errors. The same will be true for ice-filled materials, since
the conductivity of ice is 2.2 W/mK. (at about 0CC), so that the
conductivity ratio for the components will not exceed 5.
219
/
0 1 @4 of
"4 IQ
I
PIO. 29.
.o@ -,7 . pA... ramsse.
m1411Uinsge . . tarv- ap .IO pJ.Out 00111.
funWLajon av IlnvlninCfrrA.Ion. '!jgartimoplee oil
fP'aaeii peat JZ& LOPIO illato~i to d 7l'0 o",
sdaturation.
",, 1 6AI-.... -
1.140
I 11--! i,PIA..
.. V .. . t .. . n:, . .. 1. L 1...
,I , roatta lima
0 1 vu
' 4 &Lou :j ,j. ae o'f ial ti'r1 1 on.
11A
4this
materials, computed from the geometric mean equation.
results in values which are some-what low for dry materials,
since the effects of thermal radiation are neglected.
As expected,
However,
the trend as function of porosity will be the same.
220
iI
r,, 1.0-
nj. PoJ-o6%1116 ,
Pl.3 1. IIfIammannm,, mlj l~oiJo Ivininwiic~vuri oir. n '1litet~of i~v
One finds, that within the normal rancte of-poroalties, the conduc-
tivity of' water saturated materials will vary less than 3 percent
around the mean. For frozen and s~aturated materials, the correspon-
ding variation is leus than +10 percent. In both cases, the conduc-
tivity tends to increase with increasing porosity since the organic
component has lower conductivity. The variations for dry materials
are larger percentage-wise and the trend is the opposite.
-(un-±ozen Sr soils) 15
221
Ao = 0.05 .v/.I:
= 0.55 W/mK (un-frozen)
= 1.30 .lmI (frozen)
REFERKICES, CHAPTER V
222
S• • ,-
13. N. Wakao, K. Kato: Effective Thermal Conductivity of Packed
Beds. J. Chemical Engineering of Japan, 2, (1) 1969, pp. 24-33.
14. L, Barden at al: .The Collapse Mechanism in Partly Saturated
Soil. Engineering Geology, 7 (1), 1S73, pp. 49-60.
15. P.W. Kasameyer et al: Layers of High Thermal Conductivity in
the North Atlantic. J. Geophys. Res. 77 (17), 1972, pp. 3162-
3167.
16. E.H. Ratcliffe: The Thermal Conductiviites of Ocean Sediments.
J. Geophys. Res. 65, (5) 1960, pp. 1535-1541.
17. Y. Nakano, J. Brown: Mathematical Modeling and Validation
of the Thermal Regimes in Tundra Soils. Arctic and Alpine
Research (1) 1972, pp. 19-38.
18. A. Missenarii Conductivite Thermique. iditions Eyrolles,
"Paris 1965,"p. 356.
.1.
---. •..........
H 223
METHODS FOR COMPUTING TIHE:!AL CCNDUCTIVITY
OF SOIL XATERIALS
1, EXISTING METHOD3
A. Analytical methods
224
4.
an analytical expression for conductivity, with di3tance Yo as
parameter. This parameter defines the distance from the center
of each sphere to the plane intersecting the lower "meridian" boun-
dary of the meniscus, normalized to the radius. The relation be-
tween relative water volume and Yo can be deri.ved from geometrical
considerations and is
This equation is only valid when adjacent menisci are not in con-
tact with each other. The expression for conductivity has the form
-(f+g) (f-g YO o h +, h
X
a g (f-g) (f+g yo) fh2 h - h Yo
*~~~Q~j
~J~I~iTN.O
225 •
Gemant carried out a numerical example for ;ihich the particle con-
ductivity was assumed to be 4.2 '3/mi,. Figure 2 compares the func-
tional dependence of these results with Kersten's experimental va-
lues for a silt/sand mixture (Ramsey sandy loam), which has about
the same particle conductivity. The lowest density in this test
run was about 1600 kg/m3 and thus some-what higher than that given
by the mathematical model.
*,*, I= , , t-- t --
W/ni •ItUn
0. 04
II.
aA.Ii 6
(Qoytlll
2 1
226
piaical. .akozki an. uohlinszy have developed a relatively nim-
ple equation# based on curves fitted to data obtained from Gemant's
modified model. This equation will be treated in next section.
Alberts (1966) has proposed a similar model, where the soil par-
ticles are represented by tiohtly packed spheres in a cubic lattice,
while the water is partly spread out as a film on the surface of
each particle and partly contained in a meniscus around each con-
tact point (3). Part of the calculations were carriea out numeri-
cally but the results are of little value due to the very low con-
ductivity assumed for the particles (x a 1.0 W/mK), as well as some
approximations which are introduced based on this assumption about
the particle conductivity.
L nipi
1-o
where ni i4 relative volume of coniponent "il
± is conductivity for component "i"
-.1
F .. + (!i/x 0 - J , j albc
227
a ÷÷ c -= 1.0
gl x a .10 , g3 .
C.80 (Dry soils) 5
228
-A
. .4 4/, I
060
06 t
1
O.J2 0 O&O
"),.
4''-Ac
~
~r ia
4I~14or t I
pro444 4 v ,..
These graphs show that the two equations deviate from each other
by less than 5 percent within the ranges of particle conductivi-
ties and porosities that are typical (may occur) in saturated soil
materials.
De Vries suggested that the equation also could be used for calcu-
lating conductivity in moist soil materials. He separated out se-
veral degreesof saturation ranges in which the three components,
particleswater and air, were treated differently with respect to
geometrical conditions. For low degrees of saturation, the partic-
les were considered to be enclosed in water, while the air was
considered to be the continuous medium. For high degrees of satu-
229
ration, the water was considered as the continuous medium, in
which either the grains of sand or air bubbles are contained in the
form of particlesI)
This discussion has shown that analytical methods are useful only
if they can be fitted to experimental data. As will be evident
from the following section, empirical methods also have their li-
witations, due to the large number of variables which occur.
B. Empirical methods
230
.A
these fractions. When these equations were developed, results
obtained for two of the coarser materials (Northway sand and
fine sand) were not included due to the special mineral composi-
tion of these materials (Kersten, page 84), i.e. due to their low
quartz content. This gives an indication that the chosen separa-
tion (classification) according to texture in reality means a
grouping according to particle conductivity.
where c and B are empirical factors which are different for coarse
and fine soil materials:
a b a d
Fine grain soils 0.0014 1.4 0.012 0.05
231
..... I.....! I•
"ii "1 |..........................- .~"=-"..-.i -----. • , " -- f---.-.•--• .'= ,. ... : ,
-• • , , I If I I I II I I I . I
10'ky,mnl' .
\ \
t-'a,'.A ' tO I
4'\ $q
I 01
1..4
~111
0~ ~111' 1 i % 20 1 oAS 4
32~"US ' ~ tPK54 VANN
b)ni Lo
it-n av f vK".-dl-*f
232PC~VN
The previous analysis of 1Kersten's experimental results showed
that the thermal cor.ductivity of saturated soil materials can
be determined from the geometric mean equation, when quartz con-
tent and volumetric compositions are known. Figure 5 shows va-
lues for saturated materials calculated from Kersten's equations,
compared with data obtained from the geometric mean equation.
i0 j"
nod,- • .;"
VCisW,IACI ý "411L1
0Bi ',/r #
I 11i
z . t
-'a 1.)6k v -'di vfCfoa fo UIo -09'1a
110. U.#!rd.
4. 14lo iic d envIdr kgr vtnnmvtt.d$ Joret&v.i. iom -
This comparison clearly shows that the two selected soil material
groups represent differences in particle conductivity. The equa-
tion for coarse materials gives an average particle conductivity
of about 5 "//mK, while the equation for fine gqrain soil gives
about 3 73/mlt. From the previously proposed mathematical model for
particle conductivity one finds that these values correspond to
quartz contents of 60 and 30 percent, respectively. An evaluation
of reported values for the quartz content gives a mean of 57 per-
cent for the five coarse grain soils considered. The correspon-
ding average for the four fine grain soil materials at hand is
uncertain due to incomplete information about quartz content for
two of these soils.
233
------------------------------------------------------ i'--
value for quartz content in coarse grain soils that is obtained
here from Kersten's equations is not representative for Norwegian
conditions. Figure 6a shows a cumulative distribution curve based
on quartz contents in coarser materials determined by Sveian.
- to 4 4 47i - 0
4 ,i i I lllS 0 ll)0 I l *
fIl' o l b) c. - .lr*
• . 1 j-1rQ,''+vd
234
this definition, the relative contents of the fractionz d <2 irm,
2< d< 20 um and d )- 20 ' m can vary within aide limits, Since these
three fractions show such considerable differences in their dis-
stribution curves for quartz content, the definition for coarse
grain soils implies no definite classification according to quartz
content.
235
"......A
tribution function obtained for coarse materials gave a possible
variation in particle conductivity from about 2.2 to 5 t/i/mK, if
one neglects the most extreme 10 percent of the results.
w/mK
,
,• "M+
wsm 'W,. ,
11 U11111t
a .II , tI i . • .
o 010 0. 0C Oe
O 1.0
AWN,11 FV•IUIUPl
N FOA114i6l,,
One also finds variations around mean curves for the two regions
of between 20 and 30 percent for coarse grain soils and from 10
to 15 percent for fine grain soil materials. The mean curves
represent particle conductivities of 3.5 and 2.6 1V/mK for coarse
236
- +•..,• • , ......--- --•" ',,'r ..........:w .... ... • ' - , ,_• • +...4M'I
and fine Grain soils, respectively. This corresponds to quartz
contents of 40 and 20 percent, respectively.
4.0 i
03V -
F K9. . I I.....
j r . '; ',99
1[,
0 ., 4' 1 ,- ' VO.. %4.JL11 ",
237
particleI) conductivity with a broken line, representing an un-
frozen water content of 10 percent. (Nu = 0.10)
238
IWA
I "",,, I //
_I 239
should thus include possibilities to account for varying quartz
content in coarse grain materials, while the effects of un-frozen
water should be included for fine grain soils. As mentioned be-
fore, it may also be desirable to change the method for classi-
fication of the two groups of soil.
The clay content was said to affect the conductivity of soil par-
ticles. In clay-free materials, the particle conductivity was
assumed to be 5.84 W/mK while a certain clay content c (percent)
was assumed to reduce this value according to the following formu-
la.
X2 a 5.84 - 0.033c 10
240
sand, and for pure clay determined from the equation given by Makowski & Mochlinsky.
Both these curves are close to that derived from the geometric mean
equation on the basis of assumed particle conductivities.
One finds that the curves for clay-free materials on the average
follow the empirical relation for coarse grain materials, while
for pure clay the curve falls above the empirical relation for
fine grain materials.
:'ii.
AATK I LYOLIMA'NU1L
241
OW
0 ;.2 l . CIL
A .4
//o
S0 ', I. . "
0.1 0 ,3 C2 06 01 10
W 41t4I•'4IRA0
242
...
-------------------------------------.---
For this reason, they only treated un-frozen soils.
r - A 1i0 + a (cm/oC4) 11
s - sI - s 2 Yd 12
This in principle agrees with what has been pointed out earlier:
The thermal conductivity of water saturated soil material depends
primarily on relative particle volume and particle conductivity,
while the sensitivity to differences in microstructure is small.
However, calculations using this equation to determine conductivity
of saturated materials with different quartz contents show that
243
van Rooyen's equation tend to overemphasize the effect of in-
creasing dry density, particularly for low quartz contents. This
is illustrated by Figure 13, which shows calculated conductivity
for three values of quartz content and data obtained from the geo-
metric mean equation for comparison.
%-C
2.0
'.01-2.0
.40 '. ' i. 2..
10
Oaa.
k~Ir tiaIllilld. if
Vatn il-Oql
fI'~~~ $if~,'n I),"9
Q4r. t
*ju1i.... ' 1 "'h' • uVIJ "'..
0.21 4
I, M / M'
%.0 0.' 02 04086 t0
244F Sf.
r = A+ S 13
245
In both these cases, the conductivity variations are larger than
those found for dry materials during the study discussed earlier
in the last chapter. In addition, that study gave no basis for grou-
ping conductivities in soil materials according to differences in
texture, except for the proposed separation between natural soil
materials and crushed rocks.
C. Conclusion
kI 246
Another draw-back inherent in the Kersten method is the mixture
of factors related to texture and mineral content implicit in
his equations. An improved method should treat the effects of
these factors separately, e.g. so that both coarse grain soils
having low quartz content and fine grain soil materials with
high quartz content can be correctly related.
247
2. NEV METHODS FOR DETERMINING THERMAL
CONDUCTIVITY IN SOIL MATERIALS
A. Introduction
249
BEST AVAILABLE CORY
Tabpll 11Z. nrciningsrodo11 for Vkr-m.l9dni~evng.I V rv
,rnva.l4e joo r'ur.
"?L
th Trt calculatitg thermal conductivity of irnr soia.
::" ';,,4E T S 5r ..
t . .. . ............... (f
S.. . .
i•.'-'..• 2,2n .0, -nO z6 u ...... (h) : ' :a r' :J ' y
COARSE 2
250
14
mainly from rocks containing quartz, such as granite and gneiss.
Soils originating from alkaline rocks, such as gabbro, will con-
tain small quantities of quartz, as well as considerable amounts
of other minerals whose conductivities fall between the two impor-
tant groups mentioned. Coarse grain materials with low quartz
contents should consequently be treated separately (Eq (J) in Tab-
le 11). For fine grain materials it is proposed that Eq (i) be
used also for low quartz contents, since minerals in the interme-
diate groups seem to occur in small quantities in the finest frac-
tions.
251.
4 ~'i
'HE-MAL CZOJOUCtIvTY OF A'RU OIL
-b
W,,I,
.3
a CRUI
0 20 4.0 50 80 100 12 1.4. 16 le 2.0 01 02 0.4 0.6 10
1. Yvartsinnhold 7. utzoen
~.Orov/firn ~ OZiS/~
.ov
v".k
4.
"L e.JKC-
n i~a
vonvk
f a,:.,
1.
t
4c
*PAVE :77:;.....
VIA V"O
25Ir .
!Z5 ...
1.IE3'.kL. CCVl4D'.-CTVITY OF 'CkC7EN 9YILS
(WIT0* S0
2.0
-~10
- . 2
..
80
2040K0 r .4..
1.6 IIa
i010/f 20p 0 02 04 06 06
Q 20 0 B lU 1.2 1,1
BEST As aILABL
1 rtvtj Sf Inn h' . ~ ,wr o~
I-.~a of e~fJ
E~.n
C~*.0 AVEL
. .
* ,-.j j.~.....,..t.t ij.*2
Ex 2 b '-
-.It'v f~kAY
17-
e.-
253
3. EFFECTS OF VARIATIONS IN SOIL
PARA14ETERS
The two cases Just mentioned can be said to represent two informa-
tion levels, each having characteristic error limits for the con-
ductivity estimate. In chapter III, such levels of information
were given in a block diagram containing the different parameters
which affect conductivity in a soil material. This block diagram
is reproduced in Figure 16.
254
L.,
* •tL'
IO•, I '.'a"%
101 '
Thermal Cohductivity
!P1 .7 , rCNuat v, tIOr' tl. r *e; '.4V
.h' r. wqevrin t1'tr sond to,
255
For sand materials, these results show a near normal distribution,
with a mean value of about 1.90 W/mK and a standard deviation of
0.60 W/mK. 90 percent of the results fall between the limits
0.90 and 2.9 W/mK. The results for clay materials deviate sig-
nificantly from a normal distribution. The mean value was deter-
mined to be 1.80 WN/mK and 50 percent of the results fell within
the limits 1.60 and 2.0 W/AK. However, the rest of the data showed
a wider variation, with the 90 percent limits at 1.20 and 3.20
W/mK.
.1 I I'/ '
•0 r' . ...
" ... . . .....,. ... ..I.. .. !
257
.10
' 0
4~
20
* VANOIVNNPýO00 VE(T %,
d'~
.9iJsu
IJr w t 01 )atvP IJji uj 1
258
Little data exist about the distribution of dry density in road
beds (pavements). It is proposed here to assume a distribution
having a mean of 1900 kg/m3 and a standard deviation of 75 kg/m3 ,
which implies that 95 percent of the cases fall between 1750 and
2050 kg/M3 .
259
' .. ,**,-*
... i. . ... ....
Table IV. Assumed (estimated) distributions for key soil para-
meters, to be used when calculating thermal conducti-
vity. Mean values and standard deviations.
260 F
Table V. Error (confidence) limits in calculated conductivity
of soil materials. Un-frozen sand and clay.
Level 1 2 3
1. Texture Sand clay Sand clay Sand clay
2. Dry density 1770 1400 1770 1400 1770 1400
A d 340 180 -0 -0 -0 -0
3. quartz content 0.37 0.20 0.37 0.20 0.37 0.20
Aq 0.30 0.16 0.30 0.16 -0 -0
4. Deg. of satura- 0.50 0.90 0.50 0.90 0.50 0.90
tion A Sr 0.40 0.10 -0 -0 -0 -0
5. Conductivity 1.50 1.15 1.50 1.15 1.50 1.15
(x/;0) L o 0.056 0.005 0.03 0.004 0.03 0.004
6, 1
(x/aX ) X1 0.498 0.387 0.292 0.271 0.097 0.127
(0X/aKe) AKe 0.393 0.187 0.150 0.103 0.150 0.105
8. Relative error
percent 42.0 37.5 22.0 25.2 12.0 14.3
One also notes that the errors are about 25 percent in cases where
the quartz content is not known, while dry density and degree of
saturation are given. This level corresponds approximately to the
level of Kersten's empirical equations, for which the uncertainty
also was found to be t25 percent.
261
Table VI. Error (confidence) limits for different degrees of
saturation. Un-frozen sand.
From Table VI and Figure 20 one can see that the improvement in
error limits due to knowledge of the quartz content is largest
for the highest degrees of saturation.
Error limits for the conductivity were calculated for the varia-
tions in dry density and degree of saturation estimated previously
(see Table IV) and four different cases, see Table VI. For the two
first cases, dry density and degree of saturation are assumed to
deviate from their means by one a# while in the other two cases
these parameters are assumed to be known. For cases 1 and 3, the
quartz content deviates from the mean value by one • while it is
assumed known in cases 2 and 4.
262
Y-
3
- I' "gm . - f.D - (KVAR-SNNHOL. 37%/,
-QUA KJ IGON r r,'NT
e ,i: 10 --
0 ,, I I
iJ | vrun
enlwvd
In W
Io n a •s U'
t re,'ii.i n,i
I~ J.
r t j ,I1i
.' :wr,
•
re.+ivr, iiv
W , P C. 0 "1
}0.
$ 40 1 t
,..rI
. I. ?' a t4C u a
Case 4
263
when data for degree of saturation and density can be given exact-
ly. Thus, it is more important to obtain reliable data on den-
sity and moisture content in road surface layers than to accurate-
ly determine the quartz content.
2.0.
w.f .'"
f0 - - •. ,
Dry Density
264
dominated by the uncertainty in dry conductivity. This is also
shown in Figure 22, where error limits fq r the case of lacking par-
ticle shape information is compared to those obtained when all
relevant information about the gravel ma'terial at hand is available.
S'1' \ ,~
\I C
I- N
7! .13 i
I._ _ _ _ _ _ _ _ _
,
S.... . ..... . . .. ... ....... . ,•........ .. , ? - . , . •: ;=+.,+.•.++.•f•..._.•. . ..... ..... ...
S• ~
: ,: 265
- '20 S, ..
B. Conclusions
S1The in the previous examples show that
error analysis performed
In cases where one lacks relevant soil data one must resort to uaing
average values for the conductivity, with inherent uncertainties.
266
Such values are listed in Table IX. Average values and variations
for the soil parameters on which these conductivity calculations
are based are also included. In all cases, the variations repre-
sent one 0- limit in a normal distribution and will thus corre-
spond to 95 percent of possible cases (occurrences). This will
then also hold for the calculated ranges of conductivity variations.
The results are represented in terms of histograms in Figure 24.
267
Free
ground Upper I
sub.-soil -
W r
0+0 0
+I+ 000 H +
SI I0 t0 I +
HS 9+
0 + C+
0, I+0
C~0 0 .J I+ C.1+J
, C . .1 1
1+~0
1+ o+ o+ 10 1 !
I 03
. ýj1 %J1j 0 0 0
0 0 0
ci.000 v0
268
* -A
VA
FI
VA, c [:\ rq'F1 N([ ].\, bJU / rb•,L ' .'-•L ,.." ,.. v.v,' Y ',.'),/ <
0 0,5 10 15 20 2S 30 35
. 1' - Ir T*N. . V.
NATUPGRUS I __.
5",N / FILTERLAG FM.
f R.H
I Ft.
S ILT
. L il
toL', ; V tt V), ," J;,cnl ;' t .. r to $ . E ( . *'.7 V)
v•L i.
269
, il
4. EFFECTS OF CONDUCTIVITY VARIATIONS
ON TTIE1MAL CALCULATIONS
A. Introduction
270 "U
A
Effects of variations in the soil parameters at hand or uncertain-
ties in estimated conductivities will in the following be trea-
ted for some' selected cases, limited to road foundations where
freezing of the ground is to be prevented.
271
-rmi
while conductivities of sub-soil and the frost preventing layers
themselves have relatively small effect. The latter is particu-
larly true for bark layers. The large variations due to upper
layer conductivities indicates the importance of obtaining accu-
rate data for dry density and water content of materials used in
these layers when cold accumulating materials (layers) are used for
frost protection.
272
vestigato this case, calculations were performed for a structure
a
consisting of a 5 cm layer of insulation mater'ial (plates) and
* pavement (upper structure) identical to that assumed for earlier
* calculations. According to those calculations, this kind of struc-
* ture can withstand a frost quantity of about 30,000 hOC with no
frost penetrating the sub-soil,
L W A
P 30 000 hC
It,.,
I"
C
440'
:lot
is 0 64
ItOS )
It'ICI~~m4Wi4 -30 30 OS 2 4
P
NIORS. 'MA
I IiUU
"11 ,i S OM1091 02
It 1
VAMIa00
9.70 0
34it~s
A0 ol roL1 C. L
*s oil 'A 4 J
4 4..~-t 0.0* ItOJ 0 & I'iIm4ij
WVo
IJAOCRM273'
If the structure is subjected to a frost quantity of 40,000 h"C,
the sub-soil will consequently freeze to a certain depth, depending
on how large a portion of the available water actually freezes.
Figure 26 shows calculated values for the thickness of frozen 1
sub-soil, as a function of percentage un-frozen water in the sub-
soil, for this case. It is here assumed that all water "available
for freezing" has frozen at 0C.
IlOl60 . . . . . til
a.10
0 10. ...... . .. . . .-
" 0
O 10, ~~~......... 1•0
...... . ..
* kflQEy VANN V06.i
• • 1" .h
1 41 %U'1
274
:
For example, the effect of the temperature distribution in the
sub-soil on the un-frozen water content is uncertain. As men-
tioned befor% estimates of the contribution from (stored) soil
heat is another uncertain quantity. However, calculations of
the type indicated here should be performed as soon as experimen-
tal data are reduced and the computer program being developed by
the Institute for Cold Technology becomes operable, in order to
clarify the effects of under-designed frost prevention layers.
C. Conclusion
275
REFERENCES, CHAPTER VI.
276
14. S. Skaven-Haug: Frostfundamenters dimensjonering. Frysevarme
og jordvarme. Frost i jord, (3) 1911, pp. 9-27.
15. A.F. Knutson: Theory and Experience Regarding Frost Penetration
and Frost Heaving, Symp. on Frost Action on Roads, I,
Paris 1973, pp. 223-233.
16. A.K. Flemings Applications of a Computer Program to Freezing
Processes. Proceedings of the XII International
Congress of Refrigeration, Washington, D.C. 1971.
277 I
I
THERMAL CONDUCTIVITY OF SOILS
SUMMARY
278
zY~-
be said to give exact results because of the problem of
incorporating the microgeometry of natural materials in
an analytical solution of the heat conduction problem.
Such bounds as those developed by Hashin and Shtrikman
Fig. 1, I1 (1962) shows that the conductivity is very sensitive for
variations in microgeometry when the component-conductivities
are widely different, while the sensitivity is small when
the component-conductivities are of the same order of
magnitude. It follows that the conductivity may be
calculated analytically if the component-conductivities
are of the same order, while only empirical methods can
be used under other circumstances.
279
EQ. 5, V On the basis of the experimental results, two empirical
EQ. 6,V relations between conductivity of dry soils and porosity
were derived, with different relations for crushed rocks
and natural soils.
S280 t
A
temperature below freezing, was found to have an appre-
ciable influence on the thermal conductivity. This influence
is incorporated by a correction in the saturated conduc-
tivity, on the basis of volume fractions of unfrozen water.
2
281 I•,
An extreme situation is present when no relevant
soil data is available. Under such circumstances average
thermal conductivities for the given soil group must
be used. Principally, such averages may be derived on
the basis of statistics on the distribution of the
relevant soil parameters for different soil groups. Such
data are scarce. However, in conjunction with an analysis
of expected variations in thermal conductivity of different
soil groups, such distributions are sketched partly on
the basis of available data.
282
-VAN
give similar improvements when dry density and degree
of saturation is given in subgrade materials with high
degrees of saturation.
283
-
POST-SCRIPT
284
In order for the model to be used also in cases where one lacks
information about materials used in the structure, it is neces-
sary to have relevant data on expected mean values and variations
of soil parameters for various types of materials. In this report,
such values have been used to determine expected variations in
conductivity for different types of soils. However, in some ca-
see these values had to be estimated, due to lack of information
about water content and density in sub-soils and upper layers of
road structures. A study initiated by the Institute for Cold
Technology, NTH, may supply the required corrections for these da-
ta. When such results become available, a new error analysis
should be performed for the model, using revised mean values and
variations for the soil parameters.
These concluding remarks show that the work reported here is part
of an extensive project, which is by no means finished by comple-
tion of this report. As indicated, several on-going invesgations
may form the basis for necessary improvements in the result given
285
here. These results will be reported in detail when the project
"Frost in the Ground" is concluded at the end of this year (1975),
e.g. in reports issued by the "Group for Analysis of Frost in the
Ground", Institute for Cold Technology, NTH.
ae8 1 le
ma +
;,
e (R,t) - eo sin (Wt) 2
To solve this problem one can use results from text-books such
as Caralaw and Jaeger (1) and mathematl.cal handbooks (2), as well
as general control theory (3).
286
From control theory it is known that a system subjected to an
applied, stationary and periodic forced input will give a eta-
tionary response which is periodic and damped (attenuated).
For the problem at hand this implies that the temperature vari-
ation on the axis of the cylinder (r m 0) can be written (4):
G(s) - •os
S(R, s) 1 0-_(R V-3-7a)
-argjG(iw)j
since (1 + i)2• 2
This follows from the rule that the absolute value (magnitude)
of a complex number l/z, z = Xe i, equals I/X, while its argu-
ment is-
The same method can be used for finding attenuation and phaseshift
between two arbitrary radii within such a sample. For this case
the transfer function can be expressed
GC,2(a) "----
ýW-•le "o(ri_ •,M ,
(8) 42(s) IC(r2 ý712
289 *
2 2
(beix) + (berx
* (beix) + (beix)
AI,24*~. . - - Z
, i'
wher'e x, r, =w'. X2.. =r=
I~
.............. I--,
-.-........... .-.--
*g
4c *i.,__................ . .,+++
+ .Wll4I
, M. - P- 1 16,. t. • 4p 6.1 .1
F1 l "I
29,0.
41 y
, i fl.o.,•
rL Ie oo
'. '@l :4l ,1 ; I I !.I •+
r, , I , i"jI
BESTAVAILABLE COPY
29 0 t
REFERENCES
291