1. For the actone(1)/methanol(2)/water(3) system, based on Eq. (10.

5) and the Wilson

equiation make the following calculations:

(a) BUBL P . T = 338.15 K (65C), x1 = 0.3, x2 = 0.4.

(b) DEW P . T = 338.15 K (65C), y1 = 0.3, y2= 0.4.

(c) P . Tflash: T = 338.15 K (65C), P = (Pbubble + Pdew), z1 = 0.3, z2 = 0.4.

Work the preceding problem for the NRTL equation.

Peny :

Molar volume & Antoine coefficients:

Antoine coefficients:

74.05 14.3145 2756.22 228.060

V := (40.73) A := (16.5785) B := (3638.27) C := (239.500)
18.07 16.3872 3885.70 230.170

T:= (65 + 273.15)K Psat(i,T) := exp[ ]-kPa
( 273.15)+

NRTL parameters:

0.3084 0.5343 0 184.70 631.05

a =(0.3084 0 0.2994) b=( 222.64 0 253.88).
0.5343 0.2994 0 1197.41 845.21 0

i:= 1..3 j := 1..3

,
1 := 1..3 k := 1..3 i, j := . Gi, j := exp()-i,j+ i,j

(a) BUBL P calculation: Now iteration required

X1:= 0.3 X2 :,0.4 X3 :=1-X1-X2

(),,

1()1, 1
(i,x,T) := exp[ ]
, () , ,

+ [ ()1, 1[ , ()1, 1
1 1
 (,,)(,)
Pbubl := () (, , ) (, ) yi :=

0.525
y = ( 0.37 ) Pbubl = 115.3kPa Ans.
0.105

(b) DEW P calculation

y1:= 0.3 y2 := 0.4 y3 := 1 -y1 -y2

GUESS: x1 := 0.05 x2 := 0.2 x3 := 1 -x1 -x2 -x3 P := Pbubl

Given

Py1 = x1( 1,x,T)Psat(1,T) Py2 = x2(2,x,T)Psat(2,T)

Py3 = x3( 3,x,T)Psat(3,T) = 1

1
( 2 ) ()1, 2, 3,
3

0.038
X = (0.192) Pdew = 68.9kPa Ans.
0.77

+
(c). P.T -flash calculation: p := T = 338.15 K
2

z1 := 0.3 z2 := 0.4 z3 := 1 z1 z2

Guess: v := 0.5 Use x from DEW P and y from BUBL P as initial guess

Given py1(1,x,T)Psat(1,T) x1(1 V) + y1V = z1

Py2 = x2(2,x,T)Psat(2,T) X2(1 V) + y2V = z2

Py3 = x3(3,x,T)Psat(3,T) x3(1 V) + y3V = z3

= 1 = 1
 1
2
3
1 := Find()x1,x2,x3,y1,y2,y3, V
2
3
()

0.118 0.391
x = (0.347) y = (0.426) V = 0.667
0.534 0.183

2. For the acetone(1)/methanol(2)/water(3) system, based on Eq. (10.5) and the

Wilson equation make the following calculations:
(a) BUBL T: P = 101.33 kPa, x1 = 0.3, x2 = 0.4.
(b) DEW T: P = 101.33 kPa, y1 = 0.3, y2 = 0.4.
(c) P . T -flash: P = 101.33 kPa, T = (Tbubble + Tdew), z1 = 0.3, z2 = 0.2.
Peny :
Molar volumes & Antoine cofficients :

74.05 14.3145 2756.22 228.060

V:= (40.73) A:= (16.5785) B:= (3638.27) C:= (239.500)
18.07 16.3872 3885.70 230.170

Psat(i,T) := exp[ ]-kPa P := 101.33 kPa
( 273.15)+

161.88 291.27

Wilson parameters: a :=( 583.11 0 107.08).
1448.01 469.55 0

Vj
A(i,j,T) := Vi . exp i := 1..3 j := 1..3 p := 1..3

(a). BUBL T calculation:

X1:= 0.3 X2 :=0.4 X3 :=1 -X1 -X2

 (I,x,T) :=exp1 [() (, , )]

(,,)
[+ ()(,,)]

Guess: T := 300 K y1 := 0.3 y2 := 0.3 y3 := 1-y1-y2

Given

Py1 = x1 ( 1,x,T)Psat(1,T) Py2 = x2(2,x,T)Psat(2,T)

Py3 = x3 ( 3,x,T)Psat(3,T) P = . (i,x,T) Psat(i,T)

1
2
( ) ()1, 2, 3,
3

0.536
y = (0.361) Tbulb = 334.08 K Ans.
0.102

(b). DEW T calculation :

y1:= 0.3 y2 := 0.4 y3 := 1 -y1 -y2

GUESS: x1 := 0.05 x2 := 0.2 x3 := 1 -x1 -x2 T := Tbubl

Given

Py1 = x1( 1,x,T)Psat(1,T) Py2 = x2(2,x,T)Psat(2,T)

Py3 = x3( 3,x,T)Psat(3,T) = 1

1
( 2 ) ()1, 2, 3,
3

 0.043
X = (0.204) Tdew = 347.4K Ans.
0.753
+
(c). P.T -flash calculation: p := T = 340.75K
2

z1 := 0.3 z2 := 0.4 z3 := 1 z1 z2

Guess: v := 0.5 Use x from DEW P and y from BUBL P as initial guess

Given py1(1,x,T)Psat(1,T) x1(1 V) + y1V = z1

Py2 = x2(2,x,T)Psat(2,T) x2(1 V) + y2V = z2

Py3 = x3(3,x,T)Psat(3,T) x3(1 V) + y3V = z3

= 1 = 1

1
2
3
1 := Find()x1,x2,x3,y1,y2,y3, V
2
3
()

0.125 0.536
x = ( 0.17 ) y = (0.241) V = 0.426 Ans.
0.705 0.223

3. At 298.15 K (25C) and atmospheric pressure the volume change of mixing of

binary liquip mixtures of species 1 and 2 is given by the equation:
V = x1x2(45x1 + 25x2)
Where V is in cm3 mol1. At these conditions, V1 = 110 and V2 = 90 cm mol1.
Determine the partial molar volumes 1 and 2 in a mixture containing 40 m0l-%
of species l at the given conditions.
Peny :

x1 := 0.4 x2 := 1 -x1 V1 := 110 V2 := 90

VE()x1,x2 := x1,x2 := x1x2()45x1 + 25x2 VE()x1,x2 = 7.92

By Eq. (12.27):V()x1 , x2 := VE()x1,x2 + x1V1 + x2V2

V()x1,x2 = 105.92

By Eqs. (11.15) & (11.16):

3
Vbar1 := V()x1,x2 + x21 ( )2 Vbar1 = 190.28

Ans.

Vbar2 := V()x1,x2 x1(1 ()1, 2) Vbar2 = 48.68

Check by Eq. (11.11):

V := x1 Vbar1 + x2Vbar2 V = 105.92