Water Chemistry for Reactor Tech
Water Chemistry for Reactor Tech
F RT ln K
the activities or fugacities
Thus, for the reaction
(3.1)
of
the
phases
H2O H2 + O2
the equilibrium constant is
at
equilibrium
can
be
calculated.
Ke
f O12/ 2 f H 2
f H 2O
For the very small degrees of dissociation of pure water. which occur at the temperatures of
interest, we write
f H 2 PH2 2 f O2 2PO2
where the fugacity is then equal to the partial pressure, for the gases. Water does not behave
as an ideal gas, and a correction is required. Define as f/P, me ratio of fugacity-to-partial
pressure for pure water.
Thus,
f H 2O PH2O
(3.2)
Broch2 has computed for saturated water vapor (and therefore the saturated liquid) up to the
critical point, and his data may be used in the calculations. He finds log (sat)-1 = 300P0.618,
with P in pounds per square inch.
Thus, Eq. (3.2) is written
PO2 2 PO2
Ke
(3.3)
Ke
PsatH2O
2
(3.4)
PsatH2O
or
PO22
3
and
K
PO2 e PsatH2O
2
(3.5)
for the dissociation of pure water. Although this treatment is not entirely exact, further
refinements are neither warranted nor possible with existing data on mixtures of steam and
gases.
For metal reactions, similar equations are derived. Thus, for the reaction
NiO = Ni + O2
Ke
a Ni f O2
a NiO
(3.6)
PO2
(3.7)
since the solid phases have unit activity, and the fugacity of the gas maybe taken equal to its
partial pressure.
Thus,
2
PO2 K e2
(3.8)
If the equilibrium partial pressure of oxygen from NiO [Eq. (3.8)] is greater than the
equilibrium partial pressure of 0 1 from pure water [Eq. (3.5)] at a given temperature, then Ni
would be stable in pure water at that temperature. This implies that the ratio of water-to-nickel
is very large; i.e., the composition does not change with a small increment of reaction. For
such comparisons, curves of the equilibrium partial pressure of oxygen for pure water, water
containing added hydrogen, water containing added oxygen, and various oxides of interest,
have been computed and are plotted in Fig. 3.1. The curves for heavy water (D2O) would be
little different from that of light water (H2O) and are not included on the chart.
Note from Fig. 3.1 that silver is thermodynamically stable with respect to pure water over the
whole temperature range covered, and, except at relatively low temperature ( <200 F), is
stable even in water containing dissolved oxygen at 0.5 cm3/kg.*Copper is not
thermodynamically stable with respect to pure water, but is stable with respect to water
containing 25 cm3/kg; Ni is not stable io pure water but is stable in the hydrogenated water;
Fe is not stable, even in the hydrogenated water.
Volumes of gases will always refer to 0C, 1 atm, unless otherwise noted.
effect of the oxygen predominates. Moreover, the chart assumes that the composition of the
liquid phase is not changed by the reactions proceeding toward equilibrium.
H+ + B
Proton Base
(3.9)
The acid and base are conjugate to one another. Since free protons cannot be expected to exist
in solution to any significant extent, the acidic or basic properties of a solute cannot be
realized unless the solvent itself possesses acid or basic properties, i.e., can accept or donate
protons. Thus, the general case for significant acid or basic solutions must be written
A1 B2 B1 A 2
(3.10)
where Al and Bl and A2 and B2 are the conjugate acid and bases of the solute and solvent,
respectively.
Water can act as both a proton acceptor and a proton donor (amphiprotic) and therefore its
solutions can show both acidic and basic properties. In heterogeneous water reactor systems,
interest is confined largely to pure water and basic solutions, and these will be discussed in
detail.
Water dissociates to solvated protons and hydroxyl ions
2H 2O H3O OH
(3.11)
with similar reaction for heavy water. The equilibrium for this reaction is written
aH O aOH
3
Kw
aH2 2O
(3.12)
For dilute solutions, up to 10-4 M, the activities are nearly equal to the molal concentrations m,
and the water activity is essentially constant. We may write, therefore.
4
mH+ mOH- = Kw
(3.13)
and the constant K.._will apply for all conditions consistent with the above qualifications.
Measurements of K"" have been made by a variety of methods. Agreement at low
temperatures is quite good, including measurements based on the conductance of water. Very
few measurements are available at temperatures> 200 F. Those of Noyes et al.3 are most
definitive. even though they are over 50 years old. Table 3.1 and Fig. 3.2 present three sets of
data for the molal ion product of pure
Table 3.1
Molal Ion Product of Water
Temp., F
Kw 1014
Jones4
Temp., F
Kw 1014
Noyes et al.3
Kw 1014
Harned and
Owen5
0.57
54.6
222
403
308
75
0.90
64.4
0.64
200
43
212
55
300
208
313
268
400
500
424
645
480
645
583
340
540
628
water including a compilation from various sources by Jones4 and the data of Noyes et al.3
Harned and Owen5 present three relationshipsfor the variation of Kw, with temperature, and
the values for their Eq. (15-3-7a), shown here as Eq. (3.14), are included in Table 3.1 and Fig.
3.2.
log K w
4470.99
6.0875 0.0170607T
T
(3.14)
The bases of interest in water reactor technology are KOH, LiOR, and NH4OH. NaOH is not
ordinarily utilized because of problems resulting from activation of 23Na to 24Na. NH,OH is a
weak base. LiOH has been shown experimentally6 to be associated at high concentrations, and
the same is probably true for KOH at high temperatures. At the normal concentrations of
interest (10-4 M), however, only the dissociation of NH4OH need be taken into account in
computing solution pH. The dissociation constants for NH4OH and LiOH are presented in
Table 3.2 and Fig. 3.3, respectively. Table 3.2 also presents KA for the HSO4- ion. The
calculated pH at high temperature of pure water, dilute solutions of ammonia, and strong
bases, are plotted in Fig. 3.4, where SB is molal concentration of added base. In Fig. 3.4, the
point for K2SO4 , at 560 F, indicates that a 0.5 10-4 M solution of K2SO4 is more alkaline at
that temperature than a solution of NH3 which has a pH at 75F of 10.0{11.1 ppm NH3).
Table.3.2
Dissociation Constants: KB-NH4OH and LiOH; KA-HSO4-: Molal Units
NH4OH
LiOH
HSO4Temp., F
KB, Wright5
KB 106 Jones4
KB 106 Wright6
KA 106, a
75
17.9
120
21.9
0.135
160
19.1
0.0731
200
15.0
16.1
0.0804
240
12.4
0.0502
880
280
10.2
0.0348
300
7.80
300
320
7.37
0.0500
360
5.18
0.0467
400
2.83
3.42
0.0502
42.0
440
2.19
0.0320
480
1.13
1.30
0.0305
10.0
520
0.84
0.0234
540
0.43
560
0.35
2.00
a
Data from W. L Marshall and E. V. Jones. Reactor Chemistry Division Annual Progress
Report For Period Ending Jan. 31, 1965, ORNL-3789, p. 148, Oak Ridge National
Laboratory; also J. Phys. Chem,. 70, 12(December 1966).
The dissociation equilibria are calculated from the following equations:
[M+] + [H+] = [OH-]
(Condition of Electroneutrality)
(3.15)
[H+][OH-] = Kw
(Dissociation of Water)
(3.16)
(3.17)
M OH K
MOH
(3.18)
(Dissociation of Base)
Fig. 3,4. pH of solutions of strong bases (SB), NH3, and pure water.
A value of pH is assumed, defining [H+]. From this, one calculates [OH-] from Eq. (3.16), and
[M+] then from Eq. (3.15). For completely dissociated bases, [M] = [B]. For incompletely
dissociated bases, [MOH] is then calculated by substituting the values of [M+] + [OH-] in Eq.
(3.18); [B] is then calculated from Eq. (3.17). For univalent weak bases,
Kw
Kw
H 1
H
K B H
(3.19)
The calculations involve only simple arithmetic and are used to generate data for cross
plotting to give the results in more convenient form. More detailed calculations have been
reported (Meek, 13 Chap. 6).
concentration of metals in equilibrium with corroding Ni-CrFe alloy 600. His results, as
tabulated by Cheng,8 are presented in Table 3.3. An electropolished Incoloy-600 autoclave
was operated for 1500 h at 550 F (288C) at the starting coolant conditions shown in the lefthand column. For the data of the other temperatures shown, the autoclave was equilibriated
for one day at the temperature of interest. An additional test was run (see Table 3.4) to
determine the effect of short-term changes in water chemistry after previous exposure for
1500 h to a different chemistry, all at 340 F. The composition of the Incoloy-600 used for the
autoclaves is given in Table 3.5. Jones7 found that the "solubility" decreases with time, and
generally reaches a fairly steady state after 1500 h. The various metals are not found in the
water in proportion to their composition in the alloy, and the solubilities of the various
elements vary with temperature and solution pH. The short-term tests, Table 3.4, indicate that
the measured solubility is primarily dependent on the final water chemistry, not that of the
previous period of exposure. Although these data are of limited validity as solubility data,
they provide some idea of the effect of alkaline additives and of temperature on the solubility
of the various metals, which should be useful in considering ionic transport of corrosion
products.
The data as tabulated by Cheng8 represent the temperature effect quite well, but include a
great deal of scatter for the effect of pH. A better representation of the solution composition at
550 F is presented in Fig. 3.5 from the original data7 for samples after about 1000 h of
constant-temperature operation. The pH is that of the solution at 550 F.
Table 3.3
Solubilities of Incoloy-600 Corrosion Producs, Corrosion, Ref. 8)
pH
8.5
8.5
8.5
8.5
8.5
Element
ppb
Iron
Nickel
Chromium
Manganese
Cobalt
9.5
9.5
9.5
9.5
9.5
75(24)
1.0
14.4
<0.1
0.7
<0.1
16.3
200(93)
9.0
10.9
<0.1
2.6
0.3
22.9
Iron
Nickel
Chromium
Manganese
Cobalt
2.2
8.7
<0.1
0.4
0.1
11.5
2.9
4.2
<0.1
0.6
0.1
7.9
3.3
1.6
<0.1
0.7
<0.1
5.8
3.5
0.8
<0.1
0.6
<0.1
5.1
3.5
0.8
<0.1
0.3
<0.1
4.8
4.1
0.6
<0.1
0.3
<0.1
5.2
10.0
10.0
10.0
10.0
10.0
Iron
Nickel
Chromium
Manganese
Cobalt
0.8
0.2
<0.1
0.7
<0.1
1.9
1.1
0.4
<0.1
0.7
0.1
2.4
1.6
0.5
<0.1
0.8
0.1
3.1
1.0
0.7
<0.1
0.8
0.1
2.7
1.1
0.4
<0.1
0.7
0.1
2.4
2.2
0.3
<0.1
0.8
0.1
2.6
10.5
10.5
10.5
10.5
10.5
Iron
Nickel
Chromium
Manganese
Cobalt
1.8
0.6
0.1
0.2
<0.1
2.8
2.4
0.8
0.1
0.2
<0.1
3.6
2.7
0.9
0.1
0.4
<0.1
4.2
2.8
0.8
0.1
0.3
<0.1
4.1
3.0
0.7
0.1
0.2
<0.1
4.1
3.8
0.7
0.1
0.3
<0.1
5.0
10
650/343)
1.5
1.0
<0.1
1.3
0.2
4.1
Water
Lithiated
and
Hydrogen
Saturated at
Room
Temperature
pH
10.5
10.5
10.5
10.5
10.5
Element
ppb
Iron
Nickel
Chromium
Manganese
Cobalt
75(24)
0.6
3.6
<0.1
0.6
<0.1
5.0
200(93)
1.2
0.8
<0.1
0.6
<0.1
2.8
650/343)
1.7
0.3
<0.1
0.2
<0.1
2.4
Table 3.4
Comparison of Solubilities of Incoloy-600 Corro,Ion Products
(Corrosion, Ref. 8)
[Hydrogen saturated and ammoniated water at room temperature, one-day run after
~1500 h al 340 F (171 C) Static water]
Solubility ppb Fe
pH 8.5
pH 10.0
21.4
1.8
27.7
2.0
Solubility ppb Ni
pH 8.5
pH 10.0
11.0
0.5
4.2
0.4
Table 3.5
Composition of Incoloy-600*
Element
Nickel
77.49
Chromium
15.14
Iron
6.78
Silicon
0.18
Manganese
0.22
Copper
0.10
Carbon
0.030
Sulfur
0.007
*Trademark, The International Nickel Co.
Although some scatter is still evident, the trends are reasonably clear. The solubilities of
nickel and manganese decrease substantially with increasing pH. The solubilities of iron and
chromium decrease with increasing pH but not as sharply as nickel and manganese. The
selective effect of the corrosion process is evident. The solution is enriched in iron and
manganese, and depleted in chromium and nickel, compared to the base metal.
A more classical approach to solubility measurements was taken by Mingulina et al.,9 in the
case of cobalt, which because of its longlived activity (60Co) is of extreme importance in
water reactor technology. The solubility of cobalt was studied in a boiling water apparatus at
185 atm (2720 psia, ~680F) as a function of room temperature pH and chloride content of the
water. The results of this investigation are shown in Fig. 3.6.
11
12
1
1
4
Fe3O 4 (s) (2 b)H H 2 (g) Fe(OH)(2b b) ( b)H 2O
3
3
3
They have interpreted their data in terms of three values of b corresponding to three soluble
forms for the iron; namely,
Fe2+
Fe(OH)+
Fe(OH)3-
b=0
b=1
b=3
13
Fe(OH)
H P
2 b
2 b
b
1
H2
log K
A / T B(ln T 1) D
2.303R
14
1000
where C is the concentration of the solute in equivalents per liter. In all cases, the observed
conductance is corrected by subtracting the conductance of the pure solvent. Experimentally, it is
observed that the equivalent conductance is a function of concentration. The value of at infinite
dilution, termed 0 is a unique property of the solute; it is an additive property of the conductance at
infinite dilution of the individual ions (Kohlrausch's law of the independent migration of ions). Table
15
3.6 lists the ion conductances at infinite dilution at 25C (o+, o-) of a number of ions of interest, from
which 0 for a variety of salts can be obtained.
Table 3.6
Ion Conductance. at Infinite Dilution at 25 C, ohms-1, cm2 equiv-1
Cation
o+
Anion
o-
H+
Tl+
K+
NH4+
Na+
Li+
Ba2+
Ca2+
Sr2+
Mg2+
349.82
7.4.4
73.52
73.4
50.11
38.69
63.64
59.50
59.96
53.06
OHBrClNO3HCO3SO42-
197.6
78.4
76.8
71.44
44.5
79.8
0 0 C
(3.19)
Table 3.7 presents values of and for water as a function of temperature. 0 must be determined
experimentally and is a function of both pressure and temperature. The pressure dependence at
ordinary temperatures is quite small up to 5000 psi and need not be considered further for reactor
applications.
Table 3.7
Debye-Hckel-Onsager Conductivity Equation Parameters
0.2261
0.2289
0.2406
0.2549
0.2723
0.3016.
0.3139.
0.3318
0.4254
0.5272
0.6013
0.7269
Temperature, F
64.4
77.0
122.
167
212
262.4
284
312.8
424.4
500
537.8
582.8
50.50
60.10
99.75
146.5
201.0
266.6
298.0
340.0
523.5
650.8
737.3
870.0
Johnston11 showed that the variation of 0 with temperature was fairly well represented by the
equation
0 k n
(3.20)
16
where k and n are empirical constants and the solution (solvent) viscosity. Values for n for
some electrolytes in water are
NaCl
NH4OH
KCl
LiOH
= 0.94
= 0.71
= 0.86
= 0.75
In general, electrolytes with a high value of 0 at room temperature have lower values of n
than electrolytes with lower room temperature values of 0. Figure 3.7 shows the variation of
0 for a number of electrolytes over the temperature range of interest, including H-OH. These
data were obtained from the summary of Wright et al.6 and were extrapolated in part by using
the 0- relationship. Shown also is the calculated specific conductance of H2O. Quist ct al.12
have, in a similar fashion, prepared a compilation of limiting ionic conductance at high
temperatures, presented in Table 3.8.
Table 3.8
Limiting Equivalent Conductances of Several Ions
at Temperatures to 400 C (Ref. 12)
ohms-1 cm2 equiv-1
Temperature, C
H+
Li+
Na+
K+
NH4+
OH-
Cl-
100 (1 atm)
200 (0.865 g/cm3)
300 (0.7125 g/cm3)
400 (0.8 g/cm3)
364
824
894
945
156
329
562
151
304
459
440
195
364
504
455
206
394
579
447
701
821
211
391
561
520
A useful application of these data is the calculation of the diffusion coefficient for the
infinitely dilute solution. From the Nernst Formula,
Z Z2
D (cm 2 /sec) 8.93 10-10T 1
Z1 Z2
0 0
0 0
For the dilute solutions normally of interest in reactor technology, the difference between
and 0 is not ordinarily greater than about 1.5 %, as expressed by the Debye-Hckel-Onsager
Equation. However, at the higher temperatures, ion-pair formation can be significant for some
electrolytes. Bjerrum13 provided a method of analysis of the formulation of ion pairs in
completely ionized electrolytes. The association is expressed as an equivalent dissociation
constant KDI defined by
K DI
C A
C A
(3.21)
Experimental data are available for a number of electrolytes of interest. These are summarized
in Table 3.9 taken from Wright et al.6
17
Table 3.9
Ion Association-KDI (Ref. 6)
Electrolyte
Temperature, C
NaCl
KCl
HCk
NaOH
281
306
281
306
260
306
218
KDI
0.22
0.093
0.49
0.12
0.12
0.035
0.55
P1 P10 N1
(3.22)
(3.23)
18
and
P10 P1
N2
P10
(3.24)
That is, the relative vapor-pressure lowering is proportional to the mol fraction of solute. For
completely dissociated solutes in water, the mol fraction of ions must be used for N2 in the
ideal case, but departures from ideality are to be expected.
Considering evaporation as an isobaric process, it will be limited by concentration so that both
the thermal (temperature) and material (vapor pressure) exchanges cease. Thus, the required
conditions for equilibrium are that the solution concentrates until its vapor pressure at the
prescribed temperature P1T is equal to the total applied pressure , assuming that T is greater
than the saturation temperature of the pure solvent at the pressure . The solution vaporpressure curve is limited by solubility or critical phenomena, of which the former is of
primary concern in current water reactors.
Table 3.10
Vapor-Pressure Depression Data for Lithium and Sodium Hydroxide Solutions of Equivalent Concentrationsa
LiOH
NaOH
LiOH
NaOH
LiOH
NaOH
Concentration
Grams per 100 Grams H2O
5
5
8.4
10
10
15.6
15b
15
25
Weight Percent
4.76
4.76
7.83
8.59
8.59
13.5
13.04
13.04
20.0
Molality
2.09
1.25
2.09
3.92
2.35
3.92
6.25
3.75
6.25
Vapor-Pressure Depression (P), psi
248 F
1
1.2
1.95
1
2.25
4
4
3.6
6.7
320 F
4.5
3.6
6.3
3
7.15
11.5
7
11.2
18.8
392 F
12.5
9
14.5
13
16
27
18
26
44
482 F
35.5
25
38
34
42.5
69
50
67
107
572 F
73
52
83
74
93
147
90
144
226
662 F
120
98
158
140
177
280
137
275
433
a
NaOH data calculated from vapor-pressure data in the International Critical Tables; LiOH data calculated from smooth curves in Figs. 3, 4 and
5 of Report BMI-1329.18
b
Saturated, approximately.
20
22
Fig. 3.11. pH of sodium orthophosphate solutions having various mol ratios of sodium-tophosphate (Combustion. Ref. 22).
Even the use of the modified form of coordinated phosphate control does not ensure the
complete absence of free caustic. Concentrated phosphate solutions can react with iron oxide
to form sodium iron phosphate, leaving behind a more alkaline solution than that originally
present.
3-5.1.1 NH3
Because of its volatility, NH3 is attractive as a source of alkalinity in reactor and boiler
technology. In a simple (single-stage) evaporation process the limiting concentration ratio will
not exceed the reciprocal of the distribution constant, vapor to liquid. NH3 is quite stable,
thermally, at the temperature of concern in reactor technology, but is subject to radiolysis (see
Chap. 4). Morpholine and cyclohexylamine are also used in conventional boiler technology as
a source of basicity in the condensing portion of the steam cycle. Jones23 has reported
experimental determinations of the distribution coefficient for NH3 at various concentrations
and high temperatures.
The equilibrium is between un-ionized ammonia in the liquid and NH3 in the vapor. Thus,
[NH3]l [NH3]v
(3.25)
23
and
NH3 v
NH3 l
KD
The relationship between [NH3]l, identical to [NH4OH], and NH3, the total ammonia
concentration in solution, is readily derived by the method of Eqs. (3.15) to (3.18). Write
[NH4+] as NH3 and [NH4OH] as (1 -1) NH3. The relationship is
1
K
1 B H
KW
(3.26)
NH3 v
K
1 - NH3 D
(3.27)
NH3 v
NH3
K D 1 - K D'
(3.28)
Figure 3.12 shows KD for NH3 as a function of temperature for water at saturation pressure.
3-5.1.2 Iodine
Radioisotopes of iodine of considerable biological significance are formed in nuclear fission
in high yield. The chemical behavior of iodine, and particularly its volatility, are matters of
considerable practical importance in nuclear power plant design and operation. The chemistry
of iodine-water systems is quite complex. Thus, iodine reacts in water (hydrolyzes) to form
hydriodic and hypoiodous acid,
I2 + H2O HI + HIO
(3.29)
Hypoiodous acid, in turn, can disproportionate to iodic acid and iodine. Thus,
5HIO HIO3 + 2I2 + 2H2O
(3.30)
The hydrolysis reaction increases with temperature and pH, and is inversely proportional to
the concentration of iodine. Dilute solutions, less than -10-5 M, are extensively hydrolyzed
even in neutral solution. Volatility is contributed by both I2 and HIO (a weak acid), since the
degree of dissociation of HIO will decrease with increasing pressure. The apparent KD
defined as
HIO I 2 v
HIO HI I 2 l
is a complex function of temperature and pH.
Styrikovich et al.24 have published the results of a series of determinations of the distribution
of iodine compounds between steam and water, as a function of pressure (temperature) and
low temperature pH (pH0), for saturated steam.
24
Fig. 3.12. Distribution coefficient KD of undissociated ammonia, saturated water, and steam
vs temperature (J. Phys. Chem., Ref. 23).
The measurements were made over a range of total iodine concentrations of 1 to 100 mg/kg,
but are reported for concentrations of 10-4 M (25.4 ppm). Figure 3.13 shows the effect of
pressure and pH, on the apparent distribution coefficient KD at low and intermediate
pressures. The curves reflect the hydrolysis with increasing pH, decreasing the volatilization
of the I2 and increasing volatility at higher pressures in the high pH0 region, reflecting the
increased volatilization from HIO. Figure 3.14 shows the apparent distribution coefficient KD
for HIO at 425 and 1000 psi. Here, the coefficients show maxima at a value of pH0 of about
8.5, reflecting the amphoteric nature of HIO, which is the major species present. The true
distribution coefficient of HIO would be twice the value of KD at 1000 psi, assuming that
hydrolysis of I2 was practically complete, and dissociation of HIO negligible at that
temperature.
Note from Fig. 3.13 that at low pressure, such as in evaporators, and low pH, neutral or lower,
the apparent distribution coefficient of iodine is extremely high, of the order of 610-3 or
greater, and that this is reduced to about 410-5 at a pH 10. Ordinarily, iodine derived from
high-temperature water reactors is present largely as iodide, which will have a very low
volatility except at very low pH values.
25
Fig. 3.14. KD values of hypoiodous acid as a function of pH0, 71 kg/cm2 pressure (o) and at
30 kg/cm2 pressure (0) (Ref. 24).
26
equilibrium with the saturated solution A-B. The isobars G-H are the compositions in
equilibrium with solutions between the saturation line A-B and the critical line C-D in the
transcritical region.
28
Fig. 3.17. Phase equilibrium diagrams for system NaCl-H2O (Ref. 33)
29
(3.31)
representing a solute molecule reacting with m water molecules to form a complex soluble in
the vapor phase. The equilibrium constant for this reaction is defined by the equations
F 0 H 0 TS 0
(3.32a)
F 0 RT ln Ke
(3.32b)
Ke
a xmH 2 O v
(3.33c)
a x amH 2O v
30
If we define the activity of the vapor as proportional to its density , that of the dissolved
solute complex (x mH2O) as equal to the concentration of the solute in the vapor phase, and
that of the pure solid phase as unity, then
ln Cv m ln v
H v0
Constant
R
(3.33)
ax
a xmH 2O
(3.34)
K e aHm2O
Since the two phases are in equilibrium, ax must be the same. Thus,
a xmH 2O a xmH 2O
K e aHm O K e aHm O
2
2
v
l
(3.35)
As above, assuming m is the same for both phases, we can then write
ln
0
Cv
H l ,v S 0 (l , v)
m ln v
Cl
l
RT
R
(3.36a)
At the critical point, the enthalphy and entropy terms vanish; the compositions must be equal.
For a restricted range of temperatures and densities, not too far from the critical, we can write
Cv v
Cl l
(3.36b)
This equation has been used by Styrikovich and Martynova25 to correlate the distribution of
solutes between steam and water, with application to boiling water reactors. Figure 3.19 from
Styrikovich and Martynova25 indicates that, for many solutes, the simplified equation is
obeyed over a substantial range of pressures. In general, those materials that are weak
electrolytes in water -Al2O3, B2O3, SiO2- have high distribution coefficients and low values of
m of the order of 1 to 2. Strong electrolytes, such as NaOH and NaCl, have low distribution
coefficients and high values of m of the order of 4.
At high concentrations of solute there are effects due to the amphoteric nature of some of the
weak electrolytes, and departures from strong electrolyte behavior.
31
Fig. 3.19. KD, distribution of solutes between water and steam (Ref. 25)
The effect of pH is illustrated by the data for Co2+ at 680F as shown in Fig. 3.20. The
volatility data reflect the solubility behavior as shown in Fig. 3.6. Cobalt occurs only in the
form Of Co(OH)2 at high pH. The solubility is low, approximately 5 g/kg as Co2+, but KD =
0.3. For larger amounts of Co2+ in the water, the concentration in the steam will be
approximately constant at 1.5 g/kg and the value of KD will decrease. At lower pH, the total
concentration of the Co2+ dissolved in the water is able to increase greatly as a result of the
higher solubility of the basic cobalt salts formed in this range, of the type CoOHCl.
Amphoteric compounds, such as Al2O3, show maxima in the relationship of KD as a function
of pH.
If solubility of the ionic form in steam is negligible, then we can write for electrolytes
Cv
undissociated K D
Cl
(3.37)
32
-4
Apparent KD
1.3 10
4.3 10-4
(Degree of Dissociation of NaCl)
0.923
0.62
1-
0.077
0.38
Undissociated Na, ppm
1.88
87.5
-3
True KD
1.69 10
1.14 10-3
The other region of interest, that of solute solubility in steam from the saturated liquid to the
superheat and transcritical areas, is most readily examined over the whole temperature range
for constant pressure conditions. For supercritical pressures, the isobar represents a
continuum, with rapid but continuous variations in the pseudocritical region. Styrikovich et
al.28 have presented a set of solubility curves for some materials of interest at pressures of 255
and 300 atm. Figure 3.21 shows the solubility plotted as a function of temperatures for a
variety of materials.
33
Figure 3.22 is a plot of the solubilities of some of the more important solutes at 255 atm, with
fluid enthalpy as the independent variable. Comparison of Figs. 3.21 and 3.22 demonstrates
the advantage of the enthalpy plot in spreading out the plots in the pseudocritical region. The
change in solubility with temperature in this region closely parallels the variation of density
with temperature, as shown on Fig. 3.23.
The data for iron as Fe3O4 and copper as CuO are of major interest in these figures. Fe3O4
data are not specifically referenced. The copper data are those of Pocock and Stewart29 and
Deeva.30 Much of the Russian data were obtained in an experimental once-through boiler in
which the solute in the water was determined at the various parts of the cycle for various feed
concentrations. The feed concentration above which the point concentration did not increase
was taken as the solubility at that point. In the ferritic system employed, no effect of iron feed
was found at the high temperature parts of the cycle. At the mass velocities and heat transfer
rates employed, the system was therefore in equilibrium with magnetite in the absence of
additions.
The data in Figs. 3.21 and 3.22 permit an assessment of equilibrium processes in the heat
source region. The applicability to processes at high mass and heat transfer rates has been
questioned.28
35
36
However, tests in a supercritical loop16 at high heat transfer rates (300,000 Btu (h ft2)] have
indicated that chemical processes proceed according to the slate of the bulk fluid, not
according to the state corresponding to the surface temperature.
Evaluation of processes in the turbine require, in addition to the data shown in Figs. 3.20 and
3.21, data in the superheat region as shown for NaCl and SiO2 in Figs. 3.17 and 3.18. Similar
data for copper and its oxides are available only at a very limited range of temperatures and
pressures. A compilation of data has been presented by Styrikovich et al.28 and is shown in
Fig. 3.24 for the state line of a particular supercritical turbine cycle. The pressure and
temperalure curves for this turbine cycle are included in the diagram.
3-6.2 Entrainment
At low pressures and moderate temperatures, solute volatility is a minor source of transport
compared to physical carryover or entrainment. Consider a droplet of water, containing
sodium chloride, carried from the steam-generating to the superheating section of a reactor.
From Fig. 3.10 it is evident that solutions of NaCl are stable up to 590 F at 1000 psi.
The potential therefore exists for accelerated attack by chloride solutions, and carryover and
chlorides must be maintained at low values. To prevent the formation of concentrated
solutions, it would be necessary to maintain all surfaces above 590 F at all times. Droplets
carried in the steam without being deposited will dry out and the salt will ultimately volatilize
into the steam if the bulk concentration is less than the minimum solubility for the
superheater.
As previously noted, except for LiOH, caustic solutions will not dry out. These solutions are
aggressive to almost all desirable materials of construction so that free caustic in the steam
generating source should be avoided in superheat systems.
37
38