Intercambiadores
Intercambiadores
Comisión 3
Resumen
Se varió el caudal del fluido caliente y se trabajó con un único caudal de fluido frío.
También se tomaron los valores de la temperatura en la entrada, en un punto intermedio y
en la salida del intercambiador. Utilizando estos datos, se analizaron los efectos que genera
un cambio en la configuración de flujo sobre los perfiles térmicos, el calor transferido y las
eficiencias. Se calcularon los coeficientes de transferencia de calor experimentales, y se
compararon con los obtenidos teóricamente mediante dos métodos: el método NTU y de
resistencias térmicas. Se pudo ver que el coeficiente global experimental se asemeja al
calculado por el método NTU, pero se diferencia de los calculados por el método de
resistencias.
Abstract
1
Trabajo Práctico N°3
Comisión 3
Introducción
Los intercambiadores de calor son equipos que facilitan el intercambio de calor entre
dos fluidos que se encuentran a temperaturas diferentes y evitan al mismo tiempo que se
mezclen entre sí. En la práctica, los intercambiadores de calor son de uso común en
diversas aplicaciones ingenieriles.
Un intercambiador de calor está relacionado con dos fluidos que fluyen separados
por una pared sólida. En primer lugar, el calor se transfiere del fluido caliente hacia la pared
por convección, después a través de la pared por conducción y, por último, de la pared
hacia el fluido frío de nuevo por convección. Éstos son habitualmente clasificados de
acuerdo al arreglo de flujo y tipo de construcción [1].
2
Trabajo Práctico N°3
Comisión 3
(a) (b)
Entonces, el calor transferido desde cada corriente puede ser calculado en función
de sus temperaturas utilizando la ecuación 2 y 3 [1].
3
Trabajo Práctico N°3
Comisión 3
temperaturas (K).
𝑄𝑝 = 𝑄𝑐 − 𝑄𝑓 Ec. 4
Una vez obtenidos los calores intercambiados por ambas corrientes, se pueden
calcular la eficiencia global (𝜂𝑡), la eficiencia de temperatura media (𝜂m) y la eficiencia de
temperatura para el fluido caliente (𝜂𝑐) y frío (𝜂𝑓) mediante las Ecuaciones 5,6,7,8 [1].
𝑄𝑓
𝜂𝑡 = 𝑄𝑐
. 100% Ec. 5
𝑇𝑐,𝑒 −𝑇𝑐,𝑠
𝜂𝑐 = . 100% Ec. 6
𝑇𝑐,𝑒−𝑇𝑓,𝑒
𝑇𝑓,𝑠−𝑇𝑓,𝑒
𝜂𝑓 = . 100% Ec. 7
𝑇𝑐,𝑒 −𝑇𝑓,𝑒
η𝑐 +η𝑓 Ec. 8
𝜂𝑚 = 2
. 100%
4
Trabajo Práctico N°3
Comisión 3
𝐴 = 𝜋. 𝑑𝑚 . 𝐿 Ec. 11
𝑑𝑜+𝑑𝑖 Ec. 12
𝑑𝑚 = 2
𝑄𝑐
𝑈 = Ec. 13
𝐴.𝛥𝑇𝑚𝑙
𝑑
𝐿𝑛( 𝑜 )
𝑖
𝑅𝑝𝑎𝑟𝑒𝑑 = 2.𝜋.𝐾.𝐿
𝑑 Ec. 14
𝑑
𝐿𝑛( 𝑜 )
1 𝑑𝑖 1
𝑅 = 𝑅𝑡𝑜𝑡𝑎𝑙 = 𝑅𝑖 + 𝑅𝑝𝑎𝑟𝑒𝑑 + 𝑅𝑜 = + + Ec. 15
ℎ𝑖.𝐴𝑖 2.𝜋.𝐾.𝐿 ℎ𝑜.𝐴𝑜
𝛥𝑇
𝑄 = 𝑈𝐴𝑠 𝛥𝑇 = 𝑈𝑖 𝐴𝑖 𝛥𝑇 = 𝑈𝑜 𝐴𝑜 𝛥𝑇 = Ec. 16
𝑅𝑡𝑜𝑡𝑎𝑙
Ec. 17
5
Trabajo Práctico N°3
Comisión 3
1
𝑈𝑡𝑒ó𝑟𝑖𝑐𝑜 = 𝐴.𝑅𝑡𝑜𝑡𝑎𝑙
Para obtener el Uteórico, se necesita tener los coeficientes convectivos de cada fluido.
Éstos se pueden determinar sabiendo la condición de flujo mediante el Reynolds (Re) [1].
𝑝.𝑣.𝐷
𝑅𝑒 = μ
Ec. 18
mi Ec. 19
𝑣=
At . 𝑝
ℎ𝑛 . 𝐷𝑛 Ec. 20
𝑁𝑢 =
𝐾
1
𝑁𝑢 = 𝑅𝑒 0,8 𝑃𝑟 3 Ec. 21
10000< Re
6
Trabajo Práctico N°3
Comisión 3
𝑓
(𝑅𝑒𝐷 −1000)𝑃𝑟
𝑁𝑢 = 8
1
Ec. 22
2
𝑓 2
1+12,7( ) (𝑃𝑟 3−1)
8
dónde las propiedades son evaluadas a una temperatura media y el factor de fricción
f para tubos lisos puede determinarse a partir de una relación apropiada como la que se
muestra a continuación
𝑄
ε=𝑄 Ec. 24
𝑚𝑎𝑥
7
Trabajo Práctico N°3
Comisión 3
𝐶𝑖 = 𝑚𝑖 . 𝐶𝑝𝑖 Ec. 27
transferencia, NTU. Para valores específicos de U y C mín, el valor del NTU es una medida
del área superficial de transferencia de calor, As. Por ende, entre mayor sea el NTU, más
grande es el intercambiador de calor. En el análisis de los intercambiadores de calor
también resulta conveniente definir otra cantidad adimensional llamada relación de
capacidades c como:
𝐶𝑚𝑖𝑛
𝐶= Ec. 28
𝐶𝑚𝑎𝑥
Co - corriente
−𝑙𝑛(1−𝑒(1+𝑐))
𝑁𝑇𝑈 = Ec. 29
1+𝑐
Contracorriente
1 ε−1
𝑁𝑇𝑈 = 𝑐−1
(𝑙𝑛( ε.c−1) Ec. 30
8
Trabajo Práctico N°3
Comisión 3
Experimentos y resultados.
9
Trabajo Práctico N°3
Comisión 3
Luego, se pesó el contenido de agua del recipiente y dividiendo este valor por el
tiempo medido se obtuvo el caudal de operación. Lo mismo se realizó para medir el caudal
de la corriente fría, con la diferencia que, al trabajar siempre con el mismo caudal, solo se lo
midió una vez por configuración. Cabe aclarar que este procedimiento se repitió hasta
obtener tres valores de cada caudal que tuviesen un error menor al 5% y, cuando al llegar a
esos tres valores, se calculó el promedio entre ambos. Dichos valores se encuentran en la
Tabla A1 del anexo.
En las figuras 5 y 6 se representan los datos obtenidos de temperatura a lo largo del equipo
para cada configuración.
10
Trabajo Práctico N°3
Comisión 3
11
Trabajo Práctico N°3
Comisión 3
con la ecuación 4 se calculó el calor ganado o pérfido por la corriente caliente, el signo del
calor nos especifica si absorbió calor (negativo) o si lo emitió (positivo).
Caudal Qf [ W ] Qc [ W ] Qp [ W ] Qf [ W ] Qc [ W ] Qp [ W ]
𝑾
Caudal 𝜼𝒄 % 𝜼𝒇 % 𝜼𝒎 % 𝜼% 𝜟𝑻𝒎𝒍 Uexp(𝒎𝟐𝑲)
𝑾
Caudal 𝜼𝒄 % 𝜼𝒇 % 𝜼𝒎 % 𝜼% 𝜟𝑻𝒎𝒍 Uexp( )
𝒎𝟐 𝑲
12
Trabajo Práctico N°3
Comisión 3
Caudal 𝑾 𝑾 𝑾 𝑾 𝑾 𝑾
Uexp(𝒎𝟐𝑲) UNTU(𝒎𝟐𝑲) UResisT (𝒎𝟐𝑲) Uexp(𝒎𝟐𝑲) UNTU(𝒎𝟐𝑲) UResisT (𝒎𝟐𝑲)
13
Trabajo Práctico N°3
Comisión 3
Discusión de resultados
Con los calores calculados se pudo observar, en primera instancia, que los calores
de la corriente fría y caliente no eran iguales, esto era esperado ya que los calores serían
iguales sólo si el intercambiador estaría perfectamente aislado del exterior, en ese caso solo
habría transferencia de calor desde la corriente caliente a la fría. También se observó que al
aumentar los caudales aumentaba el calor transferido en ambos regímenes de flujo,
notando una mayor transferencia en el flujo contracorriente salvo por el caudal intermedio,
ya que en ese caudal el flujo co-corriente presenta una mayor transferencia de calor.
En cuanto a las eficiencias (Tabla 6 y 7), se notaron diferencias entre las eficiencias
del flujo contracorriente y el co-corriente, siendo mayores las del contracorriente, esto era
de esperarse. En ambos casos, la eficiencia global disminuye al aumentar el caudal de
corriente caliente.
De los valores de los coeficientes globales obtenidos (Tabla 8) se puede ver que
estos aumentan al aumentar el caudal, siendo más grandes en el flujo co-corriente.
También se pudo ver que los obtenidos experimentalmente son muy similares a los
obtenidos por el método NTU, pero existen diferencias significativas comparándolo con el
Método de Resistencias Térmicas.
14
Trabajo Práctico N°3
Comisión 3
Conclusiones y recomendaciones.
Para futuras experiencias, se podría añadir al equipo un buen aislante que permita
disminuir las pérdidas de calor. También resultaría interesante comparar también el efecto
que produce el aumento el número de tubos.
15
Trabajo Práctico N°3
Comisión 3
Referencias Bibliográficas.
[3] Incropera, F.P., Bergman T.L. (2011), “Fundamentos de Transferencia de Calor y Masa”.
John Willey and Sons, 7ma edición, apéndice A.
16
Trabajo Práctico N°3
Comisión 3
Anexo.
Anexo I
Tabla 1. Temperaturas obtenidas en distintos puntos del intercambiador para flujo co-
corriente.
17
Trabajo Práctico N°3
Comisión 3
Anexo II
Cálculo del coeficiente de transferencia de calor (U) para flujo en co-corriente de forma
experimental.
𝑔
Para un caudal másico de flujo caliente de ṁc = 11,6 𝑠 , se tienen los siguientes
datos:
𝑇𝑐,e = 63,2 °C
𝑇𝑐,s = 53,3 °C
𝑔
ṁ𝑓= 17,4 𝑠
𝑇𝑓,e = 15,1 °C
𝑇𝑓,s = 24,7 °C
𝐽
El valor de Cp,c es de 4,185 𝑔𝐾
, el mismo fue buscado a la temperatura media del
fluido caliente, se realizó lo mismo para el fluido frío, obteniéndose un valor de C p,f =4,1816
𝐽
𝑔𝐾
.
El calor de la corriente fría y el calor de la corriente caliente son calculados con las
Ecuaciones 2 y 3, respectivamente:
𝑔 𝐽 𝐽
𝑞c = 11,6 𝑠 . 4,185 𝑔𝐾 . (63,2 – 53,3)𝐾 = 481,02 𝑠
𝑔 𝐽 𝐽
𝑞f = 17,4 𝑠 . 4,181 𝑔𝐾 . (24,7 – 15,1)𝐾 = 700,24 𝑠
𝐽 𝐽 𝐽
𝑞𝑝 = 481,02 𝑠 - 700,2 𝑠 = 13,07 𝑠
18
Trabajo Práctico N°3
Comisión 3
𝐽
481,02 𝑠 𝑊
𝑈 = 2
= 694,94 2
0,018 𝑚 . 37,5 𝐾 𝑚 𝐾
𝐽
700,24 𝑠
𝜂𝑡 = . 100% = 145,57%
𝐽
481,02 𝑠
63,2°𝐶 − 53,3°𝐶
𝜂𝑐 = . 100% = 20,58%
63,2°𝐶 − 15,1°𝐶
24,7 °𝐶 − 15,1°𝐶
𝜂𝑓 = . 100% = 19,95%
63,2°𝐶 − 15,1°𝐶
20,58% + 19,95%
𝜂𝑚 = = 20,27%
2
Cálculo del coeficiente de transferencia de calor por medio del método NTU:
𝑔 𝐽 𝐽
𝐶𝑐,1 = 11,6 . 4,185 = 48,58
𝑠 𝑔𝐾 𝑠𝐾
𝑔 𝐽 𝐽
𝐶𝑓,1 = 17,4 . 4,1816 = 72,24
𝑠 𝑔𝐾 𝑔𝐾
𝐽
48,58 𝑠 𝐾
𝐶= = 0,66
𝐽
72,24 𝑔𝐾
19
Trabajo Práctico N°3
Comisión 3
𝐽 𝐽
𝑄𝑚á𝑥 = 48,58 (63,2𝐾 – 15,1𝐾) = 2337,09
𝑠𝐾 𝑠
𝐽
481,02
ε= 𝑠 = 0,205
𝐽
2337,09
𝑠
𝐽
48,58 𝑠 𝐾 0,252 𝑊
𝑈𝑁𝑇𝑈 = 2
= 663,78 2
0,018𝑚 𝑚 𝐾
Los datos que se necesitaron de viscosidad, densidad, conductividad tanto del agua
como del acero, se extrajeron de la tabla A-9 de (1) y de la tabla A.1 de (3).
𝑘𝑔 𝑤 𝑤
Siendo 𝜇𝑐 = 4,67 𝐸 − 4 𝑃𝑎. 𝑠,𝑝𝑐 = 983,3 , 𝑘𝑐 = 0,654 , 𝑘𝑎𝑐𝑒𝑟𝑜 = 15,1 , 𝑃𝑟𝑐 =
𝑚3 𝑚𝐾 𝑚𝐾
𝑔
2,99, ṁc = 11,6 𝑠
𝑔
mi 11,6 𝑠 1 kg m
𝑣= = . = 0,218
At . 𝑝 π . (0,0083m)2 . 983,3 𝑘𝑔 1000 g s
4 𝑚3
20
Trabajo Práctico N°3
Comisión 3
Para obtener el Uteórico, se necesita tener los coeficientes convectivos de cada fluido.
Éstos se determinaron sabiendo la condición de flujo mediante el Reynolds (Re) el cual fue
calculado por la ecuación 18.
𝑘𝑔 m
𝑝.𝑣.𝐷 983,3 .0,218 .0,0083 m
𝑚3 s
𝑅𝑒 = μ
= 4,67 𝐸−4 𝑃𝑎.𝑠
= 3813,74
𝑓
(3813,74 − 1000)2,99
𝑁𝑢 = 8 = 22,179
1
0,042 2 2
1 + 12,7 ( 8 ) (2,99 − 1) 3
𝑤
22,179.0,654 𝑚𝐾 𝑊
ℎ𝑖 = = 1747,6
0,0083 𝑚 𝑚2 𝐾
Los datos necesarios de la corriente fría, para realizar los cálculos antes
𝑘𝑔 𝑤 𝑤
mencionados son: 𝜇𝑐 = 9,84 𝐸 − 4 𝑃𝑎. 𝑠, 𝑝𝑐 = 997,84 , 𝑘𝑐 = 0,599 , 𝑘𝑎𝑐𝑒𝑟𝑜 = 15,1 ,
𝑚3 𝑚𝐾 𝑚𝐾
𝑔
ṁ𝑓= 17,4 𝑠
𝑚
Se calculó la velocidad y se llegó a un valor de 0,414 𝑠
, el cual permitió estimar el
Reynolds, este tuvo un valor de 1049,79 lo que asegura que es flujo laminar. Sabiendo esto,
se buscó en la bibliografía (1) el 𝑁𝑢 = 3,66
𝑤
3,66 . 0,599 𝑚𝐾 𝑤
ℎ𝑜 = = 877,58 2
0,0025𝑚 𝑚 𝐾
21
Trabajo Práctico N°3
Comisión 3
9,5
1 𝐿𝑛(8,3)
𝑅𝑡𝑜𝑡𝑎𝑙 = + 𝑤
𝑊
1747,6 2 . π. 0,0083m. 0,66m 2. 𝜋. 15,1 𝑚𝐾 . 0,66𝑚
𝑚 𝐾
1 𝐾
+ 𝑤 = 623,09
877,58 2 . π. 0,0095m. 0,66m 𝑊
𝑚 𝐾
1 𝑤
𝑈𝑡𝑒ó𝑟𝑖𝑐𝑜 = = 623,09 2
𝐾 𝑚 𝐾
623,09 𝑊 . π. 0,0083m. 0,66m
22