Givens

In vapor-liquid equilibria, if one phase composition is given, there are basically four types of problems, characterized by those variables which are specified and those which are to be calculated. Let T stand for temperature, P for total pressure, for the mole fraction of component i in the liquid phase, and y for the mole fraction of component i in the vapor phase. For a mixture containing m components, the four types can be organized in this way ... 

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. 

Extensive data and bibliography on hydrocarbons a few predictive correlations are also given. 

Comprehensive data collection for more than 6000 binary and multicomponent mixtures at moderate pressures. Data correlation and consistency tests are given for each data set. 

Discusses the thermodynamic basis for computer calculations for vapor-liquid equilibria computer programs are given. Now out of date. 

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

It is strictly for convenience that certain conventions have been adopted in the choice of a standard-state fugacity. These conventions, in turn, result from two important considerations (a) the necessity for an unambiguous thermodynamic treatment of noncondensable components in liquid solutions, and (b) the relation between activity coefficients given by the Gibbs-Duhem equation. The first of these considerations leads to a normalization for activity coefficients for nonoondensable components which is different from that used for condensable components, and the second leads to the definition and use of adjusted or pressure-independent activity coefficients. These considerations and their consequences are discussed in the following paragraphs. 

For a noncondensable component, therefore, it is convenient to use a normalization different from that given by Equation (13) in its place we use... 

Normally, Henry s constant for solute 2 in solvent 1 is determined experimentally at the solvent vapor pressure Pj. The effect of pressure on Henry s constant is given by... 

Because of the approximation given by Equation (22), we obtain a convenient method for determining f for a noncondensable... 

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

Figure 1 shows second virial coefficients for four pure fluids as a function of temperature. Second virial coefficients for typical fluids are negative and increasingly so as the temperature falls only at the Boyle point, when the temperature is about 2.5 times the critical, does the second virial coefficient become positive. At a given temperature below the Boyle point, the magnitude of the second virial coefficient increases with... 

As shown elsewhere (Nothnagel et al., 1973), the fugacity coefficient of component i is given by... 

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

Since parameters for many fluids of interest are not given in this monograph, it may be necessary to estimate the required parameters T, P, R, y, and n. ... 

Values of Rj, probably close to the required accuracy, can be estimated from the parachor, P the parachor can be calculated from a group-contribution method given by Reid et al. The... 

Details for calculating fugacity coefficients are given in Appendix A. 

When a condensable solute is present, the activity coefficient of a solvent is given by Equation (15) provided that all composition variables (x, 9, and ) are taicen on an (all) solute-free basis. Composition variables 9 and 4 are automatically on a solute-free basis by setting q = q = r = 0 for every solute. 

Since we make the simplifying assumption that the partial molar volumes are functions only of temperature, we assume that, for our purposes, pressure has no effect on liquid-liquid equilibria. Therefore, in Equation (23), pressure is not a variable. The activity coefficients depend only on temperature and composition. As for vapor-liquid equilibria, the activity coefficients used here are given by the UNIQUAC equation. Equation (15). ... 

To illustrate, predictions were first made for a ternary system of type II, using binary data only. Figure 14 compares calculated and experimental phase behavior for the system 2,2,4-trimethylpentane-furfural-cyclohexane. UNIQUAC parameters are given in Table 4. As expected for a type II system, agreement is good. 

The estimated true values must satisfy the appropriate equilibrium constraints. For points 1 through L, there are two constraints given by Equation (2-4) one each for components 1 and 2. For points L+1 through M the same equilibrium relations apply however, now they apply to components 2 and 3. The constraints for the tie-line points, M+1 through N, are given by Equation (2-6), applied to each of the three components. 

For a real vapor mixture, there is a deviation from the ideal enthalpy that can be calculated from an equation of state. The enthalpy of the real vapor is given by... 

Since attention is here confined to moderate pressures, the last term in Equation (15) can be neglected. The first term in Equation (15) is given by Equation (5), with x s replacing y s. 

There is a significant difference between the results shown in Figure 2 and calculated results given in Brit. Chem. Eng. Proc. Tech. 16 1036 (1971). We believe the latter to be in error. 

Two generally accepted models for the vapor phase were discussed in Chapter 3 and one particular model for the liquid phase (UNIQUAC) was discussed in Chapter 4. Unfortunately, these, and all other presently available models, are only approximate when used to calculate equilibrium properties of dense fluid mixtures. Therefore, any such model must contain a number of adjustable parameters, which can only be obtained from experimental measurements. The predictions of the model may be sensitive to the values selected for model parameters, and the data available may contain significant measurement errors. Thus, it is of major importance that serious consideration be given to the proper treatment of experimental measurements for mixtures to obtain the most appropriate values for parameters in models such as UNIQUAC. 

The method used here is based on a general application of the maximum-likelihood principle. A rigorous discussion is given by Bard (1974) on nonlinear-parameter estimation based on the maximum-likelihood principle. The most important feature of this method is that it attempts properly to account for all measurement errors. A discussion of the background of this method and details of its implementation are given by Anderson et al. (1978). 

This reduces the calculation at each step to solution of a set of linear equations. The program description and listing are given in Appendix H. 

An apparent systematic error may be due to an erroneous value of one or both of the pure-component vapor pressures as discussed by several authors (Van Ness et al., 1973 Fabries and Renon, 1975 Abbott and Van Ness, 1977). In some cases, highly inaccurate estimates of binary parameters may occur. Fabries and Renon recommend that when no pure-component vapor-pressure data are given, or if the given values appear to be of doubtful validity, then the unknown vapor pressure should be included as one of the adjustable parameters. If, after making these corrections, the residuals again display a nonrandom pattern, then it is likely that there is systematic error present in the measurements. ... 

---