Universidad Nacional de Ingeniería

Facultad de Ingeniería Química a) Componentes puros, para un gas ideal

Area Académica de Ingeniería Química
Termodinámica-II (PI-217)
 G ig 
V ig    dGigi =RTd lnP (T = cte )
 P T , x

Fugacidad y igi  P,T   iO  Po,T   RT  ln  P 

igi  io
Coeficiente Presión es una medida
de Fugacidad del potencial químico

Profesor del curso: (o) Estado estándar

Ing. Ubaldo Apaza Po = 1 bar, T

Lima 26/06/20
2
1

FUGACIDAD : medida de la tendencia al escape Gi  i (T )  RT ln  f i 

f 
Gi  Giig  RT ln  i 
dGi =RT  d lnfi  (T = cte ) P
Giig  i (T )  RT ln P
Gi  P, T   GiO  Po, T   RT  ln  fi 
Energía de GiR  RT ln i Coefficiente
i 
fi
Gi  P, T   i T   RT  ln  fi  Gibbs Residual de Fugacidad
P

Tiene unidades de presión

• Potencial Químico: μi P dP
– Provee criterio fundamental para el equilibrio de fases. ln i   ( Zi  1) (T  constante)
– Sin embargo, la energía Gibbs ( μi) se define con relación a la
0 P
energía interna y la entropía (valores absolutos desconocidos).
 Z 1
v
ln i   Z i  1  ln  Z i     i
v 
• Fugacidad : fi dv

- Una cantidad que toma el lugar de μi 4
3
1
 Ecuaciones Viriales Ecuaciones Cúbicas
EOS M. ABBOTT EOS Orbey-Vera Ecuación de Soave

B C Z 3  Z 2   A  B  B2  Z  A  B  0
Z  1
BP BPc
 B o  wB1 Z  1 
RTc v v2
RT
A  B
0.422 CPc2 ln    Z  1  ln  Z  B   ln  1  
B o  0.083  B  0.0193 
1 0.172  C o  wC 1 B  Z
Tr1.6 Tr 4.2 R 2Tc2
0.2432 0.0013
C o  0.01407  
Tr Tr 10.5 Ecuación de Peng Robinson

PB C1  0.02676 
0.05539 0.00242
 Z 3   B  1 Z 2   A  2B  3B 2  Z  A  B  B 2  B 3  0
ln    Tr 2.7 Tr 10.5

RT
B C ln    Z  1  ln  Z  B  
A 

 Z  1 2 B 


ln    Z  1  ln  Z    2
ln
v 2v  
21.5 B  Z  1  2 B 
 
5 6

VLE para especies puras

• Vapor Saturado: Giv  i (T )  RT ln f i v
Fugacidad para Líquido Puro
• La fugacidad de una especie pura “i” como líquido comprimido:
• Líquido Saturado: Gil  i (T )  RT ln f i l
fi v
Giv  Gil  RT ln fi P
fil Gi  Gisat  RT ln sat
Gi  Gisat   viL dP (T  cte)
fi Pi sat
VLE
v
fi
Giv  Gil  RT ln 0 fi 1 P
 
Vi no es función de P
fil ln sat sat
viL dP
fi RT Pi

f i v  f i l  f i sat iv  il  isat vil ( P  Pi sat )

fi v l ( P  Pi sat ) f i   P exp
sat sat
Para las especies puras donde coexistien las fases liquido y vapor, ln  i i i
f i sat RT RT
están en equilibrio cuando ellos tienen la misma temperatura, presion,
f i sat  isat Pi sat 8
fugacidad y coeficiente de fugacidad.
7
2
 Solución Analítica CEOS p2 2 p3 pq
Z 3  pZ 2  qZ  r  0 mq n  r
Tc 3 27 3
• Estado de fi Liquido vil ( P  Pi sat )
Saturación ln  D
n 2 m3
 M  3 
n
 D n
P fi sat RT 4 27 2 N3  D
• Líquido 2
Comprimido fi liq = fi vap  il  sat  iv  sat Caso I : D > 0. Una raíz real y dos raíces imaginarias conjugadas
M N p 3 i  ( 1)
a
2

3
b M  N 
L T1<Tc 2
CEOS
Z1  M  N 
p Z3  a  b * i Z2  a  b * i
Zv: Raíz mayor 3
vapor saturado
V Caso II : D = 0. Tres raíces reales, donde dos raíces son iguales
ZL: Raíz menor
líquido saturado T2<Tc  n p n p
Z1  2  3    Z 2  Z3  3  
 2 3 2 3
L-V
Z : Raíz no Caso III : D < 0. Tres raíces reales y diferentes.
significativa Z 3  Z 2   A  B  B2  Z  A  B  0  m  2  p k = 0, 1, 2
Z K  2    cos   k 
 3 3 3  3
v
 Z 1 A  B (-) cuando n>0
EOS RKS ln i   Zi  1  ln  Zi     i  dv ln    Z  1  ln  Z  B   ln  1   n2
B  Z cos     4 Se usa el signo : (+) cuando n< 0

v  3
9 m
27 10

Problema ._ Estime el valor del coeficiente de fugacidad para vapor 1-buteno a 200°C y 70 bar.
1-buteno T, °K 473.15
Método de Van Ness w 0.191 P, bar 70.00
PB
ln   
Tc, °K 420.00 Pr 1.7314
Pc, bar 40.43 Tr 1.1265
RT a(T ) RT
P 
v  b  v   b  v   b  EOS M. ABBOTT EOS SOAVE
B (cc/gmol) -223.795 w 0.90720
i 0.67150 A 0.52906
Fase líquida , valor inicial : Z = B B 0.13316
EOS ORBEY – VERA
1 B  Z
Z  B   Z   B  Z   B  C0 0.03476 Zinicial 0.72844
A C1 0.01270 Zfinal 0.54212
Fase vapor , valor inicial : Z = 1 C(cc/gmol)2 27740.449 Convergencia -5.7918E-07
V (cc/gmol) 321.645
Z B
Z 1 B  A Convergencia 4.359E-07 Z 3  Z 2   A  B  B2  Z  A  B  0
 Z   B  Z   B  Z 0.57236
i 0.64962 i 0.64635
Z 3       B  B  1 Z 2   A       B        B 2  Z B C A  B
ln    Z  1  ln  Z    2 ln( )  Z  1  ln  Z  B   ln 1  
v 2v B  Z
 A  B   ( B  B )  0
2 3

11 12
3
 1-buteno
0.10
EOS- SRK PTC Mathcad Prime 3 [2013] 1-Buteno
f(Z)
0.08

0.06 Z B
Z  1 B  A
Z Z  B Zinicial = 0.72844
0.04

0.02 Z 3  Z 2   A  B  B 2  Z  AB  0
Z
0.00
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 EOS-PR
-0.02
Z = 0.54212
-0.04

-0.06
13
HYSYS [MBWR] Z=0.747603
14

Correlaciones Generalizadas para el cálculo

del Coeficiente de Fugacidad de un gas puro
Tr 1.127
dP
ln i    Z  Z  Z
Pr
( Zi  1) r (Tr  cte) 0 1
0 Pr Pr 1.731

dPr dP
ln     ( Z 0  1)
Pr Pr
  Z 1 r (Tr  const.)
0 Pr 0 Pr
 0  0.627
o
Pr dPr Pr dPr
ln  0   ( Z 0  1) ln  1   Z 1 Fig 11.3
ln   ln  0   ln  1 con 0 Pr 0 Pr Van Ness
4ta edición
Tabla E1:E4 o Tabla E13:E16
15
164
 1
Tr

 1  1.096

Fig 11.5 Tr 1.127

Van Ness
4ta edición  0  0.627
Pr Pr 1.731
18
17

EOS Lee Kesler

c
1
B

C

D

4 
  
  exp   

2  2
Z
Vr 2 5 3 2
Vr Vr Tr  Vr  Vr   Vr 

 
ln i Z  1  ln( Z) 
Vr
B

C
2

D
5
E
2 Vr 5 Vr

b b b c c
2 3 4 2 3
B b    C c  
 1  1.096 1 Tr
Tr
2
Tr
3 1 Tr 3
Tr

d c
2 4 
   1     1 
   exp    

d 
2  2 
D E
1
  0.638
w
   0  1  f   P  (0.638)(70)  44.7 bar Tr 2 Tr  
3
  Vr   Vr  
20 5
19
 EOS Lee Kesler EOS Lee Kesler
1-Buteno Pr  1.7314 Tr  1.1265
Fluido Fluido Fluido Fluido
constante constante
simple referencia simple referencia Fluido simple Z (o)  0.46444 ln( (o) )  0.47420
b1 0.1181193 0.2026579 c3 0 0.0169010
b2 0.2657280 0.3315110 c4 0.042724 0.0415770
b3 0.1547900 0.0276550 d1x10 4
0.155488 0.4873600
Fluido referencia Z ( r )  0.53295 ln( (r) )  0.43670
b4 0.0303230 0.2034880 d2x104 0.623689 0.0740336
 0  0.62238 1  1.09885
c1 0.0236744 0.0313385  0.653920 1.2260000
c2 0.0186984 0.0503618  0.060167 0.0375400 wi 
ln i   ln    ln    ln   
(o) (r) (o)
(r)  
w
Fluido referencia : n-octano
21   0.633689 f    P  44.3582 bar
22

Tc, °K Pc, bar w Tr Pr  Z  B 

1-Buteno n-C4 425,12 37,96 0,2 1,176 1,581 ln( )  Z  1  ln  Z  B  
A
ln 
21.5 B  Z   B 


T(K) 500 M 0,6723

P(bar) 60 a 0,8896 0.10
i fi bar A
B
0,4648
0,1046 0.08
0,63501753
Z 1 f(z) 0.06
M ABBOTT 0.67150 47.005 Valor inicial, Z
Valor final , Z
0,9251
0,63501753
EXCEL
0.04
Función converg -2,7005E-09
Orbey Vera 0.64962 45.474 Coef. de Fugacidad 0,6921 0.02
0.00
LK grafico 0.63800 44.660 3 2 2 3 2
Z + (B - 1)*Z + (A - 2*B - 3*B )*Z + B + B - A*B=0 -0.02
0 0.2 0.4 0.6 0.8 1

Z3 + p*Z2 + q*Z + r=0

-0.04
LK tablas 0.63807 44.665 p -0,8954 m -0,0444
-0.06
q 0,222888 N -0,0232
r -0,03652 Z
LK analítico 0.63369 44.358 Discriminante n2 m3
D 0,00013104 D
Caso I 4 27
RKS 0.64635 45.245 M
N
0,28453243
0,05200179 Calcular la fugacidad para n-C4
a 500 °K y 60 bar
PR 0,61591 43,113 Z, solucion = 0,63501754
24
23 6
 Z( P  T)  A  CA( P  T)
B  CB( P  T) Fugacidad y Coeficiente de Fugacidad
Tb, °k =
n-Btutano
272,66
Tr =
Pr =
0,9409
0,6628
 A  B  B2  B3 
 
de Componente “i” en Solución
Tc (°k) = 425,12 m= 0,6723 2
V   A  3 B  2 B 
Pc(bar) = 37,96 alfa = 1,0407  B1  Para componente “i” en una mezcla de gases reales o en una
w= 0,2 A= 0,3563  
T, K = 400 B= 0,0548  1  solución líquida:
ANTB 2764,805 vapor  0,711721 ZZ  polyroots ( V)
Pini, bar = 25,232 liquido  0,711721 for i  0  2
ZZ i  0 if ImZZ i  0 iig  i (T )  RT ln  yi P 
Pfinal, bar = 25,15948 Fun. covergencia = 1,2064E-12 ZZ  sort ( ZZ )
solución de la EOS
Z3 -(1-B)Z2 + (A-2B-3B2)Z -(AB-B2-B3) =0 ZZ  ZZ  ZZ  10 5 
if Fugacidad del componente
i  i (T )  RT ln  fˆi 
cubica : 0  0
2 
3 2
Z -kZ + mZ -n =0 “i“ en la solución
ZZ  ZZ if  ZZ  10 
5
k m n p =m/3-(k/3)2 q=(k/3)*(p+m/6)-n/2 2 0  2 
0,9452 0,2376 0,0164 -0,02005 -0,00202 ZZ

Caso II , h < 0 h cos (theta) = 0,7104 V(P)=L(P) Fases multiples a la misma T y P están en equilibrio cuando la
-3,994E-06 theta = 0,7807 P= 25.159482 bar fugacidad de cada especie es igual en toda las fases:
pi = 3,14159 A= 0.356257
k Zk B= 0.054803

fˆi  fˆi  ...  fˆi

Zvapor 0 0,588745 a  
Z  f, bar
Zliquido 1 0,115113 Liquido 0.115113 0.711721 17.90654
2 0,241339 Vapor 0.588745 0.711721 17.90654
26
25

Propiedad residual: M R  M  M ig Relación Fundamental de la Propiedad Residual

Propiedad residual parcial: M iR  M i  M iig Gi  Gi  Gi
R ig
 nG  1 nG
d  d (nG)  dT
 RT  RT RT 2
i  i (T )  RT ln  fˆi 
d (nG)  (nV )dP  (nS )dT   i dni G  H  TS
 (nG)    i (T )  RT ln  yi P 
ig i
i     Gi i
 ni  P ,T ,n  nG  nV
d  dP 
nH G
dT   i dni
nG
 f ( P, T , ni )
j
 fˆ   RT  RT RT 2
i    RT ln  i 
ig
i
i RT RT
 yi P 
G/RT como una función de sus variables canónicas permite la
fˆ fˆi  ˆi yi P
evaluación de todas las otras propiedades termodinámicas.
Gi R  0 Para gas ideal
Gi R  RT ln ˆi ˆi  i
yi P  nG R  nV R nH R
R
 nG R  nV R nH R
dT   ln ˆi dni
G
d    dP  dT   i dni d    dP 
fˆ ig
fˆi ig  yi P RT 2
2
ˆiig  i  1 Coeficiente fugacidad de la  RT  RT RT i RT  RT  RT i
yi P especie “i “ en la solución
28
27 7
  nG R  nV R R Desarrollo de una ecuación general para el calculo de ln ˆi a
nH R Gi
d    dP  dT   dni partir de datos del factor de compresibilidad.
 RT  RT RT 2 i RT  (nG R / RT ) 
ln ˆi   
 nG R  nV R nH R  ni  P ,T ,n j
d    dP  2
dT   ln ˆi dni
 RT  RT RT nG R P dP
RT 0
i
 (nZ  n)
P
V R   (G R / RT ) 
Fijo T y composión:  
RT  P T , x
P   ( nZ  n) 
ln ˆi    
dP
0
 ni  P ,T ,n j P
Fijo P y composición: HR   (G R / RT )  n  (nZ )
 T   1  Zi
RT  T  P, x ni ni
 (nG R / RT ) 
Fijo T y P: ln ˆi    ln ˆi   ( Z i  1)
P dP Integración a Temperatura
 ni  P ,T ,n j 0 P
y composión constante
29 30

Coeficiente de Fugacidad con la EOS Virial

1   (nB)  P   (nB) 
ln ˆ1 
P
• Ecuación virial:   dP 
BP
Z  1   
RT RT 0
 n1  T ,n2 RT  n1  T ,n
2

– Segundo coeficiente virial para mezclas: B   yi y j Bij

i j B   yi y j Bij  y1  y1 B11  y2 B12   y2  y1 B21  y2 B22 
i j
– Mezcla binaria: B  y1 y1B11  y1 y2 B12  y2 y1B21  y2 y2 B22 Mezcla
B  y1 1  y2  B11  y2 1  y1  B22  2y1 y2 B12
nBP
Binaria
• “n” mol de mezcla de gases: nZ  n  B  y1 B11  y2 B22  y1 y2  2B12  B11  B22 
RT 12
  (nZ )  P   (nB) 
Z1     1   B  y1B11  y2 B22  y1 y212 12  2B12  B11  B22
 n1  P ,T ,n2 RT  n1  T ,n
2

1   (nB)  P   (nB)   
ln ˆ1 
P
 B11  y2212 
RT 
Similarmente  
ln ˆ2 
P
 B22  y1212 
RT 
ln ˆ1 
P

RT 
0 
 n
 dP 
1  T , n2 RT
 
 n1  T ,n2 32
31 8
 B   yi y j Bij  y1  y1 B11  y2 B12   y2  y1 B21  y2 B22 
1 2 i j

B11 + B22 EOS Virial : B  y1 1  y2  B11  y2 1  y1  B22  2y1 y2 B12 Mezcla

1 2 Binaria
B  y1 B11  y2 B22  y1 y2  2B12  B11  B22  12
Bmix
P  2 
1 1  
ln ˆ1   2 y j Bij  B  
RT  j 1 
P
RT
2y1 B11  2y2 B12   y1B11  y2 B22  y1 y212 
2 2
Bmix   yi y j Bi , j   1  y2  B11  2y2 B12  y2 B22  1  y2  y212
i 1 j 1 0
 B11  y2  2B12  B11  B22  12   y2212
2 2
 B11  y2212

Bmix   y1 B11   y1 B11   2  y1 y2 B12   y2 B22   y2 B22 

33
 
ln ˆ1 
P
 B11  y2212 
RT   
ln ˆ2 
P
 B22  y1212 
RT 
34

Estime ˆ1 y ˆ2 para una mezcla equimolar de metiletilcetona (1) /

Para mezcla de gases: multicomponentes
tolueno (2) a 50°C y 25 kPa. Asumir todo los kij = 0.
P    
ln ˆk  
RT 
1
 Bkk   yi y j (2 ik   ij ) 
2 i j 

 
ln ˆi 
P
RT
 2 y j Bij  B 
 j 
i,j
11
Tci,j /K
535.5
Pci,j /bar
41.5
Vci,j /cc molg-1
267
Zci,j
0.249
Wi,j
0.323
22 591.8 41.1 316 0.264 0.262
RTcij i   j  ik  2Bik  Bii  Bkk 12 562.947 41.281 290.812 0.257 0.293
Bij  ( B 0  ij B1 ) ij 
Pcij 2
i,j Tri,j B0 B1 Bi,j /cc mol-1 ˆi
Parámetro empírico de interaccion
11 0.603 -0.864 -1.296 -1375.909 0.9823
22 0.546 -1.028 -2.045 -1872.186 0.9760
Prausnitz et al. 1986 12 0.574 -0.943 -1.631 -1609.836
Tcij  (TciTcj ) (1  kij )
12 28.4235 ln ˆ1 
P
 B11  y2212  ˆ1  0.987
Z cij RTcij Z ci  Z cj V 1/ 3
V 1/ 3

3 RT
Pcij  Z cij  Vcij  
ci cj

Vcij 2  2  12  2B12  B11  B22 ln ˆ2 
P
 B22  y1212  ˆ2  0.983
  RT 9
35 36
 b1 b2
1 2 EOS2 cúbica
a1 a2 bmix   yi bi , j
1 2 i 1

bmix  y1b1  y2b2

1 1 2 2
amix   yi y j ai , j
i 1 j 1

Parámetro empírico de interacción

2 2

a1,2  a1a2 1  k1,2 

amix   y1 a1   y1 a1   2   y1 a1   y2 a2  1  k1,2    y2 a2   y2 a2 

37

38

EOS RKS
Componentes Puros
mi =0.48508  1.55171  wi  0.15613  wi
2
mi =0.48  1.574  wi  0.176  wi (API)
2

Hidrógeno
 
2
Pri Pri a i  1  mi 1  Tri 
Ai  0.42748a i Bi  0.08664  a =1.202  exp  0.30288Tr 
Tri 2 Tri

Mezclas A    xi x j Ai , j Ai , j  Ai A j 1  ki , j  B   xi Bi
i

Z 3  Z 2   A  B  B 2  Z  AB  0

A  B
Puros ln( )  Z  1  ln  Z  B   ln 1  
B  Z
Bi A2 B   B 10
Mezclas ln(ˆi )   Z  1  ln  Z  B     x j Ai , j  i  ln 1  
B BA i B   Z
39 40
 Para una mezcla gas de i-C4(1) y i -j Tc, °K Pc, bar
Vc,
Zc w
EOS Peng Robinson CO2(2)a 220 °F y 830 psia (57.21 cm3/molg
1-1 407,85 36,4 262,7 0,278 0,186
bar), con y1=0.5851 y y2=0.4149. 2-2 304,12 73,74 94,07 0,274 0,225
Puros  wi  0.49 
Calcular el coeficiente de fugacidad
0.37464 + 1.54226  wi  0.26992  wi 1-2 352,187 49,209 164,228 0,276 0,2055
2
 PR76
mi    utilizando: (a) EOS virial y (b) EOS
0.37964  1.48503  wi  0.16442  wi  0.016666  wi
 wi  0.49 
2 3
D. B. Robinson  PR78 Ding-Yu Peng de Peng Robinson, i -j Tr Pr B0 B1 B
1-1 0,926 1,572 -0,3944 -0,0988 -384,51
Twu
   
2
a i  1  mi 1  Tri 
Pri Pri L 1Tr N M • 2-2 1,242 0,776 -0,2155 0,0697 -68,52
Ai  0.45724a i
a  Tr N  M 1 e
Bi  0.0778 VEOS M. ABBOTT
Tri 2 Tri  1-2 1,072 1,163 -0,2945 0,0106 -173,93

Bmix   yi  y j  Bi , j  - 227.54 cm3 / molg

Mezclas A   xi x j Ai , j Ai , j  Ai Aj 1  ki , j  B   xi Bi ˆi fˆi , bar
i j
i
i
 P  
ˆi  exp   2   y j  Bi , j  Bmix   fˆi  ˆi  yi  P
Z 3   B  1 Z 2   A  2B  3B 2  Z  AB  B 2  B 3  0  R  T  j  
1 0.5068 16.942
2 0.9421 22.423

ln( )  Z  1  ln  Z  B  
A
ln 

 Z  1 2 B 
  EOS Peng – Robinson: k1,2=0.14
Puros
 
2 B  Z  1 2 B 
1.5

  mi Tri Pri ai Ai Bi
Mezclas
B A 2 B
ln(ˆi )   Z  1 i  ln  Z  B   1.5   x j Ai , j  i  ln 

 Z  1 2 B 
  1 0.6522 0.9258 1.5722 1.0499 0.8805 0.1321
B 2 BA i
 
B   Z  1 2 B 
  2 0.7080 1.2416 0.7761 0.8447 0.1944 0.0486
41 42

Z 3   B  1  Z 2   A  2 B  3 B 2   Z  B 3  B 2  A  B  0 Determine el coeficiente de fugacidad para el nitrógeno (1) y n-

Decano (2) a 280 °F y 800 psia para una mezcla con x1=0.0825
k1,2=0.14 k1,2=0 e y1=0.99365.

Av  0.50694 Bv  0.09739 Av  0.53509 Bv  0.09739 i -j Tr B0 B1 B yi ˆi fˆi bar

1-1 3.256 0.019 0.138 7.53
Zv  0.48143 Zv  0.35523 2-2 0.665 -0.727 -0.814 -2744.12
0.99365 1.012 55.488
1-2 1.472 -0.144 0.105 -116.75 0.00635 0.644 0.226
i ˆi fˆi bar i ˆi fˆi bar
1 0.4196 14.0265 1 0.3696 12.3542 EOS Peng Robinson EOS RKS
AV = 0.03312 Av= 0.02385
2 1.1352 27.0184 2 1.1867 28.2433 k12 = 0.155 k12 =0.15 Bv = 0.04505
Bv = 0.04046
Zv= 1.0113 Zv 1.0237
EOS Soave: k1,2=0.10 Zv= 0.46552 k1,2=0 Zv= 0.37061
i ˆi fˆi bar yi ˆi fˆi bar
yi ˆi fˆi bar ˆi fˆi,bar
1 1.0110 55.4106 1 0.99365 1.0237 58.2116
1 0.5841 0.4228 14.1314 1 0.3831 12.8064
2 1.2330 29.3453 2 0.7901 0.2767 2 0.00635 0.8807 0.3200 11
2 0.4159 1.1809 28.1068
43 44
 Para el sistema N2, C1 y n-C6, a 1500 psia y 100 F,
x1=0.113, x2=0.083, y1=0.787, y2=0.202, calcular la
fugacidad. Resolver EOS RKS,
mi Tri alfai Pri Ai Bi Kannap:
N2 0,5438 2,4638 0,4764 3,0418 0,1021 0,1070
C1 0,5037 1,6313 0,7402 2,2488 0,2674 0,1194 k12 =0.0278
n-C6 0,9380 0,6125 1,4493 3,4189 5,6453 0,4836 k13 =0.14
Fase Líquida Fase vapor
k23 =0.0374
AL 3,9318 Av 0,1424
BL 0,4108 Bv 0,1136
Zinicial 0,497 Zinicial 1,000
ZL 0,5206 Zv 1,00030
EOS PR Conv. -7,94E-11 Conv. -3,73E-08

xi yi ˆiV ˆiL fˆi V fˆi L Kcal kexp % Error

N2 0,11 0,79 1,0206 7,2245 83,069 84,430 7,0787 6,9646 1,64

C1 0,08 0,20 0,9013 2,4142 18,829 20,724 2,6787 2,4337 10,06

C6 0,80 0,01 0,4285 0,0058 0,487 0,486 0,0136 0,0137 -0,33

45 46

Tc, K Pc, bar wi Zci Wa Wb Wc

N2 126,2 34 0,038 0,289 0,4413 0,0824 0,0370
CH4 190,6 45,99 0,012 0,286 0,4435 0,0818 0,0426
n-C6 507,6 30,25 0,301 0,266 0,4588 0,0776 0,0800

Fase Líquida Fase Vapor

AL 4,0280 Av 0,1542
BL 0,3690 Bv 0,1076 EOS PTV
CL 0,3691 Cv 0,0527
Zinicial 0,437 Zinicial 0,988 k12 0,035
ZL 0,4636 Zv 0,98727
k13 0,15 Liquid Vapor
Conv -7,56E-09 Conv -2,44E-07
k23 0 phase phase
N 0,5219 N 0,1100
Nitrogen 0,1130 0,1130 0,7873
M -0,1529 M -0,0298
Methane 0,0830 0,0830 0,2018
Q 0,8910 Q 0,1901
Hexane 0,8040 0,8040 0,0109
xi yi ˆiV ˆiL fˆi V fˆi L Kcal kexp % Error
N2 0,11 0,79 1,0069 7,7885 81,955 91,021 7,7351 6,9646 11,06
C1 0,08 0,20 0,8856 2,2763 18,502 19,540 2,5703 2,4337 5,61 Presión de
C6 0,80 0,01 0,3896 0,0057 0,443 0,471 0,0145 0,0137 6,20 burbuja
47
12
P= 1629 psia
48
 Vapor phase Liquid phase
Nitrogen 0,1130 0,1130 0,000067
Methane 0,0830 0,0830 0,000171
Hexane 0,8040 0,8040 0,999762

Presión de 50
Rocío
P= 6.04 psia
49

13