Universidad Industrial de Santander
Escuela de Ingeniería Química
Laboratorio de procesos I código: 27542
Profesor: Humberto Escalante Hernández
Práctica N° 5: Evaporador en Termosifón.
ANÁLISIS Y DETERMINACIÓN DE LA EFICIENCIA DEL PROCESO DE
EVAPORACIÓN DE LA SOLUCIÓN SACAROSA-AGUA MEDIANTE
EVAPORADOR EN TERMOSIFÓN.
Grupo N° 3
Estudiantes:
David Leonardo Porras Albarracín 2162024
Ferney Sepúlveda Avila 2165183
Oscar Alberto Linares Velosa 2154614
Bucaramanga, marzo, 2020
1
1. INTRODUCCIÓN
Uno de los procesos más utilizados en la industria es la evaporación que es el método por el
cual se elimina una parte del agua contenida en un fluido a través de la transferencia de calor
obtenida por una condensación, con el fin de tener una solución más concentrada mediante la
ebullición. La solución a concentrar va a determinar el tipo de evaporación, temperatura y
presión a las que se realizará el proceso para obtener una mayor eficiencia. Los mecanismos
de transferencia de calor presentes en el fenómeno de evaporación son la conducción y
convección, siendo esta última la predominante.
Para esta práctica se trabajó el equipo de evaporador de efecto termosifón, se preparó una
solución agua-sacarosa depositando en la cámara de vaporización la solución a concentrar
correspondiente a 4 L de agua con 150 g de azúcar y se colocó en el puerto de llenado una
termocupla que permitió registrar la temperatura de la solución trabajada. Luego de ello, se
abrió el flujo de vapor que provenía de la caldera y se ajustó la presión con el valor de 8 psi
para el experimento.
Seguidamente, se inició con la recolección de 10 muestras de la solución cada 5 minutos a
través de un puerto del equipo para calcular la concentración de azúcar y se hizo un registro de
datos para las temperaturas marcadas por el sensor que posee la unidad de evaporación.
Al final de la práctica se determinan las eficiencias y pérdidas de calor en el condensador y
evaporador; además se realiza el análisis del cambio de concentración de soluto en una solución
a través del tiempo, cálculo para el coeficiente global de transferencia de calor y la
determinación de los grados BRIX de una solución.
2
2. OBJETIVOS
2.1 Objetivo general.
2.1.1 Evaluar y analizar el funcionamiento del evaporador de efecto simple
mediante un proceso de concentración de sacarosa
2.2 Objetivos Específicos
2.2.1 Calcular las pérdidas de calor en el evaporador y en el condensador.
2.2.2 Determinar las eficiencias del evaporador y el condensador.
2.2.3 Determinar el coeficiente global de transferencia de calor.
2.2.4 Estimar la relación de recirculación del líquido debido al efecto termosifón.
2.2.5 Analizar el cambio de la concentración de soluto a través del tiempo.
3
3 REGISTRO DE DATOS
3.1 Datos de la curva de calibración
Concentración Grados Brix
50 0,35
25 0,2
12,5 0,1
6,25 0,055
Tabla 1. Registro de datos de concentración y grados de Brix durante la práctica.
3.2 Datos obtenidos del Evaporador
Datos tomados de temperatura y volumenes de agua condensada
Vol. Vapor
Vol. Vapor
t (min) T1 (°C) T2 (°C) T3 (°C) T4 (°C) T5 (°C) T6 (°C) condensado
Caldera (mL)
(mL)
5 31 30 32 116 115 87 0 1050
10 31 35 33 117 115 97 58 820
15 31 33 33 116 115 98 200 950
20 31 37 33 116 115 98 230 950
25 31 36 33 116 116 98 290 1080
30 31 37 33 116 113 98 230 850
35 31 37 33 116 115 98 220 960
40 31 37 33 116 115 98 240 100
45 31 37 33 116 115 98 210 980
50 31 37 33 116 115 98 232 960
Promedio 31 35,6 32,9 116,1 114,9 96,8 Suma 1910 Suma 8700
Tabla 2. Registro de datos de temperaturas y volúmenes del Agua condensada durante la práctica.
3.3 Concentraciones salientes
Muestra Tiempo (min) Grados Brix
1 5 4
2 10 4,5
3 15 5
4 20 5
5 25 6
6 30 6
7 35 6,5
8 40 7
9 45 7,5
10 50 8
Tabla 3. Registro de datos de Concentraciones salientes y grados Brix durante la práctica.
4
4 MODELO DE CÁLCULO
4.1 Cálculo de las eficiencias
4.1.1 Eficiencia en el evaporador (𝒏𝑒 ):
𝑄𝑠
𝑛𝑒 (%) = ∗ 100 𝐸𝑐𝑢𝑎𝑐𝑖ó𝑛 (1)
𝑄𝑚𝑎𝑥
𝑄𝑠 = (𝑉𝑖𝑛𝑖𝑐𝑖𝑎𝑙 ∗ 𝜌𝑠𝑙𝑛 ∗ 𝐶𝑝𝑠𝑙𝑛 ∗ (𝑇𝑓 − 𝑇𝑖 ) + (𝑉𝑖𝑛𝑖𝑐𝑖𝑎𝑙 − 𝑉𝑓𝑖𝑛𝑎𝑙 ) ∗ 𝜌𝑎𝑔𝑢𝑎 ∗ λ𝑠(98°𝐶) ) 𝐸𝑐𝑢𝑎𝑐𝑖ó𝑛 (2)
𝐶𝑝𝑠𝑙𝑛 = (1 − 0.0056 ∗ °𝐵𝑟𝑖𝑥 ) ∗ 𝐶𝑝𝐻2𝑂 𝐸𝑐𝑢𝑎𝑐𝑖ó𝑛 (3)
𝑄𝑚𝑎𝑥 = 𝐹2 ∗ λ𝑣 𝐸𝑐𝑢𝑎𝑐𝑖ó𝑛 (4)
𝐹2 = 𝑉3 ∗ 𝜌𝐻2𝑂 𝐸𝑐𝑢𝑎𝑐𝑖ó𝑛 (5)
Sea:
λ𝑣 = Calor latente de evaporación del agua.
𝑉𝑖𝑛𝑖𝑐𝑖𝑎𝑙 = Volumen inicial de solución.
𝑉𝑓𝑖𝑛𝑎𝑙 = Diferencia entre el volumen inicial introducido al evaporador y el volumen de agua
condensada.
𝑇𝑓 = Temperatura promedio T6.
𝑇𝑖 = Temperatura ambiente.
Otra forma para calcular la eficiencia del evaporador es mediante el cálculo de la economía de
la evaporación esta se calcula a partir de la masa de vapor consumida proveniente de la caldera.
𝑉1 ∗ 𝜌𝐻2𝑂
%𝐸𝑐𝑜 = ∗ 100
𝑉2 ∗ 𝜌𝐻2𝑂
4.1.2 Eficiencia en el condensador (𝒏𝑐 )
𝑄𝑐
𝑛𝑐 (%) = ∗ 100 𝐸𝑐𝑢𝑎𝑐𝑖ó𝑛 (6)
𝑄𝑚𝑎𝑥𝑐𝑜𝑛𝑑
𝑄𝑐 = 𝑞 ∗ 𝜌𝐻2 𝑂 ∗ 𝐶𝑝𝐻2 𝑂 ∗ (𝑇2 − 𝑇1) ∗ 𝑡 𝐸𝑐𝑢𝑎𝑐𝑖ó𝑛 (7)
𝑄𝑚𝑎𝑥𝑐𝑜𝑛𝑑 = 𝑉1 ∗ 𝜌𝐻2𝑂 ∗ (λ𝑠 + 𝐶𝑝𝐻2𝑂 ∗ (𝑇6 − 𝑇3 )) 𝐸𝑐𝑢𝑎𝑐𝑖ó𝑛 (8)
Sea:
𝑞 = Caudal de agua 𝑡 = Tiempo de 𝜌𝐻2𝑂 = Densidad
de enfriamiento operación del agua en T2
4.2 Cálculo de las pérdidas de calor
4.2.1 Perdidas de calor en el evaporador 𝑸𝒑𝒆
𝑄𝑝𝑒 = (𝑄𝑚𝑎𝑥 − 𝑄𝑠 ) 𝐸𝑐𝑢𝑎𝑐𝑖ó𝑛 (9)
4.2.2 Pérdidas de calor en el condensador 𝑸𝒑𝒄
𝑄𝑝𝑐 = (𝑄𝑚𝑎𝑥𝑐𝑜𝑛𝑑 − 𝑄𝑐 ) 𝐸𝑐𝑢𝑎𝑐𝑖ó𝑛 (10)
5
4.3 Coeficiente global de transferencia de calor en (Q)
𝑄
𝑈= 𝐸𝑐𝑢𝑎𝑐𝑖ó𝑛 (11)
𝐴 ∗ ∆𝑇𝑐
∆𝑇1 − ∆𝑇2
∆𝑇𝑐 = 𝐸𝑐𝑢𝑎𝑐𝑖ó𝑛 (12)
∆𝑇
𝐿𝑛(∆𝑇1 )
2
Sea:
𝑄 = Velocidad de transferencia de calor en el evaporador.
𝐴 =Area de transferencia de calor en el evaporador. Que toma un valor de 0,087776 𝑚 2, obtenido del
diámetro externo de 2 in con una longitud de 55 cm, así el área es 0,087776 𝑚 2
∆𝑇1 = Diferencia de temperaturas entre el vapor de la caldera entrando al evaporador T4 y la solución
a la salida del evaporador T6.
∆𝑇2 = Diferencia de temperaturas entre el vapor de la caldera saliendo del evaporador T5 y la
solución a la entrada del evaporador T6.
4.4 Relación de recirculación (R)
𝑉1 ∗𝜌𝐻2𝑂 −𝑉2 ∗𝜌𝑠𝑙𝑛
𝑅= 𝐸𝑐𝑢𝑎𝑐𝑖ó𝑛 (13)
𝑉1 ∗𝜌𝐻2𝑂
𝐷𝑅
𝑉𝑟 = 𝐸𝑐𝑢𝑎𝑐𝑖ó𝑛 (14)
𝑇
Sea:
𝑉1= Volumen de agua condensado.
𝑉2 =Volumen del vapor condensado proveniente de la caldera.
𝐷𝑅 =Distancia recorrida por el pigmento.
𝑇 =Tiempo que tarda el pigmento en recorrer el trayecto 𝐷𝑅 .
4.5 Cálculo de la concentración.
Para el cálculo de la concentración de la solución durante el tiempo de operación se debe hacer
previamente una curva patrón con soluciones de azúcar establecidas. Estas concentraciones se
calculan titulando una solución madre de 150g en 100 mL a mitad, hasta tener 4 titulaciones
las cuales se grafican, en el eje x grados Brix y en el eje y concentraciones (g/mL).
𝑉1 ∗ 𝐶1 = 𝑉2 ∗ 𝐶2 𝐸𝑐𝑢𝑎𝑐𝑖ó𝑛 (15)
6
5 RESULTADOS.
5.1.Cálculo de las eficiencias.
5.1.1 Cálculo de la eficiencia en el evaporador.
Para el cálculo de la eficiencia en el evaporador es necesario la relación entre el calor
aprovechado y el calor suministrado al equipo dada por la ecuación 1.
𝑄𝑠
𝑛𝑒 (%) = ∗ 100
𝑄𝑚𝑎𝑥
El calor suministrado se calculó con la ecuación 3, además debemos usar la ecuación 4 para
el cálculo del 𝐶𝑝𝑠𝑙𝑛 pues no se conoce.
𝑄𝑠 = (𝑉𝑖𝑛𝑖𝑐𝑖𝑎𝑙 ∗ 𝜌𝑠𝑙𝑛 ∗ 𝐶𝑝𝑠𝑙𝑛 ∗ (𝑇𝑓 − 𝑇𝑖 ) + (𝑉𝑖𝑛𝑖𝑐𝑖𝑎𝑙 − 𝑉𝑓𝑖𝑛𝑎𝑙 ) ∗ 𝜌𝑎𝑔𝑢𝑎 ∗ λ𝑠(98°𝐶) )
𝑘𝐽
𝐶𝑝𝑠𝑙𝑛 = (1 − 0.0056 ∗ °𝐵𝑟𝑖𝑥) ∗ 𝐶𝑝𝐻2𝑂 = (1 − 0.0056 ∗ 4) ∗ 4.186 = 4.092
𝑘𝑔°𝐶
𝑘𝑔 𝑘𝐽 𝑘𝑔 𝑘𝐽
𝑄𝑠 = (0.004𝑚3 ∗ 1000 ∗ 4.092 ∗ (96.8 − 29)°𝐶 + (0.004 − 0.00209)𝑚3 ∗ 960.342 3 ∗ 2260 = 5255.23 𝑘𝐽
𝑚3 𝑘𝑔°𝐶 𝑚 𝐾𝑔
𝑘𝑔 𝑘𝐽
𝑄𝑚𝑎𝑥 = 𝑉3 ∗ 𝜌𝐻2𝑂 ∗ λ𝑣 = 0.001050𝑚 3 ∗ 960.342 ∗ 2260
𝑚3 𝑘𝑔
𝑄𝑚𝑎𝑥 = 2278.89(𝑘𝐽)
Eficiencia en el evaporador
Grados
Muestra Tiempo (min) Brix 𝑪𝒑𝒔𝒍𝒏 Qs (KJ) Qmax (kJ) 𝑛𝑒 (%)
1 5 4 4,0922336 5255,22603 2278,89157 -
2 10 4,5 4,0805128 5252,04735 4058,59736 -
3 15 5 4,068792 5248,86867 6120,45163 85,7594992
4 20 5 4,068792 5248,86867 8182,30591 64,1490153
5 25 6 4,0453504 5242,51131 10526,3087 49,8038911
6 30 6 4,0453504 5242,51131 12371,1256 42,3769951
7 35 6,5 4,0336296 5239,33262 14454,6836 36,2466087
8 40 7 4,0219088 5236,15394 14671,7209 35,6887509
9 45 7,5 4,010188 5232,97526 16798,6864 31,151098
10 50 8 3,9984672 5229,79658 18882,2444 27,6969012
Tabla 4. Resultados de cálculo de la eficiencia en el evaporador.
7
5.1.2 Cálculo de la eficiencia del condensador.
De manera similar al evaporador, la eficiencia energética del condensador se logra determinar
por medio de la relación de calor adquirido por el fluido de enfriamiento con el calor máximo
que cede al vapor de la solución, como se logra observar en la ecuación 6.
Para realizar la ejemplificación se utiliza la muestra número 3.
332,827186
𝑛𝑐 (%) = ∗ 100 = 51,2613887%
649,274618
Para calcular el calor adquirido y el calor máximo se usaron las ecuaciones 7 y 8
respectivamente, de tal manera que;
𝑄𝑐 = 0,008 ∗ 993,87 ∗ 4,186 ∗ (33 − 31) ∗ 5 = 332,827186 [𝑘𝐽]
(58 + 200)
𝑄𝑚𝑎𝑥𝑐𝑜𝑛𝑑 = ∗ 993,87 ∗ (2260 + 4,186 ∗ (98 − 33)) = 649,274618 [𝑘𝐽]
1000000
Eficiencia en el condensador
Tiempo Grados
Muestra (min) Brix Qc (kJ) Qmax ne (%)
1 5 4 0 0 -
2 10 4,5 665,654371 145,719661 -
3 15 5 332,827186 649,274618 51,2613887
4 20 5 998,481557 1228,08532 81,3039238
5 25 6 832,067964 1957,89013 42,498195
6 30 6 998,481557 2536,70083 39,3614234
7 35 6,5 998,481557 3090,34586 32,3097026
8 40 7 998,481557 3694,32225 27,0274624
9 45 7,5 998,481557 4222,80159 23,6450029
10 50 8 998,481557 4806,64543 20,7729397
Tabla 5. Resultados de cálculos de la eficiencia en el condensador.
5.2 Cálculo de Pérdidas de calor.
En el cálculo de las pérdidas de calor para el caso del evaporador como del condensador, es la
resta del calor máximo con el calor adquirido por la solución como el calor adquirido por el
fluido de enfriamiento, respectivamente. Para la ejemplificación se tomaron los datos obtenidos
en la muestra número 3.
8
Para el evaporador se tiene la ecuación 9, así
𝑄𝑝𝑒 = (6120,45163 − 5279,455077)[𝑘𝐽] = 843,996557 [𝑘𝐽]
De igual manera para el condensador se tiene la ecuación 10,
𝑄𝑝𝑐 = (649,274618 − 332,8271856)[𝑘𝐽] = 316,447433[𝑘𝐽]
5.2.1 Cálculo de las pérdidas de calor en el evaporador.
Pérdidas de calor en el evaporador
Tiempo (min) Qs (kJ) Qmax (kJ) Qpe
5 5282,971373 2278,89157 -
10 5279,713225 4058,59736 -
15 5276,455077 6120,45163 843,996557
20 5276,455077 8182,30591 2905,85083
25 5269,938781 10526,3087 5256,36988
30 5269,938781 12371,1256 7101,18686
35 5266,680633 14454,6836 9188,00301
40 5263,422485 14671,7209 9408,29845
45 5260,164337 16798,6864 11538,5221
50 5256,906189 18882,2444 13625,3382
Tabla 6. Resultados de cálculo de las pérdidas de calor en el evaporador.
5.2.2 Cálculo de las pérdidas en el condensador.
Pérdidas de calor en el condensador
Tiempo (min) Qc (kJ) Qmax (kJ) Qpc
5 0 0 0
10 665,6543712 145,719661 -
15 332,8271856 649,274618 316,447433
20 998,4815568 1228,08532 229,603768
25 832,067964 1957,89013 1125,82216
30 998,4815568 2536,70083 1538,21928
35 998,4815568 3090,34586 2091,8643
40 998,4815568 3694,32225 2695,84069
45 998,4815568 4222,80159 3224,32003
50 998,4815568 4806,64543 3808,16387
Tabla 7. Resultados de cálculo de las pérdidas de calor en el condensador.
9
16000
14000
Calor perdido (KJ)
12000
10000
8000
Qevap
6000
Qcond
4000
2000
0
0 10 20 30 40 50 60
Tiempo (min)
Figura 1. Comparación de calores perdidos en el evaporador y condensador durante los 50 min.
5.3 Coeficiente global de transferencia de calor.
Para el coeficiente de trasferencia de calor usamos las ecuaciones 11 y 12, así:
19,3 − 18.1
∆𝑇𝑐 = = 18,683581 °𝐶
0,064193158
5242,8292 [𝑘𝐽] 𝑘𝐽
𝑈= 2
= 3195,1953 2
0,087776 𝑚 ∗ 18,683581°𝐶 𝑚 °𝐶
5.4 Relación de recirculación.
𝑉1 ∗ 𝜌𝐻2𝑂 − 𝑉2 ∗ 𝜌𝑠𝑙𝑛
𝑅=
𝑉1 ∗ 𝜌𝐻2𝑂
𝑘𝑔 𝑘𝑔
0.004𝑚3 ∗ 997.12 − 0 𝑚3 ∗ 1000 3
𝑅= 𝑚3 𝑚
𝑘𝑔
0.004𝑚3 ∗ 997.12 3
𝑚
𝑅=1
Cálculo de la relación de recirculación
V1 (m3) V2 (m3)ρH2O(T1) R
[kg/m3]
0,004 0 997,12 1
0,004 0,000058 997,12 0,985458119
0,004 0,000258 997,12 0,935313703
0,004 0,000488 997,12 0,877647625
0,004 0,000778 997,12 0,804938222
0,004 0,001008 997,12 0,747272144
0,004 0,001228 997,12 0,692113286
0,004 0,001468 997,12 0,631939987
0,004 0,001678 997,12 0,57928835
0,004 0,00191 997,12 0,521120828
Tabla 8. Resultados de cálculos de la relación de recirculación.
10
5.5 Análisis del soluto a través del tiempo
Grados
Concentración Brix
50 0,35
25 0,2
12,5 0,1
6,25 0,055
Tabla 9. Especificación de la concentración de la solución y sus grados Brix.
Curva de Calibración
0.4
0.35
0.3 y = 0.0067x + 0.0183
Grados Brix
0.25 R² = 0.9954
0.2
0.15
0.1
0.05
0
0 10 20 30 40 50 60
Concentración (g/mL)
Figura 2. Relación de Grados Brix vs Concentración de la solución.
𝐶𝑜𝑛𝑐𝑒𝑛𝑡𝑟𝑎𝑐𝑖ó𝑛 = 0.0067 ∗ 𝐺° + 0.0183
𝑔
𝐶1 = 0.00674 ∗ 4 + 0.01826 = 0.008703
𝑚𝐿
Concentración
Muestra Grados Brix
[g/mL]
1 0,04 0,008702609
2 0,045 0,012073043
3 0,05 0,015443478
4 0,05 0,015443478
5 0,06 0,022184348
6 0,06 0,022184348
7 0,065 0,025554783
8 0,07 0,028925217
9 0,075 0,032295652
10 0,08 0,035666087
Tabla 10. Especificación de la concentración de la solución de las 10 muestras y sus grados Brix.
11
Figura 3. Relación de Concentración de la solución vs tiempo.
6 CONCLUSIONES
Con base en lo anterior analizado se puede concluir lo siguiente:
• Se puede observar en las tablas 4 y 5 que, el cálculo de las eficiencias en el
evaporador y del condensador obtenidos fueron menores esto se debe a la pérdida
de energía que ocurre en el equipo evaporador.
• La medición de los grados Brix es una medida directa de la concentración de las
soluciones y facilita la observación del fenómeno de la evaporación ya que el uso
del refractómetro es sencillo y rápido, este es un método más práctico que otros
para la determinación de concentración de las soluciones.
• En la figura 1 se puede observar que la mayor cantidad de pérdidas de calor se
presenta en el evaporador a diferencia del condensador, pues se está liberando
calor hacia el ambiente, para esto se debe implementar un asilamiento que impida
las pérdidas de calor y de esta forma se pueda aprovechar la energía suministrada
y así aumente la eficiencia tanto en el evaporador y el condensador aumente.
• En la figura 3 se puede observar que a medida que el agua se va evaporando en el
termosifón la concentración de la solución aumenta con respecto al tiempo de
operación.
12
• Se determinó para el caso de estudio, que la relación de recirculación presenta un
comportamiento decreciente ya que el volumen final varía con la temperatura T1.
En los primeros 5 minutos la relación de recirculación es 1 ya que no se evapora
nada de la solución y a medida que pasa el tiempo esta relación disminuye porque
el volumen de la solución disminuye.
13
ANEXOS
REGISTRO DE DATOS
14
