IK CAPE Equations
IK CAPE Equations
IK CAPE Equations
7KHUPRG\QDPLFV
3DFNDJH
IRU
&$3($SSOLFDWLRQV
Model Description
(Translation by DECHEMA e.V.)
Table of contents
Page
INTRODUCTION
1-1
2-1
2.1
Equations
2-1
2.2
Extrapolations
2-3
CALCULATION OF AVERAGES
3-1
ACTIVITY COEFFICIENTS
4-1
4.1
NRTL model
4-1
4.2
UNIQUAC model
4-1
4.3
Wilson model
4-2
4.4
Flory-Huggins model
4-2
4.5
UNIFAC model
4-2
POYNTING CORRECTION
5-1
HENRY CORRECTION
6-1
EQUATIONS OF STATE
7-1
7.1
Redlich-Kwong-Soave equation
7-1
7.2
Peng-Robinson equation
7-2
FUGACITY COEFFICIENTS
8-1
8.1
Redlich-Kwong-Soave equation
8-1
8.2
Peng-Robinson equation
8-1
8.3
Virial equation
8-1
8.4
8-2
ENTHALPY
9-1
9.1
9-1
9.2
9-1
9.3
9-2
9.4
9-2
9.5
9-2
9.6
The excess enthalpy
9.6.1
NRTL equation
9.6.2
UNIQUAC equation
9.6.3
Wilson equation
9.6.4
Flory-Huggins equation
9.6.5
Redlich-Kister equation
9-3
9-3
9-4
9-4
9-5
9-5
9.7
The isothermal pressure dependency of the enthalpy in vapor phase
9.7.1
Redlich-Kwong-Soave equation
9.7.2
Peng-Robinson equation
9-5
9-5
9-5
9.8
9-6
10
CHEMICAL REACTIONS
10-1
10.1
Equilibrium reactions
10-1
10.2
10-2
10.3
10-3
1 Introduction
Inside of this handbook all models included inside the IK-CAPE Thermodynamics Package are documented.
Common abbreviations:
R
T
P
Tc
Pc
v
xi
yi
i
i
1-1
POLY : Polynom
f ( T ) = a0 + a1 T + a2 T 2 + ............. + a9 T 9
f (T ) = 10 a0 +a1T +a2T
+ ....+ a9T 9
f ( T ) = a0 ( a2 T ) a1 + a3
ln( f ( T )) = a0
ln f (T ) = a 0 +
a1
T + a2
a1
+ a3T + a 4 ln T + a5T a6
T + a2
ln( f ( T )) = a0
a1
+ a2 ln T
T
2-1
f (T ) =
a0 T
a
1+ 1
T
Tr =
1
( a2 + a3 1.5 + a 4 3 + a5 6 )
Tr
T
a0
= 1 Tr
a0 = Tc
a1 = Pc
f ( T ) = a0 + a1T + a2 T 2 + a3 T 3 +
a4
T2
f ( T ) = a0 e T + a 2
a5
T
2-2
a0
T
a11+(1 a2 )
a3
ln( f ( T )) = a0 +
a1
+ a2 ln T + a3 T a4
T
a2
a4
f (T ) = a 0 + a1 T + a3 T
sinh a 2
cosh a4
T
T
f (T ) = a1 (1 Tr ) h
with
h = a 2 + a3Tr + a 4Tr2 + a5Tr3
Tr =
T
a0
a 0 = Tc
f (T ) =
a 0 T a1
a
a
1 + 2 + 32
T T
2.2 Extrapolations
To avoid nonsensical results as well as program crashes, one can use the possibility to extrapolate the above given
functions beyond the proofed temperature range using one of the given extrapolation methods.
Linear extrapolation
For linear extrapolation inside the f(T)-diagram applies for the case
2-3
f (Tupper )
f (T ) = f (Tupper ) B exp
(T Tupper ) + B
f (Tupper ) B
f (Tupper )
f (T ) = f (Tupper ) exp
(T Tupper )
f (Tupper ) B
For
f (Tlower )
f (T ) = [ f (Tlower ) B ]exp
(T Tlower ) + B
f (Tlower ) B
f (Tlower )
f (T ) = f (Tlower ) exp
(T Tlower )
f (Tlower ) B
d ln f (T )
ln f (T ) =
(T Tupper ) + ln f (Tupper )
dT T =Tupper
For
f (Tupper )
(T Tupper ) + ln f (Tupper )
f (T ) = exp
f (Tupper )
f (T ) = f (T )
2-4
f (Tupper )
f (Tupper )
a 0 = Tc
a1 = Pc
Above the critical point the Wagner equation is extrapolated linear inside the ln(p) vs. 1/T diagram
1 1
a 2a 0
T a0
f ( T ) = a1 e
a a
f ( T) = 2 2 0 f ( T)
T
at T > Tc the vapor pressure represents a virtual value.
2-5
3 Calculation of averages
x value
average =
x
mass ,i
mass ,i
mass ,i
ln(valuei )
mass ,i
m = 0 . 5 x i i +
xi
x M
=
x M
i
3-1
m =
1
x
ii
i
average =
i
y i valuei
y j Fi , j
j
Fi , j =
valuei
1 +
value j
Mol j
Moli
Moli
81 +
Mol
j
average =
i
yi valuei
y j Fi, j
j
i Mol j
4
1 +
j Moli
Fi , j =
Moli
81 +
Mol
j
average = i
i xiVoli
3-2
average =
i
xvi =
2 xvi xv j
1
1
+
valuei value j
xiVol i
x jVol j
j
3-3
4 Activity coefficients
4.1 NRTL model
G x
xG
=
+
G x
G
j ,i
ln i
j ,i
i, j
k ,i
G j ,i = e
i, j
k , j xk
G x
G x
l, j
l, j
k, j
S j ,i j ,i
j ,i = a j ,i +
b j ,i
T
+ e j ,i ln T + f j ,i T
S j ,i = c j ,i + d j ,i (T 27315
. )
Vi =
ri x i
rj x j
j
Fi =
qi x i
and F i =
qj xj
j
q i xi
xj
ln i = ln iC + ln iR
with
ln iC = ln
Vi KZ
F
V
+
q i ln i + li i
xi
2
Vi
xi
x l
j j
and
ln
R
i
F j i , j
= q i 1 ln F j j ,i
j
j F k k , j
li = 5 (ri q i ) (ri 1)
4-1
j ,i = e
a j ,i +
b j ,i
T
+ c j ,i ln T + d j ,i T
x k k ,i
ln i = 1 ln x j i, j
j
k x jk, j
j
i, j = e
ai , j +
bi , j
T
+ ci , j ln T + d i , j T
1 1
+ (1 1 ) 2 i 1, 2 i
r r
2i
ln 2 k = 1 + ln
r2 k r2 k
+ 1r2 k 1, 2 k 2 i 1, 2i
r
r
2i
with
x1
r
r x
2i = 2 i 2 i
r
1 =
r = x1 + x 2i r2i
2i
1
T
Vi =
ri
rj x j
j
ri = (ki ) Rk
4-2
Fi =
qi
qj xj
j
qi = (ki ) Qk
For the activity coefficient applies:
ln i = ln iC + ln iR
with the combinatorial part
V
V
ln ic = 1 Vi + ln Vi 5qi 1 i + ln i
Fi
Fi
and the residual part, that describes the interactions between the different functional groups.
ln iR = k( i ) ( ln k ln k( i ) )
m km
ln k = Qk 1 ln m km
m
n n nm
m =
Qm X m
Qn X n
n
Xm =
( j)
m
xj
The parameter
( j)
n
xj
j ,i = e
a j ,i
T
4-3
5 Poynting correction
ln Fp ,i =
v iL ( P Pi s )
RT
5-1
6 Henry correction
ln H i =
x
j
ln Hi, j
whereas means:
i
Index of the Henry components
j
Index of the non Henry components
ln Hi , j = a +
b
+ c ln T + d T
T
6-1
7 Equations of state
7.1 Redlich-Kwong-Soave equation
The fugacity coefficients can be calculated using the following equation:
Bi
vm
a as B v + bm
z m 1) ln z m + ln
m 2 i i ln m
(
bm
v m bm bm RT a m bm
vm
ln i =
a m = x i x j ai, j
i
asi = x j ai, j
j
and
Bi = 2 bsi bm
bm = x i x j bi, j
i
bsi = x j bi, j
j
Herein means:
ai, j =
bi, j =
Ai A j 1 ki , j
bi + b j
(1 k )
bi , j
Ai ( T ) = ai i
i = 1 + mi 1 Tr
)]
mi = 0.48 + 1574
. i 0176
. i2
R 2Tc2
a i = 0.42748
Pc
7-1
bi = 0.0867
RTc
Pc
P=
RT
a(T )
v b v ( v + b) + b(v b)
Normal form of the cubic equation for the volume explicit representation:
Pb3 + RTb 2 ab
RT 2 3Pb 2 + 2 RTb a
+
v3 + b
v
v
=0
P
P
P
Calculation of compressibility:
zm =
Pvm
RT
zm=
vm
a m vm
vm bm RT vm ( vm + bm ) + bm ( vm bm )
am = xi x j ai , j
i
bm = xibi
i
Herein means:
ai , j =
Ai ( T ) A j ( T ) (1 ki , j )
Ai ( T ) = ai i
i = 1 + mi 1
Tci
mi = 0.37464 + 154226
.
i 0.26992 i 2
7-2
R 2 Tci 2
Pci
RTci
bi = 0.07780
Pci
ai = 0.45724
7-3
8 Fugacity coefficients
8.1 Redlich-Kwong-Soave equation
i =
i0
iv
ln i0 = z i 1 ln
Pi s
A (T ) v i + bi
vi bi ) i
ln
(
RT
bi RT v i
and
ln iv =
Bi
vm
a as B v + bm
z m 1) ln z m + ln
m 2 i i ln m
(
bm
v m bm bm RT a m bm
vm
The mixture rules can be find in the previous chapter Equations of state.
i =
i 0
i v
with
(
(
)
)
v i + 1 + 2 bi
PS
Ai (T )
ln i 0 = z i 1 ln i ( v i bi )
ln
RT
2 2bi RT v i + 1 2 bi
for a pure component at saturated vapor pressure
and
ln i v =
2
bi
am
P
( z m 1) ln
( v m bm )
bm
RT
2 2bm RT a m
x a
j
i, j
(
(
)
)
bi v m + 1 + 2 bm
ln
bm v m + 1 2 bm
i =
i0
iv
8-1
i0 =
2
Bi , i
v
1+
RT
v=
2 Pi 0
2
v
iv =
v=
Bi,i
v
4P0 B
1 + 1 + i i ,i
RT
y j Bi , j
Bm
v
1+
RT
4 PBm
1 + 1 +
2P
RT
Bm = yi y j Bi, j
i
zm = 1 +
Bm
v
with
Notation:
"true" concentration
z in
i,j ... components
n
association degree n=1,2,4,6 for i 5
zi2 = zi4 = zi6 = 0 for i > 5
i,j 5
Equilibrium constants:
For each association reaction
8-2
n (i) (i)n
K in =
z in
z p n 1
n
i1
K Mij =
z Mij
z i1 z j1 p
ln K in = Ain +
Bin
T
K Mij = 2 K i 2 K j 2
Set of equations to calculate the true monomer concentrations:
Material balance:
K
i
z p n 1 = 1
n
in i1
nK
z p n 1 +
n
in i1
nK
n
1n 11
z p
n 1
z z p
Mij i1 j1
j i, j 5
z z p
M1 j 11 j1
yi
y1
j i, j 5
i =
z i1
yi
8-3
9 Enthalpy
The following abbreviations are used:
h io
h il
h iv
ri
cp i
h v
dp
p
hl = h E + xi hil
i
with
T
h =h +
l
i
0
i
cl dT
i
T0 ,i
hl = h E + xi hil
i
If the temperature for the phase change is set equal to the systems temperature, the enthalpy of the pure component can
be calculated using:
hiv
dp ri (T )
h = h + cpi dT
p
T0 ,i
ps (T )
T
l
i
0
i
If the temperature for the phase change is predetermined, the enthalpy of the pure component has to be calculated using:
Tui
h = h + cpi dT
l
i
0
i
T0,i
pis ( Tui
T
hiv
dp ri (Tui ) + cli dT
p
Tui
)
9-1
h
hv = v dp + yi hiv
p
i
0
P
If the temperature for the phase change is set equal to the systems temperature, the enthalpy of the pure component can
be calculated using:
hiv
dp
h = h + cli dT + ri (T ) +
p
T0,i
ps (T )
T
v
i
0
i
If the temperature for the phase change is predetermined, the enthalpy of the pure component has to be calculated using:
Tui
h = h + cli dT + ri (Tui ) +
0
i
v
i
T0,i
pis ( Tui
T
hiv
dp + cpi dT
p
Tui
)
h
hv = v dp + yi hiv
p
i
0
P
with
T
h = h + cpi dT
v
i
0
i
T0, i
hs = si his
i
with
T
his = hi0 +
cs dT
i
T0,i
9-2
h E = T 2
gE
1 i
ln i
T
= RT 2 xi
= RT 2 xi
T
T
i T
i
i
x G
j
j ,i
j ,i
x G
j
j ,i
h E = T 2 R xi
i
Ai Bi Bi Ai
Bi2
with
Ai = x j G j ,i j ,i
j
Ai
= Ai = x j G j ,i j ,i + G j ,i j ,i
T
j
Bi = G j ,i x j
j
Bi
= Bi = x j G j ,i
T
j
G j ,i = e
S j ,i j ,i
S j ,i = c j ,i + d j ,i (T 27315
. )
j ,i = a j ,i +
b j ,i
T
+ e j ,i ln T + f j ,i T
G j ,i
j ,i
S j ,i
j ,i +
S j ,i
= G j ,i = G j ,i
T
T
T
j ,i
T
S j ,i
T
= j ,i =
b j ,i
T2
e j ,i
T
+ f j ,i
= d j ,i
9-3
A j i , j Aj
A
T
h E = RT 2 q i xi i + F j
A
Aj 2
i
j
with
Ai = F j j ,i
j
j ,i
Ai = F j
F i =
q i xi
q
j
xj
and
j ,i = e
a j ,i +
b j ,i
T
+ c j ,i ln T + d j ,i T
j ,i j j ,i j
h E = RT 2 xi i + x j
2j
i
i
j
i
= i = x j i , j
T
j
i = x j i , j
j
i , j
bi , j ci , j
= i , j = i , j 2 +
+ di, j
T
T
T
i, j = e
9-4
ai , j +
bi , j
T
+ ci , j ln T + d i , j T
gE
= x1 ln 1 + x 2i ln 2i + x1 2i 1,2i
RT
x1
x2i
2i
2i
h E = RTx1 2i 1,2i
2i
h E = 0,5 hiE, j
i
hiE, j =
xi x j
xi + x j
( A(T ) + B(T ) x
+ C (T ) x d2 + D(T ) x d3 + E (T ) x d4 + F (T ) x d5 )
x d = xi x j
The temperature dependency of the coefficients in most cases can be characterized well by a polynom:
A( T ) = a0 + a1T + a2 T 2 ...
hi = RT ( z 1) +
dai
ai b p s
dT
ln 1 + i i
bi
zRT
hm = RT ( z 1) +
dam
am
b p
dT
ln 1 + m
zRT
bm
9-5
hi = RT ( zi 1)
Ai (T ) vi + (1 + 2 )bi
1
Ai (T ) T
ln
T vi + (1 2 )bi
2 2bi
with
Ai (T ) = ai 1 + mi 1
Tci
Ai (T )
T ai mi
= 1 + mi 1
T
Tci Tci T
hm = RT ( zm 1)
a m v m + (1 + 2 )bm
1
am T
ln
T v m + (1 2 )bm
2 2bm
with
A j (T )
a m
Ai (T )
1
Ai (T )
= xi x j (1 k i , j )
+ A j (T )
T
T
T
2 Ai (T ) A j (T )
i
j
RE =
z
i
in
nz
i
in
+ 2 z Mij
hin = RBin
R
hMij = Bi 2 + B j 2
2
9-6
j i
j i
10 Chemical reactions
10.1 Equilibrium reactions
ln f ( T ) = a0 +
a1
+ a2 ln T + a3 T + a 4 T 2 + a5 T 3
T
EQLM
= f (T )
EQVM
= f (T )
EQLC
i
x
i Vi = f (T )
xi
i
V =
i
EQVC
i
y
i Vi = f (T )
V=
z RT
P
EQLA
( x )
= f (T )
EQVP
( Py )
= f (T )
EQLF
10-1
i0 Pi s i Fpi i
i f 0 xi = f (T )
EQVF
iv P i
i f 0 yi = f (T )
KILM
= V f k x i 1 + l i,l x i
k
l
i
i
i ,k
KIVM
= V f k y i 1 + l i,l y i
k
l
i
i
i ,k
KILC
fk
= V
k
Vs
x
Vs = i
i i
i,k
KIVC
fk
= V
k
Vs i ,k
Vs =
i , k
x
i i l 1 + Vsl i i,l xi
zRT
P
KILW
10-2
y
i
i , k
i
1 + l
i ,l y i
l Vs i
fk
= VmVs
k
Vs i , k
i , k
i
1 + l
i ,l x i
l Vs i
Vs = x i Mwi
i
KIVW
fk
= VmVs
k
Vs
i ,k
i , k
y
i i l 1 + Vsl i i,l yi
Vs = y i Mwi
i
m =
zRT
P
COOR
= value
CONV
STAT
10-3