DESTILACIÓN FLASH

Andrés Benavides Prada, Ing. Qco.

Profesor Cátedra

Operaciones Unitarias II
Escuela de Ingeniería Química
Universidad Industrial de Santander
Bucaramanga, 2021 I
Destilación
Operación líquido/vapor que se fundamenta en la diferencia de puntos
de ebullición entre los compuestos de la mezcla a separar.

Conlleva un gran gasto de energía térmica. No se recomienda cuando la

volatilidad relativa de los compuestos es menor a 1,5. Para soluciones
ideales:
𝛼 = 𝑃𝑖𝑠𝑎𝑡 /𝑃𝑗𝑠𝑎𝑡

➢ Destilación flash.

➢ Destilación continua con reflujo (binaria).

➢ Destilación continua con reflujo (multicomponente).

Destilación flash
Parte de una mezcla líquida se vaporiza debido a una rápida disminución de
presión. La fase vapor es rica en compuesto volátil mientras que la fase
líquida es rica en compuesto pesado.

Flash isotérmico (𝑇𝐹 = 𝑇) Flash adiabático (𝑇𝐹 ≠ 𝑇)

𝑽, 𝒚𝒊 𝑽, 𝒚𝒊

𝑭, 𝒛𝒊 𝑭, 𝒛𝒊
𝑷, 𝑻 𝑷, 𝑻
𝑷𝑭 , 𝑻𝑭 𝑷𝑭 , 𝑻 𝑭
𝑳, 𝒙𝒊 𝑳, 𝒙𝒊

Para que la separación se lleve a cabo 𝑃𝑟𝑜𝑐í𝑜 < 𝑃 < 𝑃𝑏𝑢𝑟𝑏𝑢𝑗𝑎 .

Las fases vapor (𝑉, 𝑦𝑖 ) y líquido (𝐿, 𝑥𝑖 ) están en equilibrio termodinámico.

Equilibrio termodinámico
Si el sistema es ideal, es válida la ley de Raoult.
𝑦𝑖 𝑃 = 𝑃𝑖𝑠𝑎𝑡 𝑥𝑖 → 𝑦𝑖 = 𝐾𝑖 𝑥𝑖 𝐾𝑖 = 𝑃𝑖𝑠𝑎𝑡 /𝑃

Cuando la fase líquida es real pero la fase vapor se comporta idealmente, se

emplea la ley de Raoult modificada. Donde, 𝛾𝑖 es función de 𝑇 y 𝑥𝑖 .

𝑦𝑖 𝑃 = 𝛾𝑖 𝑃𝑖𝑠𝑎𝑡 𝑥𝑖

Para mezclas multicomponentes de hidrocarburos ligeros, los valores 𝐾𝑖 se

determinan por las gráficas de DePriester, las cuales fueron ajustadas a la
siguiente ecuación:
𝑎 𝑇1 𝑎 𝑇2 𝑎𝑃2 𝑎𝑃3
ln 𝐾 = 2 + + 𝑎 𝑇6 + 𝑎𝑃1 ln 𝑃 + 2 +
𝑇 𝑇 𝑃 𝑃
Donde, 𝑇 debe ser evaluada en R y 𝑃 en psia. Ecuación válida desde -70 hasta
200 °C, y desde 101 hasta 6000 kPa. La tabla que se muestra en la siguiente
diapositiva resume las constantes.
Constantes para el ajuste de valores de 𝐾. Fuente: Wankat, 2008.
 Flash isotérmico (𝑇𝐹 = 𝑇) Flash adiabático (𝑇𝐹 ≠ 𝑇)
𝑽, 𝒚𝒊 𝑽, 𝒚𝒊

𝑭, 𝒛𝒊 𝑭, 𝒛𝒊
𝑷, 𝑻 𝑷, 𝑻
𝑷𝑭 , 𝑻𝑭 𝑷𝑭 , 𝑻𝑭
𝑳, 𝒙𝒊 𝑳, 𝒙𝒊

Balance de materia

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

Agrupando:
𝑉 𝐹−𝑉
𝑧𝑖 = 𝑦𝑖 + 𝑥𝑖 → 𝑧𝑖 = ϕ𝑦𝑖 + (1 − ϕ)𝑥𝑖
𝐹 𝐹
Donde, ϕ es la fracción de mezcla vaporizada 𝑉/𝐹.
Considerando el equilibrio
líquido/vapor (𝑦𝑖 = 𝐾𝑖 𝑥𝑖 ):

𝑧𝑖
𝑧𝑖 = ϕ𝐾𝑖 𝑥𝑖 + (1 − ϕ)𝑥𝑖 𝑥𝑖 =
ϕ𝐾𝑖 + (1 − ϕ)

Por el principio de la conservación

de la masa:
𝒛𝒊
෍ 𝑥𝑖 = 1 ෍ =𝟏
𝝓𝑲𝒊 + (𝟏 − 𝝓)

La presión del flash es una variable de diseño, así como la alimentación al

proceso. Hay dos tipos de problemas:

❖ Se fija 𝜙 y se calcula 𝑇.

❖ Se fija 𝑇 y se calcula 𝜙.
➢ Se fija 𝜙.

➢ Las 𝐾𝑖 son desconocidas.

➢ Se calcula 𝑇 con la ecuación de balance.

➢ Con 𝑇 se determinan 𝐾𝑖 y se encuentran 𝑉, 𝑦𝑖 , 𝐿 y 𝑥𝑖 .

➢ Se fija 𝑇.

➢ Las 𝐾𝑖 son conocidas.

➢ Se calcula 𝜙 con la ecuación de balance.

➢ Se encuentran 𝑉, 𝑦𝑖 , 𝐿 y 𝑥𝑖 .
Balance de energía
𝑽, 𝒚𝒊

𝑭, 𝒛𝒊
𝐻𝐹 = ෍ 𝑧𝑖 𝐶𝑃𝑖 𝑇𝐹 − 𝑇𝑟𝑒𝑓
𝑷, 𝑻
𝑷𝑭 , 𝑻𝑭
𝑳, 𝒙𝒊
𝐻𝑣 = ෍ 𝑦𝑖 λ𝑖 + 𝐶𝑃𝑦𝑖 𝑇 − 𝑇𝑟𝑒𝑓
𝑭𝑯𝑭 + 𝑸 = 𝑽𝑯𝒗 + 𝑳𝑯𝑳

𝐻𝐿 = ෍ 𝑥𝑖 𝐶𝑃𝑖 𝑇 − 𝑇𝑟𝑒𝑓
𝑽, 𝒚𝒊

𝑭, 𝒛𝒊 Donde, 𝐶𝑃 y 𝐶𝑃𝑦 son las capacidades

𝑷, 𝑻 caloríficas promedio en fase líquida y
𝑷𝑭 , 𝑻𝑭
𝑳, 𝒙𝒊
vapor, respectivamente. λ es el calor
latente de vaporización a la temperatura
de referencia.
𝑭𝑯𝑭 = 𝑽𝑯𝒗 + 𝑳𝑯𝑳
Para mezclas binarias

Entalpía de la alimentación:

𝐻𝐹 = 𝑧𝐴 𝐶𝑃𝐴 𝑇𝐹 − 𝑇𝑟𝑒𝑓 + (1 − 𝑧𝐴 )𝐶𝑃𝐵 𝑇𝐹 − 𝑇𝑟𝑒𝑓

Entalpía del vapor:

𝐻𝑉 = 𝑦𝐴 λ𝐴 + 𝐶𝑃𝑦𝐴 𝑇 − 𝑇𝑟𝑒𝑓 + (1 − 𝑦𝐴 ) λ𝐵 + 𝐶𝑃𝑦𝐵 𝑇 − 𝑇𝑟𝑒𝑓

Entalpía del líquido:

𝐻𝐿 = 𝑥𝐴 𝐶𝑃𝐴 𝑇 − 𝑇𝑟𝑒𝑓 + (1 − 𝑥𝐴 )𝐶𝑃𝐵 𝑇 − 𝑇𝑟𝑒𝑓

Si la temperatura del flash (𝑇) se fija como referencia, la entalpía del

líquido es cero (𝐻𝐿 = 0) y la entalpía del vapor se simplifica.

𝐻𝑉 = 𝑦𝐴 λ𝐴 + (1 − 𝑦𝐴 )λ𝐵
Flash adiabático (mezclas binarias)
Ejemplo 1. Una mezcla con
𝑪𝑷 𝑻𝑭 − 𝑻 = 𝝓 𝒚𝑨 𝝀𝑨 + 𝟏 − 𝒚𝑨 𝝀𝑩
equimolar benceno/tolueno se
somete a una destilación
Flash isotérmico (mezclas binarias) instantánea a 1 atm. Calcule las
composiciones del vapor y del
𝑸 = 𝑽 𝒚𝑨 𝝀 𝑨 + 𝟏 − 𝒚 𝑨 𝝀 𝑩
líquido que abandonan el proceso
si la temperatura se fija en 95 °C.
Benceno Tolueno ¿Cuál es la temperatura de
𝜆 (𝑐𝑎𝑙/𝑚𝑜𝑙) 7360 7960 alimentación requerida para una
𝐶𝑃 (𝑐𝑎𝑙/𝑚𝑜𝑙°𝐶) 33 40 separación adiabática?

2788,5 3096,52
ln 𝑃𝐴𝑠𝑎𝑡 = 15,9 − ln 𝑃𝐵𝑠𝑎𝑡 = 16,04 −
𝑇 − 52,36 𝑇 − 53,67

Donde, 𝑃 está en mmHg y 𝑇 en K.

Se debe comprobar que la presión de operación sea adecuada.

2788,5
ln 𝑃𝐴𝑠𝑎𝑡 = 15,9 − → 𝑃𝐴𝑠𝑎𝑡 = 1171 𝑚𝑚𝐻𝑔
368 − 52,36

3096,52
ln 𝑃𝐵𝑠𝑎𝑡 = 16,04 − → 𝑃𝐵𝑠𝑎𝑡 = 487 𝑚𝑚𝐻𝑔
368 − 53,67

Punto de rocío Punto de burbuja

1
𝑃𝑟𝑜𝑐í𝑜 = 𝑃𝑏𝑢𝑟𝑏𝑢𝑗𝑎 = ෍ 𝑧𝑖 𝑃𝑖𝑠𝑎𝑡
σ 𝑧𝑖 /𝑃𝑖𝑠𝑎𝑡

1 𝑃𝑏𝑢𝑟𝑏 = 0,5 1171 + 0,5 487

𝑃𝑟𝑜𝑐í𝑜 =
0,5 0,5
+
1171 487 𝑃𝑏𝑢𝑟𝑏 = 829 𝑚𝑚𝐻𝑔

𝑃𝑟𝑜𝑐í𝑜 = 688 𝑚𝑚𝐻𝑔

Constantes de equilibrio

1171 487
𝐾𝐴 = = 1,54 𝐾𝐵 = = 0,64
760 760

Balance de materia

0,5 0,5
+ =1
𝜙(1,54) + (1 − 𝜙) 𝜙(0,64) + (1 − 𝜙)

𝝓 = 𝟎, 𝟒𝟔𝟑
Composición en las corrientes de salida

0,5
𝑥𝐴 =
(0,463)(1,54) + (1 − 0,463) 𝒙𝑨 = 𝟎, 𝟒

𝑥𝐵 = 1 − 𝑥𝐴 𝒙𝑩 = 𝟎, 𝟔
 𝑦𝐴 = (1,54)(0,4) 𝒚𝑨 = 𝟎, 𝟔𝟐

𝑦𝐵 = 1 − 𝑦𝐴 𝒚𝑩 = 𝟎, 𝟑𝟖

Balance de energía
𝐶𝑃 = 0,5 33 + 0,5 40 = 36,5 𝑐𝑎𝑙/𝑚𝑜𝑙°𝐶

(36,5) 𝑇𝐹 − 95 = (0,463) (0,62) 7360 + 0,38 (7960)

𝑻𝑭 = 𝟏𝟗𝟏, 𝟑 °𝑪

Diagrama de equilibrio
benceno/tolueno. Fuente:
McCabe at al., 2007.
Ejemplo 2. Una corriente de 100 kmol/h tiene la siguiente composición de
hidrocarburos ligeros:

10% Propano
20% n-Butano
30% n-Pentano
40% n-Hexano

La mezcla se somete a una destilación flash a 100 psia, donde se evapora un

12,2 % de la alimentación. Determine la temperatura a la cual opera el
equipo, las composiciones de las corrientes de vapor y líquido, y el calor
suministrado si la operación es isotérmica. Para ajustar el calor latente de
vaporización a la temperatura de referencia, emplee la siguiente ecuación
(el subíndice 𝑏 indica punto de ebullición):

𝜆𝑖 = 𝐶𝑃𝑖 𝑇𝑏𝑖 − 𝑇𝑟𝑒𝑓 + 𝜆𝑏𝑖 − 𝐶𝑃𝑦𝑖 𝑇𝑏𝑖 − 𝑇𝑟𝑒𝑓

Ecuación de balance
𝑧𝑖
෍ =1
𝜙𝐾𝑖 + (1 − 𝜙)

0,1 0,2 0,3 0,4

+ + + =1
0,122𝐾𝑖 + 0,878 0,122𝐾𝑖 + 0,878 0,122𝐾𝑖 + 0,878 0,122𝐾𝑖 + 0,878

𝑎 𝑇1 𝑎 𝑇2 𝑎𝑃2 𝑎𝑃3
ln 𝐾𝑖 = 2 + + 𝑎 𝑇6 + 𝑎𝑃1 ln 100 + 2
+
𝑇 𝑇 (100) (100)
 Algoritmo de solución

1. Suponer 𝑇. 3. Comprobar: σ 𝑥𝑖 = 1.

2. Calcular 𝐾𝑖 . 4. Volver a suponer 𝑇 si la sumatoria es

diferente de 1.

Para una temperatura de 660 R

(94 °C):
𝐾𝐶3 = 4,26 𝐾𝐶4 = 1,77

𝐾𝐶5 = 0,76 𝐾𝐶6 = 0,36

0,1 0,2 0,3 0,4

𝑓= + + +
(0,122)(4,26) + 0,878 (0,122)(1,77) + 0,878 (0,122)(0,76) + 0,878 (0,122)(0,36) + 0,878

𝒇≈𝟏 V= 𝟏𝟐, 𝟐 𝒌𝒎𝒐𝒍/𝒉 L= 𝟖𝟕, 𝟖 𝒌𝒎𝒐𝒍/𝒉

Composiciones de salida

0,1
𝑥𝐶3 = 𝑦𝐶3 = 4,26 0,072
(0,122)(4,26) + 0,878

𝒙𝑪𝟑 = 𝟎, 𝟎𝟕𝟐 𝒚𝑪𝟑 = 𝟎, 𝟑𝟎𝟒

0,2
𝑥𝐶4 = 𝑦𝐶4 = 1,77 0,183
(0,122)(1,77) + 0,878

𝒙𝑪𝟒 = 𝟎, 𝟏𝟖𝟑 𝒚𝑪𝟒 = 𝟎, 𝟑𝟐𝟑

0,3
𝑥𝐶5 = 𝑦𝐶5 = 0,76 0,309
(0,122)(0,76) + 0,878

𝒙𝑪𝟓 = 𝟎, 𝟑𝟎𝟗 𝒚𝑪𝟓 = 𝟎, 𝟐𝟑𝟓

𝒙𝑪𝟔 = 𝟎, 𝟒𝟑𝟔 𝒚𝑪𝟔 = 𝟎, 𝟏𝟑𝟖

Propiedades térmicas

𝑪𝑷 (𝒌𝑱/𝒌𝒎𝒐𝒍°𝑪) 𝝀 (𝒌𝑱/𝒌𝒎𝒐𝒍) 𝑻𝒃 (°𝑪)

𝑪𝟑 99 15700 -42
𝑪𝟒 134 22440 0
𝑪𝟓 168 25790 36
𝑪𝟔 193 28850 69

෣𝐶 = 𝑅 1,213 + 28,785𝑥10−3 𝑇 − 8,824𝑥10−6 𝑇 2

𝐶𝑝𝑦 3

𝑇ത = 299 𝐾 → ෣𝐶 = 75,12 𝑘𝐽/𝑘𝑚𝑜𝑙°𝐶

𝐶𝑝𝑦 3

෣𝐶 = 𝑅 1,935 + 36,915𝑥10−3 𝑇 − 11,402𝑥10−6 𝑇 2

𝐶𝑝𝑦 4

𝑇ത = 320 𝐾 → ෣𝐶 = 104,6 𝑘𝐽/𝑘𝑚𝑜𝑙°𝐶

𝐶𝑝𝑦 4
 ෣𝐶 = 𝑅 2,464 + 45,351𝑥10−3 𝑇 − 14,111𝑥10−6 𝑇 2
𝐶𝑝𝑦 5

𝑇ത = 338 𝐾 → ෣𝐶 = 134,5 𝑘𝐽/𝑘𝑚𝑜𝑙°𝐶

𝐶𝑝𝑦 5

෣𝐶 = 𝑅 3,025 + 53,722𝑥10−3 𝑇 − 16,791𝑥10−6 𝑇 2

𝐶𝑝𝑦 6

𝑇ത = 354,5 𝐾 → ෣𝐶 = 165,9 𝑘𝐽/𝑘𝑚𝑜𝑙°𝐶

𝐶𝑝𝑦 6

Calor latente (a 𝑇𝑟𝑒𝑓 ):

𝜆𝐶3 = 99 −42 − 94 + 15700 − 75,12 −42 − 94

𝜆𝐶3 = 12452 𝑘𝐽/𝑘𝑚𝑜𝑙

𝜆𝐶4 = 134 0 − 94 + 22440 − 104,6 0 − 94

𝜆𝐶4 = 19679 𝑘𝐽/𝑘𝑚𝑜𝑙

 𝜆𝐶5 = 168 36 − 94 + 25790 − 134,5 36 − 94

𝜆𝐶5 = 23851 𝑘𝐽/𝑘𝑚𝑜𝑙

𝜆𝐶6 = 193 69 − 94 + 28850 − 165,9 69 − 94

𝜆𝐶6 = 28174 𝑘𝐽/𝑘𝑚𝑜𝑙

Balance de energía

𝑄 = (12,2) (0,304)(12452) + (0,323)(19679) + (0,235)(23851) + (0,138)(28174)

𝑸 = 𝟐𝟑𝟗𝟓𝟒𝟒 𝒌𝑱/𝒉