Supersaturation Cause CC4 Around 45oc
Supersaturation Cause CC4 Around 45oc
00
Printed in the U.S.A. Copyright ~ 1988 Pergamon Press pie
A. BRANDSTE'ITER
Department of Engineering Physics, Research School of Physical Sciences, Institute of Advanced
Studies, The Australian National University, Canberra, A.C.T., Australia
Abstract~The long-term stability of calcium chloride hexahydrate phase change material, as a heat
storage medium for solar and similar applications, was assessed in a 400-g nonstirred sample of for-
mulated material. Accurate monitoring of enthalpy and temperature profiles over extended periods of
daily melt/freeze cycling was achieved by a unique method of continuous thermal analysis. Results lead
to the conclusion that in the presence of supercooling suppressants, such as strontium or barium salts
and possibly sodium chloride, the formation of calcium chloride tetrahydrate can be inhibited totally by
means of extra water above the hexahydrate composition. Heats of fusion and crystallization in excess
of 200 kJ/kg have been routinely registered.
I~0 / Temperature T
°C
weightration
(a)
weight %
(b)
degree of
hydration
(c)
o120 : L, '
~I00
J,
0 0.595 37.3 10.35
5 0.62 38.3 9.94
.80 10
15
20
0.650
0.69
0.745
39.4
40.8
42.7
9.48
8.93
8.27
60 25 0.83 45.4 7.42
28 0.93 48.2 6.62
30 1.02 50.5 6.04
20
I I I
10 20 30 /+0 50 60 rium state will as a rule feature a phase of "super-
TEMPERATURE ('C) natant" liquid. This supernatant--and also intersti-
Fig. 1. Calcium chloride solubility curve, after Pauling[l]. tiM--liquid actually is saturated solution in equilibrium
with the solidified hexahydrate and as such its com-
position, at temperatures below 30°C, is as given in
30°C, and our interpolations for 5, 15, 25, and 28°C Table 1, independent of the relative amounts of solid
based on Fig. 1. Added in Table 1 for convenient and liquid present. This was verified by examining a
reference are (2) weight-percent values in grams CaCl., number of samples of varying water content--both
per 100 gram solution, and (3) the hydration values in excess and deficient relative to the hexahydrate
in mols of water per mol CaCI2 corresponding to the composition--which all were found to have very
weight ratios. nearly identical liquid phase compositions.
The point at 30°C, for concentration 50.5 wt%, is In a system of solidified CCH and supernatant liq-
just beyond the peritectic point, P (Fig. 2), which uid, weight fractions (s) and (1), with CaC12 weight
according to M e l l o r [ l l ] , citing Roozeboom[12], is fraction X, in the solid and Xt in the liquid phase, the
located at 29.8°C and has the weight-percent con- overall CaCI2 weight fraction X is given by
centration 50.25, as compared with 50.66 for the pre-
cise 6-hydrate composition. (Lane[7] now quotes x = (s).X, + (1)% (1)
29.45°C, at 49.62 wt%, for this peritectic point and
29.6°(] for the CCH melting point). It should be em- The weight fraction of supernatant liquid is thus
phasized that all the values in Table 1 specifically given by
refer to CaC12"6HeO as the solid in equilibrium with
(l) = (X - X,)/(X, - X,), (2a)
CaCI2 solution at the given concentration-tempera-
ture points. the ubiquitous "lever rule," as for example quoted
When calcium chloride solution at or near the by Abhat[3] and Kimura and Kai[4].
hexahydrate composition is cooled below freezing Xt is as obtained from Table 1, (b), and Xs is the
temperature, for example to 20°C, the final equilib- weight fraction of CaC12 in CaCI_,-6H20 (i.e., 0.5066).
In terms of the degree of hydration H ((c) in Table
1) eqn. (2a) becomes
6.1606
0.5066
H+6.1606
, ~ 6C 2H/ (l) = (2b)
Xt- 0.5066
,., .=-t~H
Figure 3 is a graph of eqn. (2b) for a number of
I 30 29.5 29 ~CLI00 /
100 30"C
8o
g 50 , , ,
6.02 6.06 6~o 6.~c. ~o zo . 30
~9 ~'C " * " ~7"C TEMPERATURE
(C)
"7 /29C /28C
~0 /
20
Fig. 3. Weight percent supernatant CaCI2 solution vs. "hydration number" (tools H20 per mol CaCI:).
For temperatures in the 0 to 300C range. Plot shows calculated curves (not straight lines) and experi-
mental points. Inset: CaCIz-6H20 solubility curve used for these calculations.
(1), Moreover, even though the weighing may be at ture of 330C at which the points determining this
room temperature, 20 to 25°C, the amount of super- density relation were obtained was selected for prac-
natant liquid found in a sample largely reflects its tical reasons: since our samples of CCH usually were
amount at the time of its formation following freez- stored in a "hot box" at approximately this temper-
ing, at which time the temperature may have been ature, this curve provided for fast determination of
much lower. Our data points close to the 0°C and concentration by simple weighing (to 0.01 g) of a
10°C curves suggest this is a likely explanation, since given volume (approx. 250 ml) of any sample from
some of the samples were, in fact, allowed to freeze this hot box. Actual states of hydration, for example,
outdoors. to establish the density vs. H relation of Fig. 4 were
A quantity of interest directly related to the frac- determined by a dehydration method, in which the
tion of supernatant liquid is the excess of water caus- residual solid was subjected to 350 to 400°C for 15
ing it. Up to the peritectic concentration this quantity min following the initial boiling off. The curve in
is simply given by Fig. 4 represents a best linear fit to the experimental
points shown.
Weight percent excess water Note that the material used, although essentially
( H - 6) x 18.016 x 100 "technical grade," containing in particular close to
(3) 0.5 wt% NaCI, is relatively pure CaCI2 aqueous so-
219.086 + (H - 6) x 18.016
lution, as evidenced by the transition temperature in
Figure 4 presents this weight-percent excess water excess of 29°C of its CCH form. The presented den-
function for convenient reference. Figure 4 also of- sity vs. H relation should, therefore, represent fairly
fers an empirical relation of the density of CaCI2 so- closely pure CaCI2 solution.
lution as a function of concentration. The tempera-
3. TETRAHYDRATEFORMATIONAND HOW TO
P R E V E N T IT
% 1. ~ O expt't points (33 C) ' For correct evaluations in the region about the ct-
{,0F "w j tetrahydrate peritectic point it is clearly helpful to re-
draw this region, Fig. 2, on a magnified scale such
as in Fig. 5. In this drawing points P and Q are placed
in accordance with Fig. 1 in G. A. Lane's =Technical
Note" [ 13], purportedly after Roozeboom[ 12], except
Y I I
that point P here has been located at exactly 29.8°C
6.0 6.5 7.0 7.5 to comply with Mellor[11], citing Roozeboom, as al-
tool Hz0/tool CaCI2 ready mentioned.
Note that very accurate location of melting points,
Fig. 4. Calcium chloride solution specific gravity vs. "hy- etc., is here of little interest since we refer to tech-
dration number" (tools H20 per tool CaCl2), empirically
determined for "technical" material near the CaCI2"6H20 nical CCH; and second, that the 13 (and T) tetrahy-
composition. drate species (see Fig. 2) can safely be ignored as
186 A. B RANDSTETTER
alistic" geometry of the sample: an upright cylinder respectively, as given by Lane[7]. Figure 6 to 9 in-
72 mm in diameter filled with the P C M to about 70% dicate how these beat-of-fusion values were calcu-
o f its 90ram height. Some of the cycles had partic- lated. A lower limit on the "true" heat of fusion of
ularly fast or particularly slow heating or cooling CaCI2-6H20 is estimated in each case as the observed
phases within a range of about 15 to 43 °. value scaled upward according to the amount of su-
Large amounts of data of enthalpy and tempera- pematant liquid corresponding to the sample's degree
ture profiles have been collected over these 18 months of hydration and that plot's base temperature.
on magnetic disk and partly retrieved in the form of It should be borne in mind that the heat of fusion
printouts and graphic displays. A selection of plots exhibited by a cycled sample is limited by the heat
showing temperature profiles and associated enthalpy evolved during solidification, since only up to this
curves vs. time, for 24 and 12 h melt/freeze cycles amount can subsequently be absorbed upon melting.
for different ambient regimes is shown in Figs. 6 to Therefore, if (congruent-) melting proceeds to com-
9. A very slow traverse through the crystallization pletion, as normally expected, this heat of fusion must
zone is shown in Fig. 10; and Fig. 11 is a graph of equal the heat evolved during the previous solidifi-
enthalpy vs. temperature obtained from the same data cation. The variations in heats of fusion exhibited in
as Fig. 6 and is used to derive heats of fusion by the Figs. 6 to 9 probably have their origin in the varying
method described by Feilcbenfeld et al.[17]. As seen CaCI2 concentration in the trapped interstitial liquid
in Fig. 11 a correct value of the heat of fusion could as explained in the previous section. A lower limit
be obtained by this method only by way of imposing for the heat of fusion is, consequently, given by the
specific heat slopes adjusted to their assumed values upper limit of this variable found in a large number
of 1.4 and 2.2 J/g°C for the solid and liquid states, of cycles.
respectively. Cycle No. 106, Fig. 6, which ran on 7 January,
The plots shown in Figs. 6 to 9 exhibit existing 1986, is a candidate for exhibiting this upper limit,
variations in a number of characteristics, including both generally and in comparison with its close
enthalpy of fusion and crystallization, degree of sub- neighbors. The sample was inspected visually at times
cooling and subcooling time intervals. As to the en- before and after this cycle and appeared totally liquid
thalpy characteristics, their analysis has led to values when molten and near solid when frozen. A value of
for the heat of fusion which are 10% higher than re- 220 J/g is indicated in Fig. 6 as a lower limit for the
cently quoted by Lane[7]. Similar values were ob- heat of fusion of CaCI2"6H20 obtained in that cycle,
tained also by Feilchenfeld and co-workers[17], and in good agreement with the value obtained by Feil-
we feel confident that our measurements are accurate chenfeld and co-workers[17]. Many similar values
within about 2%. The values obtained depend to some were obtained in other cycles, both before and after
minor degree on the specific heats of the liquid and cycle No. 106, and cycle No. 542, shown in Fig. 9,
the solid CCH, assumed to be 2.1 and 1.4 J/g°C, presents one such example.
.~3z,
r~
~28,
N 2~,
20'
16
0 2b I00 120 I/:0I(~0180 200 2~ 2~0
TIME POINTS
2500(30
~ i n p u t : 1+.6x2.1÷10.1x1-¢,
/ -~1 Latenfheat° N usi°n : 231.-221'=k210/kkgJ/kg
..._1
I,--
,.z,
-50000
Fig. 6. Daily cycle No. 106 (7 Feb. 1986). Temperatures measured at point locations as approximately
indicated in the inset. Each "time point" represents 6 min. Cycle 106 was found to have the highest
measured heat of fusion: 220 J/g, corrected for excess water.
188 A. B RANDSTETrER
32
~3O
~. 2B
26
/Y/o+++ •++t\
22
20
18 I I "-'~- t I I' I I I I I I I I } •
0 5 10 25 3O 35 I~I ~+5 50 55 60 0 ½ 10 15 20 25
TIME POINTS' ' TIME POINTS
125000 2500O0
Specific heat estimates: :~
Solid: 7"0 + 5.0:1"~
--3 Liquid : 103 + S0: 2'IkJ/kg'C
)-
O-
.J
,(C
T I c:, T
k,-
i,J
-25000 100000
Fig. 7. Expanded-scale portions of the temperature and enthalpy curves of cycle 106, showing the sen-
sible-heat gain and loss patterns and derived specific heat values in the solid and liquid state.
250000
-S0000
P i g . 8. C y c l e nos. 41 and 42, o f the f+~st r e c o r d e d s u i t a b l e f o r p r e s e n t a t i o n . N o t e the l o n g f r e e z i n g t i m e
due to the smaller downward temperature excursion, which also is in part responsible for the relatively
lower heat o f fusion.
Thermal storage systems 189
~:28
TIME POINTS
250000 =,z-,,,- - - T . . . .
Sensible heat
/
.1 • 10.5 x 14,
\ = 37kJ/kg /
/
'~ Latent heat o~fusion: ZU,-37~=207kJ/kg /
-50000
Fig. 9. Two consecutive cycles on day 542 (March 20, 1987). Note (a) the variation in the degree and
time of subcooling (compared with cycle 41, Fig. 8); and (b) the heat of fusion values comparable to
those obtained a year earlier.
.~28 . . . . ~"~ - -
I I,~ ,,~,,,,P.,,~l ~ Starting on day NO. 578 Terminated on day NO. 581
0 20 &O 60 80 100 120 I~,0 160 180 200 220 2&O 260 280 300 320 3/+0 360
TIME POINTS
2SC)C)O0 {5 time points = 1 hour}
A
Z
W
-5(~X)0
Fig. 10. Enthalpy variation during very slow traverses through the crystallization region, after which
the PCM was visually inspected and no tetrahydrat¢ deposit was found.
190 A. B RANDSTETTER
2~0
this slope . . . . . _ _ ~ , ~ t - - ~ ' ~ F e d line
200 odjusted to F" ~" R2=0"53
2.2 ~ °c
A
t, .,/a_E :3.z
~160 /'~ ," AT
40
• slope adjusted tol.4~°C
- ....----, ~ "~ , . J , ,
20 22 2/+ 26 28 30 32 3t, 36
TEMPERATURE {°C)
Fig. 11. Enthalpy vs. temperature "scatter plot" for cycle 106 during heat-up--only minimum and
maximum temperature points are shown for each plotted enthalpy value.
ment and Demonstration Program which is administered by Change Material from Industrial Waste for Thermal
the Commonwealth Department of Resources and Energy. Storage, final report to National Energy Research, De-
Substantial contribution by Prof. S. Kaneff and helpful dis- velopment and Demonstration Council (NERDDC),
cussion with P. O. Carden and B. K. Selinger also are Grant No. 84/4203, May (1987).
gratefully acknowledged. 10. J. A. Dean (ed.), Lange's Handbook of Chemistry, llth
ed. p. 10. McGraw-Hill, New York (1973).
REFERENCES 11. J. W. Mellor, Inorganic and Theoretical Chemistry.
Vol. IIl, p. 703. Longmans, Green and Co., City (1946).
1. L. Pauling, General Chemistry 2nd edition, p. 338. W. 12. H. W. B. Roozeboom, Z. Phys. Chem. 4, 31 (1889).
H. Freeman and Co., San Francisco (1958). 13. G. A. Lane, Solar Energy 27, 73-75 (1981).
2. S. Furbo, Sunworld 6, (2) (1982). 14. M. Telkes, Thermal storage for solar heating and cool-
3. A. Abhat, Solar Energy 30, 313-332 (1983). ing. Proc. Workshop on Solar Energy Storage Sub-
4. H. Kirnura and J. Kai, Solar Energy 33, 49-55 (1984). systems for the Heating and Cooling of Buildings,
5. G. A. Lane (ed.), Solar Heat Storage : Latent Heat Charlotteville, Va., The University Press of Virginia,
Materials, Vol. I, p. 115. CRC Press Inc., Boca Ra- 17-23 (1975).
ton, Florida (1983). 15. A. Brandstetter, Heat storage in calcium chloride hexa-
6. P. Carden and A. Brandstetter, Modelling the nuclea- hydrate. Proc. Solar World Congress 1983. Pergamon
tion and subcooling behavior of calcium chloride hexa- Press, Sydney, 1845-1849 (1984).
hydrate phase change material, to be published. 16. A. Brandstetter, Phase change storage applications. Solar
7. G. A. Lane (ed.), Solar Heat Storage: Latent Heat 86 Proceedings 13-15 Nov. 1986, Adelaide, South
Materials, Vol. H, Technology, p. 14. CRC Press Inc, Australia, 244-255 (1986),
Boca Raton, Florida (1986). 17. H. Feilchenfeld, J. Fuchs, and S. Sarig, Solar Energy
8. B. Carlsson, H. Stymne and G. Wettermark, Solar En- 6, 779-784 (1984).
ergy 23, 343-350 (1979). 18. H. Feilchenfeld, J. Fuchs, F. Kahana, and S. Sarig,
9. A. Brandstetter and S. Kaneff, Evaluation of Phase Solar Energy 34, 199-201 (1985).