Fugacity phases

Chapters 2-5 deal with chemical engineering problems that are expressed as algebraic equations - usually sets of nonlinear equations, perhaps thousands of them to be solved together. In Chapter 2 you can study equations of state that are more complicated than the perfect gas law. This is especially important because the equation of state provides the thermodynamic basis for not only volume, but also fugacity (phase equilibrium) and enthalpy (departure from ideal gas enthalpy). Chapter 3 covers vapor-liquid equilibrium, and Chapter 4 covers chemical reaction equilibrium. All these topics are combined in simple process simulation in Chapter 5. This means that you must solve many equations together. These four chapters make extensive use of programming languages in Excel and MATE AB. 

The fugacity function is central to the calculation of phase equilibrium. This should be apparent from the earlier discussion of this chapter and from the calculations of Sec. 7.5, which established that once we had the pure fluid fugacity, phase behavior in a pure fluid could be predicted. Consequently, for the remainder of this chapter we will be concerned with estimating the fugacity of species in gaseous, liquid, and solid mixtures. 

Liquid-state activity-coefficient models, phase equilibria simulation, 448-452 Liquid-state fugacities, phase equilibria simulation, 446 Loans, 260... 

Pure species fugacity (phase equiiibria) at the saturation pressure (b)... 

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

The activity coefficient y relates the liquid-phase fugacity... 

When the same standard-state fugacity is used in both phases. Equation (5) can be rewritten... 

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

The standard-state fugacity of any component must be evaluated at the same temperature as that of the solution, regardless of whether the symmetric or unsymmetric convention is used for activity-coefficient normalization. But what about the pressure At low pressures, the effect of pressure on the thermodynamic properties of condensed phases is negligible and under such con-... 

We can now consider the most convenient form for writing the liquid-phase fugacity of component i. First we consider a condensable component and write... 

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. 

As discussed in Chapter 2, for noncondensable components, the unsymmetric convention is used to normalize activity coefficients. For a noncondensable component i in a multicomponent mixture, we write the fugacity in the liquid phase... 

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. 

We have repeatedly observed that the slowly converging variables in liquid-liquid calculations following the isothermal flash procedure are the mole fractions of the two solvent components in the conjugate liquid phases. In addition, we have found that the mole fractions of these components, as well as those of the other components, follow roughly linear relationships with certain measures of deviation from equilibrium, such as the differences in component activities (or fugacities) in the extract and the raffinate. 

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

At the user s option, one of two methods can be used to calculate the liquid-phase reference fugacity (1) An empirical... 

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

---