Reporte FlashAdiabatico

Modelado de Flash adiabático empleando la
ecuacion de estado Soave-Redlich-Kwong

Reporte #2
M. A. Martínez-Moreno a, J. D. Ponce-Rocha a, J. A. Sánchez-Gómez a, O. Sotelo-Pineda a.

Universidad de Guanajuato
División de Ciencias Naturales y Exactas
Maestría en Ingeniería Química
Procesos de Separación (IIMA06031)
Departamento de Ingeniería Química
Guanajuato, Guanajuato, México
Resumen
Los ingenieros químicos dentro de las industrias petroleras, farmacéuticas, alimentos, polímeros entre
otras los destiladores Flash tienen una posición importante dentro de los procesos químicos, estos son
esenciales para la mejora de los procesos ya existentes y para diseño de otros nuevos. El flash es una
etapa de separación de dos fases que se encuentran en equilibrio. La composición de cada fase está
determinada por las condiciones de presión y temperatura. El flash adiabático es un equipo de
separación en donde la carga térmica es cero. En el presente trabajo se analiza la separación de una
mezcla de hidrocarburos de la industria petroquímica. El modelo utilizado para la solución es Soave-
Redlich-Kwong. Se implementó el programa Python para la solución del problema y la obtención de
los resultados.
Keywords—Flash Adiabático, Modelo termodinámico, Soave-Redlich-Kwong, Rachford-Rice,

Método de Newton;
Introducción
Desde que Van der Waals propuso la primera ecuación de estado (EOS, por sus siglas en inglés) hace
más de un siglo, se han propuesto diversas modificaciones en la literatura para la mejora en las
predicciones del equilibrio de fases, predicciones volumétricas y propiedades termodinámicas.
Aunque la ecuación de Van der Waals no es precisa para la mayoría de las aplicaciones, puede
considerarse una gran contribución en el área de la termodinámica para la descripción del
comportamiento PVT.
De manera habitual el uso de ecuaciones de estado cubicas predicen de manera adecuada el equilibrio
líquido-vapor de ciertas mezclas multicomponentes, tal es el caso de la ecuación cubica de estado
descrita por Soave-Redlich-Kwong.
Inicialmente Redlich-Kwong, habían desarrollado una ecuación de estado capaz de describir de
manera correcta sistemas de baja y alta densidad, si bien la ecuación de estado desarrolla no tenía una
base teórica sólida, demostró dar buenos resultados para muchos sistemas gaseosos. Posteriormente
Wilson realizo una contribución importante a los intentos por generalizar la ecuación de Redlich-
Kwong, de este modo, logro considerar las variaciones en el comportamiento de diferentes fluidos a
la misma temperatura y presión reducida introduciendo el factor acéntrico de Pietzer en el término de
atracción.
Posterior a la propuesta de Wilson, dicho modelo fue ignorado durante la década de 1960 y hasta
1972, momento en el cual Soave postula una versión mejorada al modelo de Wilson, redefiniendo la
función α (𝑇𝑅 , ω) y manteniendo la funcionalidad de volumen del modelo de Redlich-Kwong. De este
modo la ecuación de estado Soave-Redlich-Kwong fue defina por:
𝐑𝐓 𝒂(𝑻)
𝐏= − (𝟏)
𝐕 − 𝐛 𝑽(𝑽 + 𝒃)
𝟎. 𝟒𝟐𝟕𝟒𝟖𝐑𝟐 𝑻𝟐𝒄
𝐚(𝐓) = 𝛂(𝑻) (𝟐)
𝑷𝒄
𝐑𝑻𝒄
𝐛 = 𝟎. 𝟎𝟖𝟔𝟔𝟒 (𝟑)
𝑷𝒄
𝑻
√𝛂(𝑻) = 𝟏 + 𝛋 (𝟏 − √ ) (𝟒)
𝑻𝒄
𝛋 = 𝟎. 𝟒𝟖𝟎 + 𝟏. 𝟓𝟕𝟒𝛚 − 𝟎. 𝟏𝟕𝟔𝛚𝟐 (𝟓)
En la actualidad el uso de la ecuación de estado desarrollada por Soave-Redlich-Kwong es empleada

en diversos procesos industriales, enfocados principalmente en sistemas de baja densidad, tales como
procesos de separación de aire, procesamiento de gas, refinación de petroquímicos en equilibrio
líquido-vapor.
Flash adiabático
Un flash consiste en una sola etapa de destilación en la que la alimentación se vaporiza parcialmente
para producir un vapor más rico en el componente más volátil. Si el equipo esta adecuadamente
diseñado, el líquido y el vapor que salen del recipiente están en equilibrio. Cuando la presión de una
corriente líquida de composición, velocidad de flujo y temperatura (o entalpía) conocidas, se reduce
adiabáticamente a través de una válvula, tal como se muestra en la figura 1, es posible realizar un cálculo
de flash adiabático para determinar la temperatura, composiciones y velocidades de flujo que resultan
para las corrientes de vapor y líquido en equilibrio para una presión especiada del tanque flash.
Figura 1. Tanque flash

Excepto para casos en los que la volatilidad relativa sea muy grande, el grado de separación de los
componentes que se puede alcanzar en una sola etapa es bajo. Por lo tanto, la operación flash es
generalmente una operación auxiliar para la preparación de corrientes de alimentación que sufrirán
un tratamiento posterior. No obstante, los métodos utilizados en el cálculo de una sola etapa son de
gran importancia.
Cuando la presión de una corriente líquida, donde se conoce la composición, velocidad de flujo y
temperatura (o la entalpía), se reduce adiabáticamente a través de una válvula, se puede realizar un
cálculo de flash adiabático para determinar la temperatura, composiciones y velocidades de flujo que
resultan para las corrientes de líquido y vapor en equilibrio para una presión del tanque flash. Como
en la figura 2., se elige la temperatura TV como variable a iterar, cuyo valor se supone. Se determinan
después ψ, x, y y L como si fuera un flash isotérmico. El valor que se supone de T V (igual a TL) se
contrasta mediante un balance de entalpia para Q=0 para obtener:
𝛙𝑯𝑽 + (𝟏 − 𝛙)𝑯𝑳 − 𝑯𝑭
𝒇(𝑻𝑽 ) = =𝟎 (𝟔)
𝟏𝟎𝟎𝟎
La división es por el orden de magnitud. El método de falsa posición se puede utilizar para iterar
|𝑇𝑉𝑘+1 −𝑇𝑉𝑘 |
sobre TV hasta que se cumpla < 0.0001. El procedimiento es más fácil cuando las mezclas
𝑇𝑉𝑘
forman soluciones ideales, o cuando los valores de K se pueden estimar a partir de 𝐾 = 𝑣𝐿
Caso de estudio
Como se vio anteriormente, el uso de la separación flash permite enriquecer las corrientes de salida
para cada fase del componente de interés, de manera habitual es común encontrar este tipo de
operaciones unitarias en la industria petroquímica, donde a menudo se presentan problemas del
siguiente tipo:
Una corriente formada por seis hidrocarburos que se indica después se somete a flash adiabático. La
entalpía de la alimentación HF es 13,210 Btu/h y la presión del tanque flash es de 300 psia.
Tabla 1. Composición de alimentación a tanque flash.
FRACCIÓN
COMPONENTE
MOL
C2H4 0.60
C2H6 0.05
C3H6 0.10
C3H8 0.03
NC4 0.02
IC4 0.20
Dicha corriente tiene un flujo molar de alimentación de 100 lbmol/hr.
Calcule la temperatura, composiciones, flujos en la corriente de líquido y vapor del tanque flash
empleando la ecuación de estado Soave-Redlich-Kwong.
Metodología.
Figura 2. Algoritmo para cálculo de flash adiabático de mezclas con amplio intervalo de
ebullición.
El algoritmo de la figura 2, es conveniente cuando la ecuación de Rachford-Rice no es sensible a TV.

Es el caso de las mezclas con amplio intervalo de ebullición, por lo tanto, el algoritmo puede fallar si
la ecuación de Rachford-Rice es extremadamente sensible al valor de TV. En este caso es preferible
seleccionar ψ como variable de tanteo y obtener TV de forma iterativa con la ecuación de Rachford-
Rice.
𝑪
𝒛𝒊 (𝟏 − 𝑲𝒊 )
𝒇(𝑻𝑽 ) = ∑ =𝟎 (𝟕)
𝟏 + 𝛙(𝑲𝒊 − 𝟏)
𝒊=𝟏
Después se obtienen x e y a partir de las ecuaciones (8) y (9) respectivamente.

𝒛𝒊
𝒙𝒊 = (𝟖)
𝟏 + 𝛙(𝑲𝒊 − 𝟏)
𝒛𝒊 𝑲𝒊
𝒚𝒊 = = 𝒙𝒊 𝑲𝒊 (𝟗)
𝟏 + 𝛙(𝑲𝒊 − 𝟏)
Para luego obtener directamente ψ de (6) teniendo en cuenta que:
𝑪
𝛙𝑯𝑽 + (𝟏 − 𝛙)𝑯𝑳 − 𝑯𝑭
𝒇(𝛙) = ∑ =𝟎 (𝟏𝟎)
𝟏𝟎𝟎𝟎
𝒊=𝟏
A partir de la cual:
𝑯𝑭 − 𝑯𝑳
𝛙= (𝟏𝟏)
𝑯𝑽 − 𝑯𝑳
Si ψ se obtiene a partir de la ecuación anterior no es igual valor de ψ que se supuso al inicio para
resolver (7), ahora el nuevo valor obtenido de ψ se introduce en (7).
En estos casos, la ecuación de Rachford-Rice y (6) no pueden ser muy sensibles a ψ y TV, y ninguno
de los anteriores procedimientos de tanteo pueden converger. Por lo tanto, es necesario combinar (8)
y (9) con (6) para obtener:
𝛙𝑯𝑽 {𝑻𝑽, 𝛙}+(𝟏 − 𝛙)𝑯𝑳 {𝑻𝑽 , 𝛙} − 𝑯𝑭

𝒇𝟐 (𝛙, 𝑻𝑽 ) = =𝟎 (𝟏𝟐)
𝟏𝟎𝟎𝟎
Y resolver esta ecuación no lineal simultáneamente con la ecuación de Rachford-Rice, también no

lineal en la forma:
𝑪
𝒛𝒊 [𝟏 − 𝑲𝒊 {𝑻𝑽 }]
𝒇𝟏 {𝛙, 𝑻𝑽 } = ∑ =𝟎 (𝟏𝟑)
𝟏 + 𝛙[𝑲𝒊 {𝑻𝑽 } − 𝟏]
𝒊=𝟏
Resultados
Tabla 2. Resultados Temperatura y relación vapor sobre alimentación (ψ).
TEMPERATURA Ψ
15.5590 °C 0.1895
Tabla 3. Composiciones de líquido y vapor.
COMPONENTE X Y
C2H4 0.55709228 0.78348172
C2H6 0.0420511 0.08399106
C3H6 0.11079473 0.05383966
C3H8 0.03183728 0.02214342
NC4 0.02365858 0.00435521
IC4 0.23456603 0.05218893
Tabla 4. Flujos de la corriente líquido y vapor.
LÍQUIDO (lbmol/hr) 81.0469

VAPOR (lbmol/hr) 18.9531
Conclusiones
Dados los resultados anteriores mostrados en las tablas, podemos determinar que el modelado de un
tanque de separación flash adiabático a partir del modelo termodinámico que propusieron Soave-
Redlich-Kwong, el cual predice de manera acertada la distribución de los componentes en las fases
líquido y vapor.
De esta manera se puede determinar que el modelo termodinámico usado predice de manera adecuada
las mezclas multicomponentes, en especial aquellas que contienen hidrocarburos. También queda
clara la ventaja del uso de modelos matemáticos fundamentados en principios termodinámicos y para
la descripción de operaciones unitarias.
Referencias
Seader, J., Henley, E. and Roper, D. (2010). Separation process principles: chemical and biochemical
operations. 3rd ed. Hoboken, NJ: John Wiley & Sons, Inc., pp. 312-315.
Grace F.Chou, J. M. (1986). Adiabatic flash calculations for continuous or semicontinuous mixtures
using an equation of state. Fluid Phase Equilibria, 75-82.
Doherty, M. F. and J. D. Perkins. "On the Dynamics of Distillation Processes I. The Simple
Distillation of Multicomponent, Non-Reacting, Homogeneous Liquid Mixtures." Chemical
Engineering Science.

Reporte FlashAdiabatico

Cargado por

Copyright:

Formatos disponibles

Reporte FlashAdiabatico

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Reporte FlashAdiabatico

Cargado por

Copyright:

Formatos disponibles

Modelado de Flash adiabático empleando la

ecuacion de estado Soave-Redlich-Kwong

M. A. Martínez-Moreno a, J. D. Ponce-Rocha a, J. A. Sánchez-Gómez a, O. Sotelo-Pineda a.

Keywords—Flash Adiabático, Modelo termodinámico, Soave-Redlich-Kwong, Rachford-Rice,

𝛋 = 𝟎. 𝟒𝟖𝟎 + 𝟏. 𝟓𝟕𝟒𝛚 − 𝟎. 𝟏𝟕𝟔𝛚𝟐 (𝟓)

En la actualidad el uso de la ecuación de estado desarrollada por Soave-Redlich-Kwong es empleada

Figura 1. Tanque flash

El algoritmo de la figura 2, es conveniente cuando la ecuación de Rachford-Rice no es sensible a TV.

Después se obtienen x e y a partir de las ecuaciones (8) y (9) respectivamente.

𝛙𝑯𝑽 {𝑻𝑽, 𝛙}+(𝟏 − 𝛙)𝑯𝑳 {𝑻𝑽 , 𝛙} − 𝑯𝑭

Y resolver esta ecuación no lineal simultáneamente con la ecuación de Rachford-Rice, también no

Tabla 2. Resultados Temperatura y relación vapor sobre alimentación (ψ).

Tabla 3. Composiciones de líquido y vapor.

Tabla 4. Flujos de la corriente líquido y vapor.

LÍQUIDO (lbmol/hr) 81.0469

También podría gustarte