Nonideal Behavior
Nonideal Behavior
Nonideal Behavior
Nonideal Behavior
P-x-y
1 atm
t-x-y diagram
Effect of Nonideality:
Chloroform(1)/Tetrahydrofuran(2)
30o
C
P-x-y
1 atm
t-x-y diagram
Effect of Nonideality:
Furan(1)/Carbontetrachloride(2)
30o
C
P-x-y
1 atm
t-x-y diagram
Effect of Nonideality:
Ethanol(1)/Toluene(2)
65o
C
P-x-y
1 atm
t-x-y diagram
Effect of Nonideality:
x y Diagram at Constant P = 1 atm
a. Tetrahydrofuran(1)/Carbon-tetrachloride
Chloroform(1)/Tetrahydrofuran(2
c. Furan(1)/Carbontetrachloride(2)
d. Ethanol(1)/Toluene(2)
M i U i , H i , Si , Vi , or G i , etc.
x i Mi
M x i Mi
M
NOT
:
M i U i , H i , Si , Vi , or G i , etc.
Partial
Properties:
.are
properties of component i in the state of
mixtures, which, in general different from
that in the state of pure species
x i Mi
nM M (T, P, n1 , n 2 , n 3 ,..........)
nM
d nM
nM
n
1
nM
dT
P,n
T , P, n 2 , n 3
nM
dn1
n
2
,...
nM
nM
d nM
dT
dP
T P,n
P T,n
NOT
x i Mi
dP
T,n
dn 2 .....
T ,P, n1 , n 3 ,...
nM
Mi
n
i
T , P,n j
nM
dn i .
n i T,P,n j
x i Mi
Tg T
g
P P
i g i
(nG )
i
n i T,P,n j
nG
nG
d nG
dT
dP
T P,n
P T,n
nG
dn i .
n i T,P,n j
nG
n
i
Gi
T , P, n j
(nG )
n
i
i Gi
T ,P, n j
x i Mi
M1 M x 2
M
x1
M 2 M x1
M
x1
Exercise
A group of students came across an unsuspected supply of laboratory
alcohol, containing 96 mass-percent ethanol and 4 mass-percent water.
As an experiment they decided to convert 2 L of this material into vodka,
having a composition of 56 mass-percent ethanol and 44 mass-percent
water. Wishing to perform the experiment carefully, they search the
literature and found the following partial-specific volume data for
ethanol water mixtures at 25 oC and 101.3 kPa.
In 96% ethanol
In vodka
V H 2O , L/kg
0.816
0.953
V EtOH , L/kg
1.273
1.243
Fugacity, f
dG S dT V dP
At constant
T
Ideal gas :
dG ig R T d ln P
Real gas
dG R T d ln f
dG
R T d ln
Residual Gibbs
energy :
G R G G ig
Fugacity
coefficient :
f
P
GR
ln
R T
Residual
RProperty
V V V ig
V R ( Z 1)
RT
P
dG i R Vi R dP
Vi R ( Zi 1)
RT
P
dG i R
RT
Vi R
dP
RT
GiR
dP
( Zi 1)
RT
P
GiR
ln i
R T
ln i ( Zi 1)
0
dP
P
P,1a
t00m1.0
Z69 P,6a
t00m1.4
Z31
2
0
0
1
.
1
3
8
7
0
0
1
.
5
0
4
3
0
0
1
.
2
0
9
8
0
0
1
.
5
7
7
4
0
0
1
.
2
8
3
9
0
0
1
.
6
4
9
500 1.356 1000 1.720
B 0.083
B1 0.139
0.422
Tr1.6
0.172
Tr 4.2
Bi P
ln i
RT
Pr 0
ln i
B B1
Tr
Tr
T
Tc
Pr
P
Pc
bi P
ZR T
ai
1
h
Z
1 h b i R T1.5 1 h
0.08664 R Tci
Pci
bi
to be solved
simultaneously
ln Z 1 ln (1 h ) Z
1.5
bR T
i i
i f Pr , Tr
0
P
Pr
Pc
i1
i1 f Pr , Tr
T
Tr
Tc
TPc
Evaluation
of Pure
Component
Fugacity, fi :
Pitzers
Generalized
Correlation
Tr
T
Tc
Evaluation
of Pure
Component
Fugacity, fi :
Pitzers
Generalized
Correlation
i 0
Pr
P
Pc
TPc
Evaluation
of Pure
Component
Fugacity, fi :
Pitzers
Generalized
Correlation
Tr
i1
T
Tc
Evaluation
of Pure
Component
Fugacity, fi :
Pitzers
Generalized
Correlation
Pr
P
Pc
Evaluation of Liquid
Pure Component Fugacity, fi
fi
1
ln
sat
R T
fi
Vi dP
Pi sat
sat
Vi P Pi
sat
sat
f i i Pi exp
R T
Poynting
Fugasity of saturated vapor, factor
calculated exactly as calculating gas
Examples of
Evaluation of Liquid Pure Component
Fugacity, fi
11.5
Estimate the fugacity of liquid acetone at
110oC and 275 bar.
At 110oC the vapor pressure of acetone is
4.36 bar and the molar volume of
saturated-liquid acetone is 73 cm3.mol-1
11.6
Estimate the fugacity of liquid n-butane at
120oC and 34 bar.
At 120oC the vapor pressure of n-butane is
22.38 bar and the molar volume of
saturated-liquid n-butane is 137 cm3.mol-1
Examples of
Evaluation of Liquid Pure Component
Fugacity, fi
11.10
The normal boiling point of nbutane is 0.5oC.
Estimate the fugacity of liquid nbutane at this temperature and
11.11
200 bar.
The normal boiling point of 1pentene is 30.0oC.
Estimate the fugacity of liquid 1pentene at this temperature and
11.12
350 bar.
The normal boiling point of isobutane
is -11.8oC.
Estimate the fugacity of liquid
Examples of
Evaluation of Gas & Liquid Pure Component Fugacity, fi
11.13
Prepare plots of f vs P and vs P for isopropanol
at 200oC for the pressure range from 0 to 50
bar. For the vapor phase, values of Z are given
by:
Z 1 9.86 10 3 P 11.41 10 5 P 2
Where P is in bars. The vapor pressures of
isopropanol at 200oC is 31.92 bar, and the
liquid-phase isothermal compressibility at
200oC is 0.3.10-3 bar-1, independent of P.
Hint:
1 V
V P
Examples of
Evaluation of Gas & Liquid Pure Component Fugacity, fi
11.14
Prepare plots of f vs P and vs P for 1,3butadiene at 40oC for the pressure range from
0 to 10 bar. At 40oC The vapor pressures of
1,3-butadiene is 4.287 bar.
Assume virial model to be valid for the vapor
phase.
The molar volume of saturated liquid 1,3butadiene at 40oC is 90.45 cm3.mol-1
P*
(Si Si )
P* : lowest value of P in
steam table
At P >= Pisat, i.e. liquid phase
water:
V P P sat
sat
sat
i
f i i Pi exp i
R T
Fugacity of Mixtures
Virial
:B P
Are formulated exactly as calculation for pure
ln i
yi y j Bij
i
Z cij
R Tcij
Pcij
Z ci Z cj
2
B0 ij B1
1
3
ij
1
3
Vcij
V V
ci
cj
i j
Z cij R Tcij
Vcij
(1 k ij )
RT
Fugacity of Components in
Mixture
f i
i
yi P
GiR
ln i
R T
Thus:
Gi R
ln i
R T
n ln
n
i
T, P, n j
ln i
Virial, binary
mixtures:
ln(i )
P
(B11 y 2 2 12 )
RT
P
ln 2
( B 22 y12 12 )
RT
ln 1
11.18
Estimate the fugacity and fugacity coefficient of
ethylene (1) and propylene (2) for a binary mixture
of 25% ethylene as a gas at 200oC and 20 bar.
Set all kij = 0
Mixing
B Rules
yi :y j Bij
i
Bij B ji
P
1
B kk
ln k
R T
2
i
R Tcij
Pcij
B 0 ij B1
yi y (2 ik i)
ik 2 B ik Bii B kk
i 2 B i Bii B
ii 0
kk 0
ki ik
yi y j a ij
i
bP
ZR T
1
a
h
Z
1 h b R T1.5 1 h
yi bi
i
to be solved
simultaneously
ln Z 1 ln (1 h ) Z
1.5
bR T
1.5
bR T
b1
b1
a
ln 1 ( Z 1) ln Z (1 h )
1
.
5
b
b
b RT
ln 2
ln Z 1 ln (1 h ) Z
b2
b
a
( Z 1) ln Z (1 h )
2
b
b RT 1.5 b
x k a1k
ln(1 h )
x k a 2k
ln(1 h )
UTS 1
G E G G id
Excess
Property
E
id
Partial
Property of the Excess
Gi Gi Gi
Property
G E G G id
Excess
Property
E
id
Partial
Property of the Excess
Gi Gi Gi
Property
Definition of
fugacity:
dG R T d ln f
Activity Coefficient
Integrati
fi on
G i G i R T ln
fi
G i id G i R T ln( x i )
G i G i id R T ln
n G E / RT
n i
fi
x i fi
(Ideal
solution)
fi
Gi E
ln
x i fi
R T
ln i
T , P, n j
GE
R T
x i ln i
Margules(1856 1920)
GE
A 21 x1 A12 x 2
x1 x 2 RT
n G E / RT
ln i
n i
ln 1 x 2 2 A12 2 (A 21 A12 ) x1
ln 2 x12 A 21 2 (A12 A 21 ) x 2
T , P, n j
van Laar
A12 ' A 21'
GE
n G E / RT
ln i
n i
'
ln 1 A12
1
ln 2 A '21 1
'
A12
x1
'
A 21 x 2
A '21 x 2
'
A12
x1
T ,P,n j
Wilson
V2
a12
12
exp
V1
RT
21
V1
a
exp 21
V2
RT
Wilson
GE
RT
x i ln x j ij
i
ln i 1 ln
x j ij
x k ki
x j kj
k
a ij
ij
exp
Vi
RT
ij 1
Vj
(i j)
(i j)
Models for
Multicomponent
Mixtures Activity
Coefficient:
UNIversal
QUAsi
Chemical
(UNIQUAC)
(Abrams & Prausnitz)
ln i ln i C ln i R
UNIquac
Functional
group
Activity
Coefficient
(UNIFAC)
(Aa Fredenslund,
Rl Jones & JM
Prausnitz)
UNIFAC: Rk & Qk
Example