Compressibility factor (gases)

Main Article
Discussion
Related Articles  [?]
Bibliography  [?]
Citable Version  [?]

This editable Main Article has an approved citable version (see its Citable Version subpage). While we have done conscientious work, we cannot guarantee that this Main Article, or its citable version, is wholly free of mistakes. By helping to improve this editable Main Article, you will help the process of generating a new, improved citable version.
(CC) Image: Milton Beychok
Figure 1: Example graphs of gas compressibility factors

The compressibility factor (Z) is a useful thermodynamic property for modifying the ideal gas law to account for behavior of real gases.[1][2][3][4][5] It is a measure of how much the thermodynamic properties of a real gas deviate from those expected of an ideal gas. It may be thought of as the ratio of the actual volume of a real gas to the volume predicted by the ideal gas at the same temperature and pressure as the actual volume.

For an ideal gas, Z always has a value of 1. For real gases, the value may deviate positively or negatively, depending on the effect of the intermolecular forces of the gas. The closer a real gas is to its critical point or to its saturation point, the larger are the deviations of the gas from ideal behavior.

The upper graph in Figure 1 illustrates how the compressibility factor varies for different gases at the same temperature and pressure. The lower graph illustrates how the compressibility factor of a gas (for example, methane) at a given pressure varies with temperature.[1]

This article deals only with the compressibility factor of gases and does not delve into the compressibility of liquids or vapor-liquid mixtures.

Determination of gas compressibility values

The ideal gas law is defined as:

${\displaystyle P\,V_{\mathrm {m} }=R\,T}$

and the ideal gas law corrected for non-ideality is defined is:

${\displaystyle P\,V_{\mathrm {m} }=ZR\,T}$
${\displaystyle P}$ where: = pressure = molar volume of the gas = compressibility factor = Universal gas constant = temperature

and thus:

${\displaystyle Z={\frac {\;P\,V_{\mathrm {m} }}{R\,T}}}$

which is the simplest and most widely used real gas equation of state (EOS). The major limitation of this equation of state is that the gas compressibility factor, Z, is not a constant but varies from one gas to another as well as with the temperature and pressure of the gas under consideration. It must be determined experimentally.

Where experimental data is available for specific gases, that data may be used to produce graphs (such as in Figure 1) of Z versus pressure at a constant temperature or of Z versus pressure for various temperatures for those specific gases. Such graphs are useful for readily obtaining interpolated values of Z between the experimentally determined values.

The compressibility factor, as mentioned earlier, may also be expressed as:

${\displaystyle Z={\frac {V_{\mathrm {actual} }}{V_{\mathrm {ideal} }}}}$

There are three regimes that affect the compressibility factor:

1. the value of Z tends toward 1 as the gas pressure approaches 0, where all gases tend toward ideal behavior
2. the value of Z is less than 1 at intermediate pressures because the intermolecular forces of attraction cause the actual volumes to be less than the ideal values
3. the value of Z is greater than 1 and ultimately tends toward infinity at high pressures because the intermolecular repulsive forces cause the actual volumes to be greater than the ideal values

Generalized compressibility factor graphs for pure gases

The unique relationship between the compressibility factor and the reduced temperature, ${\displaystyle T_{r}}$, and the reduced pressure, ${\displaystyle P_{r}}$, was first recognized by van der Waals in 1873 and is known as the two-parameter principle of corresponding states. The principle of corresponding states expresses the generalization that the properties of a gas which are dependent on intermolecular forces are related to the critical properties of the gas in a universal way. That provides a most important basis for developing correlations of molecular properties.[6]

As for the compressibility of gases, the principle of corresponding states indicates that any pure gas at the same reduced temperature, ${\displaystyle T_{r}}$, and reduced pressure, ${\displaystyle P_{r}}$, should have the same compressibility factor. [3]

The reduced temperature and pressure are defined as:

${\displaystyle T_{r}={\frac {T}{T_{c}}}}$   and   ${\displaystyle P_{r}={\frac {P}{P_{c}}}}$

${\displaystyle T_{c}}$ and ${\displaystyle P_{c}}$ are known as the critical temperature and critical pressure of a gas. They are characteristics of each specific gas with ${\displaystyle T_{c}}$ being the temperature above which it is not possible to liquify a given gas and ${\displaystyle P_{c}}$ is the minimum pressure required to liquify a given gas at its critical temperature.[7] Together they define the critical point of a fluid above which distinct liquid and gas phases of a given fluid do not exist.

The pressure-volume-temperature (PVT) data for real gases varies from one pure gas to another. However, when the compressibility factors of various single-component gases are graphed versus pressure along with temperature isotherms (as shown for methane in the lower graph of Figure 1), many of the graphs exhibit similar isotherm shapes.

(CC) Image: Milton Beychok
Figure 2: Example of a generalized compressibility factor graph[8]

In order to obtain a generalized graph that can be used for many different gases, the reduced pressure and temperature, ${\displaystyle P_{r}}$ and ${\displaystyle T_{r}}$, are used to normalize the compressibility factor data. Figure 2 is an example of a generalized compressibility factor graph derived from hundreds of experimental P-V-T data points of 10 pure gases, namely methane, ethane, ethylene, propane, n-butane, i-pentane, n-hexane, nitrogen, carbon dioxide and steam.[2][4][9][10]

There are more detailed generalized compressibility factor graphs based on as many as 25 or more different pure gases, such as the Nelson-Obert graphs.[11][12][13] Such graphs are said to have an accuracy within 1-2 percent for Z values greater than 0.6 and within 4-6 percent for Z values of 0.3-0.6.

The generalized compressibility factor graphs may be considerably in error for strongly polar gases which are gases for which the centers of positive and negative charge do not coincide. In such cases the estimate for Z may be in error by as much as 15-20 percent.

The quantum gases hydrogen, helium, and neon do not conform to the corresponding-states behavior and the reduced pressure and temperature for those three gases should be redefined in the following manner to improve the accuracy of predicting their compressibility factors when using the generalized graphs:[4][10]

${\displaystyle T_{r}={\frac {T}{T_{c}+8}}}$   and   ${\displaystyle P_{r}={\frac {P}{P_{c}+8}}}$

where the temperatures are in degrees Kelvin and the pressures are in atmospheres.

Generalized compressibility factor graphs for gas mixtures

The generalized compressibility factor graphs for pure gases (such as exemplified in Figure 2 or the Nelson-Obert graphs may also be used for gas mixtures by using Kay's Rule proposed by W. B. Kay in 1936.[14] That rule predicts the P-V-T behavior of a gas mixture by determining the compressibility factor for a gas mixture at the reduced temperature and reduced pressure defined in terms of the pseudocritical temperature, ${\displaystyle T_{pc}}$ and the pseudocritical pressure, ${\displaystyle P_{pc}}$ which are defined as:[15][16]

${\displaystyle T_{pc}=y_{1}\,T_{c_{1}}+y_{2}\,T_{c_{2}}+y_{3}\,T_{c_{3}}\,\dots }$
${\displaystyle P_{pc}=y_{1}\,P_{c_{1}}+y_{2}\,P_{c_{2}}+y_{3}\,P_{c_{3}}\,\dots }$
${\displaystyle T_{pc}}$ where: = the pseudocritical temperature of the gas mixture = the critical temperature of a component of the gas mixture = the pseudocritical pressure of the gas mixture = the critical temperature of a component of the gas = the mole (or volume) fraction of a component of the gas

The pseudoreduced temperature, ${\displaystyle T_{pr}}$, and pseudoreduced pressure, ${\displaystyle P_{pr}}$, to use with the generalized compressibility factor graphs for pure gases are then:

${\displaystyle T_{pr}={\frac {T}{T_{pc}}}}$   and   ${\displaystyle P_{pr}={\frac {P}{P_{pc}}}}$

A compressibility factor graph for natural gas (which is a mixture of hydrocarbon and other gases) using pseudoreduced temperatures and pressures is available on the Internet.[17]

The van der Waals equation

The van der Waals equation was developed in 1873 and may be expressed as:[1][2][3][4][5]

${\displaystyle \left(P+{\frac {a}{V_{\mathrm {m} }^{2}}}\right)\left(V_{\mathrm {m} }-b\right)=R\,T}$

where:

${\displaystyle a}$ is a measure of the strength of attraction between the gas molecules
${\displaystyle b}$ accounts for the volume occupied by gas molecules, which decreases the available open volume

The van der Waals equation may be re-arranged as:

${\displaystyle P={\frac {R\,T}{V_{\mathrm {m} }-b}}-{\frac {a}{V_{\mathrm {m} }^{2}}}}$

and the compressibility factor can be written as:

${\displaystyle Z={\frac {\;P\,V_{\mathrm {m} }}{R\,T}}={\frac {V_{\mathrm {m} }}{R\,T}}\left({\frac {R\,T}{V_{\mathrm {m} }-b}}-{\frac {a}{V_{\mathrm {m} }^{2}}}\right)}$

and we now have an equation for determining ${\displaystyle Z}$ by using the van der Waals parameters ${\displaystyle a}$ and ${\displaystyle b}$:

${\displaystyle Z={\frac {V_{\mathrm {m} }}{V_{\mathrm {m} }-b}}-{\frac {a}{R\,T\,V_{\mathrm {m} }}}}$

Although ${\displaystyle a}$ and ${\displaystyle b}$ are referred to as "the Van der Waals constants", they are not truly constants because they vary from one gas to another; they are, however, independent of P, V and T. In other words, they are constant for the specific gas being considered. Given the critical temperature and pressure for a specific gas, ${\displaystyle a}$ and ${\displaystyle b}$ can be obtained for that specific gas from these equations:[2][3][4][9]

${\displaystyle a={\frac {27\,R^{2}\,T_{c}^{2}}{64\,P_{c}}}}$   and   ${\displaystyle b={\frac {R\,T_{c}}{8\,P_{c}}}}$

For example, the critical pressure of ethane is 4.884 MPa and the critical temperature is 305.43 K (32.3˚C).[4] Using the universal gas constant, ${\displaystyle R}$, of 8.3145 Pa m3 K-1 mol-1, the above equations yield ${\displaystyle a}$ = 0.557 m6 Pa mol-2 and ${\displaystyle b}$ = 6.500×10-5 m3 mol-1.

Tabulated values of ${\displaystyle a}$ and ${\displaystyle b}$ for various gases are available in these references[1][4][18] as well as many other sources in the technical literature.

Modified versions of the van der Waals equation

Determining compressibility factors by using the van der Waals equation of state is an improvement over using the generalized compressibility charts. However, there are modified versions of the van der Waals equation of state which provide even more improvement. Perhaps the most commonly used equations of state by engineers working in petroleum refining, petrochemical production, natural gas processing, cryogenic distillation and related industries are: the Redlich-Kwong equation developed in 1949,[19][20] the Soave-Redlich-Kwong equation developed in 1972[19][20][21] and the Peng-Robinson equation developed in 1976.[19][20] Those three equations are essentially modified versions of the van der Waals equation. The comparative equations are:

Redlich-Kwong
${\displaystyle P={\frac {R\,T}{V_{\mathrm {m} }-b}}-{\frac {a}{V_{\mathrm {m} }\left(V_{\mathrm {m} }+b\right)T^{\,1/2}}}}$
${\displaystyle a={\frac {R^{2}\,T_{c}^{\,5/2}}{2.3393\,P_{c}}}}$   and   ${\displaystyle b={\frac {R\,T_{c}}{11.5420\,P_{c}}}}$
${\displaystyle Z={\frac {V_{\mathrm {m} }}{V_{\mathrm {m} }-b}}-{\frac {a}{R\,T\left(V_{\mathrm {m} }+b\right)T^{\,1/2}}}}$
Soave-Redlich-Kwong
${\displaystyle P={\frac {R\,T}{V_{\mathrm {m} }-b}}-{\frac {a}{V_{\mathrm {m} }\left(V_{\mathrm {m} }+b\right)}}}$
${\displaystyle a={\frac {R^{2}\,T_{c}^{2}}{2.3393\,P_{c}}}\,{\Big (}1+m\,{\Big \{}1-T_{r}^{\,1/2}{\Big \}}{\Big )}^{2}}$   and   ${\displaystyle b={\frac {R\,T_{c}}{11.5420\,P_{c}}}}$
${\displaystyle Z={\frac {V_{\mathrm {m} }}{V_{\mathrm {m} }-b}}-{\frac {a}{R\,T\left(V_{\mathrm {m} }+b\right)}}}$
${\displaystyle m=0.480+1.574\,\omega -0.176\,\omega ^{2}}$

The acentric factor, ${\displaystyle \omega }$, used as a parameter in the Soave-Redlich-Kwong and the Peng-Robinson equations was first introduced by Pitzer et al in 1955.[22] It can be expressed as:[23][24]

${\displaystyle \omega =\log _{10}\left({\frac {P_{c}}{P_{\mathrm {vp} }}}\right)-1}$

where ${\displaystyle P_{\mathrm {vp} }}$ is the saturated vapor pressure of the gas at a temperature of ${\displaystyle T=0.7\,T_{c}}$.

Peng-Robinson
${\displaystyle P={\frac {R\,T}{V_{\mathrm {m} }-b}}-{\frac {a}{V_{\mathrm {m} }^{2}-2\,V_{\mathrm {m} }\,b-b^{2}}}}$
${\displaystyle a={\frac {R^{2}\,T_{c}^{2}}{2.1870\,P_{c}}}\,{\Big (}1+k\,{\Big \{}1-T_{r}^{1/2}{\Big \}}{\Big )}^{2}}$   and   ${\displaystyle b={\frac {R\,T_{c}}{12.8535\,P_{c}}}}$
${\displaystyle Z={\frac {V_{\mathrm {m} }}{V_{\mathrm {m} }-b}}-{\frac {a}{R\,T\,V_{\mathrm {m} }+2\,R\,T\,b-\left(R\,T\,b^{2}/V_{\mathrm {m} }\right)}}}$
${\displaystyle k=0.37464+1.54226\,\omega -0.26922\,\omega ^{2}}$

Virial equation of state

The most fundamental equation of state , in that it has the most sound theoretical foundation, is the virial equation of state:[1][2][3]

${\displaystyle {\frac {P\,V_{\mathrm {m} }}{R\,T}}=Z=1+{\frac {B}{V_{\mathrm {m} }}}+{\frac {C}{V_{\mathrm {m} }^{2}}}+{\frac {D}{V_{\mathrm {m} }^{3}}}+\dots }$

It was first proposed on an empirical basis by Thiesen in 1885[25] and further developed in 1901 by Kamerlingh-Onnes.[26] A statistical-mechanical expression for B was first derived by Ornstein in his Ph.D Thesis.[27] The equation was subsequently re-developed by H.D. Ursell in 1927 on a theoretical basis derived from a statistical-mechanical (cluster) analysis of intermolecular forces.[28]

The parameters ${\displaystyle B}$, ${\displaystyle C}$ and ${\displaystyle D}$ are referred to as the second, third and fourth virial coefficients, respectively. The coefficients are not constants since they vary from one gas to another as well as the temperature of the gas under consideration. They are sometimes written as ${\displaystyle B(T)}$, ${\displaystyle C(T)}$ and ${\displaystyle D(T)}$ to denote that they are functions of temperature. The numerical value of the coefficients must be determined experimentally.

The second virial coefficient provides the largest part of the correction for the non-ideal behavior of a gas. For that reason, the virial equation is sometimes truncated after the term containing the second coefficient. However, when the third coefficient is available, the equation is usually truncated after the term containing the third coefficient. The fourth coefficient is rarely available.

An extensive compilation of virial coefficients for pure gases and gas mixtures was prepared from the virial coefficient database at the Thermodynamics Research Center, formerly at Texas A&M University, College Station, Texas and now located at the National Institute of Standards and Technology (NIST), Boulder, Colorado. The latest version of that compilation was published in 2003.[29] NIST also offers for purchase a software program, NIST Reference Fluid Thermodynamic and Transport Properties Database (REFPROP), that includes second and third virial coefficients as well as a host of other thermodynamic properties.[30]

Other equations of state

There are a number of other equations of state for predicting the behavior of real gases. Perhaps, the best known are the Beattie-Bridgeman equation developed in 1929[31] and the Benedict-Rubin-Webb equation developed in 1940.[32] In general, they are more complex than the equations discussed above.

Notation for the gas constant

The gas constant used in this article is the universal gas constant, ${\displaystyle R}$, that applies to any gas. There is also a specific gas constant, which can be denoted as ${\displaystyle R_{s}}$, and is defined as ${\displaystyle R_{s}=R/M}$ where ${\displaystyle M}$ is the molecular weight.

Unfortunately, many authors in the technical literature sometimes use ${\displaystyle R}$ as the specific gas constant without denoting it as such or stating that it is the specific gas constant. This can and does lead to confusion.

References

1. McQuarrie, Donald A. and Simon, John D. (1999). Molecular Thermodynamics. University Science Books. ISBN 1-891389-05-X.
2. Yunus A. Çengal and Robert H. Turner (2004). Fundamentals of Thermal-Fluid Sciences, 2nd Edition. McGraw-Hill Professional. ISBN 0-07-245426-1.
3. Smith, J.M. et al. (2005). Introduction to Chemical Engineering Thermodynamics, Seventh Edition. McGraw Hill. ISBN 0-07-310445-0.
4. Y.V.C. Rao (1997). Chemical Engineering Thermodynamics. Universities Press (India). ISBN 81-7371-048-1.
5. Gas Laws (Lecture by Professor M.J. Mombourquette at Queens University, Canada)
6. Hong Wei Xiang (2005). The Corresponding States Principle and Its Practice: Thermodynamic, Transport and Surface Properties of Fluids. Elsevier. ISBN 0-444-52062-7.
7. IUPAC Goldbook
8. Gouq-Jen Su (1946). "Modified Law of Corresponding States". Industrial & Engineering Chemistry 38 (8): 803-806. ISSN 0096-4484.
9. Equations of State (Professor Lyes Kadem, Concordia University, Canada)
10. Equations of State (Lecture material used by Professor Sanford Klein at the University of Wisconsin-Madison, U.S.)
11. Richard W. Miller (1996). Flow Measurement Engineering Handbook, 3rd Edition. McGraw-Hill. ISBN 0-07-042366-0.  (See Figure G-1 on page G.145)
12. Kai-Fui Vincent Wong (2000). Thermodynamics for Engineers, 1st Edition. CRC Press. ISBN 0-8493-0232-3.  (See Figure A-2, Page 278}
13. Nelson-Obert Generalized Compressibility Charts (Charts available on the Internet, from the website of San Diego State University, U.S.)
14. W.B.Kay (1936). "Gases and Vapors At High Temperature and Pressure - Density of Hydrocarbon". Industrial Engineering & Chemistry 28:1014. ISSN 0096-4484.
15. Yunus A. Çengal and Michael A. Bowes (2005). Thermodynamics: An Engineering Approach, 5th Edition. McGraw-Hill. ISBN 0-07-2884959.
16. Real Gases (Lecture material used by Professor R.M. Price at Christian Brothers University, U.S.A.)
17. Properties of Natural Gases (Professor Jon Steiner Gudmandsson, Norwegian University of Science and Technology)
18. Robert C. Weast (Editor) (1975). Handbook of Chemistry and Physics, 56th Edition. CRC Press. ISBN 0-87819-455-X.
19. Chapter 3: Equations of State (pages 6-8 of 65 pdf pages) (Part of PhD thesis by Ji Lin Wang, 2002, Swinbourne University of Technology, Australia)
20. Jean Vidal (2003). Thermodynamics: Applications in Chemical Engineering and the Petroleum Industry. Editions Technip. ISBN 2-7108-0800-5.
21. Zakia Nasri and Housam Binous (2007). "Applications of the Soave-Redlich-Kwong Equation of State Using Mathematica". Journal of Chemical Engineering of Japan 40 (6): 534-538. Retrieved from here.
22. K.S. Pitzer et al (1955). "The Volumetric and Thermodynamic Properties of Fluids". Journal of the American Chemical Society 77: 3427-3440.
23. Francis S. Manning and Richard E. Thompson (1991). Oilfield Processing of Petroleum. Pennwell Books. ISBN 0-87814-343-2.
24. Acentric Factor Definition (Physical Property Sources Index, PPSI, of the Swiss Federal Institute of Technology Zurich)
25. Max F. Thiesen (1885). "Investigations of the Equation of State". Annalen der Physik und Chemie (24): 467-492.
26. Heike Kamerlingh-Onnes (1901). "Expression of the Equation of State of Gases and Liquids by Means of Series". Commun. Phys. Lab. Leiden (71).
27. S. L. Ornstein (1908) Toepassing der statistische mechanica van Gibbs op molekulair-theoretische vraagstukken (Application of the statistical mechanics of Gibbs on molecular-theoretical problems), Leiden
28. H. D. Ursell (1927). "The Evaluation of Gibbs' Phase Integral for Imperfect Gases" (23): 685.
29. J.D. Dymond and R.C. Wilhoit (Authors), K.N. Marsch and M. Fenkel (Editors) (2003). Virial Coefficients of Pure gases and Mixtures, 1st Edition. Springer. ISBN 3-540-44340-1.
30. NIST Standard Reference Database 23
31. J.A. Beattie and O.C. Bridgeman (1929). "A New Equation of State for Fluids". Proc. Amer. Acad. Arts Sci. (63): 229-305.
32. M. Benedict, G.B. Webb and L.C. Rubin (1940). "An Empirical Equation for Thermodynamic Properties of Light Hydrocarbons and Their Mixtures". J. Chem. Phys. 8: 334-335.