Proyecto - Destilación Binaria
Proyecto - Destilación Binaria
Proyecto - Destilación Binaria
UNIVERSIDAD DEL VALLE ESCUELA DE INGENIERA QUMICA SANTIAGO DE CALI - VALLE 2011 DESTILACIN BINARIA DEL SISTEMA CICLOHEXANONA-CICLOHEXANOL, CARACTERIZACIN DEL PROCESO Y DE LA MEZCLA A SEPARAR: 1
UNIVERSIDAD DEL VALLE ESCUELA DE INGENIERA QUMICA SANTIAGO DE CALI - VALLE 2011
Contenido
I. II. 1 2 ndice de figuras .......................................................................................................................... 5 ndice de Tablas ........................................................................................................................... 7 Introduccin ................................................................................................................................ 8 Objetivos ................................................................................................................................... 10 2.1 2.2 3 Objetivo General ............................................................................................................... 10 Objetivos Especficos ......................................................................................................... 10
4 Caracterizacin del proceso, Balances de materia, presin y temperatura de cada flujo del proceso. ............................................................................................................................................. 12 5 6 7 Descripcin del equipo seleccionado y condiciones de operacin. .......................................... 13 Equilibrio Lquido-Vapor para la mezcla Ciclohexanona /Ciclohexanol. ................................... 15 Mtodos grficos para el clculo de etapas ideales ................................................................. 21 7.1 7.2 8 Mtodo de McCABE-THIELE .............................................................................................. 21 Mtodo de PONCHON-SAVARIT ........................................................................................ 22
Elaboracin del cdigo y clculo de una torre .......................................................................... 23 8.1 8.2 8.3 McCABE-THIELE ................................................................................................................. 23 PONCHON SAVARIT ........................................................................................................... 32 Eficiencia Global de la torre .............................................................................................. 36
Variacin de condiciones de operacin .................................................................................... 37 9.1 McCABE-THIELE ................................................................................................................. 38 Lquido Saturado Lquido-Vapor Lquido-Vapor Vapor Saturado ........................................................................................... 38 ........................................................................................... 40 ........................................................................................... 42 ............................................................................................. 44
PONCHON-SAVARIT........................................................................................................... 46 Lquido Saturado Lquido-Vapor Lquido-Vapor Vapor Saturado ........................................................................................... 46 ........................................................................................... 48 ........................................................................................... 50 ............................................................................................. 52 3
9.3 10 11
I.
ndice de figuras
Figura 1. Grupo de reacciones presentes en el proceso de hidrogenacin de fenol. .......................... 9 Figura 2. Diagrama del Proceso de destilacin de Ciclohexanona. ............................................... 11 Figura 3. Torre de destilacin T-2. ................................................................................................. 13 Figura 4. Del ajuste (a) Diagrama T-x-y y (b) Diagrama x-y. Los puntos son los datos experimentales y las lneas los valores calculados. ........................................................................... 18 Figura 5. Diagrama T-x-y del sistema Ciclohexanona (1) Ciclohecanol (2) a P= 200 mmHg. ......... 19 Figura 6. Diagrama x-y del sistema Ciclohexanona (1) Ciclohecanol (2) a P= 222.5 mmHg. ......... 19 Figura 7. Diagrama H-x-y del sistema Ciclohexanona (1) Ciclohecanol (2) a P= 200 mmHg, con una Temperatura de referencia Tr = 300 K. ............................................................................................. 20 Figura 8. Divisin de la torre de destilacin en secciones ................................................................ 21 Figura 9. Lnea de Alimentacin y lnea de reflujo mnimo. ............................................................. 26 Figura 10. Lnea de rectificacin ....................................................................................................... 28 Figura 11. Lnea de agotamiento ...................................................................................................... 29 Figura 12. Resultado final del cdigo para McCabe-Thiele............................................................... 30 Figura 13. Perfil de temperatura - McCabe-Thiele ........................................................................... 31 Figura 14. Perfiles de composicin - McCABE - THIELE..................................................................... 31 Figura 15. Reflujo Mnimo................................................................................................................. 32 Figura 16. Resultado final del cdigo para Ponchon-Savarit. ........................................................... 35 Figura 17. Perfil de temperatura, mtodo de Ponchon-Savarit........................................................ 35 Figura 18. Perfil de composiciones de Ponchon-Savarit ................................................................... 36 Figura 19. f=0 nRo = 1.1 ................................................................................................................ 38 Figura 20. f=0 nRo = 1.2 ................................................................................................................ 38 Figura 21. f=0 nRo = 1.5 ................................................................................................................ 39 Figura 22. f=0 nRo = 2.0 ................................................................................................................ 39 Figura 23. f=1/3 nRo = 1.1 ............................................................................................................ 40 Figura 24. f=1/3 nRo = 1.2 ............................................................................................................ 40 Figura 25. f=1/3 nRo = 1.5 ............................................................................................................ 41 Figura 26. f=1/3 nRo = 2.0 ............................................................................................................ 41 Figura 27. f=2/3 nRo = 1.1 ............................................................................................................ 42 Figura 28. f=2/3 nRo = 1.2 ............................................................................................................ 42 Figura 29. f=2/3 nRo = 1.5 ............................................................................................................ 43 Figura 30. f=2/3 nRo = 2.0 ............................................................................................................ 43 Figura 31. f=1 nRo = 1.1 ................................................................................................................ 44 Figura 32. f=1 nRo = 1.2 ................................................................................................................ 44 Figura 33. f=1 nRo = 1.5 ................................................................................................................ 45 Figura 34. f=1 nRo = 2.0 ................................................................................................................ 45 Figura 35. f=0 nRo = 1.1 ................................................................................................................ 46 Figura 36. f=0 nRo = 1.2 ................................................................................................................ 46 Figura 37. f=0 nRo = 1.5 ................................................................................................................ 47 5
Figura 38. f=0 nRo = 2.0 ................................................................................................................ 47 Figura 39. f=1/3 nRo = 1.1 ............................................................................................................ 48 Figura 40. f=1/3 nRo = 1.2 ............................................................................................................ 48 Figura 41. f=1/3 nRo = 1.5 ............................................................................................................ 49 Figura 42. f=1/3 nRo = 2.0 ............................................................................................................ 49 Figura 43. f=2/3 nRo = 1.1 ............................................................................................................ 50 Figura 44. f=2/3 nRo = 1.2 ............................................................................................................ 50 Figura 45. f=2/3 nRo = 1.5 ............................................................................................................ 51 Figura 46. f=2/3 nRo = 2.0 ............................................................................................................ 51 Figura 47. f=1 nRo = 1.1 ................................................................................................................ 52 Figura 48. f=1 nRo = 1.2 ................................................................................................................ 52 Figura 49. f=1 nRo = 1.5 ................................................................................................................ 53 Figura 50. f=1 nRo = 2.0 ................................................................................................................ 53 Figura 51. Reflujo Mnimo contra la fraccin de vapor .................................................................... 54 Figura 52. Nmero de etapas ideales en funcin de la fraccin de vapor y el factor de escalamiento del reflujo .......................................................................................................................................... 55 Figura 54. Nmero de platos ideales en funcin de la fraccin de vapor y el factor de escalamiento del reflujo. ......................................................................................................................................... 56 Figura 53. Fraccin de vapor contra el reflujo mnimo. .................................................................... 56
II.
ndice de Tablas
Tabla 1. Balance de materia.............................................................................................................. 12 Tabla 2. Composiciones y temperaturas de los flujos de la torre ...................................................... 13 Tabla 3. Constantes de la ecuacin de Antoine ................................................................................ 16 Tabla 4. Resultados del ajuste de datos experimentales.................................................................. 17 Tabla 5. Constantes para Cp_liq ....................................................................................................... 20 Tabla 6. Constantes para Clculo de la entalpa de vaporizacin (Hvap) ......................................... 20 Tabla 7. Funciones complementarias del cdigo desarrollado ........................................................ 23 Tabla 8. Resultados de los clculos por ambos mtodos. ................................................................ 36 Tabla 9. Reflujo Mnimo en funcin de la fraccin de vapor ............................................................ 54 Tabla 10. Nmeo de etapas ideales mtodo de McCabe-Thiele ...................................................... 54 Tabla 9. Reflujo Mnimo en funcin de la fraccin de vapor ............................................................ 55 Tabla 10. Nmero de etapas ideales Mtodode Ponchon Savarit ................................................... 56 Tabla 9. Desviaciones del reflujo mnimo mtodo aproximado ....................................................... 57 Tabla 10. Desviaciones del Nmero de etapas ideales del mtodo aproximado ............................. 57
1 Introduccin
La ciclohexanona tiene un nmero significativo de usos industriales, entre los que se encuentran su uso como un solvente y diluyente de barnices, especialmente los que contienen nitrocelulosa o cloruro de vinilo. Es un excelente solvente de insecticidas y de muchos otros materiales, tambin es el componente bsico en la sntesis de varios compuestos orgnicos, como lo son farmacuticos, insecticidas y herbicidas. Tambin es utilizada en la fabricacin de cintas magnticas y de vdeo. La gran importancia de la ciclohexanona a nivel industrial, se debe a que es un producto intermedio en la produccin de cido adpico, que puede ser empleado en la fabricacin de nylon-6, 6 un polmero de cido adpico y hexametilendiamina. La ciclohexanona tambin se utiliza como producto intermedio en la produccin de caprolactama, que posteriormente se utilizar en la fabricacin de nylon-6. Otro uso importante es en la fabricacin de steres para el uso de plastificantes como los ftalatos ciclohexilo y di-ciclohexilo. Como estabilizadores de jabones y detergentes sintticos y emulsiones. como disolvente y medio de contraste asistente kierbelling para la industria textil. como disolvente en pinturas, barnices y otras industrias qumicas [2]. Para la produccin industrial de Ciclohexanona, los caminos ms conocidos y utilizados son: Hidrogenacin cataltica de fenol Deshidratacin oxidativa de ciclohexanol Oxidacin cataltica de ciclohexano Donde la ruta ms eficiente entre los procesos enunciados para la sntesis de ciclohexanona es la hidrogenacin cataltica de fenol [3]. La reaccin se lleva acabo en un catalizador de paladio, de manera altamente selectiva hacia la ciclohexanona, pero en este proceso la ciclohexanona es un producto intermedio de algunas reacciones en serie. Como se puede observar en la figura 1, esta no es la nica reaccin que se da en el sistema, una pequea parte de la ciclohexanona sigue hidrogenandose hasta volverse Ciclohexanol, el fenol se puede hidrogenar directamente hasta ciclohexanol, o parte del fenol puede reducirse a benceno, entre otras reacciones de menor trasendencia en el proceso, Aunque como se observa en la figura las reacciones ms importantes son las primeras tres mensionadas. Deacuerdo a lo anterior es claro que la corriente de salida de un reactor de hidrogenacin cataltica de fenol, tiene una mezcla de ciclohexanona y ciclohexanol, junto con el fenol y el hidrogeno que no reaccionaron, ms algunas impureza como hidrocarburos ms pesados. La separacin de estas corrientes se puede realizar por un proceso de destilacin, pero es necesario hacerlo de manera eficiente, o el gasto de energa ser muy alto y la separacin ser muy costosa.
Figura 1. Grupo de reacciones presentes en el proceso de hidrogenacin de fenol. En busca de hacer ms eficiente esta operacin, en 1981, Murthy [1] presenta algunas alternativas para el mejoramiento del proceso de separacin, con un sistema de varias columnas operando de forma continua. Para este trabajo se ha seleccionado una torre de destilacin de uno de los arreglos propuestos por Murthy et al, a la cual se le realiza el clculo de las etapas de equilibrio, con los datos experimentales de una corrida a ciertas condiciones de operacin, informacin reportada por el mismo autor en la patente desarrollada, y a partir de estas etapas ideales y el nmero de platos reales, determinar la eficiencia de la torre. Para el proceso seleccionado, la pureza del destilado es muy alta, la composicin molar de Ciclohexanona en el producto de cabeza es del 99.7 %, composicin a las que resulta poco viable realizar un clculo de la torre por mtodo grfico, pues no se tendra la precisin necesaria, de all que se ha desarrollado un cdigo en la plataforma MATLAB, el cual realiza los clculo equivalentes a cada mtodo grfico, y reporta los resultados de la torre, junto con el grfico que se obtendra si se realizar el clculo de manera manual.
2 Objetivos
2.1 Objetivo General
Realizar el calculo por etapas de equilibrio, de una torre de destilacin para separar una mezcla de Ciclohexanona-Ciclohexanol, y comparar los resultados obtenidos por dos mtodos diferentes de clculo.
10
La corriente que entra al proceso (1) es una mezcla de Ciclohexanona, Ciclohexanol, Fenol, y una muy baja cantidad de algunos compuestos de alto punto de ebullicin (pesados), esta corriente entra en la torre T-1, y es separada en tres corrientes, la primera es el producto de cabeza (2), el cual es prcticamente ciclohexanona pura con un pequeo contenido de ciclohexanol, el destilador de esta torre opera a un reflujo de 3.65, a una presin de 70 torr y retira un calor de 4862 kcal/h. La corriente intermedia de esta torre (3) contiene solamente ciclohexanona y ciclohexanol, y pasa a un proceso de separacin posterior en la torre T-2. El tercer flujo es el producto de fondo (6), contiene todo el fenol y los compuestos pesados que ingresaron en la torre, ms una pequea parte de ciclohexanona y ciclohexanol, esta corriente es alimentada a la torre T-3, donde se somete a otro proceso de separacin del que no se tienen datos especficos. El reboiler en T-1 opera a 17.1 de reflujo, una presin de 362 torr y se introduce un calor de 3118 kcal/h. El condensador de la torre T-2 opera a un reflujo de 4.25, una presin de 70 torr, y retira un calor de 2938 kcal/h, mientras el reboiler opera a un reflujo de 8.0, una presin de 375 torr, e introduce un calor de 1444 kcal/h.
11
4 Caracterizacin del proceso, Balances de materia, presin y temperatura de cada flujo del proceso.
Balance de Materia Tomando como base de clculo 1 hora de operacin Tabla 1. Balance de materia flujo Ciclohexanona (mol/h) Ciclohexanol (mol/h) Fenol (mol/h) Compuestos pesados (mol/h) Total Temperatura (C) flujo Ciclohexanona (mol/h) Ciclohexanol (mol/h) Fenol (mol/h) Compuestos pesados (mol/h) Total Temperatura (C) 1 187,97 8,08 10,63 0,70 207,38 120,00 3 82,51 7,49 0,00 0,00 90,00 114,00 2 100,70 0,30 0,00 0,00 101,00 82,00 4 53,84 0,16 0,00 0,00 54,00 82,00 2a 468,26 1,40 0,00 0,00 469,65 82,00 4 282,66 0,80 0,00 0,00 283,46 82,00 2b 367,56 1,10 0,00 0,00 368,65 82,00 4b 228,82 0,64 0,00 0,00 229,46 82,00 6 4,76 0,29 10,63 0,70 16,38 161,00 5 28,67 7,33 0,00 0,00 36,00 131,00 6 86,16 2,03 61,23 0,90 150,32 161,00 5 258,03 43,98 0,00 0,00 302,01 131,00 6b 81,40 1,74 50,60 0,20 133,94 161,00 5b 229,36 36,65 0,00 0,00 266,01 131,00
12
Figura 3. Torre de destilacin T-2. Respecto a la torre se tiene la siguiente informacin: Tabla 2. Composiciones y temperaturas de los flujos de la torre Alimento Destilado Fondo (F) (D) (F) 3 4 5 Flujo Fraccin molar de Ciclohexanona (1) 0,9168 0,9970 0,7964 Fraccin molar de Ciclohexanol (2) 0,0832 0,0030 0,2036 114,00 82,00 Temperatura (C) 131,00 Condensador Reflujo de 4.25 Presin de 70 torr (mmHg) Calor retirado de 2938 kcal/h Reboiler Reflujo de 8.0 Presin de 375 torr (mmHg) Calor suministrado de 1444 kcal/h.
Nmero de Etapas y plato de alimentacin La torre de destilacin T-2 cuenta con 60 platos, y se alimenta por el plato 20 de la columna. Presin de operacin La presin de operacin de la torre no esta definida, probablemente porque la torre no es isobrica, pero el clculo se puede realizar asumiendo una presin promedio constante, que debe estar entre la presin del condensador y la del reboiler. La presin de operacin de la torre se fija en un valor medio de P, es decir
13
Como se conocen el nmero de platos reales, se va a calcular el nmero de platos ideales por el mtodo de McCabe-Thiele y por el mtodo de Ponchon-Savarit. Se determina la desviacin del mtodo simplificado (McCabe-Thiele) respecto al riguroso (Ponchon-Savarit) y despus la eficiencia de la torre respecto al mtodo riguroso. Este procedimiento se va realizar ahora cambiando el reflujo y para diferentes estados de alimentacin para elaborar un anlisis de sensibilidad en el nmero de platos ideales obtenido por cada uno de los dos mtodos empleados. Por otra parte, no se conoce el grado de sobrecalentamiento en el Rehervidor, por tando se va a asumir un evaporador, el vapor recirculado sale como vapor saturado y se va a dejar de tratar como un rehervidor.
14
: Composicin del componente i en la fase de vapor P: Presin del sistema : Composicin del componente i en la fase lquida : Coeficiente de actividad del componente i. : Presin de saturacin del componente i a la temperatura del sistema. Las presiones de saturacin se calculan con la ecuacin de antoine, y los valores de las constantes son tomadas del artculo de Goodwin y Newsham[6]. Con Psat en mmHg y T en C. 15
Tabla 3. Constantes de la ecuacin de Antoine A B C Ciclohexanona (1) 7.4280 1777.7 236.12 Ciclohexanol (2) 7.0385 1318.5 156.60 Los coeficientes de actividad se determinan con la ecuacin NRTL[8], que para un sistema binario se * ( ) ( )+
)+
El ajuste de datos se realiza calculando valores de y de T a partir de los valores de yP experimentales, mediante un algoritmo burbuja T. Estos valores calculados se comparan con los valores experimentales, se determina la diferencia, y se define las siguientes funciones para evaluar los errores en el clculo. ( ( ) )
Siendo N el nmero total de datos. Finalmente, se plantea la siguiente funcin objetivo, que se debe minimizar para obtener el valor ptimo de los parmetros.
Durante el ajuste de datos, para las presiones 100, 200 y 395 mmHg, se obtiene un buen ajuste con bajos errores en los valores calculado, pero para los datos reportados a presiones de 30 y 750 mmHg, el ajusteno es tan bueno y lo errores son bastante apreciables.
16
Tomando en cuenta lo anterior, y sabiendo que las presiones de operacin son 70 mmHg en el Rehervidor (que se va a asumir como evaporador), 222.5 mmHg en la torre y 375 mmHg en el condensador, se observa claramente que para los calculos que se van a realizar, los datos experimentales a 100, 200 y 395 mmHg, son suficiente para el rango de trabajo, y por tanto, los datos experimentales a 30 y a 750 mmHg no son tomados en cuenta para determinar el valor de los parmetrosde ajuste. Para minimizar la funcin objetivo y obtener el valor ptimo de los parmetros, se utiliza la funcin del paquete de computo MATLAB, los resultados obtenidos del ajuste son
En trminos generales, los resultados del ajuste presentan los siguientes errores promedio y mximo en el conjunto de datos experimentales. Tabla 4. Resultados del ajuste de datos experimentales P(mmHg) 395 0,108 0,228 0,450 0,590 0,831 0,105 0,235 0,438 0,593 0,644 0,807 0,107 0,247 0,434 0,598 0,740 0,813 0,147 0,294 0,519 0,654 0,856 0,165 0,337 0,550 0,681 0,722 0,855 0,185 0,394 0,598 0,713 0,827 0,875 138,8 137,7 135,4 134,3 132,6 119,6 117,7 115,1 113,4 113,0 111,6 102,6 99,7 96,3 94,3 92,8 92,2 0,147 0,292 0,523 0,651 0,858 0,165 0,336 0,551 0,687 0,729 0,856 0,198 0,397 0,596 0,732 0,833 0,881 139,0 137,6 135,5 134,3 132,7 119,9 118,0 115,4 113,7 113,2 111,7 103,0 100,4 97,4 95,2 93,6 92,8 0,000 0,002 -0,004 0,003 -0,002 0,000 0,001 -0,001 -0,006 -0,007 -0,001 -0,013 -0,003 0,002 -0,019 -0,006 -0,006 -0,2 0,1 -0,1 0,0 -0,1 -0,3 -0,3 -0,3 -0,3 -0,2 -0,1 -0,4 -0,7 -1,1 -0,9 -0,8 -0,6
200
100
17
En la Figura 4 se puede observar de manera grfica los resultados del ajuste. En las figuras 5, 6 y 7 se presentan los diferentes diagramas de equilibrio, calculados con estos parmetros a la presin de operacin de la torre.
(a)
(b) Figura 4. Del ajuste (a) Diagrama T-x-y y (b) Diagrama x-y. Los puntos son los datos experimentales y las lneas los valores calculados.
18
Figura 5. Diagrama T-x-y del sistema Ciclohexanona (1) Ciclohecanol (2) a P= 200 mmHg.
Figura 6. Diagrama x-y del sistema Ciclohexanona (1) Ciclohecanol (2) a P= 222.5 mmHg. 19
Figura 7. Diagrama H-x-y del sistema Ciclohexanona (1) Ciclohecanol (2) a P= 200 mmHg, con una Temperatura de referencia Tr = 300 K. Durante el clculo del diagrama H-x-y se utilizaron las siguientes ecuaciones, para y para respectivamente Capacidad Calorfica de lquido [9] ( A Ciclohexanona (1) Ciclohexanol (2) Calor de vaporizacin [9] ( Donde ( )) ( ) Tabla 5. Constantes para Cp_liq B C D 243 298 566 563 )
es la temperatura crtica, y A y n son parmetros para cada especie. Tabla 6. Constantes para Clculo de la entalpa de vaporizacin (Hvap) A (K) N 629.15 Ciclohexanona (1) 53.720 629.15 0.285 242.00 85.741 625.15 0.527 296.60 625.15 Ciclohexanol (2) 20
Esta aproximacin es buena en la medida de que los calores de vaporizacin de las sustancias a separar sean cercanos. Este mtodo usa el diagrama - para determinar el nmero de etapas en
21
la torre, el plato ptimo de alimentacin y las concentraciones por etapa. Si se apoya en el diagrama - - , tambin es posible conocer la temperatura por plato.
22
La ltima condicin es opcional, simplemente es necesario un dato adicional para caracterizar el alimento, lo cual se puede hacer fijando la temperatura, o fijando la fraccin de vapor; una vez fijado uno de estos valores, el otro queda determinado automaticamente; el procedimiento para este calculo se explica un poco ms adelante.
23
El cdigo elaborado esta diseado para una torre con una sola alimentacin, sin corrientes laterales, y considerando condensador total y evaporador, aunque puede ser fcilmente modificado, y adaptado al caso especfico de clculo. Para ir mostrando el desarrollo de los clculos se va a utilizar la informacin presentada en la seccin de descripcin del equipo y condiciones de operacin. En esta informacin, el reflujo de operacin tiene un valor dado y no esta en trminos del reflujo mnimo, por tanto no es necesario el clculo de este ltimo para determinar el nmero de etapas ideales, pero de igual forma se va a mostrar, para indicar el funcionamiento completo del cdigo
El cdigo viene separado por secciones para un manejo ms sencillo, y se va a presentar de igual manera.
Balance de materia
Como en el clculo de etapas de equilibrio, toda las variables involucradas son intensivas, la base de clculo seleccionada no ha de afectar en el nmero de etapas obtenidas. Acorde a esto se selecciona una base de clculo de F = 1, teniendo encuenta que es un flujo molar, despus se toman las composiciones de cada una de las corrientes, y utilizando los balances de materia global y por componente se determina la fraccin de destilado obtenido, y si se obtiene un resultado tal que (( )( )), se muestra un mensaje de Error- balance de masa inconsistente. Del balance de masa global se puede expresar como
24
El estimado inicial se realiza calculando las temperaturas de burbuja y de rocio a la presin de la torre, y ponderandolas de acuerdo a la fraccin de vapor especificada.
Para el calculo del siguiente valor, la derivada de la funcin objetivo se evalua por el mtodo de la secante ( ) ( ) ( ) ( ( Las iteraciones se mantienen hasta que termina el ciclo y se guarda el valor de . ( ) ) ) , cuando se cumple la condicin se
Un algoritmo BURB_T, prmite determinar la fraccin molar del vapor en equilibrio, y se obtiene
Con estos valores y , ya se conoce el intercepto de la curva de equilibrio con la lnea de alimentacin, este resultado se grafica, y se traza la lnea de rectificacin a reflujo mnimo. 25
Entonces
26
Como es conocidoe el valor real del reflujo de operacin, y el reflujo mnimo tambin, se determina
Se grafica la nueva lnea de operacin en la figura 10. Ahora , como , la lnea de alimento es una recta vertical, el valor de , donde se cortan la lnea de operacin real y la de alimento, es
En caso de que no ea igual a , se deben igualar las lneas de operacin y de alimentacin, despejar el valor de , y con esa valor hallar con la ecuacin (8.8)
*(
)(
)+
27
Figura 10. Lnea de rectificacin Con el valor de los interceptos se puede evaluar y graficar la lnea de agotamiento, pues se conocen dos puntos de esta recta
28
La composicin del siguiente plato se determina con la ecuacin de operacin, empezando con la lnea de rectificacin 29
Se calcula la composicin , y esto se realiza hasta que En el momento que esta condicin se deja de cumplir, hay cambio de seccin en la torre, se guarda el valor del plato de alimentacin, y se mantiene la misma secuencia, solo que ahora utilizando la ecuacin de la zona de agotamiento.
Los ciclos se mantienen, y se realizan hasta quela composicin del lquido en el plato sea menor que la composicin requerida de producto de fondo. En este momento se detiene el clculo, y como generalmente la composicin final es diferente de la composicin requerida, se determina la fraccin de esa etapa necesaria, es decir, la etapa final se no se cuenta como una etapa completa, si no como una fraccin
Figura 12. Resultado final del cdigo para McCabe-Thiele Finalmente, los resultados obtenidos para el sistema seleccionado por el mtodo de McCabeThiele son, 38.0758 Etapas ideales, plato de alimentacin en la etapa 32. 30
Tambin es posible observar los perfiles formados dentro de la torre, tanto los perfiles de composicin como los de temperatura.
31
Estos puntos muestran la lnea de equilibrio para el alimento, as que solo es generar la ecuacin de una recta con esos puntos, y hallar para que valor de entalpa corta el destilado ( ). En la figura 13 se muestra (morado), y de operacin (Negro) para un reflujo de 4.25
32
Ahora teniendo dos puntos de la recta de alimentacin, solamente falta determinar el pivote en la seccin de agotamiento, lo cual se realiza con los puntos ya conocidos ( ( ) )
Con los pivotes determinados, se obtiene la recta mostrada en la figura 13, en color negro, donde una recta pasa por los dos pivotes y por el punto de alimento.
En el mtodo grfico, las lneas de operacin toman el pivote como referencia y ascienden hasta la curva de equilibrio de la fase de vapor, en el diagrama entalpa composicin. Para lograr reproducir esto de manera analtica, cada vez que se va a calcular el nuevo valor de la fraccin de vapor en la etapa siguiente
Es necesario aplicar un clculo iterativo. Previamente se ha creado una funcin Calc_H_v que determina la entalpa de vapor saturado dada una composicin, el programa determina la temperatura por un BURB_T, y luego pasa al clculo de la entalpa con esa temperatura. Para el clculo de la etapa n, conozco la etapa n-1, as que con el pivote y el punto ( ), y se obtiene una ecuacin de recta con la forma
33
Esta ecuacin es la que se iguala con la funcin para el clculo de la entalpa, se resuelve de forma iterativa para , y finalmente se continua con los ciclos. Es necesario tambin determinar el valor de , para el cual se debe hacer el cambio de seccin. Este valor de x corresponde al corte de la recta de alimentacin con la curva de lquido saturado. Por lo cual se plantea un algoritmo similar al anterior, donde ahora la variable a encontrar no ser , sino . ( )
Tambin se dispone de una funcin que calcula la entalpa del lquido saturado Calc_h_l, esta funcin junto con la ecuacin (8.17), se resuelve de forma iterativa y permite determina el intercepto exacto con la curva de equilibrio. Una vez que la composicin del lquido es menor a este valor, se cambia de seccin en la torre, tomando como pivote el producto de fondo, en este caso la ecuacin (8.16) se cambia por ( )
Pero la estructura de solucin se mantiene. Los ciclos se realizan hasta que la composicin del lquido en el plato sea menor que la composicin requerida en el fondo ( ). Los resultados obtenidos con este algoritmo son los mostrados en la figura 16. Finalmente, los resultados obtenidos para el sistema seleccionado por el mtodo de Ponchon-Savarit son, 38.7145 Etapas ideales, plato de alimentacin en la etapa 33. En las figuras 17 y 18 se pueden observar los perfiles de temperatura y composicin respectivamente, en los que se observa una lnea suave y continua, lo que indica que el alimento est entrando por un lugar apropiado. Tambin se observa que el perfil de temperatura no cambia ms de 2 C en toda la torre, este pequeo rango sugiere que la separacin es muy sensible a la temperatura, y que un pequeo error de clculo en la temperatura puede modificar significativamente el diseo de la torre, llevando a resultados fuera de la realidad. En el perfil de composiciones se observa esta misma continuidad, y en ambos grficas se observa que los cambios ms significativos en temperatura y composicin vienen en la seccin de agotamiento, despus de introducir el alimento a la torre, lo que es de esperar, pues al alimento ser lquido, baja completamente por la seccin de agotamiento, generando flujos molares mayores, y sumado a esto, que la seccin de rectificacin es de muy alta pureza, y en este sector los cambios son mucho ms pequeos respecto a otras partes de la torre.
34
35
32
33
20
53,3
55,0
Lo primero que se observa es que los resultados entre cada uno de los mtodos no vara mucho, lo que significa, que para el sistema propuesto, trabajar con el mtodo de McCabe-Thiele no genera errores significativos en los clculos. Tambin se observa que la eficiencia global de la torre esta en valores cercanos al 50 %, se necesita aproximadamente el doble de las etapas ideales calculadas para lograr la separacin deseada.
36
Fraccin de vapor
Para la fraccin de vapor se toman 4 datos a igual distancia entre 0 y 1.
37
nRo = 1.1
nRo = 1.2
38
nRo = 1.5
nRo = 2.0
39
9.1.2 Lquido-Vapor
nRo = 1.1
nRo = 1.2
40
nRo = 1.5
nRo = 2.0
41
9.1.3 Lquido-Vapor
nRo = 1.1
nRo = 1.2
42
nRo = 1.5
nRo = 2.0
43
nRo = 1.1
nRo = 1.2
44
nRo = 1.5
nRo = 2.0
45
nRo = 1.1
nRo = 1.2
46
nRo = 1.5
nRo = 2.0
47
9.2.2 Lquido-Vapor
nRo = 1.1
nRo = 1.2
48
nRo = 1.5
nRo = 2.0
49
9.2.3 Lquido-Vapor
nRo = 1.1
nRo = 1.2
50
nRo = 1.5
nRo = 2.0
51
nRo = 1.1
nRo = 1.2
52
nRo = 1.5
nRo = 2.0
53
McCabe-Thiele
Tabla 9. Reflujo Mnimo en funcin de la fraccin de vapor Fraccin Reflujo de vapor Mnimo Ro 0 3,07577128 1/3 3,08491566 2/3 3,09408756 1 3,10349066 Los datos obtenidos muestran una dependencia completamente lineal del reflujo mnimo respecto a la fraccin de vapor.
McCabe-Thiele
3.105 3.1 Reflujo Mnimo 3.095 3.09 3.085 3.08 3.075 3.07 0 0.2 0.4 0.6 Fraccin de vapor 0.8 1 y = 0.0277x + 3.0757 R = 1
Tabla 10. Nmeo de etapas ideales mtodo de McCabe-Thiele F\nRo 0 1/3 2/3 1 1.1 72,4978308 74,3028128 76,3313115 78,4620909 1.2 51,5281426 52,4209917 53,3979158 54,4156373 1.5 33,8365627 34,1129385 34,380382 34,7346202 2.0 25,9813406 26,0578111 26,1567861 26,2553902
En los resultados obtenidos se puede observar que el aumento en la fraccin de vapor, trae consigo un aumento en el reflujo mnimo, lo que es desventajoso para el proceso, pues de cierta forma va a aumentar un poco los gastos de agua de enfriamiento y de vapor para el evaporador, pues a mayor reflujo, mayor cantidad de materia a evaporar y a condensar. En otro aspecto, cuando observamos el nmero de etapas ideales, es claro que al aumentar el reflujo, se disminuyen el nmero de etapas. En la figura 52 se puede observar que este efecto positivo de aumentar el reflujo se mantiene en el rango estudiado, pero se observa que a partir de 1.5 a 2, la pendiente de la curva se es menos pronunciada, es decir, que cada vez que aumentemos el reflujo, ser menor el beneficio recibido. Respecto al efecto de la fraccin de vapor en el nmero de etapas, el efecto no resulta ser muy significativo en el rango de estudio, pero se puede observar que un alimento como lquido saturado necesita menos etapas que un alimento como vapor saturado. Y a mayor reflujo, menor es la diferencia entre los platos obtenidos para diferentes alimentos.
54
McCABE-THIELE
80
50
30
Figura 52. Nmero de etapas ideales en funcin de la fraccin de vapor y el factor de escalamiento del reflujo
Ponchon-Savarit.
Reflujo Mnimo
Tabla 11. Reflujo Mnimo en funcin de la fraccin de vapor Fraccin Reflujo de vapor Mnimo Ro 0 3,07577128 1/3 3,08491566 2/3 3,09408756 1 3,10349066 Aunque de una forma un poco ms leve, los datos marcan una dependencia prcticamente lineal del reflujo mnimo respecto a la fraccin de vapor.
Ponchon-Savarit
3.2 3.19 3.18 3.17 3.16 3.15 3.14 3.13 0 0.2 0.4 0.6 Fraccin de vapor 0.8 1 y = 0.049x + 3.1391 R = 0.9985
55
Tabla 12. Nmero de etapas ideales Mtodode Ponchon Savarit F\nRo 1.1 1.2 1.5 2.0 0 69,5237004 50,2561063 33,3758539 25,7869438 1/3 70,8557952 50,9435103 33,5974403 25,8436349 2/3 72,2021257 51,5779222 33,8379327 25,9117651 1 73,5589772 52,2715347 34,0524878 26,0040705
Ponchon-Savarit
80
50
30
Figura 54. Nmero de platos ideales en funcin de la fraccin de vapor y el factor de escalamiento del reflujo.
56
Tabla 13. Desviaciones del reflujo mnimo mtodo aproximado respecto al riguroso Fraccin Reflujo Mnimo de vapor Ro 2,04% 0 2,22% 1/3 2,42% 2/3 2,68% 1 Tabla 14. Desviaciones del Nmero de etapas ideales del mtodo aproximado respecto al riguroso F\nRo 1.1 1.2 1.5 2.0 0 4,28% 2,53% 1,38% 0,75% 1/3 4,86% 2,90% 1,53% 0,83% 2/3 5,72% 3,53% 1,60% 0,95% 1 6,67% 4,10% 2,00% 0,97% Se observa que las desviaciones del mtodo de McCabe-Thiele respecto al de Ponchon-Savarit aumentan al aumentar la fraccin de vapor, y se hacen ms grandes a menor reflujo. Es de esperar que a menor reflujo el error sea mayor, pues a menor reflujo mayor nmero de etapas, y un mayor nmero de etapas implica un mayor nmero de clculos, y el efecto de propagacin del error hace que esta magnitud crezca significativamente al aumentar en esas proporciones el nmero de clculos a desarrollar. Es tambin se puede ver reflejado en el estado lquido, pues el nmero de etapas es siempre ligeramente menor para lquido que para vapor. Tambin se debe recordar que el mtodo aproximado se aproxima al riguroso en la medida que los calores de vaporizacin son iguales. Y cuando involucramos el clculo de entalpa de vapor, es cuando ms significativas se hacen las correcciones que el mtodo de Poncho-Savarit realiza al mtodo de McCabe-Thiele, es decir, probablemente a mayor cantidad de vapor, ms perjudiciales se vuelven las suposiciones del mtodo de McCabe-Thiele, y se vuelve mayor el error en sus clculos.
57
10 Nomenclatura
Fraccin de Destilado respecto al alimento. Fraccin de Fondo respecto al alimento. Composicin global del componente ms volatil. Presin de operacin de la torre. Fraccin Molar del componente ms volatil en el lquido. Fraccin Molar del componente ms volatil en el vapor. Temperatura, (K) El factor para escalar el reflujo real en trminos del Reflujo mnimo. Reflujo de operacin en la cabeza. Reflujo de operacin en el Evaporador. Fraccin de vapor de la corriente de alimento Coeficiente de actividad Entalpa devapor o global Entalpa de lquido
Subndices
Corriente de Alimento Corriente de Destilado Corriente de Fonfo Plato N en la seccin de rectificacin Plato M en la seccin de agotamiento Punto de interseccin de las lneas de operacin y la lnea de alimento. Condiciones mnimas de operacin. Plato n en cualquier seccin de la torre
Superndices
Pivote para el clculo en los diagramas de Entalpa composicin.
58
11 Bibliografa
1. Murthy et al, Cyclohexanone Distillation Process, U.S. patent N 4.306.944, Morris County New Jersey, The United State, 1981. 2. http://www.indiamart.com/gsfcltd/products.html Visitado en Octubre de 2011 3. KIRK-OTHMER, ENCYCLOPEDIA OF CHEMICAL TECHNOLOGY, vol 7, JONH WILLEY & SONS, The United Stated, 1982, 851p-858p. 4. Cova, D.R., J. Chem. Eng. Data, 5, 282 (1960). 5. Engelmann, K., Bittrich, H.J., J. Prakt. Chem., 19, 106 (1963). 6. Goodwin S.R., Newsham D.M.T., J. Chem. Eng. Data 19, 363 (1974). 7. Von Weber, U., J. Prakt. Chem., 312, 214 (1970). 8. Renon H., Prausnitz J.M., AIChe J., 14, 135 (1968). 9. Yaws C; Chemical Properties Handbook, Mc Graw-Hill, 1992. 10. Treybal R.E., Operaciones De Transferencia De Masa, 2 Edicin, McGraw-Hill.
59