High Pressure VLE
High Pressure VLE
High Pressure VLE
yiˆiV P = xiˆiL P
y ˆV = x ˆ L
i i i i
yi / xi = ˆiL / ˆiV = Ki (T , P, yi , xi )
Proof: For a real mixture: dGi = RTd ln fˆi (1)
For a pure component: dGi = RTd ln f i (2)
From (1) and (2): Gi − Gi = RT ln( fˆi / fi ) (3)
Giid = Gi + RT ln yi Giid − Gi = RT ln yi (4)
From (3) and (4): Gi − Gi = RT ln( fˆi / yi f i )
id
fˆiV (T , P, yi ) = yi fiV (T , P)
ˆf V / y P = y f V / y P
i i i i i
ˆf V / y P = f V / P
i i i
Thus, ˆid =
i i
Contd…
With : ˆiid = i
y i ˆiL i (T , P ) fi L (T , P) / P
L
Ki = = V = V =
xi ˆi i (T , P ) iV (T , P )
Ki = yi / xi = f (T , P)
y = K x
i
i
i
i i =1
fˆiV = fˆi L
y ˆV P = x ˆ L P
i i i i
yiˆiV = xiˆiL
ˆ ˆ
yi / xi = i / i = Ki (T , P, yi , xi )
V L
Bubble Pressure Calculation at high pressures using Cubic EOS (Note that
this may be rewritten for Bubble T, and Dew T & P with suitable changes)
Read T, xi
physical constants, Print P, yi
estimates of P and yi
If ΣKixi > 1
then Yes
P (new)=P(old)/ ΣKixi
Estimate species No
fugacity coeff. If ΣKixi < 1 Is ΣKixi=1
for both phases and Ki then
P(new)=P (old)ΣKixi
Calculate ΣKixi No
{iV }