CALCULO DE PROPIEDADES CON LA ECUACIÓN DE ESTADO VIRIAL

DATOS DE LOS COMPONENTES PUROS

Componente Tc (K) Pc (bar) ω Vc (cm3/mol) Zc "The Properties of Gases and
(1) Metano 190.56 45.99 0.011 98.60 0.286 Liquids" 5th Edition
(2) n-Hexano 507.60 30.25 0.300 368.00 0.264 B.E. Poling, J.M. Prausnitz, J.P.
(3) CO2 304.12 73.74 0.225 94.07 0.274 O'Connell
¿ ¿ ¿
Δh form Δg form Δsform Cp/R = a + b T + c T2 + d T3 + e T4 ; T en K

Componente kJ/mol kJ/mol J/mol K a 103 b 105 c 108 d 1011 e

(1) Metano -74.52 -50.45 -80.731 4.568 -8.975 3.631 -3.407 1.091
(2) n-Hexano -166.92 0.15 -560.36 8.831 -0.166 14.302 -18.314 7.124
(3) CO2 -393.51 -394.38 2.918 3.259 1.356 1.502 -2.374 1.056

PARÁMETROS BINARIOS
k12 = 0.080 k13 = 0.050 k23 = 0.200
y1 = 0.250
DATOS DE LA CORRIENTE T (K) = 323.15 y2 = 0.420
DE PROCESO P (bar) = 20.00 y3 = 0.330

A. CÁLCULO DE LOS SEGUNDOS COEFICIENTES VIRIALES

R Tc
B( o ) = 0 . 083 − 0 . 422/T 1.r 6 B= ( B( o ) + ω B( 1))
Pc

∂ B R 0 . 6752 0 . 7224
B( 1) = 0 . 139 − 0 . 172/T 4r .2 =
∂ T Pc T r 2 . 6 (
+ω
T 5.r 2 )
R= 83.14 cm3 bar / mol K
Bij dBij / dT

Pareja i-j Tc (K) Tr B(0) B(1) ω Vc (cm3/mol) Zc Pc (bar) cm3/mol cm3/mol K

1-1 190.56 1.69579 -0.098266 0.120286 0.0110 45.990 -33.40 0.53308
2-2 507.60 0.63662 -0.786161 -1.007098 0.3000 30.250 -1518.3 6.06066
3-3 304.12 1.06257 -0.299946 0.005703 0.2250 73.740 -102.41 0.90140
1-2 286.13 1.12938 -0.264351 0.035818 0.1555 204.64 0.2750 31.967 -192.58 1.82985
1-3 228.70 1.41300 -0.159708 0.098735 0.1180 96.32 0.2800 55.275 -50.93 1.18727
2-3 314.32 1.02809 -0.320704 -0.014108 0.2625 200.92 0.2690 34.987 -242.31 2.01344

nc nc

B M = ∑ ∑ y i y j Bi j = y 21 B11 + y 22 B22 + y 23 B33 +2 y 1 y 2 B 12 + 2 y1 y 3 B 13 + 2 y 2 y 3 B 23

i=1 j=1

n n
∂ BM c c
∂ Bi j ∂ B11 ∂ B 22 2 ∂ B33 ∂ B12 ∂ B 13 ∂ B 23
= ∑ ∑ yi y j = y 21 + y 22 + y3 +2 y1 y 2 + 2 y1 y3 + 2 y2 y3
∂T i=1 j=1 ∂T ∂T ∂T ∂T ∂T ∂T ∂T

BM = -397.08 cm3/mol dBM / dT = 2.3389 cm3/mol K

 B. CÁLCULO DE LAS PROPIEDADES DEL GAS IDEAL

¿ ¿
[
hi ( T , P ) = hi ( T o , Po ) + R ai ( T −T o ) +
bi
2
( T 2−T 2o) + 3i ( T 3−T 3o )+
c di
4
(T 4−T 4o) + 5i (T 5−T 5o )
e
]
¿ ¿
s i ( T , P ) = si ( T o , Po ) + R ai ln
nc
[ ( ) T
To
ci di ei
+ bi ( T −T o ) + (T 2 −T 2o ) + (T 3 −T 3o ) + (T 4 −T 4o ) − R ln
2 3 4
P
Po ] ( )
h M ( T , P )= ∑ y i hi ( T , P ) = y 1 h1 + y 2 h2 + y 3 h3
¿ ¿ ¿ ¿ ¿

n c i=1 nc

sM ( T , P ) = ∑ yi s ( T , P ) − R ∑
¿ ¿ ¿ ¿ ¿
i y i ln ( y i) = y 1 s 1 + y 2 s2 + y 3 s3 − R [ y 1 ln ( y 1 ) + y 2 ln ( y 2 ) + y 3 ln ( y3 ) ]
i=1 i=1

Estado de Referencia: hi#(To,Po) si#(To,Po)

To (K) = 298.15 J/mol J/mol K
Po (bar) = 1.00 (1) Metano -74,520 -80.731
(2) n-Hexano -166,920 -560.356
(3) CO2 -393,510 2.918
R= 8.314 J / mol K
¿
Cp P
()
T T
∫T Cp d T
o
¿
hi#(T,P) ∫T o T
d T − R ln
Po
si#(T,P)

Componente yi J/mol J/mol J/mol K J/mol K J/mol K

(1) Metano 0.250 907.21 -73,612.79 2.92138 -24.90652 -102.7163
(2) n-Hexano 0.420 3690.42 -163,229.58 11.88112 -24.90652 -573.3809
(3) CO2 0.330 938.72 -392,571.28 3.02285 -24.90652 -18.9657

h#M (T,P) = -216,508.14 J/mol s#M (T,P) = -263.8054 J/mol K

C. CÁLCULO DE LAS PROPIEDADES RESIDUALES

∂ BM
Δ h res = hM ( T , P ) − h¿M ( T , P ) = B M − T
∂T
P ( ) Δ h res = -2305.77 J/mol

∂ BM
¿
Δ sres = s M ( T , P ) − s M ( T , P ) = −
∂T
P ( ) Δ sres = -4.6778 J/mol K

D. CÁLCULO DE LAS PROPIEDADES REALES DE LA CORRIENTE DE PROCESO

¿
h M ( T , P ) = h M ( T , P ) + Δ h res
P ^
[ ] RT ¿
1
= 2( y
1
B
11
+ y
2
B
12
+y
3
B
13
)−B
M

ln ϕ
¿
¿
s M ( T , P ) = s M ( T , P ) + Δ s res ^
¿
2 [
= 2( y
1
B
21
+ y
2
B
22
+y
3
B
23
)−B
M ]
P
RT

P
ln ϕ
¿
Z =1+B
RT
nc
^
¿
i
=
( 2 ∑
j =1
y
j
B
ij
−B
M ) P
RT
^
¿
3 [
= 2( y
1
B
31
+ y
2
B
32
+y
3
B
33
)−B
M ] RT
P

v=
RT
+B ln ϕ ln ϕ
¿
P ¿

h= -218,813.9 J / mol ϕ^ 1 = 1.1476513 f^ 1 = 5.73826 bar

z= 0.70441 s= -268.4831 J / mol K ϕ^ 2 = 0.4297584 f^ =

2
3.60997 bar

v= 946.258 cm3/mol g= -132,053.6 J / mol ϕ^ 3 = 1.0776688 f̂ 3  7.11261 bar

pseudo Vc
210.3 cm3/mol El cálculo está en la región de validez de la ecuación Virial
=
 B ( cm3 / mol )
2° Coeficiente Virial - Metano 2° Coeficiente Virial - CO2
50.00
50.00
B (cm3 / mol )

0.00
0.00
-50.00

-100.00 -50.00

-150.00
-100.00 exp
-200.00 calc1
-150.00
-250.00 exp
calc -200.00
-300.00

-350.00 -250.00
-400.00
100.0 200.0 300.0
Temperatura (400.0
K) 500.0 600.0 -300.00

-350.00
200.00 300.00 400.00 500.00 600.00
Temperatura ( K700.00
) 800.00 900.00
-500.00
2° Coeficiente Virial - nHexano
B ( cm3 / mol )

-700.00

-900.00

-1100.00

-1300.00

-1500.00

exp
-1700.00 calc

-1900.00

300.0 320.0 340.0 360.0 ( K ) 380.0

Temperatura 400.0 420.0
 CALCULO DE UN COMPRESOR CON LA ECUACIÓN DE ESTADO VIRIAL
DATOS DE LOS COMPONENTES PUROS
Componente Tc (K) Pc (bar) ω Vc (cm3/mol) Zc "The Properties of Gases and
(1) Metano 190.56 45.99 0.011 98.60 0.286 Liquids" 5th Edition
(2) n-Hexano 507.60 30.25 0.300 368.00 0.264 B.E. Poling, J.M. Prausnitz, J.P.
(3) CO2 304.12 73.74 0.225 94.07 0.274 O'Connell

¿ ¿ ¿
Δh form Δg form Δs form Cp/R = a + b T + c T2 + d T3 + e T4 ; T en K

Componente kJ/mol kJ/mol J/mol K a 103 b 105 c 108 d 1011 e

(1) Metano -74.52 -50.45 -80.731 4.568 -8.975 3.631 -3.407 1.091
(2) n-Hexano -166.92 0.15 -560.36 8.831 -0.166 14.302 -18.314 7.124
(3) CO2 -393.51 -394.38 2.918 3.259 1.356 1.502 -2.374 1.056
PARÁMETROS BINARIOS
k12 = 0.080 k13 = 0.050 k23 = 0.200
ALIMENTACIÓN SALIDA y1 = 0.250
DATOS DE LAS
CORRIENTES DE T1 (K) = 323.15 P2 (bar) = 50.00 y2 = 0.420
PROCESO P1 (bar) = 10.00 y3 = 0.330
Eficiencia Adiabática del Compresor (%) = 80.00
A. CÁLCULO DE LOS SEGUNDOS COEFICIENTES VIRIALES A LA TEMPERATURA DE ALIMENTACIÓN
R= 83.14 cm3 bar / mol K Bij dBij / dT

Pareja i-j Tc (K) Tr B(0) B(1) ω Vc (cm3/mol) Zc Pc (bar) cm3/mol cm3/mol K

1-1 190.56 1.69579 -0.098266 0.120286 0.0110 45.990 -33.40 0.53308
2-2 507.60 0.63662 -0.786161 -1.007098 0.3000 30.250 -1518.3 6.06066
3-3 304.12 1.06257 -0.299946 0.005703 0.2250 73.740 -102.41 0.90140
1-2 286.13 1.12938 -0.264351 0.035818 0.1555 204.64 0.2750 31.967 -192.58 1.82985
1-3 228.70 1.41300 -0.159708 0.098735 0.1180 96.32 0.2800 55.275 -50.93 1.18727
2-3 314.32 1.02809 -0.320704 -0.014108 0.2625 200.92 0.2690 34.987 -242.31 2.01344
BM = -397.08 cm /mol 3 dBM / dT = 2.33888 cm3/mol K
B. CÁLCULO DE LAS PROPIEDADES DEL GAS IDEAL A LA TEMPERATURA Y PRESIÓN DE LA ALIMENTACIÓN

Estado de Referencia: hi#(To,Po) si#(To,Po)

To (K) = 298.15 J/mol J/mol K
Po (bar) = 1.00 (1) Metano -74,520 -80.731
(2) n-Hexano -166,920 -560.356
(3) CO2 -393,510 2.918
R= 8.314 J / mol K
T C ¿p P
∫ T C¿p d T
o
hi#(T,P) ∫T
T

o T
d T − R ln ()
P o
si#(T,P)

Componente yi J/mol J/mol J/mol K J/mol K J/mol K

(1) Metano 0.250 907.21 -73,612.79 2.92138 -19.14369 -96.9535
(2) n-Hexano 0.420 3690.42 -163,229.58 11.88112 -19.14369 -567.6181
(3) CO2 0.330 938.72 -392,571.28 3.02285 -19.14369 -13.2029
h#M (T1,P1) = -216,508.14 J/mol s#M (T1,P1) = -258.0425 J/mol K
C. CÁLCULO DE LAS PROPIEDADES RESIDUALES A LA TEMPERATURA Y PRESIÓN DE LA ALIMENTACIÓN

 h res  -1152.88 J/mol Δ sres = -2.3389 J/mol K

D. CÁLCULO DE LAS PROPIEDADES REALES DE LA CORRIENTE DE ALIMENTACIÓN

h1= -217,661.0 J / mol ϕ^ 1 = 1.0712849 f^ 1 = 2.67821 bar

z= 0.852205 s1 = -260.3814 J / mol K ϕ^ 2 = 0.6555596 f^ =

2
2.75335 bar
v= 2289.593 cm3/mol g1 = -133,518.8 J / mol ϕ^ 3 = 1.0381083 f̂ 3  3.42576 bar
 CÁLCULO DEL COMPRESOR ISOENTRÓPICO
SUPONER: T2isoent (K) =
418.53 Esta temperatura debe ser tal que: Error=S2-S1=0
E. CÁLCULO DE LOS SEGUNDOS COEFICIENTES VIRIALES A LA TEMPERATURA DE SALIDA (Reversible)
R= 83.14 cm3 bar / mol K Bij dBij / dT
Pareja i-j Tc (K) Tr B(0) B(1) ω Vc (cm3/mol) Zc Pc (bar) cm3/mol cm3/mol K
1-1 190.56 2.19630 -0.036841 0.132684 0.0110 45.990 -12.19 1.01706
2-2 507.60 0.82452 -0.491633 -0.247797 0.3000 30.250 -789.6 3.28313
3-3 304.12 1.37619 -0.170179 0.094014 0.2250 73.740 -51.10 1.29613
1-2 286.13 1.46271 -0.146646 0.104177 0.1555 204.64 0.2750 31.967 -97.07 2.76406
1-3 228.70 1.83004 -0.077463 0.125411 0.1180 96.32 0.2800 55.275 -21.56 3.18091
2-3 314.32 1.33153 -0.183903 0.087327 0.2625 200.92 0.2690 34.987 -120.24 2.75938
BM = -202.88 cm /mol 3 dBM / dT = 2.65406 cm3/mol K
F. CÁLCULO DE LAS PROPIEDADES DEL GAS IDEAL A LA TEMPERATURA Y PRESIÓN DE SALIDA (Reversible)

C ¿p P
()
T T
∫T Cp d T
o
¿ h (T,P)
i
#
∫T T o
d T − R ln P
o
si#(T,P)

Componente yi J/mol J/mol J/mol K J/mol K J/mol K

(1) Metano 0.250 4640.12 -69,879.88 13.01670 -32.52456 -100.2390
(2) n-Hexano 0.420 19901.53 -147,018.47 55.64158 -32.52456 -537.2385
(3) CO2 0.330 4756.36 -388,753.64 13.35338 -32.52456 -16.2532
h#M (T2,P2) = -207,506.43 J/mol s#M (T2,P2) = -247.1111 J/mol K
G. CÁLCULO DE LAS PROPIEDADES RESIDUALES A LA TEMPERATURA Y PRESIÓN DE SALIDA (Reversible)
 h res  -6,568.39
J/mol Δ sres = -13.2703 J/mol K
H. CÁLCULO DE LAS PROPIEDADES REALES DE LA CORRIENTE DE SALIDA (Reversible)
h2 = -214,074.8 J / mol ^ =
ϕ 1
1.1562232 f^ 1 = 14.45279 bar
z= 0.7084706 s2
= -260.3814 J / mol K ^
ϕ = 0.4294277 f^ = 2 2
9.01798 bar
v= 493.043 cm3/mol g2 = -105,098.3 J / mol ^ =
ϕ 3
1.0859418 f̂ 3  17.91804 bar
Error = s2 - s1 = 0.00000 Wreal = Wrev / Efic = 4,482.8 J / mol
Wrev = h2 - h1 = 3,586.2 J / mol h2real = h1 + Wreal = -213,178.3 J / mol

CÁLCULO DE LA CORRIENTE REAL DE SALIDA

SUPONER: T2REAL (K) = 433.29 Esta temperatura debe ser tal que: Error=h2calc -h2 real=0
I. CÁLCULO DE LOS SEGUNDOS COEFICIENTES VIRIALES A LA TEMPERATURA DE SALIDA (Real)
R= 83.14 cm3 bar / mol K Bij dBij / dT
Pareja i-j Tc (K) Tr B(0) B(1) ω Vc (cm3/mol) Zc Pc (bar) cm3/mol cm3/mol K
1-1 190.56 2.27380 -0.030373 0.133540 0.0110 45.990 -9.96 1.17324
2-2 507.60 0.85361 -0.460618 -0.195371 0.3000 30.250 -724.4 3.06205
3-3 304.12 1.42475 -0.156514 0.100111 0.2250 73.740 -45.94 1.45806
1-2 286.13 1.51432 -0.134251 0.108897 0.1555 204.64 0.2750 31.967 -87.30 3.12483
1-3 228.70 1.89462 -0.068802 0.127253 0.1180 96.32 0.2800 55.275 -18.50 3.74959
2-3 314.32 1.37851 -0.169497 0.094331 0.2625 200.92 0.2690 34.987 -108.11 3.08834
BM = -184.76 cm3/mol dBM / dT = 2.90324 cm3/mol K
J. CÁLCULO DE LAS PROPIEDADES DEL GAS IDEAL A LA TEMPERATURA Y PRESIÓN DE SALIDA (Real)
¿
Cp P
()
T T
∫T Cp d T
o
¿ h (T,P)
i
#
∫T T o
d T − R ln P
o
si#(T,P)

Componente yi J/mol J/mol J/mol K J/mol K J/mol K

(1) Metano 0.250 5261.28 -69,258.72 14.47517 -32.52456 -98.7806
(2) n-Hexano 0.420 22713.11 -144,206.89 62.24303 -32.52456 -530.6371
(3) CO2 0.330 5379.41 -388,130.59 14.81635 -32.52456 -14.7902
h#M (T2,P2) = -205,964.67 J/mol s#M (T2,P2) = -243.4911 J/mol K
K. CÁLCULO DE LAS PROPIEDADES RESIDUALES A LA TEMPERATURA Y PRESIÓN DE SALIDA (Real)
 h res  -7,213.59 J/mol Δ sres = -14.5162 J/mol K
L. CÁLCULO DE LAS PROPIEDADES REALES DE LA CORRIENTE DE SALIDA (Real)
h2 = -213,178.3 J / mol ^ =
ϕ 1
1.139733 f^ 1 = 14.24666 bar
z= 0.7084706 s2 = -258.0073 J / mol K ^ =
ϕ 0.473444 2 f^ = 2 9.94232 bar
v= 493.043 cm /mol
3 g2 = -101,385.1 J / mol ^ =
ϕ 3 1.078393 f̂ 3  17.79348 bar
Error = h2 calc - h2 real = 0.000 s gen = s2 - s1 = 2.3741 J / mol K
 CALCULO DE EQUILIBRIO (GAS)-(SÓLIDO PURO) CON LA ECUACIÓN DE ESTADO VIRIAL

DATOS DE LOS COMPONENTES PUROS

Componente Tc (K) Pc (bar) ω Vc (cm3/mol) Zc "The Properties of Gases and
(1) Metano 190.56 45.99 0.011 98.60 0.286 Liquids" 5th Edition (2001)
(2) Nitrógeno 126.20 33.98 0.037 90.10 0.289 B.E. Poling, J.M. Prausnitz, J.P.
(3) Benceno 562.05 48.95 0.210 256.00 0.268 O'Connell

PARÁMETROS BINARIOS
k12 = 0.030 k13 = 0.080 k23 = 0.046

ð Calcular la presión a la cual comienza a condensar el benceno sólido

T
DATOS: T(K) = 253.15 y1 = 0.356
y2 = 0.643 P=?
y3 = 0.001 Gas con
composición
y1, y2, y3
P°3 (bar)= 0.0066670
DATOS DEL BENCENO SÓLIDO PURO A
LA TEMPERATURA DE EQUILIBRIO vs3 (cm3/mol)= 73.98 Sólido 3
PURO

^
f^ G3 = P y 3 ϕ^ G3
P
f^ G3 = f 3S
s
f =
3 ϕ3o Po3 exp (v s3 ( P − P3o )
RT )
¿
3 [
= 2( y
1
B
31
+ y
2
B
32
+y
3
B
33
)−B
M ] RT

B33 Po3
ln ϕ
¿ ϕ o3 = exp
RT( )
nc nc

B M = ∑ ∑ y i y j Bi j = y 21 B11 + y 22 B22 + y 23 B33 +2 y 1 y 2 B 12 + 2 y1 y 3 B 13 + 2 y 2 y 3 B 23

i=1 j=1

A. CÁLCULO DE LOS SEGUNDOS COEFICIENTES VIRIALES A LA TEMPERATURA DE EQUILIBRIO

R= 83.14 cm3 bar / mol K Bij

Pareja i-j Tc (K) Tr B(0) B(1) ω Vc (cm3/mol) Zc Pc (bar) cm3/mol

1-1 190.56 1.32845 -0.184891 0.086824 0.0110 45.990 -63.36
2-2 126.20 2.00594 -0.055549 0.129757 0.0370 33.980 -15.67
3-3 562.05 0.45040 -1.428986 -4.763255 0.2100 48.950 -2319.04
1-2 150.42 1.68291 -0.100491 0.119677 0.0240 94.29 0.2875 38.134 -32.01
1-3 301.09 0.84079 -0.473947 -0.217323 0.1105 164.99 0.2770 42.026 -296.60
2-3 254.08 0.99635 -0.341475 -0.035661 0.1235 158.90 0.2785 37.024 -197.34
F u g a c id a d d e l B e n c e n o ( b a r) .

o
ϕ3 = 0.999266 BM = -29.63 cm3/mol
B. CÁLCULO DE LA PRESIÓN (Iterativo)
s
f3 ϕ^ 3
G
f^ G3
P
bar
exp ( v3s ( P − Po3 )
RT ) bar bar

2.0000 1.0070 0.006709 0.9590536 0.001918 0.010000

0.009000
4.0000 1.0141 0.006756 0.9197839 0.003679
0.008000

6.0000 1.0213 0.006804 0.8821221 0.005293 0.007000

8.0000 1.0285 0.006852 0.8460024 0.006768 0.006000

0.005000
10.0000 1.0358 0.006900 0.8113617 0.008114 0.004000 Sóli
do
12.0000 1.0431 0.006949 0.7781393 0.009338 0.003000 Gas

0.002000
8.1243 1.0289 0.006855 0.8438063 0.006855
0.001000
^
Error = f 3 − f 3S =
G
0.000000
0.0000 2.0000 Presión
4.0000 (bar)
6.0000 8.0000 10.0000 12.0000
 0.004000
do
0.003000 Gas

0.002000

0.001000
Error = ^f G3 − f 3S = 0.000000 0.000000
Presión (bar)
0.0000 2.0000 4.0000 6.0000 8.0000 10.0000 12.0000
 CALCULO DE EQUILIBRIO (GAS)-(SÓLIDO PURO) CON LA ECUACIÓN DE ESTADO VIRIAL

DATOS DE LOS COMPONENTES PUROS

Componente Tc (K) Pc (bar) ω Vc (cm3/mol) Zc "The Properties of Gases and
(1) Metano 190.56 45.99 0.011 98.60 0.286 Liquids" 5th Edition (2001)
(2) Nitrógeno 126.20 33.98 0.037 90.10 0.289 B.E. Poling, J.M. Prausnitz, J.P.
(3) Benceno 562.05 48.95 0.210 256.00 0.268 O'Connell

PARÁMETROS BINARIOS
k12 = 0.030 k13 = 0.080 k23 = 0.046

ð Calcular la solubilidad del benceno sólido en el gas

T
DATOS: T(K) = 253.15 z1 = 0.400
P (bar) = 60.00 z2 = 0.600 Gas con
composición
P
z3 = 0.000 Inicial: z1, z2
Final: y1, y2, y3
P°3 (bar)= 0.0066670
DATOS DEL BENCENO SÓLIDO PURO A
LA TEMPERATURA DE EQUILIBRIO vs3 (cm3/mol)= 73.98 Sólido 3
PURO

f^ G3 = P y 3 ϕ^ G3 
 v s P  P3o
f 3s   3o P3o exp  3
 

f^ G3 = f 3S
^ P  RT 
¿
3 [
= 2( y
1
B
31
+ y
2
B
32
+y
3
B
33
)−B
M ] RT
ϕ B33 Po3
ln

nc nc
¿ ϕ o3 = exp (RT )
B M = ∑ ∑ y i y j Bi j = y 1 B11 + y 2 B22 + y 3 B33 +2 y 1 y 2 B 12 + 2 y1 y 3 B 13 + 2 y 2 y 3 B 23
2 2 2

i=1 j=1

A. CÁLCULO DE LOS SEGUNDOS COEFICIENTES VIRIALES A LA TEMPERATURA DE EQUILIBRIO

R= 83.14 cm3 bar / mol K Bij

Pareja i-j Tc (K) Tr B(0) B(1) ω Vc (cm3/mol) Zc Pc (bar) cm3/mol

1-1 190.56 1.32845 -0.184891 0.086824 0.0110 45.990 -63.36
2-2 126.20 2.00594 -0.055549 0.129757 0.0370 33.980 -15.67
3-3 562.05 0.45040 -1.428986 -4.763255 0.2100 48.950 -2319.04
1-2 150.42 1.68291 -0.100491 0.119677 0.0240 94.29 0.2875 38.134 -32.01
1-3 301.09 0.84079 -0.473947 -0.217323 0.1105 164.99 0.2770 42.026 -296.60
2-3 254.08 0.99635 -0.341475 -0.035661 0.1235 158.90 0.2785 37.024 -197.34

B. CÁLCULO DE LA FUGACIDAD DEL SÓLIDO s

f
ϕ o3 = 0.999266 exp ( v3s ( P − Po3 )
RT ) = 1.234766 3 =0.0082261 bar

C. BALANCE DE MATERIA: y1 = z1 (1 - y3) y2 = z2 (1 - y3)

D. CÁLCULO DE LA SOLUBILIDAD DEL BENCENO SÓLIDO EN EL GAS (Iterativo)

s
f
Inicialización como gas ideal:
y3= 3 = 0.0001371
P f3
s
BM ^G
ϕ f^ G3 y 3=
3
Iteración y3 y1 y2 cm3/mol bar P ϕ^ v 3
1 0.0001371024 0.399945 0.599918 -31.20 0.282463 0.0023236 0.0004853810
2 0.0004853810 0.399806 0.599709 -31.35 0.281413 0.0081956 0.0004871920
3 0.0004871920 0.399805 0.599708 -31.35 0.281408 0.0082260 0.0004872015
4 0.0004872015 0.399805 0.599708 -31.35 0.281408 0.0082261 0.0004872015
5 0.0004872015 0.399805 0.599708 -31.35 0.281408 0.0082261 0.0004872015
 CALCULO DE EQUILIBRIO GAS-SÓLIDO PURO CON LA ECUACIÓN DE ESTADO VIRIAL
DATOS DE LOS COMPONENTES PUROS

Componente Tc (K) Pc (bar) w Vc (cm3/mol) Zc

(1) Nitrógeno 126.2 33.9 0.040 89.5 0.290
(2) Metano 190.6 46.0 0.008 99.0 0.288
(3) Benceno 562.1 48.9 0.212 259.0 0.271

PARÁMETROS BINARIOS
k12 = 0.030 k13 = 0.046 k23 = 0.080

T (K) = 253.15 R= 83.14 cm bar / mol K

3

P (bar) = 60.794

CALCULO DE LA FUGACIDAD DEL SÓLIDO

DATOS INICIALES DATOS DE LA FASE SOLIDA ( B33 o
P3 / RT )
FASE VAPOR
(y1)inic= P°3 (bar)=
ϕ o3 = e = 0.9992618
0.500 0.0066670 s
v3 ( P−P3 )/RT
o

(y2) =
inic
0.500 v 3 (cm /mol)=
s 3

f s3 (bar)=
73.98
0.0082491
f s3 =ϕ o3 P o3 e = 0.0082491001

CÁLCULO DE LOS SEGUNDOS COEFICIENTES VIRIALES A LA TEMPERATURA DEL SISTEMA

Bij
Pareja i-j Tc (K) Tr B(0) B(1) w Vc (cm3/mol) Zc Pc (bar) (cm3/mol)
1-1 126.20 2.00594 -0.055549 0.129757 0.0400 33.900 -15.59
2-2 190.60 1.32817 -0.184981 0.086778 0.0080 46.000 -63.48
3-3 562.10 0.45036 -1.429201 -4.765087 0.2120 48.900 -2331.30
1-2 150.44 1.68273 -0.100522 0.119669 0.0240 94.17 0.2890 38.385 -31.82
1-3 254.09 0.99631 -0.341506 -0.035694 0.1260 159.55 0.2805 37.140 -196.80
2-3 301.13 0.84066 -0.474080 -0.217547 0.1100 166.39 0.2795 42.055 -296.48

CÁLCULO DE LA SOLUBILIDAD EN LA FASE VAPOR COMO GAS IDEAL

f s3
y3= = 0.000136

PROCEDIMIENTO ITERATIVO PARA EL CALCULO DE LA SOLUBILIDAD (SUBSTITUCIONES SUCESIVAS)

y i=( y i )inic ( 1− y 3 ) ; i=1, 2

P ^ P ^ P
^
¿
1
=
[ 2( y 1 B 11+ y
2
B
12
+ y
3
B
13
)− B
M ] RT
¿
2 [
= 2( y
1
B
21
+ y
2
B
22
+y
3
B
23
)− B
M ] RT
¿
3
=
[ 2( y 1 B31 + y
2
B
32
+ y
3
B
33
)− B
M ] RT
ln ϕ ln ϕ ln ϕ
¿ ¿ ¿

ϕ^ 1 ϕ^ 2 ϕ^ 3
y3 y1 y2
V V V BM f s3
ITERACION
(cm3/mol)
z y 3= ERROR
P ϕ^ v
3
1 0.000110000 0.499945 0.499945 0.966714 0.841764 0.266341 -35.7238 0.8968118 0.0005094566
2 0.000509457 0.499745 0.499745 0.966328 0.841281 0.265063 -35.8928 0.8963237 0.0005119133 3.631423541
3 0.000511913 0.499744 0.499744 0.966326 0.841278 0.265055 -35.8938 0.8963207 0.0005119285 0.004822276
4 0.000511928 0.499744 0.499744 0.966326 0.841278 0.265055 -35.8938 0.8963206 0.0005119286 0.000029587
5 0.000511929 0.499744 0.499744 0.966326 0.841278 0.265055 -35.8938 0.8963206 0.0005119286 0.000000182
6 0.000511929 0.499744 0.499744 0.966326 0.841278 0.265055 -35.8938 0.8963206 0.0005119286 0.000000001
7 0.000511929 0.499744 0.499744 0.966326 0.841278 0.265055 -35.8938 0.8963206 0.0005119286 0.000000000
SISTEMA : CO2(1)-CH4(2)

Capacidad calorífica:
Cp=a + b T+ c T2 + d T3 (Cp en J/mol K; T en K)
a b c d
CO2 19.80 7.34E-02 -5.60E-05 1.72E-08
CH4 19.25 5.213E-02 1.197E-05 -1.132E-08

Segundo coeficiente virial:

Bij=c1+c2 T-1 +c3 T-2 (B en cm3/mol, T en K)
Interacción c1 c2 c3
CO2-CO2 42.5 -1.68E+04 -2.51E+06
CH4-CH4 40.4 -2.54E+04 -6.87E+06
CO2-CH4 41.4 -1.95E+04 -3.73E+06

Estado de referencia:
CO2:
To(°C)= 298.15 (h)ref= 0 J/mol
Po(bar)= 1 (s)ref= 0 J/mol K

R(J/molK)= 8.314 CH4:

(h)ref= 0 J/mol
(s)ref= 0 J/mol K

Cálculos:

T(K)= 323.15 323.15 323.15 323.15 323.15 323.15 323.15 323.15 323.15 323.15 323.15
P(bar)= 100 100 100 100 100 100 100 100 100 100 100
X1= 0.001 0.005 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09
X2= 0.999 0.995 0.99 0.98 0.97 0.96 0.95 0.94 0.93 0.92 0.91
Gas ideal:
1 puro:
hid1(J/mol)= 943.0016 943.0016 943.0016 943.0016 943.0016 943.0016 943.0016 943.0016 943.0016 943.0016 943.0016
sid1(J/mol K)= -35.2508 -35.2508 -35.2508 -35.2508 -35.2508 -35.2508 -35.2508 -35.2508 -35.2508 -35.2508 -35.2508
2 puro:
h 2(J/mol)=
id
906.5011 906.5011 906.5011 906.5011 906.5011 906.5011 906.5011 906.5011 906.5011 906.5011 906.5011
sid2(J/mol K)= -35.3685 -35.3685 -35.3685 -35.3685 -35.3685 -35.3685 -35.3685 -35.3685 -35.3685 -35.3685 -35.3685
Mezcla:
hid(J/mol)= 906.5376 906.6836 906.8661 907.2311 907.5961 907.9611 908.3261 908.6911 909.0562 909.4212 909.7862
sid(J/mol K)= -35.3026 -35.1062 -34.9017 -34.5510 -34.2447 -33.9675 -33.7121 -33.4744 -33.2515 -33.0414 -32.8426
Coef. Virial:
B11: -33.3218 -33.3218 -33.3218 -33.3218 -33.3218 -33.3218 -33.3218 -33.3218 -33.3218 -33.3218 -33.3218
B22: -103.9586 -103.9586 -103.9586 -103.9586 -103.9586 -103.9586 -103.9586 -103.9586 -103.9586 -103.9586 -103.9586
B12=B21: -54.6626 -54.6626 -54.6626 -54.6626 -54.6626 -54.6626 -54.6626 -54.6626 -54.6626 -54.6626 -54.6626
Bm: -103.8600 -103.4663 -102.9755 -101.9979 -101.0260 -100.0596 -99.0989 -98.1437 -97.1941 -96.2501 -95.3117
dB11/dT: 0.3089 0.3089 0.3089 0.3089 0.3089 0.3089 0.3089 0.3089 0.3089 0.3089 0.3089
dB22/dT: 0.6503 0.6503 0.6503 0.6503 0.6503 0.6503 0.6503 0.6503 0.6503 0.6503 0.6503
dB12/dT: 0.4078 0.4078 0.4078 0.4078 0.4078 0.4078 0.4078 0.4078 0.4078 0.4078 0.4078
dBm/dT: 0.6498 0.6479 0.6455 0.6407 0.6359 0.6311 0.6264 0.6217 0.6171 0.6124 0.6078
Residual:
1 puro:
h-hid1(J/mol): -1,331.32 -1,331.32 -1,331.32 -1,331.32 -1,331.32 -1,331.32 -1,331.32 -1,331.32 -1,331.32 -1,331.32 -1,331.32
s-sid1(J/mol K): -3.09 -3.09 -3.09 -3.09 -3.09 -3.09 -3.09 -3.09 -3.09 -3.09 -3.09
2 puro:
h-hid2(J/mol) -3,141.05 -3,141.05 -3,141.05 -3,141.05 -3,141.05 -3,141.05 -3,141.05 -3,141.05 -3,141.05 -3,141.05 -3,141.05
s-sid2(J/mol K): -6.50 -6.50 -6.50 -6.50 -6.50 -6.50 -6.50 -6.50 -6.50 -6.50 -6.50
Mezcla:
h-h (J/mol):
id
-3,138.50 -3,128.31 -3,115.60 -3,090.29 -3,065.13 -3,040.12 -3,015.25 -2,990.54 -2,965.97 -2,941.55 -2,917.29
s-sid(J/mol K): -6.50 -6.48 -6.45 -6.41 -6.36 -6.31 -6.26 -6.22 -6.17 -6.12 -6.08
Gas real:
1 puro:
h1(J/mol): -388.32 -388.32 -388.32 -388.32 -388.32 -388.32 -388.32 -388.32 -388.32 -388.32 -388.32
s1(J/mol K): -38.34 -38.34 -38.34 -38.34 -38.34 -38.34 -38.34 -38.34 -38.34 -38.34 -38.34
2 puro:
h2(J/mol) -2,234.55 -2,234.55 -2,234.55 -2,234.55 -2,234.55 -2,234.55 -2,234.55 -2,234.55 -2,234.55 -2,234.55 -2,234.55
s2(J/mol K): -41.87 -41.87 -41.87 -41.87 -41.87 -41.87 -41.87 -41.87 -41.87 -41.87 -41.87
Mezcla:
h(J/mol): -2,231.96 -2,221.62 -2,208.73 -2,183.06 -2,157.53 -2,132.15 -2,106.93 -2,081.85 -2,056.92 -2,032.13 -2,007.50
s(J/mol K): -41.80 -41.59 -41.36 -40.96 -40.60 -40.28 -39.98 -39.69 -39.42 -39.17 -38.92
Z: 0.61342 0.61489 0.61672 0.62036 0.62397 0.62757 0.63115 0.63470 0.63824 0.64175 0.64524
v(cm3/mol): 164.8069 165.2006 165.6914 166.6690 167.6409 168.6073 169.5680 170.5232 171.4728 172.4168 173.3552
Coeficiente de fugacidad:
1 puro:
f1: 0.88336 0.88336 0.88336 0.88336 0.88336 0.88336 0.88336 0.88336 0.88336 0.88336 0.88336
2 puro:
f2: 0.67913 0.67913 0.67913 0.67913 0.67913 0.67913 0.67913 0.67913 0.67913 0.67913 0.67913
Mezcla:
f1: 0.98002 0.97921 0.97820 0.97619 0.97421 0.97226 0.97033 0.96842 0.96654 0.96468 0.96285
f2: 0.67913 0.67913 0.67914 0.67916 0.67919 0.67924 0.67931 0.67938 0.67948 0.67958 0.67970
Energía de Gibbs:
1 puro:
g1(J/mol): 12,001 12,001 12,001 12,001 12,001 12,001 12,001 12,001 12,001 12,001 12,001
2 puro:
g2(J/mol): 11,296 11,296 11,296 11,296 11,296 11,296 11,296 11,296 11,296 11,296 11,296
Mezcla:
gm(J/mol): 11,276 11,217 11,156 11,052 10,964 10,884 10,811 10,745 10,682 10,624 10,570
propiedades de mezclado:
Dh(J/mol): 0.7427 3.6989 7.3605 14.5724 21.6355 28.5499 35.3157 41.9327 48.4011 54.7207 60.8916
Dg(J/mol): -21.0 -83.2 -147.7 -257.9 -353.9 -440.5 -520.1 -594.0 -663.2 -728.4 -789.9
Ds(J/mol K): 0.0672 0.2689 0.4798 0.8432 1.1620 1.4514 1.7186 1.9680 2.2022 2.4234 2.6329
Dg id(J/mol): -21.24418 -84.57383 -150.45759 -263.39866 -362.00762 -451.21035 -533.34477 -609.78663 -681.44392 -748.96106 -812.81902
 Mezcla CO2-CH4, T= 323.15 K Mezcla CO2-CH4, T= 323.15 K

200 7.0000

180

s J/mol K
6.0000
h J/mol

160

140 5.0000

120
4.0000
100
3.0000
80

60
2.0000
40
1.0000
20

0 0.0000
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Fracción mol de CO2 Fracción mol de CO2
P= 100 bar P= 50 bar P= 100 bar Ideal

Mezcla CO2-CH4, T= 323.15 K Mezcla CO2-CH4, T= 323.15 K

8.00000 0.00000

-200.00000
v cm3/mol

7.00000
g J/mol

-400.00000
6.00000
-600.00000
5.00000
-800.00000

4.00000 -1000.00000

3.00000
-1200.00000

-1400.00000
2.00000
-1600.00000
1.00000
-1800.00000

0.00000 -2000.00000
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Fracción mol de CO2 Fracción mol de CO2
A cualquier presión P= 100 bar Ideal
Coeficiente de fugacidad .

Mezcla CO2-CH4, T= 323.15 K

0.95

0.9

0.85

0.8

0.75

0.7

0.65

0.6
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Fracción mol de CO2
P= 100 bar (CO2) P= 100 bar (CH4)
P= 50 bar (CO2) P= 50 bar (CH4)