Problem Statement The process design division of your company has asked you to mathematically model to describe vapor-liquid equilibria to enable them to design large scale separation columns for binary systems. In particular, your company wants to separate ethanol from water. They need a guideline for relating xi, yi, T and p necessary for process simulation. Alternatively, they wish to be able to set up T-x-y, p-x-y and x-y diagrams.

Procedure

The procedure we will follow to calculate and if possible can you upload excel sheet

1. We will first find experimental data for vapor liquid equilibria for ethanol-water system. Do the internet search and locate some experimental data T-x-y at a pressure of 1 atm. Since most distillation processes are performed at atmospheric pressure, experimental T-x-y data at 1 atm pressure will be ideal for our work.

2. We will use these data to calculate experimental activity coefficients y; and y; at different xs and ys. Remember Y1P(T.p,x1,X3) Y1 (TIPIXLXZ) = lelsat(T) Y2 D(T,p.%1,%3) TP xp,Xy) =22L727 V2 (T,P,x1,X2) 2P, S8(T) In the above equations, we assume that the vapor phase behaves as an ideal gas.

3. Use just determined y's to calculate experimental molal excess Gibbs Free Energy GE (T,P,x1,%2) =RT [X1 In y1 (T,p,x1,X2)] + %2 In y2 (T,P,x1,%2)]

4. In addition to tabulating the values of ys, x, y and T, you should also tabulate the values of GE/RTxix, and RTxix/GE.

5. Now you have experimental values of GE. You wish to fit G to three mathematical models. These are 1) one parameter Margules model ( GE/RTxixz = A), 2) two-parameter Margules model [GE/RTxyx; = A + B(x1 -x2)] and 3) two-parameter van Laar model [R'|'X1X2/GE = A +B'(x1x2)

6. Now perform the least square fits and determine the best values of A for one parameter Margules model, A and B for two parameter Margules model and A and B for van Laar model.

7. The above models for GE do not contain pressure and temperature explicitly. In general, the effect of pressure is negligible. However, temperature does have a significant impact. The parameters in the above models are indeed temperature dependent. Fortunately, the temperature effect on vapor liquid equilibria primarily comes through the temperature dependence of the pure component vapor pressures in the above equations. Typically, a temperature increase of 10C raises the vapor pressure by a factor of 2. However, the activity coefficients over the same temperature change by only a few percentage points which is within the experimental uncertainty. As a result, the temperature dependence of GE is neglected too. Please keep in mind that this will not apply to liquid-liquid equilibria as no vapor pressures would be involved there.

8. Now you are ready to predict vapor liquid equilibria for a wide range of conditions. Please remember, you cannot expect very good results for conditions that are far removed from the conditions of underlying experimental data. Also please remember, these models do not explicitly contain T and P.

9. The procedure is outline below. a. Choose a mathematical model for GE (T,p,x1,X2). b. Determine the parameters of that model from the limited experimental data as follows. c. Use modified Raoults law to determine the numerical values of y; and y.. d. Use Equation (1) to calculate numerical values of G at different experimental values of liquid compositions. e. Try and fit these values to three models: 1) Margules one constant model, 2) Margules two-constant model and 3) two-constant van Laar model. f. Determine the associated parameters from these linear fits or least squares analysis. g. Construct p-x-y and T-x-y diagrams.

Mathematical models for G a. One parameter Margules Model The model contains only one parameter, A. GE RTxle = Iny; = Ax; Iny, = Ax? b. Two parameter MargulesModel This model has two parameters A and B. GE =A+B(x x3) (1) RT%%, RTisz =A(x + x2) + B(x x2) (2) RTifxz =(@+Bx+(A-B)x; (3) RTisz = (A20)x; + (Ar2)x2 (4) Inyy = x%(A12 + 2(A21 A12)%1) (5) Iny, = xZ(Az1 + 2(A12 Az21)%2) (6) c. Two parameter van Laar Model The model has two parameters A and B'. B%% 4+ By x) @ "T Ao+ ) +B @ - x) ( M = (4 +B)x + (A BYx) (9) RTxlxz _ 1 TE T an M +ar (10) Iny, = 4, {1+ "zl} (11) 2 2 _ar A21%2 Iy, = A, {1+ A_ile} (12)

11. Fit the GE/RTx1x, data to all three models and determine the best values of the parameters for each case. Then you can use the corresponding equations to evaluate ys over the entire composition range. This will enable you to calculate T-x-y and P-xy diagrams over the full composition range.

12. Adopt the usual convention that component 1 is the more volatile component than 2. It implies that it has a higher vapor pressure at a given temperature than the component 2. You can calculate the vapor pressure of any chemical specie from Antoine equation B; Ci+T Here the values of A, B, C can be looked up from any reference book. Just be certain to use appropriate units as specified by the source.

13. In cases where the binary mixture forms an azeotrope, one may be able to use the azeotrope information (T, P, azeotropic composition) to ones advantage to calculate the model parameters. At azeotropic composition, X1 =y; and x; =y,. Hence P Y{(T, P, %1, %;) = psa 1 P YZaz(T' P,xl,xz) = @) Hence if you know the azeotropic composition and T and P of the azeotrope, one can easily calculate the activity coefficients of both components at the azeotropic conditions. One can then use these Ys to calculate the parameters. One parameter Margules Model In az A= 1'21 or X2 n az "=t x1 If the As turn out to be different, which one would you choose? Since the formation of an azeotrope is a highly non-ideal effect, and the activity coefficient is a measure of non-ideality, one may assume that the larger value of |Iny;|carries more significant information. Hence choose A coming from that component. Two parameter_Margules Model Here the model parameters can be easily computed from the azeotropic point information. A12=(2) 1\ Iny 2lnyf Y1 + Y2 X2/ X2 X1 + X1 X2 ( 1 ) Iny$% 2lny Once A1, and A; are known, you can use Equations 5 and 6 to calculate ys at any composition. Two parameter van Laar Model Here the model parameters can be easily computed from the azeotropic point information. 2 X, lnyfz} Ay =Iny 1+ 2 n { 1 2 X1 lnyl"} Ay 21 =1In Y2 z{l+ SIny Once these parameters are known, you can use Equations (11) and (12) to compute ys at any value of x.

14. Construction of T-x-y diagram at a fixed P and P-x-y diagram at a fixed T

15. Check for thermodynamic consistency of the experimental The measured experimental data must conform to data. The Gibbs-Duhem Equation has to be satisfied. Z x;[dIny; (T,P,xy,x;)] at constant T and P = 0 i Of course, for the binary system, the above equation reduces to X1 [d In y1 (T,P,x3,%2)] +x2 [d Iny2 (T,P,x1,x2)] =0 The activity coefficients y1 and y; in their common mixture are not completely independent of one another. You can integrate the above equation numerically to check if the sum of the two terms is equal to zero or close to zero.