Métodos aproximados para el diseño de columnas de destilación multicomponente

1. TRANSFERENCIA DE MASA IIMétodos aproximados y rigurosos para el diseño de columnas de destilación de sistemas multicomponentesMaría del Rosario Caicedo Realpe, Ing. Química, M. Sc.e-mail: rcaicedo@udea.edu.co, Of. 18 - 411

2. Métodos aproximados para el diseño de columnas de destilación multicomponente 1. Método de grupo de Smith - Brinkley (SB) 2. Método de Fenske - Underwood - Gilliland (FUG) 3. Método de grupo de Kremser

3. 1. Método de grupo de Smith - Brinkley (SB) Características - Útil para evaluar el comportamiento de columnas ya existentes - Funciona bien cuando se tiene un rehervidor parcial o hay una entrada de calor por el fondo - Se debe conocer el número de platos totales - Sirve para encontrar el plato óptimo de alimentación cuando se especifican los requerimientos de recuperación - Sirve para calcular las composiciones de los productos cuando se han especificado el plato de alimentación y el flujo de los productos - El empleo de valores de k eficaces es muy útil cuando los coeficientes de distribución varían mucho a lo largo de la columna

4. Algoritmo para el caso en que se conoce el plato de alimentación y se requieren las composiciones de los productos Datos del problema 1. Número de platos totales de la columna, N (si el condensador es parcial se le suma 1) 2. Número total de platos de la zona de despojamiento (incluyendo el rehervidor parcial), M 3. Relación de reflujo, R 4. Relación de reflujo interna para la zona de enriquecimiento, L / V * 5. Flujo del destilado, D 6. Especificación total del alimento (Pf, Tf, zf, F) * L y V son valores eficaces que se van a considerar constantes, al igual que L’ y V’ para la zona de despojamiento

5. Procedimiento 1. Con R y L / V se calculan L y V 2. Con L y V se calculan L’ y V’ (cómo?) 3. Se supone una composición para D y W y se calculan TD y TW (TD es la Tr del destilado, y TW es la temperatura de burbuja de los fondos) 4. Se calculan las temperaturas eficaces, Tn y Tm, para las zonas de enriquecimiento y despojamiento, respectivamente: Tn = (T1 + Tpf)/2 y Tm = (TN+1 + Tpf)/2 (Tpf: Temp. aliment.) 5. Se hallan ki = f(Tn) y ki’ = f(Tm) 6. Se calculan los factores de separación, Sni y Smi, para las zonas de enriquecimiento y despojamiento, respectivamente: Sni = kiV/L Smi = ki’V’/L’

6. 7. Si la alimentación es básicamente líquido: hi = ((ki’ L) / (ki L’)) (1 - Sni ) / (1 - Smi) Si la alimentación es básicamente vapor: hi = (L / L’) (1 - Sni ) / (1 - Smi) 8. Se determina la fracción del componente i que sale por los fondos: fi = (W xWi) / (FzFi) 9. Se compara el valor de fi calculado en (8) con:: ((1 + SniN-M) + R (1 - Sni)) / ((1 - SniN-M) + R (1 - Sni) + hi SniN-m (1 - SmiM+1)) Si los dos valores no son iguales se regresa al paso 3 y se repite el procedimiento hasta que los valores se han iguales en la cuarta cifra decimal.

7. Conceptos importantes usados por los métodos de cálculo de columnas de destilación multicomponente Clave pesado (hk) Es el componente de menor volatilidad relativa entre los compuestos que salen por el destilado y que se halla en una cantidad considerable (mayor al 2%) Clave ligero (lk) Es el componente de mayor volatilidad relativa entre los compuestos que salen por el fondo y que se encuentra en una cantidad considerable (mayor al 2%) Componente de referencia El componente clave pesado es el componente de referencia Cómo es la volatilidad relativa de los componentes más pesados que el clave pesado?, y la de los componentes más ligeros que el clave ligero?

8. Verificación de la elección de los componentes clave A partir del hecho de que se requiere un número infinito de etapas para excluir totalmente del destilado los componentes más pesados que el clave pesado, y así mismo, de los fondos los más ligeros que el clave ligero, Shiras et al (Ind. Eng. Chem., 42, 871 (1950)) mostró que para una relación de reflujo mínima Rm se debe cumplir: Ai = yDi / zFi = ((ai - 1)/(alk - 1)) (yDlk / zFlk) + ((alk - ai)/(alk - 1)) (yDhk / zFhk) - Si A < -0.01 o A > 1.01 probablemente el componente i no se distribuirá - Si 0.01 < A < 0.99 seguramente el componente i se distribuirá Los componentes clave se distribuyen entre el destilado y los fondos Qué sucede con los demás componentes?

9. 2. Método de Fenske - Underwood - Gilliland (FUG) Características - Útil fundamentalmente para el diseño de columnas, aunque también sirve para evaluar el comportamiento de columnas ya existentes - Parte del conocimiento de la localización óptima del plato de alimentación - Supone flujos molares constantes (el Rm calculado es mucho menor que el Rm real) - Sirve para calcular los flujos y las composiciones de los productos cuando se han especificado los requerimientos de recuperación de uno o máximo dos componentes - Se pueden determinar las temperaturas del tope y del fondo de la columna

10. Ecuación de Fenske Determina el valor de Nm (a reflujo total). Nm incluye el rehervidor parcial, si la columna tiene un condensador parcial se le suma 1. (xi / xhk) = aiNm (xWi / xWhk) ó Nm = log((DyDlk / WxWlk) (WxWhk / DyDhk)) / logalk alkes un valor eficaz alk = (aNlk aN-1lk ...a2lk a1lk)1/N ó alk = (aDlk apflk aWlk)1/3 ó alk = (aDlk aWlk)1/2 Como una aproximación se puede tomar Nm como el 40 o 60% de N y suponiendo separaciones de los componentes iterativamente se pude llegar a estimar los valores iniciales de las temperaturas en los extremos

11. Ecuaciones de Underwood Es útil para determinar Rm, suponiendo ai eficaz y L/V constantes. Se deben resolver las siguientes ecuaciones: S ((ai zFi)/(ai - f)) = 1 - q (1) i son todos los componentes S ((aj yDj)/(aj - f)) = Rm + 1 (2) j son el hk y los componentes más ligeros que el pesado (?) De (1) con ahk < f < alk se resuelve (2).Si existen B componentes entre los claves el número de valores para f es B+1. (2) va a permitir construir un sistema de B+1 ecuaciones con B+1 incógnitas, (siempre y cuando se hayan especificado las composiciones de los componentes clave en el destilado), cuáles? Si yDj < 0 ó yDj > 1, los componentes clave elegidos son incorrectos.

12. Con las ecuaciones anteriores también se pueden estimar las temperaturas iniciales para el destilado y los fondos, así: 1. ai eficaz se puede asumir una temperatura promedio de la columna o tomar la temperatura del alimento para iniciar el procedimiento 2. Con las ecuaciones de Underwood se calculan los yDi’s con ai 3. Se encuentra la temperatura de rocío, T0, para el destilado con los valores de yDi’s 4. Con balances de materia se hallan los valores de xWi’s (La relación D/F debe tenerse como dato de entrada) 5. Se calcula la temperatura de burbuja, TN+1, para los fondos con los valores de xWi’s 6. Se compara el valor promedio de temperatura asumido en 1 con (TN+1 + T0)/2, sí son diferentes se vuelve al punto 1 con este nuevo valor promedio

13. 1.0 (N - Nm)/(N+1) 0. 1 Ec. Molokanov 0.01 1.0 0.1 0.01 (R - Rm)/(R+1) Correlación empírica de Gilliland Conocidos Nm y Rm, con la correlación de Gilliland, para un valor de R especificado se puede calcular N (o lo contrario) Ec. Molokanov (Int. Chem. Eng., 12(2), 209 (1972)) G = 1 - exp((1+54.4H)(H-1)/((11+117.2H)H0.5)) Donde: G = (N - Nm)/(N + 1) y H = (R - Rm)/(R+1)

14. Y ahora qué?

15. Localización del plato de alimentación Se considera solamente el sistema binario formado por los componentes clave. Partiendo del Balance del componente más liviano - zona de enriquecimiento (presentación de McCabe-Thiele) yn+1 V = xn L + yD D Re-escribiendo (2) para lk y hk y despejando x se tiene (sin considerar plato): xlk = ylk V/L - yDlk D/L xhk = yhk V/L - yDhk D/L Eliminando L de las ecuaciones anteriores y despejando ylk se tiene: ylk = (xlk/xhk)(yhk - yDhkD/V) + yDlkD/V (1)

16. Partiendo del Balance del componente más liviano - zona de despojamiento (presentación de McCabe-Thiele) xN-3 L’ = yN-2 V’ + xW W y aplicando el mismo procedimiento anterior se tiene: ylk = (xlk/xhk)(yhk + xWhkW/V’) - xWlkW/V’ (2) En el punto de alimentación las ecuaciones (1) y (2) son iguales, entonces: (xlk/xhk)int. = (yDlkD/V + xWhkW/V’)/(yDhkD/V + xWlkW/V’) Considerando R y q, la ecuación anterior se puede re-escribir como: (xlk/xhk)int. = (zFlk - yDlk(1-q)/(R+1))/(zFhk - yDhk(1-q)/(R+1)) La razón (xlk/xhk) se debe calcular plato a plato y la localización del plato de alimentación se obtiene cuando: (xlk/xhk)f-1> (xlk/xhk)int.> (xlk/xhk)f

17. Actividad Revisar el método de grupo de Kremser

19. Métodos rigurrosos para el diseño de columnas de destilación multicomponentes 1. Método de Lewis y Matheson 2. Método de Thiele y Geddes Actividad Revisar el ejemplo 9.13 de Treybal.

20. En general, para el diseño preliminar de una columna de destilación, Información que siempre se debe tener: 1. Datos de equilibrio 2. Eficiencia de los platos 3. Condiciones del alimento 4. Presión en el destilado (considerando fluido de enfriamiento usado en el condensador) De los siguientes parámetros el diseñador sólo especificará tres: 1. N 2. R 3. R’ (=VN+1 / W) 4. Concentración de un componente o máximo 2 en un producto 5. Separación de un componente o máximo 2 (DyDi / WxWi) 6. D / W.