Determination of vapor-liquid equilibria

Very important for the practical determination of vapor-liquid equilibria are collections of experimental vapor-hquid-equilibrium (VLE) data. The collection of (Gmehling and Onken 197711) consists of about 20 volumes containing vapor-liquid equilibria of nearly all binary, ternary, and quaternary mixtures pirbhshed so far. A collection of azeotropic data (Gmehling et al. 2004) lists data of 18,000 systems involving approximately 1,700 compounds. 

Fig. 1-17. LABODEST apparatus for the determination of vapor-liquid equilibria. Vacmun, ambient and over pressure operation at temperatures up to 250 C. Presentation according to data of Fischer Labor- und Verfahrenstechnik, Meckenheim near Bonn [1.61].

D]mamic Flow Method. Another method that has been widely used (Refs. 10, 21, 26, 37) for the determination of vapor-liquid equilibria is one in which a vapor is passed through a series of vessels containing liquids of a suitable composition. The vapor entering the first vessel may be of a composition somewhat different from the equilibrium vapor, but as it passes through the nstem it tends to approach equilibrium. If all the vessels have approximately the same liquid composition, the vapor will more nearly approach equilibrium as it passes through the unit. The number of vessels employed should be such that the vapor entering the last unit is of essentially equilibrium composition. 

Reichl, A. Daiminger, U. Schmidt, A. Davies, M. Hoffmann, U. Brinkmeier, C. Reder, C. Mar-quardt, W. A non-recycle flow still for the experimental determination of vapor-liquid equilibria in reactive systems. Fluid Phase Equilib. 1998, 153, 113-134. 

Yoshikawa, Y Takagi, A. Kato, M. Indirect determination of vapor-liquid equilibria by a small ebulli-ometer. Tetrahydrofuran - alcohol binary systems. J. Chem. Eng. Data 1980, 25, 344-346. 

The most important aspect of the simulation is that the thermodynamic data of the chemicals be modeled correctly. It is necessary to decide what equation of state to use for the vapor phase (ideal gas, Redlich-Kwong-Soave, Peng-Robinson, etc.) and what model to use for liquid activity coefficients [ideal solutions, solubility parameters, Wilson equation, nonrandom two liquid (NRTL), UNIFAC, etc.]. See Sec. 4, Thermodynamics. It is necessary to consider mixtures of chemicals, and the interaction parameters must be predictable. The best case is to determine them from data, and the next-best case is to use correlations based on the molecular weight, structure, and normal boiling point. To validate the model, the computer results of vapor-liquid equilibria could be checked against experimental data to ensure their validity before the data are used in more complicated computer calculations. 

The proper design of distillation and absorption columns depends on knowledge of vapor—liquid equilibrium, as do flash calculations used to determine the physical state of streams at given conditions of temperature, pressure, and composition. Detailed treatments of vapor—liquid equilibria are available (6,7). 

In this study, a thermodynamic framework has been presented for the calculation of vapor-liquid equilibria for binary solvents containing nonvolatile salts. From an appropriate definition of a pseudobinary system, infinite dilution activity coefficients for the salt-containing system may be estimated from a knowledge of vapor pressure lowering, salt-free infinite dilution activity coefficients, and a single system-dependent constant. Parameters for the Wilson equation may be determined from the infinite dilution activity coefficients. 

In order to evaluate each of the derivatives, such quantities as (V" — V-), (S l — Sj), and (dfi t/x t)T P need to be evaluated. The difference in the partial molar volumes of a component between the two phases presents no problem the dependence of the molar volume of a phase on the mole fraction must be known from experiment or from an equation of state for a gas phase. In order to determine the difference in the partial molar entropies, not only must the dependence of the molar entropy of a phase on the mole fraction be known, but also the difference in the molar entropy of the component in the two standard states must be known or calculable. If the two standard states are the same, there is no problem. If the two standard states are the pure component in the two phases at the temperature and pressure at which the derivative is to be evaluated, the difference can be calculated by methods similar to that discussed in Sections 10.10 and 10.12. In the case of vapor-liquid equilibria in which the reference state of a solute is taken as the infinitely dilute solution, the difference between the molar entropy of the solute in its two standard states may be determined from the temperature dependence of the Henry s law constant. Finally, the expression used for fii in evaluating (dx Jdx l)TtP must be appropriate for the particular phase of interest. This phase is dictated by the particular choice of the mole fraction variables. 

There are three independent variables in coexisting equilibrium vapor/liquid systems, namely temperature, pressure, and fraction liquid (or vapor). If two of these are specified in a problem, the third is determined by the phase behavior of the system. There are seven types of vapor/liquid equilibria calculations in our program, as in Figure 1 under "Single Stage Calculation."... 

Consider the separation of 100 kmol/hr of an equimolar stream of tetrahydrofuran (THF) and water using pressure-swing distillation, as shown in Figure 7.37. The tower T1 operates at 1 bar, with the pressure of the tower T2 increased to 10 bar. As shown in the T-x-y diagrams in Figure 7.38, the binary azeotrope shifts from 19 mol% water at 1 bar to 33 mol% water at 10 bar. Assume that the bottoms product from T1 contains pure water and that from D2 contains pure THF. Also, assume that the distillates from T1 and T2 are at their azeotropic compositions. Determine the unknown flow rates of the product and internal streams. Note that data for the calculation of vapor-liquid equilibria are provided in Table 7.5. 

It is not easy to determine which factors play the greatest role in obtaining good accuracy and precision. One must consider the assumptions inherent in the theory as well as the chemical, mechanical, and instrumental parameters. In general, gas chromatographic methods agree within 1-5% with other physicochemical methods. For example, Hussam and Carr (22) showed that in the measurement of vapor/liquid equilibria via headspace GC, complex thermodynamic and analytical correction factors were needed. These often came from other experimental measurements that were not necessarily accurately known. Another source of significant error can be in determination of the mass of stationary phase contained within the column (59). Other sources of error include measurement of holdup time (60), flowrate, sample mass, response factors, peak area, or baseline fidelity. 

In Equation (24), a is the estimated standard deviation for each of the measured variables, i.e. pressure, temperature, and liquid-phase and vapor-phase compositions. The values assigned to a determine the relative weighting between the tieline data and the vapor-liquid equilibrium data this weighting determines how well the ternary system is represented. This weighting depends first, on the estimated accuracy of the ternary data, relative to that of the binary vapor-liquid data and second, on how remote the temperature of the binary data is from that of the ternary data and finally, on how important in a design the liquid-liquid equilibria are relative to the vapor-liquid equilibria. Typical values which we use in data reduction are Op = 1 mm Hg, = 0.05°C, = 0.001, and = 0.003... 

In modern separation design, a significant part of many phase-equilibrium calculations is the mathematical representation of pure-component and mixture enthalpies. Enthalpy estimates are important not only for determination of heat loads, but also for adiabatic flash and distillation computations. Further, mixture enthalpy data, when available, are useful for extending vapor-liquid equilibria to higher (or lower) temperatures, through the Gibbs-Helmholtz equation. ... 

Griswold, J., D. Andres and V. A. Klein, Determination of High Pressure Vapor-Liquid Equilibria. The Vapor-Liquid Equilibrium of Benzene-Toluene," Trans. Amer. Inst. Chem. Engrs. 39, 223 (1943). 

Binary (vapor + liquid) equilibria studies involve the determination of / as a function of composition. the mole fraction in the liquid phase. Of special interest is the dependence of/ on composition in the limit of infinite dilution. In the examples which follow, equilibrium vapor pressures, p,. are measured and described. These vapor pressures can be corrected to vapor fugacities using the techniques described in the previous section. As stated earlier, at the low pressures involved in most experiments, the difference between p, and / is very small, and we will ignore it unless a specific application requires a differentiation between the two. 

There are many types of phase diagrams in addition to the two cases presented here these are summarized in detail by Zief and Wilcox (op. cit., p. 21). Solid-liquid phase equilibria must be determined experimentally for most binary and multicomponent systems. Predictive methods are based mostly on ideal phase behavior and have limited accuracy near eutectics. A predictive technique based on extracting liquid-phase activity coefficients from vapor-liquid equilibria that is useful for estimating nonideal binary or multicomponent solid-liquid phase behavior has been reported by Muir (Pap. 71f, 73d ann. meet., AIChE, Chicago, 1980). 

Kieckbusch, T.G. and King, C.J. An improved method of determining vapor-liquid equilibria for dilute organics in aqneons solution, J. Chromatogr. Sci., 17 273-276, 1979. 

The displacement of the azeotropic composition by progressive addition of magnesium nitrate has been shown in Table I above. Vapor-liquid equilibria have been determined (3, 4). Figure 5 depicts equilibrium vapor compositions in the ternary system at the boiling point, while Figure 6 shows boiling points in the system at 760 mm Hg. 

Three different uses of the Gibbs-Duhem equation associated with the integral method are discussed in this section (A) the calculation of the excess chemical potential of one component when that of the other component is known (B) the determination of the minimum number of intensive variables that must be measured in a study of isothermal vapor-liquid equilibria and the calculation of the values of other variables and (C) the study of the thermodynamic consistency of the data when the data are redundant. 

Puck, T.T., Wise, H. (1946) Studies in vapor-liquid equilibria. I. A new dynamic method for the determination of vapor pressures of liquids. J. Phys. Chem. 50, 329-339. 

Table 3. Approximate Costs of Binary Vapor-Liquid Equilibria Determination"...

A high-pressure circulation-type apparatus was designed and constructed to investigate the vapor-liquid equilibria (VLE) of systems containing limonene, linalool and supercritical carbon dioxide. VLE data of binary and ternary systems of these compounds can be determined in the ranges of pressure and temperature of interest for the deterpenation of cold-pressed orange oil. The preliminary results obtained for the binaries CC -linalool and C02-limonene were compared to data already published with acceptable accuracy and well correlated by a modified Soave-Redlich-Kwong (SRK) equation of state. 

In this work, the basis for a proper answer to the selectivity problem between the two cited compounds is set. For this purpose, a high-pressure vapor-liquid equilibria circulation-type apparatus was designed and constructed. Some vapor-liquid equilibria (VLE) data of the binary systems C02-limonene and C02-linalool were determined and compared with data available in the literature. The results obtained were accurately correlated by a modified SRK equation of state that avoids the use of critical constants [16]. 

Momentary mass of a sample may be derived from momentary weight only if the density of the gas, and thus the buoyancy of the sample, are known. Volumetric data of pure gases are calculated from precise equations of state, as they exist for carbon dioxide and other gases, or taken from tables. Cubic equations of state are used to calculate densities of gas mixtures. We have always employed van-der-Waals mixing rules and fitted the interaction parameter to vapor-liquid equilibria determined by ourselves or taken from literature. 

---