Co2 Equation of State PDF
Co2 Equation of State PDF
This work reviews the available data on thermodynamic properties of carbon dioxide
and presents a new equation of state in the form of a fundamental equation explicit in the
Helmholtz free energy. The function for the residual part of the Helmholtz free energy
was fitted to selected data of the following properties: (a) thermal properties of the
single-phase region (ppT) and (b) of the liquid-vapor saturation curve (Ps> p', p")
including the Maxwell criterion, (c) speed of sound w and (d) specific isobaric heat
capacity c p of the single phase region and of the saturation curve, (e) specific isochoric
heat capacity c v' (f) specific enthalpy h, (g) specific internal energy u, and (h) Joule-
Thomson coefficient /-t. By applying modern strategies for the optimization of the math-
ematical form of the equation of state and for the simultaneous nonlinear fit to the data of
all these properties, the resulting formulation is able to represent even the most accurate
data to within their experimental uncertainty. In the technically most important region up
to pressures of 30 MPa and up to temperatures of 523 K, the estimated uncertainty of the
equation ranges from ± 0.03% to ± 0.05% in the density, ± 0.03% to ± 1% in the speed
of sound, and ± 0.15% to ± 1.5% in the isobaric heat capacity. Special interest has been
focused on the description of the critical region and the extrapolation behavior of the
formulation. Without a complex coupling to a scaled equation of state, the new formu-
lation yields a reasonable description even of the caloric properties in the immediate
vicinity of the critical point. At least for the basic properties such as pressure, fugacity,
and enthalpy, the equation can be extrapolated up to the limits of the chemical stability of
carbon dioxide. Independent equations for the vapor pressure and for the pressure on the
sublimation and melting curve, for the saturated liquid and vapor densities, and for the
isobaric ideal gas heat capacity are also included. Property tables calculated from the
equation of state are given in the appendix. © 1996 American Institute of Physics and
American Chemical Society.
Key words: carbon dioxide; correlation; critical region; data evaluation; equation of state; extrapolation;
fundamental equation; melting line; property tables; sublimation line; thermal and caloric properties; and
vapor-liquid coexistence curve.
Contents
0047-2689/96/25(6)/1509/88/$20.00 1509 J. Phys. Chern. Ref. Data, Vol. 25, No.6, 1996
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
1510 R. SPAN AND W. WAGNER
Capacity in the Ideal-Gas State. . . . . .. 1529 5. Selected data for the critical point of carbon
4.3 Specific Isochoric Heat Capacity. . . . . . . . . .. 1530 dioxide ................................. " 1520
4.4 Speed of Sound. . . . . . . . . . . . . . . . . . . . . . . .. 1532 6. Summary of the data sets for the sublimation
4.5 Enthalpy............................... 1532 pressure of carbon dioxide. . . . . . . . . . . . . . . . . .. 1521
4.6 Internal Energy. . . . . . . . . . . . . . . . . . . . . . . .. 1532 7. Summary of the data sets for the vapor pressure
4.7 Joule-Thomson Coefficient. ............. , 1532 of carbon dioxide. . . . . . . . . . . . . . . . . . . . . . . . .. 1522
4.8 Virial Coefficients. . . . . . . . . . . . . . . . . . . . . .. 1533 8. Summary of the data sets for the saturated liquid
4.9 Liquid-Vapor Equilibrium. . . . . . . . . . . . . . .. 1533 density of carbon dioxide. . . . . . . . . . . . . . . . . . .. 1523
4.10 Adjustment of Data. . . . . . . . . . . . . . . . . . . .. 1534
9. Summary of the data sets· for the saturated vapor
4.10.1 Adjustment of Data Sets Describing
density of carbon dioxide .................. " 1524
the Critical Region. . . . . . . . . . . . . .. 1534
10. Summary of the data sets for caloric properties
4.10.2 Adjustment of ppT Data. . . . . . . . . .. 1534
on the liquid-vapor phase boundary of carbon
4.10.3 Correction of Isobaric Heat
dioxide. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. 1526
Capacities. . . . . . . . . . . . . . . . . . . . .. 1535
5. Description of Thermodynamic Propert.ies in the 11 _ Sllmm~ry of the dat~ sets available for the ppT
Critical Region. . . . . . . . . . . . . . . . . . . . . . . . . . . .. 1535 relation of carbon dioxide ................... , 1527
5.1 Limitations of Analytical Equations of State. 1536 12. Summary of selected ppT data for carbon
5.2 Use of Nonanalytic Terms as an Integral dioxide; detailed information is given on
Component in an Empirical Wide-Range the uncertainty values estimated by the authors
Equation of State. . . . . . . . . . . . . . . . . . . . . . .. 1537 and those estimated by ourselves and used
6. New Equation of State. . . . . . . . . . . . . . . . . . . . . .. 1541 in the weighting procedure .................. , 1529
6.1 Ideal-Gas Part of the Helmholtz Energy. . . .. 1541 13. Summary of the data sets available for the
6.2 Residual Part of the Helmholtz Energy. . .. .. 1543 specific isobaric heat capacity of carbon dioxide.. 1530
7. Comparisons of the New Equation of State with 14. Summary of selected data for the specific
Experimental Data and Other Equations of State_ ] 54() isobaric heat capacity of carbon dioxide; detailed
7.1 Liquid-Vapor Boundary .................. 1546 information is given on the uncertainty values
7.1.1 Thermal Properties on the Coexistence estimated by the authors and those estimated by
Curve. . . . . . . . . . . . . . . . . . . . . . . . . . .. 1546 ourselves and used in the weighting procedure. .. 1530
7.1.2 Caloric PropertIes on the CoeXIstence 15. Data sets for the isobaric heat capacity in the
Curve ............................ 1547 ideal-gas state of carbon dioxide calculated by
7.2 Single-Phase Region. . . . . . . . . . . . . . . . . . . .. 1547 theoretical approaches. . . . . . . . . . . . . . . . . . . . . .. 1530
7.2.1 Thermal Properties in the Single-Phase 16. Summary of the data sets available for the
Region. . . . . . . . . . . . . . . . . . . . . . . . . .. 1547
specific isochoric heat capacity of carbon
7.2.2 Caloric Properties in the Single-Phase
dioxide. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. 1531
Region. . . . . . . . . . . . . . . . . . . . . . . . . .. 1553
17. Summary of selected data for the specific
7.3 Extrapolation Behavior of the New
isochoric heat capaci~y of carbon dioxide; detailed
Fundamental Equation. . . . . . . . . . . . . . . . . .. 1556
information is given on the uncertainty values
7.3.1 Extrapolation Beyond the Range of
Primary Data. . . . . . . . . . . . . . . . . . . . .. 1557 estimated by the authors and those estimated by
7.3.2 Representation of "Ideal Curves"..... 1557 ourselves and used in the weighting procedure. .. 1531
8. Uncertainty of the New Fundamental Equation... 1559 18. Summary of the data sets available for the speed
9. Conclusions................................ 1559 of sound of carbon dioxide. . . . . . . . . . . . . . . . . .. 1531
10. Appendix: Thermodynamic Properties of Carbon 19. Summary of selected data for the speed of sound
Dioxide. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. 1559 of carbon dioxide; detailed information is given
11. Acknowledgments. . . . . . . . . . . . . . . . . . . . . . . . .. 1593 on the uncertainty values estimated by the
12. References.... . . . . . . . . . . . . . . . . . . . . . . . . . . . .. 1593 authors and those estimated by ourselves and
used in the weighting procedure. . . . . . . . . . . . . .. 1532
20. Summary of the data sets available for enthalpy
List of Tables
1. Available wide range equations of state for differences of carbon dioxide. . . . . . . . . . . . . . . .. 1532
carbon dioxide. . . . . . . . . . . . . . . . . . . . . . . . . . . .. 1515 21. Summary of the data sets available for differences
2. Selected scaled equations of state for carbon of the internal energy of carbon dioxide. . . . . . . .. 1532
dioxide. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. 1515 22. Summary of the data sets available for the Joule-
3. Relations of thermodynamic properties to the Thomson coefficient of carbon dioxide. . . . . . . . .. 1533
dimensionless Helmholtz function cP consisting 23. Summary of data available for the second and
of cPo and d/; see Eq. (2.1). . . . . . . . . . . . . . . . . .. 1517 third virial coefficient of carbon dioxide. For
4. Selected data for the triple point of carbon reasons explained in the text, no data were
dioxide. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. 1520 assigned to Group 1. . . . . . . . . . . . . . . . . . . . . . . .. 1533
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
A NEW EQUATION OF STATE FOR CARBON DIOXIDE 1511
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
1512 R. SPAN AND W. WAGNER
15. For temperatures between Tc and Tc + 1 K Eq. 24. Relative density deviations of selected ppT data
(6.0 results in an oscillating plot for the from values calculated from Eq. (6.0. Values
isobaric heat capacity around the critical density.. 1539 calculated from the wide-range equations of Ely
16. Relative deviations of c~ data from values et al. 15 and Angus et ai. 3 are plotted for
calculated from Eq. (6.2). The upper diagram comparison.. . . . . . . . . . .. . . . . . . . . . . . . . . . . . . .. 1551
shows data calculated from statistical 25. Relative density deviations of selected ppT data
at high temperatures from values calculated from
thermodynamics (see Table 15) and the lower
Eq. (6.1). Values calculated from the
diagram shows data extrapolated from
wide-range equations of Ely et ai. 15 and Angus
experimental results (see Table l3). Values of c~ 3
et al. are plotted for comparison.. . . . . . . . . . . . .. 1552
calculated from the corresponding equations
26. Relative density deviations of selected ppT data
of Angus et ai. 3 and Ely et ai. 15 are plotted for at high pressures from values calculated from
comparison.. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. 1540 Eq. (6.1). Values calculated from the wide-range
17. Distribution of the experimental data used for equations of Ely et ai. 15 and Angus et ai. 3 are
the establishment of the residual part of the plotted for comparison~ in this pressure range,
new fundamental equation, Eq. (6.5), in apT these two equations of state are at least partly
diagram.. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. 1542 extrapolated (see Table 0.. . . . . . . . . . . . . . . . . . .. 1552
18. Relative deviations 1OOLly = 100(Yexp -Yca1c)/Yexp 27. Relative deviations of selected isobaric heat
(y = p s ,p I , p") of the experimental saturation capacity data from values calculated from
data of Duschek et ai. 58 from values calculated Eq. (6.1). Values calculated from the wide-range
from Eq. (6.1). Values calculated from the equations of Ely et al. 15 and Angus et al. 3 are
auxiliary equations presented in Sec. 3, the plotted for comparison ..................... " 1553
equation of state of Ely et al. 15 and the auxiliary 28. Representation of the isobaric heat capacity on
equations of Angus et al. 3 are plotted for isobars in the gas region and in states on the
sublimation curve (saturated vapor). Values
comparison ................................ , 1546
calculated from the wide-range equations of Ely
19. Relative deviations 100Lly= 100(Yexp -Yca1c)IYexp
et al. 15 and Angus et al. 3 are plotted for
(Y=Ps' p', p") of the near critical experimental
comparison.. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. 1553
saturation data of Duschek et ai.58 from
29. Relative deviations of isobaric heat capacity data
values calculated from Eq. (6.1). Values
at high temperatures from values calculated from
calculated from the auxiliary equations presented Eq. (6.1). Values calculated from the
in Sec. 1 ann from the cro~~()ver equation of wide-range equations of Ely et al. 15 and Angus
Chen et ai.27 are plotted for comparison.. . . . . . .. 1546 et ai. 3 are plotted for comparison.. . . . . . . . . . . . .. 1554
20. Relative deviations 100Lly=100(Yexp-Ycalc)IYexp 30. Relative deviations of selected isochoric heat
(y = w" , w', c;) of experimental caloric data at capacity data from values calculated from Eq.
saturation from values calculated from Eq. (6.1). (6.1). Values calculated from the wide-range
Data calculated from the wide-range equations equations-bf Ely et cit. IS and Angus et ai. 3 are
of Ely et al. 15 and Angus et ai. 3 are plotted for plotted for comparison.. . . . . . . . . . . . . . . . . . . . .. 1554
comparison ................................ , 1547 31. Representatiun uf the isuchurk heat L:apadty UIl
21. Relative density deviations of very accurate ppT high-density isochores. For each of the isochores
data at subcritical temperatures from values the plotted pressure range starts at the
calculated from Eq. (6.1). Values calculated from corresponding vapor pressure. Values calculated
the wide-range equations of Ely et al. 15 and from the wide-range equations of Ely et al. 15
Angus et al. 3 are plotted for comparison ........ , 1548 and Angus et al. 3 are plotted for comparison.. . .. 1554
32. Relative deviations of speed of sound data at
22. Relative pressure deviations of very accurate
supercritical temperatures from values
ppT data in the extended critical region from
calculated from Eq. (6.1). Values calculated
values calculated from Eq. (6.1). Values
from the wide-range equations of Ely et al. 15 and
calculated from the wide-range equations of Ely of Pitzer and Schreiber,16 ::InrI, in the range of
et al. 15 and Pitzer and Schreiber 16 and -from
its validity, from the crossover equation of Chen
the crossover equation of Chen et al. 27 are plotted et aZY are plotted for comparison.. . . . . . . . . . . .. 1555
for comparison ............................. , 1549 33. Relative deviations of speed of sound data at
23. Relative density deviatiuns uf velY aL:L:uUtte fJpT high pressures from values calculated from
data at supercritical temperatures from values Eq. (6.1). Values calculated from the wide-range
calculated from Eq. (6.0. Values calculated from equations of Ely et al. 15 and Angus et al. 3 are
the wide-range equations of Ely et 0/ 15 ::Ind plotted for comparison; in this pressure range,
Angus et ai. 3 are plotted for comparison.. . . . . . .. 1550 t~ese two equations of state are at least partly
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
A NEW EQUATION OF STATE FOR CARBON DIOXIDE 1513
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
1514 R. SPAN AND W. WAGNER
Rm Molar gas constant: R m=(8.314 510±0.000 210a) stance for theoretical approaches dealing with the critical re-
J/(mol K); see Ref. 2 gion of pure fluids. Almost every physical model for the
R Specific gas constant: R=R mIM=(0.188 9241 description of the critical region has been tested for carbon
± 0.000011 6) kJ/(kg K) dioxide.
Tc Critical temperature: Tc =(304.1282±0.015)K; see In 1965, an international research project on the thermo-
Sec. 3.2 dynamic properties of carbon dioxide was established at the
Pc Critical pressure: Pc=(7.377 3±0.003 O)MPa; see Sec. instigation of the International Union of Pure and Applied
3.2 Chemistry (IUPAC). In 1976, Angus et al. 3 published a
Pc Critical density: Pc=(467.6±0.6) kg/m3 ; see Sec. 3.2 monograph which reviewed the experimental data available
Tt Triple-point temperature: T t =(216.592± 0.002) K; see up to 1973 and presented extensive tables of the thermal and
Sec. 3.1 caloric properties derived from the selected equation of state.
Pt Triple-point pressure: Pt=(0.517 95± 0.000 10) MPa; Nevertheless, knowledge of the thermodynamic properties of
see Sec. 3.1 carbon dioxide remained unsatisfactory. Thus, since 1973
To Reference temperature: T o=298.15 K numerous experiments including state-of-the-art experiments
Po Reference pressure: Po=0.101 325 MPa with significantly improved accuracy have been performed in
hg Reference enthalpy in the ideal gas state at To: order to improve the quality of the entire data set; a third of
hg=O kJ/kg the data sets available today belongs to this group.
sg Reference entropy in the ideal gas state at To, Po: In addition to the increased amount of experimental data,
s3=0 kJ/(kg K) correlation techniques have significantly improved during
the last decade. Sophisticated "multi-property" fitting
procedures4 ,5 and a new strategy for optimizing the structure
1. Introduction of empirical correlation equations6 have resulted in a new
basis for the development of an empirical equation of state.
1.1 Background
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
A NEW EQUATION OF STATE FOR CARBON DIOXIDE 1515
"Benedic(-Webb-Rubin.
bNo information published.
cTerms developed by Ewers and Wagner l3 for the description of the critical region were also used.
dBesides some experimental data, mainly data from the IUPAC tables 3 were used as input data.
CRp~inp<.: mmp pxpprimpntal nata, mainly clata calculated from the equation of Ely et al. 17 were used as input data.
et a1. 15 ,17 These equations take into account almost all of the • Extrapolation to temperatures and pressures outside the
published experimental results. However, some of the most range of validity yields unreasonable results.
important state-of-the-art experiments on carbon dioxide • The temperature values used do not correspond to the cur-
were not yet available at that time. A detailed comparison21 rent International Temperature Scale of 1990 (ITS-90).
shows that the: N]natinn givf':n in Re:f~ 14 :lno 1'i lI~ine thf> Tht: ~4ualiull uf Pilz~l ami Sl~m~I lX~ L:aIl mainly be L:UIl-
form developed by Schmidt and Wagner4 is superior to the sidered to be an example for an empirical equation of state
equation given in Ref. 17 which uses the fonn developed by with reasonable extrapolation behavior; the authors do not
Jacobsen and Stewart. 18 claim that this equation improves the description of thermo-
Tests of these correlations have YIelded that all of the ex- dynamic properties in the region where accurate data exist.
isting equations of state, independent of their different qual- The problems related to the description of properties
ity, show the following limitations: within the critical region can essentially be solved by using
• State-of-the-art data for the ppT relation are not repre- scaled equations of state; Table 2 shows selected correlations
sented to within their experimental uncertainty. with such scaling approaches for the description of the criti-
• State-or the-art data describing the liquid-vapor phase cal region of carbon dioxide. However, the relatively small
equilibrium are not represented to within their experimen- range of validity and the mathematically complex structure
tal uncertainty. of these equations limit their use.
• Within the critical region, the calculation of caloric prop- An exception with regard to the range of validity is the
erties yields unreasonable results. empirical equation of Erickson et at. 26 This equation uses an
• Unreasonable behavior can be observed in regions wIth a Improved form of the transtormatIOn procedure, originally
poor data situation. developed by Fox,29 to achieve a reasonable but not an as-
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
1516 R. SPAN AND W. WAGNER
ymptotically correct description of the caloric properties in 2. Basic Elements of the Development of
the critical region. Like other parametric approaches, this Equations of State in Form of a
equation still requires complicated iteration procedures, but
Fundamental Equation
it is valid in a range which is sufficiently wide for engineer-
ing application~. Tests cdllieu uut by Erkbun e( al. '"
showed that outside the critical region the transformed equa- The new equation of state, developed for the representa-
tion based on the form developed by Jacobsen and Stewart 18 tion of the thermodynamic properties of carbon dioxide, is an
is inferior to Ely's equation. 14,lS empirical representation of the fundamental equation explicit
The purpose of this article is to present a new equation of in the Helmholtz energy. Since the application of optimiza-
state for carbon dioxide explicit in the Helmholtz energy, tion procedures. and the technique of mUlti-property fitting
which is designed to overcome the disadvantages of the ex- were extensively discussed by Setzmann and Wagner6 ,30 and
isting correlations. The consistent use of sophisticated Saul and Wagner,S this section is restricted to some basic
fitting 30 and optimization6 procedures allows the representa- facts giving a rough understanding of the procedures used for
tion of even the most accurate data within their experimental the development of the new equation. However, those facts
uncertainty. Extensive investigations of the ci::lt~ set prevent which have not previously been published are outlined in
the equation from unreasonable behavior within the fitting more detail.
region, and new procedures are introduced to guarantee a 2.1 Helmholtz Function
reasonable extrapolation behavior of the new equation. The
representatIOn ot calonc properties in the critical region has The fundamental equation given in this article IS expressed
been improved by an empirical approach which yields rea- in form of the Helmholtz energy A with the two independent
sonable results even in the immediate vicinity of the critical variables density p and temperature T. The dimensionless
point. Helmholtz energy ¢J=A/CRT) is commonly split into a part
All the correlation equations given in this article corre- depending on the ideal-gas behavior ¢J0 and a part which
spond to the ITS-90 temperature scale. takes into account the residual fluid behavior ¢Jr, namely
r/J( 8, 'T) = r/J 0 ( 8, 'T) + tll( 8, 'T), (2.1)
1.3 Organization of the Article where o=p/Pc is the reduced density and 7=TcIT is the
inverse reduced temperature. Both the density p and the tem-
In Sec. 2 we give a brief review of the application of perature T are reduced with their critical vallles, Pc ~mci To;,
empirical equations of state explicit in the Helmholtz energy respectively.
as a function of density and temperature; this combination of Since the Helmholtz energy as a function of density and
variables is one form of a fundamental equation. The tech- temperature is one form of a fundamental equation, all the
niques used for the development of the new equation (multi- thermodynamic properties of a pure substance can be ob-
property fitting, optimization procedure) are described tained by combining derivatives of Eq. (2.1). Table 3 gives
briefly. Section 3 deals with the liquid-vapor, solid-vapor, the relations between Eq. (2.1) and the properties considered
and solid-liquid phase equilibria of carbon dioxide. The in this article. The vapor pressure and the densities of satu-
available data are discussed. In addition to the equation of rated liquid and saturated. vapor can be determined from an
state, short supplementary equations have been developed equation explidt in ¢J by simultaneous solution of the equa-
for the temperature dependence of the melting pressure tions
Pm' the sublimation pressure P sub' the vapor pressure p s ,
the saturated liquid density p', and the saturated vapor den-
sity p". The experimental information available for the
R~~' = I + 0' ¢J~( 8',7), (2.2a)
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
A NEW EQUATION OF STATE FOR CARBON DIOXIDE 1517
TABLE 3. Relations of thermodynamic properties to the dimensionless Helmholtz function 4> consisting of cpo and cpr; see Eq. (2.1)
(l+&/J~-8rdJ~~)-:~ ?;]
. (dPs/dT)/(apIJvhlu=v'
.[
Speed of sound: 2
w ( 8, r) _ r 2 r (1 + 84>~- 8r4>~".)2
w(T,p) = ~(Jplap)s ~ -1 +284>0+ 84>00- ~(4)~T+ 4>~".)
J oule-l homson coefficient: - (OcfJ8T 0 2 cfJ80 T orcfJ8T)
/-t(T,p)=(aTIJph
Fugacity:
The enthalpy hOof the ideal gas is a function of temperature stance. Thus, the equation for the residual part has to be
only, and the entropy SO of the ideal gas depends on tempera- determined in an empirical way by optimizing its functional
ture and density. Both properties can be derived from an form and fitting its coefficients to experimental results. b
equation for the ideal-gas heat capacity c~. When c~ is in- Since the Helmholtz energy is not accessible to direct mea-
serted into the expression for hO(T) and SO(p, T) in Eq. (2.3) surements, a procedure has been developed which allows the
one obtains establishment of an empirical equation for the residual part
of the Helmholtz energy by using data of different proper-
ties.
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
1518 R. SPAN AND W. WAGNER
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
A NEW EQUATION OF STATE FOR CARBON DIOXIDE 1519
set by suitable data selection if wt> O. of the decimal point were used. Data corresponding to the
• After correcting systematic deviations, a data set yields IPTS-48 temperature scale were converted to IPTS-68 ac-
results much better than expected from its original experi- cording to the 'procedure given by Bedford and Kirby.39 For
mental uncertainty if wt> 1). temperatures between 90 K and 900 K, the temperature
• In a special region, enlarged, but difficult assignable un- scales ITS-27 and IPTS-48 do not deviate from each other
certainties of the selected data are expected if wt < 1). and data older than 1927 were not used.
The algorithm used for the conversion from IPTS-68 to
To achieve as much transparency with regard to the data
the ITS-90 temperature scale 38 causes an additional uncer-
set as possible, the tables presenting the selected data (see
tainty of ± l.5 mK for temperatures below 273.15 K and
Sec. 4) contain additional information on the uncertainties
± 1 mK for temperatures above 273.15 K. This additional
used in the weighting procedure as well as on the mean value
uncertainty is not considered in the uncertainties given in the
of the weighting factor. Weighting factors deviating signifi-
tables of this section, since these uncertainties were mainly
cantly from fwt= I are discussed.
used for consistency tests between data of different authors.
In this case, the uncertainty in the absolute temperature,
3. Phase Equilibria of Carbon Dioxide which is influenced by the uncertainty of the conversion, is
less important. The comparison between two very similar
temperature values is not influenced by the uncertainty of the
An accurate description of phase equilibria by auxiliary
t::Ljualium; i~ an important precondition for the development
conversion if both temperatures are converted with the same
procedure.
of a wide-range equation of state and is also helpful for users
who are only interested in phase equilibria. Therefore, all 3.1 Triple Point
available experimental information on the triple point, the
critical point, the melting pressure, the sublimation pressure, Uunng the last YU years, the trIple-pomt temperature ot
the vapor pressure, the densities of saturated liquid and va- carbon dioxide has been determined by numerous authors
por, and on caloric properties on the liquid-vapor phase but the data situation for the triple-point pressure is rather
boundary have been reviewed. With the exception of the poor. Fortunately, the few measurements available are very
caloric properties, simple correlation equations have been consistent with each other. Table 4 shows selected data of
developed for the temperature dependency of these quanti- the triple point of carbon dioxide. After a comprehensive
ties. review of the existing data, we have chosen the following:
To condense the description of the data situation, the char- T t = (216.592± 0.003)K, (3.1)
acteristic information on the single data sets is summarized
in tables for the corresponding property. The data sets have Pt=(O.51795± O.OOOlO)MPa. (3.2)
been divided into three groups. The assignment considers the
critically assessed uncertainty of the data, size of the data set, Data of the density of the saturated liquid and the satu-
and covered temperature range. In addition, attention is paid rated vapor at the triple point are not available, but the evalu-
to the data situation for the respective property. Group I ation of the corresponding correlation equations given in
contains the data sets used for the development of the corre- Sees. 3.6 and 3.7 yields
sponding correlation equation. Group 2 contains data sets p;'=(1178.53± O.18)kg/m 3 , (3.3)
suitable for comparisons. Compared with group I data, these
data decline at least under one of the three aspects mentioned p;' = (13.76l4::!:: O.0034)kg/m 3 • (3.4)
above. Group 3 contains very small data sets and data sets
with rather high uncertainty. Consideration of these data is 3.2 Critical Point
not rt':lson:lhle on the level of accuracy aspired to here. Nev-
ertheless this means no devaluation of these data sets - the Altogether, data of the critical point of carbon dioxide are
whole ranking is intluenced more by the quality in relation to given in 75 papers. Table 5 shows selected values of the
the best available reference data than by an absolute level of critical temperature, the critical pressure, and the critical den-
uncertainty; for other purposes even group J data sets may SIty. The values found for the CrItIcal denSIty agree well
be very useful. within the expected uncertainties, but the values for the criti-
Since the new correlation equations and all temperature cal temperature show significant differences, far beyond the
values in this article correspond to the ITS-90 temperature uncertainties given by some of the authors. Essentially, the
scale. 36 the temperature values of the available data, based differences in the critical pressures can be explained by the
on older temperature scales, were converted to ITS-90. The variation of the vapor pressure with the assumed critical tem-
conversion from the IPTS-68 temperature scale to ITS-90 perature.
temperatures was carried out based on the internationally In the course of this work, no new attempt has been made
agreed procedure of Preston-Thomas et al.,37 explained by to determine the critical parameters of carbon dioxide, but
Rusby.38 The number of digits of the converted values was the evaluation of Duschek et al. 58 was tested under different
increased by one digit in order to guarantee numerically con- aspects. No reason for altering the data became obvious.
sistent reconversion. but not more than four digits to the right Thus, we have used
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
1520 R. SPAN AND W. WAGNER
aAll temperatures were converted to ITS-90. Up to an uncertainty of 0.1 mK an additional digit is added to
consistent reconversion to IPTS-68 temperatures.
uncertainty of the conversion from IPTS-68 to ITS-90 is not considered here.
Tc= (304.1282± O.015)K, temperature assumed Michels et al. 6o and Clusius et al. 61
deviate from the rel.:t:.Ill valul;; [eI. Dg. (J.l»), llll their tem-
Pc= (7.3773± O.0030)MPa, (3.6) perature values were corrected according to
Pc=(467.6± (3.7)
K, (3.8)
(The additional digit in Tc is given to guarantee numerical
consistency in any reconversion to the IPTS-68 temperature K (3.9)
scale.)
after their conversion to ITS-90 temperatures.
3.3 Pressure The corrected data are consistent within the expected un-
certainties and were used to establish a simple correlation
There are only two available sets of measurements de- equation for the melting pressure:
sen bing the metring pressure uf cwbun dioxide. In 1942
Michels et al. 60 measured 25 melting pressures at tempera- Pm
tures between 217 K and 266 K, to pressures
(3.10)
Pt
between 0.9 MPa and 284 MPa. In Clusius et al. 61
published 21 data, covering the range from 2t7 K to 222 K, with 1~=:Z16.592 K, pt=O.51795 MPa, 01 = 1955.5390, and
corresponding to 0.5 MPa to 24 MPa. a2=2055.4593. The equation is constrained to the tripIe-
Unfortunately, these data sets are inconsistent with each point pressure by its functional form. The representation of
other and with recent data of the triple point pressure. In both both the corrected and uncorrected data is shown in Fig. L.
experiments, the thennometers were calibrated at the triple- The uncertainty of Eq. (3.10) is estimated to be
point of carbon dioxide. Since the values of the triple-point tiPmIPm~± 1.5% for Tt~T~225 K and :!:O.S% for
'All temperatures converted to ITS-90. Up to four digits beyond the decimal point an additional digit is
added to guarantee consistent reconversion to IPTS-68 terntper'atures.
trrhe uncertainty of the conversion from [PIS-68 to not considered here.
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
A NEW EQUATION OF STATE FOR CARBON DIOXIDE 1521
-I::l...
Tt
~ -lOL---~mH~~~--~--~--~--~--~
In (
PSUb)Tt ( (
Pt = T' all - ~
T)
+ a2 ( 1 - ~
T) L9
-:~
~ +a3 1-~ (3.12)
I::l...
..
:! :
with Tt =216.592 K, pt=0.51795 MPa, al = -14.740846,
!
~
0
! a2=2.4327015, and a3= -5.3061778. The equation is con-
I
strained to the triple-point pressure by its functional form.
Figure 2 compares group 1 and group 2 data with values
220 230 240 250 260 270 calculated from Eq. (3.12). The uncertainty of the new equa-
Temperature T / K tiul1 is csLilllaLcll Lo uc ~Psub:::S -'- 2.50 Pa fur 18.5 K::::::: T~ Tt ,
61 61 ± 100 Pa for 170 K::::::; T::::::; 185 K, and dPsub::::::; ± 50 Pa
dPsub::::::;
o Clusius et al. • Clusius et al. (corr. data)
6O for T::::::; 170 K. Bedford et al. 47 concluded that the most recent
o Michels et al. 60 • Michels et al. (corr. data)
value for the sublimation temperature at normal pressure
FIG. 1. Relative deviations of experimental melting pressure data from val- (P0=0.101325 MPa) is T= 194.6857 K±0.0030 K~ iteration
ues calculated from the melting pressure equation, Eq. (3.10). In this figure, with Eq. (3.12) yields T=194.6855 K.
both the corrected and the uncorrected data are plotted (see Sec. 3.3).
Equation (3.12) was fitted to data at temperatures above
1.54 K. The low temperature data of Bryson er ai. 75 are not
represented to within the uncertainty given by the authors,
225 K < T~ 270 K. The simple form of the equation should
but the deviations do not exceed 0.01 Pa for T >85 K and
ensure reasonable extrapolation to higher temperatures.
0.0001 Pa for T::::::; 85 K, respectively.
TABLE 6. Summary of the data sets for the sublimation pressure of carhon dioxid~
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
1522 R. SPAN AND W. WAGNER
500
~
0...
........
.0
;:l 0
~
<'J
-500
150 160 170 180 190 200 210 220
Temperature T / K
o Ambrose 72 • Bedford et al. 47
o Giauque and Egan 70 o Fernandez-Fassnacht and del Rio 76
FIG. 2. Absolute deviations ~Psub= (Psub.exp - Psub.calc) of selected experimental sublimation pressure data from values calculated from the sublimation pressure
equation (3.12). In this figure. both the corrected and the uncorrected data of Rilkarli pf 0174 ~re plotteti (see Se~ ~4)
TABLE 7. Summary of the data sets for the vapor pressure of carbon dioxide
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
A NEW EQUATION OF STATE FOR CARBON DIOXIDE 1523
0.05~--~~~--~~~----------------------------~~--------~~~~~~
0.
x
0
0
-..~
u
(;l
u. -0.05
,~
I
, p"
x
l!- 0.05
~ I> <I
'-"'
0
a
-"
0
• Duschek et al. S8
3
- - - Angus et al. (aux. eq.)
FIG. 3. Relative deviations of selected experimental vapor pressure data from values calculated from the vapor pressure equation, Eq. (3.13). Vapor pressures
calculated from the corresponding equation of Angus et al. 3 are plotted for comparison.
TABLE 8. Summary of the data sets for the saturated liquid density of. carbon dioxide
"The data sets of Esperlll and Holste et al. 103 are different evaluations of the same measurements.
h At temperatures above 300 K. the uncertainty rises to :!: 0.2«<.
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
1524 R. SPAN AND W. WAGNER
0.1~------------------------------------------~~----------------~
----~~£>
Q..
<1
o
<>
,/
,,--0- <>
<><>
<>
<><>
<> •
o
~
-0.1~~~~~~~~~~--~~~~--~~~~--~~~~
220 240 260 280 30
Temperature T / K
lll 3
<l Holste et at. 103 t> Esper - - - Angus et at. (aux. eq.)
llO 58
<> IIayncs • Duschek et al.
FIG. 4. Relative deviations 100 dp'=100 (p~xp-P~alc)lp~xp of selected experimental saturated liquid density data from values calculated from Eq. (3.14).
Saturated liquid densities calculated from the corresponding equation of Angus et at. 3 are plotted for comparison.
zation procedure yields the same mathematical form used by the experimemal uncenaimy of (hese data, (he uncenaimy of
Duschek et al. 58 to describe this data set in the IPTS-68 tem- Eq. (3.13) is estimated to be APs~::t 0.012% for the whole
perature scale; only the coefficients ai have changed. The temperature range.
vapor pressure equation can be written as The dashed line in Fig. 3 corresponds to values calculated
from the vapor pressure equation of Angus et al. 3 after con-
version to ITS-90. The results of this correlation are remark-
(3.13) ably good but, of course, the most accurate data available
today are not represented to within their experimental uncer-
with Tc=304.1282 K, Pc=7.3773 MPa, al=-7.060 2087, tainty.
a2= 1.9391218, a3= -1.6463597, a4 = -3.2995634, tl = 1.0,
t 2=1.5, t3=2.0, and t4=4.0. 3.6 Saturated Liquid Density
Figure 3 shows both the good agreement between the data
sets marked as group 1 and group 2 data and the representa- Table 8 shows information on the 17 data sets of the satu-
tion of the data measured by Duschek et al. 58 Considering rated liquid density of carbon dioxide. Again, only the data
TABLE 9. Summary of the data sets for the saturated vapor density of carbon dioxide
aT he data sets of Esperlll and Holste et al. 103 are different evaluations of the same measurements.
b At temperatures above 295 K, the uncertainty rises to ± 0.25%.
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
A NEW EQUATION OF STATE FOR CARBON DIOXIDE 1525
0.5~---4RH~~r-----~~------------~~-------------t
e.. 0
:: ~
~
.........
,-..... Tt
u
:] -0.5
,~
..e..
220 240 260 280 300
:: ~ 0.05
~
'-'
0
0
,...-t
0
-0.05
Tt •
220 240 260 280 300
Temperature TIK
<l Holste et al. 103
0 BaadeS4 D.. Straub109
I> Esper
l11
•• Duschek et al. 58
Michels et at. 80
94
X Kholodov et al. <> Abdulagatov et al. S9
3
_ _ _ Angus et at. (aux. eq.)
FIG.5. Relative deviations of selected experimental saturated vapor density data from values calculated from Eq. (3.15). Saturated vapor densities calculated
from the corresponding equation of Angus et ai. 3 are plotted for comparison.
set of Duschek et al. 58 is used to develop a simple saturated Figure 4 shows the poor situation regarding group I and
liquid density equation. Group 2 is restricted to data sets group 2 data and the highly accurate representation of the
which deliver uncertainties comparable with common uncer- group 1 data. The dashed line corresponds to values calcu-
tamtIes in the single phase region, that means lated from an auxiliary equation given by Angus et al. 3 for
Llp'::S±0.2% outside the critical region. Unfortunately, this the saturated liquid" density. .
criterion is only met by the data sets of Haynes, 110 Esper, III
and Holste et al.; 103 all the other data proved to be less ac- 3.7 Saturated Vapor Density
curate.
For the whole temperature range, the simple correlation Data on the saturated vapor density of carbon dioxide are
given in 15 sources, which are listed in Table 9. If the same
In ( ~') =? 4
ai
( T;
1-- t·
1 (3.14) criterion as mentioned in Sec. 3.6 is used to distinguish be-
Pc 1=1 Tc tween group 2 and group 3 data, only the data of Michels
et al. 80 can be assigned to group 2. While numerous authors
with Tc=304.1282 K, Pc=467.6 kg/m3, al = 1.9245108 1
a2=-0.62385555, a3=-0.32731127, a4=0.39245142, published accurate ppT data for the homogeneous gas re-
gion, even recent data sets of the saturated vapor density
t l =0.34, t2=~' t3=~' and t 4=Jt describes the data of Du-
schek et al. 58 within their experimental uncertainty. For this scatter by up to ± 0.5% in density.
equation, both the coefficients and the exponents differ from The new correlation equation for the saturated vapor den-
the formulation given by Duschek et al.,58 which is valid for sity,
temperatures on IPTS-68. The changed functional form of
the correlation resulted in slightly improved results in the In ( ~") =;.?
5
ai
( T)
1 - - f·" (3.15)
Pc ,-I Tc
critical region. According to the uncertainty of the selected
data, the uncertainty of Eq. (3.14) is estimated to be with Tc=304.1282 K, Pc=467.6 kg/m3, al = -1.7074879,
Llp' ::s±0.015% for Tt ::ST::S295 K, Llp'::S±0.04% for a2=-0.82274670, a3=-4.6008549, a4=-10.111178,
295 K< T::S 303 K, and Llp'::S ± 1% for 303 K < T::S Tc. a5= -29.742252, t[ =0.340, t2=~' t3= 1, t4 =~, and t5= ¥,
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
1526 R. SPAN AND W. WAGNER
TABLE 10. Summary of the data sets for caloric properties on the liquid-vapor phase boundary of carbon
dioxide
aThe specific isochoric heat capacity on the phase boundary cannot be measured directly; such data can only be
determined by extrapolation from the homogeneous regions.
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
A NEW EQUATION OF STATE FOR CARBON DIOXIDE 1527
authors in the past. The uncertainties which we reviewed uncertainties. Nevertheless, in the region of overlapping
resulted from extensive comparisons and should be at least a (mainly at pressures between 8 MPa and 9 MPa) the data sets
reasonable estimation. These values were used in the weight- of Duschek et(1l. 154 and Gilgen et ai. 159 show small system-
ing procedure. atic deviations of up to about ± 0.01 % in density.
Section 4.9 contains the corresponding information on all On seven isotherms between 233 K and 523 K the two-
data used for the description of the liquid-vapor equilibrium. sinker measurements are supplemented by recent
measurements l60 ,160b with a new "single-sinker"
4.1 Thermal Properties
apparatus l60,16Oe which allows measurements up to pressures
Since 1903, the temperature and pressure dependence of of 30 MPa. On principle, for densities above 100 kg/m 3, the
the density of carbon dioxide has been investigated in 59 accuracy of this apparatus is comparable with the accuracy
papers, coveririg the single-phase region with a total of 5508 of the two-sinker apparatus, but since the measurements for
data points. Table 11 presents information on the available carbon dioxide were partly made in the test phase of the new
data sets including their classifica~ion. apparatus, slightly enlarged uncertainties had to be assumed.
For temperatures of up to 360 K and pressures of up to At pressures above 13 MPa, the adjusted data sets of
13 MPa, the description of the ppT relation is based mainly Lau,149 Kiril1in pt oL 93 ,137,140 ::mcl Mlchp.ls et 01 123 (see
on the data sets of Duschek et al. 154 and Gilgen et al. 159 Sec. 4.9.2) yield a smooth continuation of the course given
These data sets are supported by the Ps' p', p" data mea- by the data of Gilgen et ai. 159 Up to 30 MPa the adjusted
sured by Duschek et al.,58 by ppT data of the gas phase sup- data are generally consistent with the new single-sinker data
plied by Guo et al. 157 and by ppT data on the 313 K isotherm within ± 0.05% in density, but the adjusted data reach up to
measured by Nowak et al. 160d All these measurements were significantly higher pressures. Due to their high quality,
performed using the "two-sinker" buoyancy method,161 these adjusted data sets are used with substantially increased
which probably provides the most accurate ppT data avail- weighting factors; details are given in Table 12. For tempera-
able today. The data from the different two-sinker appara- tures between 523 K and 698 K and pressures up to 34 MPa,
tuses are consistent with each other far within their estimated the recent data of Fenghour et ai. 160a improve the data situ-
TABLE 11. Summary of the data sets available for the p p T relation of carbon dioxide
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
1528 R. SPAN AND W. WAGNER
TABLE 11. summary of me Ul1l:1 ))Cl:s :1vailauk fvr the I'pT relation of eurbQn diOJ'i.ide- Continued
ation significantly. The data set of Straub 109 was not included 4.2.1 Experimental Results for the Specific Isobaric Heat
in the data set used for the development of the final equation, Capacity
but it was used as a sensitive test for the description of the Today, calorimetric measurements performed with flow
critical region. apparatuses provide accurate data of the specific isobaric
heat capacity over wide ranges of temperature and pressure.
In the low density region, these results are usually more ac-
curate than isochoric heat capacity measurements.
4.2 Specific Isobaric Heat Capacity Since state-of-the-art data on the speed of sound for car-
bon dioxide are available only within a very limited tempera-
As regards the specific isobaric heat capacity, the descrip- ture range and only up to pressures of 0.9 MPa, the data sets
tion of the data situation has to be split up into data sets given in Table 13 and especially the selected data sets shown
which describe the caloric behavior of the ideal gas and data in Table 14 represent the most important source of informa-
sets which in addition contain the residual behavior. While tion on the caloric behavior of carbon dioxide. While the
the data situation for the caloric properties of the ideal gas is data set of Bender et al. l77 provides an accurate description
dominated by theoretical approaches, the description of the of the low density region, the data of Ernst and Hochberg I78
real gas behavior, where both the ideal and the residual part and Ernst et al. 179 allow a precise description of caloric prop-
of the heat capacity are considered, is virtually restricted to erties up to 90 MPa. At least for subcritical pressures, the
experimental investigations. In order to account for these dif- accuracy of these data sets is improved by suitable correc-
ferent situations, this section is divided into two subsections. tions; see Sec. 4.9.3.
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
A NEW EQUATION OF STATE FOR CARBON DIOXIDE 1529
TABLE 12. Summary of selected ppT data for carbon dioxide; detailed information is given on the uncertainty values estimated by the authors and those
estimated by ourselves and used in the weighting procedure
Ylichels 9 N.R.E.a
and Michels 122 4.00 dp= 1000 Pa, dp=0.2% Lau l49.c 42 dp=O.I%
Michels et af. 123.c 75 dpb=0.05% 3.71 dT=lO mK, dp=O.02%, dp=O.I%
1.00 dp=0.02%, dp=O.l% Esperlll 45 dp=0.03-0.l1 %
WentorrB s.c 87 dT=l mK, dp=132 Pa, dp=0.02% 1.00 dT=lO mK, dp=0.015%, dp=0.05%
0.09 dT=2 mK, dp=500 Pa, dp=0.09% Holste et al. 103 80 dT=lO mK, dp=0.015%,
Vukalovich 44 dp=0.15-0.35% 0.80 dT=lO mK, dp=0.015%, dp=0.05%
and Altunin l28 1.07 dp=O.4% Jaeschke l50 245 N.R.E.
Vukalovich 134 N.R.E. 0.69 dT=5 mK, dp=O.I%, dp=O.I%
and Altunin 129 1.00 dp=0.2% Jaeschke lSI 27 NKE.
Vukalovich et af. uu 22 N.R.E. 1.14 dT=5 mK, dp=0.05%, dp=0.05%
1.00 dp=0.2% Magee 9 dT=30 mK, dp=O.Ol%, dp=O.l%
Vukalovich et al. 131 36 XR.E. and Elyl52 0.16 dT=30 mK, dp=O.Ol%, dp=O.I%
0.64 dp=O.l % Ely et af,'? 52 dT=30 mK, dp=O.Ol%, Llp=O.1-0.15%
Juza et af. IJ2 ~:2 Ap=:2 MPa, Ap=U.j% l.UU Al=jU mK, Ap=U.Ul%, Ap=U.I-U.lS%
1.00 LlT=0.2 K, Llp=2 MPa, Ap=O.4% Hoinkis l45 72 LlT=lO mK, Llp=O.Ol%, Llp=0.16%
GoJovskii 108 LlT=80 mK, Llp=0.08%, 6.p=0.15% 1.83 6.T=lO mK, dp=O.OI%, 6.p=0.16%
and Tsymamyi 136 2.25 dT=80 mK, Llp=0.08%, dp=O.15% Duschek et al. 154 362 dT=3 mK, dp=0.006%, dp=0.015%
Kililliu t<[ ul.1"l7c 17 Af}~O.I% 1.88 AT-3 lllK, Ap-O.006%, Ap-O.015%
9.65 dp=0.05%, dp=O.2% Duschek et al. 58 87 Ap=0.03-0.25%d
Kirillin et af. 93.c 29 N.R.E. 2.25 6.p=0.03-0.25%
4.00 Llp=0.05%, dp=0.2% Iaeschke et af. 155 268 N.R.E.
0.75 6.T=5 mK, Llp=0.05-0.1%,
Kirillin et af. 138 66 6.p=0.15-0.2% Guo et af. 157 40 6.p=0.05-0.01 %
1.00 6.p=0.05%,6.p=0.2% 2.25 dT=3 mK, dp=0.005%, dp=O.OI%
Tsiklis et af. 139 49 Llp=O.3% Gilgen et aZ. 159 264 6.T=3 mK, Llp=0.005%, 6.p=O.Ol %
1.00 6. T=0.5 K, 6.p = I MPa, 6.p= 1.0% 2.18 AT=1.5 mK, dp=0.006%, dp=0.015%
Kirillin et al. 14O.c 24 Llp=O.l% Brachthauser l60 29 LlT=1.5 mK, dp=0.006%, dp=0.015%
4.00 Llp=0.05%,6.p=0.1% 1.00 6.p=0.016-0.051 %
Popov 73 N.R.E. Fenghour et af. 16Oa 120 6.p=0.016-0.051 %
and Sayapovl~1 1.00 LlT=30 mK, Llp=0.05%, Llp=0.15% 1.00 Llp=0.05-0.1O%
Vukalovich et af. 142 86 Llp =0.08-0.15% Klimeck et af. 16(Jb 60 Llp=O.lO%
1.00 AT= 15 mK, Ap =0.005%, 1.00 dT=6 mK, Llp=0.006%, Ap=0.02-0.04%
Ap=0.08-0.15% 6. T=6 mK, Llp =0.006%, Llp =0.02-0.04%
Shmonov 60 dp=0.25-0.5% Nowak et af. 160d 21 dT=1.5 mK, dp=0.006%, dp=0.015%
and Shmulovich 148 1.44 AT=O.4 K, Llp =0.2%, Llp= 1.0% 1.00 AT=1.5 mK, Ap=0.006%, dp=0.015%
The data sets of Rivkin and GUkov l 71.174 describe the su- this data set was classified in group 3. The data set given in
percritical region of carbon dioxide. Since empirical equa- Ref. 171 was used only with a reduced weighting factor.
tions of state may produce misleading results regarding ca-
loric properties within this region, these data sets were
considered to be of great interest for the development of the
new fundamental equation. Unfortunately, the data proved to 4.2.2 Results for the Specific Isobaric Heat Capacity in the
be inconsistent with each other and with state-of-the-art ppT Ideal-Gas State
data. Since the data given in Ref. 174 were deduced from According to Eq. (2.4), knowledge of the specific isobaric
measurements of a mixture with high carbon dioxide content, heat capacity in the ideal-gas state forms the basis for the
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
1530 R. SPAN AND W. WAGNER
TABLE 13. Summary of the data sets available for the specific isobaric heat capacity of carbon dioxide
aThe paper also contains extrapolated values for the specific isobaric heat capacity of carbon dioxide in the ideal
gas stare. .
bData calculated from measurements of the p p T relation.
cMeasuremems of a mixture with high carbon dioxide content; the results considered here were extrapolated
from the mixture experiments to pure carbon dioxide.
description of the ideal-gas part of the Helmholtz energy, only a few papers contain reliable data on c~, where correc-
AO(p,T). Results for c~ obtained by the extrapolation of cp tions to the simple rigid rotator, harmonic-oscillator model
measuremeIllS (0 zero pressure (cf. Table 13) are neirher ac- were taken into account in order to achieve high accuracy.
curate enough for this application nor do they cover a suffi-
Table 15 summarizes these data sets, of which the data set of
ciently wide range of temperature. Therefore, theoretical ap-
Chao l84 is the most recent one. Chao considers first order
proaches which describe the caloric behavior of carbon
corrections LO Lhe rigid roraror, hannonlc-osc1l1aror model,
dioxide in the ideal-gas state were reviewed.
which were developed by Pennington and Kobe; 185 he esti-
Information on the fundamental frequencies of the carbon
dioxide molecule can be found in various compilations, but mates that the uncertainty of his results is less than
±O.02%.
TABLE 14. Summary of selected data for the specific isobaric heat capacity
of carbon dioxide; detailed information is given on the uncertainty values
estimated by the authors and those estimated by ourselves and used in the
weighting procedure
4.3 Specific Isochoric Heat Capacity
Uncertainty estimated
by the authors Seven data sets are available for the isochoric heat capac-
No. of data Uncertainty estimated ity of carbon dioxide. The information on these data sets is
Source Meanf~l by ourselves
summarized in Table 16. More detailed information on the
Kistiakowsky ~cp==0.3o/c
and Rice l63 1.00 ~T==30 mK, Ac p =0.3o/c
Masi 12 ~T==30 mK, ACp=O.lo/c
and Petkof l65 1.00 ~T==30 mK, ~cp=O.I%
Rivkin 220 AT= 10 mK: Ap=0.OSC70,Acp==2% TABLE 15. Data sets for the isobaric heat capacity in the ideal-gas state of
and GUkOV l71 0.17 I:1T==25 mK, I:1p=0.05%, AC p ==2% carbon dioxide calculated by theoretical approaches
Bender e1 aJ. I 77.o 60 ~C p ==0.1 %-0. 15%
1.00 I:1T=lO mK, Ap==O.OI%, AC p ==0.12% Source Year Temp. range, T (K)
Ernst 9 ~cp=G.2o/c-O.99C
and Hochberg l78 .b 2.25 AT===20 mK, 6.p===O.I'K,Ac p ==0.2-0.9% WOOl ey l80 1954 50-5000
Ernst eC ai. 179 .h 55 6.c p:=:0.2'K-0.9% Baehr el al. 181 1968 10-6000
1.92 IJ.T==20 mK, IJ.p=O.l9'c, AC p ===0.2-0.99'c Gurvich lS2 1979 10-10 GOO
Chao l83 1983 50-5 000
"The consistency is given instead of the uncertainty. Chao l8 -l 1986 10- t 500
bCorrected values were used in the final data set, see Sec. 4.10.
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
A NEW EQUATION OF STATE FOR CARBON DIOXIDE 1531
TABLE 16. Summary of the data sets available for the specific isochoric heat capacity of carbon dioxide
selected data is given in Table 17. For carbon dioxide, mea- within less than 0.5% by the new equation of state. Based on
surements of the isochoric heat capacity are of great impor- these facts, we concluded that deviations up to 8% in the
tance in the two regions discussed below. Isochonc heat capaCIty were probably due to shortcommgs III
The data set of Magee and Elyl16 yields the only compre- the data set. Thus, the weighting factors were strongly re-
hensive description of caloric properties within the high den- duced for the low density measurements.
sity region. At liquid densities, this data set is limited to Within the extended critical region, the data sets of
pressures above the critical pressure, leaving a gap between Edwards 55 and of Abdulagatov et al. 59,59a.59b had great influ-
the available caloric data in the homogeneous region and on ence on the development of the new equation of state. (Ed-
the saturated liquid curve (see Sec. 7.2.2). At densities below wards's thesis 55 only contains graphical illustrations of the
600 kg/m 3, these data can only be represented with system- results. The data used here were provided by Sengers.189)
atic deviations clearly outside the uncertainty estimated by The reduced weighting factors given in Table 17 were cho-
the authors. On the other hand, in this region several accurate sen because of the large number of data points within a small
equations of state 3,15-17 agree with each other and with the and very sensitive region. The use of these data sets with
new equation of state within about 2% for the isochoric heat f wt= 1 would have resulted in an overemphasis of the critical
capacity, and data of the isobaric heat capacity are met region; details are given in Sec. 5.
The recent data of Abdulagatov et al. 59 were first pre-
sented in the Proceedings of the 11 th Symposium on Ther-
TABLE 17. Summary of selected data for the specific isochoric heat capacity
mophysical Properties in Boulder in 1991. When the corre-
of carbon dioxide; detailed information is given on the uncertainty values
estimated by the authors and those estimated by ourselves and used in the sponding paper 9a was published three years later, it
weighting procedure contained only 88 data points selected from the 331 data
points given in the proceedings. In the same year, another
Uncertainty estimated paper59b was published which contains the 88 points pub-
by the authors
No. of data Uncertainty estimated
lished before, 10 additional points on one isochore in the
Source Mean f~.! by ourselves homogeneous phase, and 132 additional points in the two
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
1532 R. SPAN AND W. WAGNER
TABLE 19. Summary of selected data for the speeds of sound of carbon TABLE 21. Summary of the data sets for differences of the internal energy of
dioxide; detailed information is given on the uncertainty values estimated by carbon dioxide
the authors and those estimated by ourselves and used in the weighting
procedure No. of Temp. range, Density range,
Source Year data T(K) p (kg/m3) Group
Uncertainty estimated
by the authors
Kriiger204 1964 454a 274-325 190-870 1-2
Baeh~o5 1968 47 291-316 466
No. of data Uncertainty estimated
Source Meanl"wt by ourselves
aSome of the measured temperature intervals partly reach into the two-phase
Novikov 234 N.R.E.a region.
and Trelin l92 0.77 LlT=O.1 K, Llp=O.I%, Llw=0.5%
Pitaevskaya 144 N.R.E.a
and Bilevich 193 2.25 LlT=O.1 K, Llp=0.2%, Llw=2%
Lemming 194 44 LlT=8 mK, Llp=4.2 mbar, Llw=0.015%
Besides the data presented in Table 18, five papersl95-199
0.06 LlT=8 mK, Llp=4.2 mbar, Llw=0.OI5% are available which deal with the speed of sound in the im-
mediate vicinity of the critical point. None of these papers,
aN.R.E.: No reasonable estimation given by the author.
which were published between 1951 and 1970, contains nu-
merical results. We did not try to evaluate the graphical. rep-
resentations given in these papers since the special problems
phase region. The uncertainty of the data of Abdulagatov
connected with experimental work in the critical region were
et ai. was estimated to be at least as big as the scatter visible
not adequately considered at that time.
in the original data set. 59
The papers of Lipa et ai., 190,52 which are often cited in
literature on carbon dioxide, do not contain numerical re-
4.5 Enthalpy
sults. Since the dutu of Edwards 55 cover the sume region,
were measured with the same equipment, and are less influ- The literature sources which deal with enthalpy measure-
enced by impurities,55 the data of Lipa et al. were not con- ments are summarized in Table.20. The data of Vukalovich
sidered in this work. and Masalov 201 ,202 and of Moller et ai. 203 were considered in
this work, but only the data of Moller et al. were used with
4.4 Speed of Sound the uncertainties given by the authors (10 data selected,
mean value of f~t=2.78, b"h=0.7% or b"h=0.6 kJlkg,
Three of the four data sets available for the speed of sound whichever is greater). These data provide important informa-
of carbon dioxide were used for the development of the new tion at high densities where the data situation for caloric
equation of state. Table 18 contains the available data sets properties is poor. The data of Vukalovich and
and Table 19 gives additional information on the selected Masalov 2ol ,202 were used to test the quality of the new equa-
data. The data sets of Novikov and Trelin 192 and of Pitaev- tion with regard to caloric properties at high temperatures.
skaya and Bilevich 193 provide important information on the
extended critical region and the high pressure region, respec-
tively. Though the quality of the data published by 4.6 Internal Energy
Lemrning 194 is superior to the other data sets, these data were
used only with reduced weighting factors. In the gas region, Table 21 lists the general information on the two data sets
the residual part of the speed of sound is less than 2% of the which are available for differences of internal energy. The
total speed of sound at pressures below 1 MPa. Since the data measured on the critical isochore by Baehr20J were not
uncertainty of measurements is related to the total value of used since the temperature displacement of the (aul aT) L'
the speed of sound and not only to the residual part, the maximum observed by Baehr could not be explained. The
contribution of thc:o.c data fOl the Lievdut'llleIll u[ all empiri- 150 selected data of Kriiger,204 which do not reach into the
cal correlation describing the residual fluid behavior is very two-phase region, were used with estimated uncertainties of
limited. b"T=40 mK, b"p=0.2%, b"u=2%-4%, and a mean weight-
ing factor of l;'t=0.88. The uncertainties given by the author
(~T= 10 mK, b"p =0.03%, b"u= 1.4%-3%) are too small to
TABLE 20. Summary of the data sets available for enthalpy differences of explain the scatter of the data.
carbon dioxide
Temp. Pressure
No. of range, range,
4.7 Joule-Thomson Coefficient
Source Year data T (K) p (MPa) Group
Since data of the Joule-Thomson coefficient can only be
Maas and Barnes 20n 1926 4 212-298 0.1 used in the nonlinear fitting procedure and not in the linear
Koppel and Smith l6i 1960 102 291-322 7.2-8.3 algorithms of the optimization procedure, their influence on
Vukalovich and Masalov 201 1964 54 423-523 2.5-9.8
the development of an empirical fundamental equation is
Vukalovich and Masalov 202 1964 68 573-773 2.5-9.8
Moller et al. 203 1993 10 233-358 15.5 fairly small. The data of Bender et al. 177 are the only ones
which were included in the data set used for the nonlinear fit
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
A NEW EQUATION OF STATE FOR CARBON DIOXIDE 1533
TABLE 22. Summary of the data sets available for the Joule-Thomson coefficient of carbon dioxide
(34 tiala Wt;It; selected with 6.T= 10 mK, 6.p=O.Ol %, Recently, Span225 pointed out that it was not very useful to
~p,=0.4%, and f~t=0.36). The data of Vukalovich et al. 206 include experimental second and third virial coefficients into
proved to be a sensitive test for the consistent description of the data set used for the development of a wide-range equa-
the B:l~ region. Tahle 22 gives a survey of the data sets deal- tion. Most of the data sets descend from ppT measurements
ing with the Joule-Thomson coefficient. which have been evaluated in order to determine virial coef-
ficients. Thus, the use of the original measurements yields
much better access to the desired experimental information.
4.8 Virial Coefficients Furthermore, fitting an equation of state to virial coefficients
is only useful if the terms of the equation corresponding to
The literature on carbon dioxide contains many papers in the virial coefficients are independent of each other, as is the
which information on the second and third vitial coefficient case in a simple virial expansion. However, if an equation of
is given. Table 23 summarizes the corresponding sources. state contains exponential functions, this condition is no
longer met. Thus, in this article, the available values for
virial coefficients were only used for comparison; none of
TABLE 23. Summary of the data sets available for the second and third virial
coefficient of carbon dioxide. For reasons explained in the text, no data were
the data sets was assigned to group 1.
assigned to Group I At temperatures below 220 K, however, the whole contri-
bution of the third virial coefficient [p / (pR T)
No. of data Temp. range, = ... + C p2+ ... ] is smaller than the uncertainty of the
Source Year S/C T(K) Group
most recent ppT data l54 throughout the gas region. Therefore,
Michels and Michels122 1935 13/13 273-423 3/3 at low temperatures a physically unreasonable representation
Schafe?09 1937 16/· .. 203-273 3/··· of the third virial coefficient calculated from an equation of
MacCormack and Schneide?1O 1950 9/9 273-873 3/3 state may occur if the equation was fitted only to ppT data.
CoHre! and Hamilton 211 1955 7/· .. 303-333 3/ .. ·
Pfeff':!rJe et al. 212 1955 212
Since reliable value~. of the th~rd virial coefficient cannot be
303 3/3
Cottrell et al. 213 1956 3/ .. · 303-363 3/ .. · established by an evaluation of experimental ppT data in this
Cook214 1957 6/· .. 213-303 3/ .. · temperature range, 13 values of the third virial coefficient
Masia and Pena llS 1958 6/ .. · 298-423 3/· .. were calculated from a simple polynomial equation which
Butcher and Dadson 216 1963 13/13 263-473 3/3 describes all the selected data of the different thermody-
Huff and Reed 217 1963 10/ .. · 298-510 3/···
Vukalovich and Masalov 218 1966 18/18 423-773 2/3
namic properties in the gas region within their experimental
Dadson et al. 219 1967 9/· .. 263-398 3/··· uncenainty and yields a reasonable plot of the lhird virial
Ku and Dodge 133 1967 III 373 3/3 coefficient at low temperatures. These "artificial" data of
Sass et af.I:I... a 1967 3/··· 348-398 3/·· . the third virial coefficient were then used during the devel-
Timoshenko et al. 22o 1970 9/· .. 224-313 3/· .. opment of the new equation of state.
Yukalovich et af. J4l 1970 10110 253-343 3/3
Vukalovich et al. 206 1970 14/ .. · 238-308 3/···
Waxman et al. 221 1973 6/ .. · 273-423 2/· ..
Bender et af.177 1981 4/··· 233-263 3/·· .
Ohgaki et al. 222 1984 2/ .. · 423-473 3/· .. 4.9 Liquid-Vapor Equilibrium
Holste et af. 1~3 1987 18/16 217-448 2/2
Mallu et al. 223 1987 3/··· 323-423 3/ .. · During the procedure of optimizing the structure of the
Hoinkisl~5 1989 4/4 298-423 2/2 new equation of state, the liquid-vapor equilibrium was used
Mallu et af. 22~.h 1989 .. ·/3 323-423 · .. /3 in a linearized way.32,30 The data set used for this purpose
Me Elroy et al.15:' 1989 4/4 303-333 3/3 consists of values of Ps(1), pi (T), and pl/(1) at 205 tempera-
Duschek et af. 15" 1990 7/4 220-340 2/2
tures which were calculated from Eqs. (3.13) to (3.15). These
"Yirial expansion developed in pressure. data cover the whole liquid-vapor phase boundary with tem-
hThis paper additionally includes the values already published!23 in 1987. perature intervals which decrease when approaching the
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
1534 R. SPAN AND W. WAGNER
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
A NEW EQUATION OF STATE FOR CARBON DIOXIDE 1535
TABLE 25. Temperature dependent corrections of isobaric heat capacity data 4.10.3 Correction of Isobaric Heat Capacities
c p.cor. - c p.exp.
Temperature, Since the low pressure limit of the isobaric heat capacity,
T(K) cp.e,p. c/' is known very well for carbon dioxide (see Sec. 6.1),
Bender et at. 177 temperature dependent errors of measured isobaric heat ca-
233 +0.40% pacities can be determined easily. Figure 7 shows systematic
243 +0.25% deviations of data measured by Bender et al. 177 and by Ernst
273 +0.13% et al. 179 Again, two different isotherms are combined in a
323 +0.05%
+0.31%
single deviation plot. A systematic error of the new equation
363
393 +0.l4% for the residual part of the Helmholtz energy, Eq. (6.5),
423- +0.l0% would result in deviations which increase with pressure since
473 -0.17% the residual contribution to the isobaric heat capacity in-
l78
Ernst and Hochberg creases.
303 +0.36%
To compensate for these temperature-dependent devia-
Ernst et al. 179
JJJ 0.18% tions, the data sets of Bender P.t a/., 177 Ernst ~nrl
363 :::1::0.00% Hochberg,178 and Ernst et al. 179 were corrected according to
393 -0.37% the values given in Table 25.
TABLE 26. Examples for power laws describing thermodynamic properties along certain paths throughout the critical region
Densities at saturation (p' - p")~(Tc- T)f3 phase boundary f3 0.326:::1::0.002 0.5 0.25
Isothermal compressibility KT~IT-Tcl-Y crit. isochore ')' 1.239:::1::0.002 1 J
Pressure ;p-Pcl-Ip-Pcl" crit. isotherm Ii 4.80:::1::0.02 3 5
Isochoric heat capacity cu-IT-Tcl- a crit. isochore a 0.110:::1::0.003 0 0
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
1536 R. SPAN AND W. WAGNER
;005 304.1282 K 6
Chen etal.
27
E
02-type equation
1
5 '1-'-'- CH4-type equation
g 0 _,._
0 1
...... g ~ 1 0
1 D
Edwards
55
FIG. 8. Relative deviations 100 !::.P = 100 (Pexp - Pcalc)/P exp of experimen-
1'------'---'---'-----'---'---'----'---'-----''----'
tal ppT data on the critical isotherm from values calculated from Eq. (6.1). 300 302 304 306 308 310
Values calculated from the crossover equation of Chen et aZY and from
refitted equations using the CH4 - and OTform (see Sec. 5.1) are plotted for Temperature T / K
comparison.
FIG. 9. Representation of representative isochoric heat capacity data in the
critical region. The plotted curves correspond to values calculated from the
rrm:<:nvpr f'qu::Ition of Ch~n I?t (I/. 27 and from refitted equations using the
CH4 - and 02-form (see Sec. 5.1).
5.1 Limitations of Analytical Equations of State
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
A NEW EQUATION OF STATE FOR CARBON DIOXIDE 1537
capacity. From a theoretical point of view, the speed of scription of the ppT surface within the critical region should
sound is expected to vanish at the critical point, whereas also yield reliable values of the specific isobaric heat capac-
analytic equations yield only a finite minimum. Again, the ity.
temperature interval, where the analytic equations fail, ex- The situation is different for the isochoric specific heat
tends over approximately ±O.5 K and ±2 K, respectively, capacity which is given by
around the critical temperature.
(5.2)
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
1538 R. SPAN AND W. WAGNER
1.80 6
T=305 K
~ 1.78 _._0_.-'-'_'_°,0-.-.-
I
0-._,_.-
Chen et at. 27
(0
N~ 1.76 I Eq. (6.1)
e 1.74 5 CH4 -type equation
"'~ 1.72 ---
---
~ 0 Edwards 55
~ 1.70
0.0 4 0 Abdulagatov et at. S9
~
.:
......
I '-"
~ .
I 3
I
~ I P=P c I
2
467.60 467.65
Density p I (kg/m3)
Eq. (6.1) -.-.- Preliminary eq.
Temperature T I K
FIG. 10. While preliminary equations of state showed a discontinuous plot
of the third density derivative of the reduced Helmholtz energy, the new FIG. 11. Representation of representative isochoric heat capacity data in the
equation of state, Eq. (6.1), yields continuous plots for the third derivatives. critical region. The plotted curves correspond to values calculated from Eq.
(6.1), the crossover equation of Chen et aI., 27 and a refitted equation using
the CH4-form (see Sec. 5.1).
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
A NEW EQUATION OF STATE FOR CARBON DIOXIDE 1539
20~~-~--------------------~
6....----------------;.....- - - - - - - - - ,
10
I .
0.5
3 1..----------, 1....
0·-:-6...........L..I..I.IJ1w,0-:·s:-'--'....L.LIJJ.J1O"-.-:-4.L....L.J...L.LIlwO-:·3:-'---'L...L.LJWl1L0·-2...L......L..J...LLWJ1O' I
T= Tc + 0.5 K IT - Tel I Te
, " ',,"I! "II""! II",," ",,'"'' II''''''
0.0003 0.003 0.03 0.3 3 30
2 IT - T) I K
55
Eq. (6.1) o Edwards
27
Chen et al. o Abdulagatov et ai.•u
Pc 02-type equation
1~~~-L-~~-~~~~~~-~~
The final equation of state for carbon dioxide, which is - - Eq. (6.1)
presented in detail in Sec. 6, has been developed by the use
of nonanalytic terms corresponding to Eq. (5.4). The proce- 2000
dure optimizing the structure of the equation of state was ~
-
OJ)
restricted to a maximum of four nonanalytic terms within a 1000
C
total number of 42 terms. Preliminary correlations with more -.
than four nonanalytic terms tended to have an unreasonable
behavior concerning the dependence of Cu and w on density
in the critical region; moreover, when using more than four
-
~
0
\.)~
500
of such complex terms, the numerical expense would have 'uro 200
0..
been increased without a significant improvement in the ro
u
quality of the equation. Anticipating Sec. 6 with regard to the 100
~
Q)
critical region, Eq. (6.1) will be discussed in this section to ...0
U
avoid repeating a discussion of the representation of caloric 50
properties in the critical region later in this paper.
Figure II shows the plot of C u on the critical isochore; for
.~
on
~
0 )IillJ
analytic equations of state this plot was already shown in 20
Fig. 9. In this figure, however, the solid line corresponds to Pc
Eq. (6.1). In contrast to the analytic CH 4-type equation, the to
350 400 450 500 550 600
new formulation is able to follow the strong curvature of the 3
Cz; plot in the immediate vicinity of the critical temperature
Density pi (kg/m )
and yields an infinite value for the specific isochoric heat FIG. 15. For temperatures between Tc and T, + 1 K Eq. (6.1) results in an
capacity at the critical point. Accordingly, the evaluation of oscillating plot for the isobaric heat capacity around the critical density.
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
1540 R. SPAN AND W. WAGNER
TABLE 27. Coefficients of the correlation equations, Eq. (6.2) and Eq. (6.3), cal point and which is still explicit in the physical variables T
for c~ and cPo, respectively
and p. In Sec. 7 it will be shown that this nonanalytic be-
havior does not affect the quality of the equation with respect
to the representation of other properties anywhere outside the
8.373 044 56 0.62105248 6.11190 critical region.
-3.70454304 6 0.411 95293 6.777 08
Nevertheless, empirical equations of state containing
2.50000000 7 1.040 289 22 11.323 84
4 1.994 27042 3.151 63 0.083 27678 27.087 92 non analytic terms have certain limits if an exact fulfillment
of the asymptotic power laws is required. The nonanalytic
terms in Eq. (6.1) do not replace the contribution of the
analytic terms in the surrounding of the critical point but
similar values for the speed of sound. However, on ap- they fill the increasing gap between the analytic and nonana-
proaching the critical temperature, both the crossover equa- lytic behavior. Therefore, efficient values of the exponent
tion and the nonanalytic wide-range equation, Eq. (6.1), de- hi are smaller than theoretically expected and result in a criti-
velop a sharp minimum in the speed of sound, whereas the cal exponent a which is too large from an asymptotic point
speed of sound calculated from the analytic equation does of view [see Eq. (5.5)].
not change significantly when approaching the critical tem- Figure 14 shows that Eq. (6.1) represents the measured
perature. Within the same density range, Fig. 13 shows the and theoretically predicted values of Cv without significant
corresponding plots of the specific isochoric heat capacity deviations for temperatures in the region IT-Tel> 0.2 K and
calculated from the same set of equations. withm about :i: 10% tor U.2 K~(1'-1~)~O.3 mK in the ho-
Thus, Eq. (6.1) is the first wide-range equation which mogeneous phase and for -0.2 K::::;(T-TJ::::; -10 mK in the
yields a nonanalytic behavior of the isochoric heat capacity two-phase region. The plotted Cv courses correspond to the
and the speed of sound in the immerliilte vicinity of the criti- critical isochore in the one- (T> TJ and t,vo-phase (T< Tc)
0.1 r-----------------~----~----------------------_r------------_.
I \ I
I I [J
c
c x x x
0..
*
><
o Ill.
'i..Jr::.. • •
---
".-...,
~ -0.1L---~--~~~---L--~--~--~--~~--~--~--~--~---L--~--~
0'i..J~
U 0 300 600 900 1200 1500
10.. 0.5 ~----------------------------------------------------------~
><
<I.)
0'i..J~ -A
'-" [!]
oo r> .-.-.-.-.
.,.-.--.-
~ - - - - -
[!]
[!]
<l <l
c;]<1 <l
v
-0.5 L....-----J--+4------L_-----'-_---l..._---L.._---l...._--'-_---L...._--I....-_---L..-_---'--_
FIG. 16. Relative de\'iations of c~ data from values calculated from Eq. (6.2). The upper diagram shows data calculated from statistical thermodynamics (see
Table 15) and the lower diagram shows data extrapolated from experimental results (see Table 13). Values of c~ calculated from the corresponding equations
of Angus et af.·l and Ely et af.ls are plotted for comparison.
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
A NEW EQUATION OF STATE FOR CARBON DIOXIDE 1541
TABLE 28. The ideal-gas part of the dimensionless Helmholtz function c/J0 and its derivatives a
8
c/J0 In 8 + aOI + a~r + a~ In(r) + L a~ln(1-e-e:T)
i=4
c/J~ 1/8 + 0 + 0 + 0 + 0
c/J 8S -lIU + 0 + 0 + 0 + 0
0
c/J~T 0 + 0 + 0 + 0 +
8
region. However, within a region of about +0.1 mK (outside where 8 = pI Pc and r = TclT with Pc = 467.6 kg/m3 and
the temperature range shown in Fig. 14) and - 10 mK around Tc = 304.1282K.
the critical temperature, Eq. (6.1) yields specific isochoric The formulations which describe the ideal-gas part of the
heat capacities which are significantly larger than the values Helmholtz energy, Eq. (6.3), and the residual part of the
predicted by the crossover equation of state. Helmholtz energy, Eq. (6.5), are introduced in this section.
Figure 15 shows a problem which is related to a small
oscillation in the derivative (op/ op h close to the critical
point. Since the numerical value of this derivative is small 6.1 Ideal-Gas Part of the Helmholtz Energy
close to the critical point, even a very small oscillation in this
According to Eq. (2.4), the ideal-gas part of the Helmholtz
derivative results in a significant oscillation of properties re-
energy can be easily obtained if the function c~(T) is known.
lated to its reciprocal value like the isobaric heat capacity
A correlation equation for c~(T) has been established by
[see Eq. (5.1)] or the isothermal compressibility. Derivatives
means of a nonlinear fitting routine using the 150 c~ data of
of these properties [e.g., (ocp/oph] should not be used in
the range 440 kg/m3:::;p:::;500 kg/m3 and (T-Tc):::;2 K. Os- Chao l84 as input values (cf. Sec. 4.2.2). With the coefficients
given in Table 27, the obtained equation
cillations are also observed in the derivative (oc v / oT) p
for approximately 420 kg/m3:::;p:::;550 kg/m3 and
(T-Tc):::; 15 K. TABLE 29. Summary of selected data used for the linear optimization pro-
cedure and for the nonlinear fit
The quality of Eq. (6.1) should be sufficient for any appli-
cation which requires numerical values of thermodynamic Number of data
properties in the critical region. Equation (6.1), however, is Details of the
data set are Linear Nonlinear
not suitable to investigate the asymptotic behavio{ of ther-
Property given in. optimization fit
modynamic properties and the derivatives mentioned above
should not be used in the regions indicated. p(T,p) Table 12 2824 2824
cu(T,p) Table 14 553 553
a
cp(T,p) Table 17 359 359
c;(T) Table 24 7Y 73
6. The New Equation of State w(T,p) Table 19 422" 422
w'(T) Table 24 2~" ?l
w"(T) Table 24 18" 18
As discussed in Sec. 2.1 the new equation of state for !1h(T,p) Sec. 4.5 10" 10
carbon dioxide is a fundamental equation expressed in form !1u(T,p) Sec. 4.6 ISO ISO
of the Helmholtz energy ~" J.L(T,p) Sec. 4.7 34
p(T,p') 205 b
A (p, T)/(RT) = ¢( 8, r) = ¢O( 8, r) + ¢r( 8, r), (6.1) p(T,p") 20S b
Maxw.-crit. 205 b
p,(T) Table 7 109
CThe recent experimental results of Wagner et al. 230 for the "thermal" criti- p'(T) Table 8 50
cal exponents f3. y. and 8. derived from ppT measurements in the imme- p"(T) Table 9 42
diate vicinity of the critical point show clear differences from the corre-
sponding values predicted by the renormalization group theory. These
suprising results might be caused by an extended validity range of the Total number of data points 5047 4667
so-called explicit influence of gravity (the implicit influence of gravity, e.g.,
the averaging errors based on density stratifications. was taken into account aLinearized data used in the optimization procedure; see Setzmann and
when evaluating the experimental ppT data). Thus. we think that the Wagner. 30
"true" asymtotic behavior at the gas-liquid critical point of a pure fluid is bLinearized solution of the Maxwell criterion when usin2: data calculated
not yet finally clarified. from the auxiliary Eqs. (3.13) to (3.15); see Wagner. 32 ~
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
1542 R. SPAN AND W. WAGNER
103 .-------------~~----------------------------------~
cPo~
d'0o
~ I ~ ~ gig ~ ~ In~L~g III
0
0
§
g e 0
0 g e
0
0
0
0 0
0 0 0 0 0 0 0 0
0
0
0
0
0
0
0
~OO~8 ~H~!~
g§gO§8ggg
t
00
<§D
00
~
()
0
00000.0.0
0
888000
0 0 0.0 0
o § 0
o 0 0
o § 0
o
o ppTdata
10- 1 ~_________0~~~____~____~____~____L -__~__L-~_____
200 300 400 500 600 800 1000
103 =---------------~--------------------------------------~
-.....
o c, data
• u data
o cp data
QQ o o..tl o 0 0 o
• h data
\1 f-1 data
QO o 0 o 0 0 o wdata L::..
o o 0 o
10- 1 L -_ _~~_G~~~~~~&__£-L----~---L---L--~~----~
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
A NEW EQUATION OF STATE FOR CARBON DIOXIDE 1543
cO(T) 8 exp( efr) discussed in Sec. 2.3. The selected data which form the ex-
-p--= 1 +a~+ ~ af(efr)2 [ (eO )-1]2 perimental basis of the new equation of state for carbon di-
R i=4 exp i r
oxide have been presented in Sec. 4. Table 29 gives a brief
(6.2)
summary of the data used and refers to the corresponding
represents Chao's data with deviations less than ± 0.005% tables, where detailed information is given. Figure 17 shows
for 10 K~T<21O K and less than ±0.002% for the distribution of the experimental data used. In addition,
210 K~ T~ 1500 K. This means that the uncertainty of Eq. data have been calculated from auxiliary equations and from
(6.2) is almost equal to the uncertainty of the data which was preliminary equations of state in order to guarantee reason-
estimated by Chao to be less than ± 0.02%. Extrapolation of able behavior in regions where the existing measurements
Eq. (6.2) yields results which are reliable to ± 1 % for yield insufficient information. In detail, these are
T~ 4000 K and to ± 2 % for 4000 K < T~ 5000 K. Figure 16
shows deviations between values calculated from the
relations given by Angus et al. 3 and by Ely et at., I data
cor- c; • 13 values of the third virial coefficient, which have been
calculated from an auxiliary equation (see Sec. 4.8),
• 23 values of the specific isochoric heat capacity, which
obtained by theoretical approaches and by extrapolation of have been calculated from the crossover equation of Chen
ex Pp.ri mp.n t::l 1 rp..;mlts, ::Inn v::Ill1es c::Ilcl11ated from Eq. (6.2). et al. 27 in order to guarantee a reasonable dependence of
which correspond to the zero line. c v on density in the critical region,
The expression for cpO( 8, r), which has been derived from • 27 Tp points describing the course of the Joule curve,
Eq. (6.2) by integration, is which have been determined by graphical extrapolation
</J
0
( 8, r) = In( 8) + a~ + a~ r+ a3 In( r) (see Sec. 7.3.2),
8 • and 70 ppT data, which have been calculated from a pre-
+L
;=4
a?ln[l-exp(-re?)]. (6.3) liminary equation of state with an exceptionally good ex-
LrapulaLiull ut;haviul ill ullkr Lv prt;~t;lVt; Lhi~ prugrt;~~ (~t;t;
The coefficients a? and e? are given in Table 27. The Sec. 7.3.1).
coefficients a~ and a~ were adjusted to give zero for the ideal All these "artificial" data are neither considered in Table 29
gas enthalpy at To = 298.15 K and the ideal gas entropy at nor in Fig. 17. The consequences of the absence of caloric
T o=298.15 K and Po = 0.101 325 MPa. In Table 28, all de- data within the liquid region and at temperatures above 500
rivatives of the ideal-gas part cpo required for the calculation K are discussed in Sec. 7.2.2.
of thermodynamic properties are explicitly given. The bank of terms which was used in the optimization of
the final equation of state (see Sec. 2.3.2) contained a total of
6.2 Residual Part of the Helmholtz Energy
860 terms. As far as parameter ranges are predetermined by
The formulation for the residual part of the Helmholtz this general form, these ranges have been established in ex-
energy has been developed with the help of the procedure tensive tests. This bank of terms can be written as
4 20 6 10 8 8 8 16 10 12
cpr=?
1=1 J=O
4 n;,j 8 i r jl4+ 4 4 1=1 J=O
n;,j 8 i r jI2e-8+ 44
1=1 J=O
n;,j 8 i r j e-
82
+~ 4 !l;,j 8 i r j e-. 83 + ~
I-I J=O
~
1-1 J-O
ni,j 8 i~je-04
10 16 15 16 48
+" "
£..; £..;
n· .
I,J
8 i~je-b·5 + "~ "~ n· .
I.J
8 i~je-b<5 + "~ n.
I
8 dirie-CXi(O-fi)2_f3i(T-Yi)2
;= 1 j=5 1=8 J=5 1= 1
2 2 3 3
(h4)
3.00 0.875 0.30 10.00 225.0 0.700 0.300 used for the parameters of the 108 nonanalytic terms in Eq.
3.50 0.925 1.00 12.50 250.0 (6.1) are given in Table 30. The values used for b j' C I' and
4.00 15.00 275.0
D In resulted from a nonlinear fit of a preliminary formula-
"Predetermined from a simultaneous fit to saturated liquid and vapor densi- tion.
ties in the critical region. Proceeding from the bank of terms defined by Eq. (6.4),
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
1544 R. SPAN AND W. WAGNER
t;
t;
1
8 0.216589615432 20X 10 1.50
9 0.158417351097 24X 10 1 2 1.50
10 -0.231327054055 03x 10° 4 2.50
11 0.581 169164 314 36X 10- 1 5 0.00
12 -0.553691372 053 82X 10° 1.50
13 0.489 466 159 094 22 x 100 5 2.00
14 -0.24275739843501 x 10- 1 6 0.00
15 0.624947905016 78X 10- 1 6 1.00
10 -U.11.1 "/)) {)U:l'.l)140X IUo b 2.UU
17 -0.370556852700 86x 10° 1 3.00
18 -0.167758797004 26x 10- 1 1 6.00
19 -0.119607366379 87X 10° 4 3.00 2
20 -0.456193625087 78x 10- 1 4 6.00 2
21 0.35612789270346XlO- 1 4 8.00 2
22 -0.744 277271 320 52x 10- 2 7 6.00 2
23 -0.173957049024 32X 10- 2 0.00
24 -0.21810121289527XIO- 1 2 7.00
25 0.243321 665592 36X 10- 1 3 12.00
26 -0.374401 334 234 63 x 10- 1 3 16.00
27 0.l43 387 l57568 78x 10° 22.00 4
28 -0.134919690832 86x 10° 24.00 4
29 -0.231 512250 53480xlO- 1 16.00 4
30 0.123631 25492901 x 10- 1 7 24.00 4
31 0.210583219729 40X 10- 2 8 8.00 4
32 -0.33958519026368xlO- 3 10 2.00 4
33 0.559936517715 92X 10- 2 4 28.00 5
34 -0.303351 180556 46X 1O-3 8 14.00 6
3
35 -0.213654886883 20X 10 2 1.00 25 325 1.16 1.00
36 0.266415691492 nx L0 5 0.00 25 300 1.19 1.00
37 -0.240272 [22045 57X 105 l.00 25 300 1.19 1.00
38 -0.283416034239 99x 103 3.00 IS 275 1.25 1.00
39 0.212472 844 001 79X 103 3.00 20 275 1.22 1.00
11; bi f3i A; Bi Ci Dj
40 -0.666422 765 4Q7 51 x 10° 3.500 0.875 0.300 0.700 0.3 10.0 275
41 0.726086323498 97x 10° 3.500 0,925 0.300 0.700 0.3 10.0 275
42 0.55068668612842xlO- 1 3.000 0.875 0.300 0.700 1.0 12.5 275
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
A NEW EQUATION OF STATE FOR CARBON DIOXIDE 1545
TABLE 32. The residual part of the dimensionless Helmholtz energy (f/ and its derivatives'
7 34 39 42
4l= ~ n;oJiri+ ~ n;oJirie-~i+ ~ n;;oJirie-ai(li-Ei)2_f3i(T-y)2 + ~ njAbiOW- with A= (P+B;[(8-1)2]a i
;= 1 ;=8 ;=35 ;=40
¢~= ±
i= I
n;d;oJ;-lri+ f
i=8
n;e-~'[oJi-lri(d;-c;OCi)]+ f n;oJ;rie-ai(li-<'i)2-f3i(T-Yi)2[~-2a;(8- EJ]
;=35 8
+ ~
i=40
n;[Abi('l'+ 8 a'l') + aAb;
a8 a8
ow-]
7 34 39
¢~li= ~ n;~;(d;-l)oJi-2r;+ ~ nje-~'[oJi-2r;«dj-cjOCi)(dj-l-cjOC;)-cJOCi)]+
;= 1 j=8
2: njrie-ai(li-<'i)2_f3i(T-Yi)2
;=35
~
+.LJ n·
;=40'
[ ab ,.( -
fh
2
'l')
a'l'+ 8a- - + 8
otJor
bi
aa - -
b
a'l'+ -aa ; ( 'l'+8-
otJ or or
a'l') + -
otJ
2 b
a a-; 8'l' 1
otJor
a'l'
a8 = -2C (8-1)'l'j
a2'l'
iAb; . a82 = [2C j ( S 1)2 1 ]2C j '¥
aT =2b abi - 1+4ffb;(b -l )ab;-2
j j a'l'
. a; = -2D;(,r-l)'l'
ia b; 2 aa a2'l'
-alXJr = - A.b. - Ab;-1(8-1 )[(8-1 )2]1/(2f3;)-1_2Bb (b -1 )Ab i - 2 _ -;;-:r =[2D,( T-l )2-1 ]?DiW
' , {3j ( ( a8
a2'l'
a8ar =4C;D;C8-1)(r-l)'l'
ment was achieved by refitting the b i' C i ' and D i of the 216 K:%;T~ 1100 K and 0 MPa:%;p:%; 800 MPa.
nonanalytic tenns.
The new fundamental equation for carbon dioxide,
Eq. (6.1), in combination with the fonnulations according to Estimations for the uncertainty of Eq. (6.1) are given in Sec.
Eqs. (6.3) and (6.5), is constrained to the critical parameters 8; the quality of the new equation of state is discussed in Sec.
given in Sec. 3.2. It is valid for the entire fluid region cov- 7 Th~ n~c.~"::lry rl~riv::ltives of Eq (6 S) are given in Table
ered by reliable data, namely for 32.
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
1546 R. SPAN AND W. WAGNER
0.02
~
<]
0 0
~
Tt
-0.02
00';
Q..
<],
0 0
0
......
T\
I •
-0.05
0.05
~ _:"}.-.-~-.-.~-.:.~.~?~~,
~
<J
0 0
0
Tl
I
-0.05
220 240 260 280 30 . 300 301 302 303 304 305
Temperature T J K Temperature T J K
58 58
o Duschek et al. o Duschek et al.
Aux. eqs. (see Sec. 3) Aux. eqs. (see Sec. 3)
15
-.-.- Ely et al. -. - . - Chen et al. 27
3
- .. - .. - Angus et at. , aux. eqs.
FIG. 19. Relative deviations 100Ay= 100 (Yexp-Ycaic)IYexp
FIG. 18. Relative deviations 100Ay= 100 (Yexp-Ycalc)IYexp(Y=Ps,P',p") (y = P s' p' ,p") of the near critical experimental saturation data of Duschek
of the experimental saturation data of Duschek at al. 58 from values calcu- et al. 58 from values calculated from Eq. (6.1). Values calculated from the
lated from Eq. (6.1). Values calculated from auxiliary equations presented in auxiliary equations presented in Sec. 3 and from the crossover of Chen
Sec. 3, the equation of state of Ely et al. 15 and the auxiliary equations of et al. 27 are plotted for comparison.
Angus et at. 3 are plotted for comparison.
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
A NEW EQUATION OF STATE FOR CARBON DIOXIDE 1547
-~
<:] -.- .--- .-'
-.-- j
sity can be considered to be within the experimental uncer-
tainty .
0
s::
0
- -~ - - - - - -- Cb@
The crossover equation of Chen et ai. 27 is able to repro-
duce the vapor pressure data within their experimental un-
Tc certainty but it cannot represent the saturated vapor and liq-
-2 300 310
270 280 290 uid densities within the uncertainty of these state-of-the-art
2 .j measurements. The auxiliary equations presented in Sec.3
\ !/ reproduce the experimental data slightly better than Eq. (6.1)
~~
<:]
0
-' \- .... 0
_Cb-ro /
.f
does. If one is interested in thermodynamically consistent
values for all properties on the phase boundary, the evalua-
§ 000 \. 00
tion of Eq. (6.1) in combination with Eq. (2.2) should be
Tt 0060 Tc
"- used and not the simpler auxiliary equations given in Sec. 3.
-2 300
220 240 260 280
2 7.1.2 Caloric Properties on the Coexistence Curve
\
~cJ
<:]
Figure 20 shows the representation of group 1 data pro-
0 0 viding information on the caloric behavior on the phase
9 buundary. Tht; spt;t;d uf suund dala uf lilt; salulalt;d vapul
state measured by Novikov and Trelin 113 are reproduced
-2~2~2-0~--2~40~L-2-6~0~L-2~8-0~--3~0~
within ±0.5% for T<301 K and to ±3% for T>301 K.
Temperature T / K Only the last point 30 mK below the critical temperature
1s exceeds this limit, but since the temperature scale of No-
o Novikov and Trelin 113
llS
- - - Elyetal.
3 vikov and Trelin is uncertain by more than 30 mK (see
o Pecceu and Van Dael -.-.- Angus et al.
<:) Magee and Ely116 Sec. 4.10), we did not rely on this data point. The analytic
equations of state of Angus et ai. 3 and Ely et ai. 15 are not
FIG. 20. Relative deviations 100Lly= 100 (Yexp - Yca\c)IYexp (y=w",w',c;) able to follow the decreasing speed of sound when approach-
of experimental caloric data at saturation from values calculated from ing the critical point.
Eq. (6.1). Data calculated from the wide-range equations of Ely et al. 15 and Equation (6.1) represents both the speed of sound and the
Angus et al. 3 are plotted for comparison.
specific isobaric heat capacity in the saturated liquid state to
within ::t 1 %; the greater deviations of the w' data of Pecceu
and Van Dael l15 at temperatures between 249 K and 261 K
presented in this work (see Sec.3), to values calculated from
are due to a systematic offset in the data set. The equation of
the equation of state of Ely et al. 15 by using the phase equi-
Ely et ai. 15 yields results which are at least similar to those
librium condition and to values calculated from auxiliary
calculated from Eq. (6.1), but the formulation of Angus
equations for the saturation properties established by Angus
et ai. 3 fails in describing caloric properties in the saturated
et al. 3 Angus et al. recommend the use of their auxiliary
liquid state (see also the representation of c v in the homoge-
equations for phase equilibrium calculations since they yield
neous liquid region,' Sec. 7.2.2and Fig. 31).
better results than the IUPAC equation of state.
Figure 18 shows that Eq. (6.1) represents the very accurate
The new equation of state was fitted to the c;
data of
Magee and Elyl16 up to temperatures of 295 K since the
vapor pressure data of Duschek et aI. 58 to within ± 0.006%
in pressure. Saturated liquid and vapor densities are repre-
conversion of c c;
(J" to seems to be uncertain for higher tem-
peratures (see Sec. 3.8). Above 295 K the deviations between
sented to within ± 0.01 % in density up to 303.6 K and the converted data and the new fundamental equation remain
297.0 K, respectively. Approaching the critical point the de-
viations in density increase, but these deviations are still
within about ± 2 %; keeping in mind that c;
is strongly di-
verging when approaching the critical point, this result un-
within the experimental uncertainty of the data; thus, all the
derlines the consistent description of the critical region by
data are represented to within their experimental uncertainty.
Eq.(6.1).
None of the equations hitherto known is able to reproduce
the data at least roughly within their uncertainty (the uncer-
tainty values are given in Tables 7-9). 7.2 Single-Phase Region
Figure 19 shows the representation of the saturation prop- 7.2.1 Thermal Properties in the Single-Phase Region
erties in the critical region where a larger deviation scale is
used. The very good representation of the vapor pressure For carbon dioxide, the region where ppT data in refer-
data is practically not affected when approaching the critical ence quality are available extends up to pressures of 13 MPa
temperature and the density deviations increase only slightly, at temperatures up to 360 K. Within this region, the data sets
except for the last point measured for the saturated liquid of Duschek et al.,58.154 Gilgen et al., 159 Nowak et aI., 160d and
density. Duschek et al. 58 do not give any estimation for the Guo et al. 157 describe the ppT surface with an uncertainty of
uncertainty of this point. but only 36 mK below the critical approximately ± 0.02% in density and in the extended criti-
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
1548 R. SPAN AND W. WAGNER
217 K - 220 K
0.051"'""""""'"-..;:::---------------------------.
il
/ 1
o~~~~--~~~~~~--~~~~~--------------~
1 -. D .. _.~- .. -O- .. 8·-
··-·· _.. - .. - .. - .. _.. _.. - .. - .. - .. - .. .::::-..,.:.:.... -.. -
-0.05
r
'-
260K
0.05 I
1
~ -0.05
I
0..
;.<
290 - 292 K
0.05
~ / /
.- .....
.-' -.-
'-" / /
0
0
,.....-I
D D
0
,/
\ ,/
,/
\
-0.05
299 - 301 K
0.05 ./
.. ---
-0.05
0 2 4 6 8 10 12 14
Pressure p / MPa
15O
0 Duschek et al. 154 t> Jaeschke et al. Phase boundary
155 J
0 Duschek et at. 511 'l Jaeschke et al. (refr. index) Angus et al.
e Guo et a1.
1S7
~
l5S
Jaeschke et ai. (Burnet) Ely et al.
1S
159 l60b
0 Gi1gen et ai. tsl K1imeck
FIG. 21. Relative density deviations of very accurate ppT data at subcritical temperatures from values calculated from Eq. (6.1). Values calculated from the
wide-range equations of Ely et at. 15 and Angus et at. 3 are plotted for comparison.
cal region of ~ 0.02% in pressure. In an extended region ~ 0.05o/c in density. None of the existing equations of state
reaching up to pressures of 30 MPa and temperatures of had access to these data sets because the data have been
523 K, the data of Brachthiiuser l60 and Klimeck et al. 160b published since 1990. Equation (6.1) is able to reproduce
describe the ppT surface with an uncertainty of ~0.02% to these data within their experimental uncertainty (for the nu-
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
A NEW EQUATION OF STATE FOR CARBON DIOXIDE 1549
304K
0.05
-0.05
305 K
0.05
o
%
0..
>< 0 ....... .' - - .. -0 . _- -0 - . \ .. ··--··0·.
a,)
~ : 00 '.
-.... '-. >- •
,;
\<l
.-....u
~
u -0.05
~
0.. 313 K
><
a,) 0.05
~
---
0
0
.,......(
.. .............. ~.--
/
-0.05
318 K - 320 K
0.05
.. - .. -."-'.
<l
o ~~~~~~~~~~~~~~~D~D~~____~____~
QlDD
<J
'. <l \
\
merical values of the uncertainty see Table 12). Figures 21 to in the gas and supercritical region. The equation of Ely
23 show the representation of some reference data on typical et al. 1S yields a suitable description of the gas region at low
isotherms in order to illustrate this statement. temperatures but it has problems for temperatures above
Figure 21 additionally contains data of Jaeschke l50 ,155 about 250 K. The equation of Angus et al. 3 is not able to
which also give a high-quality description of the ppT surface reproduce the state-of-the-art data in the gas region. In the
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
1550 R. SPAN AND W. WAGNER
337 K - 344 K
0.1
.' f
'W •/ I:
I :
'.
x x
X
-0.1
521 K- 540 K
0.1
_.' ......
o ~.,.. •• ,.., •• "*
o
~ o ~- .. - .. -"- .. -" .,..."
IZI
~.-.-.-.-.-.-.-.-.-'-.-.-.
f-
I-
I
'-
-0.1 I I I
681 K - 698 K
0.1
.. w··-··-
I-
I I I I l I
-0.1
5 6 7 8 9 10 20 30
Pressure pI MPa
103 l6011 3
In. Holste et al. W Fenghour et al. Angus et al.
15
X Hoinkis
145
lSl KIimeck 160b Ely et al.
160
o Gilgen et at. 159 IZl Brachthauser
154
o Duschek et ai.
FfG. 23. Relative density deviations of very accurate ppT data at supercritical temperatures from values calculated from Eq. (6.1). Values calculated from the
wide-range equations of Ely et af. 15 and Angus et af. 3 are plotted for comparison.
liquid region, none of the existing equations of state is able High-quality p p T data in the extended critical region are
to represent the reference data of Duschek et al. 154 and shown in Fig. 22. The data of Duschek et al. 58, 154 and Gilgen
8
Gilgen et al. 159 at least roughly to within their experimental et al. 159 are supplemented by the selected data of Wentort8
17
uncertainty. and by the data of Ely et al. which are consistent with the
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
A NEW EQUATION OF STATE FOR CARBON DIOXIDE 1551
.e
~:<
.
-
\
FIG. 24. Relative density deviations of selected ppT data from values calculated from Eq. (6.1). Values calculated from the wide-range equations of Ely
t:{ u/.I'i and AIlt;u:\ t:{ uf.' alC 1l1ultc;d fur cU1l11l<l1i:'UIl.
reference data by approximately ± 0.03% in pressure. Even resentation of the p p T data in the surrounding of the critical
in this region, Eq. (6.1) represents the reference data to isochore. However, at lower and higher densities, but clearly
within their experimental uncertainty. The deviations be- within the range of its validity, the crossover equation does
tween the data and values calculated from the equation of not reproduce the reference data within their experimental
Ely et al. 15 increase up to 0.1 % in pressure, which corre- uncertainty.
sponds to five times the uncertainty of the data, while the Figure 23 shows the representation of high-quality data at
equation of Pitzer and Schreiber l6 deviates by up to 0.2%. higher temperatures and pressures. Up to 523 K and 30 MPa,
The equation of Chen et al.,27 which is especially designed the p p T surface is defined by the data of Brachthauser l60 and
for the description of this region, yields a very suitable rep- Klimeck et al. 160b with an uncertainty of less than ± 0.05%
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
1552 R. SPAN AND W. WAGNER
0.3
o ..
773K 973K
0.3r-----------------------------~--~ 0.3 r-------------------------------...,.,
./ X X X ./
...... /
)...
X X ~-.-.-.-. /
1-'-'-'-'- )... ):-~
)...)..../
._)...
-0.3 L....-_....I....-----L...---L....l-J.....1-L....L.L.._----J_--'--....I....-.L.....L.....L.J...J.J
10 100 1 10 100
Pressure p / MPa Pressure p / MPa
l60a
'W Fenghour et al. )... Vukalovich and Altunin 129 Angus eta!.3
x Yukalovich et at. uo Ely et at.tt;
FIG. 25. Relative density deviations of selected ppT data at high temperatures from values calculated from Eq. (6.1). Values calculated from the wide-range
equations of Ely et al. 15 and Angus et al. 3 are plotted for comparison.
in density. For temperatures up to 698 K, the recent data of 200 MPa is given in Fig. 24. Generally, the new equation of
Fenghour et ai. 160a improve the uncertainty of the ppT sur- state describes the reliable data better than the older equa-
face to less than ± 0.1 % in density. tions. The comparison with the data sets of Duschek et ai. 154
A representative view on the complete group 1 set of the and Gilgen et ai. 1S9 shows that the use of adjusted data (see
ppT data up to temperatures of -173 K and pressures of Sec. 4.10) was reasonable.
~t·~:~~·;:·~·;:~~-~-,,,~.
o
o
,I
100 300 500 700 900
Pressure p / MPa Pressure p / MPa
132 3
[!] Michels et al. 123, a 6. Juza et al. -.-.- Angus et at.
15
o Tsiklis et at. 139 o Shmonov and Shmulovich 148 _.. - .. - Ely et al.
a Corrected data, see Sec. 4.10
FIG. 26. Relative density deviations of selected ppT data at high pressures from values calculated from Eq. (6.1), Values calculated from the wide-range
eqllationc; of Fly pf 11/ 15 ann Angm pf 111 3 ar~ pl()tt~n for ('ompari~on: in thic; prec;c;lIr~ ranee, thec;e two eqllations of state are at least partly extrapolated (see
Table I),
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
A NEW EQUATION OF STArE FOR CARBON DIOXIDE 1553
303 K
2
.,/~
~ •• - •• -
.. '" 0'1
-.: / •• ":..!. •• ......... _.
-'.J~ 0
- ........
, -..,
u
]
(,,)~
-2
Q.,
333 K 393 K
:l 2
~
0
2! ., . _1
~: ,~,,::".~'::",~~::~~
0
...... 0
-2
0.1 10 100
-: 0.1 1 10 100
Pressure p / MPa Pressure p / MPa
o Bender et al. 17'T, M <> Masi and Petkof 165 - .. - .. - Ely et al.
a Corrected data, see Sec. 4.10
FIG. 27. Relative deviations of selected isobaric heat capacity data from values calculated from Eq. (6.1). Values calculated from the wide-range equations of
Ely et at. 15 and Angus et at. 3 are plotted for comparison.
At temperatures above 698 K, Eq. (6.1) is essentially Sec. 4.10. The excellent representation of the data within the
based on the same data sets as the equations of Ely et al. 15 gas region, where the well-known contribution of c~ is pre-
and Angus et al. 3 Nevertheless, the three equations show sig- dominant, justifies this correction. In the region of the super-
nificantly different courses for pressures above 15 MPa. fig- critical maximum of c p' the deviations increase to about
ure 25 shows that the equation of Angus et al. 3 yields the ± 1% (+ 1.3% for a single data point) but the authors esti-
best representation of the data of Vukalovich and mate that the uncertainty of their data also increases to
Altunin. 128 ,129 However, in order to achieve a consistent de- ± 0.9% in this region.
scription of the caloric and thermal properties in other re- In Fig. 28, absolute values of the specific isobaric heat
gions, we had to assume that this data set exhibits a system- capacity in the low temperature gas region are plotted. When
atic error. The data of Fenghour et al. 160a support our approaching the sublimation curve, Eq. (6.1) yields a reason-
interpretation of the p p T surface for temperatures up to
698 K.
Figure 26 shows how the p p T data at very high pressures
are represented. Since the range of validity is limited to
100 MPa for the equation of Angus et al. 3 and to 300 MPa
for the equation of Ely et aI., 15 these equations are already
extrapolated when plotting values calculated from these 0.78
equations in Fig. 26. With an estimated uncertainty ranging ,-..,
--
OJ)
Shmonov and Shmulovich 148 represent the transition to the ~
'-"
extrapolation range of Eg. (6.1) which is discussed in Sec. ~
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
1554 R. SPAN AND W. WAGNER
o
_<>. _<>. - . - . - . -
~o
- . - . - . -<)
0.95
~ _.------
bl) ...-
...!:.d ...-
~
0.90 /'
...-
- --
...- ...-
'-'
./
cJ .,
0.85 \bc\)'f.~Itt'- - <> Magee and Ely116
\/
/ Eq. (6.1) 3
/ Angus et at.
500 550 600 650 700 r- Ely et ai.
IS
Temperature T / K 0.80
0 5 10 15 20 25 30
3
o Bender et al. 177, a Angus et at. Pressure plMPa
15
<> Schrock 166 Ely et al.
a Corrected data, see Sec. 4.10 FIG. 31. Representation of the isochoric heat capacity on high-density iso-
chores. For each of the isochores the plotted pressure range starts at the
FIG. 29. Relative deviations of isobaric heat capacity data at high tempera- corresponding vapor pressure. Values calculated from the wide-range equa-
tures from values calculated from Eq. (6.1). Values calculated from the tions of Ely et al. 15 and Angus et a1. 3 are plotted for comparison.
wide-range equarions of Ely er at. 15 and Angus et at. 3 are plOtted for com-
parison.
baric heat capacity on the sublimation line at temperatures
below T t . Due to the high triple-point pressure of CO 2
able extrapolation of the plotted isobars. The equation of (Pt = 0.51795 MPa) and the widespread use of dry ice, an
Angus et at. 3 also shows a reasonable extrapolation behav- unreasonable behavior in this region is less acceptable for
ior, but comparisons with recent data l77 at a temperature of carbon dioxide than for other substances.
233 K indicate that the resulting values are too high. In con- At ten1peratures above 400 K, the values calculated from
trast to this, the isobars calculated from the equation of Ely the different equations of state agree much better with each
et al. 15 intersect with each other and with the specific iso- other than with the data of Schrock l66 which are plotted in
3
265 kg/m
10 ! 000
Ii 00000
o
o
I; .. _.. -,,-,,_ .. _,,_,,_,,_
O~--------~I~'~~~'~~'~~~·~~
1/
-1 0 '----'----'-----'(~._ " - - - - - - - - - L _ - ' - - _ L - - _
280 300 320 340 280 300 320 340
Temperature T / K Temperature T / K
/:). Amirkhanov and Polikhronidi 188 o Magee and Ely116 Phase boundary
187 3
"V Amirkh::mnv pt nl Angus et a1.
I5
Ely et al.
FIG. 30. Relati\'e de\'iations of selected isochoric heat capacity data from values calculated from Eq. (6.1). Values calculated from the wide-range equations
of Ely el (//. I" and Angus ef al.·' are plotted for comparison.
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
A NEW EQUATION OF STATE FOR CARBON DIOXIDE 1555
308K 4~_________
29_8_K_-_3_2~3_K~______- ,
2~----------~~----~/~\-.----'
a ; \
a a a a CLO-O-fl..0 r!: '"' 'fa a
.. - .. - ..-.:.:.:::Y ::=,'t '-'r:'::"'.-. . .....
Y! "
~~ _2L---~--~--~--~---L--~~~--~
i '--
~: 373 K
io
o.......
4 6 8 10
Pressure p / MPa
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
1556 R. SPAN AND W. WAGNER
o
DODO
O~------------------------------~
.. - .. - .. - .. - .. _.. - .. - .. - .. - "
FIG.34. Relative deviations of experimentally determined Joule-Thomson coefficients from values calculated from Eq. (6.1). Values calculated from the
wide-range equations of Ely et al. 15 and Angus et al. 3 are plotted for comparison. The data of Vukalovich et al. 206 were not used when establishing Eq. (6.1).
450 MPa. These data are represented by Eq. (6.1) to within stances and general approaches, they will be discussed
their expected uncertainty of approximately::!:::: 2 %. Figure 33 elsewhere. 35 Here, the considerations are restricted to the
shows representative deviation plots covering the lowest and new fundamental equation for carbon dioxide; the following
highest temperatures of this data set. None of the previous subsections will give only a brief survey on the extrapolation
equations of state is able to yield a reasonable representation behavior of Eq. (6.0.
of these data. The equation of Angus et al. 3 is only valid up
to 100 MPa, and extrapolated values of the speed of sound
are expected to be uncertain. But the equation of Ely et aZ. 15
should yield reliable results at least within the range of its
validity, namely up to 300 MPa.
At pressures up to 1.5 MPa, the Joule Thomson measure
ments of Bender et al. 177 were used in the nonlinear fit but
these data could not be represented without systematic,
slightly temperature dependent deviations. After work on the
new equation of state had been completed, the isothermal
Joule-Thomson coefficients measured by Vukalovich ~ 200
et af. 206 were converted into differential J oule-Thomson co-
efficients by the use of the specific isobaric heat capacity ~
'-"
calculated from Eq. (6.1). The new equation represents the ~
converted data of Vukalovich et at. without systematic de- ..s o 1200 K
viations (see Fig. 34); the very accurate ppT data set in the Eo-. 150 o 1-100 K
~
gas region prevented the representation of the measurements • 1600 K
of Bender et al. which deviate from the correlation by ap- Eq. (6.1) 3
proximately 0.5% to 1%. Angus et at.
15
Ely et al.
7.3 Extrapolation Behavior of the New Fundamenta.1
Equation 100 "---1oo...l...-0-l.--2oo-L--- . L . . - -....1.- --'----400--'--0-
0 3oo0
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
A NEW EQUATION OF STATE FOR CARBON DIOXIDE 1557
-~
4)
I-c
::s
15000
I
I
o
CI)
CI)
~ 10000 I
0..
"0 1.5r
5000 / 0
C,)
~
.L_._. "0
!U
O~~~I-~~~-L-L~~~~~~ cz:::
1300 1600 1900 2200 2500 1.0r-
Density p / (kg/m3)
FIG. 36. Representation of experimental data describing the Hugoniot curve
of carbon dioxide. Hugoniot curves calculated from the wide-range equa- 0.5
tions of Ely et af. 15 and Angus et af. 3 are plotted for comparison.
0
5 10 20
7.3.1 Extrapolation Beyond the Range of Primary Data
Reduced temperature T / 1'c
The range of validity of Eq. (6.1) is based on the range FIG. 37. The so-called ideal curves calculated from Eq. (6.1) and plotted in
where reliable data of thermodynamical properties exist. a piPe' log(T1Tc) diagram. The louIe-Thomson inversion curve and the
louIe inversion curve exceed the temperature range in which Eq. (6.1) is
However, two kinds of data exceed this range.
fitted to experimental data.
Figure 35 shows the representation of fugacities reaching
up to approximately 3600 MPa at temperatures between
1200 K and 1600 K. The data of Haselton et at. 231 were not h - hOH= 0.5(p - POH)' (lIPOH+ lip), (7.1)
used when developing Eq. (6.1) since it is difficult to esti- where h is the enthalpy, p the pressure, and P the density
mate the uncertainty of data originating from measurements after releasing the shock wave and hOH , POH' and POH are the
of chemical equilibria and since the logarithmic structure of initial values. Even though it is not yet clear whether these
the dependency between the fugacity 'P and the reduced measurements describe eqUilibrium states at all, comparisons
Helmholtz energy (see Table 3) prevents an inclusion into with these data are the only source of experimental informa-
the linear optimization proct:llurt:. Nevertheless, Eq. (6.1) tion on (he extrapulaliun behavior of an equation of state at
follows the course of the measurements, whereas both the very high pressures: Figure 36 shows Hugoniot plots calcu-
equation of Angus et ai. 3 and the equation of Ely et al. IS lated from Eq. (6.1) and from the other two equations of state
yield fugacities which are significantly too large. At least at considered here, compared with the data of Nellis et al. 232
low pressures, the remaining systematic deviations cannot be and Schott. 233 At approximately 34000 MPa, Nellis et al.
explained by a misbehavior of the equation of state. All observed a kink in the course of the Hugoniot curve which is
equations of state which have been investigated result in very interpreted as an indication of a spontaneous disintegration
similar fugacities at pressures below 1000 MPa. In this re- reaction. So it can be seen that Eq. (6.1) yields a reasonable
gion, the experimental results are inconsistent with p p T data description of the Hugoniot curve of carbon dioxide up to the
at lower temperatures. limits of chemical stability. A Tp plot, which is not given
At even highp.r pressures, shock wave measurements result here, shows that the Hugoniot curves calculated from the
in data for the Hugoniot relation equations of Ely et ai. 15 and Angus et al. 3 run into low tem-
perature solutions corresponding to solid states above densi-
ties of about 1400 kg/m 3.
TABLE 33. The definition of the zeroth- and first-order ideal curves of the
compression factor 7.3.2 Representation of "Ideal Curves"
Name of the ideal curve Definition Various authors (see Refs. 234 to 238) have discussed
plots of so-called ideal curvese as a lml"'~r"::tl hehavlOf of
Classical ideal curve Z=pv/RT= 1
Boyle curve caZlap)T=o
louie-Thomson inversion curve (aZlaT),,=o e"Ideal curves" are curves which connect all states where a special property
louie inversion curve «(')Z/(Ip),,= 0 of a fluid is equal to the corresponding property of the hypothetic ideal gas
in the same state.
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
1558 R. SPAN AND W. WAGNER
G F
: . _ _ ._., F
E
! ~.-.-.-.-.-. r'-'-'-'-'-'-'-'-'-'-'
._.J D i
._._._._._ ....L _._.-,
E t;\l
Q..
._.L._.-,
c
'-'-l~'"
. ~/
j
D i'-._._. -
~
t::l.. 101
~
D
G
('-)
. E/'
/. /'./'
/'
/'
/././
FlU. 38. Tult:ram;e diagram fur densiLies calculared from Eq. (6.1). In region FlU. 40. Tolt:rance diagram fur isobaric heal capaciLies cai<.:Uiared [rum
B the uncertainty in pressure is given. Eq. (6.1). In the immediate vicinity of the critical point (region G) it is
difficult to estimate an uncertainty in c p because of the growing influence of
uncertainties in temperature and pressure measurement.
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
A NEW EQUATION OF STATE FOR CARBON DIOXIDE 1559
peratures beyond the limits of the chemical stability of car- the liquid-vapor phase boundary within their experimental
bon dioxide. uncertainty. The consideration of state-of-the-art data has re-
sulted in a previously unequalled accuracy for the regions of
8. Uncertainty of the New Equation of major technical interest. Intensive work on the consistency of
State the data set used has led to reasonable results in regions with
a poor data situation.
Estimates for the uncertainty of an empirical equation of Special interest has been focused on the behavior of ca-
state have to be guided by comparisons with experimental loric properties in the critical region and on the extrapolation
data. In regions where no data are available, comparisons behavior of empirical equations of state. The introduction of
with existing ~quations of state can be used as a substitute. A nonanalytic terms enables the new wide-range equation of
conservative estimation of the uncertainty of ppT, w, and state to represent the isochoric heat capacity and the speed of
c p values calculated from Eq. (6.1) is illustrated in the toler- sound even in the immediate vicinity of the critical point; up
ance diagrams Figs. 38 to 40. Uncertainties in c v correspond to now, this attribute has only been a domain of scaled equa-
to the uncertainties given for c p in the liquid and gas region. tions of state which introduce iterative dependencies between
In the extended critical region, the uncertainty in c v may different sets of variables and a limited range of validity. An
exceed the uncertainty in c p' The uncertainty of 6.h or 6.u examination of the extrapolation behavior of empirical equa-
values calculated from Eq. (6.1) is less than or equal to the tions of state has resulted in new approaches which have
uncertainty in c p or c v' respectively. been used for the new formulation. For the basic properties
Outside the range of its valilliLy, Ell. (6.1) shvulll yield of carbon dioxide such as pressure, enthalpy, and fugacity,
reasonable results for the basic thermodynamic properties the new fundamental equation should yield reasonable re-
like pressure, enthalpy, and fugacity within the whole chemi- sults within the whole region of chemical stability.
cally stable region of carbon dioxide. Of course, the extrapo-
lation results have an increased uncertainty which cannot be , O. Appendix: ThermOdynamic properties
estimated. The calculation of derived properties such as the of Carbon Dioxide
speed of sound or specific heat capacities is not recom-
mended beyond the limits of validity. If such data are In order to preserve thermodynamic consistency all values
needed, the results should be checked carefully. presented in Tables 34 and 35 were calculated only from the
new equation of state, Eq. (6.1). Ideally, each entry in the
table should be given to one more significant figure than the
9. Conclusions input data warrant, but a strict adherence to this principle is
difficult, and a possible conflict has always been avoided by
Based on a comprehensive study on the experimental data including more figures than are strictly necessary. Especially
for the thermodynamic properties of carbon dioxide, a new in the extended critical region, interpolation between values
fundamental equation in the form of the Helmholtz energy given may result in uncertainties which are significantly
has been developed. This empirical formulation is valid in larger then the uncertainties of Eq. (6.1). For sophisticated
the fluid region up to temperatures of 1100 K and at pres- applications, values should be calculated directly from Eq.
sures up to 800 MPa. The equation is able to represent al- (6.1); a computer-code suitable for such applications can be
most all the reliable data in the homogeneous region and on obtained from the authors. .
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
1560 R. SPAN AND W. WAGNER
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
A NEW EQUATION OF STATE FOR CARBON DIOXIDE 1561
'Triple point.
bCritical point.
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
1562 R. SPAN AND W. WAGNER
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
A NEW EQUATION OF STATE FOR CARBON DIOXIDE 1563
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
1564 R. SPAN AND W. WAGNER
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
A NEW EQUATION OF STATE FOR CARBON DIOXIDE 1565
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
1566 R. SPAN AND W. WAGNER
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
A NEW EQUATION OF STATE FOR CARBON DIOXIDE 1567
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
1568 R. SPAN AND W. WAGNER
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
A NEW EQUATION OF STATE FOR CARBON DIOXIDE 1569
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
1570 R. SPAN AND W. WAGNER
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
A NEW EQUATION OF STATE FOR CARBON DIOXIDE 1571
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
1572 R. SPAN AND W. WAGNER
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
A NEW EQUATION OF STATE FOR CARBON DIOXIDE 1573
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
1574 R. SPAN AND W. WAGNER
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
A NEW EaUATION OF STATE FOR CARBON DIOXIDE 1575
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
1576 R. SPAN AND W. WAGNER
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
A NEW EQUATION OF STATE FOR CARBON DIOXIDE 1577
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
1578 R. SPAN AND W. WAGNER
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
A NEW EQUATION OF STATE FOR CARBON DIOXIDE 1579
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
1580 R. SPAN AND W. WAGNER
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
A NEW EQUATION OF STATE FOR CARBON DIOXIDE 1581
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
1582 R. SPAN AND W. WAGNER
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
A NEW EQUArlON OF STATE FOR CARBON DIOXIDE:. 1583
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
1584 R. SPAN AND W. WAGNER
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
A NEW EQUATION OF STATE FOR CARBON DIOXIDE 1585
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
1586 R. SPAN AND W. WAGNER
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
A NEW EQUATION OF STATE FOR CARBON DIOXIDE 1587
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
1588 R. SPAN AND W. WAGNER
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
A NEW EQUATION OF STATE FOR CARBON DIOXIDE 1589
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
1590 R. SPAN AND W. WAGNER
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
A NEW EQUATION OF STATE FOR CARBON DIOXIDE 1591
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
1592 R. SPAN ANDW. WAGNER
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
A NEW EQUATION OF STATE FOR CARBON DIOXIDE 1593
"Sublimation temperature.
bMelting temperature.
'Saturation temperature.
11. Acknowledgments 10K. E. Starling, P. N. Batdorf, and Y. C. Kwok, Hydrogen Proc. 51, 2, 86
(l '172).
1I R. Meyer-Pittroff, Dissertation, TC Miinchen (1973).
The authors are grateful to the Deutsche Forschungsge-
12F. H. Huang, M. H. Li, L. L. Lee, K. E. Starling, and F. T. H. Chung, J.
meinschaft for their financial support of this project. Chern. Eng. Jpn. 18. 490 (1985).
13 J. Ewers and W. Wagner, A Methodjor Uptimizing the Structure ofEqua-
) A. Saul and W. Wagner. 1. Phys. Chern. Ref. Data 18, 1537 (1989). 171. F. Ely. W. M. Haynes. and B. C. Bain, 1. Chern. Thermodyn. 21, 879
flU. Setzrnann and W. Wagner. lnt. J. Thermophys. 10.1103 (1989). (1989).
7 E. Bender. "Equatiolls of State ExactlY Represellling the Phase Behavior
18R. T. Jacobsen and R. B. Stewart, J. Phys. Chern. Ref. Data 2,757 (1972).
of Pure Substances." in Proceedings of the Fifth Symposium 011 711ermo- 18a K. S. Pitzer and S. M. Sterner. 1. Chern. Phys. 101, 4, 311 J (1994).
phnics Properfles. C. F. Bonila. coordinating editor (AS ME. New York, 19p. Schofield. Phys. Rev. Lett. 22. 606 (1969).
1970). 20G. A. Chapela and J. S. Rowlinson. 1. Chem. Soc .. Faraday Trans. (I) 70,
Sy. Y. Altunin and O. G. Gadetskii. Teploenergetica 18, 81 (1971). 584 (1974).
9W. A. Stein. Chern. Eng. Sci. 27. 1371 (1972). 2: J. F. Ely. 1. W. Magee. and W. M. Haynes, Tlzennophysical Properties of
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
1594 R. SPAN AND W. WAGNER
::,pecial High CO 2 Content Mixtures, Research Report RR-110 (National 60 A Michels, B. Blaisse, and J. Hoogschagen, Physica 9, 565 (1942).
Bureau of Standards, Boulder, 1987). 61K Clusius, U. Piesbergen, and E. Varde, Helv. Chim. Acta 43, 1290
22M. Vicenti-Missoni, J. M. H. Levell Sengers, and M. S. Green, J. Res. (1960).
1\'at1. Bur. Stand. A 73, 563 (1969). 62H. Bois and A. P. Wills, VerhandL der Deutschen Phys. Ges. 1, 168
23T. A. Murphy, J. V. Sengers, and J. M. H. Levelt Sengers, Analysis of the (l899).
Pressure of Gases Ncar the Critical Point ill Terms of a Scaled Equation 001. P. Kuenen and W. G. Robson, Philos. Mag. 6, 149 (1902).
of State, in Proc. 6th Symp. Therrnophys. Prop., P. E. Liley. coordinating MK. H. annes and S Weber, Com. Phys. Lab. Univ. Leidcn 137b. 8
editor (AS ME, New York, 1973). (1913).
24p. C. Albright, T.J. Edwards, Z. Y. Chen, and J. V. Sengers, J. Chern. 65H. Siemens, Ann. Phys. 42, 871 (1913).
Phys 87, 1717 (1987). 66S. Weber and H. K annes, Com. Phys. Lab. Univ. Leiden 137c, 26
25p. C. Albright, Z. Y. Chen, and J. V. Sengers, Phys. Rev., Rapid Com- (1913).
mun. 36. 877 (1987). 67F. Henning, Ann. Phys. 43, 282 (1914).
26
0. D. Erickson, T. W. Leland, and J. F. Ely, Fluid Phase Equilibria 37, 68W. Heuse and J. Otto, Ann. Phys. 9, 486 (1931).
185 (1987). wW. Heuse and 1. Otto, Ann. Phys. 14, 185 (1932).
27Z. Y. Chen, A Abbaci, S. Tang, and J. V. Sengers, Phys. Rev. A 42, 70W. F. Giauque and C. J. Egan, J. Chern. Phys. 5, 45 (1936).
4470 (1990). 71 A. W. Tickner and F. P. Losing, J. Phys. Colloid Chern. 55, 733 (1951).
28 S. n. Kiselev, 1. G. Kostyukova, and A. A Povodyrev, Int. J. Thermo- 72D. Ambrose, Trans. Faraday Soc. 52, 772 (1955).
phys. 12, 877 (1991). 73M. J. Hiza, Cryogenics 10, 106 (1970).
29R. Fox, Fluid Phase Equilibria 14,45 (1983). 74Z. Bilkadi, M. W. Lee, and J. Biegeleisen, J. Chern. Phys. 62, 2087
30U. Setzmann and W. Wagner, J. Phys. Chem. Ref. Data 20,1061 (1990). (1974).
31 J. Ahrendts and H. D. Baehr, Forseh. lng. Wes. 45, 1, 1 (1979). 75c. E. Bryson, V. Cazcarra, and L L. Levenson, J. Chern. Eng. Data 19,
l2W. Wagner, Cryogenics 12, 214 (1972). 107 (1974).
33J. Ahrendts and H. D. Baehr, Forseh. lng. Wes. 45, 2, 51 (1979). 76B. Fernandez-Fassnacht and F. del Rio, J. Chern. Thermodyn. 16, 469
34W. Wagner, VDI Fortschritt-Ber., Reihe 3, Nr. 39 (VOl-Verlag, Dussel- (1984).
dorf, 1974). 77W. H. Keesom, Com. Phys. Lab. Univ. Leiden 88, I (1903).
35 R. Span and W. Wagner, Int. 1. Therrnophys. (submitted). 78C. F. Jenkin and D. R. Pye, Philos. Trans. R. Soc. 213, 67 (1914).
36H. Preston-Thomas, \1etrologia 27,3 (1990). 79 0. C. Bridgeman, J. Am. Chern. Soc. 49,1174 (1927).
37H. Preston-Thomas, P. Bloembergen, and T. J. Quinn, Supplementary 80 A. Michels, B. Blaisse, and C. Michels, Proc. R. Soc. A160, 358 (1936).
Information for the /Illernational Temperature Scale of 1990 (Bureau 81 J. R. Roebuck, T. A Murell, and E. E. Miller, J. Am. Chern. Soc. 64, 400
lnternationale des Poids et Mesures, Sevres, 1990). (1942).
38R. L. Rusby, J. Chern. Themlodyn. 23, 1153 (1991). 82 A. Michels, T. Wassenar. Th. Zwietering, and P. Smits, Physica 16,501
39R. B. Bedford and C. G. M. Kirby, Metrologia 5, 83 (1969). (1950).
4DC. H. Meyers and M. S. Van Dusen, Bur. Stand. J. Res. 10, 381 (1933). S3H. H. Reamer, B. H. Sage, and W. N. Lacey, Ind. Eng. Chern 43, 2515
41 D. Ambrose, Brit. J. AppL Phys. 8, 32 (1957). (1951).
42D. R. Lovejoy, Nature 197, 353 (1963). 84J. A Bierlein and B. K. Webster, Ind. Eng. Chern. 45, 618 (1953).
43E. Haro, E. Fernandez-Fassnacht, and F. del Rio, Rev. Soc. Quim. Mex. 85 D. Cook, Proc. R. Soc. 219, 245 (1953).
23, 369 (1979). 86 D. Cook, Trans. Faraday Soc. 49, 716 (1953).
44L A. K. Staveley, L. Q. Lobo, and J. C. G. Calado, Cryogenics 21, 131 B7E. Schmidt and W. Thomas, Forseh. lng. Wes. 20, 161 (1954).
(1981). 8SR. H. Wentorf, J. Chern. Phys. 24, 607 (1956).
45R. Blanes-Rex, E. P. A Fernandez, and F. Guzman, Cryogenics 22, 113 89 A B. KIetskii, lng. Phys. Zeitschr. 7, 40 (1964).
(1982). 90R. G. P. Greig and R. S. Dadson, Brit. J. Appl. Phys. 17, 1633 (1966).
46F. Pavese and D. Ferri, TMCSI 5, 217 (1982). 91 J. L Edwards and D. P. Johnson, J. Res. Nat!. Bur. Stand. Sec. C 72,27
47R. E. Bedford, G. Bonnier, H. Maas, and F. Pavese, Metrologia 20, 145 (1968).
(1984). 92M. P. Vukalovich, V. P. Kobelev, and N. 1. Timoshenko, Teploenergetica
48G. Bonnier, Y. Hermier, and W. B. Qin, "Triple Point of Carbon Dioxide 15, 6, 80 (1968).
in a Multicompartment Sealed Cell," in Proceedillgs of the 2nd IMEKO 90y' A. Kirillin, S. A. Ulybin, and E. P. Zherdev, Teploenergetica 16, 6, 92
Symposium on Temperature Measuremelll in Industry and Science, B. (1969).
Suhl, coordinating editor, GDR 1984. 94E. P. Kholodov, N. I. Timoshenko, and A L. Yaminov, Teploenergetica
49W. Duschek, VDI Fortschritt-Ber., Reihe 3, Nr. 187 (VOl-Verlag, Dus- 19,4, 84 (1972).
seldorf, 1989). 95 A. Fredenslund and J. Mollerup, J. Chern. Soc. Farraday Trans. (1) 70,
50 J. M. H. Levelt Sengers and W. T. Chen, J. Chern. Phys. 56, 595 (1972). 1653 (1974).
51 M. R. Moldover, J. Chern. Phys. 61, 1766 (1974). 96R. J. Gugoni, J. W. Eldridge, V. C. Okay, and T. 1. Lee, A1.Ch.E. J. 2(1,
52 J. A. Lira, c. Edwards, and M. J. Bnckingham, Phys. Rev. A 15, 778 93 (1974).
(1977). 97 G. J. Besserer. and D. B. Robinson, J. Chern. Eng. Data 21, 81 (1975).
53M. R. Moldover, J. V. Sengers, R. W. Gammon, and R. J. Hocken, Rev. 981. Davalos, W. R. Anderson, R. B. Phelps, and A J. Kidnay, J. Chem.
Mod. Phys. 51, 79 (1979). Eng. Data 21,81 (1976).
54H. J. Baade, Dissenation, TV Brannschweig (1983). 99K. Stead and J. M. Williams, J. Chem. Thermodyn. 12, 265 (1980).
SST. J. Edwards, Ph.D. Thesis, University of Western Australia (1984). 1ooH. Kwang-Bae, N. Nagahama, and M. Hirata, J. Chern. Eng. Data 27, 25
56 J. V. Sengers and J. M. H. Level! Sengers, Annu. Rev. Phys. Chern. 37, (1982).
189 (1986) 100T. A Al-Sahhaf, A. J. Kidnay, and E. D. Sloan, Ind. Eng. Chem. Fund.
57Z. Y. Chen, P. C. Albright, and J. V. Sengers, Phys. Rev. A 41, 3161 22, 372 (1983).
(1990). 102H. G. Kratz, Dissertation, Ruhr-Universitat Bochum (1984).
5SW. Duschek, R. Kleinrahm, and W. Wagner. J. Chern. Thermodyn. 22, I03J. C. Holste, K. R. Hall, P. T. Eubank, G. Esper, M. Q. Warowny, D. M.
841 (1990). Bailey, J. G. Young, and M. T. Bellamy, J. Chern. Thermodyn. 19, 1233
59!. M. AbduJagatov, 1. M. Abdurakhamov, R. G. Batyrova, and N. G. (1987).
Polikhronidi, "Isochoric Heat Capacity and Liquid-Gas Coexistence 104T. S. Brown, A. J. Kidnay, and E. D. Sloan, Fluid Phase Equilibria 40,
Curve of Carbon Dioxide in the Region of Critical State," paper pre- 169 (1988).
sented at 11th Symp. on Therrnophys. Prop., Boulder, Colorado, 1991. 105 N. N. Shah, 1. A Zollweg, and W. B. Streett, J. Chern. Eng. Data 36, 188
59ar. M. Abdulagatov, N. G. PoIikhronidi, and R. G. Batyrova, Ber. Bun- (1991).
senges. Phys. Chem. 98, 1068 (I 994). lOS'L. Yurttas, J. C. Holste, K. R. Hall, B. E. Gammon, and K. N. Marsh, 1.
59br. M. Abdulagatov, N. G. Polikhronidi, and R. G. Batyrova, J. Chenl. Chern. Eng. Data 39, 418 (1994).
Thermodyn. 26, 1031 (1994). 106M. E.-H. Amagat, C. R. Acad. Sci. 114, 1093 (1892).
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
A NEW EQUATION OF STATE FOR CARBON DIOXIDE 1595
107U. Behn, Ann. Phys. 3, 733 (1900). 153p. J. McElroy, R. Battino, and M. K. Dowd, J. Chem. Therrnodyn. 21,
108H. H. Lowry and W. R. Erickson, J. Am. Chem. Soc. 49, 2729 (1927). 1287 (1989).
109 J. Straub, Wiirme und Stoffubertragung 5, 56 (1972). 154W. Duschek, R. Kleinrahm, and W. Wagner, J. Chem. Therrnodyn. 22,
1I0W. M. Haynes, Orthobaric Liquid Densities and Dielectric Constants of 827 (1990).
Carbon Dioxide, paper presented at the Cryog. Eng. Conf., Mass. Inst. of. 155M. Jaeschke, A. E. Humphreys, P. Van Caneghem, M. Fauveau, R.
Tech., Cambridge, 1985. Janssen-Yan Rosmalen, and Q. Pellei, GERG Technical Monograph
III G. J. Esper, Dissertation, Ruhr-Universitiit Bochum (1987). 4-The GERG Databank of High Accuracy Compressibility Factor Mea-
112 A. Eucken and F. Hauck, Z. Phys. Chem. 134, 161 (1928). surements, edited by Group Europeen de Recherches Gazieres (YDI-
113 1. 1. Novikov and Y. S. Trelin, Zur. Prik. Mech. Tech. Fiz. 2, 112 (1960). Verlag, Dusseldorf, 1990).
114K. I. Amirkhanov, N. G. Polikhronidi, B. G. Alibekov, and R. G. Baty- 156p. Nebendahl, YDI Fortschritt-Ber., Reihe 3, Nr. 212 (YDI-Yerlag, Dus-
rova, Teploenergetica 19, 1,61 (1972). seldorf, 1990).
IISW. Pecceu and W. Van Dae1, Physica 63, 154 (1972). 157X. Y. Guo, R. Kleinrahm, and W. Wagner, Experimentelle Untersuchung
1161. W. Magee and J. F. Ely, Int. J. Thermophys. 7, 1163 (1986). der systematischen Me[3fehler von Betriebsdichteaufnehmernfiir Erdgas-
117 e. F. Jenkin, Proc. R. Soc. A 98, 170 (1920). Me[3trecken - Tei! 1: Me[3ergebnisse for Stickstoff, Kohlendioxid. Argon.
1180. Maass and J. H. Mennie, Proc. R. Soc. A 110, 198 (1926). Neon. Ethan und Ethen, Report Lehrstuhl fur Thermodynamik, Ruhr-
1l9D. LeB. Cooper and O. Maass, Can. J. Res. 2, 388 (1930). Universitat Bochum, Bochum, 1992.
120D. LeB. Cooper and O. Maass, Can. J. Res. 4, 283 (1931). 158L. A. Weber, Int. J. Thermophys. 13, 1011 (1992).
121 W. Cawood and H. S. Patterson, J. Chem. Soc., 619 (1933). 159R. Gilgen, R. Kleinrahm, and W. Wagner, 1. Chem. Therrnodyn. 24, 1493
122 A. Michels and e. Michels, Proc. R. Soc. A 153, 201 llY~5). (1992).
123 A. Michels, e. Michels, and H. Wouters, Proc. R. Soc. A 153, 214 160K. Brachthauser, R. Kleinrahm, H.-W. Losch, and W. Wagner, YDI
(1935). Fortschritt-Ber., Reihe 8, Nr. 371 (YDI-Yerlag Dusseldorf, 1993).
124H. H. Reamer, R. H. Olds, B. H. Sage, and W. N. Lacey, Ind. Eng. Chem. 160a A. Fenghour, W. A. Wakeham, and J. T. R. Watson, J. Chem. Thermo-
36, 88 (1944). rlyn 27, ?lQ (lQQ'i)
l2SG. A. Bottomley, D. S. Massie, and R. Whitlaw-Gray, Proc. R. Soc. A 160b J. Klimeck, R. Kleinrahm, and W. Wagner (private communication)
200,201 (1950). (1995).
126T. Batuecas and e.G. Losa, An. Fis. Quim. B 50, 845 (1954). 160cZ. Gokmenoglu, Y. Xiong, and E. Kiran, 1. Chem. Eng. Data 41, 354
I27G. e. Kennedy. Am. 1. Scie. 252. 225 (954). (1996).
128M. P. Yukalovich and Y. Y. Altunin, Teploenergetica 6, 11,58 (1959). IOOtlp. Nowak, T. Tielkes, R. Kleinrahm, and W. Wagner (unpublished).
129M. P. Yukalovich and Y. Y. Altunin, Teploenergetica 9,5, 56 (1962). 160ew. Wagner, K. Brachthiiuser, R. KIeinrahm, and H.-W. Losch, Int. J.
130M. P. Yukalovich, Y. Y. Altunin, and N. 1. Timoshenko, Teploenergetica Therrnophys. 16, 399 (1995).
161 R. Kleinrahm and W. Wagner, YDI Fortschritt-Ber., Reihe 3, Nr. 92
10, 2, 92 (1963).
(YDI-Yerlag, Dusseldorf, 1986).
131 M. P. Yukalovich, Y. Y. Altunin, and N. 1. Timoshenko, Teploenergetica
162F. G. Keyes and S. e. Collins, Proc. N. A. S. 18,328 (1932).
10, 1, 85 (1963).
163G. B. Kistiakowsky and W. W. Rice, J. Chem. Phys. 7, 5, 281 (1939).
132 J. Juza, Y. Kmonicek, and O. Sifner, Physica 31, 1735 (1965).
I64S. R. de Groot and A. Michels, Physica 14, 218 (1948).
I33p. S. Ku and B. F. Dodge, J. Chem. Eng. Data 12, 158 (1967).
165J. F. Masi and B. Petkof, J. Res. Natl. Bur. Stand. 48, 179 (1952).
134 A. Sass, B. F. Dodge, and R. H. Bretton, J. Chem. Eng. Data 12, 168
166y. E. Schrock, Calorimetric Determination of Constant-Pressure Specific
(1967).
Heats of Carbon Dioxide at Elevated Pressures and Temperatures, Na-
135M. P. Yukalovich, Y. P. Kobelev, and N. 1. Timoshenko, Teploenergetica
tional Advisory Committee for Aeronautics, Tech. Note 2838 (1952).
15, 4, 81 (1968).
167L. B. Koppel and J. M. Smith, J. Chem. Eng. Data 5,437 (1960).
136E. A. Golovskii and Y. A. Tsymamyi, Teploenergetica 16, 1,67 (1969).
168M. P. Yukalovich, Y. Y. Altunin, and A. N. Gureev, Teploenergetica 11,
my. A. Kirillin. S. A. Ulybin, and E. P. Zherdev, Teploenergetica 16, 2, 94
9, 68 (1964).
(I969).
169M. P. Yukalovich and A. N. Gureev, Teploenergetica 11, 8, 80 (1964).
l38y. A. Kirillin. S. A. Ulybin, and E. P. Zherdev, "Experimental Investi-
170M. P. Yukalovich, Y. Y. Altunin, and A. N. Gureev, Teploenergetica 12,
gation of Carbon Dioxide Density at Temperatures from -50 to +200
7,58 (1965).
°C and Pressures up to 500 bar," in Proceedings of the First Interna- 171 S. L. Rivkin and Y. M. Gukov, Teploenergetica 15, 10, 72 (1968).
tional Conference on Calorimetry and Thermodynamics, Warsaw, 1969. l72y. Y. Altunin and D~D. Kuznetsov, Teploenergetica 16, 8, 82 (1969).
139D. S. Tsiklis, L. R. Linshits, and S. S. Tsimmermann, "Measurement and l73y. Y. Altunin and D. O. Kuznetsov, Teploenergetica 17, 11,91 (1970).
Calculation of the Molar Volume and Thermodynamic Properties of Car- I7-lS. L. Rivkin and Y. M. Gukov, Teploenergetica 18, 10,82 (1971).
bon DIOxide at High Pressures and Temperatures," in f'roceedlllgs Of the 175y' Y. Altunin and D. O. Kuznetsov, Teploenergetica 19, 6, 67 (1972).
1st International Conference on Calorimetry and Thermodynamics, War- I76S. Saegusa, T. Kobayashi, Y. Takaishi, and K. Watanabe, Bull. JSME 23,
saw. 1969. 2055 (1980).
l.lOy. A. Kirillin. S. A. Ulybin. and E. P. Zherdev, Teploenergetica 17, 5, 69 I77R. Bender, K. Bier, and G. Maurer, Ber. Bunsenges. Phys. Chern. 85. 778
(1970). (19151).
1,\1 Y. N. Popov and M. K. Sayapov. Teploenergetica 17,4.76 (1970). 178G. Ernst and U. E. Hochberg, 1. Chem. Thermodyn. 21. 407 (1989).
I·e M. P. Yukalovich. N. 1. Timoshenko. and Y. P. Kobelev, Teploenergetica 179G. Ernst, G. Maurer. and E. Wiederuh. 1. Chern. Thermodyn. 21. 53
17. 12.59 (1970). (1989).
1-1< W. Schonmann. Dissertation. TH Karlsruhe (971). 179aT JJmrbin. T _Y rm::1m. ~nrl J -P E. Colier. Rev. Sci In~trl1m. 65. 10,
I-I-IE. P. Kholodov. l\. l. Timoshenko. and A. L. Yamnov. Teploenergetica 3263 (1994).
19.3.84 (1972). 180H. W. Wooley. 1. Natl. Bur. Stand. 52. 289 (1954).
1-1 5 1. Hoinkis. Dissertation. TH Karlsruhe (1989). 181 H. D. Baehr, H. Hartmann, H.-e. Poh!. and H. Schomacker, Thermody-
1-I6 G.1. Besserer and D. B. Robinson. 1. Chem. Eng. Data 18. 137 (1973). l1amisc/ze FlIllktionen idealer Gase fiir Temperatllren bis 6000 oK, edited
1-17 D. S. Rasskazov. E. K. Petrov. G. A. Spiridonov. and E. R. Ushmaikin. by H. D. Baehr (Springer-Verlag, Berlin, 1968).
Teploenergetica 21. l. 80 (1974). 182L. V. Gurvich. Thermodynamic Properties of Characteristic Substances
l-l~ V. \1. Shmono\" and K. I. ShmulO\ich. Akad. Wouk SSSR 217. 935 (Nauka. Moskwa. 1979), Vol. II, Part II.
( 19741. 183 J. Chao. in TRC ThermodYllamic Tables Hydrocarbons. K. W. Marsh,
1-l'iW._W. R. Lau. Ph.D. Thesis. Texas A&~I University (1986). coordinating editor (Texas A&M University, Austin. 1983). Vol. YIII.
I~O \1. laeschke. ppT Data from Refractive Index Measurements (private 18-11. Chao (private communication) (1986).
communication) ( 19871. 185R. E. Pennington and K. A. Kobe. 1. Chern. Phys. 22, 1442(954)
I~I ~t Jaeschke. ppT Data from Burnett ~1easurements (private communica- 186 A. Michels and 1. Strijland. Physica 18, 613 (1952).
tion 1 ([ 9871. 187 K. I. Amirkhanov. N. G. Polikhronidi, and R. G. Batyrova, Teploener-
1~2 1. W. Magee and 1. F. Ely. rnt. 1. Thermophys. 9. 547 (J 9881. getica 17. 3. 70 (1970).
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions
1596 R. SPAN AND W. WAGNER
188K. I. Amirkhanov and N. G. Polikhronidi, Teploenergetica 18, 12, 59 2l3T. L. Cottrell, R. A. Hamilton, and R. P. Taubinger, Trans. Faraday Soc.
(1971). 52, 1310 (1956).
189 J. V. Sengers (private communication) (1986). 214D. Cook, Can. J. Chem. 35, 268 (1957).
19OJ. A. Lipa, C. Edwards, and M. J. Buckingham, Phys. Rev. Lett. 25, 1086 215 A. P. Masia and M. D. Pena, An. Fis. Quim. B 54, 661 (1958).
(1970). 216E. G. Butcher and R. S. Dadson, Proc. R. Soc. 227, 448 (1963).
191C. M. Herget, 1. Chern. Phys. 8, 537 (1940). 217 J. A. Huff and T. M. Reed, J. Chern. Eng. Data 8, 306 (1963).
1921. 1. Novikov and Y. S. Trelin, Tepioenergetica 9,2, 79 (1962). 218M. P. Vukalovich and Y. F. Masalov, Teploenergetica 13, 5, 58 (1966).
193L. L. Pitaevskaya and A. V. Bilevich, Russ. J. Phys. Chem. 47, 126 2l9R. S. Dadson, E. J. Evans, and 1. H. King, Proc. Phys. Soc. 92, 1115
(1973). (1967).
194W. Lemming, VDI Fortschritt-Ber., Reihe 19, Nr. 32 (VDI-Verlag, Dus- 220N. I. Timoshenko, V. P. Kobelev, and E. P. Kholodov, Teploenergetica
seldorf, 1989). 17,9,64(1970).
195 J. Noury, Etude de la vitesse des ultrasons dans les flu ides au voisinage 221 M. Waxman, H. Davis, and 1. R. Hastings, "A New Determination of the
du point critique (Academy of Sciences, 1951). Second Virial Coefficient of Carbon Dioxide at Temperatures Between
196y' S. Trelin and E. P. Sheludiakov, J.E.T.P. Lett. 3, 2, 63 (1966). 0° and 150°C, and an Evaluation of Its Reliability," in Proceedings of
\97R. W. Gammon, H. L. Swinney, and H. Z. Cummins, Phys. Rev. Lett. 19, the Sixth Symposium on Thermophysical Properties (ASME, New York,
1467 (1967). 1973).
198G. T. Feke, K. Fritsch, and E. F. Carome, Phys. Rev. Lett. 23, 1282 222K. Ohgaki, N. Sakai, Y. Kano, and T. Katayama, 1. Chern. Eng. Jpn. 17,
(1969) 545 (1984).
199 K. Fritsch and E. F. Carome, Behavior of Fluids in the Vicinity of the 223B. V. MaUu, G. Natarajan, and D. S. Viswanath, J. Chern. Thermodyn.
Critical Point, NASA Research Report CR-1670 (1970). 19,549 (1987).
2000. Maass and W. H. Barnes, Proc. R. Soc. A 111, 224 (1926). 224B. V. MaUu, G. Natarajan, and D. S. Viswanath, 1. Chern. Thermodyn.
201 M. P. Vukalovich and Y. F. Masalov, Teploenergetica 11, 7, 78 (1964). 21, 989 (1989).
202M. P. Vukalovich and Y. F. Masalov, Teploenergetica 11, 11, 75 (1964). 225R. Span, VDI Fortschritt-Ber., Reihe 6, Nr. 285 (VDI-Verlag, Dusseldorf,
203D. Moller, B. E. Gammon, K. N. Marsh, K. R. Hall, and J. C. Holste, J. 1993).
Chern. Thermodyn. 25, 1273 (1993). 226 J. M. H. Levelt Sengers, W. L. Greer, and J. V. Sengers, J. Phys. Chern.
204K. KrUger, VDl Fortschritt-Ber., Reihe 6, Nr. 1 (VDI-Verlag, Dusseldorf, Ref. Data 5, 1 (1976).
1964). 227 J. E. Verschaffelt, Com. Phys. Lab. Leiden 28 (1896).
205 H. D. Baehr, "Measurements of the Internal Energy and the Influence of 228K. G. Wilson, Phys. Rev. 4, 3174 (1974).
Gravity Near the Critical Point of CO 2, " in Proceedings of the Fourth 229p. Hill, J. Phys. Chem. Ref. Data 19, 1233 (1990).
Symposium on Thermophysical Properties (ASME, New York, 1968). 230W. Wagner, N. Kurzeja, and B. Pieperbeck, Fluid Phase Equilibria 79,
206M. P. Vukalovich, V. V. Altunin, K. Bulle, D. S. Rasskazov, and D. 151 (1992).
Ertel, Teploenergetica 17, 5, 60 (1970). 231 H. T. Haselton, W. E. Sharp, and R. C. Newton, Geophys. Res. Lett. 5,
207E. S. Burnett, Phys. Rev. 22, 590 (1923). 753 (1978).
208V. V. Altunin and A. N. Gureev, Teploenergetica 19, 1, S9 (1972). 232W. J. Nellis, A. C. Mitchell, F. H. Ree, N. C. Holmes, and R. J. Trainor,
208aL. Cusco, S. E. McBain, and G. Saville, 1. Chern. Thermodyn. 27, 721 J. Chern. Phys. 95, 5268 (1991).
(1995). 233 G. L. Schott, High Pres. Res. 6, 187 (1991).
209K. Schafer, Z. Phys. Chern. B 36, 85 (1937). 234E. H. Brown, Bu!. lnt. Inst. Refrig. 1, 169 (1960).
210K. E. MacCormack and W. G. Schneider, J. Chem. Phys. 18, 1269 21ST. E. Morsy, Dissertation, TH Karlsruhe (1963).
(1950). 236D. Straub, Dissertation, TH Karlsruhe (1964).
211 T. L. Cottrell and R. A. Hamilton, Trans. Faraday Soc. 51, 156 (1955). 237 A. Schaber, Dissertation. TH Karlsruhe (1965).
212W. C. Pfefferle, J. A. Goff, and J. G. Miller, J. Chern. Phys. 23, 509 238D. G.Miller, Ind. Eng. Chern. Fund. 9, 585 (1970).
(1955). 239 K. M. de Reuck (private communication) (1991).
Downloaded 12 Jan 2013 to 152.14.136.96. Redistribution subject to AIP license or copyright; see http://jpcrd.aip.org/about/rights_and_permissions