Transferencia de Calor
Transferencia de Calor
Transferencia de Calor
las mismas relaciones anteriores, excepto que todas las propiedades se Pm = 1,434 x 10V5 kg/m. s CO,0347 lb,/h. ft]
evalúan a una entalpía de referencia i* dada por
(0>0758)(918)(0,70) = 3 395 x 1o6
i * zz i
OO + 0,5(i, - i,) + 0,22(i,, - i,) (5.130) Re,, m =
1,434 x 10-5 ’
El número de Stanton se redefine como Así, se llega a la conclusión de que hay que tener en cuenta la transferen-
cia de calor, tanto en la capa límite laminar como en la capa límite
Sti = hi turbulenta. Se determinan primero las temperaturas de referencia en
(5.131) ambos regímenes y después se evalúan las propiedades a esas tempe-
PU,
raturas.
Este número de Stanton se utiliza en las Ecs. (5.125), (5.126) o (5.127)
para calcular el coeficiente de transferencia de calor. Cuando se calculan
las entalpías para utilizarlas en las relaciones anteriores, debe emplearse Zona laminar
la entalpía total; esto es, debe incluirse tanto la energía química de 7-1
disociación como la energía interna térmica. El método de la entalpía de T, = T, 1 + -M2, = (233) [l + (0,2)(3)‘] = 652 K
referencia ha resultado adecuado para calcular la transferencia de calor ( 2 >
a alta velocidad con una precisión mayor del 10 por 100.
Suponiendo un número de Prandtl alrededor de 0,7, se tiene
Parte turbulenta Éste es el coeficiente medio de transferencia de calor para la parte lami-
nar de la capa límite, y la transferencia de calor se calcula a partir de
Suponiendo Pr = 0,7 se tiene
q = h.4(Tp - T,,)
Tap - Tm Top - 233 = (56,26)(0,222)(35 - 311)
r = Pr”3 = 0,888 =
To - Tm = 652 - 233 = - 3.445 W [- 11.750 Btu/h]
Tpa = 605 K = 332°C de modo que se necesitan 3.445 W de refrigeración en la zona laminar de
la placa, por metro, en la dirección z.
T* = 233 + (0,5)(35 + 40) + (0,22)(605 - 233) = 352,3 K
La concordancia entre Pr* y el valor que se ha supuesto es suficiente- Para determinar la transferencia de calor turbulenta, se debe obtener una
mente buena. Las otras propiedades a utilizar en el análisis de la transfe- expresión para el coeficiente de transferencia de calor local, a partir de
rencia de calor turbulenta son St*
x Pr*‘13 = 0,0296 Re*-
x ‘Is
(LO132 x 105)(1/20) e integrar después, desde x = 0,222 m hasta x = 0,7 m, para determinar
p* = = 0,050l kg/m3 la transferencia de calor total
(287)(352,3)
- 115
k* = 0,0302 W/m “C c; = 1,009 kJ/kg . “C Introduciendo los valores numéricos de las propiedades se tiene
h, = 94,34x- “’
Transferencia de calor laminar El coeficiente medio de transferencia de calor en la zona turbulenta se
determina por
Se supone
h= s 072 2 hx dx
= 111,46 W/m’.“C [19,6 Btu/h.ft’.“F]
P*%X, s 09222 dx
Re&, = 5 x lo5 = ~
P Haciendo uso de este valor, se puede calcular la transferencia de calor en
la zona turbulenta de la placa plana
(5 x 105)G07 x lo-‘) = o 222 m
x, = q = hA(T, - T,,)
(0,0508)(918) ’
- = (111,46)(0,7 - 0,222)(35 - 332)
hx, = - 15.823 W [ - 54.006 Btu/h]
/2 py*1/3
Nu* =
k*
= 0,664 W&)’
La cantidad potencia total de refrigeración necesaria, es la suma de las
= (0,664)(5 x 105)“2(0,697)1’3 = 416,3 transferencias de calor en las porciones laminar y turbulenta
Refrigeración total = 3.445 + 15.823 = 19.268 W [65.761 Btu/h]
h = (41624 W’3)
= 56,25 W/m2.“C [9,91 Btu/h.ft’.“F] En estos cálculos se supone longitud de 1 m en la dirección z.
0,222
5. LOS PRINCIPIOS DE LA CONVECCIÓN 183
5.13. RESUMEN
La mayor parte d&&e cap&&o se ha ocupado del Rujo sobre placas 3, Establecer el regimen de fIujo según lo determina el n@ero
pIanas y de la transferencia de calor asociada. Para comodidad deI de Reynolds.
lector, se han agrf.qad~ las ewacicmes en 1% Tabla 5.2, junto con las 4 . Seleccionar la ecuaci& apropiada, teniendo en cti&t~ eI r&-
restricciones que eor&wm. El pr~irn~~~ gerkeräl es entonces: men de Ihjo y cualquier restricción en ias propiedades del
fluido que pudiera darse,
1. Evdnar hs propiedades del fltidu; &to será gereneralmente h 5 . Cakular @los valot,kalores del coefiicente de traaskrewia
la temperatura de. pelkxda de calor por convección y/o la transferewia de 6aIor.
2. Establecer las condiciones de contorno, esto- es, temperatura
comtante o flujo de &or constante.
174 TRANSFERENCIA DE CALOR
donde 6 es, como antes, el espesor de la capa límite. Cuando lo que se En el caso 2 se tiene la condición
desea es un análisis integral, la integral de la cantidad de movimiento
puede evaluarse con la Ec. (US), por ser la subcapa laminar tan delgada. en xcrit = 5 x lo5 .Z- (5.92)
* = 4aIn
Sin embargo, el esfuerzo cortante en la pared no puede calcularse con la u,
Ec. (5.88) porque da un valor infinito en y = 0.
Para determinar el espesor de la capa límite turbulenta, se emplea la Ahora, 6,,, se calcula a partir de la relación exacta de la Ec. (5.21~)
Ec. (5.17) para la relación integral de la cantidad de movimiento, y el
esfuerzo cortante en la pared se evalúa a partir de las relaciones empíri- 6 hlll = 5,Ox,,,,(5 x 105)-“2 (5.93)
cas para la fricción superficial presentadas anteriormente. De acuerdo
con la Ec. (5.52) La integración de la Ec. (5.89) da
1/5 5
* - Ln 4 (x 4/5 - x4/5
cnt 1 (5.94)
y entonces, para Re, < lo’, de la Ec. (5.77) se obtiene Combinando las distintas relaciones anteriores se tiene
6
- = 0,381 Re;‘j5 - 10.256 Re;’ (5.95)
X
(5.89)
Esta relación sólo es aplicable en la región 5 x lo5 < Re, < 10’.
Ahora, utilizando la ecuación integral de la cantidad de movimiento
para gradiente de presión nulo [Ec. (5.17)], junto con el perfil de veloci-
dades y el esfuerzo cortante en la pared, se obtiene EJEMPLO 5.9. ESPESOR DE LA CAPA LÍMITE TURBULENTA. Cal-
cúlese el espesor de la capa límite turbulenta al final de la placa del
Ejemplo 5.7, suponiendo que se desarrolla (a) desde el borde de ataque
2 J;[i - (;)1'7](;)1-ldy = 0,0296(&~'5 de la placa y (b) desde el punto de transición con Re,,,, = 5 x 105.
Solución. Puesto que ya se ha calculado el número de Reynolds como
Integrando y reduciendo términos se tiene Re, = 1,553 x 106, es cosa sencilla introducir este valor en las Ecs. (5.91)
y (5.95) junto con x = L = 0,75 m, para dar
1
Muchos trabajos experimentales han mostrado que el perfil de velocida-
Nu, = 1,04 Nu, _ (5.87) des de una capa límite turbulenta, fuera de la subcapa laminar, se puede
_1 1 y = coos, describir por una relación con una potencia de un séptimo
h = (CW (OW‘W
x = 4,349 W/m2 . “C La diferencia de temperaturas en el borde de salida de la placa
0,4 (x = L = 0,6 m), se obtiene de las Ecs. (5.48) y (5.50) con la constante
0,453, que dan
h = (2)(4,349) = 8,698 W/m2. “C [1,53 Btu/h. ft2. ‘F]
4 = (8,698) (0,4)(60 - 27) = 114,8 W [392 Btu/h] (T, - TJxzL = (2437;);@:795) = 365,4"C
EJEMPLO 5.5. PLACA PLANA CON FLUJO DE CALOR CONSTANTE.
Se construye un calentador de 1,O kW con una placa de vidrio que lleva Una solución alternativa se basaría en el número de Nusselt, Ec. (5.51).
una película conductora de la electricidad que proporciona un flujo de
calor constante. La placa es de 60 por 60 cm y está colocada en una EJEMPLO 5.6. CORRIENTE DE ACEITE SOBRE UNA PLACA PLANA
corriente de aire a 27 T, 1 atm y u, = 5 m/s. Calcúlese el promedio de CON CALEFACCIÓN. Sobre una placa cuadrada de 20 cm de lado, se
la diferencia de temperaturas en la placa y la diferencia de temperaturas obliga a moverse aceite de motor a 20°C a una velocidad de 1,2 mjs. La
en el borde de salida. placa se calienta hasta una temperatura uniforme de 60°C. Calcúlese el
calor perdido por la placa.
Solución. Las propiedades se deberían evaluar a la temperatura de
película, pero no se conoce la temperatura de la placa, así que para un Sohción. Primero se evalúa la temperatura de película
cálculo inicial se toman las propiedades de la corriente libre en las con-
diciones de 20 + 60
TE-= 40 "C
T, = 27 "C = 300 K 2
v = 15,69 x 10m6 m’/s Pr = 0,708 k = 0,02624 W/m . “C Las propiedades del aceite de motor son
Como el número de Prandtl es tan grande, se va a emplear la Ec. (5.51) Empleando la distribución de velocidades dada por la Ec. (5.19) se tiene
para la solución. Se ve que h, varía con x de la misma manera que en la
Ec. (5.44) esto es, h, CC X- ll’, así que se tiene la misma solución que en 3 PU
--a
la Ec. (5.45) para el promedio del coeficiente de transferencia de calor. La zp 2 6
evaluación de la Ec. (5.51) en x = 0,2 da
y haciendo uso de la relación del espesor de la capa límite
(5.53)
q = hA(T, - T,) = (219,2)(0,2)‘(60 - 20) = 350,6 W La Ec. (5.44) se puede reescribir de la manera siguiente
Nu h
-2 = 2 = ()J32pr-2’3 Re,“2
Re, Pr PC,%
5.7. RELACIÓN ENTRE LA FRICCIÓN EN EL FLUIDO
Y LA TRANSFERENCIA DE CALOR El grupo de la izquierda recibe el nombre de número de Stanton
h
Ya se ha visto que los campos de temperatura y de velocidad están st, = --sL
relacionados, Se busca ahora una expresión mediante la que se pueda PC,U,
relacionar directamente la resistencia de fricción con la transferencia de de manera que
calor.
El esfuerzo cortante en la pared puede expresarse en función de un St x Pr2/3 = O,332Re-“’
x (5.55)
coeficiente de fricción C, Comparando las Ecs. (5.54) y (.5.55), se advierte que los miembros de la
1P
derecha son iguales, excepto por una diferencia de un 3 por 100 en la
= cf!+ (5.52) constante, resultado de la naturaleza aproximada del análisis integral de
ZP la capa límite. Se admite esta aproximación y se escribe
1
au La Ec. (5.56), llamada la analogia de Reynolds-Colburn, expresa la rela-
zp=$- p ción entre la fricción en el fluido y la transferencia de calor en el flujo
5. LOS PRINCIPIOS DE LA CONVECCIÓN 165
En el caso de flujo de calor constante, se puede demostrar que el En el caso de flujo de calor constante, se cambia 0,3387 por 0,4637 y
número de Nusselt local viene dado por 0,0468 se cambia por 0,0207. Las propiedades se siguen evaluando a la
temperatura de película.
hx
Nu x = - = 0 453 Re”’ Pr’j3 (5.48)
k ’ =
EJEMPLO 5.4. PLACA PLANA ISOTERMA CALENTADA EN TODA SU
que puede expresarse en función del flujo de calor de la pared y de la LONGITUD. Considérese que la placa plana de la corriente del Ejem-
diferencia de temperaturas como plo 5.3, se calienta en toda su longitud hasta una temperatura de 60°C.
Calcúlese el calor transferido en (a) los primeros 20 cm de la placa y (b)
los primeros 40 cm de la placa.
(5.49)
Solución. Se quiere obtener la transferencia de calor total en una deter-
minada longitud de la placa; así que se necesita calcular los coeficientes
El valor medio de la diferencia de temperaturas a lo largo de la placa, de transferencia de calor medios. Con este fin, se utilizan las Ecs. (5.44) y
con la condición de flujo de calor constante, se puede obtener efectuan- (5.45), evaluándose las propiedades a la temperatura de película
do la integración
27 + 60
L T, = ~ = 43,5”C = 316,5 K [110,3”F]
Tp-Tm=h (Tp-T,)dx=; oL&+ 2
s0 s x
Las propiedades, tomadas del Apéndice A, son
q&lk (5.50) Y = 17,36 x 10b6 m2/s Cl,87 x 10m4 ft2/s]
= 0,6795 Rei’ 2Pr 1’3
0 k = 0,02749 W/m . “C CO,0159 Btu/h ft . “F]
Pr = 0,7
En estas ecuaciones, qp representa el flujo de calor por unidad de área y cp = 1,006 kJ/kg . “C [0,24 Btu/lb, . ‘F]
tendrá unidades de vatio por metro cuadrado (W/m’) en el SI, o de En x = 20 cm
unidades térmicas británicas por hora por pie cuadrado (Btu/h . ft’) en el
sistema inglés. Re, = !f%? = (2) (W)
V 17,36 x lo-6 = 23’041
Otras relaciones
Nu, = % = 0,332 Re212 Pr’13
La Ec. (5.44) se aplica a fluidos con números de Prandtl entre 0,6 y 50
aproximadamente. No sería aplicable a fluidos con números de Prandtl = (0,332)(23.041)1’2(0,7)“3 = 44,74
muy bajos, como los metales líquidos, o a fluidos con números de
Prandtl altos, como los aceites pesados o las siliconas. Churchill y Ozoe h = Nu _ (44,74)(0,02W
[9] han efectuado la correlación de una gran cantidad de datos en un x
w
intervalo muy amplio de números de Prandtl, para dar la relación si-
guiente para el flujo laminar sobre una placa plana isoterma = 6,15 W/m* . “C [ 1,083 Btu/h . ft2 . ‘F]
Ésta se puede escribir en función del número de Reynolds como Para calcular el flujo másico que entra en la capa límite, proveniente de
la corriente libre, entre x = 20 cm y x = 40 cm, simplemente se toma la
6-=-
464 diferencia entre el fhrjo másico de la capa límite en esos dos valores de x.
(5.21) El flujo másico en la capa límite para cualquier valor de x viene dado
X ReXi
por la integral
donde
6
Re, = u,x pu dq’
V s0
La solución exacta de las ecuaciones de la capa límite según se dan en el donde la velocidad está dada por la Ec. (5.19)
Apéndice B conduce a
6 590
-=__ (5.21n)
x Rei”
Evaluando la integral con esta distribución de velocidades, se tiene
EJEMPLO 5.3. FLUJO MÁSICO Y ESPESOR DE LA CAPA LÍMITE.
Sobre una placa plana circula aire a 27 “C y 1 atm, y a una velocidad de
2 m/s. Calcúlese el espesor de la capa límite a distancias de 20 y 40 cm
del bórde de ataque de la placa. Calcúlese el flujo másico que entra en la
capa límite entre x = 20 cm y x = 40 cm. La viscosidad del aire a 27 “C Así, el flujo másico que entra en la capa límite es
es 185 x 10m5 kg/m. s. Supóngase la unidad de longitud en la direc-
ción 2.
Am = $PG(~,, - d,,)
Sohción. La densidad del aire se calcula a partir de
= (2, (1,177)(2,0)(0,0079 - 0,00559)
LO132 x lo5 = 3,399 x 1O-3 kg/s [7,48 x 1O-3 Ib,&]
= 1,177 kg/m3 ‘CO,073 lb,/ft3]
p = h= (287)(300)
300
j&=U’=p=
a2 460,7 ‘@l
Los términos de interés son, por tanto Los perfiles de temperaturas de la figura adjunta pueden ser algo confu-
sos. El vidrio tiene el gradiente de temperatura más pronunciado en la
base, pero su valor de k, mucho menor, produce un menor flujo de calor.
hP
Material Li m mL
FIGURA E. JI ZMPl.0 2.7
Cobre 12,99 3,604 0,3604
LO
Acero inoxidable 294,l 17,15 1,715 Cobre,- k -= 385
Vidrio 6.250 79,06 7,906 - -W/m
- -“C-
Estos valores se pueden sustituir en la Ec. (2.33~) para calcular las tem-
peraturas en diferentes localizaciones de x a lo largo de la barra, y los
resultados se muestran en la figura adjunta. Se advierte que el vidrio se
comporta como una aleta «muy larga», y su comportamiento podría
calcularse a partir de la Ec. (2.32). Los rendimientos de las aletas se
calculan con la Ec. (2.38), usando la aproximación de la longitud corre-
gida de la Ec. (2.42). Se tiene
L, = L + i = 10 + f = 10,5 cm
Material hPkA 6
Cobre 0,190 0,3784 2 4 6 8 10
Acero inoxidable 0,0084 1,8008 x, cm
Vidrio 3,9 x 1om4 8,302
Para comparar los flujos de calor se podría, tanto calcular los valores a EJEMPLO 2.8. ALETA DE ALUMINIO RECTA. Una aleta de aluminio
partir de la Ec. (2.36) para un valor unidad de B,,, como observar que el [k = 200 W/m . “C] de 3,0 mm de espesor y 7,5 cm de longitud sobresale
rendimiento de la aleta ofrece una comparación del flujo de calor relati- de una pared, igual que en la Figura 2.9. La base se mantiene a 3OO”C, y
vo, porque la transferencia de calor máxima es la misma en todos 10s la temperatura del ambiente es 50 “C, con h = 10 W/m2 . “C. Calcúlese la
casos; esto es, se está tratando con el mismo tamaño de aleta, perfil y pérdida de calor de la aleta por unidad de anchura del material.
valor de h. Se calculan así los valores de q, con la Ec. (2.38) y los valores Solución. Se puede emplear el método de resolución aproximado alar-
de mL, anteriores. gando la aleta una longitud ficticia t/2, y calculando a continuación la
transferencia de calor desde una aleta con el extremo aislado según viene
q relativo al cobre, dado por la Ec. (2.36). Se tiene
Material 4, porcentaje L, = L + t/2 = 7,5 + 0,15 = 7,65 cm
Cobre 0,955 100
Acero inoxidable 0,526 53,l
Vidrio 0,124 12,6
36 TRANSFERENCIA DE CALOR
cuando la anchura de la aleta z » t. Así, La transferencia de calor real es, entonces, el producto del flujo de calor
= 359 W/m [373,5 Btu/h. ft] Energía que entra por la cara izquierda + Energía generada
en el elemento =
Energía que sale por la cara derecha + Energía perdida por convección
EJEMPLO 2.9. ALETA DE ALUMINIO ANULAR. Para disipar calor,
se colocan unas aletas de aluminio de 1,5 cm de ancho y 1,0 mm de
espesor sobre un tubo de 2,5 cm de diámetro. La temperatura de la
superficie del tubo es 17O”C, y la temperatura del fluido ambiente es
25 “C. Calcúlese la pérdida de calor por aleta para h = 130 W/m’ . “C.
Supóngase k = 200 W/m. “C para el aluminio.
Solución. En este ejemplo, se puede calcular la transferencia de calor
utilizando las curvas de rendimiento de aletas de la Figura 2.12. Los
parámetros necesarios son
L, = L + t/2 = 1,5 + 0,05 = 1,55 cm
r1 = 2,5/2 = 1,25 cm
1
u2 130 li2
165
CAPÍTULO 3
TABLA 3-3
Eficiencia y áreas de superficie de configuraciones comunes de aletas
Aletas rectangulares rectas
tanh mLc
m 22h/kt hfin
haleta
mLc t
Lc L t/2
Aaleta 2wLc w
L
x
Aletas triangulares rectas y = (t/2) (1– x/L)
1 l1(2mL)
m 22h/kt hfin
haleta
mL l0(2mL)
fin 2w2L (t /2)
2 2
AA
aleta t
w
L