Universidad Técnica Federico Santa María
Departamento de Ingeniería Química y Ambiental
IQA-222 Transferencia de Calor
Tarea 3
Nombre: Sebastián Ignacio Oñate Inostroza Rol: 201751008-0
Diseño
Para el estudio del fenómeno, se han usado los siguientes valores:
TUBOS
k [W/m°C] 50
espesor [mm] 2,1
METANOL carcaza
flujo [kg/h] 75000
T1 °C 95
T2 °C 45
cp [kJ//kg K] 2,85
Q [kJ/h] -10687500
densidad [kg/m3] 791,8
viscosidad [Pa*s] 0,00059
Rf [m2*°C/W] 0,0002
AGUA DE MAR tubos
flujo [kg/h] 344758,0645
T1 °C 14
T2 °C 45
cp [kJ//kg K] 1
Q [kJ/h] 10687500
densidad [kg/m3] 1000
viscosidad [Pa*s] 0,001
Rf [m2*°C/W] 0,0002
Universidad Técnica Federico Santa María
Departamento de Ingeniería Química y Ambiental
IQA-222 Transferencia de Calor
La temperatura de entrada del agua de mar se ha usado como 14 °C, correspondiente a la
temperatura del mar en la V región de Valparaíso en el mes de septiembre.
Se han usado las propiedades físicas del agua para modelar el agua de mar. Además, se estima
que la temperatura máxima que puede alcanzar dicha agua es de 45 °C app, para evitar
precipitación de sales iónicas presentes en la solución
Se ha calculado el valor de Q mediante un balance de energía del metanol de la siguiente
forma:
𝑘𝑔 𝑘𝐽 𝑘𝐽
𝑄̇ = 𝑚̇𝑐𝑝 ΔT = 75000 [ ] ∗ 2,85 [ ] ∗ (45 − 95)[𝐾] = −10687500 [ ]
ℎ 𝑘𝑔 ∗ 𝐾 ℎ
A partir del cual se puede obtener el flujo de agua de mar necesario para enfriar el metanol,
sabiendo que 𝑄𝑚𝑒𝑡𝑎𝑛𝑜𝑙 = −𝑄𝑎𝑔𝑢𝑎
−𝑄̇ 10687500 𝑘𝑔
= 𝑚̇ = = 82477,7 [ ] 𝑑𝑒 𝑎𝑔𝑢𝑎
𝐶𝑝 𝑎𝑔𝑢𝑎 ∗ Δ𝑇 4,18 ∗ (45 − 14) ℎ
Posteriormente se calcula la LMTD:
Δ𝑇2 −Δ𝑇1
Δ𝑇𝑚 = Δ𝑇2 = 39,7 [s]
ln( )
Δ𝑇1
Universidad Técnica Federico Santa María
Departamento de Ingeniería Química y Ambiental
IQA-222 Transferencia de Calor
Valor que se corrige para arreglos de más de 1 paso, conociendo los valores de P y R, para
ingresar al gráfico que se muestra a continuación:
Obteniendo valores de P=0,38 y R=1,6, para finalmente obtener un F de 0,8 app.
El valor de LMTD*F=31,8[s]
Posteriormente se averigua el valor de Ao, a partir del Q, LMTD*F, y asumiendo un valor
típico de Uo.
𝑄̇
𝐴𝑜 =
𝐿𝑀𝑇𝐷 ∗ 𝐹 ∗ 𝑈𝑜
Una vez conocido Ao, se eligen dimensiones para los tubos del intercambiador de calor. Dado
que no existen limitaciones de espacio se ha elegido un largo de 7,32 [m], y un diámetro de
50 [mm], puesto que el agua de mar es un fluido muy sucio. El espesor seleccionado es de 2,1
[mm]
Luego el número de tubos necesarios se calcula:
𝐴𝑜
𝑁𝑡 =
𝜋 ∗ 𝐷𝑖 ∗ 𝐿
Se ha elegido un arreglo triangular para la distribución de tubos, por lo que el diámetro de
haz de tubos se calcula:
Universidad Técnica Federico Santa María
Departamento de Ingeniería Química y Ambiental
IQA-222 Transferencia de Calor
1
𝑁𝑡 2,207
𝐷𝑏 = 𝐷𝑜 ∗ ( )
0,249
Se asigna un valor de pitch 𝑝𝑡 = 1,25 ∗ 𝐷𝑜 , y se calcula la claridad como 𝐶𝑙 = 𝑝𝑡 − 𝐷𝑜
̇
𝑚𝑎𝑔𝑢𝑎
Para el cálculo de hi, primero es necesario conocer el valor de 𝐺 = 𝐴
𝑡𝑟𝑎𝑛𝑠 𝑡𝑢𝑏𝑜 ∗𝑁𝑡
𝐷𝑖 𝜇
Luego se obtienen 𝑅𝑒 = 𝐺 ∗ 𝑦 𝑃𝑟 = 𝑐𝑝 ∗ 𝑘
𝜇
1
𝐷𝑖 3 𝜇 0,14
Dado que el flujo es laminar, se obtiene 𝑁𝑢 = 1,86 ∗ (𝑅𝑒 ∗ 𝑃𝑟 ∗ ) ∗ (𝜇 ) , asumiendo el
𝐿 𝑤
último término equivalente a 1.
𝑘
Posteriormente ℎ𝑖 = 𝐷 ∗ 𝑁𝑢
𝑖
En el caso de la carcasa, primero se debe definir 𝑙𝑏 correspondiente a la distancia entre los
deflectores, 𝑙𝑏 = 0,4 𝐷𝑏 𝑦 𝐷𝑠 = 𝐷𝑏 + 𝑐𝑙𝑎𝑟𝑖𝑑𝑎𝑑 𝑐𝑎𝑟𝑐𝑎𝑧𝑎, donde la claridad de la carcaza se
obtiene ingresando al siguiente gráfico interceptando Db y el tipo de cabezal(U):
Universidad Técnica Federico Santa María
Departamento de Ingeniería Química y Ambiental
IQA-222 Transferencia de Calor
𝐷
Se debe averiguar 𝐴𝑠 = (𝑝𝑡 − 𝐷𝑜 ) ∗ 𝑝𝑠 ∗ 𝑙𝑏 , el cuál será utilizado para obtener 𝐺𝑠 = ̇
𝑚𝑚𝑒𝑡𝑎𝑛𝑜𝑙
𝑡 𝐴𝑠
𝐺𝑠
Se debe obtener también 𝑣𝑠 = 𝜌
𝐷 2
0,5∗𝑝𝑡 ∗0,86∗𝑝𝑡 −0,5∗𝜋∗ 𝑜
Y un 𝐷𝑒 = 4 ∗ 4
0,5∗𝜋𝐷𝑜
𝐷𝑒 𝜇
Se obtiene el número 𝑅𝑒 = 𝐺𝑠 ∗ , y 𝑃𝑟 = 𝑐𝑝 ∗ 𝑘 , y un valor de jh, ingresando en el siguiente
𝜇
gráfico con el número de Reynolds y un porcentaje de reducción de los deflectores, en este
caso se ha elegido 25%.
1
𝜇 0,14
Luego se calcula 𝑁𝑢 = 𝑗𝐻 ∗ 𝑅𝑒 ∗ 𝑃𝑟 3 ∗ (𝜇 )
𝑤
𝑘
Para obtener ℎ𝑜 = 𝐷 ∗ 𝑁𝑢
𝑒
Con ambos valores de h, solo quda calcular el valor de Uo
Universidad Técnica Federico Santa María
Departamento de Ingeniería Química y Ambiental
IQA-222 Transferencia de Calor
1
𝑈𝑜 =
𝐷
𝐷𝑜 𝑅𝑓,𝑖 𝐷𝑜 ln ( 𝐷𝑜 )
𝐷𝑜 𝑖 1
+ 𝐷 + + 𝑅𝑓,𝑜 +
ℎ𝑖 𝐷𝑖 𝑖 2𝑘 ℎ𝑜
Para calcular las caídas de presión, se debe obtener el factor jf a partir del número de
Reynolds dentro de los tubos, usando el siguiente gráfico:
𝐿 𝜌𝑣 2 𝜇 −𝑚
Luego Δ𝑃 = 8𝑗𝑓 𝐷 ∗ ∗ (𝜇 ) , donde el último término se ha reducido a 1.
2 𝑤
Resultados
Siguiendo los pasos anteriores, y comenzando la iteración con U=650, se obtienen los
siguientes valores en 3 iteraciones.
Universidad Técnica Federico Santa María
Departamento de Ingeniería Química y Ambiental
IQA-222 Transferencia de Calor
Uo Error
Uo
v [m/s] hi [W/m2°C] vs [m/s] ho [W/m2°C]2 calculado porcent jf [-] Caída de Presión [kPa]
[W/m2°C]
[W/m2°C] ual
650 0,0322 1917,8 4,407 143189,8785 1018,85112 56,74 0,085 0,110422709
1018,851
0,0505 2227,8 5,382 174870,6504 1105,58138 8,51 0,07 0,224550002
121
1105,581
0,0548 2289,3 5,581 5077105,313 1128,16864 2,04 0,07 0,264407056
382
Dimensiones Finales
Di [mm] 47,9
Do [mm] 50
L tubo [m] 7,32
N tubos 254
Arreglo triangular
Diámetro de
haz de tubos 1,15
[m]
pt [mm] 62,5
claridad
carcaza 25
[mm]
De los resultados obtenidos es posible observar que no hay caídas de presiones relevantes, por
lo que no es necesario realizar correcciones en ese sentido, sin embargo, la velocidad del flujo
dentro de los tubos es bajísima, lo que implica que implementar un intercambiador de calor con
esas características es poco viable, pues una bomba para mover el fluido no puede funcionar a
tan bajas revoluciones y con los niveles de precisión requeridos. Se sugiere disminuir el diámetro
de los tubos, con el propósito de aumentar la velocidad del fluido dentro de los tubos.
