El principal inconveniente en el análisis de columnas de destilación multicomponente, más allá de determinar el grado de fraccionamiento que se requiere para cada uno de los compuestos que conforman la alimentación, es calcular las variables de operación (presión y temperatura) en el tope, fondo y en las corrientes límites de batería de la torre. Una secuencia normal de pasos para el diseño de éste tipo de torres implica: 1. Especificar la separación deseada 2. Fijar la presión de la columna 3. Determ Determinación inación de variables variables de diseño diseño a) Determinación del reflujo mínimo. b) Número de platos y etapa de alimentación óptima. c) Optimización del refluj reflujo. o. Éstas condiciones de diseño pueden desarrollarse a partir de métodos aproximados, tales como: a) Método Fenske – Underwood y Gilliland, el cual determina el nº de etapas y el reflujo mínimo y de operación. b) Krenser, Krenser, se emplea para torres de absorción, extracción líquido líquid o – líquido. c) Edmister, para destilación multicomponente con alimentaciones y salidas laterales múltiples. 4. Cálculo riguroso de la columna Para el cual se requiere simuladores o paquetes informáticos para determinar la distribución de cada especie en todos los platos de la columna.
1. ESPECIFICAR LA SEPARACIÓN DESEADA. Lo normal es que se disponga de especificaciones para el producto final. En una columna simple y sin extr extracciones acciones laterales, laterales, se tienen los balances de materia ya conocidos: F=B+D F·(XF) = B·(XB) + D· (XD) En el caso de una destilación multicomponente se tiene un balance global y ecuaciones de balance por componentes, tantas como componentes menos uno tenga la mezcla. Este sistema resulta indeterminado debido a que aparecen más incógnitas que ecuaciones. Se deben fijar las incógnitas necesarias para poder resolver el sistema. Normalmente una de ellas suele ser la alimentación . También es posible especificar la separación indicando el porcentaje de la alimentación que sale por cabeza o por fondo, o, el porcentaje de un determinado componente componente que sale en alguno de los productos. Lo más frecuente es dar dos especificaciones de composición: la de un componente en el destilado y la de otro en el fondo. De esta manera el resto de composiciones vienen fijadas por el balance antes expuesto. Componentes Clave: Son aquellos compuestos para los que se especifica la separación, es decir, quienes delimitan la distribución de destilado y fondo.
Componente Clave ligero (LK)
Componente Clave pesado (HK)
Es el más volátil de ambos ó el más pesado entre los componentes livianos
Es el menos volátil de ambos ó, el más volátil de los pesados
Componentes
Clave
Distribuido”:
Es aquel componente cuya volatilidad se encuentra entre la volatilidad de los componentes clave, ó, en dado caso, cuya volatilidad es muy similar a la volatilidad de uno de los componentes claves. “No
2.- . FIJAR LA PRESIÓN DE LA COLUMNA Una columna puede trabajar a presión, a vacío o a presión atmosférica. Lo más económico resultaría operar a presión atmosférica, pero esto no es siempre posible Presión de Vacío
-La temperatura temperatura produce descomposición de los componentes, ensuciamiento o polimerización. -La alta temperatura en el hervidor complica el aporte de calor.
Sobre presión
Se va a operar a presión en los casos en que resulte interesante refrigerar refrigerar con agua o con aire en el condensado
En el caso de que se tenga que trabajar a vacío se debe elegir la mayor presión posible, pues a medida que aumenta el vacío aumenta el diámetro de la columna y por tanto su coste. Si se trabaja a sobrepresión, se debe elegir la mínima posible ya que el aumento de la presión produce una disminución de la volatilidad y dificulta la separación. Un procedimiento empleado para la estimación de la presión es el de Seader y Henley. Este método también permite la elección del tipo de condensador. Lo más conveniente es que la presión en cabeza sea inferior a 2,86 MPa y la temperatura mínima 50ºC (que corresponde al uso de agua de refrigeración).
Metodología para fijar la l a presión en la columna 1. Estima Estimació ción n de las las compo composic sicion iones es de cabe cabeza za y fondo fondo.. 2. En la ccabez abezaa de la colum columna na se calcula calcula la presión presión de de burbuja burbuja del dest destilado ilado (PD) (PD) a 50ºC. 50ºC. a) Si PD < 1,5 MPa (215 Psia) se usa un condensador total y utiliza agua como refrigerante. b) Si PD > 1,5 MPa ya no se puede utilizar un condensador total y se calcula la presión de rocío del destilado. c) Si PD,r < 2,5 Mpa (365 Psia.) se usa un condensador parcial con agua como refrigerante. d) Si PD,r > 2,5 MPa se busca un refrigerante para trabajar a 2,86 MPa como condensador parcial. 3. Se estima la presión en el rehervidor teniendo en cuenta que: La pérdida de carga entre destilado y tope de la torre son 5 Psia ó 35 Kpa. La pérdida de carga en la columna son 5 Psia o 35 Kpa.
Punto de Burbuja: Punto de Rocío:
Σ Yi = Σ ki * Xi = 1
→
Σ Xi = Σ Yi/Ki = 1
→
Donde Ki = Pºi/Ptotal Pºi/Ptotal
Metodología para fijar la l a presión en la columna Si se conoce la presión en un determinado plato de la torre (data empírica) se puede determinar la presión en cualquier punto de la misma, asumiendo las siguientes consideraciones:
Columnas atmosfér atmosférica ica y de sobrepr sobrepresión esión
Columnas a vacío
Se asume una caída de presión de 0,1 Psi (0.7 KPa) por cada plato
Se asume una caída de presión de 0,05 Psia (0.35 KPa) por cada plato
4. La condición térmica de la alimentación se determina por medio de un cálculo de Flash adiabático para para una presión de plato de alimentación igual a P = PD + 7,5 Psia (50 Kpa.)
Σ XF*(1-Ki) = 0 1+ V (Ki-1)
Con esto se conoce la fracción vaporizada
3.- METODO DE FENSKE-UNDERWOOD-GILLIL FENSKE-UNDERWOOD-GILLILAND AND 1.- ECUACION ECUACION DE FENSKE Se emplea para calcular el número mínimo de etapas de equilibrio (Nmin). El Nmin aumenta a medida que la separación se hace más difícil, es decir, para flujos o % de recuperación relativamente relativamen te grandes del clave liviano o valores de volatilidad relativa relativa cercanos a uno entre los componentes. El Nmin es independiente de la condición térmica de la alimentación y sólo depende del grado de separación deseado y de la volatilidad relativa.
METODO DE FENSKE-UNDERWOOD-GILL FENSKE-UNDERWOOD-GILLILAND ILAND 1.- ECUACION ECUACION DE FENSKE La ecuación de Fenske da una rápida estimación para las Etapas Mínimas requeridas para la operación, es decir, cuando la torre trabaja a Reflujo Total, si se considera la volatilidad relativa constante:
∗ = ∝ α
(LK/HK) promedio
Volatilidad del clave ligero en comparación al pesado
∝/= ∝∗ ∝ ∝/= 3 ∝∗ ∝∗ ∝
∝ ∝ ≤0,1 ∝ ∝ ∝/ ∝ 2
Si se cumple esta relación la volatilidad relativa se considera constante y se puede utilizar la ecuación de Fenske
a la Temperatura de la alimentación
Ecuación equivalente para el cálculo del Número mínimo de etapas
= ∝
LKD → Es la fracción del clave liviano recuperado en el producto de tope con respecto al flujo de clave livia liviano no alimentado a la torre.
LKB→ Es la fracción del clave liviano recuperado en el producto de fondo de la torre con respecto al flujo de clave liviano alimentado.
HKD → Es la fracción del clave pesado en el producto de tope con respecto al flujo de clave pesado alimentado a la torre.
HKB→ Es la fracción del clave pesado recuperado en el producto de fondo de la torre con respecto al flujo de clave pesado alimentado.
Distribución de cada componente “No Clave” a Reflujo total Para determinar la distribución o la concentración de otros componentes en el destilado y en los fondos en reflujo total, la ecuación de Fenske se puede reordenar y escribirse para cualquier otro component componente e de la siguiente forma :
= 1 (∝ )
∗ ∝ = 1 (∝ )
determinadas a partir de la ecuación de Fenske Nmin → Etapas teóricas determinadas
dHK →Flujo molar del clave pesado en el destilado. bHK → Flujo molar del clave pesado en el fondo de la torre. Fi→ Flujo molar del elemento analizado en la alimentación. α(i/HK) → Volatilidad promedio de un elemento en comparación al clave pesado.
EJERCICIOS 1.-Un líquido se alimenta a una torre de destilación a 450 kPa (4.4 Atm.) absolutos. Su composición en fracciones molares es como sigue: n-butano (xA = 0.40), npentano (xB = 0.25), n-hexano (xc = 0.20) y n-heptano (x0 = 0.15). Calcule el punto de ebullición y el vapor en equilibrio con el líquido. 2.-Esta alimentación líquida de 800 mol/h en el punto de ebullición se pasa a una torre de destilación y se fracciona de modo que el 90% del n-pentano (B) se recupera en el destilado y 90% del nhexano (C) en el residuo. Calcule: a) Flujo y composición de los componentes claves en el destilado y del residuo. b) Determine La temperatura superior del destilado (rocío) y el punto de ebullición del residuo. c ) Las etapas mínimas para el reflujo total y la distribución de otros componentes en el destilado y en el residuo.
2.- ECUACION DE UNDERWOOD (Reflujo mínimo Rmin) En la destilación binaria sólo ocurre un “punto comprimido” en el cual el número de etapas se vuelve infinito, y esto suele ser el plato de alimentación. En la destilación de multicomponentes ocurren dos puntos “comprimidos” o zonas de composición constante: uno en la sección por arriba del plato de alimentación y otro por debajo del mismo. El método corto de Underwood para calcular Rmin emplea valores a promedio constantes y también supone flujos constantes en ambas secciones de la torre. Este método brinda un valor razonablemente preciso. Las dos ecuaciones para determinar la relación de reflujo mínimo son:
∝ ∗ 1= ∝ ∝ ∗ R1= ∝
αi → Volatilidad a la temperatura de la alimentación. XIF→ Concentración molar de cada elemento en la alimentación. XID→ Concentración molar de cada elemento en el Destilado.
El valor de Ѳ se encuentra entre el valor α del clave ligero y el valor α del clave pesado, que es 1.0
3.- CORRELACIÓN DE GILLILAND (NÚMERO DE ETAPAS IDEALES) La correlación de Gilliland suministra una buena aproximación para determinar el número de platos ideales que se requieren para el fraccionamiento analizado. Se puede emplear para cualquier condición termodinámica de la alimentación, no obstante, otras correlaciones (Ebbar Maddox) tiene mejores aproximaciones para alimentaciones alimentaciones de Líquido L íquido saturado. saturado. Emplea el Nmin calculado por Fenske y la relación de reflujo mínimo determinado por Underwood.
Correlaciona
1 1 F
R/Rmin
Agua
1.25
Refrigerante
1.40 – 1.50
Supertorre
1.10
3.- CORRELACIÓN DE ERBAR MADDOX (NÚMERO DE ETAPAS IDEALES) La correlación de Erbar Maddox suministra una buena aproximación para determinar el número de platos ideales que se requieren para el fraccionamiento de un líquido en su punto de burbuja (líquido saturado). Correlaciona gráficamen gráficamente: te:
1
Puede emplearse para Nmin calculado por Fenske o Winn, aunque da mejores aproximaciones para Winn; y la relación de reflujo mínimo determinado por Underwood.
Las curvas corresponden a la relación de reflujo mínimo (Rm/Rm+1). Se intercepta la curva que corresponde al valor de ésta relación, con el valor del eje Y que representa (R/R+1)
Si la volatilidad relativa es variable a lo largo de la torre, no se puede utilizar la ecuación de Fenske para determinar el número mínimo de etapas requeridas para el fraccionamiento, fraccionamiento, para ello se utiliza la ecuación de WINN:
=
−− ∗ −−
() =
LK clave ligero HK clave pesado D Destilado B Fondo de la torre
Se calcula un β tanto para el fondo como para el destilado y se igualan las expresiones. De ésta forma se consiguen las dos incógnitas. Método aproximado, se toman las constantes de equilibrio en la alimentación.
Distribución de cada componente “No Clave” Para determinar la distribución o la concentración de otros componentes en el destilado y en los fondos en reflujo total, la ecuación de Winn se puede reordenar y escribirse para cualquier otro componente componente de la siguiente siguiente forma:
=
( ) 1 ∗ −
=
PLATO ÓPTIMO DE ALIMENTACIÓN En el caso de problemas donde se especifica es pecifica o determinan los flujos de destilado y el número de etapas teóricas requeridas, es difícil especificar la localización de la etapa de alimentación que conduce al máximo grado de separación posible. Sin embargo, se han desarrollado diferentes métodos gráficos y analíticos para su identificación, ente los cuales vale destacar: 1.- Construir una gráfica modificada de McCabe-Thiele basada en los componentes claves, para determinar si la etapa de
alimentación está localizada en la etapa óptima o, si por el contrario, debería trasladarse. Para esta gráfica se calculan las fracciones molares de los componentes claves sobre la base exenta de los no clave.
2.- Ecuación de Kirkbride
, = ∗ ∗
Para esta ecuación nos queda un sistema de Para ecuaciones lineales sencillos, tomando en cuenta que:
= ó
Etapa de alimentación es igual a N
¿Qué hacer cuando estemos en presencia de un elemento e lemento “No “ No clave clave distribuido” Un componente “No clave distribuido” es cualquier elemento intermedio entre los elementos claves, es decir, que posee volatilidades intermedia entre los elementos bajo los cuales se realiza el fraccionamiento. Un ejemplo de éstos sería el caso en el cual se desea separar etano (C2) de butano (C4) y en la alimentación exista propano (C3); irremediablemente, el propano también se fraccionará fraccionará en la corriente de destilado y de fondo de la torre. Según la metodología aplicada, la existencia de un “No clave distribuido” afecta la relación de reflujo de la torre, de tal forma que debemos calcular dos θ diferentes en la ecuación de Underwood, uno que sería entre el α del clave pesado y el distribuido, y otro entre el α del distribuido y el clave liviano.
1= ∝∝∗
Se calcula dos veces ésta expresión y obtendremos dos valores de Ѳ distintos.
Una forma de determinar a ciencia cierta que un elemento es Distribuido, es mediante la ecuación de Shiras, si el valor de ésta expresión es entre 0 – 1, indica que ese elemento se distribuye en la torre.
= ∝∝ 11 ∗ ∗∗ ∝∝ ∝1 ∗ ∗∗
Por supuesto, aplica la ecuación para aquellos que usted crea que se distribuyan nada más, no
DIÁMETRO DE LA COLUMNA El diámetro de una columna de destilación multicomponente puede estimarse mediante la ecuación de Heaven, la cual indica que:
=
∗+∗ ∗ = +∗ = ∗∗ ∗ ∗
TDV es la temperatura de rocío del vapor en el condensador (K) D es el flujo molar del destilado (Kmol/hora) P Presión promedio de la columna (Atm.) ρ Densidad promedio del vapor vapor en la columna (Kg/m (Kg /m 3) R Relación de reflujo de operación.
ALTURA DE LA COLUMNA
=0,61∗ 4,27
N es el número de platos ideales Η es la eficiencia promedio de los platos de la columna
