Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

IK CAPE Equations

Download as pdf or txt
Download as pdf or txt
You are on page 1of 41



7KHUPRG\QDPLFV
3DFNDJH

IRU

&$3($SSOLFDWLRQV

Project of the IK-CAPE

Model Description
(Translation by DECHEMA e.V.)

Last Update: 28.03.02

Table of contents
Page

INTRODUCTION

2 DESCRIPTION OF THE TEMPERATURE DEPENDENCY FOR PURE


COMPONENT PROPERTY DATA

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

Association in vapor phase

8-2

ENTHALPY

9-1

9.1

Enthalpy of liquid with starting phase in liquid

9-1

9.2

Enthalpy of the liquid with starting phase in vapor

9-1

9.3

Enthalpy of a gas with starting phase in liquid

9-2

9.4

Enthalpy of a gas with starting phase in vapor

9-2

9.5

Enthalpy of a solid component

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

The enthalpy of a gas with simultaneous association

CHEMICAL REACTIONS

10-1

10.1

Equilibrium reactions

10-1

10.2

Kinetically controlled reactions

10-2

10.3

Reactions of general nature

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

gas constant in J/kmol/K


temperature in Kelvin
pressure in Pascal
critical temperature in Kelvin
critical pressure in Pascal
volume in m/kmol
mol fraction in liquid phase
mol fraction in vapor phase
activity coefficients
fugacity coefficients

1-1

2 Description of the temperature dependency for pure


component property data
2.1 Equations
To describe the temperature dependency of any pure component property at current stage the following functions can be
used:

POLY : Polynom

f ( T ) = a0 + a1 T + a2 T 2 + ............. + a9 T 9

EPOL : Polynom in exponent

f (T ) = 10 a0 +a1T +a2T

+ ....+ a9T 9

WATS : extended Watson equation

f ( T ) = a0 ( a2 T ) a1 + a3

ANTO : Antoine equation

ln( f ( T )) = a0

ANT1 : extended Antoine equation

ln f (T ) = a 0 +

a1
T + a2

a1
+ a3T + a 4 ln T + a5T a6
T + a2

KIRC : Kirchhoff equation

ln( f ( T )) = a0

a1
+ a2 ln T
T

2-1

SUTH : Sutherland equation

f (T ) =

a0 T
a
1+ 1
T

WAGN : Wagner equation


ln f ( T ) = ln a1 +

Tr =

1
( a2 + a3 1.5 + a 4 3 + a5 6 )
Tr

T
a0

= 1 Tr
a0 = Tc
a1 = Pc

CPL : Equation for the specific heat capacity of liquids

f ( T ) = a0 + a1T + a2 T 2 + a3 T 3 +

a4
T2

ICPL : Equation for the specific heat capacity of liquids

f (T ) = a0 + a1T + a2T 2 + a3T 3 + a4T 4 +

VISC : Equation for the dynamic viscosity


a1

f ( T ) = a0 e T + a 2

a5
T

RACK : Rackett equation


f (T ) =

2-2

a0
T

a11+(1 a2 )

a3

KIR1 : extended Kirchhoff equation

ln( f ( T )) = a0 +

a1
+ a2 ln T + a3 T a4
T

ALYL : Aly-Lee equation


2

a2
a4

f (T ) = a 0 + a1 T + a3 T
sinh a 2
cosh a4

T
T

DIP4 : DIPPR function for HVAP and ST

f (T ) = a1 (1 Tr ) h
with
h = a 2 + a3Tr + a 4Tr2 + a5Tr3
Tr =

T
a0

a 0 = Tc

DIP5 : DIPPR function for KVAP and VISV

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

T > Tupper and f (Tupper ) > 0

f (T ) = f (Tupper ) (T Tupper ) + f (Tupper )


f (T ) = f (Tupper )

2-3

and for the case

T < Tlower and f (Tlower ) < 0

f (T ) = f (Tlower ) (T Tlower ) + f (Tlower )


f (T ) = f (Tlower )
Extrapolation to the value B
For

T > Tupper and f (Tupper ) < 0 applies with 0 B < f (Tupper )

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

T < Tlower and f (Tlower ) > 0 gilt mit 0 B < f (Tlower )

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

linear extrapolation inside a logarithmic diagram

d ln f (T )
ln f (T ) =
(T Tupper ) + ln f (Tupper )
dT T =Tupper
For

T > Tupper and f (Tupper ) < 0 applies expressed in f(T)

f (Tupper )

(T Tupper ) + ln f (Tupper )
f (T ) = exp
f (Tupper )

f (T ) = f (T )

2-4

f (Tupper )
f (Tupper )

Extrapolation of the Wagner equation for vapor pressure

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

MOLA: average based on mole fractions


average = xi valuei
i

MASS: average based on mass fractions

x value
average =
x
mass ,i

mass ,i

MOLG: logarithmic average based on mole fractions


ln(average) = xi ln(valuei )
i

MALG: logarithmic average based on mass fractions


ln(average) =

mass ,i

ln(valuei )

mass ,i

LAMB: average for the heat conductivity of gaseous mixtures

m = 0 . 5 x i i +
xi

VISC: average for the viscosity of gaseous mixtures

x M
=
x M
i

3-1

VOLU: average for density based on volume

m =

1
x
ii
i

WILK: average from Wilke for viscosity of gaseous mixtures

average =
i

y i valuei
y j Fi , j
j

Fi , j =

valuei
1 +
value j

Mol j

Moli

Moli
81 +
Mol
j

WAMA: average from Wassiljewa,Mason,Saxena fr die heat conductivity of


gaseous mixtures

average =
i

yi valuei
y j Fi, j
j

i Mol j
4
1 +

j Moli

Fi , j =

Moli
81 +
Mol
j

= gas vis cos ity

DIST: average for the surface tension based on DIPPR


xiVol i 4 valuei

average = i
i xiVoli

3-2

DIKL: average for the heat conductivity of liquids based on DIPPR

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
. )

4.2 UNIQUAC model


The volume part is defined by:

Vi =

ri x i
rj x j
j

The surface part is defined by:

Fi =

qi x i
and F i =
qj xj
j

q i xi

xj

For the activity coefficient applies:

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

4.3 Wilson model

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

4.4 Flory-Huggins model


ln i = 1 + ln

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,2i = 10,2i + 11,2i


ri
i
1,2i

1
T

segment number of component i


volume fraction of component i
interaction parameter
solvent(1)-solved species(2i)

4.5 UNIFAC model


The volume part is defined by:

Vi =

ri
rj x j
j

ri = (ki ) Rk
4-2

The surface part is defined by:

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

Herein is m the area fraction of the functional group m

m =

Qm X m
Qn X n
n

with the mole fraction

Xm =

( j)
m

xj

The parameter

X m of the functional group m.

( j)
n

xj

ji contains the group interaction coefficient between the groups j and i:

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 =

using the following mixture rule:

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

The pure component parameters a, b can be received from:

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

7.2 Peng-Robinson equation


Pressure explicit representation:

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 )

using the following mixture rules:

am = xi x j ai , j
i

bm = xibi
i

Herein means:

ai , j =

Ai ( T ) A j ( T ) (1 ki , j )

The pure component parameters a, b can be received from:

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.

8.2 Peng-Robinson equation

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

for a component inside a mixture

8.3 Virial equation

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

Bi,i = Pure component virial coefficient


Bi, j = Cross virial coefficient

8.4 Association in vapor phase


Prerequisites:
The vapor phase behaves ideal, temperature, pressure and concentration are given
up to 5 components can build Di-, Tetra- and Hexamers as well as mixed dimers among each other
any desired number of inerts is allowed

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

Mixed dimers z Mij

i,j 5

Equilibrium constants:
For each association reaction

8-2

n (i) (i)n

(i) + (j) (ij)


equilibrium constants

K in =

z in
z p n 1
n
i1

K Mij =

are generated. For each applies:

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

Material rate for each component i:

nK

z p n 1 +

n
in i1

nK

n
1n 11

z p

n 1

Set of equations for the variable

z z p

Mij i1 j1

j i, j 5

z z p

M1 j 11 j1

yi
y1

j i, j 5

z in to be solved iterative. Afterwards all z in and z Mij can be calculated.

The fugacity coefficients can be calculated afterwards as follows:

i =

z i1
yi

8-3

9 Enthalpy
The following abbreviations are used:
h io
h il
h iv
ri
cp i

Reference enthalpy of the pure component


Enthalpy of the liquid, component i
Enthalpy of the vapor, component i
Enthalpy of vaporization, component i
specific heat capacity of the ideal gas, component i
ps (T )

h v

dp
p

9.1 Enthalpy of liquid with starting phase in liquid


The enthalpy of the liquid follows from:

hl = h E + xi hil
i

with
T

h =h +
l
i

0
i

cl dT
i

T0 ,i

9.2 Enthalpy of the liquid with starting phase in vapor


The enthalpy of the liquid follows from:

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

9.3 Enthalpy of a gas with starting phase in liquid


The enthalpy of the gas follows from:

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
)

9.4 Enthalpy of a gas with starting phase in vapor


The enthalpy of the gas follows from:

h
hv = v dp + yi hiv
p
i
0
P

with
T

h = h + cpi dT
v
i

0
i

T0, i

9.5 Enthalpy of a solid component


The enthalpy of a solid component follows from:

hs = si his
i

with
T

his = hi0 +

cs dT
i

T0,i

9-2

9.6 The excess enthalpy

h E = T 2

gE
1 i
ln i
T
= RT 2 xi
= RT 2 xi
T
T
i T
i
i

The excess enthalpy can be calculated using different methods.

9.6.1 NRTL equation


gE
= xi
RT
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

9.6.2 UNIQUAC equation


F
gE
= xi ln Vi + 5 q i xi ln i q i xi ln x j F j j ,i
RT
Vi
i
i
i
j
i , j

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

9.6.3 Wilson equation


gE
= R x i ln x j i, j
T
i
j


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

9.6.4 Flory-Huggins equation

gE
= x1 ln 1 + x 2i ln 2i + x1 2i 1,2i
RT
x1
x2i
2i
2i
h E = RTx1 2i 1,2i
2i

9.6.5 Redlich-Kister equation

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 ...

9.7 The isothermal pressure dependency of the enthalpy in vapor phase


9.7.1 Redlich-Kwong-Soave equation
The pressure dependency of the pure component is characterized as follows:

hi = RT ( z 1) +

dai
ai b p s
dT
ln 1 + i i
bi
zRT

The pressure dependency of the mixture:

hm = RT ( z 1) +

dam
am
b p

dT
ln 1 + m
zRT
bm

9.7.2 Peng-Robinson equation


The pressure dependency of the pure component is characterized as follows:

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

The pressure dependency of the mixture:

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

9.8 The enthalpy of a gas with simultaneous association


hv = RE + yi hiv
i

RE =

z
i

in

hin + z Mij hMij

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 )

10.2 Kinetically controlled reactions


a1
+ a 2 ln T + a3 T
T
a
ln ( T ) = a0 + 1 + a2 ln T + a3 T
T
ln f ( T ) = a 0 +

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

10.3 Reactions of general nature

COOR

= value

CONV

ratio * M ref = r ref

STAT

ratio * M ref = r ref (ratio 1)

10-3

You might also like