Compressibility factor explained

In thermodynamics, the compressibility factor (Z), also known as the compression factor or the gas deviation factor, describes the deviation of a real gas from ideal gas behaviour. It is simply defined as the ratio of the molar volume of a gas to the molar volume of an ideal gas at the same temperature and pressure. It is a useful thermodynamic property for modifying the ideal gas law to account for the real gas behaviour.^[1] In general, deviation from ideal behaviour becomes more significant the closer a gas is to a phase change, the lower the temperature or the larger the pressure. Compressibility factor values are usually obtained by calculation from equations of state (EOS), such as the virial equation which take compound-specific empirical constants as input. For a gas that is a mixture of two or more pure gases (air or natural gas, for example), the gas composition must be known before compressibility can be calculated.
Alternatively, the compressibility factor for specific gases can be read from generalized compressibility charts^[1] that plot

Z as a function of pressure at constant temperature.

The compressibility factor should not be confused with the compressibility (also known as coefficient of compressibility or isothermal compressibility) of a material, which is the measure of the relative volume change of a fluid or solid in response to a pressure change.

Definition and physical significance

The compressibility factor is defined in thermodynamics and engineering frequently as:

Z = (p)/(\rho R _{specific} T) , where p is the pressure, \rho is the density of the gas and $R\_\backslash text\; =\; \backslash frac$ is the specific gas constant,^[2] M being the molar mass, and the T is the absolute temperature (kelvin or Rankine scale).

In statistical mechanics the description is:

Z = (pV)/(nRT) where p is the pressure, n is the number of moles of gas, T is the absolute temperature, and R is the gas constant, V is unit Volume

For an ideal gas the compressibility factor is

Z =1 per definition. In many real world applications requirements for accuracy demand that deviations from ideal gas behaviour, i.e., real gas behaviour, be taken into account. The value of Z generally increases with pressure and decreases with temperature. At high pressures molecules are colliding more often. This allows repulsive forces between molecules to have a noticeable effect, making the molar volume of the real gas (V _{m} ) greater than the molar volume of the corresponding ideal gas ((V _{m} ) _{idealgas} = RT/p ), which causes Z to exceed one.^[3] When pressures are lower, the molecules are free to move. In this case attractive forces dominate, making Z< 1 . The closer the gas is to its critical point or its boiling point, the more Z deviates from the ideal case.

Fugacity

The compressibility factor is linked to the fugacity by the relation:

f = P \exp\left(\int (Z -1)/(P) dP\right)

Generalized compressibility factor graphs for pure gases

The unique relationship between the compressibility factor and the reduced temperature,

T _{r} , and the reduced pressure, P _{r} , was first recognized by Johannes Diderik 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.

As for the compressibility of gases, the principle of corresponding states indicates that any pure gas at the same reduced temperature,

T _{r} , and reduced pressure, P _{r} , should have the same compressibility factor.

The reduced temperature and pressure are defined by

T _{r} = (T)/(T _{c}) and P _{r} = (P)/(P _{c}) . Here T _{c} and P _{c} are known as the critical temperature and critical pressure of a gas. They are characteristics of each specific gas with T _{c} being the temperature above which it is not possible to liquify a given gas and P _{c} is the minimum pressure required to liquify a given gas at its critical temperature. 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 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,

P _{r} and 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 PVT data points of 10 pure gases, namely methane, ethane, ethylene, propane, n-butane, i-pentane, n-hexane, nitrogen, carbon dioxide and steam.

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. 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:

T _{r} = (T)/(T _{c} +8) and P _{r} = (P)/(P _{c} +8)

where the temperatures are in kelvins and the pressures are in atmospheres.^[4]

Reading a generalized compressibility chart

In order to read a compressibility chart, the reduced pressure and temperature must be known. If either the reduced pressure or temperature is unknown, the reduced specific volume must be found. Unlike the reduced pressure and temperature, the reduced specific volume is not found by using the critical volume. The reduced specific volume is defined by,

\nu _{R} = (\nu _{actual})/(RT _{cr/P} _{cr})

where

\nu _{actual} is the specific volume.^[5]

Once two of the three reduced properties are found, the compressibility chart can be used. In a compressibility chart, reduced pressure is on the x-axis and Z is on the y-axis. When given the reduced pressure and temperature, find the given pressure on the x-axis. From there, move up on the chart until the given reduced temperature is found. Z is found by looking where those two points intersect. the same process can be followed if reduced specific volume is given with either reduced pressure or temperature.

Observations made from a generalized compressibility chart

There are three observations that can be made when looking at a generalized compressibility chart. These observations are:

1. Gases behave as an ideal gas regardless of temperature when the reduced pressure is much less than one (P_R ≪ 1).
2. When reduced temperature is greater than two (T_R > 2), ideal-gas behavior can be assumed regardless of pressure, unless pressure is much greater than one (P_R ≫ 1).
3. Gases deviate from ideal-gas behavior the most in the vicinity of the critical point.^[6]

Theoretical models

The virial equation is especially useful to describe the causes of non-ideality at a molecular level (very few gases are mono-atomic) as it is derived directly from statistical mechanics:

Z = 1 + (B)/(V _{m}) + (C)/(V ^{2} _{m}) + (D)/(V ^{3} _{m}) + ...

Where the coefficients in the numerator are known as virial coefficients and are functions of temperature.

The virial coefficients account for interactions between successively larger groups of molecules. For example,

B accounts for interactions between pairs, C for interactions between three gas molecules, and so on. Because interactions between large numbers of molecules are rare, the virial equation is usually truncated after the third term.^[7]

When this truncation is assumed, the compressibility factor is linked to the intermolecular-force potential φ by:

Z = 1 + 2 \pi (N _{A})/(V _{m}) \int ^{inf ty} _{0} \left( 1 - \exp \left( (\varphi)/(kT) \right)\right) r ^{2} dr

The Real gas article features more theoretical methods to compute compressibility factors.

Physical mechanism of temperature and pressure dependence

Deviations of the compressibility factor, Z, from unity are due to attractive and repulsive intermolecular forces. At a given temperature and pressure, repulsive forces tend to make the volume larger than for an ideal gas; when these forces dominate Z is greater than unity. When attractive forces dominate, Z is less than unity. The relative importance of attractive forces decreases as temperature increases (see effect on gases).

As seen above, the behavior of Z is qualitatively similar for all gases. Molecular nitrogen, N, is used here to further describe and understand that behavior. All data used in this section were obtained from the NIST Chemistry WebBook.^[8] It is useful to note that for N the normal boiling point of the liquid is 77.4 K and the critical point is at 126.2 K and 34.0 bar.

The figure on the right shows an overview covering a wide temperature range. At low temperature (100 K), the curve has a characteristic check-mark shape, the rising portion of the curve is very nearly directly proportional to pressure. At intermediate temperature (160 K), there is a smooth curve with a broad minimum; although the high pressure portion is again nearly linear, it is no longer directly proportional to pressure. Finally, at high temperature (400 K), Z is above unity at all pressures. For all curves, Z approaches the ideal gas value of unity at low pressure and exceeds that value at very high pressure.

To better understand these curves, a closer look at the behavior for low temperature and pressure is given in the second figure. All of the curves start out with Z equal to unity at zero pressure and Z initially decreases as pressure increases. N is a gas under these conditions, so the distance between molecules is large, but becomes smaller as pressure increases. This increases the attractive interactions between molecules, pulling the molecules closer together and causing the volume to be less than for an ideal gas at the same temperature and pressure. Higher temperature reduces the effect of the attractive interactions and the gas behaves in a more nearly ideal manner.

As the pressure increases, the gas eventually reaches the gas-liquid coexistence curve, shown by the dashed line in the figure. When that happens, the attractive interactions have become strong enough to overcome the tendency of thermal motion to cause the molecules to spread out; so the gas condenses to form a liquid. Points on the vertical portions of the curves correspond to N₂ being partly gas and partly liquid. On the coexistence curve, there are then two possible values for Z, a larger one corresponding to the gas and a smaller value corresponding to the liquid. Once all the gas has been converted to liquid, the volume decreases only slightly with further increases in pressure; then Z is very nearly proportional to pressure.

As temperature and pressure increase along the coexistence curve, the gas becomes more like a liquid and the liquid becomes more like a gas. At the critical point, the two are the same. So for temperatures above the critical temperature (126.2 K), there is no phase transition; as pressure increases the gas gradually transforms into something more like a liquid. Just above the critical point there is a range of pressure for which Z drops quite rapidly (see the 130 K curve), but at higher temperatures the process is entirely gradual.

The final figures shows the behavior at temperatures well above the critical temperatures. The repulsive interactions are essentially unaffected by temperature, but the attractive interaction have less and less influence. Thus, at sufficiently high temperature, the repulsive interactions dominate at all pressures.

This can be seen in the graph showing the high temperature behavior. As temperature increases, the initial slope becomes less negative, the pressure at which Z is a minimum gets smaller, and the pressure at which repulsive interactions start to dominate, i.e. where Z goes from less than unity to greater than unity, gets smaller. At the Boyle temperature (327 K for N), the attractive and repulsive effects cancel each other at low pressure. Then Z remains at the ideal gas value of unity up to pressures of several tens of bar. Above the Boyle temperature, the compressibility factor is always greater than unity and increases slowly but steadily as pressure increases.

Experimental values

It is extremely difficult to generalize at what pressures or temperatures the deviation from the ideal gas becomes important. As a rule of thumb, the ideal gas law is reasonably accurate up to a pressure of about 2 atm, and even higher for small non-associating molecules. For example, methyl chloride, a highly polar molecule and therefore with significant intermolecular forces, the experimental value for the compressibility factor is

Z =0 . 9152 at a pressure of 10 atm and temperature of 100 °C.^[9] For air (small non-polar molecules) at approximately the same conditions, the compressibility factor is only Z =1 . 0025 (see table below for 10 bars, 400 K).

Compressibility of air

Normal air comprises in crude numbers 80 percent nitrogen and 20 percent oxygen . Both molecules are small and non-polar (and therefore non-associating). We can therefore expect that the behaviour of air within broad temperature and pressure ranges can be approximated as an ideal gas with reasonable accuracy. Experimental values for the compressibility factor confirm this.

Compressibility factor for air (experimental values)^[10]

(K)	Pressure, absolute (bar)
	1 ! 5	10 ! 20	40 ! 60	80 ! 100	150 ! 200	250 ! 300	400 ! 500
75	0.0052	0.0260	0.0519	0.1036	0.2063	0.3082	0.4094	0.5099	0.7581	1.0125
80		0.0250	0.0499	0.0995	0.1981	0.2958	0.3927	0.4887	0.7258	0.9588	1.1931	1.4139
90	0.9764	0.0236	0.0453	0.0940	0.1866	0.2781	0.3686	0.4681	0.6779	0.8929	1.1098	1.3110	1.7161	2.1105
100	0.9797	0.8872	0.0453	0.0900	0.1782	0.2635	0.3498	0.4337	0.6386	0.8377	1.0395	1.2227	1.5937	1.9536
120	0.9880	0.9373	0.8860	0.6730	0.1778	0.2557	0.3371	0.4132	0.5964	0.7720	0.9530	1.1076	1.5091	1.7366
140	0.9927	0.9614	0.9205	0.8297	0.5856	0.3313	0.3737	0.4340	0.5909	0.7699	0.9114	1.0393	1.3202	1.5903
160	0.9951	0.9748	0.9489	0.8954	0.7803	0.6603	0.5696	0.5489	0.6340	0.7564	0.8840	1.0105	1.2585	1.4970
180	0.9967	0.9832	0.9660	0.9314	0.8625	0.7977	0.7432	0.7084	0.7180	0.7986	0.9000	1.0068	1.2232	1.4361
200	0.9978	0.9886	0.9767	0.9539	0.9100	0.8701	0.8374	0.8142	0.8061	0.8549	0.9311	1.0185	1.2054	1.3944
250	0.9992	0.9957	0.9911	0.9822	0.9671	0.9549	0.9463	0.9411	0.9450	0.9713	1.0152	1.0702	1.1990	1.3392
300	0.9999	0.9987	0.9974	0.9950	0.9917	0.9901	0.9903	0.9930	1.0074	1.0326	1.0669	1.1089	1.2073	1.3163
350	1.0000	1.0002	1.0004	1.0014	1.0038	1.0075	1.0121	1.0183	1.0377	1.0635	1.0947	1.1303	1.2116	1.3015
400	1.0002	1.0012	1.0025	1.0046	1.0100	1.0159	1.0229	1.0312	1.0533	1.0795	1.1087	1.1411	1.2117	1.2890
450	1.0003	1.0016	1.0034	1.0063	1.0133	1.0210	1.0287	1.0374	1.0614	1.0913	1.1183	1.1463	1.2090	1.2778
500	1.0003	1.0020	1.0034	1.0074	1.0151	1.0234	1.0323	1.0410	1.0650	1.0913	1.1183	1.1463	1.2051	1.2667
600	1.0004	1.0022	1.0039	1.0081	1.0164	1.0253	1.0340	1.0434	1.0678	1.0920	1.1172	1.1427	1.1947	1.2475
800	1.0004	1.0020	1.0038	1.0077	1.0157	1.0240	1.0321	1.0408	1.0621	1.0844	1.1061	1.1283	1.1720	1.2150
1000	1.0004	1.0018	1.0037	1.0068	1.0142	1.0215	1.0290	1.0365	1.0556	1.0744	1.0948	1.1131	1.1515	1.1889

Z values are calculated from values of pressure, volume (or density), and temperature in Vasserman, Kazavchinskii, and Rabinovich, "Thermophysical Properties of Air and Air Components;' Moscow, Nauka, 1966, and NBS-NSF Trans. TT 70-50095, 1971: and Vasserman and Rabinovich, "Thermophysical Properties of Liquid Air and Its Component, "Moscow, 1968, and NBS-NSF Trans. 69-55092, 1970.

See also

• Fugacity
• Real gas
• Theorem of corresponding states
• Van der Waals equation

External links

• Compressibility factor (gases) A Citizendium article.
• Real Gases includes a discussion of compressibility factors.

Notes and References

1. http://iptibm1.ipt.ntnu.no/~jsg/undervisning/naturgass/parlaktuna/Chap3.pdf Properties of Natural Gases
2. Book: Zucker, Robert D. . Biblarz, Oscar. Fundamentals of Gas Dynamics. 2nd. Wiley Books . 2002. 0-471-05967-6. page 327
3. Book: McQuarrie, Donald A. . Simon, John D.. Molecular Thermodynamics. University Science Books . 1999. 1-891389-05-X. page 55
4. Book: Y.V.C. Rao. Chemical Engineering Thermodynamics. Universities Press (India). 1997. 81-7371-048-1.
5. Book: Cengel, Yunus A. . Boles, Michael A.. Thermodynamics: An Engineering Approach, Eighth Edition. McGraw-Hill Education . 2015. 978-0-07-339817-4. page 140
6. Book: Cengel, Yunus A. . Boles, Michael A.. Thermodynamics: An Engineering Approach, Eighth Edition. McGraw-Hill Education . 2015. 978-0-07-339817-4. page 139
7. Book: Smith, J.M.. Introduction to Chemical Engineering Thermodynamics. Seventh. McGraw Hill. 2005. 0-07-310445-0. etal. page73
8. http://webbook.nist.gov/chemistry/fluid/ NIST Chemistry WebBook
9. Book: Perry's chemical engineers' handbook. registration. 6. MCGraw-Hill . 1984. 0-07-049479-7. page 3-268
10. Book: Perry's chemical engineers' handbook . registration . 6 . MCGraw-Hill . 1984 . 0-07-049479-7 . 3-162 .