UNIVERSIDAD MAYOR REAL Y PONTIFICIA DE SAN FRANCISCO XAVIER DE

CHUQUISACA
FACULTAD DE TECNOLOGÍA
INGENIERÍA QUÍMICA

ASIGNATURA: PRQ204 transferencia de masa I

NOMBRE DEL ESTUDIANTE: Sánchez Plaza Edson

FECHA DE ENTREGA DE LA PRÁCTICA: 28/03/2025

Sucre- Bolivia
2024
1. Un líquido que contiene 60 moles % de tolueno y 40 moles % de benceno entra como alimentación continua a
una unidad de destilación de una sola etapa que opera a la presión atmosférica. ¿Qué porcentaje del benceno
contenido en la alimentación sale con el vapor si el 90% del tolueno que entra con la alimentación sale con el
líquido?
Supóngase una volatilidad relativa de 2,3 y obténgase la solución gráficamente.

Base de cálculo de 1 mol

𝑧𝐹𝑏 = 60 𝑚𝑜𝑙 𝑑𝑒 𝑡𝑜𝑙𝑢𝑒𝑛𝑜

𝑧𝐹𝑎 = 40 𝑚𝑜𝑙 𝑑𝑒 𝑏𝑒𝑛𝑐𝑒𝑛𝑜

𝐹 =𝐿+𝑉

𝐹𝑧𝐹𝑎 = 𝐿𝑥𝑎 + 𝑉𝑦𝑎

Conociendo los datos de alimentación podemos armar un sistema de dos ecuaciones con dos incógnitas:

𝐹 =𝐿+𝑉
{
𝐹𝑧𝐹𝑎 = 𝐿𝑥𝑎 + 𝑉𝑦𝑎

1=𝐿+𝑉
{
0,4 ∗ 1 = 𝐿 ∗ 𝑥𝑎 + 𝑉 ∗ 𝑦𝑎

1=𝐿+𝑉
{
0,4 = 𝐿 ∗ 𝑥𝑎 + 𝑉 ∗ 𝑦𝑎
Sabemos que el 90% de la alimentación de tolueno sale como líquido, por lo tanto, el restante 10% de la
alimentación es del benceno.

𝐿𝑥𝑎 = 0,1𝐹𝑧𝐹𝑎

𝐿𝑥𝑎 = 0,1 ∗ 1 ∗ 0,4

𝐿𝑥𝑎 = 0,04

0,04
𝑥𝑎 =
𝐿

𝑉𝑦𝑎 = 0,9𝐹𝑧𝐹𝑎

𝑉𝑦𝑎 = 0,9 ∗ 1 ∗ 0,4

𝑉𝑦𝑎 = 0,36

0,36
𝑦𝑎 =
𝑉
Con estas ecuaciones nos armamos el siguiente algoritmo para hallar la fracción molar en ambas fases.

Donde a través de cálculos iterativos en Excel se obtuvieron los siguientes resultados:

Xa Ya Xb Yb Ka Kb
0,246725 0,42966 0,753275 0,57034 1,74145 0,75715
Donde vemos que el benceno sale un 42,966% en el vapor.

Con estos resultados podemos calcular la fracción vaporizada.

𝑛
𝑧𝐹 (1 − 𝐾𝑖 )
∑ =0
1 + 𝑓(𝐾𝑖 − 1)
𝑖=1

𝑧𝐹𝑎 (1 − 𝐾𝑎 ) 𝑧𝐹𝑏 (1 − 𝐾𝑏 )
+ =0
1 + 𝑓(𝐾𝑎 − 1) 1 + 𝑓(𝐾𝑏 − 1)
0,4(1 − 1,74145) 0,6(1 − 0,75715)
+ =0
1 + 𝑓(1,74145 − 1) 1 + 𝑓(0,75715 − 1)
Resolviendo esta ecuación tenemos como resultado:

𝑓 = 0,855
Estos cálculos nos dan la siguiente gráfica:

Podemos observar que nuestros cálculos son correctos, ya que, al considerar una fracción vaporizada de 0,855, se
evidencia cómo este valor influye en las composiciones tanto de la fase vapor como de la fase líquida, tal y como se observa
en la figura.
2. Una mezcla líquida consistente en 100 kmol con 60 moles % de benceno, 25 moles % de tolueno y 15 moles % de
o-xileno se somete a una operación de flash a 1 atm y 100°C.
a) Calcúlense las cantidades de líquido y vapor que se forman, así como sus composiciones.
b) Repítanse los cálculos para 100°C y 2 atm.
c) Repítanse los cálculos para 105°C y 0,1 atm.
d) Repítanse los cálculos para 150°C y 1 atm.
Considérense soluciones ideales y utilícese la ecuación de Antoine.
Para el apartado a).
Con valores de constantes para la ecuación de Antoine en T [=] °C, P [=] KPa:
A B C
Benceno 6,03055 1211,033 220,790
Tolueno 6,07954 1344,800 219,482
o-xileno 6,12381 1474,679 213,686
Con la ecuación de Antoine:
𝐵
log 𝑃 = 𝐴 −
𝑇+𝐶
𝐵
𝑃 = 10𝐴−𝑇+𝐶
Para el benceno:
1211,033
6,03055−
100°𝐶+220,790
𝑃 = 10
𝑃 = 180,04945 𝐾𝑃𝑎
Para el tolueno:
1344,800
6,07954−
𝑃 = 10 100°𝐶+219,482

𝑃 = 74,16964 𝐾𝑃𝑎
Para el o-xileno:
1474,679
6,12381−
𝑃 = 10 100°𝐶+213,686

𝑃 = 26,46543 𝐾𝑃𝑎
Con la presión de operación y las presiones de saturación calculamos las constantes de equilibrio sabiendo que 1
atm es igual a 101,3 Kpa:
𝑃𝑖𝑠𝑎𝑡
𝐾𝑖 =
𝑃𝑡𝑜𝑡𝑎𝑙
Para el benceno:
180,04945 𝐾𝑃𝑎
𝐾𝑏 =
101,3 𝐾𝑃𝑎
𝐾𝑏 = 1,7773884
Para el tolueno:
74,16964 𝐾𝑃𝑎
𝐾𝑡 =
101,3 𝐾𝑃𝑎
𝐾𝑡 = 0,7321781
Para el o-xileno:
26,46543 𝐾𝑃𝑎
𝐾𝑜 =
101,3 𝐾𝑃𝑎
𝐾𝑜 = 0,2612579
Con la ecuación de Rachford-Rice:
𝑛
𝑧𝐹 (1 − 𝐾𝑖 )
∑ =0
1 + 𝑓(𝐾𝑖 − 1)
𝑖=1

𝑧𝐹𝑏 (1 − 𝐾𝑏 ) 𝑧𝐹𝑡 (1 − 𝐾𝑡 ) 𝑧𝐹𝑜 (1 − 𝐾𝑜 )

+ + =0
1 + 𝑓(𝐾𝑏 − 1) 1 + 𝑓(𝐾𝑡 − 1) 1 + 𝑓(𝐾𝑜 − 1)
0,6 ∗ (1 − 1,7773884) 0,25 ∗ (1 − 0,7321781) 0,15 ∗ (1 − 0,2612579)
+ + =0
1 + 𝑓(1,7773884 − 1) 1 + 𝑓(0,7321781 − 1) 1 + 𝑓(0,2612579 − 1)
La ecuación se vuelve de una ecuación lineal en una ecuación cuadrada, dando como resultado:
𝑓 = 0,681172
𝑓 = 2,755258
Donde tomamos como valor verdadero la primera raíz, puesto que, la fracción vaporizada no puede ser mayor que la
unidad ya que sabemos que no se puede vaporizar más de lo que se alimenta al separador flash.
Tomando como base de cálculo 100 moles de alimentación:
𝑉
𝑓=
𝐹
𝑉 = 𝐹∗𝑓
𝑉 = 100 𝑚𝑜𝑙 ∗ 0,681172
𝑉 = 68,1172 𝑚𝑜𝑙
𝐿=𝐹−𝑉
𝐿 = 100 𝑚𝑜𝑙 − 68,1172 𝑚𝑜𝑙
𝐿 = 31,8828 𝑚𝑜𝑙
Para calcular las fracciones molares en el líquido usaremos:
𝑧𝐹𝑖
𝑥𝑖 =
1 + 𝑓(𝐾𝑖 − 1)
Para el benceno:
𝑧𝐹𝑏
𝑥𝑏 =
1 + 𝑓(𝐾𝑏 − 1)
0,6
𝑥𝑏 =
1 + 0,681172 ∗ (1,7773884 − 1)
𝑥𝑏 = 0,392276
Para el tolueno:
𝑧𝐹𝑡
𝑥𝑡 =
1 + 𝑓(𝐾𝑡 − 1)
0,25
𝑥𝑡 =
1 + 0,681172 ∗ (0,7321781 − 1)
𝑥𝑡 = 0,305785
Para el o-xileno:
𝑧𝐹𝑜
𝑥𝑜 =
1 + 𝑓(𝐾𝑜 − 1)
0,15
𝑥𝑜 =
1 + 0,681172 ∗ (0,2612579 − 1)
𝑥𝑜 = 0,301939
Para calcular las fracciones molares en el vapor usaremos:
𝑦𝑖 = 𝐾𝑖 ∗ 𝑥𝑖
Par el benceno:
𝑦𝑏 = 𝐾𝑏 ∗ 𝑥𝑏
𝑦𝑏 = 1,7773884 ∗ 0,392276
𝑦𝑏 = 0,697227
Para el tolueno:
𝑦𝑡 = 𝐾𝑡 ∗ 𝑥𝑡
𝑦𝑡 = 0,7321781 ∗ 0,305785
𝑦𝑡 = 0,223889
Para el o-xileno:
𝑦𝑜 = 𝐾𝑜 ∗ 𝑥𝑜
𝑦𝑜 = 0,2612579 ∗ 0,301939
𝑦𝑜 = 0,078884
Para el inciso b).
Como trabaja a la misma temperatura, ya tenemos calculadas las presiones de saturación.
Para el benceno:
𝑃 = 180,04945 𝐾𝑃𝑎
Para el tolueno:
𝑃 = 74,16964 𝐾𝑃𝑎
Para el o-xileno:
𝑃 = 26,46543 𝐾𝑃𝑎
Con la presión de operación y las presiones de saturación calculamos las constantes de equilibrio sabiendo que 2 atm es
igual a 202,6 Kpa:
𝑃𝑖𝑠𝑎𝑡
𝐾𝑖 =
𝑃𝑡𝑜𝑡𝑎𝑙
Para el benceno:
180,04945 𝐾𝑃𝑎
𝐾𝑏 =
202,6 𝐾𝑃𝑎
𝐾𝑏 = 0,8886942
Para el tolueno:
74,16964 𝐾𝑃𝑎
𝐾𝑡 =
202,6 𝐾𝑃𝑎
𝐾𝑡 = 0,3660890
Para el o-xileno:
26,46543 𝐾𝑃𝑎
𝐾𝑜 =
202,6 𝐾𝑃𝑎
𝐾𝑜 = 0,1306290
Con la ecuación de Rachford-Rice:
𝑛
𝑧𝐹 (1 − 𝐾𝑖 )
∑ =0
1 + 𝑓(𝐾𝑖 − 1)
𝑖=1
 𝑧𝐹𝑏 (1 − 𝐾𝑏 ) 𝑧𝐹𝑡 (1 − 𝐾𝑡 ) 𝑧𝐹𝑜 (1 − 𝐾𝑜 )
+ + =0
1 + 𝑓(𝐾𝑏 − 1) 1 + 𝑓(𝐾𝑡 − 1) 1 + 𝑓(𝐾𝑜 − 1)
0,6 ∗ (1 − 0,8886942) 0,25 ∗ (1 − 0,3660890) 0,15 ∗ (1 − 0,1306290)
+ + =0
1 + 𝑓(0,8886942 − 1) 1 + 𝑓(0,3660890 − 1) 1 + 𝑓(0,1306290 − 1)
La ecuación se vuelve de una ecuación lineal en una ecuación cuadrada, dando como resultado:
𝑓 = 1,302306
𝑓 = 4,452247
Como vemos, ambas raíces son mayores a la unidad, lo que quiere decir que; la ecuación de Rachford-Rice no tiene
solución física en el intervalo 0 ≤ f ≤ 1. Esto indica que la mezcla permanece completamente en fase líquida a 100°C y 2
atm. Por lo tanto, la composición de la fase líquida será la misma que la alimentación.
Para el incido c).
Con la ecuación de Antoine:
𝐵
log 𝑃 = 𝐴 −
𝑇+𝐶
𝐵
𝑃 = 10𝐴−𝑇+𝐶
Para el benceno:
1211,033
6,03055−
𝑃 = 10 105°𝐶+220,790

𝑃 = 205,74549 𝐾𝑃𝑎
Para el tolueno:
1344,800
6,07954−
105°𝐶+219,482
𝑃 = 10
𝑃 = 86,11687 𝐾𝑃𝑎
Para el o-xileno:
1474,679
6,12381−
𝑃 = 10 105°𝐶+213,686

𝑃 = 31,36440 𝐾𝑃𝑎
Con la presión de operación y las presiones de saturación calculamos las constantes de equilibrio sabiendo que 0,1 atm es
igual a 10,13 Kpa:
𝑃𝑖𝑠𝑎𝑡
𝐾𝑖 =
𝑃𝑡𝑜𝑡𝑎𝑙
Para el benceno:
205,74549 𝐾𝑃𝑎
𝐾𝑏 =
10,13 𝐾𝑃𝑎
𝐾𝑏 = 20,3105121
Para el tolueno:
86,11687 𝐾𝑃𝑎
𝐾𝑡 =
10,13 𝐾𝑃𝑎
𝐾𝑡 = 8,5011716
Para el o-xileno:
31,36400 𝐾𝑃𝑎
𝐾𝑜 =
10,13 𝐾𝑃𝑎
𝐾𝑜 = 3,0961895
Con la ecuación de Rachford-Rice:
𝑛
𝑧𝐹 (1 − 𝐾𝑖 )
∑ =0
1 + 𝑓(𝐾𝑖 − 1)
𝑖=1

𝑧𝐹𝑏 (1 − 𝐾𝑏 ) 𝑧𝐹𝑡 (1 − 𝐾𝑡 ) 𝑧𝐹𝑜 (1 − 𝐾𝑜 )

+ + =0
1 + 𝑓(𝐾𝑏 − 1) 1 + 𝑓(𝐾𝑡 − 1) 1 + 𝑓(𝐾𝑜 − 1)

0,6 ∗ (1 − 20,3105121) 0,25 ∗ (1 − 8,5011716) 0,15 ∗ (1 − 3,0961895)

+ + =0
1 + 𝑓(20,3105121 − 1) 1 + 𝑓(8,5011716 − 1) 1 + 𝑓(3,0961895 − 1)
La ecuación se vuelve de una ecuación lineal en una ecuación cuadrada, dando como resultado:
𝑓 = −0,417534
𝑓 = −0,108662
En este caso, como la presión de operación es muy baja (0,1 atm) y la temperatura es relativamente alta (105°C), la mezcla
estará por completo en la fase de vapor. No habrá fase líquida en equilibrio, por lo que el flujo de líquido L=0 y el flujo de
vapor V=F. Por lo tanto, toda la alimentación se evapora y la composición de la fase de vapor será la misma que la
alimentación.

Para el inciso d).

Con la ecuación de Antoine:
𝐵
log 𝑃 = 𝐴 −
𝑇+𝐶
𝐵
𝑃 = 10𝐴−𝑇+𝐶
Para el benceno:
1211,033
6,03055−150°𝐶+220,790
𝑃 = 10
𝑃 = 581,38173 𝐾𝑃𝑎
Para el tolueno:
1344,800
6,07954−
𝑃 = 10 150°𝐶+219,482

𝑃 = 275,32766 𝐾𝑃𝑎
Para el o-xileno:
1474,679
6,12381−150°𝐶+213,686
𝑃 = 10
𝑃 = 117,21852 𝐾𝑃𝑎
Con la presión de operación y las presiones de saturación calculamos las constantes de equilibrio sabiendo que 1 atm es
igual a 101,3 Kpa:
𝑃𝑖𝑠𝑎𝑡
𝐾𝑖 =
𝑃𝑡𝑜𝑡𝑎𝑙

Para el benceno:

581,38173 𝐾𝑃𝑎
𝐾𝑏 =
101,3 𝐾𝑃𝑎

𝐾𝑏 = 5,7392076
Para el tolueno:

275,32766 𝐾𝑃𝑎
𝐾𝑡 =
101,3 𝐾𝑃𝑎

𝐾𝑡 = 2,7179433

Para el o-xileno:
117,21852 𝐾𝑃𝑎
𝐾𝑜 =
101,3 𝐾𝑃𝑎
𝐾𝑜 = 1,1571423
Con la ecuación de Rachford-Rice:
𝑛
𝑧𝐹 (1 − 𝐾𝑖 )
∑ =0
1 + 𝑓(𝐾𝑖 − 1)
𝑖=1

𝑧𝐹𝑏 (1 − 𝐾𝑏 ) 𝑧𝐹𝑡 (1 − 𝐾𝑡 ) 𝑧𝐹𝑜 (1 − 𝐾𝑜 )

+ + =0
1 + 𝑓(𝐾𝑏 − 1) 1 + 𝑓(𝐾𝑡 − 1) 1 + 𝑓(𝐾𝑜 − 1)
0,6 ∗ (1 − 5,7392076) 0,25 ∗ (1 − 2,7179433) 0,15 ∗ (1 − 1,1571423)
+ + =0
1 + 𝑓(5,7392076 − 1) 1 + 𝑓(2,7179433 − 1) 1 + 𝑓(1,1571423 − 1)
La ecuación se vuelve de una ecuación lineal en una ecuación cuadrada, dando como resultado:
𝑓 = −5,457993
𝑓 = −0,472088
En este caso, como la presión de operación es relativamente baja (1 atm) con respecto a la temperatura (150°C), la mezcla
estará por completo en la fase de vapor. No habrá fase líquida en equilibrio, por lo que el flujo de líquido L=0 y el flujo de
vapor V=F. Por lo tanto, toda la alimentación se evapora y la composición de la fase de vapor será la misma que la
alimentación.
 3. La mezcla que se muestra a continuación se condensa parcialmente y se separa en dos fases. Calcúlense las
cantidades y composiciones de las fases en equilibrio, V y L.

Kmol/h
H2 72,53
Ni 7,98
Benceno 0,13
Tolueno 150,00

Primero precedemos a calcular la fracción molar de cada compuesto en la alimentación, primero calculamos el flujo
total:
𝐹𝑇 = ∑ 𝐹𝑖
𝐹𝑇 = 𝐹𝐻2 + 𝐹𝑁𝑖 + 𝐹𝑏 + 𝐹𝐶6
𝐹𝑇 = 72,53 𝑘𝑚𝑜𝑙 + 7,98 𝑘𝑚𝑜𝑙 + 0,13 𝑘𝑚𝑜𝑙 + 150,0 𝑘𝑚𝑜𝑙
𝐹𝑇 = 230,64 𝑘𝑚𝑜𝑙
Ya con el flujo total calculamos las fracciones molares:
𝐹𝑖
𝑧𝑖 =
𝐹𝑇
Para el hidrogeno:
𝐹𝐻2
𝑧𝐻2 =
𝐹𝑇
72,53 𝑘𝑚𝑜𝑙
𝑧𝐻2 =
230,64 𝑘𝑚𝑜𝑙
𝑧𝐻2 = 0,314473
Para el níquel:
𝐹𝑁𝑖
𝑧𝑁𝑖 =
𝐹𝑇
7,98 𝑘𝑚𝑜𝑙
𝑧𝑁𝑖 =
230,64 𝑘𝑚𝑜𝑙
𝑧𝑁𝑖 = 0,034599
Para el benceno:
𝐹𝑏
𝑧𝑏 =
𝐹𝑇
0,13 𝑘𝑚𝑜𝑙
𝑧𝑏 =
230,64 𝑘𝑚𝑜𝑙
𝑧𝑏 = 0,000564
Para el ciclo hexano:
𝐹𝐶6
𝑧𝐶6 =
𝐹𝑇
150 𝑘𝑚𝑜𝑙
𝑧𝐶6 =
230,64 𝑘𝑚𝑜𝑙
𝑧𝐶6 = 0,650364
Convirtiendo los datos al sistema internacional de unidades tenemos que P = 2068,43 Kpa y que T = 48,89°C a la entrada
del flash.
Con valores de constantes para la ecuación de Antoine en T [=] °C, P [=] KPa:
A B C
H2 4,24710 2556,400 -67,0
Níquel 4,71570 2233,200 -38,4
Benceno 6,03055 1211,033 220,79
Ciclo hexano 5,96620 1201,531 222,647
Con la ecuación de Antoine:
𝐵
log 𝑃 = 𝐴 −
𝑇+𝐶
𝐵
𝑃 = 10𝐴−𝑇+𝐶
Para el hidrogeno:
2556,400
4,24710−
𝑃 = 10 48,89°𝐶−67,0

𝑃 = 2,55𝑥10145 𝐾𝑃𝑎
Para el níquel:
2233,200
4,71570−
𝑃 = 10 48,89°𝐶−38,4

𝑃 = 6,72𝑥10−209 𝐾𝑃𝑎
Para el benceno:
1211,033
6,0305−
𝑃 = 10 48,89°𝐶+220,79

𝑃 = 34,667328 𝐾𝑃𝑎
Para el ciclo hexano:
1201,531
5,96620−
𝑃 = 10 48,89°𝐶+222,647

𝑃 = 34,775588 𝐾𝑃𝑎
Calculamos las constantes de equilibrio:
𝑃𝑖𝑠𝑎𝑡
𝐾𝑖 =
𝑃𝑡𝑜𝑡𝑎𝑙

Para el benceno:

34,667328 𝐾𝑃𝑎
𝐾𝑏 =
2068,43 𝐾𝑃𝑎
𝐾𝑏 = 0,016760

Para el ciclo hexano:

34,775588 𝐾𝑃𝑎
𝐾𝐶6 =
2068,43 𝐾𝑃𝑎
𝐾𝐶6 = 0,016813
Para el hidrógeno, se ha tomado un valor de Ki=50, basado en datos obtenidos de la bibliografía, ya que los resultados
indican que este elemento tiende a permanecer en fase gaseosa. En cuanto al níquel, se asumió un valor de Ki=0, debido
a que, al ser un metal, requiere condiciones extremas para pasar a fase líquida y, con mayor razón, para estar en fase
gaseosa.
Ahora procedemos a calcular la fracción vaporizada con la ecuación de Rachford-Rice:
𝑧𝐹 (1 − 𝐾𝑖 )
∑ =0
1 + 𝑓(𝐾𝑖 − 1)
𝑧𝐹𝐻2 (1 − 𝐾𝐻2 ) 𝑧𝐹𝑁𝑖 (1 − 𝐾𝑁𝑖 ) 𝑧𝐹𝑏 (1 − 𝐾𝑏 ) 𝑧𝐹𝐶6 (1 − 𝐾𝐶6 )
+ + + =0
1 + 𝑓(𝐾𝐻2 − 1) 1 + 𝑓(𝐾𝑁𝑖 − 1) 1 + 𝑓(𝐾𝑏 − 1) 1 + 𝑓(𝐾𝐶6 − 1)
0,314473 ∗ (1 − 50) 0,034599 ∗ (1 − 0) 0,000564 ∗ (1 − 0,016760) 0,650364 ∗ (1 − 0,016813)
+ + + =0
1 + 𝑓(50 − 1) 1 + 𝑓(0 − 1) 1 + 𝑓(0,016760 − 1) 1 + 𝑓(0,016813 − 1)

𝑓 = 0,306
𝑓 = 1,017
𝑓 = 1,001
Al resolver la ecuación se obtuvieron tres raíces, de las cuales tomaremos la primera para hallar las fracciones molares
en cada una de las fases.

Para la fase líquida:

𝑧𝐹
𝑥𝑖 =
1 + 𝑓(𝐾𝑖 − 1)
Para el hidrogeno:
0,314473
𝑥𝐻2 =
1 + 0,306 ∗ (50 − 1)
𝑥𝐻2 = 0,019662
Para el níquel:
0,034599
𝑥𝑁𝑖 =
1 + 0,306 ∗ (0 − 1)
𝑥𝑁𝑖 = 0,049825
Para el benceno:
0,000564
𝑥𝑏 =
1 + 0,306 ∗ (0,016760 − 1)
𝑥𝑏 = 0,000305
Para el ciclo hexano:
0,650364
𝑥𝐶6 =
1 + 0,306 ∗ (0,016813 − 1)
𝑥𝐶6 = 0,930208
Para calcular las fracciones molares en fase gaseosa usaremos:
𝑦𝑖
𝐾𝑖 =
𝑥𝑖
𝑦𝑖 = 𝐾𝑖 ∗ 𝑥𝑖
Para el hidrogeno:
𝑦𝐻2 = 50 ∗ 0,019662
𝑦𝐻2 = 0,983096
Para el níquel:
𝑦𝑁𝑖 = 0 ∗ 0,049825
𝑦𝐻2 = 0,0
Para el benceno:
𝑦𝑏 = 0,016760 ∗ 0,0000305
𝑦𝑏 = 0,000914
Para el ciclo hexano:
𝑦𝐶6 = 0,016813 ∗ 0,930208
𝑦𝐶6 = 0,015990

El hidrógeno está casi completamente en la fase vapor, lo cual es consistente con su alta volatilidad. El níquel permanece
en la fase líquida, ya que no es volátil. El benceno y el ciclohexano están mayormente en la fase líquida, aunque una
pequeña cantidad de ciclohexano pasa a la fase vapor.

Para conocer que cantidad de vapor y de líquido salen, se resuelve el siguiente sistema de ecuaciones:
𝐹 =𝐿+𝑉
{
𝐹𝑧𝐻2 = 𝐿𝑥𝐻2 + 𝑉𝑦𝐻2

230,64 = 𝐿 + 𝑉
{
230,64 ∗ 0,314473 = 𝐿 ∗ 0,019662 + 𝑉 ∗ 0,983096

230,64 = 𝐿 + 𝑉
{
72,530053 = 𝐿 ∗ 0,019662 + 𝑉 ∗ 0,983096

𝐿 = 160,0641 𝑘𝑚𝑜𝑙/ℎ
{
𝑉 = 70,5759 𝑘𝑚𝑜𝑙/ℎ
 4. La siguiente mezcla se introduce en una columna de destilación como líquido saturado a 1,72 MPa. Calcúlese la
temperatura del punto de burbuja utilizando los valores K.

Kmol/h
Etano 1,5
Propano 10,0
n-butano 18,5
n-pentano 17,5
n-hexano 3,5

Primero precedemos a calcular la fracción molar de cada compuesto en la alimentación, primero calculamos el flujo
total:
𝐹𝑇 = ∑ 𝐹𝑖
𝐹𝑇 = 𝐹𝐶2 + 𝐹𝐶3 + 𝐹𝐶4 + 𝐹𝐶5 + 𝐹𝐶6
𝐹𝑇 = 1,5 𝑚𝑜𝑙 + 10 𝑚𝑜𝑙 + 18,5 𝑚𝑜𝑙 + 17,5 𝑚𝑜𝑙 + 3,5 𝑚𝑜𝑙
𝐹𝑇 = 51 𝑚𝑜𝑙
Ya con el flujo total calculamos las fracciones molares:
𝐹𝑖
𝑧𝑖 =
𝐹𝑇
Para el etano:
𝐹𝐶2
𝑧𝐶2 =
𝐹𝑇
1,5 𝑚𝑜𝑙
𝑧𝐶2 =
51 𝑚𝑜𝑙
𝑧𝐶2 = 0,02941
Para el propano:
𝐹𝐶3
𝑧𝐶3 =
𝐹𝑇
10 𝑚𝑜𝑙
𝑧𝐶3 =
51 𝑚𝑜𝑙
𝑧𝐶3 = 0,19607
Para el n-butano:
𝐹𝐶4
𝑧𝐶4 =
𝐹𝑇
18,5 𝑚𝑜𝑙
𝑧𝐶4 =
51 𝑚𝑜𝑙
𝑧𝐶4 = 0,36275
Para el n-pentano:
𝐹𝐶5
𝑧𝐶5 =
𝐹𝑇
17,5 𝑚𝑜𝑙
𝑧𝐶5 =
51 𝑚𝑜𝑙
𝑧𝐶5 = 0,34314
Para el n-hexano:
𝐹𝐶6
𝑧𝐶6 =
𝐹𝑇
3,5 𝑚𝑜𝑙
𝑧𝐶6 =
51 𝑚𝑜𝑙
𝑧𝐶6 = 0,06863
Ya que el problema nos dice que la alimentación es líquido saturado, asumimos que f=0 por lo que z i = xi.

Conociendo esta información: nos podemos plantear el siguiente algoritmo de cálculo para hallar la temperatura del punto
de burbuja:
Por medio de cálculos iterativos en Excel, se obtuvieron los siguientes resultados:

Suponiendo una temperatura de:

𝑇𝑖𝑠𝑎𝑡 = 93,16755°𝐶

Con la ecuación de Antoine y sus constantes en T [=] °C, P [=] KPa:

A B C
Etano 5,95942 663,700 256,470
Propano 5,92888 803,810 246,990
n-butano 5,93386 935,860 238,730
n-pentano 6,00122 1075,780 233,205
n-hexano 5,99514 1168,720 224,210
𝐵
log 𝑃 = 𝐴 −
𝑇+𝐶
𝐵
𝑃 = 10𝐴−𝑇+𝐶
Para el etano:
663,700
5,95942−
𝑃 = 10 93,16755°𝐶+256,470

𝑃 = 11512,469979 𝐾𝑃𝑎
Para el propano:
803,810
5,92888−
𝑃 = 10 93,16755°𝐶+246,990

𝑃 = 3679,831543 𝐾𝑃𝑎
Para el n-butano:
935,860
5,93386−
𝑃 = 10 93,16755°𝐶+238,730

𝑃 = 1300,572166 𝐾𝑃𝑎
Para el n-pentano:
1075,780
6,00122−
𝑃 = 10 93,16755°𝐶+238,205

𝑃 = 507,046964 𝐾𝑃𝑎
Para el n-hexano:
1168,720
5,99514−
𝑃 = 10 93,16755°𝐶+224,210

𝑃 = 205,452735 𝐾𝑃𝑎
Con la presión de operación y las presiones de saturación calculamos las constantes de equilibrio sabiendo que 1,72 MPa
es igual a 1720 Kpa:
𝑃𝑖𝑠𝑎𝑡
𝐾𝑖 =
𝑃𝑡𝑜𝑡𝑎𝑙

Para el etano:

11512,469979 𝐾𝑃𝑎
𝐾𝐶2 =
1720 𝐾𝑃𝑎
𝐾𝐶2 = 6,693296

Para el propano:

3679,831543 𝐾𝑃𝑎
𝐾𝐶3 =
1720 𝐾𝑃𝑎
𝐾𝐶3 = 2,139437
Para el n-butano:
1300,572166 𝐾𝑃𝑎
𝐾𝐶4 =
1720 𝐾𝑃𝑎
𝐾𝐶4 = 0,756147
Para el n-pentano:
507,046964 𝐾𝑃𝑎
𝐾𝐶5 =
1720 𝐾𝑃𝑎
𝐾𝐶5 = 0,294795
Para el n-hexano:
205,452735 𝐾𝑃𝑎
𝐾𝐶6 =
1720 𝐾𝑃𝑎
𝐾𝐶3 = 0,119449
Ahora calculamos yi sabiendo que f=0 por lo que zi = xi:
𝑦𝑖
𝐾𝑖 =
𝑥𝑖
𝑦𝑖 = 𝐾𝑖 ∗ 𝑥𝑖
𝑦𝑖 = 𝐾𝑖 ∗ 𝑧𝑖
Para el etano:
𝑦𝐶2 = 6,693296 ∗ 0,02941
𝑦𝐶2 = 0,196862
Para el propano:
𝑦𝐶3 = 2,139437 ∗ 0,19607
𝑦𝐶3 = 0,419498
Para el n-butano:
𝑦𝐶4 = 0,756147 ∗ 0,36275
𝑦𝐶4 = 0,274288
Para el n-pentano:
𝑦𝐶5 = 0,294795 ∗ 0,34314
𝑦𝐶5 = 0,101155
Para el n-hexano:
𝑦𝐶6 = 0,119449 ∗ 0,06863
𝑦𝐶2 = 0,008197
Para comprobar que los cálculos están bien hechos se verifica con la siguiente sumatoria:

∑ 𝑦𝑖 = 1

0,196862 + 0,419498 + 0,274288 + 0,101155 + 0,008197 = 1

 5. Para una mezcla que consta de 45 moles % de n-hexano, 25 moles % de n-heptano y 30 moles % de n-octano a 1
atm,
a) Hallar las temperaturas del punto de burbuja y del punto de rocío.
b) La mezcla se somete a una destilación de flash a 1 atm de forma que se vaporiza el 50% de la alimentación.
Hallar la temperatura de flash, así como la composición y las cantidades relativas de los productos líquido y
vapor.
c) Repítanse los apartados (a) y (b) a 5 atm y 0,5 atm.
d) Si se retira el 90% del hexano como vapor, ¿cuánto octano se retira como vapor?
Para el inciso a).
Para calcular la temperatura del punto de burbuja, nos planteamos el siguiente algoritmo de cálculo:

Por medio de cálculos iterativos en Excel, se obtuvieron los siguientes resultados:

Suponiendo una temperatura de:

𝑇𝑖𝑠𝑎𝑡 = 85,73805°𝐶

Con la ecuación de Antoine y sus constantes en T [=] °C, P [=] KPa:

A B C
n-hexano 5,99514 1168,72 224,210
n-heptano 6,01875 1264,37 216,636
n-octano 6,03430 1349,82 209,385
𝐵
log 𝑃 = 𝐴 −
𝑇+𝐶
𝐵
𝑃 = 10𝐴−𝑇+𝐶
Para el n-hexano:
1168,72
5,99514−
𝑃 = 10 85,73805°𝐶+224,210

𝑃 = 167,665454 𝐾𝑃𝑎
Para el n-heptano:
1264,37
6,01875−
𝑃 = 10 85,73805°𝐶+216,636

𝑃 = 68,750118 𝐾𝑃𝑎
Para el n-octano:
1349,82
6,0343−
𝑃 = 10 85,73805°𝐶+209,385

𝑃 = 28,876648 𝐾𝑃𝑎
Con la presión de operación y las presiones de saturación calculamos las constantes de equilibrio sabiendo que 1 atm es
igual a 101,3 Kpa:
𝑃𝑖𝑠𝑎𝑡
𝐾𝑖 =
𝑃𝑡𝑜𝑡𝑎𝑙

Para el n-hexano:

167,665454 𝐾𝑃𝑎
𝐾𝐶6 =
101,3 𝐾𝑃𝑎
𝐾𝐶6 = 1,655138

Para el n-heptano:

68,750118 𝐾𝑃𝑎
𝐾𝐶7 =
101,3 𝐾𝑃𝑎
𝐾𝐶7 = 0,678678

Para el n-octano:
28,876648 𝐾𝑃𝑎
𝐾𝐶8 =
101,3 𝐾𝑃𝑎
𝐾𝐶8 = 0,285061
Ahora calculamos yi sabiendo que f=0 por lo que zi = xi:
𝑦𝑖
𝐾𝑖 =
𝑥𝑖
𝑦𝑖 = 𝐾𝑖 ∗ 𝑥𝑖
𝑦𝑖 = 𝐾𝑖 ∗ 𝑧𝑖
Para el n-hexano:
𝑦𝐶6 = 1,655138 ∗ 0,45
𝑦𝐶6 = 0,744812
Para el n-heptano:
𝑦𝐶7 = 0,678678 ∗ 0,25
𝑦𝐶7 = 0,169670
Para el n-octano:
𝑦𝐶8 = 0,285061 ∗ 0,3
𝑦𝐶8 = 0,085518
Para comprobar que los cálculos están bien hechos se verifica con la siguiente sumatoria:

∑ 𝑦𝑖 = 1

0,744812 + 0,169670 + 0,085518 = 1

Para calcular la temperatura del punto de rocío, nos planteamos el siguiente algoritmo de cálculo:

Por medio de cálculos iterativos en Excel, se obtuvieron los siguientes resultados:

Suponiendo una temperatura de:

𝑇𝑖𝑠𝑎𝑡 = 102,29461°𝐶
Con la ecuación de Antoine y sus constantes en T [=] °C, P [=] KPa:

A B C
n-hexano 5,99514 1168,72 224,210
n-heptano 6,01875 1264,37 216,636
n-octano 6,03430 1349,82 209,385
𝐵
log 𝑃 = 𝐴 −
𝑇+𝐶
𝐵
𝑃 = 10𝐴−𝑇+𝐶
Para el n-hexano:
1168,72
5,99514−
𝑃 = 10 102,29461°𝐶+224,210

𝑃 = 260,405384 𝐾𝑃𝑎
Para el n-heptano:
1264,37
6,01875−
𝑃 = 10 102,29461°𝐶+216,636

𝑃 = 113,330115 𝐾𝑃𝑎
Para el n-octano:
1349,82
6,0343−
102,29461°𝐶+209,385
𝑃 = 10
𝑃 = 50,525032 𝐾𝑃𝑎
Con la presión de operación y las presiones de saturación calculamos las constantes de equilibrio sabiendo que 1 atm es
igual a 101,3 Kpa:
𝑃𝑖𝑠𝑎𝑡
𝐾𝑖 =
𝑃𝑡𝑜𝑡𝑎𝑙

Para el n-hexano:

260,405384 𝐾𝑃𝑎
𝐾𝐶6 =
101,3 𝐾𝑃𝑎
𝐾𝐶6 = 2,570636

Para el n-heptano:

113,330115 𝐾𝑃𝑎
𝐾𝐶7 =
101,3 𝐾𝑃𝑎
𝐾𝐶7 = 1,118757

Para el n-octano:
50,525032 𝐾𝑃𝑎
𝐾𝐶8 =
101,3 𝐾𝑃𝑎
𝐾𝐶8 = 0,498766
Ahora calculamos yi sabiendo que f=1 por lo que zi = yi:
𝑦𝑖
𝐾𝑖 =
𝑥𝑖
𝑦𝑖
𝑥𝑖 =
𝐾𝑖
𝑧𝑖
𝑥𝑖 =
𝐾𝑖
Para el n-hexano:

0,45
𝑥𝐶6 =
2,570636
𝑥𝐶6 = 0,175054

Para el n-heptano:

0,25
𝑥𝐶7 =
1,118757
𝑥𝐶7 = 0,223462

Para el n-octano:

0,3
𝑥𝐶8 =
0,498766
𝑥𝐶8 = 0,601484

Para comprobar que los cálculos están bien hechos se verifica con la siguiente sumatoria:

∑ 𝑥𝑖 = 1

0,175054 + 0,223462 + 0,601484 = 1

Para el inciso b).

Para calcular la temperatura del flash, nos planteamos el siguiente algoritmo de cálculo:
Por medio de cálculos iterativos en Excel, se obtuvieron los siguientes resultados:

Suponiendo una temperatura de:

𝑇𝑖𝑠𝑎𝑡 = 93,54761°𝐶

Con la ecuación de Antoine y sus constantes en T [=] °C, P [=] KPa:

A B C
n-hexano 5,99514 1168,72 224,210
n-heptano 6,01875 1264,37 216,636
n-octano 6,03430 1349,82 209,385
𝐵
log 𝑃 = 𝐴 −
𝑇+𝐶
𝐵
𝑃 = 10𝐴−𝑇+𝐶
Para el n-hexano:
1168,72
5,99514−
𝑃 = 10 93,54761°𝐶+224,210

𝑃 = 207,546954 𝐾𝑃𝑎
Para el n-heptano:
1264,37
6,01875−
𝑃 = 10 93,54761°𝐶+216,636

𝑃 = 87,609538 𝐾𝑃𝑎
Para el n-octano:
1349,82
6,0343−
𝑃 = 10 93,54761°𝐶+209,385

𝑃 = 37,884145 𝐾𝑃𝑎
Con la presión de operación y las presiones de saturación calculamos las constantes de equilibrio sabiendo que 1 atm es
igual a 101,3 Kpa:
𝑃𝑖𝑠𝑎𝑡
𝐾𝑖 =
𝑃𝑡𝑜𝑡𝑎𝑙

Para el n-hexano:

207,546954 𝐾𝑃𝑎
𝐾𝐶6 =
101,3 𝐾𝑃𝑎
𝐾𝐶6 = 2,048835

Para el n-heptano:

87,609538 𝐾𝑃𝑎
𝐾𝐶7 =
101,3 𝐾𝑃𝑎
𝐾𝐶7 = 0,864852

Para el n-octano:
37,884145 𝐾𝑃𝑎
𝐾𝐶8 =
101,3 𝐾𝑃𝑎
𝐾𝐶8 = 0,373980
Sabiendo que f=0,5, con la ecuación de Rachford-Rice comprobamos si:
𝑧𝐹 (1 − 𝐾𝑖 )
∑ =0
1 + 𝑓(𝐾𝑖 − 1)
 𝑧𝐹𝐶6 (1 − 𝐾𝐶6 ) 𝑧𝐹𝐶7 (1 − 𝐾𝐶7 ) 𝑧𝐹𝐶8 (1 − 𝐾𝐶8 )
+ + =0
1 + 𝑓(𝐾𝐶6 − 1) 1 + 𝑓(𝐾𝐶7 − 1) 1 + 𝑓(𝐾𝐶8 − 1)
0,45 ∗ (1 − 2,048835) 0,25 ∗ (1 − 0,864852) 0,3 ∗ (1 − 0,373980)
+ + =0
1 + 0,5 ∗ (2,048835 − 1) 1 + 0,5 ∗ (0,864852 − 1) 1 + 0,5 ∗ (0,373980 − 1)
0=0
Como vemos que se cumple esta igualdad, ahora pasamos a calcular las fracciones parciales en ambas fases con:
𝑧𝐹 𝑧𝐹 𝐾𝑖
𝑥𝑖 = ; 𝑦𝑖 =
1 + 𝑓(𝐾𝑖 − 1) 1 + 𝑓(𝐾𝑖 − 1)
Para la fase líquida:
Para el n-hexano:
0,45
𝑥𝐶6 =
1 + 0,5 ∗ (2,048835 − 1)

𝑥𝐶6 = 0,295195
Para el n-heptano:
0,25
𝑥𝐶7 =
1 + 0,5 ∗ (0,864852 − 1)
𝑥𝐶7 = 0,268118
Para el n-octano:
0,3
𝑥𝐶8 =
1 + 0,5 ∗ (0,373980 − 1)

𝑥𝐶8 = 0,436687
Para verificar si están bien los cálculos:

∑ 𝑥𝑖 = 1

0,295195 + 0,268118 + 0,436687 = 1

Para la fase vapor:

Para el n-hexano:
0,45 ∗ 2,048835
𝑦𝐶6 =
1 + 0,5 ∗ (2,048835 − 1)

𝑦𝐶6 = 0,604805
Para el n-heptano:
0,25 ∗ 0,864852
𝑦𝐶7 =
1 + 0,5 ∗ (0,864852 − 1)

𝑦𝐶7 = 0,231882
Para el n-octano:
0,3 ∗ 0,373980
𝑦𝐶8 =
1 + 0,5 ∗ (0,373980 − 1)

𝑦𝐶8 = 0,163313
Para verificar si están bien los cálculos:

∑ 𝑥𝑖 = 1

0,604805 + 0,231882 + 0,163313 = 1

Para conocer las cantidades relativas de líquido y vapor, tomamos el dato de la fracción vaporizada: f=0,5, que nos dice
que el 50% de la alimentación sale como vapor por lo que el líquido representa el otro 50%.

Para el inciso c).

Se repiten todos los cálculos a 0,5 atm y 5 atm, que nos dejan la siguiente tabla de resultados.

Para 5 atm:

n-hexano n-heptano n-octano

𝑃 = 5 𝑎𝑡𝑚
𝑝𝑢𝑛𝑡𝑜 𝑑𝑒 𝑏𝑢𝑟𝑏𝑢𝑗𝑎; 𝑇𝑠𝑢𝑝 = 152,15541°𝐶
𝑃𝑖𝑠𝑎𝑡
775,996831 389,345362 199,882807
𝐾𝑖
1,532077 0,768698 0,394635
𝑦𝑖
0,689434 0,192174 0,118392
𝑝𝑢𝑛𝑡𝑜 𝑑𝑒 𝑟𝑜𝑐í𝑜; 𝑇𝑠𝑢𝑝 = 167,25084°𝐶
𝑃𝑖𝑠𝑎𝑡
1022,362943 531,066869 282,106941
𝐾𝑖
2,018486 1,048503 0,556973
𝑥𝑖
0,222939 0,238435 0,538626
𝑓𝑙𝑎𝑠ℎ; 𝑇𝑠𝑢𝑝 = 159,36105 °𝐶
𝑃𝑖𝑠𝑎𝑡
887,553868 452,936632 236,446002
𝐾𝑖
1,752327 0,894248 0,466823
𝑥𝑖
0,326996 0,263957 0,409047
𝑦𝑖
0,573004 0,236043 0,190953
Para 0,5 atm:

n-hexano n-heptano n-octano

𝑃 = 0,5 𝑎𝑡𝑚
𝑝𝑢𝑛𝑡𝑜 𝑑𝑒 𝑏𝑢𝑟𝑏𝑢𝑗𝑎; 𝑇𝑠𝑢𝑝 = 63,7387°𝐶
𝑃𝑖𝑠𝑎𝑡
86,369713 32,298149 12,363566
𝐾𝑖
1,705226 0,637673 0,244098
𝑦𝑖
0,767353 0,159418 0,073229
𝑝𝑢𝑛𝑡𝑜 𝑑𝑒 𝑟𝑜𝑐í𝑜; 𝑇𝑠𝑢𝑝 = 80,8588°𝐶
𝑃𝑖𝑠𝑎𝑡
145,926888 58,707356 24,191234
𝐾𝑖
2,881084 1,159079 0,477616
𝑥𝑖
0,156192 0,215688 0,628120
𝑓𝑙𝑎𝑠ℎ; 𝑇𝑠𝑢𝑝 = 71,77041 °𝐶
𝑃𝑖𝑠𝑎𝑡
111,301314 43,128349 17,112968
𝐾𝑖
2,197459 0,851498 0,337867
𝑥𝑖
0,281473 0,270052 0,448475
𝑦𝑖
0,618527 0,229948 0,151525
Para el inciso d).

Si la alimentación tiene una fracción vaporizada de f=0,5; y en equilibrio ideal el 90% del hexano se encuentra en la fase
de vapor, entonces también se espera que el 90% del octano de la alimentación esté en la fase de vapor
 6. El sistema que se representa seguidamente se utiliza para enfriar el efluente de un reactor y separar los gases
ligeros de los hidrocarburos más pesados. Calcúlese la composición y la velocidad de flujo del vapor que sale de
la cámara de flash. ¿Influye sobre el resultado la velocidad de flujo del líquido agotado?
Valores K para los componentes a 500 psia y 100°F:

K
H2 80
CH4 10
Benceno 0,010
Tolueno 0,004

Para el inciso a).

Primero precedemos a calcular la fracción molar de cada compuesto en la alimentación, primero calculamos el flujo
total:
𝐹𝑇 = ∑ 𝐹𝑖
𝐹𝑇 = 𝐹𝐻2 + 𝐹𝐶𝐻4 + 𝐹𝐵 + 𝐹𝑇
𝐹𝑇 = 2000 𝑘𝑚𝑜𝑙 + 2000 𝑘𝑚𝑜𝑙 + 500 𝑘𝑚𝑜𝑙 + 100 𝑘𝑚𝑜𝑙
𝐹𝑇 = 4600 𝑘𝑚𝑜𝑙
Ya con el flujo total calculamos las fracciones molares:
𝐹𝑖
𝑧𝑖 =
𝐹𝑇
Para el hidrogeno:
𝐹𝐻2
𝑧𝐻2 =
𝐹𝑇
2000 𝑘𝑚𝑜𝑙
𝑧𝐻2 =
4600 𝑘𝑚𝑜𝑙
𝑧𝐻2 = 0,43478
Para el metano:
𝐹𝐶𝐻4
𝑧𝐶𝐻4 =
𝐹𝑇
2000 𝑘𝑚𝑜𝑙
𝑧𝐶𝐻4 =
4600 𝑘𝑚𝑜𝑙
𝑧𝐶𝐻4 = 0,43478
Para el benceno:
𝐹𝐵
𝑧𝐵 =
𝐹𝑇
500 𝑘𝑚𝑜𝑙
𝑧𝐵 =
4600 𝑘𝑚𝑜𝑙
 𝑧𝐵 = 0,10870
Para el ciclo hexano:
𝐹𝑇
𝑧𝑇 =
𝐹𝑇
100 𝑘𝑚𝑜𝑙
𝑧𝑇 =
4600 𝑘𝑚𝑜𝑙
𝑧𝑇 = 0,02174
Convirtiendo los datos al sistema internacional de unidades tenemos que P = 3447,38 Kpa y que T = 37,78°C a la entrada
del flash.
Con valores de constantes para la ecuación de Antoine en T [=] °C, P [=] KPa:
A B C
H2 4,24710 2556,400 -67,0
Metano 5,82051 405,420 267,777
Benceno 6,03055 1211,033 220,790
Tolueno 6,07954 1344,800 219,482
Con la ecuación de Antoine:
𝐵
log 𝑃 = 𝐴 −
𝑇+𝐶
𝐵
𝑃 = 10𝐴−𝑇+𝐶
Para el hidrogeno:
2556,400
4,24710−
𝑃 = 10 37,78°𝐶−67,0

𝑃 = 5,43𝑥1091 𝐾𝑃𝑎
Para el metano:
405,420
5,82051−
𝑃 = 10 37,78°𝐶−267,777

𝑃 = 31166,428831 𝐾𝑃𝑎
Para el benceno:
1211,033
6,03055−
𝑃 = 10 37,78°𝐶−220,790

𝑃 = 22,231619 𝐾𝑃𝑎
Para el tolueno:
1344,800
6,07954−
𝑃 = 10 37,78°𝐶−219,482

𝑃 = 7,115153 𝐾𝑃𝑎
Calculamos las constantes de equilibrio:
𝑃𝑖𝑠𝑎𝑡
𝐾𝑖 =
𝑃𝑡𝑜𝑡𝑎𝑙

Para el metano:

31166,438831 𝐾𝑃𝑎
𝐾𝐶𝐻4 =
3447,38 𝐾𝑃𝑎
𝐾𝐶𝐻4 = 9,040616

Para el benceno:

22,231619 𝐾𝑃𝑎
𝐾𝐵 =
3447,38 𝐾𝑃𝑎
 𝐾𝑏 = 0,006449

Para el tolueno:

7,115153 𝐾𝑃𝑎
𝐾𝑇 =
3447,38 𝐾𝑃𝑎
𝐾𝑇 = 0,002064

Para el hidrógeno, se ha tomado un valor de Ki=50, basado en datos obtenidos de la bibliografía, ya que los resultados
indican que este elemento tiende a permanecer en fase gaseosa.
Ahora procedemos a calcular la fracción vaporizada con la ecuación de Rachford-Rice:
𝑧𝐹 (1 − 𝐾𝑖 )
∑ =0
1 + 𝑓(𝐾𝑖 − 1)
𝑧𝐹𝐻2 (1 − 𝐾𝐻2 ) 𝑧𝐹𝐶𝐻4 (1 − 𝐾𝐶𝐻4 ) 𝑧𝐹𝐵 (1 − 𝐾𝐵 ) 𝑧𝐹𝑇 (1 − 𝐾𝑇 )
+ + + =0
1 + 𝑓(𝐾𝐻2 − 1) 1 + 𝑓(𝐾𝐶𝐻4 − 1) 1 + 𝑓(𝐾𝐵 − 1) 1 + 𝑓(𝐾𝑇 − 1)
0,43478 ∗ (1 − 50) 0,43478 ∗ (1 − 9,040616) 0,10870 ∗ (1 − 0,006449) 0,02174 ∗ (1 − 0,002064)
+ + + =0
1 + 𝑓(50 − 1) 1 + 𝑓(9,040616 − 1) 1 + 𝑓(0,006449 − 1) 1 + 𝑓(0,002064 − 1)

𝑓 = −0,0727
𝑓 = 0,8655
𝑓 = 1,0028
Al resolver la ecuación se obtuvieron tres raíces, de las cuales tomaremos la segunda para hallar las fracciones molares en
cada una de las fases.

Para la fase líquida:

𝑧𝐹
𝑥𝑖 =
1 + 𝑓(𝐾𝑖 − 1)
Para el hidrogeno:
0,43478
𝑥𝐻2 =
1 + 0,8665 ∗ (50 − 1)
𝑥𝐻2 = 0,010016
Para el metano:
0,43478
𝑥𝐶𝐻4 =
1 + 0,8665 ∗ (9,040616 − 1)
𝑥𝐶𝐻4 = 0,054527
Para el benceno:
0,10870
𝑥𝐵 =
1 + 0,8665 ∗ (0,006449 − 1)
𝑥𝐵 = 0,775946
Para el tolueno:
0,02174
𝑥𝑇 =
1 + 0,8665 ∗ (0,002064 − 1)
𝑥𝑇 = 0,159511
Para calcular las fracciones molares en fase gaseosa usaremos:
𝑦𝑖
𝐾𝑖 =
𝑥𝑖
𝑦𝑖 = 𝐾𝑖 ∗ 𝑥𝑖
Para el hidrogeno:
 𝑦𝐻2 = 50 ∗ 0,010016
𝑦𝐻2 = 0,500792
Para el metano:
𝑦𝐶𝐻4 = 9,040616 ∗ 0,054527
𝑦𝐶𝐻4 = 0,493859
Para el benceno:
𝑦𝐵 = 0,006449 ∗ 0,775946
𝑦𝐵 = 0,005004
Para el tolueno:
𝑦𝑇 = 0,002064 ∗ 0,159511
𝑦𝑇 = 0,000345
Para conocer que cantidad de vapor y de líquido salen, se resuelve el siguiente sistema de ecuaciones:
𝐹 =𝐿+𝑉
{
𝐹𝑧𝐻2 = 𝐿𝑥𝐻2 + 𝑉𝑦𝐻2

4600 = 𝐿 + 𝑉
{
4600 ∗ 0,43478 = 𝐿 ∗ 0,010016 + 𝑉 ∗ 0,500792

4600 = 𝐿 + 𝑉
{
1999,988 = 𝐿 ∗ 0,010016 + 𝑉 ∗ 0,500792

𝐿 = 618,7246 𝑘𝑚𝑜𝑙/ℎ
{
𝑉 = 3981,275 𝑘𝑚𝑜𝑙/ℎ
La velocidad del flujo de vapor es de 3981,275 kmol/h, y su composición es 0,500792 para el hidrogeno, 0,493859 para el
metano, 0,005004 para el benceno y 0,000345 para el tolueno.

Par el inciso b)

No influirá pues esta recirculación ingresa al flujo de alimentación mucho antes de que la mezcla entre al flash, además
de que los datos de alimentación son dados en el momento en el que se unen ambas corrientes,
 7. A continuación, se muestra el sistema de cabeza de una columna de destilación. Se indica la composición total del
destilado, retirándose 10 moles % del mismo como vapor. Determínese la presión en el tanque de reflujo si la
temperatura es 100°F. Utilícense los valores K que se dan seguidamente para cualquier otra presión, suponiendo
que K es inversamente proporcional a la presión.

Composición K (100°F, 200psia)

C2 2,7
C3 0,95
C4 0,34

Con el dato de que en el destilado total se retira el 10% como vapor, asumimos que f = 0,1.

Conociendo valores de Ki a presión y temperatura conocidas podemos escalarlo a la nueva presión en el tanque de reflujo.

Si la presión en el tanque de reflujo es P, los nuevos valores de K estarán relacionados con la nueva presión como:

𝐾𝑖 (100°𝐹, 200 𝑝𝑠𝑖𝑎) ∗ 200

𝐾𝑖 (100°𝐹, 𝑃) =
𝑃
Usamos la ecuación de equilibrio para la fase vapor:

𝑦𝑖
𝐾𝑖 =
𝑥𝑖
𝑦𝑖
𝑥𝑖 =
𝐾𝑖
Y como la suma de las fracciones molares en la fase líquida debe ser 1:

∑𝑥𝑖 = 1
𝑦𝐶2 𝑦𝐶3 𝑦𝐶4
+ + =1
𝐾𝐶2 𝐾𝐶3 𝐾𝐶4
Conociendo valores de Ki a presión y temperatura conocidas podemos escalarlo a la nueva presión en el tanque de reflujo.
𝑦𝐶2 𝑦𝐶3 𝑦𝐶4
+ + =1
𝐾𝐶2 ∗ 200 𝐾𝐶3 ∗ 200 𝐾𝐶4 ∗ 200
𝑃 𝑃 𝑃
0,1 ∗ 𝑃 0,2 ∗ 𝑃 0,7 ∗ 𝑃
+ + =1
2,7 ∗ 200 0,95 ∗ 200 0,34 ∗ 200
0,1 0,2 0,7
( + + )∗𝑃 = 1
540 190 68
1
𝑃=
0,1 0,2 0,7
( + + )
540 190 68
𝑃 = 86,7157 𝑝𝑠𝑖𝑎
En este problema, si no se menciona ninguna caída de presión, suponemos que la presión del tanque de reflujo es la
misma que la del flash.