Universidad Nacional de Colombia

Facultad de Ingenierı́a

Primer parcial. Transferencia de

Calor.

Prof: Gustavo Adolfo Orozco Alvarado

Presentan:
Kelly Johana Jáuregui Acevedo
Juan Felipe Mora Rios
Xiomara Lizeth Vélez Diaz

2021
 Parte 1
1. Realizar una discusión general de los mecanismos de transferencia de calor
en un evaporador.
dentro de un evaporador, dependiendo el tipo de sustancia problema a evaporar gene-
ralmente se realiza evaporación por convección libre, o forzada dependiendo del tipo
de evaporador, Un evaporador es aquella unidad de proceso qué utiliza un fluido como
medio de transferencia de calor. Este aprovecha el calor latente, o sensible, de un flui-
do para aportar calor latente a otro. es importante destacar que todo compuesto que
presente una transferencia de calor por convección también lo presenta por convección

2. Plantear el balance general de energı́a para múltiples unidades

Para plantear el balance de energı́a de múltiples unidades tomamos como ejemplo la
figura 6 es indispensable realizar un balance de materia como se propone a continuación.

Figura 1: Secuencia III-I-II numeración de lı́quido.

F= Corriente de alimentación L= Corriente de Concentrado S= Corriente de vapor

n= Número de efectos
Empezaremos por proponer el balance de materia para el primer efecto; no se presenta
generación, destrucción o reacción, por lo que el balance de materia para este caso serı́a
igual a la corriente de alimentación (lo que entra), igual a lo que sale, que para este
caso es la corriente de concentrados y la del vapor.

F = L1 + V1 (1)

1
Para el segundo efecto, tenemos en cuenta que la corriente de entrada serı́a la corriente
de concentrados que sale del primer efecto y eso seria igual a las salidas, para este caso,
la corriente de concentrados del efecto 2, más la corriente de de vapor del efecto 2.

L1 = L2 + V1 (2)

y ası́ generamos una expresión, para los balances de masa en 2 efectos, como se muestra
a continuación reemplazando la ecuación 2 en la 1

F − V1 − V2 = L2 (3)
Podemos ver una secuencia dentro de la ecuación número 3, simplificando para obtener
una ecuación general para n efectos nos queda como aparece a continuación
X
F− V i = Ln (4)
m

donde; vi =Salida de vapor desde el efecto m hasta el efecto n Li=Salida de concentrados

en el conjunto de efectos Balance de energı́a para un evaporador simple efecto
las corrientes de entrada para un evaporador simple efecto, como se mencionan en la
ecuación 5 son las corrientes de calor del evaporador, el flujo de alimentación, con su
respectiva entalpı́a y por último la entrada de vapor con su entalpı́a, las corrien

Q1 + F HFL + V1 HL1 = F HL1 + V H1V (5)

tes de salida comprenden el flujo del lı́quido de salida (Concentrados) con su entalpı́a,
y el flujo de vapor de salidaPmultiplicado por la entalpı́a, sabemos que por el balance de
masa, ecuación 4 Li = F − Vi , en este caso pretendemos formular un sistema de ecua-
ciones donde las incógnitas sean las corrientes de vapor. ası́ que buscamos simplificar
la ecuación, a continuación se mostrarán los pasos consecutivos para la simplificación.

Q1 + F HFL = F HL1 − V1 HL1 + V H1V (6)

Q1 + V1 (H1L − H1V ) = F (H1L − HFL ) (7)

Q1 + F L1 = LH1L + V H1V (8)

para realizar el balance de energı́a para un evaporador multi-efecto similarmente uti-

lizamos la ecuación numero 5, pero en lugar de los subı́ndices 1, en este caso vamos a
utilizar los subı́ndices i que se refieren a cualquier efecto de interés dentro del evapo-
rador, además tenemos en cuenta que i-1 pertenece al evaporador anterior, la corriente
en este caso será reemplazada por L1 que significa los concentrados del evaporador
anterior.
L
Qi + Li Hi−1 + Vi HiL = Li HiL + Vi HiV (9)

2
 de una manera general para más evaporadores tenemos las siguientes ecuaciones para
simplificar el balance de energı́a

L
Qi + Li Hi−1 + Vi HiL − Li HiL = Vi HiV (10)
X X
L
Qi + (F − Vm )Hi−1 + Vi HiL − Vi HiV = (F − Vm )HiL + (11)
X X
L L
Qi + F iHi−1 − VM Hi−1 + Vi HiL = Fi HiL − Vm HiV (12)
X
Qi + Vm (HiL − Hi−1
V
) + Vi [HiL − HiV ] = F (HiL − Hi−1
V
) (13)
Salvo por el ultimo efecto n, en este el balance de manera análoga corresponde a:
 
L
X
V xf L L
Qi − Hn−1 Vm − Vn Hn = F H − Hn−1 (14)
xs n

3. Hacer un estimado del coeficiente global de transferencia de calor en cada

evaporador. Determinar las resistencias involucradas y hacer un estudio de
cuales pueden ser despreciables.
El coeficiente global de transferencia de calor de un evaporador de pelı́cula descendente
se puede calcular experimentalmente o mediante la ecuación (Mackereth,1992):
1 1 Ss Sw Sf 1
= + + + + (15)
U hi ks kw kf he
Donde
U = Coeficiente global de transferencia.
hi = Coeficiente de transferencia de calor interno.
he = Coeficiente de transferencia de calor externo.
Ss = Grosor de la capa de sarro del lado del vapor de calentamiento.
Sw = Grosor de la pared
Sf = Grosor de las incrustaciones del lado por donde pasa la leche.
ks = Conductividad térmica del sarro
ks = Conductividad térmica de la pared.
ks = Conductividad térmica de las incrustaciones.

3
Figura 2: Perfil de temperatura a través de la pared del tubo de calandria (Mackereth,1992).

El hi y he son los que hacen referencia a la transferencia de calor convectivo que

existe tanto por dentro como por fuera de un tubo del evaporador, uno para la leche
evaporándose y otro para el vapor de calentamiento condensándose.
El sarro hace referencia a aquellas impurezas que pueden llegar a ingresar junto con el
vapor de calentamiento, debidas a por ejemplo la presencia de oxido en el suministro
de vapor, estas formaran una capa por la cual el calor tendrá que pasar.
Las incrustaciones hacen referencia a las proteı́nas y minerales de la leche que se quedan
en pegadas en las paredes de los tubos, también puede deberse a un mal flujo de la leche,
si es insuficiente la pelı́cula se romperá y ocasionara que se queden secciones de tubo
expuestas esto generara una disminución en la eficiencia de la transferencia de calor y
que se formen ya que si el flujo no es uniforme en todo el tubo habrán partes del flujo
donde este se haga muy concentrado, en consecuencia muy viscoso, que eventualmente
se endurece y puede bloquear el tubo, además de que puede propiciar la proliferación de
microorganismos que no se quiere en la leche producto.En la figura 3 se puede apreciar
de mejor manera la capa formada por leche entera después de 22 horas de trabajo en
un evaporador.

4
 Figura 3: Leche endurecida en un evaporador (Broome, 2005).

W
La conductividad térmica de este material esta entre 0.5-3 Km 2 la cual es mucho menor

que la del material de los tubos, razón por la cual incluso una capa muy delgada es
capaz de reducir el coeficiente global de transferencia de calor. En lo que respecta a la
W
conductividad de la pared del tubo como tal, el acero inoxidable tiene k = 13.4 Km 2,

para minimizar la resistencia y la formación de depósitos se utilizan tubos lisos de 1mm

de grosor (Mackereth,1992).
Existen distintas maneras de calcular los coeficientes de transferencia de calor tanto
internos como externos, para el presente calculo se usó la siguiente ecuación para el
coeficiente interno del efecto i que pertenece al conjunto de efectos I:
 V L V
1
v 3 3
gρi (ρ i − ρi )(κi )
hi = 0,002(ReVi )0,6 ,i ∈ I (16)
(µvi )2
donde
g = Gravedad
ReVi = Es el numero de Reynolds del vapor producido del efecto i.
µvi = Es la viscosidad del vapor producido del efecto i.
ρi = Es la densidad ya sea del vapor o de liquido de leche producidas en el efecto
i.
κvi = Conductividad térmica del vapor generado en el efecto i
Para poder estimar las propiedades que se solicitan en la leche se utilizara para la
densidad la ecuación propuesta por (articulo profe) y para la viscosidad de la leche se
utilizara la ecuacion propuesta por (Fernández-Martı́n, 1972):

Figura 4: Viscosidad dinámica de la leche.

5
Donde s es el porcentaje en masa de solidos totales presente y t corresponde a la
temperatura a la cual se pretende leer el dato. Para la información correspondiente al
agua se leerá de tablas.
Este numero de Reynolds se calcula:
4Vi
(ReVi ) = (17)
πΦi µvi
donde Φi es el diámetro del evaporador.
Para el calculo del coeficiente exterior:
 2 2 0,25
gρli κli ∆hcond
he = 0, 934 ,i ∈ I (18)
µli ∆T Lo

donde

g = Gravedad
ρli = Es la densidad del condensado que genera el vapor de calentamiento del
efecto i.
µli = Es la viscosidad del condensado que genera el vapor de calentamiento del
efecto i..
∆hcond = Es la entalpı́a de condensación para el vapor de calentamiento que entro
al efecto i .
∆T = Es la diferencia de temperatura entre el vapor de calentamiento y la pared
del tubo del evaporador
L = Es la longitud de los tubos del evaporador

Para la verificación del calculo efectuado del coeficiente de transferencia de calor glo-
bal se tendrá en cuenta las expresiones que propuso (Broome, 2005) en función de la
cantidad total de solidos presentes a ala salida del efecto :

Para leche desnatada la correlación es valida cuando solidos totales estén entre 11 y 50 %
y para la leche entera cuando estén entre 15 y 50 %. Usando las expresiones mencionadas
anteriormente, para una temperatura de presión de 269°C y presión de 0,4053 bar, con
la composición de la leche mencionada en el inciso numero 6 del presente documento
y una concentración de solidos a la salida de 29,5343 %, los coeficientes convectivos
fueron:

6
 hi = 3.5328e-10 mW2 K y he = 753.1339 mW2 K
Tomando estos valores y despreciando los términos correspondientes a las incrustaciones
tanto dentro del tubo como por fuera del mismo el coeficiente global de transferencia
de calor U corresponde a 732.5476 mW2 K este calculo fue realizado en matlab (para mas
detalles del código remitirse al anexo) las propiedades del vapor de agua fueron to-
madas de tablas disponibles en la literatura. El dato que resulta de usar la expresión
desarrollada por Broome para la leche entera es U = 1775 mW2 K debido a que el error
es de mas del 50 % se cree mas conveniente realizar los cálculos con las expresiones de-
sarrolladas por Broome. Para el calculo del coeficiente global de transferencia de calor
se tuvieron en cuenta solo estos aportes primeramente porque no se tenia información
sobre el como calcular el grosor de las capas de las incrustaciones, también se tuvo en
cuenta el método de calculo propuesto por (Madoumier,2018)

4. Discutir acerca del tipo de evaporadora usar, haciendo un cuadro compara-

tivo y basado en criterios de selección escoger el evaporador.
Para realizar esta selección nos referimos a la tabla 1 que se encuentra a continuación,
Para nuestro diseño es indispensable que el evaporador tenga una temperatura sensible,
debido a que si la leche pasa una temperatura indicada puede desproteinizarse, además
también tiene que poseer la habilidad de manejar sólidos o cristales, debido a que
la leche se le retirará agua, esto puede producir que se cristalice, o se solidifique.El
propósito de este trabajo es encontrar un sistema de evaporación que contenga una
alta capacidad, debido a que será utilizado en una fabrica de leche en polvo, ası́ este se
convierte también en un aspecto de gran importancia, gracias a esto, podemos ver que
los evaporadores de pelı́cula descendente, bandejas, y pelı́cula delgada agitada
cumplen con nuestros objetivos.
La decisión final de la elección de evaporadores conviene debido que gracias a la lite-
ratura, para la evaporación de leche, es indispensable utilizar evaporadores de pelı́cula
descendente como se muestra en la imagen ??.

7
Tipo de limpio Alta SolidosAbordaje Espuma Temperaturaviscosidad
evapora- ca- o Sensible
dor paci- cris-
dad tales

Batch x x
Tubos Ho- x
rizontales
tubos Ver- x
ticales cor-
tos
tubos Ver- x x
ticales lar-
gos
pelı́cula x x x
ascendente
y descen-
dente
pelı́cula x x x x x
descenden-
te
Circulación x x x x x
forzada
Lámina x x x x x x x
pelı́cula x x x x x x x
delgada
agitada

Cuadro 1: caracterı́sticas para la selección de criterios para evaporadores

8
 Figura 5: Evaporador pelı́cula descendente elegido para la realización del proyecto

5. Hacer una discusión acerca del modelo termodinámico.

Para la elección del modelado termodinámico tenemos en cuenta los siguientes paráme-
tros.

La leche es un compuesto a-polar

Se compone por componentes Reales (También por pseudo-componentes como la
caceı́na y la prolina, pero no se pudo modelar en Aspen, por lo tanto reelegimos el
modelo sin estas sustancias asimilando las demás como puras)

Gracias a la figura 6 evidenciar que los modelos apropiados a seguir serı́an Peng-
Robinson, Redich-Kwong-soave, Lee-Kesler-pockler

9
 Figura 6: Elección modelo termodinámico.

Para nuestra simulación terminamos por elegir el modelo de peng robinson debido a
que era el que menos errores presentó en la modelación.

6. ¿Como va a modelar el compuesto leche?

Gracias a la literatura pudimos encontrar varios modelados para la leche, simulamos
en Aspen como una solución acuosa de grasas, carbohidratos,, minerales y proteı́nas,
realmente la leche consiste en cientos de productos quı́micos pero el acercamiento es
muy preciso al verdadero valor de la leche, por ejemplo propiedades como la densidad,

10
la viscosidad dinámica, y Capacidad calorı́fica especı́fica tienen un error de desviación
de menos del 10 %

Componente Porcentaje en peso

Agua 90.95
Ácido palmı́tico 1.40
ácido oleico 1.40
Lactosa 3.30
Calcio 0.28
Potasio 0.27
Prolina 2.40

Cuadro 2: Composición de compuestos de porcentaje en peso para leche

Es importante destacar que el compuesto de proteı́na que en este caso comprendı́a la

prolina, no pudo ser modelado en aspen, por lo que para esta situación lo que hicimos
fue obtener el valor de composición de la prolina y lo dividimos en los demás menos el
agua.

Parte 2
a) Realice todas las suposiciones necesarias y partiendo de un sistema
completamente especificado (g.libertad=0), de forma manual resuelva
el proceso para un evaporador simple y uno doble efecto. Determine la
economı́a del(os) evaporador(es), y el área

Figura 7: Planteamiento del problema para un efecto.

11
 Balance de energı́a

HA + Q = HB + HC (19)
Kg Kj KJ
HB = mCp∆T = 4,00686 ∗ 3930 ∗ (95 − 269)k = −2,74E6 (20)
h Kg ∗ k h
KJ Kg KJ
HC = 2605,13 ∗ 25,9931 = 67715,4 (21)
Kg h h
Kg KJ KJ
HA = 30 ∗ −13973,4 = −419202 (22)
h Kg h

Despejamos el calor

KJ
Q = 2,25308E6
h
fm=845.392Kg h
E=28.1767
A=7.295m2
(23)

Figura 8: Planteamiento del problema para dos efectos.

12
 Figura 9: Balance para 2 efectos.

Figura 10: Cálculos de calor, flujo másico, economı́a, y área para 2 efectos.

b) Utilizando el simulador:
Compare la evaporación en simple y varios efectos en configuración paralelo
y Contracorriente. Grafique la economı́a como función del numero de efectos
cuál es el número óptimo de efectos.

13
 Efectos en configuración en paralelo
Para los efectos en configuración en paralelo, fuimos tomando el valor de
economı́a a medida de que se agregaban los efectos, decidimos graficar 5 efectos
figura debido a que si aumentamos la cantidad de efectos el porcentaje de agua
en los concentrados no disminuı́a.

Figura 11: Modelación en aspen evaporadores multiefecto con circulación paraelela.

Graficamos la economı́a Vs el número de efectos y podemos observar que el

número óptimo de efectos es 4

Figura 12: Economı́a vs cantidad de efectos para una configuración en paralelo.

Efectos en configuración en Contracorriente

Para los efectos en configuración en Contracorriente, fuimos tomando el valor
de economı́a a medida de que se agregaban los efectos, decidimos graficar 5

14
 efectos figura debido a que si aumentamos la cantidad de efectos el porcentaje
de agua en los concentrados no disminuı́a.

Figura 13: Modelación en aspen evaporadores multiefecto con circulación Contracorriente.

Graficamos la economı́a Vs el número de efectos y podemos observar que el

número óptimo de efectos es 5, por este método es mucho mayor el costo
de inversión que si se hace por evaporadores multiefecto en paralelo que en
contracorriente

Figura 14: Economı́a vs cantidad de efectos para una configuración en contracorriente.

Utilice parte del vapor de salida de la última unidad para precalentar el ali-
mento y concluya que pasa con la economı́a al hacer esto

Efectos en configuración en paralelo con precalentamiento

Para los efectos en configuración en paralelo, fuimos tomando el valor de
economı́a a medida de que se agregaban los efectos, decidimos graficar 5 efectos

15
 figura debido a que si aumentamos la cantidad de efectos el porcentaje de agua
en los concentrados no disminuı́a.

Figura 15: Modelación en aspen evaporadores multiefecto con circulación paraelela con pre-
calentamiento.

Graficamos la economı́a Vs el número de efectos y podemos observar que el

número óptimo de efectos es 1

Figura 16: Economı́a vs cantidad de efectos para una configuración en paralelo con precalen-
tamiento.

Efectos en configuración en Contracorriente con precalentamiento

Para los efectos en configuración en Contracorriente, fuimos tomando el valor
de economı́a a medida de que se agregaban los efectos, decidimos graficar 5
efectos figura debido a que si aumentamos la cantidad de efectos el porcentaje
de agua en los concentrados no disminuı́a.

16
Figura 17: Modelación en aspen evaporadores multiefecto con circulación Contracorriente
con precalentamiento.

Graficamos la economı́a Vs el número de efectos y podemos observar que el

número óptimo de efectos es 1, por este método es mucho mayor el costo
de inversión que si se hace por evaporadores multiefecto en paralelo que en
contracorriente

Figura 18: Economı́a vs cantidad de efectos para una configuración en contracorriente con
precalentamiento.

Basado en simulaciones, explique que pasa con la presion de operación de los

evaporadores, ¿ como la escogió ?, ¿Qué estrategia puede proponer para hacer
el proceso mas eficiente?
Se puede observar teóricamente que en los evaporadores multiefecto en dispo-
sición paralela para la presión de cada uno de ellos se espera que aumente a
medida que hay más evaporadores conectados, para los evaporadores multi-
efecto de configuración contracorriente, se espera que la presión aumentee de
igual manera pero no como la de la configuración en paralelo como se puede
observar en las siguientes figuras

17
 Figura 19: Cambio de las propiedades del sistema en cada evaporador sistema paralelo.

Figura 20: Cambio de las propiedades del sistema en cada evaporador sistema contracorriente

De ser posible optimice la presion en los efectos para maximizar la produccion

de vapor (existe una función en el simulador llamada optimizador)
La presión en los evaporadores dependı́a del método a usar, en esta caso usa-
mos Peng-Robinson, en la literatura se nos indicaba que no se pueden usar
presiones ni temperaturas muy altas, esto se comprobaba al cambiar la presión
de los evaporadores por unas de mayor magnitud, arrojaba en la simulación
errores para realizar los respectivos cálculos y se paraba el proceso. por otro la-
do en la literatura se anunciaba que en procesos donde se usaba como alimento
la leche, las presiones deberı́an ser menores a la atmosférica para garantizar
una mayor evaporación del agua que entraba al proceso, además de asegurar
una correcta composición de la leche simulada.
En el momento de usar la herramienta optimización, las presiones de los eva-
poradores daban demasiado similares a las establecidas inicialmente; por tanto
la generación de vapor era casi la misma, por tan razón decidimos mantener

18
 las presiones iniciales ya que no se notaba mucha diferencia. si no hubiéramos
encontrado en la literatura seguramente si habrı́a un cambio mayor al usar
esta herramienta.
Proponga estrategias adicionales y basado en simulaciones muestre que tan
efectivas son.
Como estrategia se puede emplear el pre-calentamiento de la leche para de este
modo lograr disminuir la concentración de agua en el producto final, además se
pueden instalar válvulas para disminuir mas la presión y aumentar la generación
de vapor. realizar un correcto enfriamiento del producto para conservar todas las
concentraciones. .

Referencias
c) franco, I., nd DISEÑO DE UN SIMULADOR POR COMPUTADORA DE PRO-
CESOS DE EVAPORACIÓN EN UNA LÍNEA DE EVAPORADORES DE MÚLTI-
PLE EFECTO.
d ) Ang, K., 2011. Investigación de las propiedades reológicas de la leche concentrada
y el efecto de estas propiedades en el flujo dentro de los evaporadores de pelı́cula
descendente. , K., 2011. Investigación de las propiedades reológicas de la leche
concentrada y el efecto de estas propiedades en el flujo dentro de los evaporadores
de pelı́cula descendente.
e) Broome, S. R. (2005). Liquid Distribution And Falling Film Wetting In Dairy Eva-
porators. Department of Chemical and Process Engineering. https://drive.google.com/file/d/1-
Q6ojC8CA3k4tvKjnJ7iuB7y5CgPumGd/view
f ) Ang, K., 2011. Investigación de las propiedades reológicas de la leche concentrada
y el efecto de estas propiedades en el flujo dentro de los evaporadores de pelı́cula
descendente
g) Autores, V., 2018. Evaluación de variables de operación de evaporadores verticales
mediante herramienta computacional.
h) Martial Madoumier, Catherine Azzaro-Pantel, Gaelle Tanguy-Sai, Genevieve Gesan-
Guiziou. Modelling the properties of liquid foods for use of process flowsheeting
simulators: Appli-cation to milk concentration. Journal of Food Engineering, El-
sevier, 2015,
i ) AUTORES, V., 2013. Desarrollo de componentes hipotéticos para la simulación
del proceso de la leche utilizando un simulador de proceso comercial.
j ) AUTORES, V., 2014. Modelización de un evaporador de pelı́cula descendente para
procesos lácteos
k ) Dı́az Ovalle, C. O., Galván Ángeles, E. (2013). METODOLOGÍA PARA EL
DISEÑO ÓPTIMO DE EVAPORADORES DE PELÍCULA DESCENDENTE.
Avances en Ciencias e Ingenierı́a. https://drive.google.com/file/d/197VLDT26r7059yG4HMuJtz

19
.

20