DISEÑO Y SIMULACIÓN DE EVAPORADORES MULTIEFECTO

M.T. González, E. Ugrin y M.Urbicain

PLAPIQUI (Univ. Nacional del Sur-CONICET)
(8000) Bahía Blanca – ARGENTINA
Tel +54 (291) 486 1700 – e-mail:urbicain@criba.edu.ar

RESUMEN
Se desarrolló un programa para el Diseño y Simulación de Evaporadores Multiefecto
utilizados en la manufactura de jugos de fruta concentrados. Resuelve simultáneamente el
sistema de ecuaciones de transferencia de calor y balances de materia en cada efecto, previa
linealización como se describe en Ugrin y Urbicain (1999).
El programa puede operar en dos modos:
 Diseño: Para una alimentación de caudal y concentración conocidos, calcula el área de
intercambio requerida para lograr la concentración final especificada del producto
 Simulación: Dado un equipo de área de intercambio conocida, calcula la concentración
de salida del producto, para una alimentación de caudal y concentración conocidos.
Precisamente la característica destacada del mismo es la posibilidad de efectuar ambos tipos
de cálculo, mediante el artificio de convertir un problema de Simulación en uno de Diseño.
Se contemplan algunas opciones que el usuario elige a voluntad, tales como la existencia de
ascenso ebulloscópico de la solución en cada efecto, la posibilidad de precalentamiento de la
alimentación 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 calefacción.
El programa incluye subrutinas para el cálculo de las propiedades termofísicas del jugo
correspondiente.

BASE TEÓRICA
Durante décadas, el Diseño y Simulación de evaporadores se ha llevado a cabo en forma
aproximada por el método de prueba y error. La aparición del diseño asistido por ordenador
permitió el desarrollo de técnicas numéricas para allanar las dificultades que plantea la
resolución de sistemas de ecuaciones no lineales, como las que gobiernan el fenómeno, así
como los requerimientos iterativos impuestos por características físicas del sistema, tales
como el ascenso ebulloscópico, .

Formulación del Problema

El fundamento teórico del trabajo puede encontrarse, detallado, en el trabajo de
Ugrin y Urbicain (1999), cuyos aspectos más relevantes se resumen en los siguientes
párrafos.
La siguiente es la formulación correspondiente a un sistema en contracorriente, pero
el razonamiento es completamente válido para un arreglo co-corriente, por un simple cambio
de subíndices de las variables involucradas.
Diseño y simulación de Evaporadores Multiefecto – González et al. Pag. 2 de 12

Con referencia a la Figura 1, las ecuaciones de gobierno se pueden escribir como un

sistema de 3N ecuaciones con 3N incógnitas.

V0 , T0 V1 , T1 Vi-1 , Ti-1 Vi , Ti VN-1 , TN-1 VN , TN

Efecto Efecto Efecto

L1 1 L2 Li i Li+1 LN N LA
X1 X2 Xi Xi+1 XN XA
T1 T2 Ti Ti+1 TN TA

V0 (cond.) Vi-1 (cond.) VN-1 (cond.)

Figura 1. – Esquema de un evaporador de N efectos en contracorriente

) N balances de masa
Primer efecto: L2 - L1 - V1 = 0 (1)
Efecto iésimo: Li+1 - Li - Vi = 0 ( i = 2, …, N-1) (2)
Último efecto LF- LN - VN= 0 (3)

) N balances entálpicos:
Primer efecto: H0 V0 + h2L2- h1L1- H1V1= 0 (4)
Efecto iésimo: Hi-1 Vi-1 + hi+1 Li+1 - hi Li - Hi Vi = 0 (i = 2, …, N-1) (5)
Ultimo efecto: HN-1 VN-1 + hF LF - hN LN - HN VN = 0 (6)

) N ecuaciones de transferencia de calor:

Primer efecto: H0 V0 - U1 R1 A (T0 - T1) = 0 (7)

Efecto i-ésimo: Hi-1 Vi-1 - Ui Ri A (Tsi-1 - Ti) = 0 ( i = 2, …, N-1) (8)
Ultimo efecto: HN-1 VN-1 - UNRN A (TsN-1 - TN) = 0 (9)

El problema de Diseño se resuelve sin dificultades, porque en las diferencias de

temperaturas, si se introduce el ascenso ebulloscópico, éste puede calcularse iterando sobre
la variable concentración de modo que, finalmente, se tiene una constante conocida, que
afecta la fuerza impulsora definida por las temperaturas del vapor y el líquido saturado.
Diseño y simulación de Evaporadores Multiefecto – González et al. Pag. 3 de 12

El de Simulación, en cambio, introduce una variable adicional, puesto que la

concentración de salida es desconocida. Esto hace que el modo Simulación no sea tratado con
frecuencia. En este trabajo, con un adecuado manejo para eliminar las no linealidades, el
sistema se transforma en uno de diseño mediante la adopción de aquella incógnita.

Inicialización
Dado que xp está acotada entre xF 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 número de efectos.
Si se contempla la existencia de aumento ebulloscópico, 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, V t, por el número de efectos. A su vez
Vt surge del balance global de materia una vez adoptada xp.
La solución numérica del sistema se llevó a cabo mediante la Subrutina LEQTIF de la
biblioteca IMSL, que es una versión modificada del método de Gauss (método de Crout), con
pivoteo parcial por fila. Probó ser particularmente útil para este tipo de problemas.

Convergencia
En el modo Diseño, 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
próximos. En símbolos:

k 1
V0
k
1  e
V0
(10)
donde V0k-1 y V0k son los V0 obtenidos en las iteraciones (k - 1) y k respectivamente.
En cada iteración, con el valor calculado de V0, Ts y B1 para el flujo contracorriente (o
BN, en co-corriente), el area requerida Ar se calcula con la equación (11)

 0 V0
Ar  (11)
U 1 R 1  T0   Ts1  B1  

Si Ar es menor que A, xp debe ser incrementada y viceversa. El nuevo valor es el inicial

de la siguiente iteración. Eventualmente, Ar resultará mayor que A, con lo que la solución del
sistema se dará para un valor de xp ubicado entre los dos últimos ciclos, tal que:
Diseño y simulación de Evaporadores Multiefecto – González et al. Pag. 4 de 12

k1  k
x x x
p p p
(12)

En este punto, habiendo acotado el valor de la concentración del producto, se cambia

el criterio de convergencia por otro que garantice una más rápida solución. A partir de
experimentos numéricos se ha podido establecer una relación entre la composición y el area
del tipo:

x  a.A b (13)

donde a y b son dos constantes. Por lo tanto los valores sucesivos de x p se calculan con la
ecuación (14)
1/ b
k 1  A
x k
x   (14)
 Ar 

En el modo Simulación, que se transforma en uno de Diseño, la solución 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 ecuación (15).

Ar
 1  e (15)
A
donde e es el error especificado.
En el Apéndice se describe, como ejemplo, el algoritmo utilizado para el modo Simulación

RESULTADOS
El programa se aplicó a la Simulación de un Evaporador de Jugo de Manzana, tal
como se muestra en el ejemplo presentado en el Manual del Usuario. La consistencia interna
del método quedó probada puesto que los resultados de la simulación se alimentaron como
datos del diseño, recuperándose como resultados los datos del primer problema.
Para su ejecución, se le incorporaron subrutinas para calcular propiedades
termofísicas del jugo, implementadas con correlaciones seleccionadas de la literatura. En
particular calor específico (Constenla y asoc., 1989), y ascenso ebulloscópico (Crapiste y
Lozano, 1988).
Las figuras 2, 3, 4 y 5, son impresiones de las pantallas con los datos y resultados.

CONCLUSIONES
Diseño y simulación de Evaporadores Multiefecto – González et al. Pag. 5 de 12

El método propuesto ha sido probado con numerosos ejemplos de literatura y de

instalaciones industriales, con resultados muy ssatisfactorios, tanto el modo diseño como en
simulación. Opera en forma estable y converge invariablemente
La estructura del sistema permite la incorporación de precalentamiento de la
alimentación en haces de tubos incorporados a los efectos, de modo que se calefaccionen con
el mismo vapor utilizado para la ebullición de la solución. Simplemente requiere adicionar
las ecuaciones de transferencia y de balance entálpico que representan el caslentamiento.
También admite la consideración de termocompresión del vapor efluente de algún
efecto para ser realimentado al mismo o a otro efecto precedente, si bien esa alternativa no
fué incluída en el presente trabajo.
Diseño y simulación de Evaporadores Multiefecto – González et al. Pag. 6 de 12

APENDICE
ALGORITMO PARA EL CASO SIMULACION

Datos necesarios:
 Area de referencia A, y razón de las áreas de cada efecto respecto de ella, R i
 Número de efectos N
 Temperatura o presión de saturación del efecto N
 Temperatura o presión del vapor vivo, si es saturado, o ambos en caso contrario
 Caudal masa, composición y temperatura de la corriente de alimentación.

Resultados
 Composición del producto
 Cuadales masa de las corrientes de líquido y de vapor en todos los efectos
 Consumo de vapor vivo

Algoritmo
1. Adopte una composición 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 parámetros
operativos.
5. Calcule el area requerida AT, correspondiente al valor de xp adoptado, por medio de la
ecuación (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 continúe con el paso 7.
k k 1
7. Seleccione un valor de xpdefinido por los dos últimos ciclos, x p y x p as indicated in
Eq. (12).
8. Ejecute el paso 2.
9. Ejecute los pasos 3 y 4.
10. Calcule el área requerida Ar, correspondiente al valor interpolado xp, 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 cálculo.
12. Actualice la composición xi en cada efecto or medio de la Eq. (17). Vaya al paso 9.
13. FIN DEL CALCULO
Diseño y simulación de Evaporadores Multiefecto – González et al. Pag. 7 de 12

NOMENCLATURA

A Area de intercambio calórico de referencia [m2]

AJC Jugo Concentrado de Manzana
a Par;ametro en la Eq. (16) [1/m2/b]
b Parámetro en la Eq. (16) [-]
H Entalpía del vapor saturado [J/kg]
h Entalpía del líquido a temperatura de ebullición [J/kg]
L Caudal masa del líquido [kg/h]
N Numero de efectos del evaporador
R Razón del área calculada y la de referencia
T Temperatura [C]
U Coeficiente global de transferenia de calor [W / m2.C]
V Caudal masa de vapor [kg/h]
x Fracción en peso de sólidos en la solución [-]
B Ascenso ebulloscópico [C]
H Calor latente de vaporización [J/kg]

Subíndices

F Alimentación
i Efecto I
N Efecto N
o Vapor vivo
t Total
p Producto
r Requerido
s Condiciones de saturación
1 Efecto 1

Superíndices

k k-ésima iteración
* Correspondiente a la solución numérica

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).

AGRADECIMIENTOS
Los autores agradecen al CONICET y a la Universidad Nacional del Sur el apoyo brindado
al presente trabajo.
Diseño y simulación de Evaporadores Multiefecto – González et al. Pag. 8 de 12

MANUAL DEL USUARIO

El programa MULTIEF para evaporadores de múltiple efecto lleva a cabo el cálculo

de las condiciones de salida para un área conocida (modo simulación), o el área si se
conocen las condiciones de salida (modo diseño). El programa de cálculo fue desarrollado en
Visual Basic 6.0, y está basado en el algoritmo de Ugrin y Urbicain (1999), escrito
originalmente en lenguaje FORTRAN.
Las librerías necesarias se suministran junto con el programa ejecutable
MULTIEF.EXE, y deberán ser guardadas en la misma carpeta que éste. Una vez descargados
los archivos de los dos diskettes, se puede ejecutar la aplicación. Al hacerlo, se carga un
formulario que consta de 4 solapas:
 Datos Generales
 Datos de Evaporación
 Datos de Precalentamiento
 Resultados
y una barra de menús con dos opciones:
 Archivo
 Ejecutar
A continuación se describen las variables de cada una de las solapas, y las opciones
de la barra de menús. Las palabras o cantidades en negrita itálica corresponden a valores que
se eligen de un menú desplegable.
Cabe agregar que cuando una opción elegida hace innecesario el ingreso de algún
dato, automáticamente la casilla correspondiente al mismo es deshabilitada, por lo que el
usuario no podrá ingresar dato alguno en ella.

1. Datos Generales
Tipo de Cálculo: Diseño (se conocen condiciones de las corrientes a la entrada y a la
salida, y se desea calcular el área), o Simulación (se conocen las condiciones de entrada y el
área, y se deben calcular las de salida)
Número de Efectos: el programa soporta de 2 a 5 efectos. Los efectos se numeran
siguiendo la corriente de solución a concentrar.
Precalentamiento: las opciones son Sí o No.
Ascenso Ebulloscópico: las opciones son Sí o No, y dependiendo de la opción elegida
el programa calculará, o no, el ascenso ebulloscópico debido al aumento en la concentración
de la solución.
Arreglo Corrientes: Co-corriente si el vapor vivo ingresa en el primer efecto,
Contracorriente si lo hace en el último.
Alimentación: Caudal, Concentración, Temperatura: el caudal es variable de entrada
para Diseño, y se calcula en modo Simulación; las otras dos variables son entrada en todos
los casos.
Temperatura del Vapor: Vapor vivo: se ingresará la temperatura del vapor vivo a la
presión de trabajo de la cámara de condensación del primer efecto para disposición en co-
corriente, y del último efecto para flujo en contracorriente. Temperatura de saturación: es la
Diseño y simulación de Evaporadores Multiefecto – González et al. Pag. 9 de 12

correspondiente al equilibrio del vapor de agua a la presión de trabajo de la cámara de

evaporación del último efecto.
Concentración del Producto: se debe ingresar para el modo Diseño.

2. Datos de Evaporación
Area de Evaporación: es el valor de un área de referencia, y se ingresa sólo en el
modo Simulación.
Relación de Areas: las opciones son Iguales o Diferentes. En el primer caso no debe
agregarse ningún otro dato; en el segundo deberá llenarse la grilla que se presenta a
continuación con valores entre 0 y 1 para los diferentes efectos; estas fracciones constituirán
el factor por el que se multiplicará el valor del área de referencia (dato para Simulación,
calculada para Diseño) 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 sólo en el
modo Simulación.
Temperatura de salida del precalentador: est valor se ingresa sólo en el caso de
Simulación con Precalentamiento.
Relación de Areas: las opciones son Iguales o Diferentes. En el primer caso no debe
agregarse ningún otro dato; en el segundo deberá llenarse la grilla que se presenta a
continuación con valores entre 0 y 1 para los diferentes efectos; estas fracciones constituirán
el factor por el que se multiplicará el valor del área de referencia (dato para Simulación,
calculada para Diseño) 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: Alimentación, Producto y Vapor Vivo: en modo Simulación todos son
calculados por el programa; en modo Diseño el caudal de Alimentación es dato.
Composición: Alimentación y Producto: ambos son datos para Diseño, y para
Simulación la composición de la alimentación es dato, y la del producto se calcula.
Temperatura: Alimentación, Producto y Vapor Vivo: la primera y la última son
datos, y la del producto la calcula el programa en función del valor ingresado para el
equilibrio a la presión de trabajo de la última cámara de evaporación y el ascenso
ebulloscópico de la mezcla lograda en la misma.
Area de Evaporación y Precalentamiento: ambas son la suma del área
correspondiente a los N efectos. Son calculadas tanto en modo Diseño como Simulación; 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.
Diseño y simulación de Evaporadores Multiefecto – González et al. Pag. 10 de 12

Economía 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 evaporación 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 extensión “.mtf”.
Guardar Como: permite guardar en un archivo de texto los datos del caso
actualmente en uso. Los guarda con la extensión “.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 ejecución, 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 ratón sobre el mismo se ejecuta la
rutina de diseño y simulación de evaporadores. Al terminar la ejecución, aparece un cartel
que dice “Programa ejecutado”, y los resultados pueden verse en la solapa “Resultados”, que
hasta ese momento había aparecido con sus casilleros en blanco.

Ejemplo de Aplicación

Como ejemplo, se presenta el caso de Simulación de un evaporador de 4 efectos sin

precalentamiento, que debe procesar jugo con una concentración de 0,105 kg de sólidos
solubles por kg de solución. Se dispone de vapor a 123°C, y se desea llegar a una
concentración de 0,70 kg de sólidos solubles por kg de solución a una presión tal que la
temperatura de saturación del vapor puro es 50°C. La alimentación ingresa a 45°C, el
funcionamiento es co-corriente y se calcula el ascenso ebulloscópico a medida que la solución
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 vacía porque no se ha contemplado en el ejemplo esa
opción. En la figura 5 se presentan los resultados.
Obsérvese que en todos los casos, la solapa activa es la que se destaca con letra en
negrita.
Diseño y simulación de Evaporadores Multiefecto – González et al. Pag. 11 de 12

Figura 2 – Ingreso de datos generales

Figura 3 – Ingreso de Datos del área de evaporación

Diseño y simulación de Evaporadores Multiefecto – González et al. Pag. 12 de 12

Figura 4 – Ingreso de Datos del Area de Precalentamiento (no incluída en este ejemplo)

Figura 5 – Resultados de la Simulación