UNIVERSIDAD AUTÓNOMA DEL ESTADO
DE MORELOS
FACULTAD DE CIENCIAS QUIMICAS E INGENIERÍA
TRANSFERENCIA DE CALOR
ACTIVIDAD #4
“CONVECCIÓN FORZADA INTERNA”
EQUIPO 2
INTEGRANTES:
FLORES SALAZAR HERNÁN
ORTIZ CASTILLO BRANDON JAFETH
SÁNCHEZ SÁNCHEZ HERIBERTO
VAZQUEZ TORICES LUIS ANGEL
FECHA: 09/06/2020
INDICE
Página
1. Resumen convección forzada……………………………… 1-5
2. Ejercicios
2.1 Ejercicios régimen laminar………………………....... 6-10
2.2 Ejercicios régimen en transición………………..…… 11-18
2.3 Ejercicios régimen turbulento……………………..…. 19-26
3. Bibliografías…………………………………………………... 27
CONVECCIÓN FORZADA INTERNA.
En las aplicaciones de calentamiento y enfriamiento es común el uso de un flujo de
un fluido a través de tubos o ductos el cual es forzado a desplazarse por medio de
un ventilador o bomba por un tramo de tubo o ducto lo suficientemente largo para
realizar la transferencia de calor.
Temperatura de bulto:
Es la temperatura promedio a la entrada y salida de una tubería.
𝑇𝑒𝑛𝑡𝑟𝑎𝑑𝑎 + 𝑇𝑆𝑎𝑙𝑖𝑑𝑎
𝑇𝑏 = (𝑬𝒄. 𝟏)
2
Temperatura de película:
Es la temperatura promedio entre la temperatura de la corriente del fluido y
la temperatura de la pared del tubo.
𝑇𝑓𝑙𝑢𝑖𝑑𝑜 + 𝑇𝑝𝑎𝑟𝑒𝑑
𝑇𝑏 = (𝑬𝒄. 𝟐)
2
Velocidad promedio:
Es la velocidad del fluido cuando un fluido pasa a través del área transversal
del ducto.
𝑚̇ = 𝜌𝑣𝐴 (𝑬𝒄. 𝟑) 𝑉̇ = 𝑣𝐴 (𝑬𝒄. 𝟒)
𝑚̇ 𝑉̇
𝑣= (𝑬𝒄. 𝟓) 𝑣= (𝑬𝒄. 𝟔)
𝐴𝜌 𝐴
El flujo en un tubo puede ser laminar o turbulento dependiendo de las condiciones
del mismo el flujo es laminar a velocidades bajas, pero se vuelve turbulento
conforme se incrementa la velocidad. La transición fluctúa entre laminar y
turbulento.
El régimen de flujo se determina a partir del número adimensional de Reynolds.
1
Para un fluido por el interior de un tubo de sección circular el número de Reynolds
se cuenta con la siguiente expresión:
𝜌𝑣𝑑𝑖
𝑅𝑒𝑖 = (𝑬𝒄. 𝟕)
𝜇
𝑅𝑒 ≤ 2300 𝑅𝑒𝑔𝑖𝑚𝑒𝑛 𝑙𝑎𝑚𝑖𝑛𝑎𝑟
𝑅𝑒 < 2300 < 10000 𝑅𝑒𝑔𝑖𝑚𝑒𝑛 𝑒𝑛 𝑡𝑟𝑎𝑛𝑠𝑖𝑐𝑖𝑜𝑛.
𝑅𝑒 > 10000 𝑅𝑒𝑔𝑖𝑚𝑒𝑛 𝑡𝑢𝑟𝑏𝑢𝑙𝑒𝑛𝑡𝑜. Figura.1 Vista frontal
de una tuberia.
ℎ 𝑖 𝑑𝑖
𝑁𝑢𝑖 = (𝑬𝒄. 𝟖)
𝑘
El número de Reynolds y el Nusselt se basan en el
diámetro hidráulico.
4𝐴𝑐
𝐷ℎ = ( 𝑬𝒄. 𝟗)
𝑃
𝜌𝑣𝑑ℎ ℎ𝐷ℎ
𝑅𝑒 = (𝑬𝒄. 𝟏𝟎) 𝑁𝑢 = 𝑬𝒄. (𝟏𝟏)
𝜇 𝑘
Figura 2. Diámetro hidráulico
para distintas geometrías.
Para un doble tubo un fluido fluye por el tubo exterior mientras que otro fluye por la
sección anular.
𝜌𝑣𝐷ℎ
𝑅𝑒𝑜 = (𝑬𝒄. 𝟏𝟐)
𝜇
ℎ𝑜 𝐷𝐻
𝑁𝑢𝑜 = (𝑬𝒄. 𝟏𝟑)
𝑘
𝐷ℎ = 𝐷𝑜 − 𝐷𝑖 (𝑬𝒄. 𝟏𝟒)
Figura 3. Vista frontal Figura 4. Vista frontal
4 𝑚̇
de un doble tubo 𝑅𝑒𝑜 = (𝑬𝒄. 𝟏𝟓) de un tubo exterior.
𝜋(𝐷𝑜 + 𝐷𝑖 )𝜇
2
Se utilizan ecuaciones específicas para estimar los coeficientes de transferencia
de calor para cada caso de estudio.
Régimen laminar.
Para un flujo laminar completamente desarrollado en el interior de un tubo
circular con la temperatura superficial constante (𝑇𝑠 = 𝑐𝑡𝑒):
𝑁𝑢𝑖 = 3.66 𝑪𝒐𝒏𝒅𝒊𝒄𝒊𝒐𝒏 𝟏.
Para el flujo laminar completamente desarrollado en un tubo circular sujeto a
flujo de calor constante en la superficie (𝑄̇ = 𝑐𝑡𝑒):
𝑁𝑢𝑖 = 4.36 𝑪𝒐𝒏𝒅𝒊𝒄𝒊𝒐𝒏 𝟐.
Para el flujo laminar en tubos no circulares obtenemos la siguiente tabla
donde se presenta el número de Nusselt para el flujo laminar:
Tabla 1. Numero de Nusselt y factor de fricción para flujo laminar.
Configuración Numero de Nusselt
𝑎 Factor de
geométrica del 𝑜 𝜃° 𝑇𝑠 = 𝑐𝑜𝑛𝑠t 𝑞̇ 𝑠 = 𝑐𝑜𝑛𝑠𝑡
tubo 𝑏 fricción
3.66
- 4.36 64/Re
1 2.98 3.61 56.92/Re
2 3.39 4.12 62.2/Re
3 53.96 4.79 68.36/Re
4 4.44 5.33 72.92/Re
6 5.14 6.05 78.8/Re
8
5.6 6.49 82.32/Re
∞ 7.54 8.24 96/Re
1 3.66 4.36 64/Re
2 3.74 4.56 67.28/Re
4 3.79 4.88 72.96/Re
8 3.72 5.09 76.6/Re
16 3.65 5.18 78.16/Re
10° 1.61 2.45 50.8/Re
30° 2.26 2.91 52.28/Re
60° 2.47 3.11 53.32/Re
90° 2.34 2.98 52.6/Re
120° 2 2.68 50.96/Re
3
Ecuaciones de transferencia de calor.
𝑄̇ = (𝑚̇𝑐𝑝)𝐻 (𝑇𝑠𝑎𝑙𝑖𝑑𝑎 − 𝑇𝑒𝑛𝑡𝑟𝑎𝑑𝑎 ) ( 𝑬𝒄. 𝟏𝟔)
𝑄̇ = (𝑚̇𝑐𝑝)𝐻 (𝑇𝑒𝑛𝑡𝑟𝑎𝑑𝑎 − 𝑇𝑠𝑎𝑙𝑖𝑑𝑎 ) (𝑬𝒄. 𝟏𝟕)
𝑄̇ = ℎ𝐴𝑠 (𝑇𝑝𝑎𝑟𝑒𝑑 − 𝑇𝑏𝑢𝑙𝑡𝑜 ) ( 𝑬𝒄. 𝟏𝟖)
𝑄̇ = ℎ𝐴𝑠 (𝑇𝑏𝑢𝑙𝑡𝑜 − 𝑇𝑝𝑟𝑒𝑑 ) ( 𝑬𝒄. 𝟏𝟗)
𝑄̇𝐶 = 𝑄̇𝐻 (𝑬𝒄. 𝟐𝟎)
Convección forzada interna: Transición-Turbulento.
Las ecuaciones utilizadas para la estimación de los coeficientes convectivos de
transferencia de calor son semi-empíricas.
El error estimado de las ecuaciones es de 5% a 20%.
Es importante considerar:
Balance de energía: Transferencia de calor entre fluidos
Geometría: tipo de ductos, sección interna, anular.
Régimen de flujo: Transición o completamente turbulento.
Régimen de transición
Gnielinsky propone la siguiente correlación semi-empírica, para calcular el
número de Nusselt promedio.
𝑓
(2) (𝑅𝑒 − 1000)𝑃𝑟
𝑁𝑢 = 1 (𝑬𝒄. 𝟐𝟏) 2,300 < 𝑅𝑒 10,000
𝑓 2 2
1 + 12.7 (2) (𝑃𝑟 3 −1) 0.5 < 𝑃𝑟 < 2,000
Nota: Las propiedades se calculan a la temperatura de bulto.
4
Régimen turbulento.
La ecuación de Dittus-Boelter, para flujo completamente desarrollado en
régimen turbulento:
𝑁𝑢 = 0.023𝑅𝑒 0.8 𝑃𝑟 𝑛 (𝑬𝒄. 𝟐𝟐) 𝑅𝑒 ≥ 10,000
0.6 ≤ 𝑃𝑟 ≤ 160
Nota: Las propiedades se calculan a la temperatura de bulto.
Petukov propone la siguiente ecuación:
𝑓 𝑓 = [(1.58 ∗ ln 𝑅𝑒) − 3.28)−2 ] ( 𝑬𝒄. 𝟐𝟒)
(2) 𝑅𝑒 ∗ 𝑃𝑟
𝑁𝑢 = 1 (𝑬𝒄. 𝟐𝟑)
𝑓 2
1.07 + 12.7 (2) (𝑃𝑟 2/3 −1)
10,000 < 𝑅𝑒 < 1𝑥10−6
0.5 < 𝑃𝑟 < 2,000
5
EJERCICIOS REGIMEN LAMINAR
Primer Ejercicio. Un flujo másico de 1.5 kg/min de agua líquida entra a un
intercambiador de calor tubular de 2.54 cm de diámetro interno a 20 °C; mientras que
abandona el intercambiador a 40 °C.
Figura 5.
Intercambiador de
calor tubular.
1. Calcular la Temperatura de bulto con la Ec1.
20℃ + 40℃
𝑇𝑏𝑢𝑙𝑡𝑜 = = 30 °𝐶
2
2. Las propiedades del agua a la temperatura de bulto las obtenemos de la Tabla A9.
𝑘𝑔
𝜌 = 996
𝑚3
𝑘𝑔
𝜇 = 0.798𝑥10−3
𝑚∙𝑠
𝐽
𝑐𝑝 = 4178
𝑘𝑔°𝐶
𝑊
𝑘 = 0.615
𝑚°𝐶
𝑃𝑟 = 5.42
3. El flujo volumétrico del agua que ingresa al intercambiador con la Ec.4
𝑉̇ = 𝑣𝐴 (𝑬𝒄. 𝟒)
𝑚3
𝑉̇ = 2.51𝑥10−5
𝑠
4. Calcular el Número de Reynolds a partir de la Ec.7
𝜌𝑣𝑑𝑖
𝑅𝑒𝑖 = … … … (𝐸𝑐. 7)
𝜇
6
𝑘𝑔 𝑚
(996 ) (0.04953 𝑠 ) (0.0254𝑚)
𝑚3
𝑅𝑒 = = 𝟏, 𝟓𝟕𝟎. 𝟐𝟏𝟐𝟕 ∴ 𝑭𝒍𝒖𝒋𝒐 𝑳𝒂𝒎𝒊𝒏𝒂𝒓
−3 𝑘𝑔
0.798𝑥10 𝑚 ∙ 𝑠
El fluido se encuentra en el régimen laminar debido a que:
𝑹𝒆 ≤ 𝟐, 𝟑𝟎𝟎
5. Obtener el Coeficiente convectivo de transferencia de calor con la Ec.8.
ℎ𝑖 𝑑𝑖
𝑁𝑢𝑖 = … … … (𝐸𝑐. 8)
𝑘
El número de Nusselt se obtiene a partir de la tabla 1. Sabiendo que 𝑇𝑠 = 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒
𝑁𝑢𝑖 = 3.66
Despejando y sustituyendo la Ec.8:
ℎ𝑖 𝑑𝑖
𝑁𝑢𝑖 =
𝑘
𝑊
𝑁𝑢𝑖 ∙ 𝑘 (3.66) (0.615 𝑚℃)
ℎ𝑖 = =
𝑑𝑖 0.0254𝑚
𝑾
𝒉𝒊 = 𝟖𝟖. 𝟔𝟏𝟖𝟏
𝒎𝟐 °𝑪
7
Segundo Ejercicio. Un flujo másico de 1.5 kg/min de agua líquida entra a un
intercambiador de calor de tipo rectangular a 40 °C; mientras que abandona el
intercambiador a 60 °C. Considera que la condición de transferencia de calor
se realiza de forma constante.
Figura 6.
Intercambiador de calor
de tipo rectangular
Para un ducto rectangular:
Ver Tabla 1 𝑞̇ 𝑠 = Flujo de calor por unidad de área constante
8
Procedimiento de cálculo:
1. Calcular la Temperatura de bulto con la Ec.1
40℃ + 60℃
𝑇𝑏𝑢𝑙𝑡𝑜 = = 𝟓𝟎 °𝑪
2
2. Propiedades del agua a la temperatura de bulto obtenidas de la tabla (A9)
𝑘𝑔
𝜌 = 988.1
𝑚3
𝑘𝑔
𝜇 = 0.547𝑥10−3
𝑚∙𝑠
𝐽
𝑐𝑝 = 4181
𝑘𝑔°𝐶
𝑊
𝑘 = 0.644
𝑚°𝐶
𝑃𝑟 = 3.55
3. Flujo volumétrico del agua que ingresa al intercambiador con la Ec.4
𝑚3
𝑉̇ = 2.5301𝑥10−5
𝑠
4. Calcular la Velocidad con la Ec.6 Sabiendo que: El área de sección transversal
del ducto A=a*b
𝑉̇
𝑣= (𝑬𝒄 𝟔)
𝐴
𝑚3
(2.5301𝑥10−5 𝑠 ) 𝑚
𝑣= = 0.0158
(0.08𝑚 𝑥 0.02 𝑚) 𝑠
9
5. Obtener el diámetro hidráulico con la Figura 2. Para la geometría del ducto.
2𝑎𝑏 (2)(0.08𝑚)(0.02𝑚)
𝐷ℎ = = = 0.032 𝑚
𝑎+𝑏 0.08𝑚 + 0.02𝑚
6. Ahora podemos calcular el Número de Reynolds con la Ec.10
𝜌𝑣𝐷ℎ
𝑅𝑒 = (𝑬𝒄. 𝟏𝟎)
𝜇
𝑘𝑔 𝑚
(988.1 ) (0.0158 𝑠 ) (0.032𝑚)
𝑚3
𝑅𝑒 = = 913.3151 ∴ 𝐹𝑙𝑢𝑗𝑜 𝐿𝑎𝑚𝑖𝑛𝑎𝑟
𝑘𝑔
0.547𝑥10−3 𝑚 ∙ 𝑠
7. Justifica por qué razón el régimen de flujo del fluido es laminar
El fluido se encuentra en el régimen laminar debido a que: 𝑅𝑒 ≤ 2,300
8. Coeficiente convectivo de transferencia de calor con la Ec.11
𝑁𝑢𝐷ℎ = 5.33 𝑜𝑏𝑡𝑒𝑛𝑖𝑑𝑜 𝑑𝑒 𝑙𝑎 𝑡𝑎𝑏𝑙𝑎 1
ℎ𝑖 𝐷ℎ
𝑁𝑢𝐷ℎ = (𝑬𝒄. 𝟏𝟏)
𝑘
𝑊
(𝑁𝑢𝐷ℎ )(𝑘) (5.33) (0.644 𝑚℃)
ℎ𝑖 = =
𝐷ℎ 0.032 𝑚
𝑾
𝒉𝒊 = 𝟏𝟎𝟕. 𝟐𝟔𝟔𝟐
𝒎𝟐 °𝑪
10
EJERCICIO RÉGIMEN EN TRANSICIÓN.
Se requiere enfriar un flujo de agua líquida de 25 L/min desde 80 °C hasta 50 °C.
Para el enfriamiento se dispone de 50 L/min de agua líquida fría a 20 °C.
Para el enfriamiento se dispone de un intercambiador de calor de tubos
concéntricos. El tubo interior es de 3 pulgadas (di=77.9 mm, Di= 88.9 mm) y el tubo
exterior es de 4 pulgadas (Do=114.3 mm).
El agua caliente fluye por el tubo interno; mientras que el agua de enfriamiento fluye
por la sección anular. El intercambiador de calor se encuentra aislado para evitar
las pérdidas de calor con el ambiente.
Calcula:
a) Flujo de calor requerido para el proceso de enfriamiento del agua
R = 𝟓𝟏. 𝟐𝟐𝟑𝟕 𝒌𝑾
b) Temperatura de salida del agua fría a la salida del intercambiador
R = 𝟓𝟎 °𝑪
c) Coeficiente convectivo de transferencia de calor para el tubo interno
𝑾
R = 𝟔𝟗𝟖. 𝟔𝟖𝟓𝟗 𝒎𝟐 °𝑪
d) Coeficiente convectivo de transferencia de calor para la sección anular
𝑾
R = 𝟏𝟐𝟓𝟒. 𝟓𝟖𝟐𝟏 𝒎𝟐 °𝑪
Tabla 2. Resumen de resultados
ℎ
Número
𝑚̇ 𝑊
Nusselt [ 2 ]
𝑘𝑔 𝑚 °𝐶
[ ]
𝑠
Agua caliente 0.4078 82.59125 698.6859
Agua fría 0.817 51.1499 1254.5821
11
Procedimiento de cálculo:
1. Cálculo de propiedades de los fluidos
Para el agua caliente:
𝑇𝑒𝑛𝑡𝑟𝑎𝑑𝑎 + 𝑇𝑠𝑎𝑙𝑖𝑑𝑎
𝑇𝑏,𝐻 = (𝑬𝒄. 𝟏)
2
Sustituyendo:
50 °𝐶 + 80 °C
𝑇𝑏,𝐻 = = 65°𝐶
2
Para el agua fría:
𝑇𝑒𝑛𝑡𝑟𝑎𝑑𝑎 + 𝑇𝑠𝑎𝑙𝑖𝑑𝑎
𝑇𝑏,𝐶 = (𝑬𝒄. 𝟏)
2
Sustituyendo:
20 °𝐶 + 50 °C
𝑇𝑏,𝐶 = = 35 °𝐶
2
Tabla 3. Propiedades obtenidas a las temperaturas de bulto.
𝜌 𝑐𝑝 𝜇 𝑘 𝑃𝑟
𝑘𝑔 𝐽 𝑘𝑔 𝑊
[ 3] [ ] [ ] [ ]
𝑚 𝑘𝑔°𝐶 𝑚∗𝑠 𝑚°𝐶
Agua
0.433 x
caliente 980.4 4187 0.659 2.75
10-3
(65 °C)
Agua fría 0.720 x
994.0 4178 0.623 4.83
(35 °C) 10-3
12
2. Para flujos másicos
Para el flujo másico del agua caliente:
𝐿 1 𝑚3 1 𝑚𝑖𝑛 𝑚3 𝑘𝑔 𝑘𝑔
25 ( )( ) = 0.000416 (980.4 3 ) = 0.4078
𝑚𝑖𝑛 1000 𝐿 60 𝑠 𝑠 𝑚 𝑠
Para el flujo másico del agua fría:
𝐿 1 𝑚3 1 𝑚𝑖𝑛 𝑚3 𝑘𝑔 𝑘𝑔
50 ( )( ) = 0.0008333 (994.4 3 ) = 0.8286
𝑚𝑖𝑛 1000 𝐿 60 𝑠 𝑠 𝑚 𝑠
3. Flujo de calor que cede el agua caliente, durante su enfriamiento
De la ecuación:
𝑄̇𝐻 = (𝑚̇𝑐𝑝)𝐻 (𝑇𝑒𝑛𝑡𝑟𝑎𝑑𝑎 − 𝑇𝑠𝑎𝑙𝑖𝑑𝑎 ) ( 𝑬𝒄. 𝟏𝟕)
Sustituyendo:
𝑘𝑔 𝐽
𝑄̇𝐻 = [0.4078 (4187 )] (80°𝐶 − 50°𝐶) = 𝟓𝟏𝟐𝟐𝟑. 𝟕 𝑾 = 𝟓𝟏. 𝟐𝟐𝟑𝟕 𝒌𝑾
𝑠 𝑘𝑔 °𝐶 𝐻
4. Temperatura de salida del agua de enfriamiento
Del balance de energía:
Por la ecuación:
𝑄̇𝐶 = 𝑄̇𝐻 (𝑬𝒄. 𝟐𝟎 )
Tomamos:
𝑄̇𝐶 = (𝑚̇𝑐𝑝)𝐶 (𝑇𝑠𝑎𝑙𝑖𝑑𝑎 − 𝑇𝑒𝑛𝑡𝑟𝑎𝑑𝑎 ) (𝑬𝒄. 𝟏𝟔)
13
Despejando:
𝑄̇𝐻
𝑇𝑠𝑎𝑙𝑖𝑑𝑎 = + 𝑇𝑒𝑛𝑡𝑟𝑎𝑑𝑎
(𝑚̇𝑐𝑝)𝐻
Sustituyendo en:
𝐽
51223.7 𝑠
𝑇𝑠𝑎𝑙𝑖𝑑𝑎 = + 20 °𝐶 = 𝟓𝟎 °𝑪
𝑘𝑔 𝐽
[(0.4078 𝑠 ) (4187 𝑘𝑔 °𝐶 )]
5. Coeficiente convectivo de transferencia de calor para el tubo interno
Primeramente, se determina el régimen de flujo, con el número de Reynolds para
el agua caliente en el interior del tubo interno:
Número de Reynolds
𝜌𝑣𝑑𝑖
𝑅𝑒𝑖 = (𝑬𝒄. 𝟕)
𝜇
Para 𝑣
𝑚̇𝑖 = 𝜌𝑣𝐴𝑖 (𝑬𝒄. 𝟑)
Sustituyendo en Ec 5:
𝑘𝑔
0.4078 𝑠 𝑚
𝑣= = 0.08727
𝑘𝑔 𝜋(0.0779 𝑚)2 𝑠
(980.4 3 ) ( )
𝑚 4
14
Se sustituye en Ec. 7:
𝑘𝑔 𝑚
(980.4 ) (0.08727 𝑠 ) (0.0779 𝑚)
𝑚3
𝑅𝑒𝑖 = = 15392.8075 ∴ 𝑇𝑢𝑟𝑏𝑢𝑙𝑒𝑛𝑡𝑜
𝑘𝑔
0.433 × 10−3 𝑚 ∗ 𝑠
Número de Nusselt para el agua caliente en el tubo interno
Debido a que el régimen es (turbulento o transición), se calcula el número de Nusselt
mediante la Ecuación de Petukov:
𝑓
(2) (𝑅𝑒 − 1000)𝑃𝑟
𝑁𝑢𝑖 = (𝑬𝒄. 𝟐𝟏)
𝑓 1/2
1 + 12.7 (2) (𝑃𝑟 2/3 − 1)
Para 𝑓:
2
𝑓 = [((1.58 ∗ 𝑙𝑛 𝑅𝑒) − 3.28) ] (𝑬𝒄. 𝟐𝟒)
Sustituyendo en Ec. 24:
2
𝑓 = [((1.58 ∗ 𝑙𝑛 (15392.8075)) − 3.28) ] = 0.006998
De nuevo para 𝑁𝑢 y sustituyendo en Ec. 23:
0.006998
( ) (15392.8075)2.75
𝑁𝑢 = 2 = 82.59125
0.006998 1/2 2/3
1.07 + 12.7 ( ) (2.75 − 1)
2
15
A partir del número de Nusselt se obtiene el coeficiente convectivo de transferencia
de calor, en el interior del tubo:
De la ecuación:
ℎ 𝑖 𝑑𝑖
𝑁𝑢𝑖 = (𝑬𝒄. 𝟖)
𝑘
Despejando y sustituyendo en:
𝑊
𝑁𝑢𝑖 ∗ 𝑘 (82.59125) (0.659 𝑚 °𝐶 ) 𝑾
ℎ𝑖 = = = 𝟔𝟗𝟖. 𝟔𝟖𝟓𝟗 𝟐
𝑑𝑖 0.0779 𝑚 𝒎 °𝑪
6. Coeficiente convectivo de transferencia de calor para la sección anular
Primeramente, se determina el régimen de flujo, con el número de Reynolds para el
agua fría, en la sección anular:
4𝑚̇
𝑅𝑒𝑜 = (𝑬𝒄. 𝟏𝟓)
𝜋(𝐷𝑜 + 𝐷𝑖 )𝜇
Sustituyendo en Ec. 15:
𝑘𝑔
4 (0.8286 𝑠 )
𝑅𝑒𝑜 = = 7211.05566 ∴ 𝑇𝑟𝑎𝑛𝑠𝑖𝑐𝑖ó𝑛
𝑘𝑔
𝜋(0.1143 𝑚 + 0.0889 𝑚)0.720 × 10−3 𝑚∗𝑠
16
Número de Nusselt para el agua fría en la sección anular
Debido a que el régimen es (turbulento o transición), se calcula el número de Nusselt
mediante la Ecuación de Petukov:
En la ecuación:
𝑓
(2) (𝑅𝑒 − 1000) ∗ 𝑃𝑟
𝑁𝑢𝑜 = (𝑬𝒄. 𝟐𝟏)
𝑓 1/2
1 + 12.7 (2) (𝑃𝑟 2/3 − 1)
Para 𝑓:
−2
𝑓 = [((1.58 ∗ 𝑙𝑛 𝑅𝑒) − 3.28) ] (𝑬𝒄. 𝟐𝟒)
Sustituyendo:
−2
𝑓 = [((1.58 ∗ 𝑙𝑛 (7211.05566)) − 3.28) ] = 0.008644
De nuevo para 𝑁𝑢𝑜 y sustituyendo en Ec.21:
0.008644
( ) (7211.05566 − 1000) ∗ 4.83
𝑁𝑢𝑜 = 2 = 51.1499
0.008644 1/2 2/3
1 + 12.7 ( ) (4.83 − 1)
2
A partir del número de Nusselt se obtiene el coeficiente convectivo de transferencia
de calor, en la sección anular, considerando el diámetro hidráulico:
De la ecuación:
ℎ𝑜 𝐷𝐻
𝑁𝑢𝑜 = (𝑬𝒄. 𝟏𝟑)
𝑘
17
Para 𝐷ℎ
𝐷ℎ = 𝐷𝑜 − 𝐷𝑖 (𝑬𝒄. 𝟏𝟒)
Sustituyendo en Ec 14:
𝐷ℎ = 0.1143 𝑚 − 0.0889 𝑚 = 0.0254 𝑚
De nuevo para 𝑁𝑢𝑜
ℎ𝑜 𝐷𝐻
𝑁𝑢𝑜 = (𝑬𝒄. 𝟏𝟑)
𝑘
Despejando ℎ𝑜 y sustituyendo
𝑊
𝑁𝑢𝑜 𝑘 (51.1499) (0.623 𝑚 °𝐶 ) 𝑾
ℎ𝑜 = = = 𝟏𝟐𝟓𝟒. 𝟓𝟖𝟐𝟏 𝟐
𝐷𝐻 0.0254 𝑚 𝒎 °𝑪
18
EJERCICIO RÉGIMEN TURBULENTO.
Ejercicio: Se requiere calentar un flujo de agua líquida de 5000 kg/h desde 20 °C
hasta 35 °C. Para el calentamiento se dispone de agua líquida caliente a 140 °C y
se permitirá un descenso de 15 °C.
Para el calentamiento se dispone de un intercambiador de calor de tubos
concéntricos. El tubo interior es de 2 in (di=2.067 in, Di= 2.357 in) y el tubo exterior
es de 3 in (Do=3.068 in).
El agua caliente fluye por el tubo interno; mientras que el agua de enfriamiento fluye
por la sección anular. El intercambiador de calor se encuentra aislado para evitar
las pérdidas de calor con el ambiente.
Calcula:
e) Flujo de calor requerido para el proceso de calentamiento del agua
f) Flujo másico de agua de calentamiento que ingresa al intercambiador de calor
g) Coeficiente convectivo de transferencia de calor para el tubo interno
h) Coeficiente convectivo de transferencia de calor para la sección anular.
Tabla 3. Resumen de resultados
𝑚̇ Número ℎ Número ℎ
𝑘𝑔 Nusselt 𝑊 Nusselt 𝑊
[ ] [ 2 ] [ 2 ]
𝑠 (Petukov) 𝑚 °𝐶 (Dittus-Boelter) 𝑚 °𝐶
(Dittus-Boelter)
(Petukov)
Agua caliente 1.3596 356.8561 4,647.6258 361.4820 4,707.8727
Agua fría 1.3888 115.5766 2,778.0213 93.8858 2,256.6579
19
Procedimiento de cálculo:
1. Cálculo de propiedades de los fluidos
Para el cálculo de la temperatura de bulto en ambos fluidos usaremos la Ec.1.
Para el agua caliente:
𝑻𝒆𝒏𝒕𝒓𝒂𝒅𝒂 + 𝑻𝒔𝒂𝒍𝒊𝒅𝒂 140℃ + 125℃
𝑻𝒃,𝑯 = = = 𝟏𝟑𝟐. 𝟓℃
𝟐 2
Para el agua fría:
𝑻𝒆𝒏𝒕𝒓𝒂𝒅𝒂 + 𝑻𝒔𝒂𝒍𝒊𝒅𝒂 20℃ + 35℃
𝑻𝒃,𝑪 = = = 𝟐𝟕. 𝟓℃
𝟐 2
Interpolando valores para conocer las propiedades de las temperaturas de
bulto del fluido caliente y frío nos quedan dadas como:
Tabla 4. Propiedades obtenidas a la temperatura de bulto.
𝜌 𝑐𝑝 𝜇 𝑘 𝑃𝑟
𝑘𝑔 𝐽 𝑘𝑔 𝑊
[ 3] [ ] [ ] [ ]
𝑚 𝑘𝑔°𝐶 𝑚∗𝑠 𝑚°𝐶
Agua 932.02 4268.75 0.210x10−3 0.68375 1.315
caliente
Agua 996.5 4179 0.845x10−3 0.611 5.78
fría
𝑘𝑔⁄ ℎ⁄
El flujo másico de 5000 ℎ lo convertimos a 𝑠
𝑘𝑔 1ℎ 𝑘𝑔
5000 ( ) = 1.3838
ℎ 3600 𝑠 𝑠
20
2. Calculamos e l flujo de calor que recibe el agua para su calentamiento
con la Ec. 16
𝑸̇𝑪 = (𝒎̇𝒄𝒑)𝑪 (𝑻𝒔𝒂𝒍𝒊𝒅𝒂 − 𝑻𝒆𝒏𝒕𝒓𝒂𝒅𝒂 )
𝑘𝑔 𝐽
𝑸̇𝑪 = (1.3888 ) (4179 ) (35℃ − 20℃) = 𝟖𝟕𝟎𝟓𝟔. 𝟗𝟐𝟖 𝑾 = 𝟖𝟕. 𝟎𝟓𝟔𝟗 𝒌𝑾
𝑠 𝑘ℎ℃
3. Obtener el flujo másico de agua de calentamiento que ingresa al
intercambiador de calor
Del balance de energía sabemos que:
𝑸̇𝑪 = 𝑸̇𝑯 = 𝟖𝟕𝟎𝟓𝟔. 𝟗𝟐𝟖 𝑾
Y utilizando la Ec.2.
𝑸̇𝑯 = (𝒎̇𝒄𝒑)𝑯 (𝑻𝒆𝒏𝒕𝒓𝒂𝒅𝒂 − 𝑻𝒔𝒂𝒍𝒊𝒅𝒂 )
𝑸̇𝑯
𝒎̇𝑯 =
𝑪𝒑𝑯 (𝑻𝒆𝒏𝒕𝒓𝒂𝒅𝒂 − 𝑻𝒔𝒂𝒍𝒊𝒅𝒂)
87056.928 𝑊 𝒌𝒈
𝒎̇𝑯 = = 𝟏. 𝟑𝟓𝟗𝟔
𝐽 𝒔
(4268.75 𝑘𝑔℃) (140℃ − 125℃)
4. Coeficiente convectivo de transferencia de calor para el tubo interno
Primeramente, se determina la velocidad del fluido con la Ec.3, para después
determinar el régimen de flujo, con el número de Reynolds utilizando la Ec. 7 para
el agua caliente en el interior del tubo interno:
21
𝒎̇𝒊 = 𝝆𝒗𝑨𝒊 (𝑬𝒄. 𝟑)
𝒎̇𝒊
𝑫𝒆𝒔𝒑𝒆𝒋𝒂𝒏𝒅𝒐 𝒗 → 𝒗 = (𝑬𝒄. 𝟓)
𝝆𝑨𝒊
𝑘𝑔
1.3596 𝑠
𝒗= = 𝟎. 𝟔𝟕𝟑𝟖 𝒎⁄𝒔
𝑘𝑔 𝜋(0.05250 𝑚)2
(932.02 3 ) ( )
𝑚 4
𝝆𝒗𝒅𝒊
𝑹𝒆𝒊 = (𝑬𝒄. 𝟕)
𝝁
𝑘𝑔 𝑚
(932.02 ) (0.6738 𝑠 ) (0.05250 𝑚)
𝑚3
𝑅𝑒𝑖 = = 𝟏𝟓𝟔𝟗𝟗𝟖. 𝟕𝟔𝟗 ∴ 𝑭𝒍𝒖𝒋𝒐 𝑻𝒖𝒓𝒃𝒖𝒍𝒆𝒏𝒕𝒐
𝐾𝑔
0.210𝑥10−3 ⁄𝑚 ∙ 𝑠
5. Número de Nusselt para el agua caliente en el tubo interno
Debido a que el régimen es turbulento, se calcula el número de Nusselt mediante la
Ec.23 de Petukov:
𝒇
( ) 𝑹𝒆 ∗ 𝑷𝒓
𝟐
𝑵𝒖𝒊 = (𝑬𝒄. 𝟐𝟑)
𝒇 𝟏/𝟐
𝟏. 𝟎𝟕 + 𝟏𝟐. 𝟕 (𝟐) (𝑷𝒓𝟐/𝟑 − 𝟏)
Donde 𝒇 esta dada por la Ec. 24.
𝒇 = [(𝟏. 𝟓𝟖 ∗ 𝒍𝒏 𝑹𝒆) − 𝟑. 𝟐𝟖)−𝟐 ] = [(1.58)(𝐼𝑛(156998.769)) − 3.28]−2
𝒇 = 𝟎. 𝟎𝟎𝟒𝟎𝟗𝟔 ≈ 𝟎. 𝟎𝟎𝟒𝟏
22
𝟎. 𝟎𝟎𝟒𝟏
( ( )( )
𝑵𝒖𝒊 = 𝟐 ) 𝟏𝟓𝟔𝟗𝟗𝟖. 𝟕𝟔𝟗 𝟏. 𝟑𝟏𝟓
𝟏
𝟎. 𝟎𝟎𝟒𝟏 ⁄𝟐 𝟐
𝟏. 𝟎𝟕 + 𝟏𝟐. 𝟕 + ( ) ((𝟏. 𝟑𝟏𝟓) ⁄𝟑 − 𝟏)
𝟐
𝑵𝒖𝒊 = 𝟑𝟓𝟔. 𝟖𝟓𝟔𝟏
A partir del número de Nusselt se obtiene el coeficiente convectivo de transferencia
de calor, en el interior del tubo con la Ec 8:
𝒉𝒊 𝒅𝒊
𝑵𝒖𝒊 = (𝑬𝒄. 𝟖)
𝒌
(𝑵𝒖𝒊 )(𝒌)
𝑫𝒆𝒔𝒑𝒆𝒋𝒂𝒏𝒅𝒐 𝒉𝒊 → 𝒉𝒊 =
𝒅𝒊
𝑊
(356.8561) (0.68375 )
ℎ𝑖 = 𝑚℃
0.05250 𝑚
𝑾
𝒉𝒊 = 𝟒𝟔𝟒𝟕. 𝟔𝟐𝟓𝟖
𝒎𝟐 ℃
Aplicando la ecuación de Dittus-Boelter:
𝑵𝒖 = 𝟎. 𝟎𝟐𝟑𝑹𝒆𝟎.𝟖 𝑷𝒓𝒏 (𝑬𝒄. 𝟐𝟐)
En donde n= 1/3
1
𝑵𝒖 = (0.023)(156998.769 )0.8 (1.315)3 = 𝟑𝟔𝟏. 𝟒𝟖𝟐𝟎
23
Despejando y sustituyendo nuevamente la Ec.8. El coeficiente convectivo de
transferencia de calor para el tubo interno será:
𝑊
(𝑵𝒖𝒊 )(𝒌) (361.4820) (0.68375 𝑚℃) 𝑾
𝒉𝒊 = = = 𝟒𝟕𝟎𝟕. 𝟖𝟕𝟐𝟕 𝟐
𝒅𝒊 0.05250 𝑚 𝒎 ℃
6. Conocer el coeficiente convectivo de transferencia de calor para la
sección anular
Primeramente, se determina el régimen de flujo, con el número de Reynolds para el
agua fría, en la sección anular:
𝟒𝒎̇
𝑹𝒆𝒐 = (𝑬𝒄. 𝟏𝟓)
𝝅(𝑫𝒐 + 𝑫𝒊 )𝝁
𝑘𝑔
(4) (1.3596
𝑠)
𝑅𝑒𝑜 =
𝑘𝑔
𝜋(0.07792 𝑚 + 0.05250 𝑚) (0.845𝑥10−3 𝑚 ∙ 𝑠)
𝑹𝒆𝒐 = 𝟏𝟓, 𝟕𝟎𝟕. 𝟗𝟖𝟏𝟎𝟗 ∴ 𝑭𝒍𝒖𝒋𝒐 𝑻𝒖𝒓𝒃𝒖𝒍𝒆𝒏𝒕𝒐
Aplicando el número de Nusselt para el agua fría en la sección anular
Debido a que el régimen es turbulento, se calcula el número de Nusselt mediante la
Ecuación de Petukov:
𝒇
(𝟐) 𝑹𝒆 ∗ 𝑷𝒓
𝑵𝒖𝒐 = (𝑬𝒄. 𝟐𝟑)
𝒇 𝟏/𝟐
𝟏. 𝟎𝟕 + 𝟏𝟐. 𝟕 (𝟐) (𝑷𝒓𝟐/𝟑 − 𝟏)
24
Donde 𝒇 esta dada por la Ec. 24.
−𝟐
𝒇 = [(𝟏. 𝟓𝟖 ∗ 𝒍𝒏 𝑹𝒆) − 𝟑. 𝟐𝟖)−𝟐 ] = [((𝟏. 𝟓𝟖)(𝑰𝒏(𝟏𝟓, 𝟕𝟎𝟕. 𝟗𝟖𝟏𝟎𝟗)) − 𝟑. 𝟐𝟖) ]
= 𝟎. 𝟎𝟎𝟔𝟗
Sustituyendo la (Ec.23):
0.0069
(
𝑁𝑢𝑜 = 2 ) (15,707.98109)(5.78) = 𝟏𝟏𝟓. 𝟓𝟕𝟔𝟔
1
0.0069 2 2
1.07 + 12.7 ( ((5.78)3
2 ) − 1)
A partir del número de Nusselt se obtiene el coeficiente convectivo de transferencia
de calor, en la sección anular Ec.13, considerando el diámetro hidráulico Ec. 14:
𝒉𝒐 𝑫𝑯
𝑵𝒖𝒐 = (𝑬𝒄. 𝟏𝟑)
𝒌
𝑫𝒉 = 𝑫𝒐 − 𝑫𝒊 = 0.07792 𝑚 − 0.05250 𝑚 = 0.02542 𝑚
Despejando ℎ𝑜
ℎ𝑜 𝐷𝐻 (𝑁𝑢𝑜 )(𝑘)
𝑁𝑢𝑜 = → ℎ𝑜 =
𝑘 𝐷𝐻
Entonces el coeficiente convectivo de transferencia de calor para la sección anular
será:
𝑊
(115.5766) (0.611
𝒉𝒐 = 𝑚℃) = 𝟐𝟕𝟕𝟖. 𝟎𝟐𝟏𝟑 𝑾
0.02542 𝑚 𝒎𝟐 ℃
25
Aplicando la ecuación de Dittus-Boelter:
𝑵𝒖𝒐 = 𝟎. 𝟎𝟐𝟑𝑹𝒆𝟎.𝟖 𝑷𝒓𝒏 (𝐄𝐜. 𝟐𝟐)
𝒏 = 𝟏⁄𝟑
𝟏
𝑵𝒖𝒐 = (𝟎. 𝟎𝟐𝟑)(𝟏𝟓𝟕𝟎𝟕. 𝟗𝟖𝟏𝟎𝟗)𝟎.𝟖 (𝟓. 𝟕𝟖)𝟑
𝑵𝒖𝒐 = 𝟗𝟑. 𝟖𝟖𝟓𝟖
𝑾
(𝟗𝟓. 𝟖𝟖𝟓𝟖) (𝟎. 𝟔𝟏𝟏
𝒉𝒐 = 𝒎℃) = 𝟐𝟐𝟓𝟔. 𝟔𝟓𝟕𝟗 𝑾
𝟎. 𝟎𝟐𝟓𝟒𝟐 𝒎 𝒎𝟐 ℃
26
BIBLIOGRAFÍA.
Çengel, Y. A., Cengel, Y. A., & Castellanos, J. H. P. (2007). Transferencia
de calor y masa (4.a ed.). CAPITULO 8 – CONVECCIÓN INTERNA
FORZADA New York, Estados Unidos: McGraw-Hill Educación.
27