Vapor-liquid equilibrium programming

VPLQFT is a computer program for correlating binary vapor-liquid equilibrium (VLE) data at low to moderate pressures. For such binary mixtures, the truncated virial equation of state is used to correct for vapor-phase nonidealities, except for mixtures containing organic acids where the "chemical" theory is used. The Hayden-0 Connell (1975) correlation gives either the second virial coefficients or the dimerization equilibrium constants, as required. 

The research programs on extractive distillation by salt effect and on salt effect in vapor-liquid equilibrium at the Royal Military College of Canada are supported by the Defence Research Board of Canada, Grant No. 9530-142. 

The current state-of-the-art is such that there are no reliable methods of predicting liquid-liquid equilibria of polymer-solvent systems. Thus, the recommended procedures and computer programs included in this Handbook treat only vapor-liquid equilibrium. A discussion of the correlation of LLE data is included in Chapter 2. 

Special correlations have been developed to handle the vapor-liquid equilibrium of such systems, and these are incorporated in most design and simulation programs. 

Write a Fortran computer program for batch distillation. Calculate the mass fraction composition in the pot and the mass retrieved from the condenser given the initial mass and composition in the pot, and the final mass in the pot. Use Example 6.3 for the vapor-liquid equilibrium data and the general setup. 

In using simulation software, it is important to keep in mind that the quality of the results is highly dependent upon the quahty of the liquid-liquid equilibrium (LLE) model programmed into the simulation. In most cases, an experimentally vmidated model will be needed because UNIFAC and other estimation methods are not sufficiently accurate. It also is important to recognize, as mentioned in earlier discussions, that binary interaction parameters determined by regression of vapor-liquid equilibrium (VLE) data cannot be rehed upon to accurately model the LLE behavior for the same system. On the other hand, a set of binary interaction parameters that model LLE behavior properly often will provide a reasonable VLE fit for the same system—because pure-component vapor pressures often dominate the calculation of VLE. 

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. 

Using one of these activity coefficient equations it is possible to calculate liquid-liquid equihbrium (LLE) behavior of multicomponent hquid systems. Consider, for example, the ternary system of Figure 1. A system of overall composition A splits into two liquid phases B and C. The calculation of compositions of B and C is analogous to the flash ciculation of vapor-liquid equilibrium problems. By using the UNIQUAC equations to obtain the partition coefficients, Kj, this problem can be solved for any composition A of the overall system. The calculations are lengthy but computer programs for this purpose (2) have been published. In this paper simpler approximate methods for phase equilibrium problems of environmental interest is sought. For the moment it is sufficient to note that the activity coefficients provide the means of complete liquid-liquid equihbrium computations. 

The same vapor-liquid equilibrium ratio (AT) charts were used for the rigorous solution as for the shortcut, For the rigorous, however, values of K, were combined with total pressure P, and the Kf product was treated as effective vapor pressure" in the Antoine equation. The rigoroas program was nin in a mal-and-error fashion, with constant reflux ratio of 1.722 (from the shortcut) and with iterations of DIF ratio, total plates, and feed plate location. 

The dashed lines in the figure are the predictions, at all temperatures for the acetone-water system that result from setting the binary parameter k 2 equal to zero. Note that very nonideal behavior is predicted, which shows that setting ]c 2 = 0 is not equivalent to assuming ideal solution behavior. In fact, such extreme nonideal behavior is. predicted that vapor-liquid equilibrium calculations made with the program VLMU do not even converge for acetone mole fractions less than about 0.15. 

When the data are assembled, graphs are often prepared with curves positioned to provide a good representation, especially for experimental data with scatter. Alternatively, the coefficients of equations, theoretical or empirical, are computed using regression analysis programs. This is especially common for thermophysical property data, such as the vapor pressure, P as a function of the temperature, T, and vapor-liquid equilibrium data, as discussed later in this section. 

The current commercial process simulators include special programs for calculating vapor-liquid equilibrium under different conditions. 

The determination of the quality function represents the most time-consuming phase of the parameter set optimization. The high amount of computing-time results directly from the repeatedly calculated vapor-liquid equilibrium conditions. Nevertheless, these computations of the defined quality function F can be performed in parallel mode for all sets of parameters as the optimization algorithm only requires the quahty values calculated for the selections and fresh mutations. Consequently, the use of parallel processors is possible without a great deal of programming and only a very low level of data communications being required. 

The computer program ELBT of the CD version of this book permits the retrieval of data tables using any single or combined search criterion based on chemical system (a single component or the binary mixture), the vapor-liquid equilibrium property type, and/or the reference (authors, original source of data, and/or year range)... 

The ELBT Program on the CD includes a large, almost exhaustive, bibliographical database, EVLM 2006, giving 13476 references to experimental vapor-liquid equilibrium in mixtures and solution measurements for 20937 two- to nine-component systems. These components are organic or inorganic electrolytes and nonelectrolytes, ionic liquids, alloys, fused salts, polymers and other materials. 

In the GAMS program the vapor-liquid-equilibrium of equation (2) may be defined by ... 

Several of these programs will be briefly described here. VLEFIT is a fitting program which uses vapor-liquid equilibrium (VLE) data to determine the best values for the adjustable parameters... 

Figure 6. Vapor-liquid equilibrium data fitting program (VLEFIT)...

Vapor-Liquid Equilibrium Plus Chemical Reaction. To conclude this discussion of VLE calculations, I d like to mention that we have initiated a program to allow VLE plus chemical reaction in column calculations. A modification of Frank s dynamic simulation method ( ) is being used. The initial trials with complex VLE and kinetic models have been moderately successful, but further testing is required. Other simulation techniques are also being investigated. 

Summarizing, it takes a flexibile equation of state and the proper combination of thermodynamics, common sense, patience and programming finesse to utilize equations of state effectively in computerized vapor-liquid equilibrium and other design calculations. 

The general methods of design for multicomponent distillation apply, the principal difficulty being the paucity of vapor-liquid equilibrium data for these highly nonideal mixtures [23]. Computer programs are available [4]. 

The van Krevelen et al. (1949) study is considered to be a pioneering effort in Ibe correlation of vapor-liquid equilibrium data for systems of tbis type. It has saved as the basis for several sour water stripper design procedures, including the widely used approach described by Beychok (1967) and Wild s calculator program for sour water stiippa design (1979). 

---