Case vapor-liquid equilibrium

The vapor-liquid equilibrium separation calculations considered here are for two cases, isothermal and adiabatic, both at fixed pressure. 

Data on the gas-liquid or vapor-liquid equilibrium for the system at hand. If absorption, stripping, and distillation operations are considered equilibrium-limited processes, which is the usual approach, these data are critical for determining the maximum possible separation. In some cases, the operations are are considerea rate-based (see Sec. 13) but require knowledge of eqmlibrium at the phase interface. Other data required include physical properties such as viscosity and density and thermodynamic properties such as enthalpy. Section 2 deals with sources of such data. 

By using vapor-liquid equilibrium data the above integral can be evaluated numerically. A graphical method is also possible, where a plot of l/(y - xj versus Xr is prepared and the area under the curve over the limits between the initial and fmal mole fraction is determined. However, for special cases the integration can be done analytically. If pressure is constant, the temperature change in the still is small, and the vapor-liquid equilibrium values (K-values, defined as K=y/x for each component) are independent from composition, integration of the Rayleigh equation yields ... 

Figure 9-96. Vapor-liquid equilibrium showing X and application cases referred to in the text. Used by permission, Koshy, T. D., and Rukovena, F. Jr., Hydrocarbon Processing V. 85, No. 5 (1986) p. 64 all rights reserved.

Thermodynamic energy terms (and equilibrium constants) may differ for compounds containing different isotopic species of an element. This effect is described in theoretical detail by Urey (1947), and applications to geochemistry are discussed by Broecker and Oversby (1971) and Faure (1977). A good example is the case of the vapor/liquid equilibrium for water. The vapor pressure of a lighter isotopic species, H2 0, is higher relative to that of heavier species, (or HD O), and others. 

These models are semiempirical and are based on the concept that intermolecular forces will cause nonrandom arrangement of molecules in the mixture. The models account for the arrangement of molecules of different sizes and the preferred orientation of molecules. In each case, the models are fitted to experimental binary vapor-liquid equilibrium data. This gives binary interaction parameters that can be used to predict multicomponent vapor-liquid equilibrium. In the case of the UNIQUAC equation, if experimentally determined vapor-liquid equilibrium data are not available, the Universal Quasi-chemical Functional Group Activity Coefficients (UNIFAC) method can be used to estimate UNIQUAC parameters from the molecular structures of the components in the mixture3. 

In the case of vapor-liquid equilibrium, the vapor and liquid fugacities are equal for all components at the same temperature and pressure, but how can this solution be found In any phase equilibrium calculation, some of the conditions will be fixed. For example, the temperature, pressure and overall composition might be fixed. The task is to find values for the unknown conditions that satisfy the equilibrium relationships. However, this cannot be achieved directly. First, values of the unknown variables must be guessed and checked to see if the equilibrium relationships are satisfied. If not, then the estimates must be modified in the light of the discrepancy in the equilibrium, and iteration continued until the estimates of the unknown variables satisfy the requirements of equilibrium. 

Consider first, a binary mixture of two Components A and B the vapor-liquid equilibrium exhibits only a moderate deviation from ideality, as represented in Figure 4.4a. In this case, as pure A boils at a lower temperature than pure B in the temperature-composition... 

Separation systems include in their mathematical models various vapor-liquid equilibrium (VLE) correlations that are specific to the binary or multicomponent system of interest. Such correlations are usually obtained by fitting VLE data by least squares. The nature of the data can depend on the level of sophistication of the experimental work. In some cases it is only feasible to measure the total pressure of a system as a function of the liquid phase mole fraction (no vapor phase mole fraction data are available). 

To test the validity of the extended Pitzer equation, correlations of vapor-liquid equilibrium data were carried out for three systems. Since the extended Pitzer equation reduces to the Pitzer equation for aqueous strong electrolyte systems, and is consistent with the Setschenow equation for molecular non-electrolytes in aqueous electrolyte systems, the main interest here is aqueous systems with weak electrolytes or partially dissociated electrolytes. The three systems considered are the hydrochloric acid aqueous solution at 298.15°K and concentrations up to 18 molal the NH3-CO2 aqueous solution at 293.15°K and the K2CO3-CO2 aqueous solution of the Hot Carbonate Process. In each case, the chemical equilibrium between all species has been taken into account directly as liquid phase constraints. Significant parameters in the model for each system were identified by a preliminary order of magnitude analysis and adjusted in the vapor-liquid equilibrium data correlation. Detailed discusions and values of physical constants, such as Henry s constants and chemical equilibrium constants, are given in Chen et al. (11). 

A wide variety of data for mean ionic activity coefficients, osmotic coefficients, vapor pressure depression, and vapor-liquid equilibrium of binary and ternary electrolyte systems have been correlated successfully by the local composition model. Some results are shown in Table 1 to Table 10 and Figure 3 to Figure 7. In each case, the chemical equilibrium between the species has been ignored. That is, complete dissociation of strong electrolytes has been assumed. This assumption is not required by the local composition model but has been made here in order to simplify the systems treated. 

Equation 27 is similar to the solid-liquid equilibrium relation used for non-electrolytes. As in the case of the vapor-liquid equilibrium relation for HC1, the solid-liquid equilibrium expression for NaCl is simple since the electrolyte is treated thermodynamically the same in both phases. 

An equilibrium-flash calculation (using the same equations as in case A above) is made at each point in time to find the vapor and liquid flow rates and properties immediately after the pressure letdown valve (the variables with the primes F , F l, y], x j,.. . shown in Fig. 3.8). These two streams are then fed into the vapor and liquid phases. The equations describing the two phases will be similar to Eqs. (3.40) to (3.42) and (3.44) to (3.46) with the addition of (1) a multi-component vapor-liquid equilibrium equation to calculate Pi and (2) NC — 1 component continuity equations for each phase. Controller equations relating 1 to Fi and P to F complete the model. 

Separation processes which involve non-volatile salts arise in two situations. First, as an alternative to extractive or azeotropic distillation, salts may be added to a system to alter the vapor-liquid equilibrium behavior. Second, there are cases where a salt is generated in the process before final product purification. For example, product streams from processes involving esterification, etherification, or neutralization contain salts and are often fed to separation units such as distillation or stripping columns. 

Prediction of salt effect. The procedure for calculation of the preferential solvation number S has been described above. By reversing this procedure, that is, by determining xia from S, we can estimate the salt effect using the vapor-liquid equilibrium without a salt. When the salt concentration is below saturation, the preferential solvation number S can be expressed as follows in cases where the solvation is formed with the first component. 

The salt effects of potassium bromide and a series office symmetrical tetraalkylammonium bromides on vapor-liquid equilibrium at constant pressure in various ethanol-water mixtures were determined. For these systems, the composition of the binary solvent was held constant while the dependence of the equilibrium vapor composition on salt concentration was investigated these studies were done at various fixed compositions of the mixed solvent. Good agreement with the equation of Furter and Johnson was observed for the salts exhibiting either mainly electrostrictive or mainly hydrophobic behavior however, the correlation was unsatisfactory in the case of the one salt (tetraethylammonium bromide) where these two types of solute-solvent interactions were in close competition. The transition from salting out of the ethanol to salting in, observed as the tetraalkylammonium salt series is ascended, was interpreted in terms of the solute-solvent interactions as related to physical properties of the system components, particularly solubilities and surface tensions. 

The data in Tables I-XVI (see Appendix for all tables) show the isobaric vapor-liquid equilibrium results at the boiling point for potassium, ammonium, tetramethylammonium, tetraethylammonium, tetra-n-propylammonium, and tetra-n-butylammonium bromides in various ethanol-water mixtures at fixed liquid composition ratios. The temperature, t, is the boiling temperature for all solutions in these tables. In all cases, the ethanol-water composition was held constant between 0.20 and 0.35 mole fraction ethanol since it is in this range that the most dramatic salt effects on vapor-liquid equilibrium in this particular system should be observed. That is, previous data (12-15,38) have demonstrated that a maximum displacement of the vapor-liquid equilibrium curve by salts frequently occurs in this region. In the results presented here, it should be noted that Equation 1 has been modified to... 

It is easy to see that the BET adsorption isotherm has the correct limits at very high [A] and when multilayer adsorption is negligible. First, consider the case where the pressure of A approaches the value for saturated vapor pressure of A in equilibrium with the liquid. Let the corresponding concentration be designated [A]sa/. The vapor/liquid equilibrium process is written... 

Basically, DESIGNER can use different physical property packages that are easy to interchange with commercial flowsheet simulators. For the case considered, the vapor-liquid equilibrium description is based on the UNIQUAC model. The liquid-phase binary diffusivities are determined using the method of Tyn and Calus (see Ref. 72) for the diluted mixtures, corrected by the Vignes equation (57), to account for finite concentrations. The vapor-phase diffusion coefficients are assumed constant. The reaction kinetics parameters taken from Ref. 202 are implemented directly in the DESIGNER code. 

In Eqs. (3) and (4), the vapor phase mole fractions y and the liquid phase mole fractions %i have to fulfil the vapor-liquid-liquid equilibrium conditions in case of heterogeneous liquid mixtures, and the vapor-liquid equilibrium conditions in homogeneous mixtures. For the latter case, the reaction term 0 in Eq. (5) simplifies to... 

Summing up, the influence of the kinetics of a chemical reaction on the vapor-liquid equilibrium is very complex. Physical distillation boundaries may disappear, while new kinetic stable and unstable nodes may appear. As result, the residue curve map with chemical reaction could look very different from the physical plots. As a consequence, evaluating the kinetic effects on residue curve maps is of great importance for conceptual design of reactive distillation systems. However, if the reaction rate is high enough such that the chemical equilibrium is reached quickly, the RCM simplifies considerably. But even in this case the analysis may be complicated by the occurrence of reactive azeotropes. 

For the specific case of multicomponent vapor/ liquid equilibrium, Eq. (11. becomes... 

Equation (11.22) results as a special case when this relation is applied to vapor/liquid equilibrium of pure species L... 

When liquid and gas phases are both present in an equilibrium mixture of reacting species, Eq. (11.30), a criterion of vapor/liquid equilibrium, must be satisfied along with the equation of chemical-reaction equilibrium. There is considerable choice in the method of treatment of such cases. For example, consider a reaction of gas A and water B to form an aqueous solution C. The reaction may be assumed to occur entirely in the gas phase with simultaneous transfer of material between phases to maintain phase equilibrium. In this case, the equilibrium constant is evaluated from AG° data based on standard states for the species as gases, i.e., the ideal-gas states at 1 bar and the reaction temperature. On the other hand, the reaction may be assumed to occur in the liquid phase, in which case AG° is based on standard states for the species as liquids. Alternatively, the reaction may be written... 

In order to correlate the results obtained, a modified SRK equation of state with Huron-Vidal mixing rules was used. Details about the model are reported in the paper by Soave et al. [16]. This approach is particularly adequated when experimental values of the critical temperature and pressure are not available as it was the case for limonene and linalool. Note that the flexibility of the thermodynamic model to reproduce high-pressure vapor-liquid equilibrium data is ensured by the use of the Huron-Vidal mixing rules and a NRTL activity coefficient model at infinite pressures. Calculation results are reported as continuous curves in figure 2 for the C02-linalool system and in figure 3 for C02-limonene. Note that the same parameters values were used to correlated the data of C02-limonene at 45, 50 e 60 °C. 

Related Calculations. These calculations show how to use vapor-liquid equilibrium data to obtain parameters for activity-coefficient correlations such as those of Van Laar and Wilson. (Use of liquid-liquid equilibrium data for the same purpose is shown in Example 1.20.) If the system forms an azeotrope, the parameters can be obtained from a single measurement of the azeotropic pressure and the composition of the constant boiling mixture. If the activity coefficients at infinite dilution are available, the two parameters for the Van Laar equation are given directly, and the two in the case of the Wilson equation can be solved for as shown in the example. 

In principle, the parameters can be evaluated from minimal experimental data. If vapor-liquid equilibrium data at a series of compositions are available, the parameters in a given excess-free-energy model can be found by numerical regression techniques. The goodness of fit in each case depends on the suitability of the form of the equation. If a plot of GE/X X2RT versus X is nearly linear, use the Margules equation (see Section 3). If a plot of Xi X2RT/GE is linear, then use the Van Laar equation. If neither plot approaches linearity, apply the Wilson equation or some other model with more than two parameters. 

Related Calculations. When vapor-liquid equilibrium data are taken under isobaric rather than isothermal conditions, as is often the case, the right-hand side of the preceding Gibbs-Duhem equation cannot as readily be taken to approximate zero. Instead, the equation should be taken as... 

---