Compressibility factor (gases)
From Citizendium, the Citizens' Compendium
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 vaporliquid mixtures.
Contents

Determination of gas compressibility values
The ideal gas law is defined as:
and the ideal gas law corrected for nonideality is defined is:
where:  
= pressure  
= molar volume of the gas  
= compressibility factor  
= Universal gas constant  
= temperature 
and thus:
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:
There are three regimes that affect the compressibility factor:
 the value of Z tends toward 1 as the gas pressure approaches 0, where all gases tend toward ideal behavior
 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
 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, , and the reduced pressure, , was first recognized by van der Waals in 1873 and is known as the twoparameter 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, , and reduced pressure, , should have the same compressibility factor. ^{[3]}
The reduced temperature and pressure are defined as:
 and
and are known as the critical temperature and critical pressure of a gas. They are characteristics of each specific gas with being the temperature above which it is not possible to liquify a given gas and 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 pressurevolumetemperature (PVT) data for real gases varies from one pure gas to another. However, when the compressibility factors of various singlecomponent 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.
In order to obtain a generalized graph that can be used for many different gases, the reduced pressure and temperature, and , are used to normalize the compressibility factor data. Figure 2 is an example of a generalized compressibility factor graph derived from hundreds of experimental PVT data points of 10 pure gases, namely methane, ethane, ethylene, propane, nbutane, ipentane, nhexane, 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 NelsonObert graphs.^{[11]}^{[12]}^{[13]} Such graphs are said to have an accuracy within 12 percent for Z values greater than 0.6 and within 46 percent for Z values of 0.30.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 1520 percent.
The quantum gases hydrogen, helium, and neon do not conform to the correspondingstates 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]}
 and
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 NelsonObert 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 PVT 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, and the pseudocritical pressure, which are defined as:^{[15]}^{[16]}
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, , and pseudoreduced pressure, , to use with the generalized compressibility factor graphs for pure gases are then:
 and
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]}
where:
 is a measure of the strength of attraction between the gas molecules
 accounts for the volume occupied by gas molecules, which decreases the available open volume
The van der Waals equation may be rearranged as:
and the compressibility factor can be written as:
and we now have an equation for determining by using the van der Waals parameters and :
Although and 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, and can be obtained for that specific gas from these equations:^{[2]}^{[3]}^{[4]}^{[9]}
 and
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, , of 8.3145 Pa m^{3} K^{1} mol^{1}, the above equations yield = 0.557 m^{6} Pa mol^{2} and = 6.500×10^{5} m^{3} mol^{1}.
Tabulated values of and 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 RedlichKwong equation developed in 1949,^{[19]}^{[20]} the SoaveRedlichKwong equation developed in 1972^{[19]}^{[20]}^{[21]} and the PengRobinson equation developed in 1976.^{[19]}^{[20]} Those three equations are essentially modified versions of the van der Waals equation. The comparative equations are:
RedlichKwong
 and
SoaveRedlichKwong
 and
The acentric factor, , used as a parameter in the SoaveRedlichKwong and the PengRobinson equations was first introduced by Pitzer et al in 1955.^{[22]} It can be expressed as:^{[23]}^{[24]}
where is the saturated vapor pressure of the gas at a temperature of .
PengRobinson
 and
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]}
It was first proposed on an empirical basis by Thiesen in 1885^{[25]} and further developed in 1901 by KamerlinghOnnes.^{[26]} A statisticalmechanical expression for B was first derived by Ornstein in his Ph.D Thesis.^{[27]} The equation was subsequently redeveloped by H.D. Ursell in 1927 on a theoretical basis derived from a statisticalmechanical (cluster) analysis of intermolecular forces.^{[28]}
The parameters , and 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 , and 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 nonideal 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 BeattieBridgeman equation developed in 1929^{[31]} and the BenedictRubinWebb 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, , that applies to any gas. There is also a specific gas constant, which can be denoted as , and is defined as where is the molecular weight.
Unfortunately, many authors in the technical literature sometimes use 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.0} ^{1.1} ^{1.2} ^{1.3} ^{1.4} McQuarrie, Donald A. and Simon, John D. (1999). Molecular Thermodynamics. University Science Books. ISBN 189138905X.
 ↑ ^{2.0} ^{2.1} ^{2.2} ^{2.3} ^{2.4} Yunus A. Çengal and Robert H. Turner (2004). Fundamentals of ThermalFluid Sciences, 2nd Edition. McGrawHill Professional. ISBN 0072454261.
 ↑ ^{3.0} ^{3.1} ^{3.2} ^{3.3} ^{3.4} Smith, J.M. et al. (2005). Introduction to Chemical Engineering Thermodynamics, Seventh Edition. McGraw Hill. ISBN 0073104450.
 ↑ ^{4.0} ^{4.1} ^{4.2} ^{4.3} ^{4.4} ^{4.5} ^{4.6} Y.V.C. Rao (1997). Chemical Engineering Thermodynamics. Universities Press (India). ISBN 8173710481.
 ↑ ^{5.0} ^{5.1} Gas Laws (Lecture by Professor M.J. Mombourquette at Queens University, Canada)
 ↑ Hong Wei Xiang (2005). The Corresponding States Principle and Its Practice: Thermodynamic, Transport and Surface Properties of Fluids. Elsevier. ISBN 0444520627.
 ↑ IUPAC Goldbook
 ↑ GouqJen Su (1946). "Modified Law of Corresponding States". Industrial & Engineering Chemistry 38 (8): 803806. ISSN 00964484.
 ↑ ^{9.0} ^{9.1} Equations of State (Professor Lyes Kadem, Concordia University, Canada)
 ↑ ^{10.0} ^{10.1} Equations of State (Lecture material used by Professor Sanford Klein at the University of WisconsinMadison, U.S.)
 ↑ Richard W. Miller (1996). Flow Measurement Engineering Handbook, 3rd Edition. McGrawHill. ISBN 0070423660. (See Figure G1 on page G.145)
 ↑ KaiFui Vincent Wong (2000). Thermodynamics for Engineers, 1st Edition. CRC Press. ISBN 0849302323. (See Figure A2, Page 278}
 ↑ NelsonObert Generalized Compressibility Charts (Charts available on the Internet, from the website of San Diego State University, U.S.)
 ↑ W.B.Kay (1936). "Gases and Vapors At High Temperature and Pressure  Density of Hydrocarbon". Industrial Engineering & Chemistry 28:1014. ISSN 00964484.
 ↑ Yunus A. Çengal and Michael A. Bowes (2005). Thermodynamics: An Engineering Approach, 5th Edition. McGrawHill. ISBN 0072884959.
 ↑ Real Gases (Lecture material used by Professor R.M. Price at Christian Brothers University, U.S.A.)
 ↑ Properties of Natural Gases (Professor Jon Steiner Gudmandsson, Norwegian University of Science and Technology)
 ↑ Robert C. Weast (Editor) (1975). Handbook of Chemistry and Physics, 56th Edition. CRC Press. ISBN 087819455X.
 ↑ ^{19.0} ^{19.1} ^{19.2} Chapter 3: Equations of State (pages 68 of 65 pdf pages) (Part of PhD thesis by Ji Lin Wang, 2002, Swinbourne University of Technology, Australia)
 ↑ ^{20.0} ^{20.1} ^{20.2} Jean Vidal (2003). Thermodynamics: Applications in Chemical Engineering and the Petroleum Industry. Editions Technip. ISBN 2710808005.
 ↑ Zakia Nasri and Housam Binous (2007). "Applications of the SoaveRedlichKwong Equation of State Using Mathematica". Journal of Chemical Engineering of Japan 40 (6): 534538. Retrieved from here.
 ↑ K.S. Pitzer et al (1955). "The Volumetric and Thermodynamic Properties of Fluids". Journal of the American Chemical Society 77: 34273440.
 ↑ Francis S. Manning and Richard E. Thompson (1991). Oilfield Processing of Petroleum. Pennwell Books. ISBN 0878143432.
 ↑ Acentric Factor Definition (Physical Property Sources Index, PPSI, of the Swiss Federal Institute of Technology Zurich)
 ↑ Max F. Thiesen (1885). "Investigations of the Equation of State". Annalen der Physik und Chemie (24): 467492.
 ↑ Heike KamerlinghOnnes (1901). "Expression of the Equation of State of Gases and Liquids by Means of Series". Commun. Phys. Lab. Leiden (71).
 ↑ S. L. Ornstein (1908) Toepassing der statistische mechanica van Gibbs op molekulairtheoretische vraagstukken (Application of the statistical mechanics of Gibbs on moleculartheoretical problems), Leiden
 ↑ H. D. Ursell (1927). "The Evaluation of Gibbs' Phase Integral for Imperfect Gases" (23): 685.
 ↑ 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 3540443401.
 ↑ NIST Standard Reference Database 23
 ↑ J.A. Beattie and O.C. Bridgeman (1929). "A New Equation of State for Fluids". Proc. Amer. Acad. Arts Sci. (63): 229305.
 ↑ 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: 334335.