5.1 Unidad 5. Convección Forzada (Placas Planas)
5.1 Unidad 5. Convección Forzada (Placas Planas)
Contenido
5.1.1 Convección forzada en una placa plana
5.1.2 Análisis integral de la capa límite en una placa plana (método
de von Kárman)
a) Ecuación integral de continuidad
b) Ecuación integral de cantidad de movimiento
c) Ecuación integral de energía
5.1.3 Convección forzada en una placa plana en régimen laminar
a) Perfil de velocidad y espesor de la capa límite
hidrodinámica
b) Perfil de temperatura y espesor de la capa límite térmica
c) El coeficiente de transferencia de calor
5.1.4 Analogía entre la transferencia de calor y la fricción
5.1.5 Convección forzada en una placa plana en régimen
turbulento
Bibliografía
Manrique V., J. A., Transferencia de Calor, 2ª ed, Oxford University Press, México,
2002.
Ҫengel, Yunus A., Transferencia de Calor, McGraw-Hill, México.
Incropera, F. P., DeWitt, D. P., Fundamentos de Transferencia de Calor, Prentice-
Hall, México.
1
5.1 Convección forzada
La convección es un proceso de transporte de energía que resulta como
consecuencia del movimiento de un fluido.
Las aplicaciones de la convección en la ingeniería atañen, en general, a la
transferencia de calor entre un fluido en movimiento y la superficie de un sólido con
una temperatura diferente de la del fluido.
En la convección forzada el movimiento de un fluido es provocado por algún agente
externo al proceso. Por ejemplo, una bomba, un ventilador, el viento.
El flujo de calor por convección entre la superficie de un cuerpo y el fluido que pasa
sobre ella se calcula mediante la ley de Newton de Enfriamiento. Esto es,
q hA(Ts T )
h f ( , c p , k , , u, L)
2
ciento de la velocidad de corriente libre u . Este concepto de capa límite fue
introducido por L. Prandtl y básicamente divide el campo de flujo en dos regiones:
una capa muy delgada en donde las fuerzas viscosas de corte son importantes y
una capa exterior a esta capa en donde los efectos viscosos son despreciables
( 0 ). El fluido que se mueve en la capa límite hidrodinámica experimenta
gradientes de velocidad y, por lo tanto, está sometido a esfuerzos viscosos; en el
borde de la capa límite ( y ) la velocidad del fluido se hace prácticamente igual a
la velocidad de corriente libre y el gradiente de velocidad se hace cero, por lo que
el esfuerzo viscoso prácticamente desaparece.
3
Inicialmente el desarrollo de la capa límite hidrodinámica es laminar. Es decir, el
fluido se desplaza a lo largo de láminas y las partículas de éste siguen una
secuencia ordenada y continua sin cruzarse entre sí. Sin embargo, a cierta distancia
crítica xcrit que depende del campo de flujo y de las propiedades del fluido, se
empiezan a experimentar pequeñas perturbaciones que se amplifican a medida que
la distancia aumenta. Como consecuencia de este fenómeno se presenta un
proceso de transición hasta que el flujo se hace completamente turbulento. Esta
transición de régimen laminar a régimen turbulento no es abrupta y depende de las
condiciones de rugosidad de la superficie y de los niveles de turbulencia en la
corriente libre del fluido. Para determinar si el flujo es laminar o turbulento se usa el
número de Reynolds, que constituye un cociente entre fuerzas inerciales y fuerzas
de corte. El número de Reynolds evaluado en la posición x , Reynolds local, es
u x
Re x = Número de Reynolds local =
(5.1)
Usualmente para una placa plana se establece como criterio de transición que
u xcrit (5.2)
Rec = Número de Reynolds crítico = 5x105
4
(a)
(b)
Figura 5.3. (a) Diferentes regímenes de flujo en la capa límite sobre una placa
plana. (b) Perfil de velocidad en la capa límite turbulenta.
Dada la relación física que existe entre el movimiento del fluido y la transferencia de
energía que tiene lugar, la determinación analítica del coeficiente de transferencia
de calor implica un conocimiento completo de la distribución de velocidad del fluido
en la capa límite hidrodinámica y de la distribución de temperatura del fluido en la
capa límite térmica.
5
5.1.2 Análisis integral de la capa límite en una placa plana (método de von
Kárman)
La dificultad matemática para resolver las ecuaciones de cambio (ecuación de
continuidad, ecuación de cantidad de movimiento, ecuación de energía)
diferenciales para obtener soluciones exactas conduce a la búsqueda de métodos
aproximados que conduzcan rápidamente a una solución, aunque la exactitud de
ésta sea inferior. Un método que conduce a soluciones aproximadas es el método
integral de von Kárman. Consideraciones: el flujo en estado estable e
incompresible, no hay gradientes de presión, las propiedades del fluido son
constantes.
a) Ecuación integral de continuidad
En la figura 5.4 se muestra un volumen de control dentro de la capa límite
hidrodinámica.
0
uzdy
uzdy
x
0 x x
m 0 (5.3)
6
donde m denota el flujo másico que entra por la parte superior del volumen de
control.
0
u 2zdy u zdy m u
x x
0
2
x
s xz (5.4)
s xz
7
0
u 2dy
x x
0
u 2dy
x 0
uu dy uu dy
x x
0
x
s
x x
d d
dx dx
u 2
dy uu dy s
d
dx ( uu u 2 )dy s
d u u u u
dx u u u u
u2 dy s
d u u
dx 0 u u
u2 1 dy s (5.5)
8
Figura 5.6 Volumen de control en la capa límite térmica, balance de energía.
0
uH Zdy uH Zdy
x 0
x x
mT H qs 0 (5.6)
mT uzdy uzdy
0 x x
0 x
(5.7)
Combinando las ecuaciones 5.6 y 5.7, dividiendo entre xZ , y haciendo tender
x a cero, se obtiene que,
d T
dx 0
u ( H H )dy k
y y o
9
Puesto que ( H H ) c p (T T ) ,
d T
dx 0
u (T T )dy
y y o
(5.8)
k
donde difusividad térmica
cp
d
dx 0
u ( )dy
y y o
(5.9)
u 1 y y 2
u b 2c 2 3d 3
y
10
2u 1 y
u 2c 2 6d 3
y
2
Por otra parte, la condición de continuidad en la velocidad al pasar del perfil dentro
de la capa límite a la velocidad de corriente libre requiere que,
u u en y (5.12)
2u
0 en y0 (5.14)
y 2
Sustituyendo las condiciones 5.11, 5.12, 5.13 y 5.14 en la ecuación 5.10 (perfil de
velocidad supuesto) se obtiene que a = 0, b = 3/2, c = 0, y d = -1/2. Por lo tanto,
3
u 3 y 1 y
(5.15)
u 2 2
du 3 u
De la ecuación 5.15 el gradiente de velocidad en la superficie es
dy y 0 2
11
du 3 u
s (5.16)
dy y 0 2
d 39 3 3
u
dx 280 2 2
Separando variables,
d 140
(5.17)
dx 13 u
Integrando,
2 140 x
C1
2 13 u
4.64 x / u (5.18)
4.64
(5.19)
x Re x
donde Re x u x / u x /
12
Con las ecuaciones 5.15, perfil de velocidad, y 5.19, espesor de la capa límite
hidrodinámica, queda definida completamente la distribución de velocidad en la
capa límite hidrodinámica para el flujo laminar.
T Ts
2 3
y y y
A B C D (5.20)
T Ts
T Ts , 0 en y0 (5.21)
T T , en y (5.22)
T / y 0 , / y 0 en y (5.23)
2T / y 2 0 , / y 2 0
2
en y0 (5.24)
Sustituyendo las condiciones 5.21, 5.22, 5.23 y 5.24 en la ecuación 5.20 se obtiene
que: A = 0, B = 3/2, C = 0 y D = -1/2. Por lo tanto,
T Ts 3 y 1 y
3
T Ts 2 2
(5.25)
d 3
dy 2 (5.25a)
y 0
13
Para calcular el espesor de la capa límite térmica es necesario sustituir la
distribución de velocidades (ecuación 5.15), la distribución de temperaturas
(ecuación 5.25) y el gradiente de temperatura en Y = 0 (ecuación 5.25a) en la
ecuación integral de energía (ecuación 5.9). Entonces, al hacer las sustituciones e
integrando, se obtiene que
d 3 2 3 4 3
u
dx 20 280 2
3 d 3
u 2
20 dx 2
Al derivar, queda
1 d d
u 2 2
10 dx dx
2 2 d 3 d
2 10
dx dx u
14
d 140
Pero, del análisis de la capa límite hidrodinámica se tiene que y
dx 13 u
280 x
2 , entonces,
13 u
280 x 2 d 140
2 3 10
13 u dx 13 u u
3 2 d 13
+4x dx 14
d 1 d 3
Pero viendo que 2 , se tiene la ecuación diferencial para
dx 3 dx
3 4 d 3 13
+ 3 x dx 14 (5.27)
13 x0
3/4
3
1 (5.28)
14 x
15
1/3
donde
/ cp
Pr Número de Prandtl =
k / c p k
T
hx A(Ts T ) kA
y y o
Entonces,
k T / y ) y 0
hx (5.30)
Ts T
16
donde hx denota el coeficiente local de transferencia de calor. El gradiente de
temperatura en la interfase ( y 0) está dado por la ecuación 5.25a, entonces,
k T / y ) y 0 3k 3k
hx (5.31)
Ts T 2(Ts T ) 2
3 k
hx 1/3
2 (0.9756Pr )(4.64 x Re1/2 )
k
hx 0.332 Re1/2
x Pr
1/3
(5.31a)
x
17
O, en forma adimensional,
hx x
0.332Re1/2
x Pr
1/3
(5.31b)
k
hx x
Nux 0.332(Re x )1/2 (Pr)1/3 (5.32)
k
1/2 1/2
1 L 1 L ku x k u L
h hx dx 0.332 Pr1/3 dx (0.332) Pr1/3 x 1/2 dx
L 0 L 0 x L 0
k u L
1/2
h 2 (0.332) Pr 2hx L
1/3
L
Por lo tanto,
hL
Nu 0.664Re1/2
L Pr
1/3
(5.33)
k
18
mayoría de los líquidos y gases caen dentro de esta categoría. Los metales líquidos
deben tratarse de manera diferente porque tienen números de Pr del orden de 0.01.
Dado que en el análisis integral se supuso que las propiedades del fluido son
constantes, es recomendable evaluar las propiedades del fluido a la temperatura de
película Tf, que es el promedio aritmético de la temperatura de la superficie y la del
fluido. Es decir,
Ts T
Tf
2
Ejemplo 5.1
Considere una placa de 0.1 m de longitud y 0.1 m de ancho que se encuentra a una
temperatura de 80 °C. Se hace pasar agua sobre su superficie a una velocidad de
0.10 m/s y con una temperatura de 40 °C. Calcule: a) el calor disipado y, b) el
espesor de la capa límite hidrodinámica y el espesor de la capa límite térmica en
x = 0.1 m.
Ts T 80 40
Temperatura de película: T f 60 C
2 2
Determinación de las propiedades del agua a la temperatura de 60 °C:
19
El número de Reynolds es menor que 5x105, entonces el flujo es laminar y puede
emplearse la ecuación 5.33. Según esta relación,
hL
Nu 0.664Re1/2
L Pr
1/3
0.664(2.09 104 )1/2 (3.02)1/3 138.70
k
Por consiguiente,
El flujo de calor es
4.64 x 4.64(0.1 m)
3.21 103 m = 3.21 mm
Re x 2.09 10 4
0.9756Pr 1/3 0.9756(3.02) 1/3 0.67
20
0.67(3.21 mm) 2.15 mm
u 3 u u
0.323 Re1x/ 2 (3)
y y 0 2 x
u u
s 0.323 Re1x/ 2 (4)
y y 0 x
c fx s
0.323Rx1/ 2 (5)
2 u2
21
Coeficiente local de transferencia de calor y número de Nusselt local:
hx x
Nu x 0.332 Re1x/ 2 Pr1 / 3 (6)
k
Nux hx x / k hx
St x (7)
Re x Pr ( u x / )(c p / k ) uc p
Nux
St x 0.332Rex1/2 Pr 2/3 (7)
Re x Pr
de donde se obtiene
C fx
St x Pr 2 / 3 (9)
2
22
5.1.5 Convección forzada en una placa plana en régimen turbulento
En el flujo paralelo sobre una placa plana de longitud L, figura 5.7, el desarrollo de
la capa límite laminar comienza desde x 0 y la transición al flujo turbulento ocurre
a una distancia x xc , donde se alcanza un número de Reynolds crítico
Rec u xc / . Desde este punto y hasta x L se desarrolla la capa límite
turbulenta.
Figura 5.7 Placa plana con capa límite laminar y capa límite
turbulenta.
C fx 0.0592Re x1 / 5
(2)
23
Aplicando la analogía de Reynolds del punto anterior, se obtiene que
C fx
St x Pr 2 / 3 0.0296 Re x1 / 5
2
Nu x
Pero St x , entonces
Re x Pr
hx x
Nu x 0.0296 Re x4 / 5 Pr 1 / 3
k
k
hx 0.0296 Re x4 / 5 Pr 1 / 3 (3)
x
1 xc
0 hx laminardx x hx turbulentodx
L
h (4)
L c
1 xc k u x k u x
1/ 2 4/5
h 0.332 Pr 1 / 3 dx
L
Pr 1/ 3
dx 0.0296
L 0 x xc x
24
k u
1/ 2
1 / 2
u
4/5
h Pr1 / 3 0.332
xc L
xc
1 / 5
x dx (0.0296) x dx
L 0
k u
1/ 2
u
4/5
5 4 / 5 5 4 / 5
h Pr1 / 3 0.332 (2 x 1/ 2
) (0.0296) L xc
L
c
4 4
k u x
1/ 2
u L
4/5
u x
4/5
h Pr1 / 3 0.664 c 0.037 0.037 c
L
h
k
L
Pr
1/3
1/2 4/5 4/5
0.664 Rec 0.037 ReL 0.037 Rec (5)
Nu
hL
k
Pr1/ 3 0.037 ReL4 / 5 871 (6)
Si la transición de flujo laminar a flujo turbulento ocurre en x 0 , (no hay capa límite
laminar) el número de Reynolds crítico se hace cero y, de la ecuación 5, el número
de Nusselt promedio es,
Nu
hL
k
Pr
1/3
4/5
0.037 ReL (7)
25
Para números de Reynolds entre 1 10 y 1 10 , el coeficiente local de fricción
7 9
El coeficiente de fricción promedio en una placa plana con una capa límite laminar
hasta Rec 5 105 y una capa límite turbulenta desde este punto y hasta Re L , se
puede calcular con la ecuación
0.455 1742
Cf 2.584
(log Re L ) Re L
Cf
St Pr 2/3 0.288(log Re) 2.584 871Re L1
2
Entonces, el número de Nusselt promedio es
Nu
hL
k
Pr1 / 3 0.288 Re L (log Re L ) 2.584 871 (8)
Ejemplo 5.2
Considere una placa plana de 1 m de longitud y 1 m de ancho que tiene una
temperatura superficial de 80 °C. Sobre la placa pasa agua a una velocidad de 1
m/s y con una temperatura de 40 °C. Considere que la transición de flujo laminar a
flujo turbulento ocurre cuando Rec 5 10 . a) Determine la coordenada xc donde
5
26
Ts T 80 40
Temperatura de película: T f 60 C
2 2
u xc
Rec
u L u L (1 m/s)(1 m)
Re L 2.09 106
0.478 x106 m2 / s
5x105 Re L 1x107
Nu
hL
k
Pr1/ 3 0.037 ReL4 / 5 871
hL
Nu 3.021/3 0.037(2.09 106 )4/5 871 4826.86
k
27
El coeficiente promedio de transferencia de calor es,
Nu k (4826.86)(0.651 W/mK)
h 3142.28 W/m2 K
L (1 m)
Por lo tanto,
q 125.69 kW
Ejemplo 5.3
Aire a 10°C fluye sobre una placa plana de 1.5 m de longitud y 1 m de profundidad
con una temperatura superficial de 90°C. Determine la velocidad del aire para que
el calor transferido de la placa hacia el aire sea 3 kW.
Ts T 900 C 100 C
Tf 500 C
2 2
Pr 0.703
28
q 3x103 W
h
25 W/m2 °C
A(Ts T ) (1.5 m )(90 10) C
2
c) Número de Reynolds
Entonces,
5/4 5/4
(Nu/Pr1/3 ) 871 (1339.29 / 0.7031/3 ) 871
Re L 1.0229 106
0.037 0.037
6 2.025 10
5
kg/m s
u Re L (1.0229 10 ) (1.095 kg/m3 )(1.5 m) 12.61 m/s
L
29