Must

For vapor-liquid equilibria, the equations of equilibrium which must be satisfied are of the form... 

In vapor-liquid equilibria, it is relatively easy to start the iteration because assumption of ideal behavior (Raoult s law) provides a reasonable zeroth approximation. By contrast, there is no obvious corresponding method to start the iteration calculation for liquid-liquid equilibria. Further, when two liquid phases are present, we must calculate for each component activity coefficients in two phases since these are often strongly nonlinear functions of compositions, liquid-liquid equilibrium calculations are highly sensitive to small changes in composition. In vapor-liquid equilibria at modest pressures, this sensitivity is lower because vapor-phase fugacity coefficients are usually close to unity and only weak functions of composition. For liquid-liquid equilibria, it is therefore more difficult to construct a numerical iteration procedure that converges both rapidly and consistently. 

The use of Henry s constant for a standard-state fugacity means that the standard-state fugacity for a noncondensable component depends not only on the temperature but also on the nature of the solvent. It is this feature of the unsymmetric convention which is its greatest disadvantage. As a result of this disadvantage special care must be exercised in the use of the unsymmetric convention for multicomponent solutions, as discussed in Chapter 4. 

The standard-state fugacity of any component must be evaluated at the same temperature as that of the solution, regardless of whether the symmetric or unsymmetric convention is used for activity-coefficient normalization. But what about the pressure At low pressures, the effect of pressure on the thermodynamic properties of condensed phases is negligible and under such con-... 

If we vary the composition of a liquid mixture over all possible composition values at constant temperature, the equilibrium pressure does not remain constant. Therefore, if integrated forms of the Gibbs-Duhem equation [Equation (16)] are used to correlate isothermal activity coefficient data, it is necessary that all activity coefficients be evaluated at the same pressure. Unfortunately, however, experimentally obtained isothermal activity coefficients are not all at the same pressure and therefore they must be corrected from the experimental total pressure P to the same (arbitrary) reference pressure designated P. This may be done by the rigorous thermodynamic relation at constant temperature and composition ... 

It is important to be consistent in the use of fugacity coefficients. When reducing experimental data to obtain activity coefficients, a particular method for calculating fugacity coefficients must be adopted. That same method must be employed when activity-coefficient correlations are used to generate vapor-liquid equilibria. 

Finally, values of n must be obtained when specific ("chemical") interactions can occur. These are difficult to estimate a priori but reasonable approximations can often be made by choosing a value (Appendix C) for a chemically similar system. 

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. 

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. 

For each experiment, the true values of the measured variables are related by one or more constraints. Because the number of data points exceeds the number of parameters to be estimated, all constraint equations are not exactly satisfied for all experimental measurements. Exact agreement between theory and experiment is not achieved due to random and systematic errors in the data and to "lack of fit" of the model to the data. Optimum parameters and true values corresponding to the experimental measurements must be found by satisfaction of an appropriate statistical criterion. 

If the m Equations (7-1) are considered independent, additional restrictions must be imposed as... 

Many different manipulations of these equations have been used to obtain solutions. As discussed by King (1971), many of the older approaches work in terms of V/L, which has the disadvantage of being unbounded and which, in the classical implementation, leads to poorly convergent iterative calculations. A preferable arrangement of this equation system for solution is based on the ratio V/F, which must lie between 0 and 1. If we substitute in Equation (7-1) for L from Equation (7-2) and for y from Equation (7-4), and then divide by F, we obtain... 

The equilibrium ratios are not fixed in a separation calculation and, even for an isothermal system, they are functions of the phase compositions. Further, the enthalpy balance. Equation (7-3), must be simultaneously satisfied and, unless specified, the flash temperature simultaneously determined. 

It is important to stress that unnecessary thermodynamic function evaluations must be avoided in equilibrium separation calculations. Thus, for example, in an adiabatic vapor-liquid flash, no attempt should be made iteratively to correct compositions (and K s) at current estimates of T and a before proceeding with the Newton-Raphson iteration. Similarly, in liquid-liquid separations, iterations on phase compositions at the current estimate of phase ratio (a)r or at some estimate of the conjugate phase composition, are almost always counterproductive. Each thermodynamic function evaluation (set of K ) should be used to improve estimates of all variables in the system. 

However, when carboxylic acids are present in a mixture, fugacity coefficients must be calculated using the chemical theory. Chemical theory leads to a fugacity coefficient dependent on true equilibrium concentrations, as shown by Equation (3-13). ... 

Appendix C presents the best set of constants for Equation (2). Also shown are the temperature limits of the real experimental data. Users must exercise caution when using the correlation outside the range of real data such use should, in general, be avoided. 

If only two parameters are estimated, C must be set to some value. 

If only two parameters are fit, C must set to some arbitrary value, usually one, and only Pqj and P 2) estimated from the VLE... 

If only two parameters are fit, a must be set to some predetermined value, and only and P 2) determined from the binary... 

NMD(1,J) cols 1-6 six-character identification code for component one this must match the code read in with the pure-component data. 

Multiple sets of binary VLE data may be correlated by continuing with another set of cards starting at part B. The last set of cards must be followed with a blank card to end the program. 

INITIAL CALL TO THIS SUBROUTINE MUST BE MADE WITH KEY.EQ.O... 

SHOULD BE AS IN VALIK. IF VALIK HAS NOT BEEN CALLED, KEY MUST BE 9. 

THE SUBROUTINE ACCEPTS BOTH A LIQUID FEED OF COMPOSITION XF AT TEMPERATURE TL(K) AND A VAPOR FEED OF COMPOSITION YF AT TVVAPOR FRACTION OF THE FEED BEING VF (MOL BASIS). FDR AN ISOTHERMAL FLASH THE TEMPERATURE T(K) MUST ALSO BE SUPPLIED. THE SUBROUTINE DETERMINES THE V/F RATIO A, THE LIQUID AND VAPOR PHASE COMPOSITIONS X ANO Y, AND FOR AN ADIABATIC FLASHf THE TEMPERATURE T(K). THE EQUILIBRIUM RATIOS K ARE ALSO PROVIDED. IT NORMALLY RETURNS ERF=0 BUT IF COMPONENT COMBINATIONS LACKING DATA ARE INVOLVED IT RETURNS ERF=lf ANO IF NO SOLUTION IS FOUND IT RETURNS ERF -2. FOR FLASH T.LT.TB OR T.GT.TD FLASH RETURNS ERF=3 OR 4 RESPECTIVELY, AND FOR BAD INPUT DATA IT RETURNS ERF=5. 

The temperature, T, and overall mole fractions, Z(I), of the system must be specified. 

PARIN first loads all pure component data by reading two records per component. The total number of components, M, in the library or data deck must be known beforehand. Next the associ-ation/solvation parameters are input for M components. Finally all the established UNIQUAC binary interaction parameters (or noncondensable-condensable interaction parameters) are read. 

Thus the complexity of chemical process synthesis is twofold. First, can we identify all possible structures Second, can we optimize each structure for a valid comparison When optimizing the structure, there may be many ways in which each individual task can be performed and many ways in which the individual tasks can be interconnected. This means that we must simulate and optimize... 

Most processes are catalyzed where catalysts for the reaction are known. The strategy will be to choose the catalyst, if one is to be used, and the ideal characteristics and operating conditions needed for the reaction system. Decisions must be made in terms of reactor... 

Before we can proceed with the choice of reactor and operating conditions, some general classifications must be made regarding the types of reaction systems likely to be encountered. We can classify reaction systems into five broad types ... 

The reaction is very rapid and must be stopped quickly to prevent excessive byproduct formation. 

Very often the choice is not available. For example, if reactor temperature is above the critical temperature of the chemical species, then the reactor must be gas phase. Even if the temperature can be lowered below critical, an extremely high pressure may be required to operate in the liquid phase. 

The choice of reactor temperature, pressure, arid hence phase must, in the first instance, take account of the desired equilibrium and selectivity effects. If there is still freedom to choose between gas and liquid phase, operation in the liquid phase is preferred. 

---