Problem 1) : CHEG 211 Chemical Engineering Thermodynamics Solution To Homework #3
Problem 1) : CHEG 211 Chemical Engineering Thermodynamics Solution To Homework #3
Solution to Homework #3
Problem 1)
1.1
B
ln
TC
T = 75oC = 348.15 K
P-x diagram:
. .
P x P x P x e .
. 1 x e
.
.
.
.
82.2717x 42.14241 x 42.1424 40.1293x
P-y diagram:
1 1
P
y 1 y y 1 y
. . 82.2717 42.1424
e .
.
. e .
.
.
1
0.02373 0.01157y
1.2
T-x diagram:
. .
70 x P x P x e .
. 1 x e .
.
at x1 = 0,
. .
70 0e .
. 1 0e .
.
at x1 = 1.0,
. .
70 1e .
. 1 1e .
.
T = 343.21 K = 70.06oC
T-y diagram:
Using 70
. .
. .
.
.
1.3
from (1.1),
1
P
0.02373 0.01157y
1
62.2
0.02373 0.01157y
y1 = 0.66
y 0.66
K 1.32
x 0.5
1.4
from (1.2),
1
70
y 1 y
. .
e .
. e .
.
T = 354.81 K = 81.66oC
. .
70 x e .
.. 1 x e .
.
.
x1 = 0.34
0.5
K 1.47
0.34
1.5
L+V=1
xiL + yiV = zi (N equations)
yiP = xiPi* (N equations)
y 1
z K
1
1 VK 1
where K
P = 58.88 kPa
1.6
xiL + yiV = zi
yiP = xiPi*
x1 = 0.417, y1 = 0.583
Problem 2)
2.1
P γ x P γ x P
2.2
1
P y y
γ P γ P
1 1
P
0.05 0.95 0.05 0.95
79.8γ 40.5γ
79.8e.
40.5e.
Nonlinear equations
f(x) = 79.8*x*exp(0.95*(1-x)^2)/0.05 - 1/(0.05/(79.8*exp(0.95*(1-
1
x)^2))+0.95/(40.5*exp(0.95*x^2))) = 0
1
P
0.05 0.95
79.8e.
. 40.5e.
.
P = 42.2 kPa
2.3
At the azeoptrope,
γ P γ P
x1= 0.857 = y1
x2 = 0.143 = y2
Problem 3)
3.1 V 120x 70x 15x 8x x x 120x 701 x $15x 81
x %x 1 x 70 58x x 7x
70 58x x 7x x 58 2x 21x 70 x 14x
!
3.2 V V x
substituting x1 = 0.5
V 72
dV
V V x 70 58x x 7x 1 x 58 2x 21x
dx
128 2x 20x 14x
By substituting x1 = 0.5
V 123.75
3.3
Show that
dV x dV x dV 0 is met.
dV 2 40x 42x dx
dV 2x 42x dx
x dV x dV x 2 40x 42x dx 1 x 2x 42x dx 0
3.4
By Gibbs-Duhem equation,
x dV x dV 0
Since x1 and x2 are positive, for a negative dV or a decrease in partial molar volume of
component 1, dV must be positive to meet the Gibbs-Duhem equation or the partial molar
volume of component (2) should increase.
Problem 4)
P, bar 0 10 20 40 60 80 100 200 300 400 500
Z 1.000 0.985 0.970 0.942 0.913 0.885 0.869 0.765 0.762 0.824 0.910
(Z-1)/P, -1.5* -1.5 -1.5 -1.45 -1.45 -1.435 -1.31 -1.175 -0.793 -0.44 -0.18
bar-1x103
*: (Z-1)/P at P=0 is indefinite since (1-1)/0. It was estimated by extrapolation.
Z 1
ln' ( dP
P
"
The integration is carried out by obtaining the area between the horizontal axis and
bounded by two pressures, one at P=0 and the other at the pressure at which fugacity
coefficient is to be determined.
Problem 5)
P
lnФ +B y δ /
RT
P
lnФ +B y δ /
RT
δ 2B B B 2x 466— 276— 809 153
2 bar
lnФ + 276 0. 5 x153/ 0.0164
cm • bar
83.1 x348.15 K
mol • K
Ф 0.984
2 bar
lnФ + 809 0. 5 x153/ 0.0533
cm • bar
83.1 x348.15 K
mol • K
Ф 0.948
Problem 6)
Problem 7)
7.1
Vapor pressure of methanol = system pressure at mole fraction of methanol = 1 or x1=1.
84.562 kPa
Vapor pressure of water = system pressure at mole fraction of methanol = 0 or x1=0.
19.953 kPa
7.2
P/kPa x1 y1 $ $
γ1(= ) γ2(= )
By linear interpolation:
γ1 = 1.162 + (1.097-1.162)/(0.5282-0.4461)x(0.5-0.4461) = 1.119
γ2 = 1.157 + (1.233-1.157)/(0.5282-0.4461)X(0.5-0.4461) = 1.207
7.3
G%
x lnγ x lnγ
RT
x1 γ lnγ x2 γ lnγ G%
RT
0.0000 1.0000 1.000 0.0000 0.0000
0.1686 1.572 0.4523 0.8314 1.013 0.0129 0.0870
0.2167 1.470 0.3853 0.7833 1.026 0.0257 0.1036
0.3039 1.320 0.2776 0.6961 1.075 0.0723 0.1347
0.3681 1.246 0.2199 0.6319 1.112 0.1062 0.1481
0.4461 1.162 0.1501 0.5539 1.157 0.1458 0.1477
0.5282 1.097 0.0926 0.4718 1.233 0.2095 0.1478
0.6044 1.050 0.0488 0.3956 1.311 0.2708 0.1366
0.6804 1.031 0.0305 0.3196 1.350 0.3000 0.1166
0.7255 1.021 0.0208 0.2745 1.382 0.3235 0.1039
0.7776 1.012 0.0119 0.2224 1.410 0.3436 0.0857
1.0000 1.000 0.0000 0.0000 0.0000
Problem 8)
Method A:
P γ x P γ x P
Determine the constants for the model to best fit measured system pressures
Pressure x1
39.223 0.1686
42.984 0.2167
48.852 0.3039
52.784 0.3681
56.652 0.4461
60.614 0.5282
63.998 0.6044
67.924 0.6804
70.229 0.7255
72.832 0.7776
Model: P = x*84.562*exp(A12*(1+A12*x/(A21*(1-x)))^(-2))+(1-x)*19.953*exp(A21*(1+A21*(1-
x)/(A12*x))^(-2))
By Polymath
A12 = 0.8301
A21 = 0.4677
Method B:
& &
Determine the constants for ln & predicted by the model to best fit ln & determined from
experimental data.
γ γ γ
x1 ln
γ
0.1686 1.572 1.013 -0.4394
0.2167 1.470 1.026 -0.3560
0.3039 1.320 1.075 -0.2053
0.3681 1.246 1.112 -0.1138
0.4461 1.162 1.157 -0.00431
0.5282 1.097 1.233 0.1169
0.6044 1.050 1.311 0.2220
0.6804 1.031 1.350 0.2696
0.7255 1.021 1.382 0.3027
0.7776 1.012 1.410 0.3317
Model: y = A21*(1+A21*(1-x)/(A12*x))^(-2)-A12*(1+A12*x/(A21*(1-x)))^(-2)
By polymath
A12 = 0.8283
A21 = 0.5227
Problem 9)
By Wilson model,
< <
lnγ ln; ; < x
; ; < ; ; <
< <
lnγ ln; ; < ;
; ; < ; ; <
γ ; ; < < <
ln ln
γ ; ; < ; ; < ; ; <
& &
Determine the constants for ln predicted by the model to best fit ln determined from
& &
experimental data.
γ γ γ
x1 ln
γ
0.1686 1.572 1.013 -0.4394
0.2167 1.470 1.026 -0.3560
0.3039 1.320 1.075 -0.2053
0.3681 1.246 1.112 -0.1138
0.4461 1.162 1.157 -0.00431
0.5282 1.097 1.233 0.1169
0.6044 1.050 1.311 0.2220
0.6804 1.031 1.350 0.2696
0.7255 1.021 1.382 0.3027
0.7776 1.012 1.410 0.3317
Model: y = ln((x+(1-x)*A12)/((1-x)+x*A21))-(A12/(x+(1-x)*A12)-A21/((1-x)+x*A21))
Precision
R^2 0.9970359
R^2adj 0.9966654
Rmsd 0.0045897
Variance 0.0002633
Comparison between experimental and predicted values in activity coefficient and system
pressure.
x1 pressure- pressure- γ -exp γ -cal (γ -exp - γ -exp γ -cal (γ -exp -
exp cal γ -cal)2 γ -cal)2
0.0000 19.953 1.000
0.1686 39.223 40.061 1.572 1.608 0.0013 1.013 1.033 0.0004
0.2167 42.984 44.002 1.470 1.491 0.0004 1.026 1.051 0.0006
0.3039 48.852 49.476 1.320 1.334 0.0002 1.075 1.094 0.0004
0.3681 52.784 53.156 1.246 1.250 0.0000 1.112 1.130 0.0003
0.4461 56.652 57.291 1.162 1.173 0.0001 1.157 1.180 0.0005
0.5282 60.614 61.421 1.097 1.114 0.0003 1.233 1.239 0.0000
0.6044 63.998 65.153 1.050 1.074 0.0006 1.311 1.300 0.0001
0.6804 67.924 68.830 1.031 1.045 0.0002 1.350 1.365 0.0002
0.7255 70.229 71.020 1.021 1.032 0.0001 1.382 1.407 0.0006
0.7776 72.832 73.540 1.012 1.020 0.0001 1.410 1.458 0.0023
1.0000 84.562 1.000
total 0.0033 total 0.0054
Problem 10)
Maingroup subgroup k Rk Qk νk
number
species (1) species (2)
1 CH3 1 0.9011 0.848 1 0
5 OH 15 1.0000 1.200 1 0
7 H2O 17 0.9200 1.400 0 1
lnγ lnγ lnγ
' (
The subscript C refers to combinatorial term to account for molecular size and shape,
and R is a residual property to account for molecular interaction.
J J
lnγ' 1 J lnJ 5q 1 ln
L L
#
where J ∑
#
q
L
∑* q * x *
r ν + R +
+
q ν+ Q +
+
ν+ is the number of subgroups of type k in a molecular species of k. Rk and Qk come from
β+ β+
lnγ( q +1 θ+ e+ ln /
s+ s+
+
ν Q+
e+ +
q
τ,+ exp where amk come from tabulations in the literature.
β+ e, τ,+
,
∑ x q e+
θ+
∑* x * q *
s+ θ, τ,+
,
+
r1 = 1 x 0.9011 + 1 x 1.000 = 1.9011
r2 = 1 x 0.9200 = 0.9200
(2) Calculation of Ji
#
J ∑
#
J1 = 1.9011/(0.5282 x 1.9011 + 0.4718 x 0.9200) = 1.3218
J2 = 0.9200/(0.5282 x 1.9011 + 0.4718 x 0.9200) = 0.6397
+
q1 = 1 x 0.848 + 1 x 1.200 = 2.048
q2 = 1 x 1.400 = 1.400
(4) Calculation of Li
-
L ∑
-
L1 = 2.048/(0.5282 x 2.048 + 0.4718 x 1.400) = 1.1755
L2 = 1.400/(0.5282 x 2.048 + 0.4718 x 1.400) = 0.8035
(5) Calculation of lnγ'
J J
lnγ' 1 J lnJ 5q 1 ln
L L
1.3218 1.3218
lnγ' 1 1.3218 ln1.3218 5 G 2.048 H 1 ln I
1.1755 1.1755
0.0305
0.6397 0.6397
lnγ' 1 0.6397 ln0.6397 5 G 1.400 H 1 ln I
0.8035 0.8035
0.08238
(7) Calculation of θ+
∑ x q e+
θ+
∑* x * q *
∑* x* q * 0.5282 x 2.048 0.4718 x 1.400 1.7423
0.5282 x 2.048 x 0.4141 0.4718 x 1.400 x 0
θ 0.2570
1.7423
0.5282 x 2.048 x 0.5859 0.4718 x 1.400 x 0
θ
0.3638
1.7423
0.5282 x 2.048 x 0 0.4718 x 1.400 x 1.000
θ 0.3791
1.7423
(8) Calculation of τ,+
a,+
τ,+ exp
T
From Table H.2
a1,1 = 0.00, τ, 1
a1,15 = 986.50, τ,
0.05176
a15,1 = 156.4, τ
, 0.6253
a15,15 = 0.00, τ
,
1
a1,17 = 1318, τ, 0.0191
a17,1 = 300.00, τ, 0.4064
a17,15 = -229.10, τ,
1.9891
a15,17 = 353.5, τ
, 0.3461
a17,17 = 0.00, τ, 1
(9) Calculation of s+
s+ θ, τ,+
,
s θ τ θ
τ
, θ τ,
0.2570 x 1 0.3638 x 0.6253 0.3791 x 0.4064 0.6386
s θ τ,
θ
τ
,
θ τ,
0.2570 x 0.05176 0.3638 x 1 0.3791 x 1.9891 1.1312
s θ τ, θ
τ
, θ τ,
0.2570 x 0.0191 0.3638 x 0.3461 0.3791 x 1 0.5099
Experimental UNIFAC
γ 1.097 1.148
γ 1.233 1.297