CAPÍTULO 3: CONVECCIÓN

OBJETIVOS
• Comprender los fundamentos y aplicar las leyes que gobiernan el
fenómeno de transporte de calor por convección (Ley de enfriamiento
de Newton para convección y ecuaciones de balance).
• Establecer una relación entre convección y la Ley de Fourier para
conducción pura.
• Comprender y aplicar métodos para la determinación del coeficiente
de transmisión de calor por convección. Definición de Números de
Nusselt, Reynolds, Prandtl y Grashof.
• Relaciones funcionales entre parámetros adimensionales relevantes
para la determinación de coeficientes de convección h.
• Resolver problemas de convección.
•
Q gp = h gp A p (T g , − T p , s )

LEY DE ENFRIAMIENTO DE NEWTON

(Isaac Newton; 1642-1727)
La ley establece que la tasa de transferencia de
calor por convección entre un sólido y un fluido,
que se encuentran en contacto a través de una
superficie A y a diferentes temperaturas, es
proporcional a la diferencia de temperatura
entre ambos medios.

dQ/dt=h A (Ts-T∞) (W)

Donde:
h=coeficiente de transferencia de calor por
convección (W/m2 K)
ORDENES DE MAGNITUD DE “h”
Rangos de valores del coeficiente de transferencia
de calor por convección

Valores promedios de ℎ
Convección libre 𝑊
ℎ
𝑚2 𝐾
Gases 2-25
Líquidos 50-1.000
Convección forzada 𝑊
ℎ
𝑚2 𝐾
Gases 25-250
Líquidos 50-20.000
 ORDENES DE MAGNITUD DE “h”
Coeficiente de transferencia de calor por convección
Fluido 𝑾 𝑩𝒕𝒖
𝒎𝟐 𝑲 𝒉𝒇𝒕𝟐 °𝑭
Aire, convección libre 6-30 1-5
Vapor o aire 30-300 5-50
sobrecalentado,
convección forzada
Aceite, convección forzada 60-1.800 10-300
Agua, convección forzada 300-18.000 50-3.000
Agua en ebullición 3.000-60.000 500-10.000
Vapor condensándose 6.000-120.000 1.000 -20.000

Estos órdenes de magnitud son referenciales y orientadores y, en

general, no deben ser usados en el diseño térmico de sistemas.
Coeficientes de convección en la NCh 853

El valor de Rsi = 0.12 m2 K/W corresponde a hi = 8.3 W/m2 K

El valor de Rse = 0.05 m2 K/W corresponde a he = 20 W/m2 K

Coeficientes de convección en macerado de madera

Este proceso se puede

modelar, en primera
instancia, considerando
un coeficiente de
superficie igual a

h = 500-600 W/m2 K
 VARIABLES PRINCIPALES QUE INFLUYEN EN LA CONVECCIÓN

• El coeficiente de transferencia de calor por convección depende del

fluido (en el cual se estudia el fenómeno) a través de su conductividad
térmica, y de la velocidad del flujo. Así:

h = h (kf, vf)
• El coeficiente de película h, en general presenta la siguiente relación:

h líquidos > h gases

h CF > h CN
h turbulento > h laminar
PROBLEMA DE CONVECCIÓN EN SISTEMA DE CALEFACCIÓN
PROBLEMA DE CONVECCIÓN EN INGENIERÍA ELECTRÓNICA

Aire Dispositivos
10 𝑚/𝑠 0,005 m
electrónicos
320 K

Turbulador
PROBLEMA DE CONVECCIÓN EN INGENIERÍA TÉRMICA SOLAR

0.05 m

Cubierta de cristal, 35 °C
Placa absorbente, 70 °C

Aislante
0°, 30 ° o 60°
PROBLEMA DE CONVECCIÓN EN INGENIERÍA TÉRMICA SOLAR
PROBLEMA DE CONVECCIÓN EN INGENIERÍA TÉRMICA SOLAR
PROBLEMA DE CLIMATIZACIÓN DE PISCINAS
FUNDAMENTOS FÍSICOS QUE GOBIERNAN LA CONVECCIÓN
Gradiente de la velocidad producto de la viscosidad y
gradiente térmico producto de la conductividad térmica

Flujo de fluido
𝜕𝑇
ቤ
𝜕y 𝑦=0
𝑈∞
𝑦
𝑢(𝑦) 𝑇(𝑦)
Superficie de la placa

𝑞ሶ 𝑇∞ 𝑇𝑠
Flujo de calor
 FUNDAMENTOS FÍSICOS QUE GOBIERNAN LA CONVECCIÓN
Conceptos de capa límite fluodinámica y capa límite térmica
𝑈∞

𝑞ሶ 𝑇∞ 𝑈∞
𝑇∞ 𝛿𝑇 𝑦
𝛿
𝑢(𝑦)
𝑈∞ , 𝑇∞ 𝑇(𝑦)
𝑥

𝑇𝑠 𝛿 =capa delimitadora de velocidad.

𝛿𝑇 =capa delimitadora térmica.

Desarrollo de las capas delimitadoras de velocidad y térmica en un flujo sobre una

superficie con forma arbitraria
 FUNDAMENTOS FÍSICOS QUE GOBIERNAN LA CONVECCIÓN
Así como la ley de enfriamiento de Newton establece que el flujo de calor por
unidad de área es directamente proporcional a la diferencia de temperatura
entre la superficie y el fluido, también se puede establecer, a través de la ley
de Fourier, que esa densidad de flujo es igual al producto de la conductividad
térmica del fluido por el gradiente de T° en la pared.

𝜕𝑇
𝑞ሶ = −𝑘𝑓 ቤ = ℎ (𝑇𝑠 − 𝑇∞ )
𝜕y 𝑦=0
 FUNDAMENTOS FÍSICOS QUE GOBIERNAN LA CONVECCIÓN
Definición del coeficiente h
En el planteamiento anterior, la rigurosidad en la aplicación de la ley
de Fourier, a un fluido en movimiento, se basa en que producto de la viscosidad
del fluido, la velocidad en la pared es nula. Por lo tanto aplica el concepto de
conducción pura que subyace en la ley de Fourier.

La transferencia de calor en la primera capa de fluido adyacente

a ella, es por conducción pura.

Así, la ecuación que define el coeficiente de transferencia de calor en W/m2 K es:

convección es:
𝑑𝑇
−𝑘 𝑦=0
𝑑𝑦
ℎ=
(𝑇𝑠 − 𝑇∞ )
 FUNDAMENTOS FÍSICOS QUE GOBIERNAN LA CONVECCIÓN
Definición de parámetro adimensionales
En base a la ecuación anterior se define el parámetro adimensional
de Nusselt (Wilhelm Nusselt; 1882-1957)
Coeficiente adimensional de
transferencia de calor; relación entre la
Número de Nusselt (𝑁𝑢𝐿 ) = transferencia de calor por convección y
la conducción en una capa de fluido con
espesor L.
ℎത 𝐿 𝑇𝑟𝑎𝑛𝑠𝑓𝑒𝑟𝑒𝑛𝑐𝑖𝑎 𝑑𝑒 𝑐𝑎𝑙𝑜𝑟 𝑝𝑜𝑟 𝑐𝑜𝑛𝑣𝑒𝑐𝑐𝑖ó𝑛
𝑁𝑢𝐿 = =
𝑘𝑓 𝑇𝑟𝑎𝑛𝑠𝑓𝑒𝑟𝑒𝑛𝑐𝑖𝑎 𝑑𝑒 𝑐𝑎𝑙𝑜𝑟 𝑝𝑜𝑟 𝑐𝑜𝑛𝑑𝑢𝑐𝑐𝑖ó𝑛

En la ecuación anterior, la longitud característica L es propia de cada

problema y debe ser determinada con atención especial.
 FUNDAMENTOS FÍSICOS QUE GOBIERNAN LA CONVECCIÓN
Definición de parámetro adimensionales

Nota: No confundir el N° de Nusselt con el N° de Biot

Además, se define el parámetro adimensional de Prandtl,

como propiedad física del fluido (Ludwig Prandtl; 1875-1953)

𝜈 𝑐𝑝 𝜇 Relación entra la difusividad de cantidad

Número de Prandtl (𝑃𝑟) = = de movimiento y la difusividad térmica.
𝛼 𝑘
FUNDAMENTOS FÍSICOS QUE GOBIERNAN LA CONVECCIÓN
Ludwig Prandtl
FUNDAMENTOS FÍSICOS QUE GOBIERNAN LA CONVECCIÓN
Capa límite
 FUNDAMENTOS FÍSICOS QUE GOBIERNAN LA CONVECCIÓN
𝑑𝑇 En la ecuación que define al
−𝑘 𝑦=0
𝑑𝑦 coeficiente h, la velocidad del flujo
h = h (kf, vf) ℎ=
(𝑇𝑠 − 𝑇∞ ) está implícita en el gradiente de
temperatura en la pared ya que el
campo de temperaturas en el
fluido depende de la velocidad.
Continuidad Esta dependencia se da a través de
𝜕𝑢 𝜕𝑣
𝜕𝑥
+ 𝜕𝑦=0 las ecuaciones de balance de
Cantidad de movimiento masa, cantidad de movimiento y
𝜕𝑢 𝜕𝑢 𝜕 2 𝑢 𝜕𝑝 energía, que para un flujo
𝜌 𝑢 +𝑣 =𝜇 2−
𝜕𝑥 𝜕𝑦 𝜕𝑦 𝜕𝑥 incompresible en régimen
Energía estacionario 2D se escriben como:
𝜕𝑇 𝜕𝑇 𝜕2𝑇
𝑢 +𝑣 =∝ 2
𝜕𝑥 𝜕𝑦 𝜕𝑦
 3.1 CONVECCIÓN FORZADA
CASO PLACA PLANA EN RÉGIMEN LAMINAR-SOLUCIÓN ANALÍTICA
PROBLEMA 1
PROBLEMA 2
 CONVECCIÓN FORZADA
CASO PLACA PLANA EN RÉGIMEN LAMINAR-SOLUCIÓN ANALÍTICA
B D
𝑦 𝑦𝑠
𝑇∞

𝑈∞
Espesor de la capa
Espesor de la delimitadora de
capa velocidad 𝛿
u(𝑦) delimitadora de
𝑇(𝑦) temperatura 𝛿𝑡

A 𝑇𝑠 C 𝑥
dx

Volumen de control para el análisis integral de la conservación

integral de la energía
 CONVECCIÓN FORZADA
CASO PLACA PLANA EN RÉGIMEN LAMINAR

Solución analítica de Pohlhausen (1921)

1 1
𝑘 𝑈∞ 𝑥 2 𝜈 3
ℎ𝑥 = 0,332
𝑥 𝜈 𝛼
Donde:
ℎ𝑥 = coeficiente de transferencia convectiva local (en x)
𝑈∞ 𝑥
𝑅𝑒𝑥 = = 𝑁ú𝑚𝑒𝑟𝑜 𝑑𝑒 𝑅𝑒𝑦𝑛𝑜𝑙𝑑𝑠 𝑙𝑜𝑐𝑎𝑙 𝑒𝑛 𝑥
𝜈
𝜈
𝑃𝑟 = = 𝑁𝑢𝑚𝑒𝑟𝑜 𝑑𝑒 𝑃𝑟𝑎𝑛𝑑𝑡𝑙 𝑑𝑒𝑙 𝑓𝑙𝑢𝑖𝑑𝑜
𝛼
 CONVECCIÓN FORZADA
CASO PLACA PLANA EN RÉGIMEN LAMINAR

Así, si
ℎ𝑥 𝑥 1 1
= 0,332 𝑅𝑒𝑋 2 𝑃𝑟 3 ∗
𝑘
O
1 1
𝑁𝑢𝑥 = 0,332 𝑅𝑒𝑋 2 𝑃𝑟 3
Donde
ℎ𝑥 𝑥
𝑁𝑢𝑥 = Número de Nusselt local en x
𝑘
A partir de este resultado se espera que en general
𝑁𝑢𝐿 = 𝑐 𝑅𝑒𝐿𝑚 𝑃𝑟 𝑛
 CONVECCIÓN FORZADA
CASO PLACA PLANA EN RÉGIMEN LAMINAR

O
𝑁𝑢𝐿 = 𝑓(𝑅𝑒𝐿 , 𝑃𝑟)

Que es una forma adimensional de escribir la relación

ℎ = 𝑓(𝑣, 𝑘)

Integrando la ecuación ∗ a lo largo de toda la placa de largo L,

se obtiene ℎത 𝐿
 CONVECCIÓN FORZADA
CASO PLACA PLANA EN RÉGIMEN LAMINAR

𝑥=𝐿
‫=𝑥׬‬0 ℎ𝑥 𝑑𝑥
ℎത 𝐿 =
𝐿
1
1 𝑇∞
𝑥=𝐿 2 3
1 0.332𝑅𝑒𝑥 Pr 𝑘
ℎത 𝐿 = න 𝑑𝑥 𝑦
𝛿𝑇
𝐿 𝑥=0 𝑥 𝑥

1 1 𝑥
𝑥=𝐿 21 𝑇𝑠
0.332 3 𝑈∞ 𝑥 𝐿
ℎത 𝐿 = 𝑘 Pr න 𝑑𝑥
𝐿 𝑥=0 𝜈 𝑥

1 1 𝐿
𝑘 3 𝑈∞ 𝑥 2
ℎത 𝐿 = 0.332 Pr 2 ቮ
𝐿 𝜈
𝑥=0
 CONVECCIÓN FORZADA
CASO PLACA PLANA EN RÉGIMEN LAMINAR
𝑘 1 1
ℎത 𝐿 = 0.664 𝑅𝑒𝐿 2 𝑃𝑟 3
𝐿
O en términos adimensionales
1 1 ℎത 𝐿 𝐿
𝑁𝑢𝐿 = 0.664 𝑅𝑒𝐿 2 𝑃𝑟 3 =
𝑘
Resultando válido si:

𝑅𝑒𝐿 < 𝑅𝑒𝑐 = 5 𝑥 105 , 𝑐 = 𝑐𝑟í𝑡𝑖𝑐𝑜

Las propiedades se evalúan a temperatura media 𝑇ത𝑓 : 𝑇° media de film o película.

𝑇∞ + 𝑇𝑠
𝑇ത𝑓 =
2
 CONVECCIÓN FORZADA
CASO PLACA PLANA EN RÉGIMEN LAMINAR

RESUMEN
ℎ𝑥 𝑥 1 1
𝑁𝑢𝑥 = = 0.332 𝑅𝑒𝑥 2 𝑃𝑟 3
𝑘

𝑅𝑒𝐿 < 5x105

1 1
𝑁𝑢𝐿 = 0.664 𝑅𝑒𝐿 2 𝑃𝑟 3
L < xc

Cuando el número de Reynolds supera el valor crítico de 5x105 , se produce la transición desde el
régimen laminar a un régimen turbulento, con lo cual el perfil de velocidad también se transforma
producto de velocidades transversales al flujo principal.
 CONVECCIÓN FORZADA
CASO PLACA PLANA EN RÉGIMEN TURBULENTO

Ley de Viscosidad de Newton

𝑑𝑢
𝜏=𝜇
𝑑𝑦

Debido a la viscosidad del fluido que genera el perfil de velocidades con velocidad nula en la pared
sólida, se produce un esfuerzo de corte gobernado por la Ley de Viscosidad de Newton.
 CONVECCIÓN FORZADA
CASO PLACA PLANA EN RÉGIMEN TURBULENTO
Transición desde un régimen laminar a un régimen turbulento para un flujo
paralelo sobre una placa plana
 CONVECCIÓN FORZADA
CASO PLACA PLANA EN RÉGIMEN TURBULENTO
Si 𝑅𝑒𝐿 > 𝑅𝑒𝑐 , el número de Nusselt local se puede evaluar con la siguiente expresión:

0.8
ℎ𝑥 𝑥 1 𝑈∞ 𝑥
𝑁𝑢𝑥 = = 0.0288 𝑃𝑟 3
𝑘 𝜈

Al comparar con el resultado obtenido por Pohlhausen para régimen laminar y reflejada en la
expresión siguiente:

ℎ𝑥 𝑥 1 1
𝑁𝑢𝑥 = = 0.332 𝑅𝑒𝑥 2 𝑃𝑟 3
𝑘

Se puede demostrar que en flujo turbulento el coeficiente de transferencia de calor es mayor

que en régimen laminar.
 CONVECCIÓN FORZADA
CASO PLACA PLANA EN RÉGIMEN TURBULENTO
Al integrar la ecuación del coeficiente local de transferencia de calor en régimen turbulento, sobre
una placa de largo L

1 𝐿
ℎത = න ℎ𝑐𝑥 𝑑𝑥
𝐿 0
se obtiene:
ℎത 𝐿 1
0.8
𝑁𝑢𝐿 = = 0.036 𝑃𝑟 3 𝑅𝑒𝐿
𝑘
En la expresión anterior, al integrar
desde x=0 hasta x = L, se ha
despreciado la zona de la capa límite
laminar y por lo tanto se puede
aplicar si L >> xc; la distancia xc se
obtiene de la expresión:
Rexc = 5 x10+5
 CONVECCIÓN FORZADA
CASO PLACA PLANA EN RÉGIMEN TURBULENTO
Al considerar la zona de la capa límite laminar, se obtiene una expresión que considera las zonas
laminar y turbulenta:

1
𝑁𝑢𝐿 = 0.036 𝑃𝑟 3 (𝑅𝑒𝐿0.8 − 23,200)
 CONVECCIÓN FORZADA
CASO FLUJO EN INTERIOR DE UN CILINDRO EN RÉGIMEN LAMINAR
ഥ
𝜌𝑈𝐷
𝑅𝑒𝐷 = < 2300
𝜇
BALANCE DE MOMENTUM,
Mec de Fluidos

BALANCE DE ENERGÍA,
Termodinámica y Transferencia de calor
 CONVECCIÓN FORZADA
CASO FLUJO EN INTERIOR DE UN CILINDRO EN RÉGIMEN LAMINAR
 CONVECCIÓN FORZADA
CASO FLUJO EN INTERIOR DE UN CILINDRO EN RÉGIMEN LAMINAR
SOLUCIÓN ANALÍTICA DE LAS ECUACIONES DE BALANCE.
Cap. 6.2.1 y 6.2.2 Texto de F. Kreith y M Bohn (2001)

ℎ𝑐 𝐷
𝑁𝑢𝐷 = = 3,66 𝑇𝑠 = 𝑐𝑜𝑛𝑠𝑡
𝑘

ℎ𝑐 𝐷
𝑁𝑢𝐷 = = 4,364 𝑞𝑠 ′′ = 𝑐𝑜𝑛𝑠𝑡
𝑘
 CONVECCIÓN FORZADA
DATOS EXPERIMENTALES PARA FLUJO EN TUBERÍA EN RÉGIMEN LAMINAR Y TURBULENTO
 CONVECCIÓN FORZADA
CASO FLUJO EN INTERIOR DE UN CILINDRO EN RÉGIMEN LAMINAR
EFECTO DE LA VARIACIÓN DE LA VISCOSIDAD CON LA
TEMPERATURA

Cuando de trata de flujos muy viscosos, el efecto de la T°

puede ser importante en el cálculo del coeficiente de
transferencia de calor por convección.

En la figura se muestran tres perfiles de velocidades en el

interior de una tubería:

A) perfil de velocidades obtenido en mecánica de

Fluidos en un dominio isotérmico,
B) corresponde a un perfil deformado por efecto de la
T° en un flujo que recibe calor desde el exterior.
C) corresponde a un perfil deformado por efecto de la
VARIACIÓN VISCOSIDAD T° en un flujo que rechaza calor hacia el exterior.
PERFIL DE VELOCIDADES
CINEMÁTICA CON LA T°
EN FLUJOS CON ENFRIAMIENTO
Y CALENTAMIENTO
 CONVECCIÓN FORZADA
CASO FLUJO EN INTERIOR DE DUCTOS EN RÉGIMEN LAMINAR
ECUACIÓN EMPÍRICA DE SIEDER Y TATE (1936) PARA
LÍQUIDOS (Considera efecto de T° sobre viscosidad y h)

ഥ
𝜌𝑈𝐷
0.33 0.14 𝑅𝑒𝐷 = < 2300
𝑅𝑒𝐷𝐻 𝑃𝑟𝐷𝐻 𝜇𝑏 𝜇
𝑁𝑢𝐷𝐻 = 1.86
𝐿 𝜇𝑠
0.4 < Pr < 16.7

Esta ecuación introduce el factor empírico de corrección del perfil de velocidades por efecto
de la T° en un flujo laminar.

Las propiedades, excepto la viscosidad en la superficie, se evalúan a la temperatura media

del volumen del flujo másico.

Se usa para tubos con temperatura superficial constante.

 CONVECCIÓN FORZADA
CASO FLUJO EN INTERIOR DE DUCTOS EN RÉGIMEN LAMINAR

ECUACIÓN EMPÍRICA DE HAUSEN (1983)

(Considera efecto de T° sobre viscosidad y h)

0.0668𝑅𝑒𝐷𝐻 𝑃𝑟 𝐷ൗ𝐿 𝜇𝑏
0.14
ഥ
𝑁𝑢𝐷𝐻 = 3.66 + 𝜌𝑈𝐷
0.66 𝑅𝑒𝐷 = < 2300
1 + 0.045 𝑅𝑒𝐷𝐻 𝑃𝑟 𝐷ൗ𝐿 𝜇𝑠 𝜇

Donde: 100 < 𝑅𝑒𝐷𝐻 𝑃𝑟 𝐷Τ𝐿<1500)

Esta ecuación también introduce el factor empírico de corrección del perfil de

velocidades por efecto de la T° en un flujo laminar.

Las propiedades, excepto la viscosidad en la superficie, se evalúan a la temperatura

media del volumen del flujo.
 CONVECCIÓN FORZADA
FLUJO EN INTERIOR DE DUCTOS Y TUBERÍAS CIRCULARES
RÉGIMEN TURBULENTO
ECUACIÓN DE DITTUS-BOELTER (1930) PARA FLUIDOS EN GENERAL

ℎത 𝐷
𝑁𝑢𝐷 = = 0.023 𝑅𝑒𝐷0.8 𝑃𝑟 𝑛 0.5 < 𝑃𝑟 < 120
𝑘 6000 < 𝑅𝑒𝐷 < 107
donde 𝐿
60 <
n=0.4 para calentamiento (𝑇𝑠 > 𝑇𝑏 ) 𝐷
n=0.3 para enfriamiento (𝑇𝑠 < 𝑇𝑏 )

Esta ecuación también NO introduce el factor empírico de corrección del perfil de velocidades por
efecto de la T°. Se debe usar si las diferencias de T° entre superficie y fluido son moderadas
situación que tiende a presentarse si el fluido no es muy viscoso.

Las propiedades se evalúan a la temperatura media del volumen del flujo.

 CONVECCIÓN FORZADA
FLUJO EN INTERIOR DE DUCTOS Y TUBERÍAS CIRCULARES
RÉGIMEN TURBULENTO
ECUACIÓN DE SIEDER Y TATE (1936) PARA FLUIDOS EN GENERAL
0.7 < 𝑃𝑟 < 10000
0.14 6000 < 𝑅𝑒𝐷 < 107
0.8
1 𝜇𝑏 𝐿
𝑁𝑢𝐷 = 0.027 𝑅𝑒𝐷 𝑃𝑟 3 60 <
𝜇𝑠 𝐷

La ecuación de Sieder y Tate es adecuada para situaciones con temperatura superficial uniforme.
Sin embargo, para considerar la variación de las propiedades físicas producto del gradiente axial de T°, las
temperaturas superficial y del fluido deberán tener valores medios entre la entrada y salida del ducto.

En ductos con sección transversal no circular, las ecuaciones de Dittus-Boelter y de Sieder y Tate se pueden
emplear reemplazando el diámetro D por el diámetro hidráulico.

Las propiedades, excepto la viscosidad en la superficie, se evalúan a la temperatura media del volumen
del flujo. Apta para flujos con elevadas diferencias de T° entre superficie y fluido.
 CONVECCIÓN FORZADA
FLUJO EN INTERIOR DE DUCTOS LARGOS EN RÉGIMEN TURBULENTO

ECUACIÓN DE KAYS Y LONDON (1984) PARA FLUIDOS EN GENERAL

𝑛 0.02 con temperatura superficial uniforme 𝑇𝑆

𝑇𝑏 C 0.7 < 𝑃𝑟 < 10000
𝑁𝑢𝐷𝐻 = 𝐶 𝑅𝑒𝐷0.8 𝑃𝑟 0.3 0.021 con flujo de calor uniforme 𝑞𝑠 ′′ 6000 < 𝑅𝑒𝐷 < 107
𝐻 𝑇𝑠 𝐿
0.0575 para calentamiento 60 <
C 𝐷
0.150 para enfriamiento

La ecuación de Kays y London, para considerar la variación de las propiedades físicas producto del
gradiente axial de T°, considera valores medios de las temperaturas superficial y del fluido, entre la
entrada y salida del ducto.

Las propiedades , excepto la viscosidad en la superficie, se evalúan a la temperatura media del volumen
del flujo.
 CONVECCIÓN FORZADA
FLUJO EN EXTERIOR DE CILINDROS
FLUJO POTENCIAL (INCOMPRESIBLE Y NO VISCOSO) SOBRE UN CILINDRO

Al despreciar el efecto de la viscosidad del fluido, la solución obtenida es ideal y corresponde al

denominado Flujo Potencial. No existe capa límite.
Este flujo genera un perfil de velocidades y presiones simétrico respecto de los ejes vertical y horizontal
del flujo.
Las fuerzas de arrastre viscoso y de forma son nulas.
 CONVECCIÓN FORZADA
FLUJO EN EXTERIOR DE CILINDROS

FLUJO REAL (INCOMPRESIBLE Y VISCOSO) SOBRE UN CILINDRO

Al considerar el efecto de la viscosidad del fluido, se genera la

capa límite en la zona cercana al sólido y la solución ideal sólo
permanece en el exterior de la capa límite.

Este flujo genera un perfil de velocidades y presiones

simétrico respecto del eje horizontal pero asimétrico respecto
del vertical.

Las fuerzas de arrastre viscoso y de forma son relevantes.

 CONVECCIÓN FORZADA
FLUJO EN EXTERIOR DE CILINDROS

FLUJO REAL (INCOMPRESIBLE Y VISCOSO) SOBRE UN CILINDRO

A partir de un cierto punto sobre el

borde sólidos, denominado punto de
separación, se genera una inversión
en el perfil de velocidades, dando
origen a vórtices que son propios del
régimen turbulento.

Por lo tanto, una parte del flujo se

desarrolla en régimen laminar y la
otra en régimen turbulento.
 CONVECCIÓN FORZADA
FLUJO EN EXTERIOR DE CILINDROS EN RÉGIMEN LAMINAR Y
TURBULENTO
ECUACIÓN DE ZUKAUSKAS (1972) PARA LÍQUIDOS Y GASES. Ts UNIFORME

𝑚 0.25
ℎത 𝐷 𝑈∞ 𝐷 𝑃𝑟
𝑁𝑢𝐷 = =𝐶 𝑃𝑟 𝑛
𝑘 𝜈 𝑃𝑟𝑠

n = 0.37 si Pr < 10
n = 0.36 si Pr > 10

Las propiedades, excepto Prs que se evalúa a la T° superficial, se evalúan a

la temperatura de corriente libre T∞.
 CONVECCIÓN FORZADA
FLUJO EN EXTERIOR DE CILINDROS NO PERPENDICULARES AL FLUJO
ECUACIÓN DE GROEHN (1986)

𝑁𝑢𝐷 = 0.206 𝑅𝑒𝑁0.63 𝑃𝑟 0.36 𝑅𝑒𝑁 = 𝑅𝑒𝐷 𝑠𝑒𝑛(𝜃)

Donde el número de Reynolds se basa en la

componente de flujo normal al eje del
cilindro.
Por tanto el Angulo 𝜃es el ángulo entre la
dirección del flujo y el eje del cilindro.
Ecuación válida para 2500 < ReN < Recrit
 CONVECCIÓN FORZADA
FLUJO EN EXTERIOR DE TUBOS
APLICACIÓN A INTERCAMBIADORES DE CALOR
 CONVECCIÓN FORZADA
FLUJO EN EXTERIOR DE BANCOS O HACES DE TUBOS
APLICACIÓN A INTERCAMBIADORES DE CALOR
 CONVECCIÓN FORZADA
FLUJO EN EXTERIOR DE BANCOS O HACES DE TUBOS
APLICACIÓN A INTERCAMBIADORES DE CALOR

𝑈𝑚𝑎𝑥 𝐷
𝑅𝑒𝐷 =
𝜈
 CONVECCIÓN FORZADA
FLUJO EN EXTERIOR DE BANCOS O HACES DE TUBOS
APLICACIÓN A INTERCAMBIADORES DE CALOR

En general, las ecuaciones de correlación adimensional son de la forma

0.25
𝑃𝑟
𝑁𝑢𝐷 = 𝐶𝑅𝑒𝐷𝑚 𝑃𝑟 0.36
𝑃𝑟𝑠

Donde el subíndice S indica que la propiedad Pr se debe evaluar a la

temperatura de la superficie del tubo y las demás propiedades se
evalúan a la temperatura media de film.
 CONVECCIÓN FORZADA
FLUJO EN EXTERIOR DE BANCOS DE TUBOS
APLICACIÓN A INTERCAMBIADORES DE CALOR
Tubos en línea
0.25
𝑃𝑟
𝑁𝑢𝐷 = 0.8𝑅𝑒𝐷0.4 𝑃𝑟 0.36
𝑃𝑟𝑠

Donde el subíndice S indica que

la propiedad Pr se debe evaluar
a la temperatura de la superficie
del tubo y las demás
propiedades se evalúan a la
temperatura media de film.
Válida para 10 < ReD < 100.
Régimen laminar.
 CONVECCIÓN FORZADA
FLUJO EN EXTERIOR DE BANCOS DE TUBOS
APLICACIÓN A INTERCAMBIADORES DE CALOR
Tubos en línea

0.25
𝑃𝑟
𝑁𝑢𝐷 = 0.27𝑅𝑒𝐷0.63 𝑃𝑟 0.36
𝑃𝑟𝑠

Válida para 10+3 < ReD < 10+5.

Régimen de transición
 CONVECCIÓN FORZADA
FLUJO EN EXTERIOR DE BANCOS DE TUBOS
APLICACIÓN A INTERCAMBIADORES DE CALOR
Tubos en línea

0.25
𝑃𝑟
𝑁𝑢𝐷 = 0.021 𝑅𝑒𝐷0.84 𝑃𝑟 0.36
𝑃𝑟𝑠

Válida para 2x10+5< ReD.

Régimen turbulento
 CONVECCIÓN FORZADA
FLUJO EN EXTERIOR DE BANCOS DE TUBOS
APLICACIÓN A INTERCAMBIADORES DE CALOR
Tubos alternados
0.25
𝑃𝑟
𝑁𝑢𝐷 = 0.9 𝑅𝑒𝐷0.4 𝑃𝑟 0.36
𝑃𝑟𝑠

Donde el subíndice S indica que

la propiedad Pr se debe evaluar
a la temperatura de la superficie
del tubo y las demás
propiedades se evalúan a la
temperatura media de film.
Válida para 10 < ReD < 100
 CONVECCIÓN FORZADA
FLUJO EN EXTERIOR DE BANCOS DE TUBOS
APLICACIÓN A INTERCAMBIADORES DE CALOR
Tubos alternados
0.25
𝑃𝑟
𝑁𝑢𝐷 = 0.022 𝑅𝑒𝐷0.84 𝑃𝑟 0.36 Si Pr > 1
𝑃𝑟𝑠

𝑁𝑢𝐷 = 0.019 𝑅𝑒𝐷0.84 Si Pr = 0.7

Ambas ecuaciones válidas para

2x10+5 < ReD.
Régimen turbulento
3.2 CONVECCIÓN NATURAL
FUNDAMENTOS
Desde un punto de vista de la
Mecánica de Fluidos, el régimen de
flujo (régimen laminar o régimen
turbulento) queda determinado por la
relación existente entre las fuerzas
motrices (Empuje-Peso) y las fuerzas
resistentes de origen viscoso.

Aplicación: calefactor de agua caliente

 CONVECCIÓN NATURAL. CASO PLACA PLANA VERTICAL
Número de Grashof (Franz Grashof; 1826-1893)

Por lo anterior el número de Reynolds

no es un buen predictor del tipo de
régimen de flujo (relacionada las
fuerzas de inercia con las de origen
viscoso) y por lo tanto, tampoco es un 𝑔𝛽 𝑇𝑠 − 𝑇∞ 𝐿3
parámetro relevante para predecir el 𝐺𝑟 =
𝜈2
coeficiente de transferencia de calor. En donde
𝑔=aceleración de gravedad.
𝛽= coeficiente de expansión volumétrica.
Se define, por tanto, un nuevo 𝑇𝑠 =temperatura de una superficie.
parámetro adimensional denominado 𝑇∞ =temperatura ambiente.
𝜈= viscosidad cinemática.
Número de Grashof. L= longitud característica del fenómeno.
 CONVECCIÓN NATURAL. CASO PLACA PLANA VERTICAL
Coeficiente de expansión volumétrica
En la ecuación del número de Grashof Para gases ideales y con base en la
es fundamental el coeficiente de ecuación de estado pV=mRT, es
expansión volumétrica del fluido β, posible calcular β, como:
definido como
1
𝛽=
𝑇∞

El número de Grashof representa la

relación entre las fuerzas de empuje
efectivo y las fuerzas resistentes de
1 𝜕𝜌 𝜌∞ − 𝜌 naturaleza viscosa.
𝛽=− ≅
𝜌 𝜕𝑇 𝜌 (T−𝑇∞ )
 CONVECCIÓN NATURAL
CASO PLACA PLANA VERTICAL

La figura muestra correlaciones

adimensionales entre el número de
Nusselt y el número de Grashof, en
régimen laminar y en régimen
turbulento, obtenidos
experimentalmente. En este gráfico
se puede apreciar un nuevo
parámetro adimensional como lo es
el número de Rayleigh, Ra:

Ra = Gr Pr
CONVECCIÓN NATURAL. CASO PLACA PLANA VERTICAL
Número de Rayleigh (Lord Rayleigh; 1842-1919)
𝑔𝛽
𝑅𝑎𝐿 = 𝑇𝑠 − 𝑇∞ L3 = GrL Pr
𝜈𝛼

Donde:
𝑅𝑎𝐿 : Numero de Rayleigh
Gr𝐿 : numero de Grashof
Pr: numero de Prandtl
𝑔: Aceleración de gravedad
𝐿: Longitud característica
𝑇𝑝 : Temperatura pared
𝑇∞ : Temperatura del fluido alejado de la pared o corriente libre
𝛼: difusividad térmica 𝑇𝑠 − 𝑇∞
𝑇𝑓 =
𝛽:coeficiente de expansión térmica 2
 CONVECCIÓN NATURAL. CASO PLACA PLANA VERTICAL
RÉGIMEN LAMINAR
En convección natural por tanto, la correlación adimensional será de la forma
general siguiente:

𝑁𝑢𝐿 = 𝑁𝑢𝐿 𝐺𝑟𝐿 , 𝑃𝑟

En el caso de una placa plana vertical, en régimen laminar
1
1 𝐺𝑟𝑥4 𝑘 Rango de validez:
ℎ𝑥 = 0.508 𝑃𝑟 2 1𝑥
0.952 + 𝑃𝑟 4
10 < GrL Pr < 10+8
Y el espesor de la capa límite queda dado por:
1
𝑃𝑟 + 0.56 4
𝛿 𝑥 = 4.3𝑥
𝑃𝑟 2 𝐺𝑟𝑥
 CONVECCIÓN NATURAL
CASO PLACA PLANA VERTICAL EN RÉGIMEN LAMINAR
Y el coeficiente medio de transferencia de calor para una placa de altura L,
sometida a un régimen laminar será:
1
1 𝑥=𝐿 1 𝐺𝑟𝐿4 𝑘
ℎത = න ℎ𝑥 𝑑𝑥 = 0.68 𝑃𝑟 2 1𝐿
𝐿 𝑥=𝑜
0.952 + 𝑃𝑟 4

Que en su forma adimensional se puede expresar como

(Ec. de Gryzagoridis, 1971): ℎത 𝐿
𝑁𝑢𝐿 =
𝑘
1
1 𝐺𝑟𝐿4
= 0.68 𝑃𝑟 2
1
0.952 + 𝑃𝑟 4
 CONVECCIÓN NATURAL
CASO PLACA PLANA VERTICAL EN RÉGIMEN TURBULENTO
Ecuación de McAdams, 1954:

ℎത 𝐿 1
𝑁𝑢𝐿 = = 0.13 𝐺𝑟𝐿 𝑃𝑟 3
𝑘

Rango de validez: GrL > 10+9

 CONVECCIÓN NATURAL
CASO PLACA PLANA VERTICAL INCLINADA (Ec. de Fujii e Imura, 1972)

1
𝑁𝑢𝐿 = 0.56 𝐺𝑟𝐿 Pr 𝑐𝑜𝑠𝜃 4

Rango de validez: 105 < 𝐺𝑟𝐿 Pr 𝑐𝑜𝑠𝜃 < 1011 y 0° ≤ 𝜃 ≤ 89°

CONVECCIÓN NATURAL. PLACAS HORIZONTALES
Superficie superior caliente o
superficie inferior fría
1
𝑁𝑢𝐿 = 0.54𝑅𝑎𝐿4 (105 ≤ 𝑅𝑎𝐿 ≤ 107 )
1
𝑁𝑢𝐿 = 0.15𝑅𝑎𝐿3 (107 ≤ 𝑅𝑎𝐿 ≤ 1010 )

Superficie inferior caliente o

superficie superior fría
1
𝑁𝑢𝐿 = 0.27𝑅𝑎𝐿4 (105 ≤ 𝑅𝑎𝐿 ≤ 1010 )

á𝑟𝑒𝑎 𝑑𝑒 𝑙𝑎 𝑠𝑢𝑝𝑒𝑟𝑓𝑖𝑐𝑖𝑒
𝐿=
𝑝𝑒𝑟í𝑚𝑒𝑡𝑟𝑜
CONVECCIÓN NATURAL. CILINDROS HORIZONTALES

Rango de validez: 1
10+3 < GrL < 10+9 𝑁𝑢𝐷 = 0.53 𝐺𝑟𝐷 𝑃𝑟 4

Pr > 0.5
CONVECCIÓN NATURAL. ESPACIOS CERRADOS
Ec. de Berkovsky y Polevikov

Definición de N° de Grashof

𝑔𝛽 𝑇1 − 𝑇2 𝛿 3
𝐺𝑟𝛿 =
𝜈2
CONVECCIÓN NATURAL. ESPACIOS CERRADOS

1
−4 0.28
𝐿 𝑃𝑟
𝑁𝑢𝛿 = 0.22 𝑅𝑎𝛿
𝛿 0.2 + 𝑃𝑟

Rango de validez:
𝐿
2 < < 10 𝑃𝑟 < 10 𝑦 𝑅𝑎𝛿 < 1010
𝛿
 CONVECCIÓN NATURAL. ESPACIOS CERRADOS

0.29
𝑃𝑟
𝑁𝑢𝛿 = 0.18 𝑅𝑎𝛿
0.2 + 𝑃𝑟

Rango de validez:

𝐿 𝑅𝑎𝛿 𝑃𝑟
1< <2 10−3 < 𝑃𝑟 < 105 𝑦 103 < < 1010
𝛿 0.2 + 𝑃𝑟
 CONVECCIÓN NATURAL. SUPERFICIES CON ALETAS
Ecuación de Raithby y Hollands, 1998

2.7 0.37
0.89 2.7 1
𝑅𝑎𝑠
𝑁𝑢𝑠 = + 0.62𝑅𝑎𝑠4 𝑁𝑢𝑠 = 𝐶 𝑅𝑎𝑠𝑏
18
CONVECCIÓN NATURAL o CONVECCIÓN FORZADA:
¿CÓMO DECIDIR?
El siguiente parámetro adimensional se emplea para analizar
cualitativamente el efecto que tienen en un problema dado,
la convección natural sobre la convección forzada, cuando
ambos fenómenos se encuentran presente.

Se ha demostrado que si el parámetro es cercano a 1.0

ambos fenómenos actúan con la misma intensidad, lo que
equivale a decir que las fuerzas de empuje (Gr), que
predominan en convección natural, son del mismo orden de
magnitud que las fuerzas de inercia (Re), que son las que
generan movimiento en convección forzada.

𝑔𝛽𝐿3 𝑇0 − 𝑇∞ /𝑣 2 GrL
2
= 2≈1
𝑈∞ 𝐿/𝑣 ReL
 COMPARACIÓN ENTRE CONVECCIÓN NATURAL Y CONVECCIÓN
FORZADA
CASO 1: Si GrL > (ReL)+2, entonces predominan
las fuerzas de empuje por sobre las de inercia en
el problema y se aborda como si fuese
convección natural.

CASO 2: Si GrL < (ReL)+2, entonces

predominan las fuerzas de inercia
por sobre las de empuje en el
problema y se aborda como si
fuese convección forzada.