Eberth Julian Rodriguez Redondo

Cod: 2143253
Enunciado
 Inciso A: Flash Adiabático

Código: 2143253
𝑋 5
𝐹 = 1000 ∗ = 1000 = 1666,666𝑘𝑚𝑜𝑙/ℎ
𝑌 3
𝑇𝐹 = 115°𝐶
𝑃 = 200 𝐾𝑝𝑎

Compuesto Zi Teb [°C] Teb [K] λ [kJ/kmol] Cp [kJ/kmol*K]

Propano 0,3 -42 231,15 15700 99
n-Butano 0,1 0 273,15 22440 134
n-Pentano 0,15 36 309,15 25790 168
n-Hexano 0,45 69 342,15 28850 193
P operación [kpa] 200 P operación [psi] 29,015544
1er Paso
Se debe suponer una temperatura para calcular las Psat de cada compuesto.
También tener en cuenta que a presión calculada debe estar entre
Procio<Poperación<Pburbuja.
𝑇𝑠𝑢𝑝𝑢𝑒𝑠𝑡𝑎 ≠ 𝑇𝐹
5
𝑇𝑠𝑢𝑝𝑢𝑒𝑠𝑡𝑎 = 580,43°𝑅 = 580,4 − 491,67 ∗ = 49,3°𝐶
9
𝑇𝑠𝑢𝑝𝑢𝑒𝑠𝑡𝑎 = 49,3°𝐶 = 322,45𝐾°

Con esta temperatura supuesta, usamos la ecuación

Antoine. Y calculamos la 𝑃𝑟𝑜𝑐í𝑜 𝑦 𝑃 𝑏𝑢𝑟𝑏𝑢𝑗𝑎.

Ctes. Antoine A B C P sat [bar] P sat [kpa]

Propano 3,92828 803,997 247,04 16,4130883 1641,308829
n-Butano 3,93266 935,773 238,789 4,83561795 483,5617949
n-Pentano 3,97786 1064,84 232,014 1,55820134 155,8201341
n-Hexano 4,00139 1170,875 224,317 0,52740228 52,74022765

𝐵
𝐿𝑜𝑔10 𝑃𝑠𝑎𝑡 𝑏𝑎𝑟 =𝐴−
𝐶 + 𝑇(𝐾)
2do Paso
Con los datos de Psat debemos calcular la Procio y Pburbuja para comprobar si
Pop=200Kpa se encuentra dentro de estos valores

𝑃burbuja= (0.3 ∗ 16.91) + (0.1 ∗ 4.85) + (0.15 ∗ 1.55) + (0.45 ∗ 0.527) = 6,028𝑏𝑎𝑟 = 602,84Kpa

𝑃𝑟𝑜𝑐𝑖𝑜 < 𝑃𝑜𝑝 < 𝑃𝑏𝑢𝑟𝑏𝑢𝑗𝑎

𝑃𝑟𝑜𝑐í𝑜= 0,3 1
0,1 0,15 0,45 = 1,012𝑏𝑎𝑟 = 101,202𝐾𝑝𝑎
+ + +
16,91 4,85 1,55 0,527 101.202 𝐾𝑝𝑎 < 200𝐾𝑝𝑎 < 602.84𝐾𝑝𝑎
𝑆𝑖 𝑐𝑢𝑚𝑝𝑙𝑒 𝑒𝑠𝑡𝑎 𝑐𝑜𝑛𝑑𝑖𝑐𝑖ó𝑛
3er paso
Se deben calcular los Ki a partir de la ecuación de DePriester.
Debemos tener las temperaturas en R° y P en psia.

𝑇𝑠𝑢𝑝𝑢𝑒𝑠𝑡𝑎 = 580°𝑅
1𝑝𝑠𝑖𝑎
𝑃 = 200 𝐾𝑝𝑎 ∗ = 29𝑝𝑠𝑖𝑎
6,89476𝐾𝑝𝑎
 3er

Con los valores Ki debemos

calcular la fracción
vaporizada de la mezcla.(Ф)

𝑎𝑇1 𝑎𝑇2
+𝑇
𝑎 𝑎
+𝑎 𝑇6 +𝑎 𝑃1 ln 𝑃 + 𝑃 2𝑃2+ 𝑃𝑃3 𝑑𝑜𝑛𝑑𝑒 𝑓 ≈ 1
K= 𝑒 𝑇2

aT1 aT2 aT6 aP1 aP2 aP3 Ki

Propano -970688,563 0 7,15059 -0,76984 0 6,90224 6,78243049
n-Butano -1280557 0 7,94986 -0,96455 0 0 2,4599579
n-Pentano -1524891 0 7,33129 -0,89143 0 0 0,8208097
n-Hexano -1778901 0 6,96783 -0,84634 0 0 0,31252048
 0.3 0.1 0.15 0.45
𝑓≈ + + +
Ø∗6.785+ 1−Ø Ø∗2.461+ 1−Ø Ø∗0.821+ 1−Ø Ø∗0.312+ 1−Ø

Se resuelve por medio solver en Excel o en la calculadora

Texas, donde el valor da Ø=0.5239 con 𝑓 ≈ 1 = 0,9999

4to paso 5to paso

Hallar las fracciones de la salida 𝑥𝑖 𝑦 𝑦𝑖.
Hallar los Flujos a partir del balance de masa
y de la fracción vaporizada Ø=V/F .

Por ley de Raoult se puede hallar la

fracción de vapor de cada compuesto.
V=F∗Ø = 873,23 𝑘𝑚𝑜𝑙/ℎ
Compuesto Xi Yi
𝐿 = 793,428298 𝑘𝑚𝑜𝑙/ℎ
Propano 0,07437864 0,50472418
n-butano 0,05662002 0,1393621
n-pentano 0,16550537 0,13593075
n-hexano 0,70348157 0,21999604
0,9999856 1,00001307
6to paso
Calcular ₮, 𝐶𝑝𝑦𝑖, 𝜆𝑖, 𝐻𝑣
Compuesto Zi Teb [°C] Teb [K] λ [kJ/kmol] Cp [kJ/kmol*K]
Propano 0,3 -42 231,15 15700 99
n-Butano 0,1 0 273,15 22440 134
n-Pentano 0,15 36 309,15 25790 168
n-Hexano 0,45 69 342,15 28850 193
P operación [kpa] 200 P operación [psi] 29,015544

1. Tprom

𝑇𝑠𝑢𝑝𝑢𝑒𝑠𝑡𝑎 = 𝑇 = 49.31°𝐶 = 322.46 𝐾

2. Cpy

Cap. Caloríficas A B C Tprom [K] Cpy

Propano 1,213 0,028785 -8,824E-06 276,8 70,70731632
n-Butano 1,935 0,036915 -1,1402E-05 297,8 99,07878981
n-Pentano 2,464 0,045351 -1,4111E-05 315,8 127,8573549
n-Hexano 3,025 0,053722 -1,6791E-05 332,3 158,1547539
 3. 𝜆𝑖
𝑇𝑏𝑖 = 𝑇𝑒𝑏

99𝑘𝐽 𝑘𝐽 𝑘𝐽
𝜆𝑖 = ∗ 231.15 − 322.46 𝐾 + 15700 − 70.708 231.15 − 322.46 𝐾
𝑘𝑚𝑜𝑙 ∗ 𝐾 𝑘𝑚𝑜𝑙 𝑘𝑚𝑜𝑙 ∗ 𝐾
𝑘𝐽
𝜆𝑖 = 13116.63
𝑘𝑚𝑜𝑙

4. Hv
𝑇𝑟𝑒𝑓 = 𝑇

𝐻𝑣 = 0.505 ∗ 13116.63 + 0.139 ∗ 20718.02 + 0.136 ∗ 25255.62 + 0.220 ∗ 29535.97

Hv=19438.26 kJ/kmol

Cap. Caloríficas Tprom [K] Cpy Lambda ajust. Hv

Propano 276,8 70,70731632 13116,87798 6623,20562
n-Butano 297,8 99,07878981 20718,38434 2887,71883
n-Pentano 315,8 127,8573549 25256,10282 3431,76072
n-Hexano 332,3 158,1547539 29536,45135 6492,39803
total Hv 19435,0832
 7mo paso
Calcular el Cp promedio y hacer el balance de 𝐶𝑝 = 155.15 𝐾𝐽/𝑘𝑚𝑜𝑙 ∗ 𝐾
energia.
𝐻𝑣 = 19438.26 kJ/kmol
𝐶𝑝 = 𝑧𝐴𝐶𝑝𝐴 + (1 − 𝑧𝐴)𝐶𝑝𝐵 𝑇𝐹 = 115°𝐶
Ø=0.5239
𝐶𝑝 = 0,3 ∗ 99 + 0,1 ∗ 134 + 0,12 ∗ 168 + 0,45 ∗ 193
𝐾𝐽
𝐶𝑝 = 155,15
𝑘𝑚𝑜𝑙 ∗ 𝑘
Hacer el balance de energía, para determinar que el Q=0 sea
adiabático y encontrar la temperatura hallada aplicada para ver
si es adecuada o no.
Reemplazo y calculo por medio Excel o solver de
Texas y llega a este valor de la Temperatura hallada
𝐶𝑝 𝑇𝐹 − 𝑇 = ϕ 𝑦𝐴λ𝐴 + (1 − 𝑦𝐴)λ𝐵
Flash adiabático (Q=0) 𝑇𝑠𝑢𝑝𝑢𝑒𝑠𝑡𝑎 = 𝑇 = 49.31°𝐶 = 322.46𝐾
𝐸𝑟𝑟𝑜𝑟 = 𝑇ℎ𝑎𝑙𝑙𝑎𝑑𝑎 − 𝑇𝑠𝑢𝑝𝑢𝑒𝑠𝑡𝑎 = 3 ∗ 10−4
𝐶𝑝 𝑇𝐹 − 𝑇 = ϕ ∗ 𝐻𝑣 𝐶𝑝 𝑇𝐹 − 𝑇 − ϕ ∗ 𝐻𝑣 = 0 𝑇 = 49.3123°𝐶
𝐾𝐽 Esta es la temperatura aplica para esta
𝐶𝑝 = 155,15 ; 𝐻𝑣 = 19438,26 ; 𝑇𝐹 = 115°𝐶 ; ϕ = 0,5239 operación de destilación Flash Adiabática
𝑘𝑚𝑜𝑙 ∗ 𝑘
Inciso B

Ø=0.5239
𝐹 = 1666,66𝑘𝑚𝑜𝑙/ℎ 𝑉 = 873,23 𝑘𝑚𝑜𝑙/ℎ

𝐿 = 793,428298 𝑘𝑚𝑜𝑙/ℎ

𝑇𝐹 = 115°𝐶 = 388,15𝐾 = 698,67°𝑅

Compuesto Zi Teb [°C] Teb [K] λ [kJ/kmol] Cp [kJ/kmol*K]

Propano 0,3 -42 231,15 15700 99
n-Butano 0,1 0 273,15 22440 134
n-Pentano 0,15 36 309,15 25790 168
n-Hexano 0,45 69 342,15 28850 193
P operación [kpa] 200 P operación [psi] 29,015544
 1er paso
Lo primero que se debe hacer es calcular las Psat para cada compuesto con la temperatura de operación.

𝑇𝐹 = 388,15𝐾
𝐵
𝐿𝑜𝑔10 𝑃𝑠𝑎𝑡 𝑏𝑎𝑟 =𝐴−
𝐶 + 𝑇(𝐾)

Ctes. Antoine A B C P sat [bar] P sat [kpa]

Propano 3,92828 803,997 247,04 50,9963623 5099,63623
n-Butano 3,93266 935,773 238,789 19,3935317 1939,35317
n-Pentano 3,97786 1064,84 232,014 8,11485258 811,485258
n-Hexano 4,00139 1170,875 224,317 3,55392443 355,392443

Luego hallar
𝑃𝑟𝑜𝑐í𝑜 𝑦 𝑃 𝑏𝑢𝑟𝑏𝑢𝑗𝑎.

𝑃𝑟𝑜𝑐𝑖𝑜 < 𝑃𝑜𝑝 < 𝑃𝑏𝑢𝑟𝑏𝑢𝑗𝑎

640,43𝑘𝑃𝑎 < 𝑃𝑜𝑝 < 2005,475𝑘𝑃𝑎
P rocío [kpa] 640,4330652
P burb [kpa] 2005,475574
 2do paso
Asumo una presión de operación entre los valores de Pburbuja y Procio para asegurar separación.
Determinar los términos de Ki a través de la ecuación de Depriester.

𝑃𝑜𝑝𝑒 = 10𝑏𝑎𝑟 = 1000𝐾𝑝𝑎

𝑎𝑇1 𝑎𝑇2 𝑎𝑃2 𝑎𝑃3
𝐿𝑛 𝐾 = + + 𝑎 𝑇6 + 𝐴 𝑃1 ln 𝑃 + +
𝑇2 𝑇 𝑃2 𝑃
Se toman los mismos valores para la ecuación de Depriester
del inciso A.
𝑎𝑇1 𝑎𝑇2 𝑎 𝑎
+ 𝑇 +𝑎 𝑇6 +𝑎 𝑃1 ln 𝑃 + 𝑃 2𝑃2+ 𝑃𝑃3
K= 𝑒 𝑇2

aT1 aT2 aT6 aP1 aP2 aP3 Ki

Propano -970688,5625 0 7,15059 -0,76984 0 6,90224 4,05495895
n-Butano -1280557 0 7,94986 -0,96455 0 0 1,73613115
n-Pentano -1524891 0 7,33129 -0,89143 0 0 0,81425449
n-Hexano -1778901 0 6,96783 -0,84634 0 0 0,42059105
 aT1 aT2 aT6 aP1 aP2 aP3 Ki
Propano -970688,5625 0 7,15059 -0,76984 0 6,90224 4,05495895
n-Butano -1280557 0 7,94986 -0,96455 0 0 1,73613115
n-Pentano -1524891 0 7,33129 -0,89143 0 0 0,81425449
n-Hexano -1778901 0 6,96783 -0,84634 0 0 0,42059105

𝑍𝑖
𝑓=෍
Ø 𝐾𝑖 + (1 − Ø)
𝑍𝑖
𝑓=෍
𝑃𝑠𝑎𝑡
Ø
𝑃
+ (1 − Ø)
3er paso
Calculados los Ki procedemos
En donde f=1
a calcular la Presión de
Solver para hallar Pop
operación

𝑃𝑜𝑝 = 141.241𝑝𝑠𝑖𝑎 = 973.82kPa = 9.738 bar

4to paso
Hallar las fracciones de la salida 𝑥𝑖 𝑦 𝑦𝑖.
Compuesto Xi Yi
Propano 0,11532039 0,46761945
n-Butano 0,07215584 0,125272
n-Pentano 0,16618107 0,13531369
Por ley de Raoult se puede hallar la n-Hexano 0,64630423 0,27182977
fracción de vapor de cada compuesto. 0,99996154 1,00003491

5to paso
Calcular ₮, 𝐶𝑝𝑦𝑖, 𝜆𝑖,

1.Calcular la T prom
2. Cpy

Cap. Caloríficas
A B C Tprom [K] Cpy
Propano 1,213 0,028785 -8,824E-06 309,65 77,1556042
n-Butano 1,935 0,036915 -1,1402E-05 330,65 107,203798
n-Pentano 2,464 0,045351 -1,4111E-05 348,65 137,682649
n-Hexano 3,025 0,053722 -1,6791E-05 365,15 169,628616

3. 𝜆𝑖

Compuesto Composicion Teb(C) Teb(K) Cp KJ/kmol*K λ(kJ/kmol) Lambda ajust.

propano 0.3 -42 231.15 99 15700 12270,42986
n-butano 0.1 0 273.15 134 22440 19358,43675
n-pentano 0.15 36 309.15 168 25790 23394,92925
n-hexano 0.45 69 342.15 193 28850 27774,91633

λ𝑖 = 𝐶𝑝𝑖 𝑇𝑏𝑖 − 𝑇𝑟𝑒𝑓 + λ𝑏𝑖 − 𝐶𝑝𝑦𝑖(𝑇𝑏𝑖 − 𝑇𝑟𝑒𝑓)

λ 𝑝𝑟𝑜𝑝 = 99 231,15 − 388,15 + 15700 − 77,15 ∗ (231,15 − 388,15)

𝐾𝐽
λ 𝑝𝑟𝑜𝑝 = 12270,43
𝑘𝑚𝑜𝑙
 6to paso
Hacer el balance de energía, para determinar el
Qrequerido. 𝑉 = 873,23 𝑘𝑚𝑜𝑙/ℎ

𝑄 = 𝑉[𝑦𝐴λ𝐴 + (1 − 𝑦𝐴 ) λ 𝐵 ]
𝑄 = 873,23 ∗ (0,467 ∗ 12270,43 + 0,125 ∗ 19358,44 + 0,135 ∗ 23394,93 + 0,2719 ∗ 27774,92)
𝐾𝐽
𝑄=16494057,01 ( )
ℎ

Proceso Adiabatico
Compuesto xi yi ki Cpyi(KJ/LmolK) Top(°C) Pope (Kpa) Ф=V/F Q(KJ/h)
Propano 0,07444789 0,504937656 6,782430 70,7073163
n-butano 0,05665932 0,139379538 2,459958 99,0787898
49,3 200 0,5239 0
n-pentano 0,16554199 0,135878474 0,820810 127,8573549
nhexano 0,7033449 0,219809684 0,312520 158,1547539
Proceso Isotérmico
Compuesto xi yi ki Cpyi(KJ/LmolK) Top(°C) Pope (Kpa) Ф=V/F Q(KJ/h)
Propano 0,11532039 0,46761945 4,05495895 77,15560423
n-butano 0,07215584 0,125272 1,73613115 107,2037978
n-pentano 137,6826487
115 973,82483 0,5239 16494057,01
0,16618107 0,13531369 0,81425449
nhexano 0,64630423 0,27182977 0,42059105 169,6286159
 Conclusiones

• La presión de operación en el flash adiabática es mayor que la presión

de flash isotérmica.
• La operación del flash de forma adiabática requiere de una chaqueta
(aislante térmico) que evite el intercambio de calor con el medio. A
pesar de esto el flash adiabático sigue siendo mejor que el
isotérmico.
• En este caso se recomienda el sistema adiabático debido a que el
isotérmico requiere una alta presión, además de un alto gasto
energético para mantener la temperatura constante, lo que indica un
alto costo de operación.