Equation UNIQUAC

Earlier experience with the UNIQUAC equation indicated... 

Figure 1 compares data reduction using the modified UNIQUAC equation with that using the original UNIQUAC equation. The data are those of Boublikova and Lu (1969) for ethanol and n-octane. The dashed line indicates results obtained with the original equation (q = q for ethanol) and the continuous line shows results obtained with the modified equation. The original equation predicts a liquid-liquid miscibility gap, contrary to experiment. The modified UNIQUAC equation, however, represents the alcohol/n-octane system with good accuracy. 

Figure 4-1. Effect of UNIQUAC equation modification for an alcohol-hydrocarbon system.

An adequate prediction of multicomponent vapor-liquid equilibria requires an accurate description of the phase equilibria for the binary systems. We have reduced a large body of binary data including a variety of systems containing, for example, alcohols, ethers, ketones, organic acids, water, and hydrocarbons with the UNIQUAC equation. Experience has shown it to do as well as any of the other common models. V7hen all types of mixtures are considered, including partially miscible systems, the... 

Figure 4 shows experimental and predicted phase equilibria for the acetonitrile/benzene system at 45°C. This system exhibits moderate positive deviations from Raoult s law. The high-quality data of Brown and Smith (1955) are very well represented by the UNIQUAC equation. 

Figure 5 shows the isothermal data of Edwards (1962) for n-hexane and nitroethane. This system also exhibits positive deviations from Raoult s law however, these deviations are much larger than those shown in Figure 4. At 45°C the mixture shown in Figure 5 is only 15° above its critical solution temperature. Again, representation with the UNIQUAC equation is excellent.

Figure 7 shows a fit of the UNIQUAC equation to the iso-baric data of Nakanishi et al. (1967) for the methanol-diethyl-amine system this system also exhibits strong negative deviations from Raoult s law. The UNIQUAC equation correctly re-... 

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. 

Table 3 shows results obtained from a five-component, isothermal flash calculation. In this system there are two condensable components (acetone and benzene) and three noncondensable components (hydrogen, carbon monoxide, and methane). Henry s constants for each of the noncondensables were obtained from Equations (18-22) the simplifying assumption for dilute solutions [Equation (17)] was also used for each of the noncondensables. Activity coefficients for both condensable components were calculated with the UNIQUAC equation. For that calculation, all liquid-phase composition variables are on a solute-free basis the only required binary parameters are those for the acetone-benzene system. While no experimental data are available for comparison, the calculated results are probably reliable because all simplifying assumptions are reasonable the... 

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

UNIQUAC equation with binary parameters estimated by supplementing binary VLE data with ternary tie-line data. 

Figure 15 shows results for a difficult type I system methanol-n-heptane-benzene. In this example, the two-phase region is extremely small. The dashed line (a) shows predictions using the original UNIQUAC equation with q = q. This form of the UNIQUAC equation does not adequately fit the binary vapor-liquid equilibrium data for the methanol-benzene system and therefore the ternary predictions are grossly in error. The ternary prediction is much improved with the modified UNIQUAC equation (b) since this equation fits the methanol-benzene system much better. Further improvement (c) is obtained when a few ternary data are used to fix the binary parameters. 

Unfortunately, good binary data are often not available, and no model, including the modified UNIQUAC equation, is entirely adequate. Therefore, we require a calculation method which allows utilization of some ternary data in the parameter estimation such that the ternary system is well represented. A method toward that end is described in the next section. 

Figure 4-16. Representation of ternary liquid-liquid equilibria using the UNIQUAC equation is improved by incorporating ternary tie-line data into binary-parameter estimation. Representation of binary VLB shows small loss of accuracy. ---- Binary...

The ternary diagrams shown in Figure 22 and the selectivi-ties and distribution coefficients shown in Figure 23 indicate very good correlation of the ternary data with the UNIQUAC equation. More important, however, Table 5 shows calculated and experimental quarternary tie-line compositions for five of Henty s twenty measurements. The root-mean-squared deviations for all twenty measurements show excellent agreement between calculated and predicted quarternary equilibria. 

Wittrig, T. S., "The Prediction of Liquid-Liquid Equilibria by the UNIQUAC Equation," B.S. Degree Thesis, University of Illinois, Urbana (1977). 

When the UNIQUAC equation (Chapter 4) is substituted into Equation (16), assuming all parameters and a to be inde-... 

Equations for Liquid-Phase Nonidealities A. Modified UNIQUAC Equation... 

Activity coefficients for condensable components are calculated with the UNIQUAC Equation (4-15)/ and infinite-dilution activity coefficients for noncondensable components are calculated with Equation (4-22). ... 

If the mutual solubilities of the solvents A and B are small, and the systems are dilute in C, the ratio ni can be estimated from the activity coefficients at infinite dilution. The infinite dilution activity coefficients of many organic systems have been correlated in terms of stmctural contributions (24), a method recommended by others (5). In the more general case of nondilute systems where there is significant mutual solubiUty between the two solvents, regular solution theory must be appHed. Several methods of correlation and prediction have been reviewed (23). The universal quasichemical (UNIQUAC) equation has been recommended (25), which uses binary parameters to predict multicomponent equihbria (see Eengineering, chemical DATA correlation). 

The UNIQUAC equation treats g G /RT as comprised of two additive parts, a combinatorial term g, accounting for molecular size and shape differences, and a residual term g (not a residual property), accounting for molecular interactions ... 

Outlined below are the steps required for of a X T.E calciilation of vapor-phase composition and pressure, given the liquid-phase composition and temperature. A choice must be made of an equation of state. Only the Soave/Redlich/Kwong and Peng/Robinson equations, as represented by Eqs. (4-230) and (4-231), are considered here. These two equations usually give comparable results. A choice must also be made of a two-parameter correlating expression to represent the liquid-phase composition dependence of for each pq binaiy. The Wilson, NRTL (with a fixed), and UNIQUAC equations are of general applicabihty for binary systems, the Margules and van Laar equations may also be used. The equation selected depends on evidence of its suitability to the particular system treated. Reasonable estimates of the parameters in the equation must also be known at the temperature of interest. These parameters are directly related to infinite-dilution values of the activity coefficients for each pq binaiy. 

Only the Wilson, NRTL, and UNIQUAC equations are suited to the treatment of multicomponent systems. For such systems, the parameters are determined for pairs of species exactly as for binary systems. 

For most LLE applications, the effect of pressure on the Yi < an be ignored, and thus Eq. (4-327) constitutes a set of N equations relating equilibrium compositions to each other and to temperature. For a given temperature, solution of these equations requires a single expression for the composition dependence of suitable for both liquid phases. Not all expressions for suffice, even in principle, because some cannot represent liquid/liquid phase splitting. The UNIQUAC equation is suitable, and therefore prediction is possible by the UNIFAC method. A special table of parameters for LLE calculations is given by Magnussen, et al. (Jnd E/ig Chem Process Des Dev, 20, pp. 331-339 [1981]). 

---