Transferencia de Calor

Convección forzada
David Fuentes Díaz

Escuela de Ingeniería Mecánica

Universidad Industrial de Santander
 Convección forzada

Flujo interno

En flujo interno, el perfil de temperaturas depende de la

posición en la sección transversal y a lo largo del conducto:
T=T(x,r).
TS

R T Distribución de la
r temperatura
U Tf

Distribución de temperaturas en el flujo de tuberías

Tf = Tb (temperatura de masa)

Junio 2011 - Esc. Ing Mecánica UIS Transferencia de Calor - Convección Forzada 2
 Convección forzada

Flujo interno
R
Temperatura del fluido Tf: ρ ⋅ Cp ⋅ 2 ⋅ π ⋅ r ⋅ T(r) ⋅ dr
temperatura característica para cada ∫
Tf = 0 R
sección. Corresponde a un
∫ ρ ⋅ Cp ⋅ 2 ⋅ π ⋅ r ⋅ dr
promedio energético de la sección: 0

El calor transmitido por convección en cada posición se

refiere a dicha temperatura de masa:
dQ(x)=dA(x) h·(Ts(x) –Tf(x))= P(x) dx h·(Ts(x) –Tf(x))
Las propiedades necesarias para el cálculo del coeficiente
de convección se evalúan a la temperatura de masa
promediada entre la entrada y la salida del conducto:
Tfi + Tfo
2 3
Junio 2011 - Esc. Ing Mecánica UIS Transferencia de Calor - Convección Forzada
 Convección forzada

Flujo interno

Temperatura de pared constante.

Siendo P el perímetro del tubo, tenemos:
dT f h⋅P
dQ = h ⋅ P ⋅ dx ⋅ (Ts − T f ) = m ⋅ Cp ⋅ dT f =− ⋅ dx
(T f − TS ) m ⋅ Cp
integrando
h⋅ P ⋅ x T
T f − TS −
m⋅Cp TS
=e
T f i − TS
Tf
∆T0 ∆TL

Tfi

x
Junio 2011 - Esc. Ing Mecánica UIS Transferencia de Calor - Convección Forzada 4
 Convección forzada

Flujo interno

L L − h⋅ P ⋅ x

Q = ∫ dQ = − ∫ h ⋅P ⋅ (T fi − Ts ) ⋅ e m⋅Cp
⋅ dx = h ⋅ P ⋅ L ⋅ DMLT
0 0

siendo DMLT la Diferencia Media Logarítmica de

Temperatura: ∆T − ∆T
DMLT = 0 L ≈ (T − T )
∆T S B promedio
ln 0
∆T
L

Acotada por ∆T0 y ∆TL

min(∆T0 , ∆TL ) ≤ DMLT ≤ max( ∆T0 , ∆TL )

Junio 2011 - Esc. Ing Mecánica UIS Transferencia de Calor - Convección Forzada 5
 Convección forzada

Flujo interno

Flujo de calor uniforme.

Junio 2011 - Esc. Ing Mecánica UIS Transferencia de Calor - Convección Forzada 6
 Convección forzada

Flujo interno laminar

U⋅D
FLUJO LAMINAR Re D = < 2300
µ/ ρ

-Todas las propiedades se evalúan al promedio de la temperatura del fluido.

Tubo largo (flujo desarrollado) (constantes válidas sólo para tubos

cilíndricos)

 h⋅D 
Si TS es uniforme: Nu D =   = 3.66
 kf 
 Q  es uniforme:  h⋅D 
Si   Nu D =   = 4.36
A  kf 

Junio 2011 - Esc. Ing Mecánica UIS Transferencia de Calor - Convección Forzada 7
 Convección forzada

Flujo interno laminar

FLUJO LAMINAR
(Región de entrada)

Junio 2011 - Esc. Ing Mecánica UIS Transferencia de Calor - Convección Forzada 8
 Convección forzada

Flujo interno laminar

Para un tubo:

Ap=(π/4)D2

Pw= πD

Extensión de la longitud de entrada:

Leh= (Xent/D)H = 0.05Re
Régimen Laminar
LeT= (Xent/D)T = 0.05RePr

Le= X/D = 10 – 20 Régimen Turbulento

Junio 2011 - Esc. Ing Mecánica UIS Transferencia de Calor - Convección Forzada 9
 Convección forzada

Flujo interno laminar

Para secciones no circulares. Re=(ρUDh)/µ
Nu=(hDh)/k
Dh=4Ap/Pw

Junio 2011 - Esc. Ing Mecánica UIS Transferencia de Calor - Convección Forzada 10
 Convección forzada

Flujo interno laminar

TUBOS CONCENTRICOS EN REGIMEN LAMINAR
En el régimen laminar la transferencia de calor por convección se debe
exclusivamente a la interacción molecular de una capa de fluido con otra y no de la
posibilidad de mezcla cruzada de partículas macroscópicas de fluido.
qo qo

qi qi
Pared
adiabática

q =h o (Tso -Tm )
ho Dh
Nu o =
o k
q =h i (Tsi -Tm )
h Dh
Nu i = i
i k
4 π  D 2 -D 2 
4 o i 
= (Do -Di )
4 Ap
D h= =
Pw π (Do + Di )
Junio 2011 - Esc. Ing Mecánica UIS Transferencia de Calor - Convección Forzada 11
 Convección forzada

Flujo interno laminar

o A Temp. Constante

Junio 2011 - Esc. Ing Mecánica UIS Transferencia de Calor - Convección Forzada 12
 Convección forzada

Flujo interno laminar

o Flujo de calor uniforme.
Puede existir qo y qi simultáneamente,
Nuii Nuoo
Nui =
1- ( )
qo
qi
θ
*
i
Nuo =
1- ( )
qo
qi
θ
*
o

Junio 2011 - Esc. Ing Mecánica UIS Transferencia de Calor - Convección Forzada 13
 Convección forzada

Flujo interno turbulento

Cf
= St Pr =
2/3 Nu Dh Pr 2 / 3 TUBOS EN REGIMEN TURBULENTO
2 Re Dh Pr
Cf una correlación apropiada
para para régimen turbulento
0.8 1/3
Nu = 0.023Re Pr Ecuación de Colburn
0.8 n
Nu = 0.023Re Pr Ecuación de Dittus-Boelter

Donde: 0.7≤Pr≤160
n= 0.3 para enfriamiento del fluido (Tsup<Tm) Re≥10000
n= 0.4 para calentamiento del fluido (Tsup>Tm) L/D≥10
Tent + Tsal
Las propiedades del fluido se obtienen con la T =
2
CORRECCION POR TEMPERATURA
0.14
0.8 n µ
Nu = 0.023Re Pr µs es la viscosidad a la
µs temperatura de la superficie.

Junio 2011 - Esc. Ing Mecánica UIS Transferencia de Calor - Convección Forzada 14
 Convección forzada

Flujo interno turbulento

Corrección por el efecto de la LIQUIDOS GASES
variación de propiedades entre la m m m
f  Pr  µ  T 
temperatura de la superficie y el = s  ó  s ó  s
flujo f
B  PrB   µ B   TB 
n n n
Nu  Pr  µ  T 
= s  ó  s ó  s
Nu
B  PrB   µ B   TB 
Tipo de flujo Fluido Sentido transferencia calor m n
µs
Laminar Líquidos Fluido se calienta 0.58 -0.11
µB Fluido se enfría 0.50 -0.11

Ts
Gases Calentándose o enfriándose 1 0
TB

µs
Turbulento Líquidos Fluido se calienta 0.25 -0.11
µB Fluido se enfría 0.25 -0.25
Ts
Gases Fluido se calienta -0.2 -0.55
TB Fluido se enfría -0.1 0
Junio 2011 - Esc. Ing Mecánica UIS Transferencia de Calor - Convección Forzada 15
 Convección forzada

Flujo interno turbulento

FACTOR DE FRICCIÓN
• Tubería Lisa.
f = [0.79Ln (Re - 1.64 )]
-2
Petukhov

0.316 4
f = Re < 2x10
Re 0.23
Moody
0.184 4
f = Re ≥ 2x10
Re 0.2
Las propiedades del fluido se obtienen con la
• Tubería Rugosa. Tent + Tsal
T=
1 ε/d 2.54 2
= -2 log +
f 3.7 Re f
1.11
1 6.4 ε/d
= -1.8log +
f Re 3.7

Junio 2011 - Esc. Ing Mecánica UIS Transferencia de Calor - Convección Forzada 16
 Convección forzada

Junio 2011 - Esc. Ing Mecánica UIS Transferencia de Calor - Convección Forzada 17
 Convección forzada

Flujo externo placas

(T + T )
s ∞
Propiedades a la temperatura promedio Tm =
2
El desarrollo de la capa límite en función de la distancia viene caracterizado
U∞ ⋅ x
por el valor del número de Reynolds a lo largo de la placa
υ

< 5 ⋅105 → flujo laminar

Re
> 5 ⋅105 → flujo turbulento

Expresión del Número de Nusselt: es local, dependiente de la posición,

calculándose el coeficiente de convección según la expresión:
Nu x ⋅ k
hx =
x
Junio 2011 - Esc. Ing Mecánica UIS Transferencia de Calor - Convección Forzada 18
 Convección forzada

Flujo externo placas

ZONA FLUJO LAMINAR Li

Q= ∫h x ⋅ (Ts − T∞ ) ⋅ w ⋅ dx
0.6 ≤ Pr ≤ 50
Li −1
Válida para

Nusselt local Nu x =0.332 ⋅(Re x )1/2 ⋅Pr 1/3 obteniendo h local, hx

Nusselt promedio Nu 0 −L =0.664⋅(Re x =L )1/2 ⋅Pr 1/3 obteniendo hL
L
1
(Nusselt obtenido a partir de h promediado) = ∫ h X • dx
L0
Q = hL A(Ts − T∞ )

Junio 2011 - Esc. Ing Mecánica UIS Transferencia de Calor - Convección Forzada 19
 Convección forzada

Flujo externo placas

ZONA FLUJO TURBULENTO (Válida para 0.6 ≤ Pr ≤ 50 )

Si 5·105 <Re <107 Nux=0.0288·Rex (0.8)·Pr1/3 Cfx=0.0576/Re0.2

(Hasta 108 con error inferior al 15%)

Valores promedio
[
C f = 0.072 Re −L0.2 − 0.0464( xc / L ) ]
[
Nu = 0.036Pr1/ 3 Re 0L.8 − 23200 ] No tiene en cuenta la capa límite laminar previa

Nu = 0.036Pr1/ 3 [Re 0.8

L − 850 ] Tiene en cuenta la capa límite laminar previa

Nu = hL / k Nu x = hx / k
τs
Re L = ρU ∞ L / µ Re x = ρU ∞ x / µ C fx =
ρU ∞2 / 2
Junio 2011 - Esc. Ing Mecánica UIS Transferencia de Calor - Convección Forzada 20
 Convección forzada

Flujo externo cilindros

Estela del fluido que se

desprende de la superficie

u=0

T +T
1.Determinación de las propiedades necesarias a T = ∞ S
m 2
C y m están dados en la siguiente tabla para el Re.
ReD C m
0.4 – 4 0.989 0.33
ρu D
Re = x ext
4 – 40 0,911 0,385 µ
40 – 4000 0.683 0.466
4000 – 40000 0.193 0.618 Nu =C Re m Pr1/ 3
D
40000 – 400000 0.027 0.805 ext

Junio 2011 - Esc. Ing Mecánica UIS Transferencia de Calor - Convección Forzada 21
 Convección forzada

Flujo externo cilindros

CORRELACION DE ZHUKAUSKAS

 Cilindros.
1/4
m Pr n
ReD C m
N u D ext = CRe Pr 1-40 0.75 0.4
Prs
40-103 0.51 0.5
Propiedades a T∞ y Prs a Tsup 103-2 105 0.26 0.6
2 105-106 0.076 0.7
0.7 ≤Pr < 500
6
1 < Re < 10
n=0.36

 Esfera. 1/4
0.71 ≤Pr < 380
µ
[ 1/2
N u = 2 + 0.4Re + 0.06Re
D
2/3
D ]Pr 0.4

µs 3.5 < Re < 76x10 4

µ
1 ≤ ≤3.2
Propiedades a T∞ y µs a Tsup µs

Junio 2011 - Esc. Ing Mecánica UIS Transferencia de Calor - Convección Forzada 22
 Convección forzada

Flujo externo secciones no circulares

Junio 2011 - Esc. Ing Mecánica UIS Transferencia de Calor - Convección Forzada 23
 Convección forzada

Flujo externo bancos de tubos

Fluido en flujo cruzado

• Este tipo de arreglo es muy sobre un banco de tubos
utilizado en la industria, Flujo interno a
traves del tubo
siendo uno de los tantos
tipos de intercambiadores de
calor. Consiste en que un
fluido fluye dentro de los
tubos, mientras un segundo
fluido, a diferente
temperatura fluye a través
de ellos.

Junio 2011 - Esc. Ing Mecánica UIS Transferencia de Calor - Convección Forzada 24
 Convección forzada

Flujo externo bancos de tubos

SL SD
SL
D

D
ST
A1 ST
U∞ U∞ A2 A1
T∞ T∞

A2 A1

Disposición de los tubos en un banco. (a) Alineados (b) Alternados

SL= Separación Longitudinal. Nf= Numero de filas.

ST= Separación Transversal. ND= Numero de tubos por fila.
D= Diámetro del tubo.

Junio 2011 - Esc. Ing Mecánica UIS Transferencia de Calor - Convección Forzada 25
 Convección forzada

Flujo externo bancos de tubos

Nomenclatura usada en los bancos de tubos

• SL.: Distancia entre dos hileras consecutivas. Todas las correlaciones que rigen al flujo
externo a través de un banco de
• ST.: Paso transversal; distancia, transversal al flujo, entre
tubos están basadas en el ReDmax.
dos cilindros consecutivos.
• SD.: Paso diagonal (sólo para disposición triangular).
• NL.: Número de tubos en dirección longitudinal.
• NT.: Número de tubos en dirección transversal.
Donde :
D: Diámetro de los tubos
ν: Viscosidad cinemática
Umax: Velocidad máxima.

Para los efectos de transferencia de calor la disposición triangular es más eficiente que la
disposición alineada, pero la caída de presión en la primera disposición es mucho mayor. A la
hora de diseñar un banco de tubos se busca un compromiso, una transferencia de calor eficiente,
y una caída de presión baja; de lo contrario se necesitará de una bomba más eficiente para
transportar el flujo, lo que puede implicar un incremento en los costos.

Junio 2011 - Esc. Ing Mecánica UIS Transferencia de Calor - Convección Forzada 26
 Convección forzada

Flujo externo bancos de tubos

Correlación de GRIMISON.
NL ≥ 10
m ρVmaxD
Nu D = C 1 Re D,max 2000 < Re D,max < 40000 ReD,max ≡
µ
Pr = 0.7

Donde C1 y m se presentan en la siguiente tabla

Junio 2011 - Esc. Ing Mecánica UIS Transferencia de Calor - Convección Forzada 27
 Convección forzada

Flujo externo bancos de tubos

• De lo anterior se ha desarrollado otra Si NL<10 se toma en cuenta el

expresión para utilizarla con diferentes siguiente factor de corrección:
fluidos

Restricciones: Donde C2 viene dado en la

siguiente tabla:

Las propiedades se evalúan a la

temperatura promedio entre la superficie
y el fluido.
C1 y m se obtienen de la Tabla anterior
Junio 2011 - Esc. Ing Mecánica UIS Transferencia de Calor - Convección Forzada 28
 Convección forzada

Flujo externo bancos de tubos

La correlación mas reciente es la

propuesta por Zhukauskas:

Restricciones:

Las propiedades se evalúan en la media

aritmetica de las temperaturas del fluido y Si NL<20 se aplica el siguiente factor de
de la superficie, excepto Prs que es corrección
evaluado a la temperatura de la superficie.
C se obtiene de la siguiente Tabla:

Junio 2011 - Esc. Ing Mecánica UIS Transferencia de Calor - Convección Forzada 29
 Convección forzada

Flujo externo bancos de tubos

CALCULO DE LA VELOCIDAD MÁXIMA

ST Vmax
Para el arreglo alineado, Vmax = V
ST - D

Para el arreglo escalonado,

Vmax
2(SD - D)< (ST - D)
1/2
2
ST
2 ST + D
SD = S + L <
2 2
ST
Vmax = V
2(SD - D)

Junio 2011 - Esc. Ing Mecánica UIS Transferencia de Calor - Convección Forzada 30
 Convección forzada

Problema 1

De cual de las resistencias que se muestran en la figura se debe

suministrar más calor en la placa para que Ts sea constante. Cada
placa es de 1x1 [m2] y cada elemento de 4 cm.

Junio 2011 - Esc. Ing Mecánica UIS Transferencia de Calor - Convección Forzada 31
 Convección forzada

Problema 2

Dos habitaciones se encuentran ubicadas en la parte superior de un inmueble.

La longitud del techo de ambas es de 10 m , teniendo la pared correspondiente
un espesor de 0.2 m y una conductividad térmica equivalente, k = 0.8 W/m K. El
coeficiente de convección interior se estima en 10 W/m2 K. Dentro de las
habitaciones se mantiene una temperatura de 20 ºC. Considérese que en el
exterior la temperatura es de -10 ºC y que sopla un viento paralelo a la pared
del techo de 80 km/h. Calcular las pérdidas de calor para cada una de las
habitaciones.

Datos del problema:

Viento exterior 80 km/h y 263 K Temperatura del ambiente interior: Ti = 293 K
Temperatura del ambiente exterior: Te =263 K
Velocidad del aire: Ua =80 km h-1
Presión del aire ambiente: Pa =100000 Pa
Habitación 1 Habitación 2
Longitud de cada uno de los techos: L =10 m
Espesor de la pared del techo: ep =0.2 m
hi = 10 W m-2 K-1 Coeficiente de convección interior: hi = 10 W m-2 K-1
Ti = 293 K Conductividad equivalente de la pared: k pared = 0.8 W m-1 K-1

Junio 2011 - Esc. Ing Mecánica UIS Transferencia de Calor - Convección Forzada 32
 Convección forzada

Problema 3

La anemometría de hilo caliente es una técnica utilizada para determinar la

velocidad puntual de un fluido. El método se basa en medir la intensidad eléctrica
que debe pasar por un filamento de material conductor inmerso en la corriente de
fluido, para mantener su temperatura constante, compensando el enfriamiento
producido por la circulación del fluido en torno al conductor. Suponer un filamento
de 0.2 mm de diámetro que se quiere mantener a una temperatura constante de 77
ºC en una corriente de aire a 27 ºC. Determinar la corriente eléctrica requerida si la
velocidad del aire es de 10 m/s y la resistividad del platino que compone el
filamento de 17.1*10-5 ohmios*m.

Ieléctrica Datos del problema:

Temperatura del aire Ta = 300 K
Temperatura del filamento Tf= 350 K
Velocidad del aire Ua =10 m s-1
Diámetro del filamento Df = 0.2 10-3 m
Tf Ua Ta Resistividad del filamento: ρ= 17.1 10-5 Ω m

Junio 2011 - Esc. Ing Mecánica UIS Transferencia de Calor - Convección Forzada 33
 Convección forzada

Problema 4

Un termopar T1 se inserta en un conducto por el que circula aire caliente. Dicho termopar se
encuentra soldado en el extremo de una vaina metálica (conductividad 70 W/m K), de
longitud L = 0.03 m y diámetros interior y exterior respectivamente Di = 2 mm y Do = 3 mm. La
parte central del termopar está rellena de un material que puede considerarse perfectamente
aislante. Un segundo termopar T2 se utiliza para determinar la temperatura de la pared del
conducto. Considérese unas condiciones en las que el aire que circula por el conducto tiene
una velocidad de 20 m/s y los dos termopares registran temperaturas T1 = 450 K y T2 = 375 K.
Suponiendo que la temperatura indicada por el termopar T2 es la que existe en la base de la
vaina, y despreciando el calor transmitido por radiación, evaluar la temperatura del aire.

T1 T2
Datos del problema:
Longitud de la vaina del termopar: Lv =0.03 m
Diámetro interior: Di = 2 10-3 m
Ua Ta Diámetro exterior: De = 3 10-3 m
Vaina metálica Velocidad del aire: U aire = 20 m s-1
Conductividad de la vaina: kvaina = 70 W m-1 K-1
Temperatura del extremo: T1 = 450 K
Temperatura en la base de la vaina:T2 =375 K

Junio 2011 - Esc. Ing Mecánica UIS Transferencia de Calor - Convección Forzada 34
 Convección forzada

Problema 5
Un transformador de potencia eléctrica de 300mm de diámetro y 800mm
de altura, disipa 1000w. Se desea mantener su temperatura superficial a
47oc, suministrando glicerina a 24oc, por una tubería de pared delgada de
20mm de diámetro, soldada a la superficie lateral del transformador. Se
supone que se transfiere el calor a la glicerina. Suponiendo que ∆Tmax para la
glicerina es de 6oc. Determine el flujo másico de refrigerante que se requiere,
el espaciado lateral S entre cada vuelta.

Junio 2011 - Esc. Ing Mecánica UIS Transferencia de Calor - Convección Forzada 35
 Convección forzada

Problema 6
En un calentador de aire tubular (de dos pasos por el aire) de una caldera de vapor, se debe calentar un flujo másico de
aire ma= 21.5 kg/s desde una temperatura de te,aire=30 °C a una temperatura de salida ts,aire=260 °C.
Los gases de combustión fluyen a razón de mgases=19.6 kg/s y circulan por el interior de unos tubos de acero con
conductividad térmica k=46.5 W/m K, y con diámetro externo e interno de d2/d1=53/50 mm a una velocidad media de
Vgases=14 m/s. La temperatura de los gases a la entrada del calentador es de Te,gases=380 °C. Considerar los gases de
combustión como aire.
El aire se mueve transversalmente a través de los tubos con una velocidad media en la sección frontal de 0.8 m/s. Los
tubos están dispuestos en orden alternado con una separación longitudinal St=2.0 d2 y una separación longitudinal
Sl=1.5 d2.
Calcular:
El área de transferencia de calor requerida, la longitud de los tubos, y el número de tubos necesarios para lograr el
calentamiento del aire.

Aire

Gases

Te,gase
s=380º
C

Ltubo Ltubo

Te,aire=30ºC Ts,aire=260ºC

Junio 2011 - Esc. Ing Mecánica UIS Transferencia de Calor - Convección Forzada 36
 Convección forzada

Problema 7
Un calentador solar simple se compone de una placa plana de vidrio debajo de la cual hay una
bandeja poco profunda. En el espacio que queda entre la placa de vidrio y la bandeja circula agua
a una velocidad de 0.1 m/s. La radiación solar, que incide sobre la placa, pasa a través del vidrio a
razón de 500 W/m2. La placa de vidrio tiene una emisividad ε=0.85, con una absortividad α=0.1 y
una reflectividad r=0.05. El agua se encuentra en la entrada a una temperatura de 23°C, mientras que el aire
se encuentra a 23 °C y la temperatura del cielo es igual a -40°C. Si por el exterior del colector solar incide
aire a 0.5 m/s, calcule el calor ganado por el agua, así como la temperatura de salida del agua, si se asume que
la parte inferior de la bandeja se encuentra perfectamente aislada y que la cubierta de vidrio tiene una
temperatura constante en toda su área. El calentador solar tiene una profundidad (respecto del corte mostrado
en el papel) de 2 m.
Cubierta de
vidrio

0.07 m Agua

0.8 m Bandeja

NOTA: Para el cálculo del diámetro hidráulico para el cálculo del número Reynolds utilice el perímetro
mojado por el agua. Para el caso del cálculo del número de Nusselt, utilice el perímetro calentado.

Junio 2011 - Esc. Ing Mecánica UIS Transferencia de Calor - Convección Forzada 37
 Convección forzada

Problema 8
Se desea estimar las perdidas de calor a través de una chimenea. La chimenea tiene una altura
de 10m, diámetro 0.6m y un espesor de 1cm. Esta fabricado en fundición de k= 50w/mk. Se tiene
un flujo de gases procedentes de los hornos a razón de 1Kg/s entrando a la chimenea a 500oc.
Las propiedades de los gases se evalúan como las del aire asumiendo que la Tchim es constante y
circula aire alrededor de la chimenea a 20oc con una velocidad perpendicular a la chimenea de
80Km/h, calcular la temperatura de salida de los gases y la temperatura de la chimenea.
Tsal

Aire a 20oc
V=50 km/h

Gases de combustión a 500oc

Junio 2011 - Esc. Ing Mecánica UIS Transferencia de Calor - Convección Forzada 38