Fugacity vapor-phase

In vapor-liquid equilibria, it is relatively easy to start the iteration because assumption of ideal behavior (Raoult s law) provides a reasonable zeroth approximation. By contrast, there is no obvious corresponding method to start the iteration calculation for liquid-liquid equilibria. Further, when two liquid phases are present, we must calculate for each component activity coefficients in two phases since these are often strongly nonlinear functions of compositions, liquid-liquid equilibrium calculations are highly sensitive to small changes in composition. In vapor-liquid equilibria at modest pressures, this sensitivity is lower because vapor-phase fugacity coefficients are usually close to unity and only weak functions of composition. For liquid-liquid equilibria, it is therefore more difficult to construct a numerical iteration procedure that converges both rapidly and consistently. 

In Chapter 2 we discuss briefly the thermodynamic functions whereby the abstract fugacities are related to the measurable, real quantities temperature, pressure, and composition. This formulation is then given more completely in Chapters 3 and 4, which present detailed material on vapor-phase and liquid-phase fugacities, respectively. 

The calculation of vapor and liquid fugacities in multi-component systems has been implemented by a set of computer programs in the form of FORTRAN IV subroutines. These are applicable to systems of up to twenty components, and operate on a thermodynamic data base including parameters for 92 compounds. The set includes subroutines for evaluation of vapor-phase fugacity... 

Equation (1) is of little practical use unless the fuga-cities can be related to the experimentally accessible quantities X, y, T, and P, where x stands for the composition (expressed in mole fraction) of the liquid phase, y for the composition (also expressed in mole fraction) of the vapor phase, T for the absolute temperature, and P for the total pressure, assumed to be the same for both phases. The desired relationship between fugacities and experimentally accessible quantities is facilitated by two auxiliary functions which are given the symbols (f... 

A rigorous relation exists between the fugacity of a component in a vapor phase and the volumetric properties of that phase these properties are conveniently expressed in the form of an equation of state. There are two common types of equations of state one of these expresses the volume as a function of... 

In the calculation of vapor-liquid equilibria, it is necessary to calculate separately the fugacity of each component in each of the two phases. The liquid and vapor phases require different techniques in this chapter we consider calculations for the vapor phase. 

At pressures to a few bars, the vapor phase is at a relatively low density, i.e., on the average, the molecules interact with one another less strongly than do the molecules in the much denser liquid phase. It is therefore a common simplification to assume that all the nonideality in vapor-liquid systems exist in the liquid phase and that the vapor phase can be treated as an ideal gas. This leads to the simple result that the fugacity of component i is given by its partial pressure, i.e. the product of y, the mole fraction of i in the vapor, and P, the total pressure. A somewhat less restrictive simplification is the Lewis fugacity rule which sets the fugacity of i in the vapor mixture proportional to its mole fraction in the vapor phase the constant of proportionality is the fugacity of pure i vapor at the temperature and pressure of the mixture. These simplifications are attractive because they make the calculation of vapor-liquid equilibria much easier the K factors = i i ... 

The fugacity fT of a component i in the vapor phase is related to its mole fraction y in the vapor phase and the total pressure P by the fugacity coefficient ... 

The fugacity coefficient is a function of temperature, total pressure, and composition of the vapor phase it can be calculated from volumetric data for the vapor mixture. For a mixture containing m components, such data are often expressed in the form of an equation of state explicit in the pressure... 

Two additional illustrations are given in Figures 6 and 7 which show fugacity coefficients for two binary systems along the vapor-liquid saturation curve at a total pressure of 1 atm. These results are based on the chemical theory of vapor-phase imperfection and on experimental vapor-liquid equilibrium data for the binary systems. In the system formic acid (1) - acetic acid (2), <() (for y = 1) is lower than formic acid at 100.5°C has a stronger tendency to dimerize than does acetic acid at 118.2°C. Since strong dimerization occurs between all three possible pairs, (fij and not... 

As discussed in Chapter 3, at moderate pressures, vapor-phase nonideality is usually small in comparison to liquid-phase nonideality. However, when associating carboxylic acids are present, vapor-phase nonideality may dominate. These acids dimerize appreciably in the vapor phase even at low pressures fugacity coefficients are well removed from unity. To illustrate. Figures 8 and 9 show observed and calculated vapor-liquid equilibria for two systems containing an associating component. 

To illustrate calculations for a binary system containing a supercritical, condensable component. Figure 12 shows isobaric equilibria for ethane-n-heptane. Using the virial equation for vapor-phase fugacity coefficients, and the UNIQUAC equation for liquid-phase activity coefficients, calculated results give an excellent representation of the data of Kay (1938). In this case,the total pressure is not large and therefore, the mixture is at all times remote from critical conditions. For this binary system, the particular method of calculation used here would not be successful at appreciably higher pressures. 

Figure 13 presents results for a binary where one of the components is a supercritical, noncondensable component. Vapor-phase fugacity coefficients were calculated with the virial... 

Enthalpies are referred to the ideal vapor. The enthalpy of the real vapor is found from zero-pressure heat capacities and from the virial equation of state for non-associated species or, for vapors containing highly dimerized vapors (e.g. organic acids), from the chemical theory of vapor imperfections, as discussed in Chapter 3. For pure components, liquid-phase enthalpies (relative to the ideal vapor) are found from differentiation of the zero-pressure standard-state fugacities these, in turn, are determined from vapor-pressure data, from vapor-phase corrections and liquid-phase densities. If good experimental data are used to determine the standard-state fugacity, the derivative gives enthalpies of liquids to nearly the same precision as that obtained with calorimetric data, and provides reliable heats of vaporization. 

As discussed in Chapter 3, the virial equation is suitable for describing vapor-phase nonidealities of nonassociating (or weakly associating) fluids at moderate densities. Equation (1) gives the second virial coefficient which is used directly in Equation (3-lOb) to calculate the fugacity coefficients. 

These were converted from vapor pressure P to fugacity using the vapor-phase corrections (for pure components), discussed in Chapter 3 then the Poynting correction was applied to adjust to zero pressure ... 

The computer subroutines for calculation of vapor-phase and liquid-phase fugacity (activity) coefficients, reference fugac-ities, and molar enthalpies, as well as vapor-liquid and liquid-liquid equilibrium ratios, are described and listed in this Appendix. These are source routines written in American National Standard FORTRAN (FORTRAN IV), ANSI X3.9-1978, and, as such, should be compatible with most computer systems with FORTRAN IV compilers. Approximate storage requirements and CDC 6400 execution times for these subroutines are given in Appendix J. 

PHIS calculates vapor-phase fugacity coefficients, PHI, for each component in a mixture of N components (N 5. 20) at specified temperature, pressure, and vapor composition. 

PHIS CALCULATES VAPOR PHASE FUGACITY COEFFICIENTS PHI, FOR ALL N... 

CALCULATE VAPOR PHASE FUGACITY COEFFICIENTS FOR ACTUAL COMPOSITION OF... 

In the calculation of vapor phase partial fugacities the use of an equation of state is always justified. In regard to the liquid phase fugacities, there is a choice between two paths ... 

In the case of three-phase equilibria, it is also necessary to account for the solubility of hydrocarbon gases in water. This solubility is proportional to the partial pressure of the hydrocarbon or, more precisely, to its partial fugacity in the vapor phase. The relation which ties the solubility expressed in mole fraction to the fugacity is the following ... 

Effective use of this general equation requires expHcit introduction of the compositions of the phases. This is done either through the activity coefficient, y, or the fugacity coefficient, ( ) Two procedures are in common use. By the gamma—phi approach, activity coefficients for the Hquid phase enter by equation 202 and fugacity coefficients for the vapor phase by equation 164 equation 220 then becomes equation 221 ... 

The second common procedure for VLE calculations is the equation-of-state approach. Here, fugacity coefficients replace the fugacities for both Hquid and vapor phases, and equation 220 becomes equation 226 ... 

Eor nonideal vapor-phase behavior, the fugacity coefficient for component i in the mixture must be determined ... 

The fugacity coefficient departure from nonideaHty in the vapor phase can be evaluated from equations of state or, for approximate work, from fugacity/compressibiHty estimation charts. References 11, 14, and 27 provide valuable insights into this matter. 

Gamma/Phi Approach For many XT E systems of interest the pressure is low enough that a relatively simple equation of state, such as the two-term virial equation, is satisfactoiy for the vapor phase. Liquid-phase behavior, on the other hand, may be conveniently described by an equation for the excess Gibbs energy, from which activity coefficients are derived. The fugacity of species i in the liquid phase is then given by Eq. (4-102), written... 

When Eq. (4-282) is applied to XT E for which the vapor phase is an ideal gas and the liquid phase is an ideal solution, it reduces to a veiy simple expression. For ideal gases, fugacity coefficients and are unity, and the right-hand side of Eq. (4-283) reduces to the Poynting factor. For the systems of interest here this factor is always veiy close to unity, and for practical purposes <1 = 1. For ideal solutions, the activity coefficients are also unity. Equation (4-282) therefore reduces to... 

A generally apphcable alternative to the gamma/phi approach results when both the hquid and vapor phases are described by the same equation of state. The defining equation for the fugacity coefficient, Eq. (4-79), may be applied to each phase ... 

The limits of the Lewis fugacity rule are not determined by pressure but by composition the Lewis rule becomes exact at any pressure in the limit as y( - 1, and therefore it always provides a good approximation for any component i which is present in excess. However, for a component with small mole fraction in the vapor phase, the Lewis rule can sometimes lead to very large errors (P5, R3, RIO). 

Fig. 4. Vapor-phase solubility of naphthalene in ethylene. Data points from G. A. M. Diepen and F. E. C. Scheffer, J. Am. Chem. Soc. 70, 4085 (1948) vapor-phase fugacities from (---) Redlich-Kwong equation (-) Ideal gas law.

Vapor-phase fugacity coefficients are needed not only in high-pressure phase equilibria, but are also of interest in high-pressure chemical equilibria (D6, K7, S4). The equilibrium yield of a chemical reaction can sometimes be strongly influenced by vapor-phase nonideality, especially if reactants and products have small concentrations due to the presence in excess of a suitably chosen nonreactive gaseous solvent (S4). 

In Section II, we discussed the fugacity coefficient, which relates the vapor-phase fugacity to the total pressure and to the composition. The fugacity coefficient can be calculated exactly from an equation of state and, therefore, the problem of calculating vapor-phase fugacities reduces to the problem of... 

---