DISEO Y SIMULACIN DE EVAPORADORES MULTIEFECTO RESUMEN Se desarroll un programa para el Diseo y Simulacin de Evaporadores Multiefecto utilizados en la manufactura de jugos

de fruta concentrados. Resuelve simultneamente el sistema de ecuaciones de transferencia de calor y balances de materia en cada efecto, previa linealizacin como se describe en Ugrin y Urbicain (1999). El programa puede operar en dos modos: Diseo: Para una alimentacin de caudal y concentracin conocidos, calcula el rea de intercambio requerida para lograr la concentracin final especificada del producto Simulacin: Dado un equipo de rea de intercambio conocida, calcula la concentracin de salida del producto, para una alimentacin de caudal y concentracin conocidos.

Precisamente la caracterstica destacada del mismo es la posibilidad de efectuar ambos tipos de clculo, mediante el artificio de convertir un problema de Simulacin en uno de Diseo. Se contemplan algunas opciones que el usuario elige a voluntad, tales como la existencia de ascenso ebulloscpico de la solucin en cada efecto, la posibilidad de precalentamiento de la alimentacin en haces de tubos independientes pero alojados en los efectos y el manejo de reas y coeficientes de transferencia de calor distintos para cada efecto. Opera indistintamente con arreglo de flujo en cocorriente o en contracorriente, esto es con la corriente de proceso en el mismo, u opuesto sentido respectivamente, que el de los vapores de calefaccin. El programa incluye subrutinas para el clculo de las propiedades termofsicas del jugo correspondiente. BASE TERICA Durante dcadas, el Diseo y Simulacin de evaporadores se ha llevado a cabo en forma aproximada por el mtodo de prueba y error. La aparicin del diseo asistido por ordenador permiti el desarrollo de tcnicas numricas para allanar las dificultades que plantea la resolucin de sistemas de ecuaciones no lineales, como las que gobiernan el fenmeno, as como los requerimientos iterativos impuestos por caractersticas fsicas del sistema, tales como el ascenso ebulloscpico, . Formulacin del Problema El fundamento terico del trabajo puede encontrarse, detallado, en el trabajo de Ugrin y Urbicain (1999), cuyos aspectos ms relevantes se resumen en los siguientes prrafos.

Pag. 2 de 12

La siguiente es la formulacin correspondiente a un sistema en contracorriente, pero el razonamiento es completamente vlido para un arreglo co-corriente, por un simple cambio de subndices de las variables involucradas. Con referencia a la Figura 1, las ecuaciones de gobierno se pueden escribir como un sistema de 3N ecuaciones con 3N incgnitas.

V0 , T0

V1 , T1

Vi-1 , Ti-1

Vi , Ti

VN-1 , TN-1

VN , TN

L1 X1 T1

Efecto 1 L2 X2 T2 V0 (cond.)

Li Xi Ti

Efecto i

Li+1 Xi+1 Ti+1

LN XN TN

Efecto N LA XA TA VN-1 (cond.)

Vi-1 (cond.)

Figura 1. Esquema de un evaporador de N efectos en contracorriente ) N balances de masa Primer efecto: Efecto isimo: N-1) (2) L 2 - L1 - V 1 = 0 Li+1 - Li - Vi = 0 ltimo efecto (3) (1) ( i = 2, , LF- LN - VN= 0

) N balances entlpicos: Primer efecto: Efecto isimo: (5) Ultimo efecto: (6) H0 V0 + h2L2- h1L1- H1V1= 0 (4) Hi-1 Vi-1 + hi+1 Li+1 - hi Li - Hi Vi = 0 (i = 2, , NHN-1 VN-1 + hF LF - hN LN - HN VN = 0

1)

) N ecuaciones de transferencia de calor: Primer efecto: (7) H0 V0 - U1 R1 A (T0 - T1) = 0

Pag. 3 de 12

Efecto i-simo: (8) Ultimo efecto: (9)

Hi-1 Vi-1 - Ui Ri A (Tsi-1 - Ti) = 0

( i = 2, , N-1)

HN-1 VN-1 - UNRN A (TsN-1 - TN) = 0

El problema de Diseo se resuelve sin dificultades, porque en las diferencias de temperaturas, si se introduce el ascenso ebulloscpico, ste puede calcularse iterando sobre la variable concentracin de modo que, finalmente, se tiene una constante conocida, que afecta la fuerza impulsora definida por las temperaturas del vapor y el lquido saturado. El de Simulacin, en cambio, introduce una variable adicional, puesto que la concentracin de salida es desconocida. Esto hace que el modo Simulacin no sea tratado con frecuencia. En este trabajo, con un adecuado manejo para eliminar las no linealidades, el sistema se transforma en uno de diseo mediante la adopcin de aquella incgnita. Inicializacin Dado que xp est acotada entre x F y 1, se adopta como valor inicial x p = xF + 0.1. La temperatura inicial en cada efecto se calcula distribuyendo el salto total de temperaturas, Tt, en forma inversamente proporcional a las resistencias a la transferencia en cada uno. Si todos los coeficientes son iguales, resulta de dividir Tt por el nmero de efectos. Si se contempla la existencia de aumento ebulloscpico, se calculan aproximadamente las concentraciones en cada efecto. Para ello se admite que todos producen la misma cantidad de vapor, resultado de dividir el vapor total producido, Vt, por el nmero de efectos. A su vez V t surge del balance global de materia una vez adoptada xp. La solucin numrica del sistema se llev a cabo mediante la Subrutina LEQTIF de la biblioteca IMSL, que es una versin modificada del mtodo de Gauss (mtodo de Crout), con pivoteo parcial por fila. Prob ser particularmente til para este tipo de problemas. Convergencia En el modo Diseo, el criterio de convergencia est basado en el caudal del vapor vivo requerido, V0, de tal modo que dos valores sucesivos de la variable sean aceptablemente prximos. En smbolos:

V0

k 1 k

V0
(10)

1 e

donde V0k-1 y V0k son los V0 obtenidos en las iteraciones (k - 1) y respectivamente.

Pag. 4 de 12

En cada iteracin, con el valor calculado de V 0, Ts y B1 para el flujo contracorriente (o BN, en co-corriente), el area requerida A r se calcula con la equacin (11)

Ar =

0 V0 U 1 R 1 ( T0 ( Ts1 + B1 ) )

(11)

Si Ar es menor que A, x p debe ser incrementada y viceversa. El nuevo valor es el inicial de la siguiente iteracin. Eventualmente, A r resultar mayor que A, con lo que la solucin del sistema se dar para un valor de x p ubicado entre los dos ltimos ciclos, tal que:

x x x
(12)

k 1 k p p p

En este punto, habiendo acotado el valor de la concentracin del producto, se cambia el criterio de convergencia por otro que garantice una ms rpida solucin. A partir de experimentos numricos se ha podido establecer una relacin entre la composicin y el area del tipo:
1 b

x = a.A
(13)

donde a y b son dos constantes. Por lo tanto los valores sucesivos de x p se calculan con la ecuacin (14)
x
k +1

A =x Ar
k

1/ b

(14) En el modo Simulacin, que se transforma en uno de Diseo, la solucin del problema se obtiene cuando el rea calculada es igual a la disponible. En rigor esa igualdad se da con un cierto error admitido, tal como se expresa en la ecuacin (15).

Ar 1 e A donde e es el error especificado.

(15)

Pag. 5 de 12

En el Apndice se describe, como ejemplo, el algoritmo utilizado para el modo Simulacin RESULTADOS El programa se aplic a la Simulacin de un Evaporador de Jugo de Manzana, tal como se muestra en el ejemplo presentado en el Manual del Usuario. La consistencia interna del mtodo qued probada puesto que los resultados de la simulacin se alimentaron como datos del diseo, recuperndose como resultados los datos del primer problema. Para su ejecucin, se le incorporaron subrutinas para calcular propiedades termofsicas del jugo, implementadas con correlaciones seleccionadas de la literatura. En particular calor especfico (Constenla y asoc., 1989), y ascenso ebulloscpico (Crapiste y Lozano, 1988). Las figuras 2, 3, 4 y 5, son impresiones de las pantallas con los datos y resultados. CONCLUSIONES El mtodo propuesto ha sido probado con numerosos ejemplos de literatura y de instalaciones industriales, con resultados muy ssatisfactorios, tanto el modo diseo como en simulacin. Opera en forma estable y converge invariablemente La estructura del sistema permite la incorporacin de precalentamiento de la alimentacin en haces de tubos incorporados a los efectos, de modo que se calefaccionen con el mismo vapor utilizado para la ebullicin de la solucin. Simplemente requiere adicionar las ecuaciones de transferencia y de balance entlpico que representan el caslentamiento. Tambin admite la consideracin de termocompresin del vapor efluente de algn efecto para ser realimentado al mismo o a otro efecto precedente, si bien esa alternativa no fu includa en el presente trabajo.

Pag. 6 de 12

APENDICE ALGORITMO PARA EL CASO SIMULACION Datos necesarios: Area de referencia A, y razn de las reas de cada efecto respecto de ella, Ri Nmero de efectos N Temperatura o presin de saturacin del efecto N Temperatura o presin del vapor vivo, si es saturado, o ambos en caso contrario Caudal masa, composicin y temperatura de la corriente de alimentacin.

Resultados Composicin del producto Cuadales masa de las corrientes de lquido y de vapor en todos los efectos Consumo de vapor vivo

Algoritmo 1. Adopte una composicin del producto xp, entre xF y 1. 2. Inicialice temperaturas y composiciones en todos los efectos. 3. Calcule, mediante subrutinas adecuadas, los factores involucrados en los coeficientes de las equaciones (H,U, B) 4. Resuelva el sistema de ecuaciones de (1) a (9) obteniendo todos los parmetros operativos. 5. Calcule el area requerida AT, correspondiente al valor de x p adoptado, por medio de la ecuacin (14). 6. Compare Ar con el rea calculada A. Si A r resultara menor que A, incremente x y vuelva al paso 2. De lo contrario contine con el paso 7. 7. Seleccione un valor de indicated in Eq. (12). 8. Ejecute el paso 2. 9. Ejecute los pasos 3 y 4. 10.Calcule el rea requerida A r, correspondiente al valor interpolado x p, con la Eq. (11). 11.Aplique el criterio de convergencia expresado por la desigualdad (15). Si se satisface vaya al paso 13, con lo que finaliza el clculo. 12.Actualice la composicin xi en cada efecto or medio de la Eq. (17). Vaya al paso 9. 13. FIN DEL CLCULO
k k 1 x pdefinido por los dos ltimos ciclos, x p y x p as

Pag. 7 de 12

NOMENCLATURA A AJC a b H h L N R T U V x B H Subndices F i N o t p r s 1 Superndices k * REFERENCIAS Constenla D.T. , Crapiste, G.H. y Lozano, J.E.. "Thermophysical properties of clarified apple juice as a function of concentration and temperature". J. Food Sci.(54), 663-669 (1989). Crapiste G. and Lozano, J.E.. "Effect of Concentration and Pressure on the Boiling Point Rise of Apple Juice and Related Sugar Solutions." J. Food Science (53), 3, 865 - 868, l988. Ugrin, E. y Urbicain, M. Design and Simulation of Multieffect Evaporators , Heat Transfer Engineering, Vo. 20, No. 4, 1999 (38-44). k-sima iteracin Correspondiente a la solucin numrica Alimentacin Efecto I Efecto N Vapor vivo Total Producto Requerido Condiciones de saturacin Efecto 1 Area de intercambio calrico de referencia [m2] Jugo Concentrado de Manzana Par;ametro en la Eq. (16) [1/m2/b] Parmetro en la Eq. (16) [-] Entalpa del vapor saturado [J/kg] Entalpa del lquido a temperatura de ebullicin [J/kg] Caudal masa del lquido [kg/h] Numero de efectos del evaporador Razn del rea calculada y la de referencia Temperatura [C] Coeficiente global de transferenia de calor [W / m2.C] Caudal masa de vapor [kg/h] Fraccin en peso de slidos en la solucin [-] Ascenso ebulloscpico [C] Calor latente de vaporizacin [J/kg]

AGRADECIMIENTOS Los autores agradecen al CONICET y a la Universidad Nacional del Sur el apoyo brindado al presente trabajo.

Pag. 8 de 12

MANUAL DEL USUARIO El programa MULTIEF para evaporadores de mltiple efecto lleva a cabo el clculo de las condiciones de salida para un rea conocida (modo simulacin), o el rea si se conocen las condiciones de salida (modo diseo). El programa de clculo fue desarrollado en Visual Basic 6.0, y est basado en el algoritmo de Ugrin y Urbicain (1999), escrito originalmente en lenguaje FORTRAN. Las libreras necesarias se suministran junto con el programa ejecutable MULTIEF.EXE, y debern ser guardadas en la misma carpeta que ste. Una vez descargados los archivos de los dos diskettes, se puede ejecutar la aplicacin. Al hacerlo, se carga un formulario que consta de 4 solapas: Datos Generales Datos de Evaporacin Datos de Precalentamiento Resultados Archivo Ejecutar

y una barra de mens con dos opciones:

A continuacin se describen las variables de cada una de las solapas, y las opciones de la barra de mens. Las palabras o cantidades en negrita itlica corresponden a valores que se eligen de un men desplegable. Cabe agregar que cuando una opcin elegida hace innecesario el ingreso de algn dato, automticamente la casilla correspondiente al mismo es deshabilitada, por lo que el usuario no podr ingresar dato alguno en ella. 1. Datos Generales Tipo de Clculo: Diseo (se conocen condiciones de las corrientes a la entrada y a la salida, y se desea calcular el rea), o Simulacin (se conocen las condiciones de entrada y el rea, y se deben calcular las de salida) Nmero de Efectos: el programa soporta de 2 a 5 efectos. Los efectos se numeran siguiendo la corriente de solucin a concentrar. Precalentamiento: las opciones son S o No. Ascenso Ebulloscpico: las opciones son S o No, y dependiendo de la opcin elegida el programa calcular, o no, el ascenso ebulloscpico debido al aumento en la concentracin de la solucin. Arreglo Corrientes: Co-corriente si el vapor vivo ingresa en el primer efecto, Contracorriente si lo hace en el ltimo. Alimentacin: Caudal, Concentracin, Temperatura: el caudal es variable de entrada para Diseo, y se calcula en modo Simulacin; las otras dos variables son entrada en todos los casos. Temperatura del Vapor: Vapor vivo: se ingresar la temperatura del vapor vivo a la presin de trabajo de la cmara de condensacin del primer efecto para disposicin en co-corriente, y del ltimo efecto para flujo en contracorriente.

Pag. 9 de 12

Temperatura de saturacin: es la correspondiente al equilibrio del vapor de agua a la presin de trabajo de la cmara de evaporacin del ltimo efecto. Concentracin del Producto: se debe ingresar para el modo Diseo. 2. Datos de Evaporacin Area de Evaporacin: es el valor de un rea de referencia, y se ingresa slo en el modo Simulacin. Relacin de Areas: las opciones son Iguales o Diferentes. En el primer caso no debe agregarse ningn otro dato; en el segundo deber llenarse la grilla que se presenta a continuacin con valores entre 0 y 1 para los diferentes efectos; estas fracciones constituirn el factor por el que se multiplicar el valor del rea de referencia (dato para Simulacin, calculada para Diseo) a fin de obtener el rea de cada efecto. Coeficiente Global de Transferencia de Calor: las opciones son Iguales o Diferentes. En el primer caso debe suministrarse el primer valor de U de la grilla, y en el ltimo tantos valores como efectos tenga el evaporador. 3. Datos de Precalentamiento Area de Precalentamiento: es el valor de un rea de referencia, y se ingresa slo en el modo Simulacin. Temperatura de salida del precalentador: est valor se ingresa slo en el caso de Simulacin con Precalentamiento. Relacin de Areas: las opciones son Iguales o Diferentes. En el primer caso no debe agregarse ningn otro dato; en el segundo deber llenarse la grilla que se presenta a continuacin con valores entre 0 y 1 para los diferentes efectos; estas fracciones constituirn el factor por el que se multiplicar el valor del rea de referencia (dato para Simulacin, calculada para Diseo) a fin de obtener el rea de cada efecto. Coeficiente Global de Transferencia de Calor: las opciones son Iguales o Diferentes. En el primer caso debe suministrarse el primer valor de U de la grilla, y en el ltimo tantos valores como efectos tenga el evaporador. 4. Resultados Caudal: Alimentacin, Producto y Vapor Vivo: en modo Simulacin todos son calculados por el programa; en modo Diseo el caudal de Alimentacin es dato. Composicin: Alimentacin y Producto: ambos son datos para Diseo, y para Simulacin la composicin de la alimentacin es dato, y la del producto se calcula. Temperatura: Alimentacin, Producto y Vapor Vivo: la primera y la ltima son datos, y la del producto la calcula el programa en funcin del valor ingresado para el equilibrio a la presin de trabajo de la ltima cmara de evaporacin y el ascenso ebulloscpico de la mezcla lograda en la misma. Area de Evaporacin y Precalentamiento: ambas son la suma del rea correspondiente a los N efectos. Son calculadas tanto en modo Diseo como

Pag. 10 de 12

Simulacin; en el primer caso se desconoce su valor, y en el segundo el usuario ingresa un rea de referencia y relaciones del rea de cada efecto con sta. Economa de Vapor: cociente entre cantidad de vapor generada en todos los efectos y vapor vivo ingresado al sistema. Resultados en cada efecto: se presentan en forma de tabla los valores que toman las variables a la salida de cada efecto, tanto en evaporacin como en precalentamiento, si existiera. 5. Men Archivo Nuevo: carga un nuevo formulario en blanco. Abrir: permite cargar un archivo de datos existente. Estos son archivos de texto, y tienen la extensin .mtf. Guardar Como: permite guardar en un archivo de texto los datos del caso actualmente en uso. Los guarda con la extensin .mtf. Imprimir: imprime la solapa activa, o sea que si se desea imprimir los resultados deber estar visible la solapa Resultados, etc. Salir: con este comando se termina la ejecucin, de la misma manera que pulsando la cruz del extremo superior derecho. 6. Men Ejecutar Este men no es desplegable, y al hacer Click con el ratn sobre el mismo se ejecuta la rutina de diseo y simulacin de evaporadores. Al terminar la ejecucin, aparece un cartel que dice Programa ejecutado, y los resultados pueden verse en la solapa Resultados, que hasta ese momento haba aparecido con sus casilleros en blanco.

Ejemplo de Aplicacin Como ejemplo, se presenta el caso de Simulacin de un evaporador de 4 efectos sin precalentamiento, que debe procesar jugo con una concentracin de 0,105 kg de slidos solubles por kg de solucin. Se dispone de vapor a 123C, y se desea llegar a una concentracin de 0,70 kg de slidos solubles por kg de solucin a una presin tal que la temperatura de saturacin del vapor puro es 50C. La alimentacin ingresa a 45C, el funcionamiento es co-corriente y se calcula el ascenso ebulloscpico a medida que la solucin se concentra. Las reas de los efectos son iguales, con un valor de 48 m 2, y se estima un coeficiente de transferencia de calor de 2300 en el primer efecto, 2100 en el segundo, 1800 en el tercero y 1500 en el cuarto. En las figuras 2, 3, 4 se presentan impresiones de las pantallas sobre las que el usuario vuelca los datos, para luego pulsar Ejecutar. La pantalla correspondiente a los datos de precalentamiento aparece vaca porque no se ha contemplado en el ejemplo esa opcin. En la figura 5 se presentan los resultados. Obsrvese que en todos los casos, la solapa activa es la que se destaca con letra en negrita.

Pag. 11 de 12

Figura 2 Ingreso de datos generales

Figura 3 Ingreso de Datos del rea de evaporacin

Pag. 12 de 12

Figura 4 Ingreso de Datos del Area de Precalentamiento (no includa en este ejemplo)

Figura 5 Resultados de la Simulacin