Intercambiador de Tubos y Coraza
Intercambiador de Tubos y Coraza
Intercambiador de Tubos y Coraza
Objetivos
Teoría
𝑞 𝑑𝑇
= −𝑘
𝐴 𝑑𝑥
En el diseño de los intercambiadores es necesario considerar los coeficientes
convectivos de transferencia de calor, la capa límite entre el fluido y la pared, así
como la temperatura de pared, lo cual es difícil de lograr y obliga a la
experimentación. Para dichos cálculos se encuentran expresiones aritméticas que
expresan la transferencia de calor lograda por los diferentes tipos de
intercambiadores de calor.
∆𝑇1 + ∆𝑇2
∆𝑇𝑚 =
∆𝑇
ln (∆𝑇2 )
1
𝑞 = 𝑈𝑖 𝐴𝑖 𝐹𝑇 ∆𝑇𝑚𝑙 = 𝑈0 𝐴0 𝐹𝑇 ∆𝑇𝑚𝑙
El fluido con la capacidad térmica mínima puede ser tanto el frío como el caliente,
dependiendo de los flujos másicos y los calores específicos.
Con los valores Tml es posible determinar el área superficial del intercambiador si
además se conoce el coeficiente total de transferencia de calor U. Sin embargo, si
las temperaturas de los fluidos que salen del intercambiador no se conocen y se
va a usar un intercambiador determinado, es necesario utilizar un tedioso
procedimiento por prueba y error. En estos casos se utiliza el método llamado
eficacia del intercambiador de calor (ε), que no implica ninguna de las
temperaturas de salida.
ln (𝐷𝑜 /𝐷𝑖 )
𝑅pared =
2π𝑘𝐿
En donde k es la conductividad térmica del material de la pared y L es la longitud
del tubo. Entonces la resistencia térmica total queda:
1 ln (𝐷𝑜 /𝐷𝑖 ) 1
𝑅 = 𝑅total = 𝑅𝑖 + 𝑅pared + 𝑅𝑜 = + +
ℎ𝑖 𝐴𝑖 2π𝑘𝐿 ℎ𝑜 𝐴𝑜
∆𝑇
𝑄= = 𝑈𝐴𝑠 ∆𝑇 = 𝑈𝑖𝐴𝑖 ∆𝑇 = 𝑈𝑜𝐴𝑜 ∆𝑇
𝑅
En donde U es el coeficiente total de transferencia de calor, cuya unidad es
W/m2 · °C, la cual es idéntica a la unidad del coeficiente de convección común, h.
Cancelando ∆T, la ecuación anterior se convierte en:
1 1 1 1 1
= = =𝑅= + 𝑅pared +
𝑈𝐴𝑠 𝑈𝑖 𝐴𝑖 𝑈𝑜 𝐴𝑜 ℎ𝑖 𝐴𝑖 ℎ𝑜 𝐴𝑜
OPERACIÓN y PROCEDIMIENTO
Sólo se deberán operar los equipos bajo la supervisión directa del profesor o
responsable. Antes de iniciar, cerciorarse de que las condiciones de operación a
las que trabajara, están dentro del rango de las especificaciones del
intercambiador.
Utilizar prendas de protección personal adecuadas, tales como: bata,
guantes y lentes de seguridad, ya que es posible que haya fugas de vapor o
agua caliente que pueden ocasionar lesiones severas a los usuarios.
Es necesario conocer los siguientes datos para proceder con los cálculos:
Datos
𝐷𝑒𝑥𝑡 (𝑖𝑛) 1
𝐷𝑖𝑛𝑡 (𝑖𝑛) 0.75
𝐿𝑡 (𝑖𝑛) 356.75
𝑤 45.2
𝑘𝑎𝑐𝑒𝑟𝑜 ( )
𝑚 ∗ °𝐶
𝑃𝑎𝑡𝑚 (𝑎𝑡𝑚) 0.84
𝐽 4179
𝑐𝑝,𝑎𝑔𝑢𝑎 ( )
𝑘𝑔 ∗ °𝐶
𝐴(𝑚2 ) 0.72307157
Para calcular la hliquido, primeramente, es necesario evaluar la temperatura
de película:
25 + 49
𝑇𝑝𝑒𝑙𝑖𝑐𝑢𝑙𝑎 = = 37°𝐶
2
Dicha temperatura se muestra en la primera columna de la siguiente tabla para
cada una de las corridas:
Liquido
t media Pr 𝒌𝒈 𝒌𝒈 Re 𝒘 𝒘
𝝁( ) 𝝆𝒍í𝒒𝒖𝒊𝒅𝒐 ( 𝟑 ) 𝒌𝒍 ( ) 𝒉𝑳𝒊𝒒𝒖𝒊𝒅𝒐 ( 𝟐 )
𝒎∗𝒔 𝒎 𝒎 ∗ °𝑪 𝒎 ∗ °𝑪
37 4.58 0.0006932 993.24 30082.3097 0.6262 4802.67673
37.5 4.57 0.0006865 995.05 29723.5995 0.627 4761.15022
5
36 4.72 0.0007066 993.62 29071.9875 0.6246 4710.93491
8
40 4.32 0.000653 992.1 22927.1449 0.631 3819.17986
40.5 4.27 0.000683 991.7 21920.0961 0.6318 3677.3407
9
39.5 4.37 0.0006417 992.29 27361.706 0. 4411. 630215916
1
Dónde:
Cálculo de hvapor.
Como para el caso del cálculo anterior es necesario conocer algunas propiedades del
vapor que fluye por el exterior de los tubos, las cuales se muestran en la siguiente tabla:
vapor
𝑻𝒔𝒂𝒕 𝒌𝒈 𝒌𝒈 𝒘 𝑱 𝑱 𝒌𝒈 𝒘
𝝆𝒍 ( 𝟑 ) 𝝆𝒗𝒂𝒑𝒐𝒓 ( 𝟑 ) 𝒌𝒍 ( ) 𝒉𝒇𝒈 ( ) 𝒄𝒑𝒍 ( ) 𝝁𝒍 ( ) 𝒉𝒗𝒂𝒑𝒐𝒓 ( 𝟐 )
𝒎 𝒎 𝒎 ∗ °𝑪 𝒌𝒈 𝒌𝒈 ∗ °𝑪 𝒎∗𝒔 𝒎 ∗ °𝑪
H de
evaporización
103 955.71 0.66635 0.6799 2248900 4220.6 0.0002739 9596.52
103 955.71 0.66635 0.6799 2248900 4220.6 0.0002739 9856.25
103 955.71 0.66635 0.6799 2248900 4220.6 0.0002739 10156.06
105 954.25 0.71205 0.6805 2243500 4223 0.0002685 9839.47
105 954.25 0.71205 0.6805 2243500 4223 0.0002685 9858.92
105 954.25 0.71205 0.6805 2243500 4223 0.0002685 9918.46
Dónde se observa que las condiciones como, densidad de líquido, densidad de vapor,
conductividad térmica de líquido, calor latente de vaporización y viscosidad de líquido
fueron evaluadas a la temperatura de saturación de vapor, una vez conocidas estas
propiedades se procede a calcular el hvapor, mediante la expresión:
1
4
𝑔𝜌𝑙 (𝜌𝑙 − 𝜌𝑣 )𝑘𝑙3 3
ℎ𝑣𝑎𝑝𝑜𝑟 = 0.555 (( ) ∗ (ℎ𝑓𝑔 + 𝑐𝑝𝑙 (𝑇𝑠𝑎𝑡 − 𝑇𝑠 )))
𝜇𝑙 (𝑇𝑠𝑎𝑡 − 𝑇𝑠,𝑣𝑎𝑝𝑜𝑟 )𝐷 8
𝐽 3 𝐽 𝑤
∗ (2248900 + 4220.6 (103°𝐶 − 89.3°𝐶)) = 9596.52
𝑘𝑔 8 𝑘𝑔 ∗ °𝐶 𝑚2
∗ °𝐶
)
Una vez conocidos los valores de hlíquido y hvapor es posible calcular el valor teórico
de U, mediante la siguiente expresión:
1
𝑈=
1 ∆𝑥 1
+ + + 𝑅𝑖𝑛𝑐𝑟𝑢𝑠𝑡𝑎𝑐𝑖𝑜𝑛𝑒𝑠
ℎ𝐿í𝑞𝑢𝑖𝑑𝑜 𝑘𝑎𝑐𝑒𝑟𝑜 ℎ𝑣𝑎𝑝𝑜𝑟
Factores de
incrustaciones
(de ambos
lados)
0.00042194
0.00042447
0.0006586
0.00038755
0.00033465
0.00041499
Los cuales fueron calculados mediante los valores experimentales, que se
mostrará más adelante.
Cálculo de UTeórico:
R P F ∆𝑇𝑙𝑚 𝑄(𝑤) 𝑤 𝑤 𝑤
𝑈( 2 ) ℎ𝐿𝑖𝑞𝑢𝑖𝑑𝑜 ( 2 ) ℎ𝑣𝑎𝑝𝑜𝑟 ( 2 )
𝑚 ∗ °𝐶 𝑚 ∗ °𝐶 𝑚 ∗ °𝐶
0.57083333 0.30769231 1 59.00023201 31292.352 733.505205 1803.20739 3158.90345
0.492 0.32051282 1 59.12283722 31896.2175 746.109602 1764.48466 3586.35094
0.49545455 0.28205128 1 61.38282117 28257.1443 636.649111 1776.33266 3585.25858
0.45 0.35443038 0.98 58.36175733 26210.688 633.786081 1294.61043 2876.91207
0.43103448 0.36708861 0.98 57.85840943 27146.784 662.131933 1294.61043 3003.4962
0.45185185 0.34177215 1 59.09142259 29641.2291 693.729999 1518.2775 3360.12234
𝑇1 − 𝑇2
𝑅=
𝑡2 − 𝑡1
103°𝐶 − 89.3°𝐶
𝑅= = 0.5708333
49°𝐶 − 25°𝐶
Mientras que el parámetro p es:
𝑡2 − 𝑡1 49°𝐶 − 25°𝐶
𝑃= = = 0.30769231
𝑇1 − 𝑡1 103°𝐶 − 25°𝐶
𝑄 = 𝐶𝑝 ∗ 𝑤𝐿 ∗ (𝑇𝑠 − 𝑇𝑒 )𝑝𝑟𝑢𝑒𝑏𝑎
Una vez conocido el calor es posible evaluar hlíquido, hvapor y Upráctica, mediante las
expresiones:
𝑄
ℎ𝑙𝑖𝑞𝑢𝑖𝑑𝑜 =
𝐴 ∗ (𝑇𝑠 − 𝑇𝑒 )𝑝𝑟𝑢𝑒𝑏𝑎
𝑄
ℎ𝑉𝑎𝑝𝑜𝑟 =
𝐴 ∗ (𝑇𝑠 − 𝑇𝑒 )𝑐𝑎𝑙𝑒𝑛𝑡𝑎𝑛𝑡𝑒
𝑄
𝑈𝑝𝑟𝑎𝑐𝑡𝑖𝑐𝑎 =
𝑓 ∗ ∆𝑇𝑙𝑚 ∗ 𝐴
hlíquido
31292.352 w 𝑤
ℎ𝑙𝑖𝑞𝑢𝑖𝑑𝑜 = 2
= 1803.20739 2
0.72307157𝑚 ∗ (49°𝐶 − 25°𝐶) 𝑚 ∗ °𝐶
hvapor
31292.352 w 𝑤
ℎ𝑣𝑎𝑝𝑜𝑟 = 2
= 3158.90345 2
0.72307157𝑚 ∗ (103°𝐶 − 89.3°𝐶) 𝑚 ∗ °𝐶
UPráctico
31292.352 w 𝑤
𝑈𝑝𝑟𝑎𝑐𝑡𝑖𝑐𝑎 = 2
= 733.505205 2
1 ∗ 59.00023201°C ∗ 0.72307157𝑚 𝑚 °𝐶
1 1 1 ∆𝑥
𝑅𝑖𝑛𝑐𝑟𝑢𝑠𝑡𝑎𝑐𝑖𝑜𝑛𝑒𝑠 = − − −
𝑈𝑃𝑟𝑎𝑐𝑡𝑖𝑐𝑎 ℎ𝐿í𝑞𝑢𝑖𝑑𝑜 ℎ𝑣𝑎𝑝𝑜𝑟 𝑘𝑎𝑐𝑒𝑟𝑜
1 1 1
𝑅𝑖𝑛𝑐𝑟𝑢𝑠𝑡𝑎𝑐𝑖𝑜𝑛𝑒𝑠 = 𝑤 − 𝑤 − 𝑤
733.505205 2 1803.20739 2 3158.90345 2
𝑚 °𝐶 𝑚 ∗ °𝐶 𝑚 ∗ °𝐶
(1 − 0.75) ∗ (0.0254𝑚) 𝑚2 ∗ °𝐶
− 𝑤 = 0.00042194
2 ∗ 45.2 𝑚 ∗ °𝐶 𝑤
ℎ𝑙𝑖𝑞𝑢𝑖𝑑𝑜,𝑡𝑒ó𝑟𝑖𝑐𝑜 − ℎ𝑙𝑖𝑞𝑢𝑖𝑑𝑜,𝑃𝑟á𝑐𝑡𝑖𝑐𝑜
% 𝑑𝑒 𝑑𝑒𝑠𝑣𝑖𝑎𝑐𝑖ó𝑛 = ∗ 100
ℎ𝑙𝑖𝑞𝑢𝑖𝑑𝑜,𝑡𝑒ó𝑟𝑖𝑐𝑜
𝑤 𝑤
4802.68 − 1803.20 2
% 𝑑𝑒 𝑑𝑒𝑠𝑣𝑖𝑎𝑐𝑖ó𝑛 = 𝑚2 ∗ °𝐶 𝑚 ∗ °𝐶 ∗ 100 = 62.45%
𝑤
4802.68 2
𝑚 ∗ °𝐶
Tabla con los % de desviación
% desviación
𝒉𝒍𝒊𝒒𝒖𝒊𝒅𝒐
62.45
62.94
62.29
66.10
64.79
65.58
ℎ𝑣𝑎𝑝𝑜𝑟,𝑡𝑒ó𝑟𝑖𝑐𝑜 − ℎ𝑣𝑎𝑝𝑜𝑟,𝑃𝑟á𝑐𝑡𝑖𝑐𝑜
% 𝑑𝑒 𝑑𝑒𝑠𝑣𝑖𝑎𝑐𝑖ó𝑛 = ∗ 100
ℎ𝑣𝑎𝑝𝑜𝑟,𝑡𝑒ó𝑟𝑖𝑐𝑜
𝑤 𝑤
9596.52 − 3158.90 2
% 𝑑𝑒 𝑑𝑒𝑠𝑣𝑖𝑎𝑐𝑖ó𝑛 = 𝑚2 ∗ °𝐶 𝑚 ∗ °𝐶 ∗ 100 = 67.08%
𝑤
9596.52 2
𝑚 ∗ °𝐶
% desviación
𝒉𝒗𝒂𝒑𝒐𝒓
67.08
63.61
64.70
70.76
69.54
66.12
Para U global:
𝑈 𝑡𝑒ó𝑟𝑖𝑐𝑜 − ℎ𝑣𝑎𝑝𝑜𝑟,𝑃𝑟á𝑐𝑡𝑖𝑐𝑜
% 𝑑𝑒 𝑑𝑒𝑠𝑣𝑖𝑎𝑐𝑖ó𝑛 = ∗ 100
𝑈 𝑡𝑒ó𝑟𝑖𝑐𝑜
𝑤 𝑤
1242.85 − 733.51 2
% 𝑑𝑒 𝑑𝑒𝑠𝑣𝑖𝑎𝑐𝑖ó𝑛 = 𝑚2 ∗ °𝐶 𝑚 ∗ °𝐶 ∗ 100 = 40.98%
𝑤
1242.85 2
𝑚 ∗ °𝐶
Tabla con los % de desviación
%desviación
𝑼 𝒕𝒆ó𝒓𝒊𝒄𝒐
40.98
39.85
33.82
47.95
48.47
43.62
𝐽 𝑘𝑔 𝑤
𝑞 = 4179 ∗ 0.312 ∗ (49°𝐶 − 25°𝐶) = 0.7231𝑚2 ∗ (3158.90 2 ) ∗ (103°𝐶 − 89.3°𝐶)
𝑘𝑔 ∗ 𝑘 𝑠 𝑚 ∗ °𝐶
31292.35𝑊 = 31292.35𝑊
Ahora con:
ℎ𝑣𝑎𝑝𝑜𝑟,𝑃𝑟á𝑐𝑡𝑖𝑐𝑜 ∗ ∆𝑇𝑣 = ℎ𝑙í𝑞𝑢𝑖𝑑𝑜,𝑃𝑟á𝑐𝑡𝑖𝑐𝑜 ∗ ∆𝑇𝑙
𝑤 𝑤
(3158.90 2 ) ∗ (103°𝐶 − 89.3°𝐶) = (1803.21 2 ) ∗ (49°𝐶 − 25°𝐶)
𝑚 ∗ °𝐶 𝑚 ∗ °𝐶
43277𝑊 = 43277𝑊
Conclusión
Bibliografía:
[1] (ÇENGEL, Yunus A.; GHAJAR, A. Transferencia de calor y masa. Editorial: Editorial McGraw–
Hill, 2007.)