TRANSFERENCIA DE CALOR EN FLUIDOS

Durante el calentamiento de un fluido frío en contacto con una pared caliente se desarrolla la capa
límite laminar y turbulenta en el seno del fluido. En la zona laminar se tiene la mayor resistencia al
flujo de calor y esta resistencia es mínima en la zona turbulenta donde hay un intenso mezclado del
fluido motivado por la presencia de los remolinos.
En un fluido la transmisión
En un sólido la de calor es por colisiones .
transmisión de calor es Las partículas frías y calientes que
por la vibración Tsólido están en movimiento se encuentran,
Tfluido colisionan e intercambian energía.

Red de átomos
Partícula fría

Partícula
vibración zona
Capa Laminar caliente
caliente
Región
pared turbulenta
Colisión
red vibrando

dx

1#
 ¿Cómo se transmite calor en el sólido?
Por conducción
¿Cómo se transmite
calor en el fluido?
Por convección
¿Dónde está la
mayor resistencia
al flujo de calor
en el fluido?
En la capa
laminar
¿Por qué?
- k (dT) =de la
Q /A =característicos
Porque los torbellinos - kturbulencia
(dT) se encargan deQmezclar
= h A dT
el fluido mejorando
la transmisión de calor dx dx
Siendo h el coeficiente pelicular que se debe calcular en cada caso en particular debido a que:
Analizando la fluido
k depende del pérdiday de
de calor en la capamedia
la temperatura laminar
de del fluido,
la capa desde el punto de vista de la
laminar
laconducción:
temperatura media de la capa laminar a su vez depende del espesor de la capa laminar
El espesor de la capa laminar es muy difícil de calcular porque depende de muchas variables:
viscosidad, velocidad, rugosidad de la superficie sólida, geometría, si el fluido se está enfriando, en
ebullición o si condensa, etc.etc.
2#
 CLASIFICACIÓN

•Convección natural: La diferencia de temperaturas que existe en el sistema provoca diferencias

de densidad en el seno del fluido y esto genera el movimiento del fluido.

•Convección forzada: La transmisión de calor se debe a corrientes de remolinos existentes en un

fluido con movimiento turbulento.

Según el tipo de calor intercambiado

•Convección sin cambio de estado (Calor sensible)

•Convección con cambio de estado (Calor latente)

De acuerdo a la configuración:

• Flujo en conductos

• Flujo alrededor de objetos sumergidos

3#
 FLUJO EN EL INTERIOR DE CONDUCTOS
CILINDRICOS

Q = velocidad de transferencia de calor

De acuerdo a la Ley de
enfriamiento de Newton: h = coeficiente convectivo W/m2C
A = Área de transferencia: 2  DL
D: diámetro interior del conducto
Q = h A T L: longitud del conducto
T = Diferencia de Temperatura: Hay
distintas formas de medirla.

4#
 DIFERENCIA DE TEMPERATURA PARA DISTINTOS CASOS DE FLUJO EN CONDUCTOS

FLUJO PARALELO CONTRACORRIENTE

Tho Ths Tho T
Ths co

Tcs
Tco Tcs Tho
Ths Tho
Ths
1) Intervalo de temperatura
1) Intervalo de temperatura
T1 = Tho – Tcs T2 = Ths - Tco
T1 = Ths – Tcs T2 = Tho – Tco
2) T aritmética 2) T aritmética
T1 + T2 T1 + T2
Ta = Ta = 2
2
3) T Media Logarítmica 3) T Media Logarítmica
T2 - T1 T2 - T1
TML = TML =
ln T2 ln T2
T1 T1

TML = (Tho - Tco ) - (Ths - Tcs) (Ths - Tco) - (Tho - Tcs)

(T - T ) TML =
ln (T - T )
ho co ln hs co
(Ths - Tcs) (Tho - Tcs)
5#
 SON LAS TEMPERATURAS MEDIAS DE LA CORRIENTE
Ths Tho Tcs Tco

TML ES EL MÁS CONVENIENTE PARA USAR

Si 0,5 < (T2 / T1) < 2,0 Ta solo tiene 4% de error

Donde:
TML : Diferencia media logarítmica entre los fluídos en circulación.
T ML = ( T max - T min )/ [ ln ( Tmax / Tmin)]
El término diferencia media logarítmica de temperatura se presenta debido a que la temperatura de
uno de los dos fluidos en consideración varía de acuerdo con su recorrido en la dirección de flujo.
A través de una longitud L, en un punto donde L = 0 existirá un T entre los dos fluidos y en un
punto donde L = L existirá otra T realizándose la medición de L en el sentido del flujo de uno de
los dos fluidos. Indiferente del sentido recorrido, T max. Será el T mayor y el T min. Será el T
menor en cada uno de los puntos citados.
Entre los fluidos en circulación se cumplirá la ley de conservación de la energía, de tal modo que
el calor cedido por uno de ellos al enfriarse, será tomado por otro al calentarse . Considerando por
ejemplo, que el fluido que circula por la tubería interior varía su temperatura durante el
intercambio de calor, se podrá calcular la velocidad de transferencia de calor, si se conoce la
cantidad de masa de fluido que circula durante un tiempo determinado, su calor específico
promedio y las temperaturas iniciales y finales.

6#
 DIFERENCIA DE TEMPERATURA
BALANCE ENTÁLPICO PARA DISTINTOS CASOS DE CALENTAMIENTO

1) CALENTAMIENTO CON UN VAPOR CONDENSANTE

Th = Temp. del vapor condensante

Temperatura oK

Th
T1 Tco = Temp. entrada fluido frío Tcs
T2 T
Tcs = Temp. salida fluido frío
T = Diferencia puntual de Temp.
Tco
T1 = Th - Tcs
Longitud del tubo, m T2 = Th - Tco
Tcs - Tco = T del fluido frío

Balance entálpico:

mh  = mc Cpc (Tcs - Tco) = q

mh [ + Cph (Ths - Tho) ] = mc Cpc (Tcs - Tco) = q

7#
 DIFERENCIA DE TEMPERATURA
BALANCE ENTÁLPICO PARA DISTINTOS CASOS DE CALENTAMIENTO
2) CALENTAMIENTO CON OTRO FLUIDO

FLUJO PARALELO CONTRACORRIENTE

Tho Tho

Temperatura oK
Temperatura oK

Ths Ths
T2 T T1 T T1
T2
Tcs
Tco Tco Tcs

Longitud del tubo, m

Longitud del tubo, m

Tho = Temp. entrada fluido caliente

Ths = Temp. salida fluido caliente
Tco = Temp. entrada fluido frío
Tcs = Temp. salida fluido frío
Tcs - Tco = T del fluido frío
Ths -Tho = T del fluido
caliente ; T2 = Ths - Tco
T = Diferencia puntual de
Temp.
Th 1(T
mh Cp = hsT-hsT-hoT)cs] = mc Cpc (Tcs - T; co) T
= q2 =
Tho - Tco ; T1 =
8#
Tho - Tcs
EJEMPLO: Área de intercambio para flujo en contracorriente y en paralelo.
Un caudal másico de 0,5 kg/s de aceite (Cp=2090 J/kgoK) se enfría de 375o a 350oK
intercambiando calor con agua que se calienta de 280o a 311oK. Calcular el área de intercambio
cuando el flujo es en paralelo y en contracorriente.
FLUJO PARALELO CONTRACORRIENTE
Tho=375
Tho=375 Ths=350 Ths=350

Tco=280
Tcs=311 Tco=280 Tcs=311

T1 = Ths - Tcs = 350-311 = 39oK T1 = Tho - Tcs = 375-311 = 64oK

T2 = Tho - Tco = 375-280 = 95oK T2 = Ths - Tco = 350-280 = 70oK

TML = (Tho - Tco ) - (Ths - Tcs) TML = (Ths - Tco) - (Tho - Tcs)
(T - T ) (Ths - Tco)
ho
ln co ln
(Ths - Tcs) (Tho - Tcs)

(375 - 280 ) - (350 - 311) (350 - 280 ) - (375 - 311)

TML = TML =
(375 - 280 ) (350 - 280 )
ln (350 - ln (375 -
311) 311)

TML = 63oK TML = 67oK

9#
 • Cálculo del calor entregado por el aceite
Q = m Cp T = 0,5 kg/s 2090 J/kg.oK 25 oK = 26125 W

FLUJO PARALELO CONTRACORRIENTE

Tho=375 Ths=350 Tho=375
Ths=350
Tco=280
Tcs=311 Tco=280 Tcs=311

• Cálculo del A • Cálculo

del A Q 26125 W
Q 26125 W A= =
A= = h T 250 W m-2 K-1 67 oK

h TA = 1,66 m2250 W m-2 K-1 A = 1,56 m2

63 oK

Asumimos que el coeficiente convectivo es igual para ambas configuraciones

10#
 EL COEFICIENTE CONVECTIVO h
En consecuencia si tenemos distintas formas de medir el T también tendremos distintos
valores de h:
1) h basado en la diferencia de temperaturas de entrada
Q = h1  D L (Tb1 – T01)

2) h basado en la media aritmético de diferencia de temperaturas de entrada y salida

(Tb1 – T01) + (T02 – Tb2)

Q = ha D L
2
3) h basado en la media logarítmica de diferencia de temperaturas
(T01 – Tb1 ) - (T02 – Tb2)
Q = hML  D L
ln (T01 – Tb1 )
(T02 – Tb2)
La forma más conveniente para expresar la ecuación que gobierna la transferencia de calor entre el
fluido de mayor temperatura y el de menor temperatura es utilizando TML

Q = h MLA T ML

11#
POSIBLES CAMINOS PARA DETERMINAR EL COEFICIENTE CONVECTIVO

1)Resolución analítica de las ecuaciones de continuidad, movimiento y energía. Sólo es

aplicable a geometrías sencillas.
2)Utilización de las analogías entre transferencia de calor y cantidad de movimiento.
3)Métodos semiempíricos basados en el análisis dimensional y resultados experimentales.
Podemos sistematizar y simplificar este trabajo si recurrimos al Teorema P que:
➢Requiere identificar perfectamente los parámetros aplicables a cada situación:
➢ geometría
➢ propiedades del fluido
➢ régimen de flujo
➢ diferencia de temperatura
➢ flujo interno o flujo externo
➢ convección forzada o convección libre
➢Manipularlos algebraicamente para obtener una funcionalidad de una menor cantidad
de grupos adimensionales
h = f (a1, a2, a3, …) € 1 = f (2,
3)
➢Apoyarse en datos experimentales y correlaciones existentes para encontrar la
funcionalidad entre ellos. De no encontrarlas deberemos realizar nuestras propias
experiencias en laboratorio.
12#
 CONVECCIÓN FORZADA CONVECCIÓN LIBRE

AIRE
AIRE

El calor es barrido por la corriente El calor se transmite hacia arriba por el

forzada de aire aire caliente que asciende
Las características del flujo están dadas Las características del flujo están dadas por
por una fuerza externa. el efecto de flotación del fluido caliente.

Primero se halla el perfil de velocidad, que se Los perfiles de velocidad y de temperatura

emplea para calcular el perfil de temperatura. están íntimamente relacionados.

Nu = f (Re, Pr, L/D, Br) Nu = f (Gr, Pr, L/D, Br)

Nu = Calor transferido por convección

Calor transferido si solo hubiera conducción
13#
 CONVECCIÓN FORZADA: LA SENSACIÓN TÉRMICA

La SENSACION TERMICA es un índice que combina el efecto de la temperatura con la

velocidad del viento. El cuerpo pierde calor más rápidamente cuanto más fuerte es el viento y
menor es la temperatura. Se siente como si la temperatura fuera menor de lo que en realidad es.

SENSACIÓN

Oc
TÉRMICA
AIRE A 7Oc

DÍA CALMO VIENTO 15 M/S SENSACIÓN TÉRMICA

T = 7OC T = 7OC
S.T.= 7OC S.T.: -7Oc
VELOCIDAD DEL VIENTO m/s

Es la temperatura externa que tendría que haber a una velocidad del viento de 6 km/h
(velocidad promedio a la que camina una persona) de modo que la sensación térmica sea
idéntica a la de la situación real.
La menor sensación térmica registrada fue de -92oC en Pelly Bay, Canadá, y correspondió a
una temperatura de -51oC y una velocidad del viento de 60 km/h.

14#
10 4 -1 -7 -12 -18 -23 -29 -34 -
40

15#
 CORRELACIONES PARA CALCULAR EL COEFICIENTE CONVECTIVO h

El punto de partida para poder usar la ecuación que gobierna la transferencia convectiva de calor es
poder predecir el coeficiente pelicular h

Q = h MLA T ML

El calculo del coeficiente de película h se basa en correlaciones de la forma:

Nu = f (Re, Pr, L/D, Br) Convección forzada

Nu = f (Gr, Pr, L/D, Br) Convección libre

Calor transferido por convección

Nu = =hD/k
Calor transferido si solo hubiera conducción

Estudiaremos a continuación cómo surgen estas funcionalidades y algunas correlaciones disponibles

para los casos de flujo interno o externo, convección forzada, libre o combinadas y otras situaciones
posibles

16#
 CONVECCIÓN FORZADA: NÚMEROS ADIMENSIONALES
Variables involucradas: M L T 
velocida v 0 1 -1 Variables independientes:
d  0
densidad
viscosidad  1 -3 -1
-1 0 0 D 0 1 0 0
0 k 1 1 -3 -
conductividad térmica del k 1 1 -3 -1  1 -1 -1 1 0
fluido 0 =
diámetro D 0 1 0 0 v 0 1 -1 0 /
calor específico a presión cte Cp 0 2 -2 -1
coeficiente convectivo h 1 0 -3 -1

7 variables - 4 (rango de matriz dimensional) = 3 grupos 

◆  1 = Da. kb . c. vd. 
1 = La x (M LT-3 -1)b x (M L-1T-1)c x (L T-1)d x (M L-3)

M: b+c+1=0 a=1
L: a+b-c+d-3= b=0
0 T :-3b - c - d = 0 c=-1
d=1 1=Dv  Re
 : -b = 0 

17#
◆  2 = De. kf . g. vh. Cp
2 = Le x (M LT-3 -1)f x (M L-1T-1)g x (L T-1)h x (L2 T-2 -1)

M: f+g=0 e=0
L: e +f-g+h+2= f=-1
0 T :-3f - g - h - 2 = 0 g=1
 Cp
 : -f - 1 = 0 2= = Pr
h=0 k

◆  3 = Di. kj . k. vm. h

2 = Li x (M LT-3 -1)j x (M L-1T-1)k x (L T-1)m x (M T-3 -1)

M : j+ k + 1 = 0 i=1
L: i+j-k+m=0 j=-1
T : -3j - k - m - 3 = 0 k=0
 : -j - 1 = 0 m=0

hD
3= k= Nu

18#
 CONVECCIÓN FORZADA FLUJO INTERNO
CORRELACIONES EMPÍRICAS
1) Régimen laminar Re < 2100
Para (Re Pr D/L) < 100
0,085 (Re Pr D/L)
Nu = 3,66 + (b/ 0)0,14
1 + 0,045 (Re Pr D/L)0,66
Para (Re Pr D/L) > 100
Nu = 1,86 (Re Pr D/L)1/3 (b/0)0,14

2) Régimen turbulento

Para 104<Re<105
0,6<Pr<100
Nu = 0,026 Re0,8 Pr1/3 (b/0)0,14
L/D>100

Todas las propiedades físicas están evaluadas a la temperatura media del fluido (Tf1+Tf2)/2 a
excepción de 0 que se evalúa a la temperatura media de pared (T01+T02)/2

19#
Otras correlaciones para régimen turbulento:
• Diltus y Boelter

Para: n = 0,4 para fluido calentándose

n = 0,3 para fluido enfriándose
Nu = 0,023 Re0,8 Prn
Re > 104 , 0,7<Pr<100 ,
L/D>60
propiedades físicas evaluadas a
Tf
• Colburn
Para: Re >104 , 0,7<Pr<160 , L/D>60
Re y Pr evaluados a (T f1+T )/2
f2
St = 0,023 Re-0,2 Pr-2/3
St evaluado a (Tf+Tp)/2
NSTANTON : St = Nu Re-1 Pr-1
• Sieder y Tate para Pr elevados
Para: Re >104 St = 0,023 Re-0,2 Pr-2/3 (b/0)0,14
0,7<Pr<17000
L/D>60

20#
 CONVECCIÓN FORZADA FLUJO INTERNO
CORRELACIÓN DE SIEDER Y TATE
Modificando la ecuación para régimen laminar a la forma:
h
C G Pr2/3 ( / )-0,14 = 1,62 (D/L) (DG/)
1/3 -2/3
Velocidad másica G =  v
p b 0

y la ecuación para régimen turbulento a la forma:

h
Pr 2/3
CG ( / )-0,14
= 0,023 [1+ (D/L)0,7] (DG/)-0,2
p b 0

En un gráfico doble-log Sieder y Tate grafican Y en función del Re

(Tb2-Tb1) D hD h
Pr2/3
Pr 2/3 b 0-0,14
4 L ( / ) k Re-1 Pr-1/3 ( /b ) 0 ( b/ )0-0,14
-0,14
(Tb2-T )b1 ln CpG

La ecuación para régimen laminar es una recta de pendiente -2/3

La ecuación para régimen turbulento es una recta de pendiente -0,2
y se pueden interpolar los valores para la zona de transición 2100<Re<10000
La ordenada està expresada de distintas formas equivalentes

21#
 GRÀFICO DE SIEDER - TATE

(Tb2-Tb1) D
Pr2/3 (b/ 
)-0,14
4L 0
(Tb2-Tb1)ln

h D Re-1 Pr-1/3 (b/0)-0,14

k

h
Pr2/3 (b/ ) -0,14
0
CpG

Coeficiente de transmisión de calor para flujo totalmente desarrollado en tubos lisos

Las curvas para flujo laminar solo deben utilizarse para Re Pr D/L >> 10 y Tpared cte

22#
 Nu para flujo laminar completamente desarrollado en
tubos de diferente sección transversal

Nu para la región de entrada en flujo laminar en tubos de sección circular con:

a) temperatura superficial uniforme b) flujo de calor superficial uniforme

23#
RESUMEN DE CORRELACIONES DE CONVECCIÓN FORZADA FLUJO
INTERNO

Ver fig. anterior

Ver tabla anterior

Ver tabla anterior

24#
 SECCIONES NO CIRCULARES: DIÁMETRO EQUIVALENTE

Área de la sección transversal

D e = 4 Rh Rh =
perímetro mojado

◆ Para una sección circular

De = = D

◆ Para una sección cuadrada De =

◆ Para un ánulo
De =

EFECTOS DE ENTRADA EN CONDUCTOS CON L/D < 60

Para 2<L/D<20 hcorr/h = 1+(D/L)0,7

Para 20<L/D<60 hcorr/h = 1+6D/L

25#
EFECTO DE LA RUGOSIDAD ARTIFICIAL EN EL Nu

ROSCA TRIANGULAR RESALTO RECTANGULAR

ROSCA ONDULADA RESALTO TRIANGULAR

Nufl = 0,021 Ref 0,8 Pr 0,43 (Pr /Pr )0,25 

l
fl
fl
p
rug
Para 0,7<Pr<80 (s/h)opt= 13+/-1
rug = 1,04 Pr 0.04 exp [0,85 f (s/h)]
(s/h)>(s/h) opt , f (s/h) = (s/h)opt / (s/h)

(s/h)<(s/h)opt , f (s/h) = (s/h) / (s/h) opt

26#
 EJEMPLO: TRANSMISIÓN DE CALOR EN UN TUBO
En un sistema de enfriamiento entra etilenglicol a un tubo de Cu estirado en frío de 5 m de
longitud y 100 mm de diámetro con una v = 5 m/s. Estimar el coeficiente de transferencia de
calor si la pared del tubo se mantiene a 100oC y la temperatura media del fluido es 20oC.
Solución

1) Cálculo del Pr
20C = 1,92 x 10-5 m2/s k20C = 0,294 W/m0K (b/o)-0,14 = 1,38
Pr = 204 Pr1/3 = 5,89

2) Cálculo del Re Y = 0,0035

Dv 0,1 m . 5 m/s
Re    2,6 x 104
 1,92 x 10 m /s
-5 2

3) Cálculo de Y a partir del gráfico Re = 2,6 x 104

hD
Y k Re-1 Pr-1/3 ( /b )-0,14
o = 0,0035

0,249 W/moK h = 1842 W/m2 oK

h   . 2,6 x 10 . 5,89 . 1,38
4
0,1 m

27#
 EJEMPLO: TRANSMISIÓN DE CALOR EN UN TUBO
A un tubo de 2,2 m de largo y 12 mm de diámetro ingresan 0,083 kg/s de agua a 30oC.
Determinar la temperatura de salida del agua si la temperatura interior de la pared del tubo se
mantiene a 60oC.
Solución
1) Suponiendo Ts = 40oC >>>> Tm,fl = (Te+Ts)/2 = (30+40)/2 = 35oC
35C = 7,28 x 10-4 kg/ms ; 60C = 4,5 x 10-4 kg/ms ; (b/o)-0,14 = 0,935
k35C = 0,539 kcal/hm0C ; Cp = 1 kcal/kgoC

4 . 0,083 kg/s

Re 4  D  7,24 x 10-4 kg/ms .  . 0,012m  1,2 x 104 >>Turbulento

7,28 x 10-4 kg/ms . 1 kcal/kgoC

 4,86
Pr    Cp 
k 0,539 kcal/hm0C
Y = 0,004
D/4L = (0,012 m)/(4.2,2 m) = 1,36 x 10 -3

(Tb2-Tb1)
Pr2/3 (b/o)-0,14
Y= (T01-Tb1) - (T02-Tb2)
(T01-Tb1) D 4
ln L
(T02-Tb2)

28#
 (Tb2-30)
0,004 = 1,36 x 10-3 . 4,862/3 . 0,935
(Tb2-30) € Tb2 = 50oC
(30) (Valor distinto al supuesto, se debe iterar)
ln (60-Tb2)

2) Suponiendo Ts = 50oC >>>> Tm,fl = (Te+Ts)/2 = (30+50)/2 = 40oC

40C = 7,20 x 10-4 kg/ms ; 60C = 4,5 x 10-4 kg/ms ; (b/o)-0,14 = 0,936
k40C = 0,539 kcal/hm0C ; Cp = 1 kcal/kgoC

Re 
4  = 1,22 x 104
D
Y = 0,0039
  Cp
Pr   4,81

D/4L = 1,36k x 10-3

(Tb2-30)
0,0039 = 1,36 x 10-3 . 4,812/3 . 0,936
(Tb2-30)
Tb2 = 50 oC
(30)
ln (60-Tb2)

29#
 CONVECCIÓN FORZADA FLUJO EXTERNO
CORRELACIONES EMPÍRICAS PARA DISTINTAS FORMAS GEOMÈTRICAS

La transferencia de calor en objetos sumergidos se calcula: Q = h A (T0 - Too)

Transmisión de calor por convección forzada placa plana

Para: 103 < Re <

Nu = 0,664 Re Pr1/2 1/3
5x105 Pr > 0,5

Para: 2x 105 < Re < 3x107

Nu = 0,036 (Re0,8 – 9200) Pr0,43 (b/0)-0,25 0,7 < Pr < 380
0,26 < (b/0) < 3,5

Todas las propiedades físicas están evaluadas a la temperatura media del fluido (Tf1+Tf2)/2 a
excepción de 0 que se evalúa a la temperatura de pared.

30#
 CONVECCIÓN FORZADA FLUJO EXTERNO
CORRELACIONES EMPÍRICAS PARA DISTINTAS FORMAS GEOMÈTRICAS
Transmisión de calor por convección forzada a una esfera
aislada
Nu = 2 + 0.60 Re1/2 Pr1/3

Correlaciòn de Whitaker
N = hD/k

u voo

Re Pr2/3 propiedades físicas evaluadas a la

temperatura media del fluido

Para líquidos y 1<Re<200.000 Nu = 2,0+0,60 Re1/2 Pr1/3

Para gases y 1<Re<20 St = 2,2 Re-1 + 0,48 Re-1/2 propiedades físicas evaluadas
Para gases y 20<Re<150.000 Nu = 0,37 Re0,6 Pr1/3 a la “temperatura de película”
(T0+ Too)/2
Para líquidos y 1<Re<70.000 Nu = 2,0 + 0,6 Re1/2 Pr1/3
Para aire y 20<Re<150000 Nu = 0,31 Re0,6

31#
 Transmisión de calor por convección forzada para flujo normal a un tubo

Nu = hD/k

Re

Douglas y Churchill Nu = 0,46 Re1/2 + 0,00128 Re Para aire y Re>500

Hsu Nu = 0,43 + 0,48 Re1/2 Para aire y Re<500

Para otro líquido se multiplica por (1,1 Pr1/3)

Simplificadas para Aire: Nu = 0,245 Re0,6 si Re>103

Nu = 0,49 Re1/2 si Re<103

32#
 Transmisión de calor por convección forzada para flujo normal a un tubo

Para flujo cruzado sobre un cilindro

Para flujo de aire cruzado sobre un cilindro

33#
Transmisión de calor por convección forzada para distintas configuraciones

Aire:
NuDf = h D / kf = Gg ReDf
n

Líquidos:
Nu Df = h D / kf = Gl Prf
ReDf
1/3 n

Las propiedades físicas están evaluadas

a la “temperatura de película”

(T0 + T∞ ) / 2

34#
 Convección haces de tubos

Correlación de Dittus-Boelter

calentamiento enfriamiento

Correlación para bancada v

de tubos

Alineados
Alternados

v v

35#
 TRANSFERENCIA DE ENERGÍA Y MOMENTO
ANALOGÍAS

ANALOGÍA DE REYNOLDS
• Pr = 1
h f • No hay arrastre de forma
St =
 v Cp = 2

EJEMPLOS DE APLICACIÓN:
• Flujos paralelos a superficies planas
• Flujo en conductos

ANALOGÍA DE DEISSLER
• Pr >200
• Flujo altamente turbulento
St = f 1/2
0,0789 =
Pr3/4

36#
 ANALOGÍA DE COLBURN

f • 0,5<Pr<50
jH = St Pr2/3 =
2 • Re>10.000
• No hay arrastre de forma
• Tubos lisos largos

• Las analogías de Reynolds y Colburn son coincidentes para Pr = 1

Nu Re-1Pr-1/3 (b/0)-0,14

Re

37#
 CONVECCIÓN LIBRE
La convección libre sigue el mismo modelo que la convección forzada

Nu = C (Re)b Prd

a excepción de que la velocidad ahora depende de la fuerza de flotación

= . g . V
Fflotación

experimentada por las capas de fluido que están junto a la superficie caliente y que son menos densas
que las porciones frías.

Este efecto puede expresarse a través del coeficiente de expansión volumétrico  que representa la
variación de la densidad del fluido en función de T a P cte.

= - (1/) (d/dT)P

A mayor  mayor variación de volumen con T

◆ Para gases ideales  a una temperatura T es
gas ideal = 1/T
◆ Para líquidos  puede considerarse
constante en un dado intervalo de T
líquidos = (
38#
El cambio de densidad es
    

El trabajo hecho por m3 para lograr un desplazamiento vertical L es

W = (m/V).g.L = .g.L =    g L

Si todo el trabajo se convierte en energía cinética

  g L =  v2

Siendo: Nu = C (Re);b Prd

Nu = C ( Lv2 2v=2 /

) b/2 ( C / k)
gL p d

Nu = C (2 L2   g L / ) b/2 ( Cp / k) d

Nu = C (2 L3   g / ) b/2 ( Cp / k) d

2 L3   g
Gr = Número de Grashof
2

39#
 CONVECCIÓN LIBRE: NÚMEROS ADIMENSIONALES
Variables involucradas:
M L T 
gravedad g o 1 -2 0 Variables independientes:
densidad  1 -3 0 0
viscosidad  1 -1 -1 0 D 0 1 0 0
conductividad térmica del fluido k 1 1 -3 -1 k 1 1 -3 -1 =/ 0
 1 -1 -1 0
longitud característica D 0 1 0 0 g 0 1 -2 0
calor específico a presión cte Cp 0 2 -2 -1
coeficiente convectivo h 1 0 -3 -1
coeficiente de expansión térmica  0 0 0 -1
diferencia de temperatura T 0 0 0 1

9 variables - 4 (rango de matriz dimensional) = 5 grupos 

◆  1 = Da . b. kc. gd. 
1 = La x (M L-1T-1)b x (M LT-3 -1)c x (L T-2)d x (M L-3)

M: b+c+1=0 a=3/2
L: a-b+c+d-3= b=1  1 = D3/2  g1/2
0 T :-b - 3 c - 2 d = 0 c=0
d=1/2 
 : -c = 0
40#
◆  2 = De. f . kg. gh. Cp
2 = Le x (M L-1T-1)f x (M LT-3 -1)g x (L T-2)h x (L2 T-2 -1)
M: f+g=0
L: e -f+g+h+2= e=0
f=1  Cp
0 T :-f - 3g - 2h - 2 = 0 2= = Pr
g = -1 k
 : -g - 1 = 0 h=0
◆  3 = Di. j . kk. gm. 
2 = Li x (M L-1T-1)j x (M LT-3 -1)k x (L T-2)m x (T-1) M
: j+ k = 0
L: i-j+k+m=0 i = 1/2  D1/2
T : -j - 3k - 2m - 1 = 0 j=0
k=0 3=
g1/2
 : -k = 0
m = -1/2
◆  4 = Dn. ñ . ko. gp. T
4 = Ln x (M L-1T-1)ñ x (M LT-3 -1)o x (L T-2)p x (T)
M: ñ+o= n = -1/2
T g1/2
0L : n - ñ + o + p = 0 ñ=0
o=0 4=
T : -ñ – 3o – 2p + 1 = 0 D1/2
 : -o = 0 p = 1/2
41#
◆  5 = Dq. r . ks. gt. 
5 = Lq x (M L-1T-1)r x (M LT-3 -1)s x (L T-2)t x (M T-1 1)

M: r+s+1=0 q=1
L: q-r+s+t=0 r=0 hD
T : -r – 3s – 2t - 3 = 0 s=1 5= k = Nu

 : -s - 1 = 0 t=0

◆ El 12 3 4 = Gr (Número de Grashof)

producto
 D1/2 T g1/2
= D  g  T
3 
D3 g g1/2 D1/2

◆ Nu = f (Gr, Pr)

42#
 CONVECCIÓN LIBRE
CORRELACIONES EMPÍRICAS PARA DISTINTAS GEOMETRÍAS
Convección libre en placas horizontales calientes

Superficie caliente hacia arriba Superficie caliente hacia abajo

Nu = 0.27(Gr.Pr)1/4
Nu = 0.54(Gr.Pr)1/4 2.6x104< Ra <107
3x105< Ra <3x1010
Nu = 0.15(Gr.Pr)1/3 107< Ra <1011

Convección libre en placas verticales calientes

log Nu

Nu = 0.1(Gr.Pr)1/4 104< Ra <109

Nu = 0.1(Gr.Pr)1/3 1 09< Ra
<3x1012
Las propiedades físicas están evaluadas a
log Gr Pr la “temperatura de película” (T0 + Too)/2

43#
 Convección libre en cilindros verticales pequeños

log Nu

log Gr Pr D/L

Convección libre en cilindros horizontales

log Nu

Nu = 0.4 0< Ra < 10-5

Nu = 0.53(Gr.Pr)1/4 104< Ra <109

Nu = 0.13(Gr.Pr)1/3 109< Ra <1012

log Gr Pr

Correlación de Churchill y Chu

10-5< Ra <1012
Propiedades evaluadas a Tf
Para gases  = 1/Tf

44#
Convección libre en esferas

Nu = 0.60 (Gr . Pr) ¼

Expresiones simplificadas para convección libre al aire

h = C (T/L)n

45#
Correlaciones válidas para convección libre distintas geometrías

Ver fig. anterior

Nu = C (Gr Pr) a

Las propiedades físicas están evaluadas

a la “temperatura de película”
Ver fig. anterior
(T0 + T∞ )/2

Ver fig. anterior

46#
 EMISIÓN CALORÍFICA DE METALES FUNDIDOS
CORRELACIONES EMPÍRICAS

❖ Metales que pueden emplearse:

❖ Na, K, Li, Cs, Hg, Bi
❖ aleaciones de Bi con Pb, etc
❖ Características:
❖ Poseen alta conductividad térmica
❖ Pr << 1
❖ ¿Cuándo utilizarlos?
❖ cuando se requiere una intensa evacuación de calor
❖ cuando se necesita una alta temperatura del fluido a presión baja
❖ Cuidados a tomar en cuenta
❖ Los conductos de circulación deben ser herméticos y estar llenos con un gas neutro
porque en contacto con el aire puede ocurrir una oxidación enérgica.
❖ El depósito de capas de óxido disminuye la transferencia de calor

47#
Convección libre

102<Gr<109 C=0,52 n=0,25

Nu = C Grn Pr0,4 109<Gr<1013 C=0,106
n=0,33 Las propiedades físicas están
evaluadas a la “temperatura de
película” (T0 + Too)/2

Convección Forzada

tubos limpios
Nu = 4,36 + 0,025 Pe0,8
20<Pe<104
L/D>30
corrección para L/D<30
 = 1,7(D/L)0,16
tubos oxidados, flujo turbulento
Nu = 3,3 + 0,014 Pe0,8

Las propiedades físicas están

evaluadas a la temperatura
media del fluido

48#
 FLUJO DE CALOR EN UN ESPACIO CERRADO LIMITADO
Las condiciones del movimiento dependen de:
➢ la forma y dimensiones geométricas
➢ las características del fluido
➢ la diferencia de temperatura entre los medios

Canales o hendiduras verticales Canales o hendiduras horizontales

d d
Tp1
Tp2
Tp2 d

Tp2

d
Tp1 Tp1

Tp1

49#
 Capas intermedias cilíndricas o esféricas

Tp2
Tp2
Tp1

Tp1
Tp1
Tp2

• Coeficiente equivalente de conductividad térmica: keq = Q d / A

T
• Coeficiente de convección: c = keq / k siendo c = f ( Grf , Prf )
c = 1 Grf Prf < 1000 La transferencia solo se debe a kfluido
c = 0,105 (Grf Prf)0,3 03 < Grf Prf < 106
c = 0,40 (Grf Prf)0,2 06 < Grf Prf < 1010

Para cálculos aproximados: c = 0,18 (Grf Prf)0,25 Grf Prf > 103

50#
 Otras correlaciones definen un hm promedio de transferencia de calor:
Propiedades del
q = hm A (T1 - T2)
fluido a (T1+T2)/2

1) ESPACIOS CERRADOS HORIZONTALES

Para Aire - Placa superior caliente

q = hm A (T1 - T2) = kdA (T1 - T2)

hd
Nu = =1
k No hay movimiento del fluido, la
transferencia es por conducción

Para Aire - Placa inferior

caliente
Nu = 1 Gr < ~2000 conducción > convección
Nu = 0,195 Gr 1/4 04 < Gr < 4x105 convección ~
Nu = 0,068 Gr1/3 Gr > 4x105 conducción convección
> conducción
51#
2) ESPACIOS CERRADOS VERTICALES

Capas verticales de aire con L/d > 3

Nu = 1 Gr < 2000
Nu = 0,18 Gr 1/4 (L/d)-1/9 2x104 < Gr < 2x105
Nu = 0,065 Gr1/3 (L/d)-1/9 2x105 < Gr < 1,1x107

Capas verticales de
líquido
04 < Gr Pr < 107
Nu = 0,42 (Gr Pr)1/4 Pr0,012 (L/d)-0,3  < Pr < 2x104
0 < L/d < 40

106 < Gr Pr < 109

Nu = 0,046 (Gr Pr)1/3
1 < Pr < 20
1 < L/d < 40

52#
 CONVECCIÓN FORZADA Y LIBRE COMBINADAS

Gr Fuerza de empuje Gr/Re2 < 0.1 predominio de convección forzada

= Gr/Re2 > 10 predominio de convección libre
Re2 Fuerza de
0.1 < Gr/Re2 < 10 convección forzada y libre
inercia
combinadas
Flujo asistido: El efecto de empuje tiene la misma dirección que el movimiento forzado, los
efectos de convección natural y forzada se suman mejorando la transmisión de calor.
Flujo opuesto: El efecto del empuje y el movimiento forzado tienen direcciones opuestas, los
efectos de convección natural empeoran la transmisión de calor.
Flujo transversal: El efecto del empuje es perpendicular al movimiento forzado. Ayuda al
mezclado del fluido mejorando la transferencia de calor.

Circulación transversal en un tubo horizontal

producida por la convección libre del fluido
a) calentamiento del líquido
b) enfriamiento del líquido

a) b)

53#
Regímenes de flujo en tubos verticales

54#
Regímenes de flujo en tubos horizontales

Longitud
característica: D
T entre la T de
pared y T
media
volumétrica del
fluido

Nr de Graetz
Gz = Re.Pr.D/L

55#
COEFICIENTES DE TRANSMISIÓN DE CALOR PARA CONDENSACIÓN
DE VAPORES PUROS

56#
El proceso de condensación:
➢El vapor fluye sobre la superficie debido a un pequeño gradiente de presión.
➢Algunas moléculas de vapor chocan contra la fase líquida y rebotan.
➢Otras penetran y ceden su calor latente de condensación
➢El calor fluye a través del condensado hasta alcanzar la pared refrigerada
➢A su vez el condensado escurre por gravedad por la pared
➢El condensado es la única resistencia importante para la transmisión de calor
➢El condensado forma una película continua si la superficie está limpia, condensación en película
➢Superficies con impurezas grasas ocasionan una condensación en gotas, con mayores velocidades
de transmisión de calor, pero es difícil de mantener este régimen

Q = hm A (Td - Ts)

hm: Coeficiente de transferencia de calor para la película del líquido

A: Area de la superficie sólida
Td: T de rocío del vapor, es la del líquido en la interfase L-V
Ts: T constante de la superficie refrigerante
57#
 Expresiones para calcular hm para flujo laminar sin ondulaciones
1/3
1) Tubos horizontales

hm = 0,954 (
k3 2 g L
 )
k3 g2 
)
1/4

m
h = 0,725 (  D (Td-Ts)

2) Tubos o paredes verticales 1/3

Buena aproximación hm = 4/3 (

k3 2 g  D
3
)
para tubos

( k3 g2 
)
1/4
verticales cortos hm = 2/3 2
(L<15 cm) Para  L (Td-Ts)
L<24 m puede
haber hasta un 70%
d e e ro r
Expresiones p a r a c alcular hm para flujo turbulento (Re >
350)
1) Tubos o paredes verticales 1/2

hm = 0,003 ( k
3 2g (T -T
d )L

s
)
k3 2g 
)

1/4

hm = 0,021 (   D
58#
 EN RESUMEN

❖El punto de partida en convección consiste en calcular el coeficiente de película h.

❖Para el cálculo de h se utilizan correlaciones empíricas. Sólo pocos casos permiten una
solución analítica.
❖ El coeficiente películar h se calcula a partir del número de Nusselt
Nu = ( h L ) / k € h = ( Nu k ) / L .
❖ El cálculo del Nusselt depende si estamos en presencia de:
❖ Convección forzada: Nu = f ( Re, Pr )
❖ Convección natural: Nu = f ( Gr, Pr ) = f ( Ra ), siendo Ra = Gr Pr.
❖ y de la configuración física del sistema
❖ Flujo interno
❖ Flujo Externo
❖ Y se selecciona la correlación/gráfico correspondiente para poder calcular el Nu.
❖ A partir del Un Nusselt se calcula el coeficiente de película: h = ( Nu k ) / L
❖Conocido este coeficiente, se puede calcula el flujo de calor mediante la Ley de
Enfriamiento de Newton: Q = h A ( Tt - Ts ).

59#
 BALANCE MACROSCÓPICO DE ENERGÍA
ENTRADA
“1” DE CALOR
“2”
Q

SALIDA DE
TRABAJO
MECÁNICO
W

• Los planos “1”y”2” son perpendiculares a las

paredes del tubo.
• En los planos “1” y “2” se cumple:
Velocidad de tiempo ajustado paralela a las paredes
Densidad y demás propiedades físicas constantes en
toda la sección
Velocidad de Velocidad de Velocidad Velocidad
incremento entrada-salida neta de adición neta de trabajo
de energía
de energía de calor producida
= interna + -
interna cinética y al sistema por
cinética y potencial procedente el sistema
potencial por de los sobre los
convección alrededores alrededores
60#
60
#
 ^ ^

d (Utot + Ktot + tot) =  <v1> Û1 S1 -  <v2> Û2 S2 + 1/2  <v1 > S1 - 1/2  <v2 > S2
3 3
dt
+  <v1> 1 S1 -  <v2> 2 S2 + Q – W + p <v1> S1 - p <v2> S2

W incluye el trabajo transmitido mediante partes móviles (turbinas, compresor,etc.)

Q incluye toda la energía calórica que llega al fluido a través de las superficies sólidas por conducción o radiación
Definiendo:
Etot = (Utot + Ktot + tot)
1 =  <v1> S1 y  =  <v2> S2

p <v3> ^ BALANCE MACROSCÓPICO

d tot = -  [ (Û + DE ENERGÍA
d tE  + 1/2 <v> +  ) ] + Q - W
EN ESTADO NO
ESTACIONARIO
Es el Primer Principio de la Termodinámica aplicado a un sistema de flujo

p <v3> ^ BALANCE MACROSCÓPICO

^
 (Û +  + 1/2 <v> +  ) = ^ Q- DE ENERGÍA
W EN ESTADO ESTACIONARIO

^Q:Calor suministrado por unidad de masa de fluido que circula a través del
p ^
sistema
Û + W ^
= H:Trabajo producido
es la entalpía por unidad
por unidad de masa de fluido al pasar de “1” a “2”
de masa

2012 L.Colombo-FT-Mec.Transf.Calor 61#

61
#
 Casos particulares
T2 T2
1) Gases ideales: R   dT
^ ^
p=RT y Cp - Cv = 
H = Cp dT = M
R T1
T1
2) Líquidos incompresibles: T2

^ ^
 = cte y Cp = Cv 1
Cp dT +  (p - p )
2 1
H =
T
3) Energía interna, cinética y del trabajo despreciables
1

^ ^
H = Q
4) Simplificaciones
1/2 <v3>/ <v> = 1/2 <v2> Flujo turbulento, energía cinética pequeña
  = g  h, siendo  h diferencia de altura entre planos “1” y “2”

62#
RESUMEN DE BALANCES MACR OSCÓPICOS PARA SISTEMAS ISOTÉRMICOS
DE FLUJO CONTENIE NDO UNA SOLA ESPECIE QUÍMICA

63#