The International Association For The Properties of Water and Steam
The International Association For The Properties of Water and Steam
The International Association For The Properties of Water and Steam
Kyoto, Japan
September 2017
Please cite as: International Association for the Properties of Water and Steam, IAPWS R16-17, Release
on the IAPWS Formulation 2017 for the Thermodynamic Properties of Heavy Water (2017)
This release replaces the corresponding revised release of 2005 (IAPWS R3-84(2005), which replaced the
original release of 1984), and contains 17 pages, including this cover page.
This Release has been authorized by the International Association for the Properties of Water and
Steam (IAPWS) at its meeting in Kyoto, Japan, 27 August to 1 September, 2017. The members of
IAPWS are: Britain and Ireland, Canada, the Czech Republic, Germany, Japan, New Zealand, Russia,
Scandinavia (Denmark, Finland, Norway, Sweden), and the United States, and associate members
Argentina and Brazil, Australia, China, Egypt, France, Greece, Italy, and Switzerland. The President at
the time of adoption of this document was Prof. Hans-Joachim Kretzschmar of Germany.
Summary
The formulation provided in this Release is recommended for calculating thermodynamic
properties of heavy water, which IAPWS defines as water whose hydrogen atoms are entirely the
deuterium isotope (2H or D) and whose oxygen atoms have the isotopic composition of ordinary water
[1]. Further details about the formulation can be found in the article “A New Reference Equation of State
for Heavy Water” by S. Herrig et al. [2]. This formulation provides the most accurate representation
available at the time this Release was prepared for the thermodynamic properties of the fluid phases of
heavy water over a wide range of conditions.
Further information about this Release and other documents issued by IAPWS can be obtained
from the Executive Secretary of IAPWS (Dr. R.B. Dooley, [email protected]) or from
http://www.iapws.org.
2
Contents
1 Nomenclature 2
2 Reference Constants 3
3 The Formulation 3
4 Relations of Thermodynamic Properties to the Dimensionless Helmholtz Free Energy 5
5 Range of Validity 9
6 Melting and Sublimation Curves 10
7 Estimates of Uncertainty 11
8 Computer-Program Verification 11
9 References 13
1 Nomenclature
2 Reference Constants
Tc = 643.847 K (1)
c = 17.775 55 mol dm3 (2)
R = 8.314 459 8 J mol1 K1 (3)
M = 20.027 508 g mol1 (4)
The numerical value for the critical temperature Tc is identical to that given in the IAPWS
release on the critical parameters of ordinary and heavy water [3]. The value for the critical
density c is the massic critical density given in [3] converted to molar density by the molar mass
M, which is that given in [1]. The value of the molar gas constant R is that given by the 2014
CODATA evaluation [4].
Due to the use of the molar gas constant, Eq. (5) corresponds to a mole-based
formulation. In order to convert values of molar properties to mass-based properties, a molar
mass must be used. In general, this should be the value of M given by Eq. (4). In some special
cases (such as heavy water that is enriched in heavy oxygen isotopes), specific information may
be available that the isotopic composition differs from the definition of heavy water given in [1];
in such cases it may be appropriate to use a different molar mass corresponding to that different
isotopic composition.
3 The Formulation
The formulation is a fundamental equation for the molar Helmholtz free energy f. This
equation is expressed in dimensionless form, = f/(RT), and is separated into two parts, an ideal-
gas part and a residual part r, so that:
f T ,
, o , r , , (5)
RT
where = Tc/T and = /c with Tc , c , and R given by Eqs. (1), (2), and (3).
The ideal-gas part of the dimensionless Helmholtz free energy is obtained from a fit to
the high-accuracy calculations by Simkó et al. [5] of the isobaric heat capacity in the ideal-gas
state. It is given by:
4
o , a1 a2 ln c0 1 ln vi ln 1 exp ui Tc . (6)
i 1
4
Table 1. Numerical values of the coefficients and parameters of the ideal-gas part
of the dimensionless Helmholtz free energy, Eq. (6)
Parameter Value Parameter Value
c0 4.0
a1 8.678 502 998 052 1
a2 6.963 637 340 437 6
1 0.106 33 × 101 u1 308 K
2 0.997 87 u2 1695 K
3 0.214 83 × 101 u3 3949 K
4 0.354 9 u4 10317 K
The form of the residual part r of the dimensionless Helmholtz free energy is:
6 12
r , ni d t ni d t exp l
i i i i i
i 1 i 7
24
(7)
ni exp i i i i .
di ti 2 2
i 13
All parameters (coefficients ni, temperature exponents ti, density exponents di and li, and
parameters of the Gaussian bell-shaped terms i, i, i, and i) are listed in Table 2.
The IAPWS reference-state convention is that the internal energy and the entropy of the
saturated liquid at the triple point are set equal to zero. Thus, at the triple-point temperature
Tt = 276.969 K [6],
u t 0 , s t 0. (8)
In order to meet this condition, the coefficients a1 and a2 in Eq. (6) have been adjusted
accordingly. Because of differences in computer codes (such as convergence tolerances), a
particular computer code may produce values that differ slightly from the zeros in Eq. (8). If a
user wants to reproduce Eq. (8) more precisely, the constants a1 and a2 in Eq. (6) may be
readjusted by imposing the condition u t 0, s t 0 with the desired accuracy.
5
All thermodynamic properties can be derived from Eq. (5) by using the appropriate
combinations of the ideal-gas part , Eq. (6), and the residual part r, Eq. (7), of the
dimensionless Helmholtz free energy and their derivatives. Relations between thermodynamic
properties and and r and their derivatives are summarized in Table 3. All required
derivatives of the ideal-gas part and of the residual part of the Helmholtz free energy are given in
Table 4 and Table 5, respectively.
Besides the single-phase region, the formulation also describes vapor-liquid equilibria.
For a given temperature T, solving the phase-equilibrium equations (see Table 3) yields the
saturation properties p, , and at T. With these, all other saturation properties can be
derived from Eq. (5). In this way, the properties calculated on the saturation curve are
thermodynamically consistent with the properties of the single-phase region.
6
Table 3. Relations of thermodynamic properties to the ideal-gas part o and the residual part
a
r of the dimensionless Helmholtz free energy and their derivatives
Property Relation
Pressure p ,
1 r
p 2 f / T RT
Molar internal energy u ,
u f T f /T RT
o r
Molar entropy s ,
s f /T
R
o r o r
Molar enthalpy h ,
h f T f /T f / T RT
1 o r r
Molar isochoric heat capacity cv ,
cv u /T R
2
o r
Molar isobaric heat capacity 2
c p , 1 r
r
c p h /T p
R
2
o
r
1 2r 2
r
2
Speed of sound 1 r
r
w2 , r 2 r
w M 1 2 p/
1/2 1 2
s RT M 2
o r
r 2
r r
Joule-Thomson coefficient
R
T / p h
1 r r 2 o r 1 2r 2r
2
p p
Phase-equilibrium condition 1 r , ; 1 r ,
(Maxwell criterion) RT RT
p 1 1 r
ln , ,
r
RT
r 2 r r 2 r 2 r 2 o
a r , r , r , r , r , o , o .
o
2 2 2
7
Table 4. The ideal-gas part o of the dimensionless Helmholtz free energy and its derivativesa
4
o = ln + a1 + a2 + co 1 ln + vi ln 1 eu T
i c
i 1
o = 1/ + 0 + 0 + 0 + 0
o
= 1/ 2 + 0 + 0 + 0 + 0
4
1
o
= 0 + 0 + a2 + co 1 / + vi ui Tc ui Tc
1
i 1 1 e
e i
4 u Tc
vi ui Tc
2
o
= 0 + 0 + 0 – co 1 / 2
–
i 1 1 e ui Tc
2
o
= 0 + 0 + 0 + 0 + 0
a o , o , o , o , o .
o 2 o o 2 o 2 o
2 2
8
Table 5. The residual part r of the dimensionless Helmholtz free energy and its derivativesa
6 12 24 2 2
r ni di ti ni di ti e ni di ti e i i i i
li
i 1 i 7 i 13
6 12 24 2 2
r ni di di 1 ti ni di 1 ti di li li e ni di 1 ti e i i i i
li
di 2i i
i 1 i 7 i 13
6 12 24 2 2
li i i i i
r
ni d i d i 1 di 2 ti ni di 2 ti d i li li d i 1 li li li2 li e ni ti e
i 1
i 7
i 13
6 12 24 2 2
ti
r ni ti di ti 1 ni ti di ti 1e ni di ti ei i i i
li
2 i i
i 1 i 7 i 13
6 12 24 2 2 t
2
r ni ti ti 1 di ti 2 ni ti ti 1 di ti 2e ni di ti e i i i i
li t
i 2 i i i2 2 i
i 1 i 7 i 13
6 12 24 2 2
di ti
ni di ti di 1 ti 1 ni ti di 1 ti 1 di li li e ni di ti e i i
li i i
2i i 2 i i
r
i 1 i 7 i 13
______________________________________________________________________________________________________________
r 2 r r 2 r 2 r
a
r , r
2 ,
r
, r
2 , r
.
9
5 Range of Validity
IAPWS has tested the formulation and endorses its validity in the following way:
(1) The formulation is valid in the entire stable fluid region of D2O from the melting-pressure
curve (see Section 6) to 825 K at pressures from the triple-point pressure to 1200 MPa; the
lowest temperature on the melting-pressure curve is 254.415 K (at 222.41 MPa). At
pressures below the triple-point pressure and temperatures above 100 K, it is also valid for
the stable vapor at temperatures from the sublimation curve1 up to 825 K.
In this entire region, Eq. (5) represents the experimental data available at the time the
Release was prepared (except for very few data points) to within their uncertainties.
Although Eq. (5) is also in satisfactory agreement with the limited experimental data in the
critical region, the equation, like any analytical Helmholtz energy formulation, is not able to
reproduce the nonclassical values of critical exponents in the immediate vicinity of the
critical point.
(2) In the stable fluid region, the formulation can also be extrapolated beyond the limits given
under item (1). Tests show that Eq. (5) behaves reasonably when extrapolated to pressures
up to about 100 GPa and temperatures up to about 5000 K. This holds at least for the density
and enthalpy of undissociated D2O.
(3) The formulation behaves reasonably when extrapolated into the metastable regions.
Equation (5) is in general agreement, in most cases within experimental uncertainties, with
the currently available experimental data of the subcooled liquid (metastable with respect to
solid). No reliable experimental data are available for either the superheated liquid
(metastable with respect to vapor) or the supercooled vapor (metastable with respect to
liquid), but the behavior of Eq. (5) when extrapolated into these regions is physically
reasonable.
1
Note that the equation for the sublimation curve in Section 6 should not be extrapolated below 180 K.
10
While IAPWS has not made a official recommendation for the melting and sublimation
curves of heavy water, it is necessary to have a representation of these curves in order to
compute the boundaries of the range of validity as given in Section 5. For that purpose, the
correlations given in [2] are adopted, and summarized here.
Melting-pressure equation for D2O ice Ih (temperature range from 276.969 K to 254.415 K):
Melting-pressure equation for D2O ice III (temperature range from 254.415 K to 258.661 K):
Melting-pressure equation for D2O ice V (temperature range from 258.661 K to 275.748 K):
Melting-pressure equation for D2O ice VI (temperature range from 275.748 K to 315 K):
Sublimation-pressure equation for D2O ice Ih (temperature range from 210 K to 276.969 K):
p
ln subl 0.140 626 9 102 1 1.47 0.337 529 1 102 1 1.29 , (13)
pn
with = T/Tn, Tn = 276.969 K, and pn = 0.000 661 48 MPa. For checking computer imple-
mentations, at 245 K Eq. (13) gives a sublimation pressure of 0.326 856 321 × 104 MPa.
Equation (13) extrapolates reasonably down to at least 180 K.
11
7 Estimates of Uncertainty
Estimates have been made of the uncertainty of the density, speed of sound, isobaric heat
capacity, and vapor-liquid saturation pressure and coexisting densities when calculated from the
formulation, Eq. (5). These estimates were derived from comparisons with the various sets of
experimental data together with the judgment of the Working Group on Thermophysical
Properties of Water and Steam of IAPWS. IAPWS considers these as estimates of the expanded
uncertainties with coverage factor k = 2, roughly corresponding to a 95 % confidence interval.
For the single-phase region, these values are indicated in Figs. 1–3, which give the
estimated uncertainties in various regions. With regard to the uncertainty for the speed of sound
and the specific isobaric heat capacity, see Figs. 2 and 3, it should be noted that the uncertainties
for these properties increase drastically when approaching the critical point. The statement “no
definitive uncertainty estimates possible” for some regions is based on the lack of experimental
data in the region. However, because the formulation was validated against experimental data
and with regard to correct physical behavior over its entire range of validity, the results in these
regions should be physically reasonable.
As noted in the captions of Figs. 1–3, the uncertainties of properties for the vapor at low
pressures become small. This is because the low-density vapor approaches the ideal-gas limit;
the ideal-gas behavior is known essentially exactly for the density and to a relatively small
uncertainty for the speed of sound and isobaric heat capacity.
For the saturation properties, the estimates of the expanded uncertainties of vapor
pressure, saturated liquid density, and saturated vapor density are shown in Fig. 4.
8 Computer-Program Verification
To assist the user in computer-program verification, three tables with test values are
given. Table 6 contains values of the ideal-gas part and the residual part r of the
dimensionless Helmholtz free energy together with the corresponding derivatives. Table 7 lists
values for the pressure p, the specific isochoric heat capacity cv, the speed of sound w, and the
specific entropy s calculated at selected values of temperature T and density . Table 8 gives
values for the vapor pressure p, values for the density , enthalpy h, and entropy s for the
saturated liquid, and values for the density , enthalpy h, and entropy s for the saturated
vapor. All these saturation values have been calculated with Eq. (5) by using the phase-
equilibrium condition (see the corresponding comment in Section 4).
12
Table 6. Values for the ideal-gas part , Eq. (6), and for the residual part r,
Eq. (7), of the dimensionless Helmholtz free energy together with the
corresponding derivativesa for T = 500 K and = 46.26 mol dm–3
o = 0.196 026 959 × 101 r = –0.342 285 760× 101
o = 0.384 253 134 r = –0.366 996 322
o r
= –0.147 650 471 = 0.835 582 360
o = 0.939 589 569 × 101 r = –0.589 810 764 × 101
o = –0.209 517 144 × 101 r = –0.244 654 026 × 101
o r
= 0 = –0.113 453 328 × 10
1
a For the abbreviated notation of the derivatives of and r see the footnotes of Tables 4 and 5,
respectively.
Table 7. Thermodynamic property values in the single-phase region for selected values of T and
T/K /(mol dm–3) p/MPa cv/(J mol–1 K–1) w/(m s–1) s/(J mol–1 K–1)
2 a 2 4 1
–1
300 0.551 26 × 10 0.769 726 647 × 10 0.838 021 206 × 10 0.140 382 481 × 10 0.679 005 151 × 10
2 3 2 4 1
0.600 00 × 10 0.238 232 461 × 10 0.741 687 428 × 10 0.177 139 316 × 10 0.545 863 860 × 10
2 3 2 4 1
0.650 00 × 10 0.623 823 936 × 10 0.701 406 088 × 10 0.228 172 318 × 10 0.280 312 863 × 10
1 2 3 3
500 0.500 00 × 10 0.206 052 849 0.294 270 645 × 10 0.514 481 350 × 10 0.140 941 847 × 10
1 2 3 3
0.500 00 0.188 978 723 × 10 0.365 500 382 × 10 0.489 681 127 × 10 0.120 296 840 × 10
2 1 2 4 2
0.462 60 × 10 0.863 679 942 × 10 0.626 150 807 × 10 0.118 054 395 × 10 0.496 096 271 × 10
2 3 2 4 2
0.500 00 × 10 0.107 843 032 × 10 0.616 765 862 × 10 0.148 486 405 × 10 0.469 893 566 × 10
2 3 2 4 2
0.600 00 × 10 0.718 164 137 × 10 0.577 780 154 × 10 0.239 579 704 × 10 0.394 344 288 × 10
2 2 3 3 2
643.8 0.200 00 × 10 0.216 436 864 × 10 0.101 648 572 × 10 0.250 365 829 × 10 0.817 985 464 × 10
1 –1 2 3 3
800 0.100 00 × 10 0.664 864 054 × 10 0.340 033 367 × 10 0.642 794 795 × 10 0.169 130 005 × 10
1 2 3 3
0.250 00 0.164 466 073 × 10 0.344 322 768 × 10 0.639 287 536 × 10 0.142 187 731 × 10
a In the liquid-water region at low pressures, small changes in density along an isotherm cause large changes in pressure. For
this reason, due to an accumulation of small errors, a particular computer code or a particular computer may fail to reproduce
the pressure value with nine significant figures.
13
Table 8. Thermodynamic property values in the two-phase region for selected temperaturesa
T = 280 K T = 450 K T = 625 K
p /MPa 0.822 981 146 × 10–3 0.921 135 794 0.172 108 749 × 102
/ (mol dm–3) 0.552 036 867 × 102 0.492 883 097 × 102 0.306 541 687 × 102
/ (mol dm–3) 0.353 715 935 × 10–3 0.264 016 992 0.694 095 402 × 101
h/(J mol–1) 0.260 035 930 × 103 0.145 280 634 × 105 0.324 689 025 × 105
h/(J mol –1) 0.466 283 044 × 105 0.515 237 547 × 105 0.472 606 343 × 105
s/(J mol–1 K–1) 0.933 712 599 0.407 157 082 × 102 0.731 638 526 × 102
s/(J mol–1 K–1) 0.166 534 672 × 103 0.122 928 355 × 103 0.968 306 235 × 102
a All these test values were calculated from the Helmholtz free energy, Eq. (5), by applying the phase-equilibrium
condition (Maxwell criterion).
9 References
[1] IAPWS, G5-01(2016), Guideline on the Use of Fundamental Physical Constants and Basic
Constants of Water (2001). Available from http://www.iapws.org.
[2] Herrig, S., Thol, M., Span, R., Harvey, A.H., and Lemmon, E.W., A New Reference Equation of
State for Heavy Water, J. Phys. Chem. Ref. Data, to be submitted (2018).
[3] IAPWS, R2-83(1992), Release on the Values of Temperature, Pressure, and Density of Ordinary
and Heavy Water Substances at Their Respective Critical Points (1992). Available from
http://www.iapws.org.
[4] Mohr, P.J., Newell, D.B., and Taylor, B.N., CODATA Recommended Values of the Fundamental
Physical Constants: 2014, J. Phys. Chem. Ref. Data 45, 043102 (2016) [also published in Rev. Mod.
Phys. 88, 035009 (2016)].
[5] Simkó, I., Furtenbacher, T., Hrubý, J., Zobov, N.F., Polyansky, O.L., Tennyson, J., Gamache, R.R.,
Szidarovszky, T., Dénes, N., and Császár, A.G., Recommended Ideal-gas Thermochemical Functions
for Heavy Water and its Substituent Isotopologues, J. Phys. Chem. Ref. Data 46, 023104 (2017).
[6] Markó, L., Jákli, G., and Jancsó, G., Vapour pressure of heavy water at its triple point, J. Chem.
Thermodyn. 21, 437-441 (1989).
14
Fig. 1. Expanded relative uncertainties in density, / , estimated for Eq. (5). In the enlarged critical region
(triangle), the uncertainty is given as percentage uncertainty in pressure, p /p. This region is bordered by
the two isochores 8 mol dm3 and 29 mol dm3 and by the 30 MPa isobar. The positions of the lines
separating the uncertainty regions are approximate. At low pressures for the vapor, the uncertainties become
much smaller than indicated because the vapor is nearly an ideal gas.
15
Fig. 2. Expanded relative uncertainties in speed of sound, w/w, estimated for Eq. (5). For the definition of the
triangle around the critical point, see the Fig. 1 caption. The positions of the lines separating the uncertainty
regions are approximate. At low pressures for the vapor, the uncertainty becomes small because the vapor
approaches the ideal-gas limit.
16
Fig. 3. Expanded relative uncertainties in specific isobaric heat capacity, cp /cp , estimated for Eq. (5). For the
definition of the triangle around the critical point, see the Fig. 1 caption. The positions of the lines
separating the uncertainty regions are approximate. The uncertainty in the vapor phase at low pressures
approaches the uncertainty of the ideal-gas heat capacity, which is less than 0.02 %.
17
Fig. 4. Uncertainties in vapor pressure, p /p in saturated liquid density, / , and in saturated vapor density,
/ , estimated for Eq. (5). The uncertainties for the saturated densities increase to 1.5 % at the critical
temperature.