Frank Kreith C-4 PDF

CAPÍTULO 4
Análisis de
transferencia de calor
por convección
Flujo a alta velocidad sobre el
cuerpo de un misil genérico que
Conceptos y análisis que se deben aprender
muestra una onda de impacto La transferencia de calor por convección comprende dos mecanismos
oblicua en la punta del cuerpo, de ocurrencia simultánea, difusión o conducción, acompañada con
un embudo de expansión en el
transporte macroscópico de calor hacia (o desde) un fluido en movi-
hombro y una región de aire
miento o circulando. Este modo de transferencia de calor se encuentra
inerte debida a la interacción de
en una gran variedad de aplicaciones que incluyen, entre muchas
la capa límite laminar-onda
de impacto y características de otras, efectos de sensación térmica en el invierno, enfriamiento de
flujo base a alta velocidad en la toberas de cohetes, enfriamiento de láminas microelectrónicas, recu-
parte posterior. peración de calor de gases de la combustión en un intercambiador
Fuente: Impresa con permiso de Erdem de calor, enfriamiento de álabes de turbinas de gas, calentamiento de
E., Yang L. y Kontis K. Presentada agua en un colector solar y enfriamiento de una mezcla de agua-glicol
como “Drag Reduction by Energy en un radiador automotriz. Así pues, es necesario tener una buena
Deposition in Hypersonic Flows”, en
la 16a. AIAA/DLR/DGLR International
comprensión de los conceptos y expresiones matemáticas que descri-
Space Planes and Hypersonic Systems ben la transferencia de calor por convección para el diseño ingenieril
and Technologies Conference, Bremen, de esos sistemas y dispositivos. Estos conceptos y expresiones mate-
Alemania 2009. Afiliación: The School máticas se desarrollan en este capítulo y su estudio le enseñará:
of MACE, The University of Manchester.
• Cómo modelar una capa límite en transferencia de calor por con-
vección.
• Cómo deducir las ecuaciones matemáticas para la conservación
de la masa, cantidad de movimiento y energía térmica.
• Cómo efectuar análisis dimensionales y desarrollar correlacio-
nes para transferencia de calor por convección con fluidos
diferentes en flujo laminar y turbulento.
• Cómo obtener soluciones analíticas para las ecuaciones de capa
límite de flujo laminar comunes.
• Cómo aplicar la analogía entre la cantidad de movimiento y la
transferencia de calor para resolver problemas de convección de
flujo turbulento.
4.1 Introducción
En los capítulos anteriores, se consideró la convección sólo hasta el punto en que
proporciona la condición de frontera cuando la superficie de un cuerpo está en
contacto con un fluido a una temperatura diferente. Sin embargo, a partir de los
problemas ilustrativos, puede darse cuenta que es difícil que existan problemas que
se puedan tratar sin tener conocimiento del mecanismo mediante el cual el calor se
transfiere entre la superficie de un cuerpo y el medio circundante. Por tanto, en este
capítulo se ampliará el tratamiento de la convección para adquirir una mejor com-
prensión del mecanismo y de algunos de los parámetros clave que influyen en él.
4.2 Transferencia de calor por convección

Antes de intentar calcular un coeficiente de transferencia de calor, se debe examinar
el proceso de convección con cierto detalle y relacionar la transferencia de calor
hacia o desde el fluido en movimiento. En la figura 4.1 se muestra una placa plana
caliente que se enfría por una corriente de aire que circula sobre ella. También se
muestran las distribuciones de velocidad y temperatura que representan esta situa-
ción de convección. El primer punto que se debe observar es que la velocidad dis-
minuye en la dirección hacia la superficie como resultado de fuerzas viscosas que
actúan en el fluido. Como la velocidad de una capa de fluido adyacente a la pared
es cero, la transferencia de calor entre la superficie y esta capa de fluido debe ser,
por conducción:
0T
qc¿¿ = - kf ` = hc1Ts - Tq2
0y y = 0 (4.1)
Si bien esta ecuación sugiere que el proceso se puede considerar como conducción, el
gradiente de temperatura en la superficie (0T/0y)y = 0 se determina por la tasa a la que
el fluido alejado de la pared puede transportar la energía hacia la corriente principal.
Flujo de fluido
∂T
U∞ ∂y y=0
y
u(y) Superficie
T(y)
de la placa
qc T∞ Ts
Flujo de calor
FIGURA 4.1 Distribuciones de la velocidad y la

temperatura en un flujo por convección forzada
laminar sobre una placa plana calentada a tem-
peratura Ts.
231
232 Capítulo 4 Análisis de transferencia de calor por convección
Por tanto, el gradiente de temperatura en la pared depende del campo de flujo, donde
las velocidades mayores pueden producir gradientes de temperatura mayores y tasas
mayores de transferencia de calor; cómo se explica más adelante, la transferencia de
calor por convección en un flujo turbulento a una velocidad mayor por lo general es
mayor que en un flujo laminar a una velocidad menor de los fluidos. Sin embargo,
al mismo tiempo la conductividad térmica del fluido tiene un efecto. Por ejemplo, el
valor de kf para el agua es un orden de magnitud mayor que para el aire; por tanto,
como se muestra en la tabla 1.3, el coeficiente de transferencia de calor por convección
para el agua es mayor que para el aire.
EJEMPLO 4.1 Sobre una placa plana cuya temperatura superficial es 100 °C circula aire a 20 °C.
En una cierta ubicación, la temperatura se mide como una función de la distancia
desde la superficie de la placa; los resultados están trazados en la figura 4.2. Con
estos datos, determine el coeficiente de transferencia de calor por convección en esa
ubicación.
SOLUCIÓN De la ecuación (4.1), el coeficiente de transferencia de calor se puede expresar en la

forma
- kf (0T> 0y)y = 0
hc =
Ts - Tq
La tabla 28 del apéndice 2 muestra que la conductividad térmica del aire a la tempe-
ratura promedio entre la placa y la corriente de fluido (60 °C) es 0.028 W/m K. El
gradiente de temperatura 0T/0y en la superficie se obtiene de manera gráfica trazando
una línea tangente a los datos de temperatura medidos que se muestran en la figura
4.2. De esta manera se obtiene (0T/0y)y = 0 L -66.7 K/mm. Sustituyendo este valor del
gradiente en la superficie calentada de la placa en la ecuación (4.1) se obtiene
-10.028 W/m K21- 66.7 K/mm2

hc = * 103 mm/m
(100 - 20) K
= 23.3 W/m K
2
= Mediciones experimentales
∂T (20 − 100) K
= ΔT = = −66.7 K
2 mm ∂y y=0 Δy 1.2 mm mm
Flujo
de aire
1 mm
Superficie de la placa
20 °C 50 °C 100 °C
FIGURA 4.2 Datos experimentales sobre la distribución

de la temperatura del ejemplo 4.1.
4.3 Fundamentos de la capa límite 233
u(y) y
T∞ T(y)
∂T
∂y y=0
Tg
FIGURA 4.3 Distribuciones de

temperatura y velocidad en con-
vección natural sobre una placa
caliente inclinada en un ángulo b
con respecto a la horizontal.
La situación es muy similar en convección natural, como se muestra en la

figura 4.3. Se debe recordar que en la convección natural o libre, el movimiento
del fluido lo ocasionan los efectos de flotación debidos a la inversión de la den-
sidad del fluido (lo que requiere una diferencia de temperatura), en tanto que en
convección forzada, un gradiente de presión impuesto acciona el flujo de fluido.
La diferencia principal en su distribución de velocidad respectiva es que en con-
vección forzada la velocidad se aproxima al valor de corriente libre impuesto por
una fuerza externa, en tanto que en convección natural, la velocidad primero se
incrementa al aumentar la distancia desde la placa debido a que el efecto de la vis-
cosidad disminuye rápidamente en tanto que la diferencia de densidad disminuye
lentamente. Sin embargo, con el tiempo la fuerza de flotación disminuye conforme
la densidad del fluido se aproxima al valor del fluido circundante, lo que ocasiona
que la velocidad primero alcance un máximo y después tienda a cero alejada de
la superficie calentada. Los campos de temperatura en convección natural y for-
zada tienen formas similares y en los dos casos, el mecanismo de transferencia de
calor en la interfaz fluido/sólido es de conducción.
El análisis anterior indica que el coeficiente de transferencia de calor por con-
vección depende de la densidad, viscosidad y velocidad del fluido así como de sus
propiedades térmicas. En convección forzada, la velocidad se suele imponer en el
sistema por una bomba o un ventilador y se puede especificar de manera directa.
En convección natural, la velocidad depende del efecto de flotación debido a la
diferencia de temperatura entre la superficie y el fluido, del coeficiente de dilata-
ción térmica del fluido, lo que determina el cambio de densidad por diferencia de
temperatura unitaria y del campo de fuerza, que en sistemas ubicados en la Tierra es
simplemente la fuerza gravitacional.
4.3 Fundamentos de la capa límite

Para comprender los parámetros que son significativos en la convección forzada, se
examinará el campo de flujo con más detalle. En la figura 4.4 se muestra la distri-
Región laminar Transición Turbulenta

U∞ d
y Subcapa Superficie
U∞ U∞ viscosa de la placa
U∞
x
du
x=0 dy
FIGURA 4.4 Perfiles de velocidad en capas de límite laminar,

de transición y turbulenta en flujo sobre una placa plana.
bución de velocidad a varias distancias del borde de ataque de una placa. Desde el
borde de ataque hacia adentro, se desarrolla una región en el flujo en la que las fuer-
zas viscosas ocasionan que el fluido aminore su velocidad. Estas fuerzas viscosas
dependen del esfuerzo cortante t. En flujo sobre una placa plana, la velocidad del
fluido paralela a la placa se puede utilizar para definir este esfuerzo como
du
t = m (4.2)
dy
donde du/dy es el gradiente de velocidad y la constante de proporcionalidad m se

denomina viscosidad dinámica. Si el esfuerzo cortante se expresa en newtons por
metro cuadrado y el gradiente de velocidad en (segundos)-1, entonces m tiene las
unidades de newton-segundos por metro cuadrado (N s/m2).*
La región de flujo cerca de la placa donde la velocidad del fluido disminuye por
las fuerzas viscosas se denomina capa límite. La distancia desde la placa a la que la
velocidad alcanza 99% de la velocidad de corriente libre se designa arbitrariamente
como el espesor de la capa límite y la región más allá de este punto se denomina co-
rriente libre sin perturbar o régimen de flujo potencial.
Inicialmente, el flujo en la capa límite es completamente laminar. El espesor
de la capa límite crece al aumentar la distancia desde el borde de ataque y a cierta
distancia crítica, xc, los efectos inerciales se vuelven lo suficientemente grandes
comparados con la acción de amortiguamiento viscoso que las pequeñas perturba-
ciones en el flujo comienzan a crecer. Conforme estas perturbaciones se amplifican,
la regularidad del flujo viscoso se perturba y tiene lugar una transición de flujo
laminar o flujo turbulento. En la región de flujo turbulento, trozos macroscópicos
de fluido se mueven a través de las líneas de corriente y vigorosamente transportan
energía térmica así como cantidad de movimiento. Como se muestra en mecánica
de fluidos [1], el parámetro adimensional que relaciona de manera cuantitativa las
*Observe que m también se puede expresar en las unidades de kg/m s. Con estas unidades es necesario
emplear la constante de conversión gc por consistencia dimensional y la ecuación 4.2 sería entonces
t = (m/gc) du/dy.
4.4 Ecuaciones de conservación de masa, cantidad de movimiento... 235
fuerzas viscosas e inerciales y cuyo valor determina la transición de flujo laminar a

turbulento es el número de Reynolds* Rex, que se define como
rUq x Uq x
Rex = = (4.3)
m n
donde Uq = velocidad de corriente libre

x = distancia desde el borde de ataque
Ȟ = m/r = viscosidad cinemática del fluido
r = densidad del fluido
En la figura 4.4 se muestran bosquejos de formas aproximadas de los perfiles de
velocidad en flujo laminar y turbulento. En el intervalo laminar, el perfil de la capa
límite es aproximadamente parabólico. En el intervalo turbulento, existe una
capa delgada cerca de la superficie (la subcapa viscosa), a través de la que el perfil
de velocidad es casi lineal. Fuera de esta capa el perfil de velocidad es plano com-
parado con el perfil laminar.
El valor crítico de Rex al cual ocurre la transición, Rec,x, depende de la rugosidad
superficial y del nivel de actividad turbulenta, el nivel de turbulencia, en la corriente
principal. Cuando en el flujo principal están presentes perturbaciones grandes, la
transición comienza cuando Rex = 105, pero en campos de flujo menos perturbados,
no empezará hasta Rex = 2 * 105 [1, 2]. Si el flujo está muy libre de perturbaciones,
la transición puede que no ocurra hasta que Re = 106. Por ejemplo, considere el flujo
de aire a 1 m/s y a 20 °C paralelo a una placa plana. Dado que Ȟ = 15.7 * 10-6 m2/s,
la distancia desde el borde de ataque donde ocurre la transición está dada por
m2
Rec,x n (105)a15.7 * 10-6 b
s
xc = = = 1.57 m
a1 b
Uq m
s
El régimen de transición se extiende hasta un número de Reynolds de aproximada-

mente el doble del valor al que comienza la transición y más allá de este punto la
capa límite es turbulenta.
4.4 Ecuaciones de conservación de masa, cantidad de movimiento

y energía para flujo laminar sobre una placa plana
En el enfoque clásico de la convección, se deducen ecuaciones diferenciales para
los equilibrios de cantidad de movimiento y energía en la capa límite y después se
resuelven estas ecuaciones para el gradiente de temperatura en el fluido en la interfaz
*El número de Reynolds, que describe la similitud dinámica adimensional de flujos de fluidos dada por
la relación de las fuerzas inerciales a las viscosas, se nombra así en honor del matemático y profesor
de ingeniería británico, Osborne Reynolds (1842-1912), que nació en Irlanda; estudió en la Cambridge
University; trabajó en el Owens College en Manchester, Inglaterra y realizó experimentos que desarro-
llaron la base de la similitud para determinar la transición de flujos laminares a turbulentos.
fluido/pared para evaluar el coeficiente de transferencia de calor. Un enfoque un

poco más simple, pero más útil es deducir ecuaciones integrales en vez de diferen-
ciales y utilizar un análisis aproximado para obtener una solución. En esta sección,
se deducirán las ecuaciones diferenciales que rigen el flujo de un fluido sobre una
placa plana, para ilustrar la similitud entre transferencia de calor y cantidad de mo-
vimiento y para introducir parámetros adimensionales apropiados que se relacionan
con el proceso. Después se desarrollarán y resolverán las ecuaciones integrales para
flujo sobre una superficie plana, para ilustrar un enfoque analítico que también se
utilizará para obtener los coeficientes de la capa límite de transferencia de calor en
flujo turbulento.
Para deducir la ecuación de conservación de masa o de continuidad, considere
un volumen de control dentro de la capa límite, como se muestra en la figura 4.5
y suponga que prevalecen condiciones en régimen permanente. No hay gradientes
en la dirección z (perpendicular al plano del bosquejo) y el fluido es incompresible.
Entonces las tasas de flujo másico hacia y fuera del volumen de control, respectiva-
mente, en la dirección x son
0u
ru dy y r au + dxbdy
0x
Por tanto, el flujo másico neto hacia el elemento en la dirección x es
0u
-r dx dy
0x
De manera similar, el flujo másico neto hacia el volumen de control en la dirección
y es
0v
-r dx dy
0y
U∞
dx
Superficie de
dy la placa plana
a)
∂υ
(
ρ s+
∂y )
dy
ρ μ + ∂μ dx
ρs dy ( ∂x )
dx
(x, y)
ρs
FIGURA 4.5 Volumen de control
(dx dy ∙ 1) para la conservación de masa
en una capa límite incompresible en
flujo sobre una placa plana.
4.4 Ecuaciones de conservación de masa, cantidad de movimiento... 237
Como el gasto másico neto que sale del volumen de control debe ser cero, se ob-
tiene
0u 0v
- ra + bdx dy = 0
0x 0y
de donde se deduce que en un flujo en régimen permanente bidimensional, la con-
servación de la masa requiere que
0u 0v
+ = 0 (4.4)
0x 0y
La ecuación de conservación de la cantidad de movimiento se obtiene de la
aplicación de la segunda ley del movimiento de Newton al elemento. Suponiendo
que el flujo es newtoniano, que no hay gradientes de presión en la dirección y
y que el cortante viscoso en la dirección y es insignificante, las tasas de flujo de la
cantidad de movimiento en la dirección x para el fluido circulando a través de las caras
verticales izquierda y derecha (consulte la figura 4.6) son ru2 dy y r[u + (0u/0x)
dx]2 dy, respectivamente. Sin embargo, se debe observar que el flujo a través de las
caras horizontales también contribuirá al equilibrio de la cantidad de movimiento en
la dirección x. El flujo de la cantidad de movimiento x entrante a través de la cara
inferior es ruv dx y el flujo de la cantidad de movimiento por ancho unitario saliente
a través de la cara superior es
0v 0u
r av + dyb au + dybdx
0y 0y
Flujos de cantidad de movimiento

Â (s + ∂υ dy) ( u + ∂u dy)
∂y ∂y
(Â u)u dy Â (u + ∂u dx) . ( u + ∂u dx)

∂x ∂x
dx
(x, y)
(Âs)u
Fuerzas
∂u
r y + dy = + ∂ ( ∂uy ) dy
∂y ∂y
∂p
p dy ( p + ∂x dx )
dx
(x, y)
r y = ∂u
∂y
FIGURA 4.6 Volumen de control diferencial para la con-
servación de la cantidad de movimiento en una capa límite
incompresible bidimensional.
La fuerza cortante viscosa en la cara inferior es t ƒ y = -m (0u/0y) dx y sobre la cara

superior es
0u 0 0u
ƒ t ƒ y + dy = m dxc + a bdy d
0y 0y 0y
Por tanto, el cortante viscoso neto en la dirección x es m dx(02u/0y2) dy.
La fuerza de presión sobre la cara izquierda es p dy y sobre la cara derecha es -[p
+ (0p/0x) dx] dy. Por tanto, la fuerza de presión neta en la dirección de movimiento es
-(0p/0x) dx dy. Igualando la suma de las fuerzas con la velocidad de circulación de la
cantidad de movimiento que sale del volumen de control en la dirección x da
0u 2
0v 0u
rau + dxb dy - ru2 dy + r av + dyb au + dybdx - rvu dx
0x 0y 0y
0 2u 0p
= m dx dy - dx dy
0y 2 0x
Ignorando las diferenciales de segundo orden y utilizando la ecuación de conser-
vación de la masa, la ecuación de conservación de la cantidad de movimiento se
reduce a
0u 0u 0 2u 0p
rau + v b = m 2 - (4.5)
0x 0y 0y 0x
En la figura 4.7 se muestra la tasa a la que la energía se conduce por convección
hacia dentro y hacia fuera del volumen de control. Hay cuatro términos convectivos
además de los términos conductivos deducidos en el capítulo 2. Un equilibrio de ener-
gía requiere que la tasa neta de conducción y convección sea cero. Esto produce
0 2T 0 2T 0T 0u 0u 0T
k dx dya + b - crcp au + T + dxbddx dy
0x 2
0y 2 0x 0x 0x 0x
0T 0v 0v 0T
- crcp av + T + dybd dx dy = 0
0y 0y 0y 0y
( ∂T ∂T
2
− kdx + dy ) ρ cp (s + ∂s dy ) (T + ∂T dy ) dx
∂y ∂y2
∂y ∂y
ρ cp (u + ∂u dx ) (T + ∂T dx ) dy
ρ cpuTdy ∂y ∂x
dy
∂T
− kdy
( ∂T + ∂ T dx )
2
dx − kdy
∂x
∂x ∂x2
∂T
− kdx ρ s cpTdx
∂y
FIGURA 4.7 Volumen de control diferencial para la con-

servación de energía.
4.5 Ecuaciones adimensionales de la capa límite y parámetros de similitud 239
La ecuación de la conservación de energía se dedujo con base en la suposición de

que todas las propiedades físicas son dependientes de la temperatura y que la velocidad
del flujo es suficientemente pequeña que se puede ignorar la disipación viscosa. Cuando
la disipación viscosa no se puede ignorar, la energía mecánica se convierte irreversible-
mente en energía térmica, dando origen a un término adicional en el lado izquierdo de
la ecuación. El término, denominado disipación viscosa £, está dado por
0u 0v 2 0u 2 0v 2 2 0u 0v 2
£ = mca + b + 2c a b + a b d - a + b s
0y 0x 0x 0y 3 0x 0y
para propiedades constantes. El efecto de la disipación viscosa puede ser signi-

ficativo si el fluido es muy viscoso, como en las chumaceras, o si la velocidad
de cortante en el fluido es extremadamente alta [3, 4].
Utilizando la ecuación de conservación de la masa e ignorando los términos de
segundo orden, como se hizo en la deducción de la ecuación de conservación de la
cantidad de movimiento, da la expresión siguiente para la ecuación de la energía sin
disipación:
0T 0T 0 2T 0 2T
u + v = aa 2 + 2 b (4.6)
0x 0y 0x 0y
Como una capa límite es muy delgada, en condiciones normales 0T/0y W 0T/0x.
Además, el término de presión en la ecuación de la cantidad de movimiento es cero
para flujo sobre una placa plana ya que (0Uq/0x) = 0. Entonces la similitud entre las
ecuaciones de la cantidad de movimiento y de la energía se vuelve aparente:
0u 0u 0 2u
u + v = va 2 b (4.7a)
0x 0y 0y
0T 0T 0 2T
u + v = aa 2 b (4.7b)
0x 0y 0y
En las relaciones anteriores, Ȟ es la viscosidad cinemática, que es igual a m/r y con
frecuencia se denomina difusividad de la cantidad de movimiento. La relación Ȟ/a
es igual a (m/r)/(k/pcp), que es el número de Prandtl, Pr:
cpm n
Pr = = (4.8)
k a
Si Ȟ es igual a a, entonces Pr es 1 y las ecuaciones de la cantidad de movimiento
y energía son idénticas. Para esta condición, las soluciones no dimensionales de
u(y) y T(y) son idénticas si las condiciones de frontera son similares. Por tanto, es
evidente que el número de Prandtl, que es la relación de las propiedades del fluido
(o más rigurosamente la relación de la cantidad de movimiento y las difusividades
térmicas), controla la relación entre las distribuciones de velocidad y temperatura.
4.5 Ecuaciones adimensionales de la capa límite y parámetros de similitud

Las soluciones de la ecuaciones (4.7a y b), llamadas ecuaciones de capa límite lami-
nar para convección forzada a baja velocidad, producirán los perfiles de velocidad
y temperatura. En general, estas soluciones son muy complicadas y el lector debe

consultar Schlichting [1), Van Driest [2] y Langhaar [5] para ver el planteamiento de
los procedimientos matemáticos, los cuales están fuera del alcance de este libro. Sin
embargo, se puede obtener una visión adicional considerable de los aspectos físicos
del flujo de la capa límite así como de la forma de los parámetros de similitud que
rigen los procesos de transporte al dimensionar las ecuaciones gobernantes, incluso
sin resolverlas.
En la figura 4.8 se muestra el desarrollo de las capas límite de velocidad y tér-
mica para el flujo sobre una superficie plana de forma arbitraria. Para expresar las
ecuaciones de la capa límite en forma adimensional, se definen las variables adimen-
sionales siguientes, que son similares a las definidas en la sección 2.2:
x v
x* = v* =
L Uq
y p
y* = p* =
L rqU2q
u T - Ts
u* = T* =
Uq Tq - Ts
donde L es una dimensión de longitud característica como la longitud de una placa,

Uq es la velocidad de corriente libre, Ts es la temperatura superficial, Tq es la tem-
peratura de corriente libre y rq es la densidad de corriente libre.
Sustituyendo las variables adimensionales anteriores en las ecuaciones dimensio-
nales (4.4), (4.5) y (4.7b) se obtienen las ecuaciones correspondientes de capa límite:
0u* 0v*
+ = 0 (4.9a)
0x* 0y*
0u* 0u* 0p* 1 0 2u*
u* + v* = - + (4.9b)
0x* 0y* 0x* ReL 0y*2
0T* 0T* 1 0 2T*
u* + v* = (4.9c)
0x* 0y* ReLPr 0y*2
Observe que al no dimensionar las ecuaciones de capa límite se han puesto en una
forma en la que aparecen los parámetros de similitud adimensionales ReL y Pr. Estos
parámetros de similitud permiten aplicar soluciones de un sistema a otro sistema
U∞
qʺ T∞ U∞
T∞ δT
δ y
U∞, T∞ u(y)
T(y)
x
δ = Capa límite de velocidad

Ts
δT = Capa límite térmica
FIGURA 4.8 Desarrollo de las capas límite de velocidad y térmica

en flujo sobre una superficie plana de forma arbitraria.
4.5 Ecuaciones adimensionales de la capa límite y parámetros de similitud 241
geométricamente similar, siempre que los parámetros de similitud tengan el mismo

valor en los dos. Por ejemplo, si el número de Reynolds es el mismo, las distribu-
ciones de velocidad adimensionales para aire, agua y glicerina fluyendo sobre una
placa plana serán las mismas a valores dados de x*.
Al analizar la ecuación (4.9a) se observa que v* está relacionada con u*, y* y x*:
v* = f1(u*, y*, x*) (4.10)
y que de la ecuación (4.9b) la solución para u*, como corresponde, se puede expre-
sar en la forma
0p*
u* = f2 ax*, y*, ReL, b (4.11)
0x*
La distribución de la presión sobre la superficie de un cuerpo se determina por
su forma. De la ecuación de la cantidad de movimiento y se puede demostrar que
0p*/0y* = 0 y p* sólo es una función de x*. De aquí, dp*/dx* se puede obtener de
manera independiente y representa la influencia de la forma sobre la distribución
de la velocidad en la corriente libre justo fuera de la capa límite.
4.5.1 Coeficiente de fricción

De la ecuación (4.2), el esfuerzo cortante en la superficie ts está dado por
0u mUq 0u*
ts = m ` = ` (4.12)
0y y = 0 L 0y* y* = 0
Definiendo el coeficiente de resistencia de rozamiento local Cf como
ts
Cfx = (4.13)
rU2q/2
y sustituyendo la ecuación (4.12) para ts da
2 0u*
Cfx = ` (4.14)
ReL 0y* y* = 0
De la ecuación (4.11), es aparente que el gradiente de velocidad adimensional

0u*/0y* en la superficie (y* = 0) depende sólo de x*, ReL y dp*/dx*. Pero dado que
dp*/dx* se determina por completo por la forma geométrica de un cuerpo, la ecua-
ción (4.14) se reduce a la forma
2
Cfx = f (x*, ReL) (4.15)
ReL 3
para cuerpos de forma similar. La relación anterior implica que para flujo sobre
cuerpos de forma similar el coeficiente de rozamiento local está relacionado a x* y
ReL mediante una función universal que es independiente de la velocidad del fluido
o de la velocidad de corriente libre.
_
La resistencia de rozamiento promedio sobre un cuerpo t se puede determinar _
integrando el esfuerzo cortante local t sobre la superficie del__cuerpo. De aquí, t debe
ser independiente de x* y el coeficiente de fricción promedio Cf sólo depende del valor
del número de Reynolds para el flujo sobre cuerpos geométricamente similares:
tq 2
Cqf = = f (Re ) (4.16)
2
rUq /2 ReL 4 L
EJEMPLO 4.2 Para flujo sobre una superficie ligeramente curva, el esfuerzo cortante local está
dado por la relación
rm 0.5 1.5
ts(x) = 0.3a b Uq
x
A partir de esta ecuación dimensional, obtenga relaciones no dimensionales para los

coeficientes de fricción local y promedio.
SOLUCIÓN De la ecuación (4.13), el coeficiente de fricción local es
rm 0.5 U1.5
b
ts(x) q
Cfx = = 0.6a
1 2
2 rUq
x rU2q
0.5
b =
m 0.6 0.6
= 0.6a 0.5
=
rUq x Rex (ReLx*)0.5
Integrando el valor_local y dividiendo entre el área por ancho unitario (L * 1) da el

cortante promedio t:
L
rm 0.5 rm 0.5 1.5
0.3a b U1.5 b Uq
1
tq = q dx = 0.6a
L L0 x L
y por tanto el coeficiente de fricción promedio es
tq 1.2
Cq f = 2
=
rUq/2 Re0.5
L
4.5.2 Número de Nusselt

En la transferencia de calor por convección, la incógnita clave es el coeficiente de
transferencia de calor. De la ecuación (4.1), se obtiene la ecuación siguiente en tér-
minos de los parámetros adimensionales:
kf (Tq - Ts) 0T* kf 0T*

hc = - c d ` = + ` (4.17)
L (Ts - Tq) 0y* y* = 0 L 0y* y* = 0
4.6 Evaluación de los coeficientes de transferencia de calor por convección 243
El análisis de esta ecuación sugiere que la forma adimensional apropiada del coeficiente
de transferencia de calor es el denominado número de Nusselt, Nu, definido por
0T*
`
hcL
Nu = K (4.18)
kf 0y* y * = 0
De las ecuaciones (4.9a) y (4.9c), es aparente que para una geometría preescrita el
número de Nusselt local depende sólo de x*, ReL y Pr:
Nu = f5(x*, ReL, Pr) (4.19)
Una vez que se conoce esta relación funcional ya sea a partir de un análisis o de
experimentos con un fluido particular, se puede utilizar para obtener el valor de Nu
para otros fluidos y para cualesquiera valores de Uq y L. Además, del valor local
de Nu, primero se puede obtener el valor_ local de hc y después un valor promedio___
del coeficiente de transferencia de calor hc y un número de Nusselt promedio NuL.
Como el coeficiente de transferencia de calor promedio se obtiene integrando sobre
la superficie de transferencia de calor de un cuerpo, es independiente de x* y el
número de Nusselt promedio es una función de ReL y Pr:
qhcL
NuL = = f6(ReL, Pr) (4.20)
kf
4.6 Evaluación de los coeficientes de transferencia de calor por convección

Existen cinco métodos generales para la evaluación de los coeficientes de transfe-
rencia de calor por convección:
1. Análisis dimensional combinado con experimentos
2. Soluciones matemáticas exactas de las ecuaciones de la capa límite
3. Análisis aproximados de las ecuaciones de la capa límite mediante métodos
integrales
4. La analogía entre transferencia de calor y de cantidad de movimiento
5. Métodos de análisis numérico, o modelado con dinámica de fluidos compu-
tacional (CFD)
Estas cinco técnicas han contribuido a la comprensión de la transferencia de calor
por convección. No obstante ningún método por sí solo puede resolver todos los proble-
mas, debido a que cada uno tiene limitaciones que restringen su alcance de aplicación.
El análisis dimensional es matemáticamente simple y ha encontrado un inter-
valo amplio de aplicación [5, 6]. La limitación principal de este método es que los
resultados obtenidos son incompletos y muy inútiles sin datos experimentales. El
análisis dimensional contribuye poco al entendimiento del proceso de transferencia,
pero facilita la interpretación y extiende el intervalo de datos experimentales al
correlacionarlo en términos de grupos adimensionales.
Existen dos métodos diferentes para determinar grupos adimensionales adecuados
para correlacionar datos experimentales. El primero de estos métodos, analizado en la
sección siguiente, sólo requiere una lista de las variables pertinentes a un fenómeno.
Esta técnica es simple de utilizar, pero si se omite una variable pertinente, se originan
resultados erróneos. En el segundo método, los grupos adimensionales y las condi-
ciones de similitud se deducen a partir de ecuaciones deferenciales que describen
el fenómeno. Este método se prefiere cuando el fenómeno se puede describir mate-

máticamente, pero la solución de las ecuaciones resultantes con frecuencia es dema-
siado complicada para ser práctica. Esta técnica se presentó en la sección 4.5.
Los análisis matemáticos exactos requieren la solución simultánea de las
ecuaciones que describen el movimiento de fluido y la transferencia de energía en
el fluido en movimiento [7]. El método presupone que los mecanismos físicos se
comprenden lo suficiente para describirlos en lenguaje matemático. Este reque-
rimiento preliminar limita el alcance de las soluciones exactas debido a que las
ecuaciones matemáticas completas que describen los mecanismos de flujo de fluido
y de transferencia de calor sólo se pueden escribir para flujo laminar. Incluso para
flujo laminar, las ecuaciones son muy complicadas, pero se han obtenido soluciones
para una variedad de sistemas simples como el flujo sobre una placa plana, sobre
una superficie aerodinámica o sobre un cilindro circular [7].
Las soluciones exactas son importantes dado que las suposiciones hechas en el
curso del análisis se pueden especificar con precisión y su validez se puede verificar
con experimentos. También sirven como una base de comparación y comprobación
en métodos aproximados más simples. Además, el desarrollo de computadoras de
alta velocidad ha aumentado el intervalo de problemas apropiados para su solu-
ción matemática y los resultados de cálculos para sistemas diferentes se publican
de manera continua en la bibliografía especializada.
El análisis aproximado de la capa límite evita la descripción matemática deta-
llada del flujo en la capa límite. En su lugar, se utiliza una ecuación simple, pero razo-
nable, para describir las distribuciones de velocidad y temperatura en la capa límite.
Después el problema se analiza en una base macroscópica aplicando la ecuación de
movimiento y la ecuación de energía al total de las partículas de fluido contenidas
dentro de la capa límite. Este método es relativamente simple; además, produce solu-
ciones para problemas que no se pueden tratar con un análisis matemático exacto.
En casos donde existan otras soluciones, éstas concuerdan dentro de la precisión
ingenieril con las soluciones obtenidas mediante este método aproximado. La técnica
no está limitada a flujo laminar sino que también se puede aplicar a flujo turbulento.
La analogía entre la transferencia de calor y de cantidad de movimiento es una
herramienta útil para analizar procesos de transferencia turbulenta. Nuestro cono-
cimiento de los mecanismos de intercambio turbulento es insuficiente para escribir
ecuaciones matemáticas que describan la distribución de la temperatura de manera
directa, pero el mecanismo de transferencia de calor se puede describir en términos
de un modelo simplificado. De acuerdo con un modelo ampliamente aceptado, un
movimiento mezclado en una dirección perpendicular al flujo principal explica la
transferencia de cantidad de movimiento así como de energía. El movimiento de mez-
clado se puede describir con una base estadística mediante un método similar al que
se utilizó para representar el movimiento de moléculas de gas en la teoría cinética. De
ninguna manera puede decirse que hay un acuerdo de que este modelo corresponde
a condiciones que en realidad existen en la naturaleza, pero para fines prácticos, su
uso se puede justificar por el hecho de que los resultados experimentales concuerdan
sustancialmente con predicciones analíticas basadas en el modelo hipotético.
Con los métodos numéricos es posible resolver en una forma apropiada las
ecuaciones exactas de movimiento [8, 9]. La aproximación resulta de la necesidad
de expresar las variables de campo (temperatura, velocidad y presión) en puntos dis-
cretos en el tiempo y el espacio en lugar de continuamente. Sin embargo, la solución
se puede hacer lo suficientemente precisa si se tiene el cuidado de discretizar las
4.7 Análisis dimensional 245
ecuaciones exactas. Una de las ventajas más importantes de los métodos numéricos es
que una vez que se haya programado el procedimiento de solución, es posible calcular
con facilidad soluciones para diferentes condiciones de frontera, variables de propie-
dades, etc. En general, con los métodos numéricos se pueden manejar con facilidad
condiciones de frontera complejas. Los métodos numéricos para resolver problemas
de convección se analizan en [9] y son una extensión de los métodos presentados en
el capítulo 3 para problemas de conducción.
4.7 Análisis dimensional

El análisis dimensional difiere de otros enfoques en que no produce ecuaciones que se
puedan resolver. En lugar de eso, combina diferentes variables en grupos adimensiona-
les, como el número de Nusselt, que facilitan la interpretación y extienden el intervalo
de aplicación de datos experimentales. En la práctica, los coeficientes de transferencia de
calor por convección por lo general se calculan a partir de ecuaciones empíricas obteni-
das correlacionando datos experimentales con ayuda del análisis dimensional.
La limitación más seria del análisis dimensional es que no proporciona información
acerca de la naturaleza del problema. De hecho, para aplicar el análisis dimensional es
necesario saber de antemano qué variables influyen en el fenómeno y el éxito o falla del
método depende de la selección apropiada de estas variables. Por tanto, es importante
tener al menos una teoría preliminar o una comprensión física completa de un fenómeno
antes de que se pueda efectuar un análisis dimensional. Sin embargo, una vez que se
conocen las variables pertinentes, el análisis dimensional se puede aplicar a la mayoría
de los problemas mediante un procedimiento de rutina descrito a continuación.*
4.7.1 Dimensiones primarias y fórmulas dimensionales

El primer paso es seleccionar un sistema de dimensiones primarias. La elección de
las dimensiones primarias es arbitraria, pero las fórmulas dimensionales de todas las
variables pertinentes se deben poder expresar en términos de ellas. En el sistema SI, se
utilizan la dimensiones primarias de longitud L, tiempo t, temperatura T y masa M.
La fórmula dimensional de una cantidad física proviene de definiciones o
leyes físicas. Por ejemplo, la fórmula dimensional para la longitud de una barra
es [L] por definición.† La velocidad promedio de una partícula de fluido es igual
a la distancia dividida entre el intervalo de tiempo que le toma recorrerla. Por tanto,
la fórmula dimensional de la velocidad es [L/t] o [Lt-1] (es decir, una distancia
o longitud dividida entre el tiempo). Las unidades de la velocidad se podrían expre-
sar en metros por segundo, pies por segundo o millas por hora ya que todas son una
longitud dividida entre un tiempo. Las fórmulas dimensionales y los símbolos de
cantidades físicas que se presentan con frecuencia en problemas de transferencia
de calor se dan en la tabla 4.1.
*La teoría algebraica del análisis dimensional no se desarrollará aquí. Para consultar un tratamiento
riguroso y amplio de los fundamentos matemáticos se recomiendan los capítulos 3 y 4 de Langhaar [5].
†Los paréntesis rectangulares indican que la cantidad tiene la fórmula dimensional indicada entre los
paréntesis.
TABLA 4.1 Cantidades físicas importantes de transferencia de calor, masa

y sus dimensiones
Cantidad Símbolo Dimensiones en el sistema MLtT
Longitud L, x L
Tiempo t t
Masa M M
Fuerza F ML/t2
Temperatura T T
Calor Q ML2/t2
Velocidad u, v, Uq L/t
Aceleración a, g L/t2
Trabajo W ML2/t2
Presión p M/t2L
Densidad r M/L3
Energía interna e L2/t2
Entalpía i L2/t2
Calor específico c L2/t2T
Viscosidad absoluta m M/Lt
Viscosidad cinemática v = m/r L2/t
Conductividad térmica k ML/t3T
Difusividad térmica a L2/t
Resistencia térmica R Tt3/ML2
Coeficiente de dilatación b 1/T
Tensión superficial s M/t2
Esfuerzo cortante t M/Lt2
Coeficiente de transferencia de calor h M/t3T
Gasto másico ∙
m M/t
4.7.2 Teorema p de Buckingham

Para determinar el número de grupos adimensionales independientes necesario para
obtener una relación que describa un fenómeno físico, se puede utilizar el teorema p
de Buckingham.‡ De acuerdo con esta regla, el número requerido de grupos adi-
mensionales independientes que se puede formar combinando las variables físicas
pertinentes a un problema es igual al número de estas cantidades físicas n (por ejem-
plo, densidad, viscosidad, coeficiente de transferencia de calor) menos el número de
dimensiones primarias m necesario para expresar las fórmulas dimensionales de las
n cantidades físicas. Si estos grupos se designan por p1, p2 y así sucesivamente, la
ecuación que expresa la relación entre las variables tiene una solución con la forma
F(p1, p2, p3, . . .) = 0 (4.21)
‡Una regla más rigurosa propuesta por Van Driest [26] demuestra que el teorema p es válido siempre que
el conjunto de ecuaciones simultáneas formado al igualar los exponentes de cada dimensión primaria a cero
sea linealmente independiente. Si una ecuación en el conjunto es una combinación lineal de una o más de
las otras ecuaciones (es decir, si las ecuaciones son linealmente dependientes), el número de grupos adi-
mensionales es igual al número total de variables n menos el número de ecuaciones independientes.
En un problema que comprenda cinco cantidades físicas y tres dimensiones prima-

rias, n – m es igual a dos y la solución tiene la forma
F(p1, p2) = 0 (4.22)
o bien
p1 = f (p2) (4.23)
En ese caso, los datos experimentales se pueden presentar de manera conveniente
trazando p1 contra p2. La curva empírica resultante revela la relación funcional
entre p1 y p2, que no se puede deducir a partir de un análisis dimensional.
Para un fenómeno que se pueda describir en términos de grupos adimensionales
(es decir, si n – m = 3), la ecuación (4.21) tiene la forma
F(p1, p2, p3) = 0 (4.24)
pero también se puede escribir como
p1 = f (p2, p3)
Para ese caso, los datos experimentales se pueden correlacionar trazando p1 contra
p2 para varios valores de p3. En ocasiones es posible combinar de alguna manera
dos de las p y trazar este parámetro contra la p restante en una sola curva.
4.7.3 Determinación de grupos adimensionales

Ahora se ilustrará un método simple para determinar grupos adimensionales apli-
cándolo al problema de correlacionar datos experimentales de transferencia de calor
por convección para un fluido que circula a través de un tubo caliente. Se podría
emplear exactamente el mismo enfoque para transferencia de calor en flujo a través
de un tubo o sobre una placa.
A partir de la descripción del proceso de transferencia de calor por convección,
es razonable esperar que las cantidades físicas dadas en la tabla 4.2 sean pertinentes
al problema.
Existen siete cantidades físicas y cuatro dimensiones primarias. Por tanto, se
espera que se requieran tres grupos adimensionales para correlacionar los datos.
Para determinar estos grupos adimensionales, se escribe p como un producto de las
variables, cada una elevada a una potencia desconocida:
c d e f qg
p = D ak bUq r m cp hc (4.25)
TABLA 4.2 Cantidades físicas pertinentes en transferencia de calor

por convección
Variable Símbolo Dimensiones
Diámetro del tubo D [L]
Conductividad térmica del fluido k [ML/t3T]
Velocidad de corriente libre del fluido Uq [L/t]
Densidad del fluido r [M/L3]
Viscosidad del fluido m [M/Lt]
Calor específico a presión constante cp
_ [L2/t2T]
Coeficiente de transferencia de calor hc [M/t3T]
y se sustituyen las fórmulas dimensionales
ML b L c M d M e L2 f M g
p = [L]a c d c d c 3 d c d c 2 d c 3 d (4.26)
t3T t L Lt tT tT
Para que p sea adimensional, la suma de los exponentes por separado de cada dimen-
sión primaria debe sumar cero. Igualando la suma de los exponentes de cada dimensión
primaria a cero, se obtiene el conjunto de ecuaciones siguiente:
b + d + e + g = 0 para M
a + b + c - 3d - e + 2f = 0 para L
- 3b - c - e - 2f - 3g = 0 para t
- b - f - g = 0 para T
Es evidente que cualquier conjunto de valores de a, b, c, d y e que satisfaga simul-

táneamente este conjunto de ecuaciones hará que p sea adimensional. Hay siete
incógnitas, pero sólo cuatro ecuaciones. Por tanto, se pueden elegir valores para tres
de los exponentes en cada uno de los grupos adimensionales. La única restricción
en la elección de los exponentes es que cada uno de los exponentes seleccionados
sea independiente de los otros. Un exponente es independiente si el determinante
formado con los coeficientes de los términos restantes no desaparece (es decir, no
es igual a cero)._
Dado que hc, el coeficiente de transferencia de calor por convección, es la
variable que finalmente se quiere evaluar, es conveniente igualar su exponente g
a la unidad. El mismo tiempo, se hace c = d = 0 para simplificar las manipulacio-
nes algebraicas. Al resolver las ecuaciones simultáneamente, se obtiene a = 1, b = -1,
e = f = 0. Entonces el primer grupo adimensional es
qhcD
p1 =
k
___
el cual se conoce como el número de Nusselt, NuD_.
Para p2, se selecciona g igual a cero, tal que hc no aparezca de nuevo y se hace
a = 1 y f = 0. La solución simultánea con estas elecciones da b = 0, c = d = 1, e = -1 y
UqDr
p2 =
m
Este grupo adimensional es el número de Reynolds, ReD, con el diámetro del tubo
como el parámetro de longitud.
Si se hace e = 1 y c = g = 0, se obtiene el tercer grupo adimensional,
cpm
p3 =
k
que es el número de Prandtl, Pr.
Se observa que aunque el coeficiente de transferencia de calor por convección
es una función de seis variables, con ayuda del análisis dimensional las siete varia-
bles originales se han combinado en tres grupos adimensionales. De acuerdo con la

ecuación (4.24), la relación funcional se puede escribir como
NuD = f (ReD, Pr)
y ahora los datos experimentales se pueden correlacionar en términos de tres varia-

bles en vez de las siete originales. La importancia de esta reducción en el número de
variables se vuelve aparente cuando se intenta planear experimentos y correlacionar
datos experimentales.
4.7.4 Correlación de datos experimentales

Suponga que en una serie de pruebas con aire fluyendo sobre un tubo de 25 mm de
diámetro exterior, el coeficiente de transferencia de calor se ha medido experimen-
talmente a velocidades variando de 0.15 a 30 m/s. Este intervalo de velocidades
corresponde a números de Reynolds basados en el diámetro DrUq/m que van de 250
a 50 000. Como la velocidad fue la única variable en estas pruebas, los resultados
Coeficiente de transferencia de calor, hc
210
180
150
(W / m2 K)
120
90
60
30
0
0.1 0.2 0.5 1 2 5 10 20 50
U∞, Velocidad (m/s)
a)
160
Número de Nusselt, hc D/ k
140
120
100
80
60
40
20
0
100 200 500 1000 2000 5000 10 000 20 000 50 000 100 000
Número de Reynolds, U∞ DÂ /
b)
FIGURA 4.9 Variación del número de Nusselt con el número de Reynolds para
flujo transversal de aire sobre un tubo o cilindro largo a) gráfica dimensional,
b) gráfica adimensional.
están correlacionados
_ en la figura 4.9a) trazando el coeficiente de transferencia de
calor hc contra
_ la velocidad Uq. La curva resultante permite la determinación
directa de hc en cualquier velocidad para el sistema utilizado en las pruebas, pero
no se puede emplear para determinar los coeficientes de transferencia de calor para
cilindros que sean mayores o menores que el utilizado en las pruebas. Tampoco
se podría evaluar el coeficiente de transferencia de calor si el aire estuviera a presión
y su densidad fuera diferente de la empleada en las pruebas. A menos que los datos
experimentales se pudieran correlacionar de manera más efectiva, sería necesario
efectuar experimentos separados para cada diámetro del cilindro, cada densidad, etc.
El esfuerzo sería enorme.
Sin embargo, con ayuda del análisis dimensional, los resultados de una serie
de pruebas se pueden aplicar a una variedad de otros problemas, como se ilustra en
la figura 4.9b), donde los datos de la figura 4.9a) se volvieron a trazar en términos
de grupos adimensionales pertinentes. La abscisa en la figura_ 4.9b) es el número de
Reynolds UqDr/m y la ordenada es el número de Nusselt hcD/k. Esta correlación
de los datos permite la evaluación del coeficiente de transferencia de calor para aire
fluyendo sobre cualquier tipo de tubo o alambre siempre que el número de Reynolds
del sistema se encuentre dentro del intervalo cubierto en el experimento y que los
sistemas sean geométricamente similares.
Datos experimentales obtenidos solamente con aire no revelan la dependencia
del número de Nusselt en el número de Prandtl ya que el número de Prandtl es una
combinación de propiedades físicas cuyo valor no varía de manera apreciable para
gases. Para determinar la influencia del número de Prandtl, es necesario utilizar
fluidos diferentes. De acuerdo con el análisis anterior, se necesitan datos experimen-
tales con varios fluidos cuyas propiedades físicas produzcan un intervalo amplio de
números de Prandtl para completar la correlación.
En la figura 4.10 están trazados los resultados experimentales de varias investi-
gaciones independientes para transferencia de calor entre aire, agua y aceites en flujo
transversal sobre un tubo o un alambre para un intervalo amplio de temperaturas,
tamaños de cilindros y velocidades. La ordenada en la figura 4.10 es la cantidad
10.0
La ecuación de correlación para
3 < ReD < 100 es
5.0 0.4
NuD / Pr 0.3 = 0.82 ReD
3.0
±15%
NuD / Pr0.3
2.0
1.0 Símbolo Fluido Diám. (mm) Ref

Agua 0.0254 19
Aceites 0.1-0.2 20
±15% Aire 0.2-150 21
Aire 0.023-0.15 22
0.3
0.1 1.0 10.0 100
ReD = DrU∞ /m
FIGURA 4.10 Correlación de datos experimentales de transferencia
de calor para varios fluidos en flujo transversal sobre tubos, alambres
y cilindros circulares.
___
adimensional* NuD/Pr0.3 y la abscisa es ReD. Un examen de los resultados muestra
que todos los datos siguen una sola línea razonablemente bien y por tanto se pueden
correlacionar empíricamente.
Se debe observar que los datos experimentales de las figuras 4.9 y 4.10 cubren
intervalos diferentes de números de Reynolds entre 0.1 y 100 en la figura 4.10 y entre
aproximadamente 200 y 50 000 en la figura 4.9. La extrapolación de la ecuación de corre-
lación en la figura 4.10 en un intervalo del número de Reynolds mucho mayor que 200
o 300 conduciría a errores graves. Por ejemplo, para aire (Pr = 0.71) fluyendo sobre un
cilindro en ReD = 20 000, la ecuación de correlación en la figura 4.10 anticiparía
NuD = Pr0.3 * 0.82 * Re0.4

D = 0.71
0.3
* 0.82 * 20000 0.4 = 39
en tanto que los datos experimentales en la figura 4.9 dan NuD = 85, una diferencia
sustancial. Cuando no existan datos en un intervalo encontrado en un diseño, puede ser
necesario hacer una extrapolación, pero como se muestra en el ejemplo anterior, la extra-
polación de datos más allá del intervalo de parámetros adimensionales cubiertos en
experimentos se deberá evitar si es posible. Si no hay oportunidad de conducir experi-
mentos apropiados, los resultados de una extrapolación se deben tratar con precaución.
4.7.5 Principio de similitud

El resultado extraordinario de la figura 4.10 se puede explicar mediante el principio
de similitud. De acuerdo con este principio, con frecuencia denominado ley, el com-
portamiento de dos sistemas será similar si las relaciones de sus dimensiones lineales,
fuerzas, velocidades y así sucesivamente son las mismas. En condiciones de convec-
ción forzada en sistemas geométricamente similares, los campos de velocidad serán
similares siempre que la relación de las fuerzas inerciales con relación a las fuerzas
viscosas sea la misma en los dos fluidos. El número de Reynolds es la relación de estas
fuerzas y en consecuencia, se esperan condiciones de flujo similares en convección
forzada para un valor dado del número de Reynolds. El número de Prandt es la rela-
ción de dos propiedades de transporte molecular, la viscosidad cinemática Ȟ = m/r, que
afecta la distribución de velocidad y la difusividad térmica k/rcp, que afecta el perfil de
temperatura. En otras palabras, es un grupo adimensional que relaciona la distribución
de temperatura con la distribución de velocidad. De aquí que en sistemas geométrica-
mente similares con los mismos números de Prandtl y Reynolds, las distribuciones de
temperatura serán similares. El número de Nusselt es igual a la relación del gradiente
de temperatura en una interfaz fluido a superficie a un gradiente de temperatura de refe-
rencia. Por tanto, se espera que, en sistemas con geometrías y campos de temperatura
similares, los valores numéricos de los números de Nusselt serán iguales. Esta predicción
se confirma por los resultados experimentales que aparecen en la figura 4.10.
Se han realizado análisis dimensionales para numerosos sistemas de transferen-
cia de calor y en la tabla 4.3 se resumen los grupos adimensionales más importantes
utilizados en el diseño.
*El combinar el número de Nusselt con el de Prandtl para trazar los datos es simplemente un asunto de
conveniencia. Como ya se mencionó, cualquier combinación de parámetros adimensionales es satisfac-
toria. La selección del parámetro más conveniente suele hacerse con base en la experiencia mediante
prueba y error con ayuda de resultados experimentales, aunque en ocasiones los grupos característicos
se sugieren por los resultados de las soluciones analíticas.
TABLA 4.3 Grupos adimensionales de importancia en la transferencia de calor

y flujo de fluidos
Grupo Definición Interpretación
Número de Biot (Bi) qhL Relación de la resistencia térmica interna
ks de un cuerpo sólido a su resistencia
ts térmica superficial
Coeficiente de arrastre (Cf) Relación del esfuerzo cortante superficial
rU2q/2
a la energía cinética de corriente libre
Número de Eckert (Ec) U2q Energía cinética del flujo relativa a la
cp(Ts - Tq) diferencia de entalpía de la capa límite
Número de Fourier (Fo) at Tiempo adimensional; relación de la tasa
L2 de conducción de calor a la tasa de alma-
cenamiento de energía interna en un sólido
¢p
Factor de fricción ( f ) Caída de presión adimensional para flujo
(L/D)(rU2m/2) interno a través de conductos
Número de Grashof (GrL) gb(Ts - Tq)L3 Relación entre fuerzas de flotación y viscosas
n2
Factor j de Colburn j ( jH) StPr 2/3 Coeficiente de transferencia de calor

adimensional
Número de Nusselt (NuL) hqcL Coeficiente de transferencia de calor
kf adimensional; relación de transferencia
de calor por convección a conducción
en una capa de fluido de espesor L
Número de Peclet (PeL) Re LPr Producto de los números de Reynolds
y Prandtl
cpm n
Número de Prandtl (Pr) = Relación de la difusividad de la cantidad de
k a movimiento molecular a la difusividad
térmica
Número de Rayleigh (Ra) Gr LPr Producto de los números de Grashof y Prandtl
Uq L
Número de Reynolds (ReL) Relación de fuerzas de inercia y viscosas
n
Número de Stanton (St) hqc NuL Coeficiente de transferencia de calor
=
rUqcp ReLPr adimensional
4.8* Solución analítica para el flujo laminar de capa límite sobre

una placa plana†
En la sección anterior se determinaron grupos adimensionales para correlacionar
datos experimentales para transferencia de calor por convección forzada. Se encontró
que el número de Nusselt depende del número de Reynolds y de Prandtl, es decir,
Nu = f(Re)c(Pr) (4.27)
†En el resto de este capítulo, se pueden omitir los detalles matemáticos en un curso introductorio sin
romper la continuidad de la presentación.
4.8 Solución analítica para el flujo laminar de capa límite sobre una placa plana 253
En las siguientes secciones se considerarán métodos analíticos para determinar las

relaciones funcionales en la ecuación (4.27) para flujo a baja velocidad sobre una
placa plana. Este sistema se ha seleccionado principalmente debido a que es el más
simple de analizar, pero los resultados tienen muchas aplicaciones prácticas, por
ejemplo, son buenas aproximaciones para flujo sobre las superficies de cuerpos
aerodinámicos como las alas de aviones o los álabes de turbinas.
En vista de las diferencias en las características de flujo, las fuerzas de fricción
así como la transferencia de calor están gobernadas por relaciones diferentes para
flujo laminar y turbulento. Primero se abordará la capa límite laminar, la cual se
presta para un método de solución exacto y aproximado. La capa límite turbulenta
se analiza en la sección 4.10.
Para determinar el coeficiente de transferencia de calor por convección forzada
y el coeficiente de fricción para flujo incompresible sobre una superficie plana, se
deben satisfacer las ecuaciones de continuidad, cantidad de movimiento y energía
simultáneamente. Estas relaciones se dedujeron en la sección 4.4 y por convenien-
cia se repiten a continuación.
Continuidad:
0u 0v
+ = 0 (4.4)
0x 0y
Cantidad de movimiento:
0u 0u 0 2u 0p
rau + v b = m 2 - (4.5)
0x 0y 0y 0x
Energía:
0T 0T 0 2T
u + v = a 2 (4.7b)
0x 0y 0y
4.8.1 Espesor de la capa límite y fricción en la superficie

La ecuación (4.5) se debe resolver simultáneamente con la ecuación de continuidad,
ecuación (4.4), a fin de determinar la distribución de velocidad, el espesor de la capa
límite y la fuerza de fricción en la pared. Estas ecuaciones se resuelven primero defi-
niendo una función de corriente c(x, y) que automáticamente satisface la ecuación
de continuidad, es decir,
0c 0c
u = y v = -
0y 0x
Introduciendo la nueva variable
Uq
h = y
C nx
se puede hacer
c = 1nxUq f (h)
donde f(h) denota una función de corriente adimensional. En términos de f(h), las
componentes de la velocidad son
0c 0c 0h d[ f (h)]
u = = = Uq
0y 0h 0y dh
y
0c
e h - f (h) f
1 vUq d[ f (h)]
v = - =
0x 2C x dh
Expresando 0u/0x, 0u/0y y 02u/0y2 en términos de h e insertando las expresiones resul-
tantes en la ecuación de la cantidad de movimiento se obtiene la ecuación diferencial
ordinaria, no lineal y de tercer orden
d2[ f (h)] d3[ f (h)]
f (h) + 2 = 0
dh2 dh3
la cual se debe resolver sujeta a las tres condiciones de frontera
d[ f (h)]
f(h) = 0 y = 0 en h = 0
dh
y
d[ f (h)]
= 1 en h = q
dh
La solución para esta ecuación diferencial la obtuvo numéricamente Blasius en 1908

[10]. Los resultados significativos se muestran en las figuras 4.11 y 4.12.
1.0
Pendiente = 0.332
x = 1.0 cm
x = 2.0 cm
x = 2.5 cm
u
x = 4.0 cm
U∞
x = 5.0 cm
0.5 x = 7.5 cm
x = 10.0 cm
U∞ = 8 m/s
x = 12.5 cm
Placa delgada (Núm. 1)
x = 15.0 cm
x = 17.5 cm
Blasius
0
0 1 2 3 4 5 6 7
y rU∞ x
x m
FIGURA 4.11 Perfil de velocidad en una capa límite laminar de acuerdo
con Blasius, con datos experimentales de Hansen [11].
Fuente: Cortesía del National Advisory Committee for Aeronautics, NACA TM 585.
100.0
10.0
Cf x × 103
Intervalo
de transición
1.0
0.1
10 102 103 104 105 106
Rex
FIGURA 4.12 Coeficiente de fricción local contra el número de

Reynolds basado en la distancia desde el borde de ataque para flujo
laminar sobre una placa plana.
En la figura 4.11, los perfiles de velocidad de Blasius en la capa límite laminar

en una placa plana están trazados en forma adimensional junto con datos experimen-
tales obtenidos por Hansen [11]. La ordenada es la velocidad local en la dirección
x, que es u, dividida entre la velocidad de corriente libre Uq y la abscisa es un pará-
________
metro de distancia adimensional (y/x)¥ (rUqx)/m . Se observa que una sola curva es
suficiente para correlacionar las distribuciones de velocidad en todas las estacio-
de la placa. La velocidad u alcanza 99% del valor de corriente libre
nes a lo largo________
Uq en (y/x)¥ (rUqx)/m = 5.0. Si se define el espesor de la capa límite hidrodinámica
como la distancia desde la superficie a la que la velocidad local u alcanza 99% del
valor de corriente libre Uq, el espesor de la capa límite @ se vuelve
5x
d = (4.28)
1Rex
donde Rex = rUqx/m, que es el número de Reynolds local. La ecuación (4.28) satis-
face la descripción cualitativa del crecimiento de la capa límite, @ es cero en el borde
de ataque (x = 0) y aumenta con x a lo largo de la placa. En cualquier situación,
es decir, a cualquier valor dado de x, el espesor de la capa límite es inversamente
proporcional a la raíz cuadrada del número de Reynolds local. De aquí, un aumento
en la velocidad resulta en una disminución en el espesor de la capa límite.
La fuerza cortante en la pared se puede determinar a partir del gradiente de
velocidad en y = 0 en la figura 4.11. Se observa que
0(u> Uq)
` = 0.332
0(y> x) 1Rex y = 0
y por tanto en cualquier valor especificado de x el gradiente de velocidad en la
superficie es
0u
`
Uq
= 0.332 1Rex
0y y = 0 x
Sustituyendo este gradiente de velocidad en la ecuación para el cortante, el cortante

en la pared por área unitaria ts es
0u
`
Uq
ts = m = 0.332m 1Rex (4.29)
0y y = 0 x
Se observa que el cortante en la pared cerca del borde de ataque es muy grande y
disminuye con el aumento en la distancia desde el borde de ataque.
Para una representación gráfica, es más conveniente utilizar coordenadas adi-
mensionales. Dividiendo los dos lados de la ecuación (4.29) entre la presión de la
velocidad de la corriente libre rU q
2 /2, se obtiene
ts 0.664
Cfx = =
rU2q> 2
(4.30)
1Rex
donde Cfx es el coeficiente de rozamiento local adimensional o coeficiente de

fricción. La figura 4.12 es una gráfica de Cfx contra Rex y muestra la variación del
coeficiente de fricción local de manera gráfica. El coeficiente de fricción promedio
se obtiene integrando la ecuación (4.30) entre el borde de ataque x = 0 y x = L:
L m
1
Cqf = Cfx dx = 1.33 (4.31)
L L0 C Uq rL
Por __
tanto, para flujo laminar sobre una placa plana, el coeficiente de fricción prome-
dio Cf es igual al doble del valor del coeficiente de fricción local en x = L:
4.8.2 Transferencia de calor por convección

La ecuación de conservación de la energía para una capa límite laminar es
0T 0T 0 2T
u + v = a 2 (4.7b)
0x 0y 0y
Las velocidades en la ecuación de conservación de la energía, u y Ȟ, tienen los mis-

mos valores en cualquier punto (x, y) como en la ecuación de dinámica de fluidos,
ecuación (4.5). En el caso de la placa plana, Pohlhausen [12] utilizó las velocidades
calculadas antes por Blasius [10] para obtener la solución del problema de transfe-
rencia de calor. Sin considerar los detalles de esta solución matemática, se pueden
obtener resultados significativos comparando la ecuación (4.7b) con la ecuación
(4.5), la ecuación de la cantidad de movimiento. Las dos ecuaciones son similares;
de hecho u(x, y) también es una solución para la distribución de temperatura T(x, y)
si Ȟ = a y si la temperatura de la placa Ts es constante. Esto se puede verificar con
facilidad remplazando el símbolo T en la ecuación (4.7b) por el símbolo u y obser-
vando que las condiciones límites para T y u son idénticas. Si se utiliza la temperatura
superficial como referencia y se deja que la variable en la ecuación (4.7b) sea

(T – Ts)/(Tq - Ts), entonces las condiciones de frontera son
T - Ts u
= 0 y = 0 en y = 0
Tq - Ts Uq
T - Ts u
= 1 y = 1 en y : q
Tq - Ts Uq
donde Tq es la temperatura de corriente libre.

La condición de que Ȟ = a corresponde a un número de Prandtl igual a la unidad
ya que
cpm n
Pr = =
k a
Para Pr = 1, la distribución de velocidad es, por tanto, idéntica a la distribución de

temperatura. Una interpretación en términos de procesos físicos es que la transferen-
cia de cantidad de movimiento es análoga a la transferencia de calor cuando Pr = 1.
Las propiedades físicas de la mayoría de los gases son tales que tienen números
de Prandtl que varían de 0.6 a 1.0 y por tanto, la analogía es satisfactoria. Por otro
lado, los líquidos tiene números de Prandtl considerablemente diferentes y el análisis
anterior no se puede aplicar directamente [13].
Utilizando los resultados analíticos del trabajo de Pohlhausen, la distribución de
temperatura en la capa límite laminar para Pr = 1 se puede modificar empíricamente
para incluir fluidos con números de Prandtl diferentes de la unidad. En la figura
4.13 se muestran perfiles de temperatura calculados teóricamente en la capa límite
para valores de Pr de 0.6. 0.8, 1.0, 3.0, 7.0, 15 y 50. Ahora se define un espesor de
capa límite térmica @t como la distancia desde la superficie a la que la diferencia
de temperatura entre la pared y el fluido alcanza 99% del valor de corriente libre.
Al analizar los perfiles de temperatura se tiene que la capa límite térmica es mayor
que la capa límite hidrodinámica para fluidos con Pr menores que la unidad, pero
1.0
0.8
Pr = 50 15 7 3 1 0.8 0.6
0.6
T∞ − Ts
T − Ts
0.4
0.2
0
0 0.4 0.8 1.2 1.6 2.0 2.4 2.8 3.2 3.6 4.0
y U∞ x
h=x Ȟ
FIGURA 4.13 Distribuciones de temperatura adimensionales
en un fluido que circula sobre una placa caliente para varios
números de Prandtl.
1.5
∂ (T − TS)/(T∞ − TS)
y = 0.332
∂ x √Rex Pr 1/3 y=0
1.0
T∞ − TS
T − TS
Distribución de temperatura
0.5
0
0 1.0 2.0 3.0 4.0 5.0
y
x √RexPr
1/3
FIGURA 4.14 Distribución de temperatura adimensional

para flujo laminar sobre una placa caliente a temperatura
uniforme.
son menores cuando Pr es mayor que la unidad. De acuerdo con los cálculos de
Pohlhausen, la relación entre las capas límites térmica e hidrodinámica es aproxi-
madamente
@>@t = Pr1/3 (4.32)
Utilizando el mismo factor de corrección, Pr1/3, a cualquier distancia desde la super-
ficie, las curvas de la___
figura 4.13 se vuelven a trazar en la figura 4.14. La nueva
abscisa es Pr1/3(y/x)¥ Rex y la ordenada es la temperatura adimensional (T – Ts)/
(Tq - Ts), donde T es la temperatura local del fluido, Ts es la temperatura superficial
de la placa y Tq la temperatura de corriente libre. Esta modificación de la ordenada
produce los perfiles de temperatura para un intervalo amplio de números de Prandtl
en una sola línea, que es la curva para Pr = 1.
4.8.3 Evaluación del coeficiente de transferencia

de calor por convección
Ahora se pueden determinar la tasa de transferencia de calor por convección y el
coeficiente de transferencia de calor por convección. El gradiente de temperatura
adimensional en la superficie (en y = 0) es
0[(T - Ts)> (Tq - Ts)]
` = 0.332
0[(y> x)1Rex Pr1/3] y=0
Por tanto, en cualquier valor especificado de x,

0T Re1/2 1/3
`
x Pr
= 0.332 (Tq - Ts) (4.33)
0y y = 0 x
y la tasa local de transferencia de calor por convección por área unitaria se convierte,
al sustituir 0T/0y de la ecuación (4.33), en
0T Re1/2 1/3
`
x Pr
qcœœ = - k = - 0.332k (Tq - Ts) (4.34)
0y y = 0 x
La tasa total de transferencia de calor de una placa de ancho b y longitud L, obtenida

integrando qƎc de la ecuación (4.34) entre x = 0 y x = L, es
L Pr b(Ts - Tq)
q = 0.664k Re1/2 1/3 (4.35)
El coeficiente de transferencia de calor por convección es

qcœœ k
hcx = = 0.332 Re1/2Pr1/3 (4.36)
(Ts - Tq) x x
y el número de Nusselt local correspondiente es
hcxx
Nux = = 0.332 Re1/2
x Pr
1/3 (4.37)*
k
_
El número de Nusselt promedio hcL/k se obtiene integrando el lado derecho de _la
ecuación (4.36) entre x = 0 y x = L y dividiendo_el resultado entre L para obtener hc,
que es el valor promedio de hcx; multiplicando hc por L/k da
NuL = 0.664 Re1/2

L Pr
1/3 (4.38)
El valor promedio del número de Nusselt sobre una longitud L de la placa es,
por tanto, el doble del valor local de Nux = en x = L. Se puede verificar con facilidad
que la misma relación entre el valor promedio y el local es válida también para el
coeficiente de transferencia de calor, es decir,
_
hc = 2hc(x = L) (4.39)
En la práctica, las propiedades físicas en las ecuaciones (4.32) a (4.38) varían
con la temperatura, en tanto que para fines de análisis se supuso que las propiedades
físicas son constantes. Se han determinado datos experimentales que concuerdan
satisfactoriamente con los resultados anticipados analíticamente si las propiedades se
evalúan a una temperatura media a la mitad entre la de la superficie y la temperatura
de corriente libre; esta temperatura media se denomina temperatura de película.
EJEMPLO 4.3 Un colector solar de placa plana se coloca horizontalmente sobre un techo, como
se muestra en la figura 4.15. Para determinar su eficiencia, se necesita calcular la
pérdida de calor de su superficie al entorno. El colector es una franja larga de 1 ft
de ancho. La temperatura superficial del colector es 140 °F. Si un viento a 60 °F
sopla sobre el colector a una velocidad de 10 ft/s, calcule las cantidades siguientes
en x = 1 ft y x = xc en unidades inglesas y SI:
a) el espesor de la capa límite
b) el coeficiente de fricción local
c) el coeficiente de fricción promedio
d) el rozamiento local o esfuerzo cortante debido a la fricción
*Observe que la ecuación (4.37) se dedujo con la suposición de que Pr Ú 1. Por tanto, no es válida para
valores pequeños de Pr, es decir, metales líquidos. Una ecuación empírica para los números de Nusselt
locales para metales líquidos (Pr 6 0.1) se da en la tabla 4.5 en la página 283.
Colector 1 ft
Superficie del colector a 140 °F
x Aire a
60 °F
d
Viento x
FIGURA 4.15 Colector solar de placa plana del ejemplo 4.3.
e) espesor de la capa límite térmica

f) coeficiente de transferencia de calor por convección local
g) coeficiente de transferencia de calor por convección promedio
h) tasa de transferencia de calor por convección
SOLUCIÓN Las propiedades relevantes del aire a 100 °F en unidades inglesas son
r = 0.071 lbm/ft3
cp = 0.240 Btu/lbm °F
m = 1.285 * 10-5 lbm/ft s
k = 0.0154 Btu/h ft °F
Pr = 0.72
El número local de Reynolds en x = 1 ft (0.305 m) es
Uq rx (10 ft /s)(0.071 lbm /ft3)(1 ft)

Rex = 1 = = = 55200
m 1.285 * 10-5 lbm/ft s
y en x = 9 ft es
Rex=9 = 5 * 105
Suponiendo que el número crítico de Reynolds es 5 * 105, la distancia crítica es
Recm (5 * 105)(1.285 * 10-5 lbm /ft s)

xc = = = 9.0 ft (2.8 m)
Uqr (10 ft/s)(0.071 lbm /ft3)
Las cantidades deseadas se determinan sustituyendo los valores apropiados de la

variable en las ecuaciones pertinentes. Los resultados de los cálculos se muestran en
la tabla 4.4 y se sugiere que el lector los verifique.
4.9 Análisis integral aproximado de la capa límite 261
TABLA 4.4 Resultados del ejemplo 4.3

Ecuación Unidades Resultado Resultado Unidades Resultado Resultado
Parte Símbolo utilizada inglesas (x = 1 ft) (x = 9 ft) SI (x = 0.305 m) (x = 2.8 m)
a) @ (4.28) ft 0.0213 0.0638 m 0.00648 0.0195
b) Cfx (4.30) — 0.00283 0.000942 — 0.00283 0.000942
_
c) Cf (4.31) — 0.00566 0.00189 — 0.00566 0.00189
d) ts (4.29) lbf/ft2 3.12 * 10-4 1.04 10-4 N/m2 0.0149 0.00497
e) @t (4.32) ft 0.0237 0.0712 m 0.00723 0.0217
f) h_cx (4.36) Btu/h ft2 °F 1.08 0.359 W/m2 K 6.12 2.04
g) hc (4.39) 2
Btu/h ft °F 2.18 0.718 W/m2 K 12.23 4.08
h) q (4.35) Btu/h 172 517 W 50.5 152
Una relación útil entre el número de Nusselt local, Nux y el coeficiente de fric-
ción correspondiente, Cfx, se obtiene dividiendo la ecuación (4.37) entre RexPr1/3.
Cfx
a bPr2/3 =
Nux 0.332
= (4.40)
RexPr Re1/2
x 2
La relación adimensional Nux/RexPr se conoce como el número de Stanton, Stx.

De acuerdo con la ecuación (4.40), el número de Stanton por el número de Prandtl
elevado a la potencia de dos tercios es igual a la mitad del valor del coeficiente de
fricción. Esta relación entre transferencia de calor y fricción del fluido la propuso
Colburn [14] e ilustra la interrelación de los dos procesos.
4.9* Análisis integral aproximado de la capa límite

Para evitar los problemas implicados en la solución de ecuaciones diferenciales par-
ciales de la capa límite, se puede utilizar un enfoque integral. Para ese fin, considere
un volumen de control elemental que se extiende de la pared hasta más allá del límite
de la capa límite en la dirección y, tiene un espesor dx en la dirección x y tiene un
ancho unitario en la dirección z, como se muestra en la figura 4.16. Para obtener una
relación para el flujo neto de entrada de la cantidad de movimiento y el transporte
B
D U∞
y=d U∞
Capa límite
de velocidad
u(y)
A C x
dx
FIGURA 4.16 Volumen de control para el análisis integral

de la cantidad de movimiento.
de energía neto, se procede de una manera similar a la empleada para deducir las
ecuaciones de la capa límite en la sección anterior.
El flujo de la cantidad de movimiento a través de la cara AB en la figura 4.16 será:
d
ru2 dy
L0
De manera similar, el flujo de la cantidad de movimiento a través de la cara CD será

d d
d
ru2 dy + ru2 dy dx
L0 dx L0
El fluido también entra en el volumen de control a través de la cara BD a una tasa de

d
d
ru dy dx
dx L0
Esta cantidad es la diferencia entre la tasa de flujo que sale a través de la cara CD
y la del flujo que entra a través de la cara AB. Como el fluido que entra a través
de BD tiene una componente de la velocidad en la dirección x igual a la velocidad de
corriente libre Uq, el flujo de la cantidad de movimiento x hacia el volumen
de control a través de la cara superior es
d
d
Uq ru dy dx
dx L0
Sumando las componentes de la cantidad de movimiento x da
d d d
d d d
ru2 dy dx - Uq ru dy dx = - ru(Uq - u) dy
dx L0 dx L0 dx L0
No habrá cortante a través de la cara BD ya que esta cara está fuera de la capa límite,
donde du/dy es igual a cero. Sin embargo, existe un esfuerzo cortante tw que actúa
sobre las caras AB y CD. Escribiendo las fuerzas netas que actúan sobre el volumen
de control y sumándolas se obtiene la relación
pd - ap +
dp dp
dxbd - tw dx = - d dx - tw dx (4.41)
dx dx
Para flujo sobre una placa plana, el gradiente de presión en la dirección x se puede
ignorar y entonces la ecuación de la cantidad de movimiento se puede escribir en
la forma
d
d
ru(Uq - u)dy = tw (4.42)
dx L0
La ecuación integral de la energía se puede deducir de una manera similar. Sin

embargo, en este caso en la deducción se debe utilizar un volumen de control que se
extienda más allá de los límites tanto de la temperatura como de las capas límites de
y B
ys D
T∞
U∞
Espesor de la
Espesor de la capa límite
u ( y) capa límite hidrodinámica
T ( y) térmica, dt térmica, d
A C x
Ts
dx
FIGURA 4.17 Volumen de control para el análisis integral de conservación
de la energía.
la velocidad (consulte la figura 4.17). La primera ley de la termodinámica demanda
que se considere la energía en forma de entalpía, energía cinética y calor, así como
de trabajo de corte. Sin embargo, para velocidades bajas los términos de la energía
cinética y del trabajo de corte son pequeños comparados con las otras cantidades y
se pueden ignorar; entonces la tasa a la que entra la entalpía a través de la cara AB
está dada por
ys
cpruT dy
L0
en tanto que la tasa de flujo de salida de la entalpía a través de la cara CD es

ys s y
d
cpruT dy + cpruT dy dx
L0 dx L0
La entalpía transportada hacia el volumen de control a través de la cara superior está

dada por
s y
d
cpTs ru dy dx
dx L0
Por último, el calor se conducirá a través de la interfaz entre el fluido y la superficie

del sólido a la tasa
0T
- k dxa b
0y y = 0
Sumando todas las cantidades de energía se obtiene la ecuación integral para la

conservación de la energía en la forma
y y
0T
b
s s
d d
cpTq ru dy dx - rcpTu dy dx - k dxa = 0 (4.43)
dx L0 dx L0 0y y = 0
Sin embargo, se debe observar que fuera del límite de la capa límite térmica, la
temperatura es igual a la temperatura de corriente libre, Tq, de manera que la inte-
gración sólo se necesita calcular hasta y = @t. Por tanto, la ecuación (4.43) se puede
simplificar a la forma
d
0T
(T - T )u dy - a a b
t
d
= 0 (4.44)
dx L0 q 0y y = 0
que suele conocerse como la ecuación integral de la energía de la capa límite laminar
para flujo a baja velocidad.
4.9.1 Evaluación de la transferencia de calor y de los

coeficientes de fricción en flujo laminar
En el método integral aproximado, el primer paso es suponer contornos de velocidad
y temperatura en la forma de polinomios. Después, los coeficientes de los polino-
mios se evalúan para satisfacer las condiciones de frontera. Suponiendo un polinomio
de cuatro términos para la distribución de la velocidad [15]
u(y) = a + by + cy2 + dy3 (4.45)
las constantes se evalúan aplicando condiciones de frontera
en y = 0: u = 0 por tanto a = 0
0 2u
u = v = 0 por tanto = 0
0y2
0u
y = d: u = Uq y = 0
0y
Estas condiciones proporcionan cuatro ecuaciones para la evaluación de los cua-

tro coeficientes desconocidos en términos de la velocidad de corriente libre y del
espesor de la capa límite. Se puede verificar con facilidad que los coeficientes que
satisfacen estas condiciones de frontera son
3 Uq Uq
a = 0 b = c = 0 d =-
2 d 2d3
Sustituyendo estos coeficientes en la ecuación (4.45) y dividiendo todo entre la velo-
cidad de la corriente libre Uq para no dimensionalizar el resultado, se obtiene
1 y 3
- a b
u 3 y (4.46)
=
Uq 2 d 2 d
Sustituyendo la ecuación (4.46) para la distribución de la velocidad en la ecuación
integral de la cantidad de movimiento [ecuación (4.42)] se obtiene
d
- a b d # c1 -
1 y 3 1 y 3
c
2 3 y
+ a b ddy = tw = ma b
d 3 y du (4.47)
rUq
dx L0 2 d 2 d 2 d 2 d dy y = 0
El esfuerzo cortante en la pared tw se puede obtener evaluando el gradiente de

velocidad de la ecuación (4.46) en y = 0. Sustituyendo para tw y efectuando la inte-
gración en la ecuación (4.47) da
arUq b = m
d 2 39d 3 Uq (4.48)
dx 280 2 d
La ecuación (4.48) se puede reacomodar e integrar para obtener el espesor de la
capa límite en términos de la viscosidad, de la distancia desde el borde de ataque y
distribución de la velocidad de la corriente libre:
d2 140nx (4.49)
= + C
2 13Uq
Como @ = 0 en el borde de ataque (es decir, x = 0), el coeficiente C en la relación

anterior debe ser igual a 0 y
280nx
d2 =
13Uq
o
d 4.64 (4.50)
=
x Re1/2
x
Para evaluar el coeficiente de fricción, se sustituye la ecuación (4.46) en la ecuación

(4.47):
`
du 3 Uq
tw = m = m
dy y = 0 2 d
Sustituyendo para @ de la ecuación (4.50) da
3 mUq 1/2
tw = Rex
9.28 x
y el coeficiente de fricción Cfx es
tw 0.647
Cfx = = (4.51)
1 2
2 rUq
Re1/2
x
A continuación se recurre a la ecuación de energía y se propone una distri-

bución de temperatura en la capa límite de la misma forma que la distribución de
velocidad:
T(y) = e + fy + gy2 + hy3 (4.52)
Las condiciones de frontera para la distribución de la temperatura son en y = 0,
T = Ts; en y = @t (el espesor de la capa límite térmica), T = Tq y dT/dy = 0. Además,
de la ecuación (4.7b), d2T/dy2 en y = 0 debe ser cero debido a que u y v son cero en
la interfaz. A partir de estas condiciones, se deduce que las constantes son:
3 (Tq - Ts) (Tq - Ts)
e = Ts f = g = 0 h =
2 dt 2d3t
Si la variable en la ecuación de energía se toma como la temperatura en el fluido

menos la temperatura en la pared, la distribución de temperatura se puede escribir
en la forma adimensional
T - Ts 1 y 3
- a b
3 y (4.53)
=
Tq - Ts 2 dt 2 dt
Utilizando las ecuaciones (4.53) y (4.46) para T - Ts y u, respectivamente, la integral

en la ecuación (4.44) se puede escribir como
dt dt
(Tq - T)u dy = [(Tq - Ts) - (T - Ts)]u dy
L0 L0
dt
1 y 3 3 y 1 y 3
c1 - + a b dc - a b d dy
3 y
= (Tq - Ts)Uq
L0 2 dt 2 dt 2 d 2 d
Efectuando la multiplicación bajo el signo de integral, se obtiene
dt
a y6bdy
3 9 2 3 4 1 3 3 1
(Tq - Ts)Uq y - y + 3
y - 3
y + 3
y4 -
L0 2d 4ddt 4ddt 2d 4dtd 4d3t d3
que produce, después de integrar,
3 d2t 3 d2t 3 d2t 1 d4t 3 d4t 1 d4t

(Tq - Ts)Uq a - + - + - b
4 d 4 d 20 d 8 d3 20 d3 28 d3
Si se hace z = @t /@, la expresión anterior se puede escribir como
z b
3 2 3 4
(Tq - Ts)Uqda z -
20 280
Para fluidos que tienen un número de Prandtl igual a o mayor que la unidad, z es
igual a o menor que la unidad y el segundo término entre paréntesis se puede ignorar
comparado con el primero.* Sustituyendo esta forma aproximada para la integral en
la ecuación (4.44), se obtiene
0d 0T 3 Tq - Ts
`
3
Uq(Tq - Ts)z2 = a = a
20 0x 0y y = 0 2 dz
o
1 0d
U z3d = a
10 q 0x
*Para metales líquidos, que tienen Pr V 1, z 7 1 y el segundo término se puede ignorar.

4.10 Analogía entre la cantidad de movimiento y la transferencia de calor... 267
De la ecuación (4.50), se obtiene

0d n
d = 10.75
0x Uq
y con esta expresión se obtiene

10 a
z3 =
10.75 n
o
@t = 0.976@ Pr-1/3 (4.54)
Excepto para la constante numérica (0.976 comparada con 1.0), el resultado anterior
concuerda con el cálculo exacto de Pohlhausen [12].
La tasa de flujo de calor por convección de la placa por área unitaria es, de las
ecuaciones (4.1) y (4.53),
0T
`
3 k
qcœœ = - k = - (T - Ts)
0y y = 0 2 dt q
Sustituyendo en las ecuaciones (4.50) y (4.54) para @ y @t se obtiene
3 k Pr1/3Re1/2
x k
qœœ = - (Tq - Ts) = 0.33 Re1/2Pr1/3(Ts - Tq) (4.55)
2 x (0.976)(4.64) x x
y el número de Nusselt local, Nux, es
hcxx qcœœ x
Nux = = = 0.33Re1/2
x Pr
1/3
(4.56)
k (Ts - Tq) k
Este resultado concuerda muy bien con la ecuación (4.37), el resultado de un análisis
exacto de Pohlhausen [12].
El ejemplo anterior ilustra la utilidad del análisis aproximado de la capa límite.
Con un poco de visión e intuición física como guía, esta técnica produce resultados
satisfactorios sin la complicación matemática inherente en las ecuaciones exactas de
la capa límite. El método aproximado se ha aplicado a muchos otros problemas y los
resultados están disponibles en obras de consulta.
4.10* Analogía entre la cantidad de movimiento y la transferencia

de calor en flujo turbulento sobre una superficie plana
En la mayoría de aplicaciones prácticas, el flujo en la capa límite es turbulento en
lugar de laminar. De forma cuantitativa el mecanismo de intercambio en flujo tur-
bulento se puede representar como una magnificación del intercambio molecular
en flujo laminar. En flujo laminar en régimen permanente, las partículas de fluido
siguen líneas de corriente bien definidas. El calor y la cantidad de movimiento se trans-
fieren a través de las líneas de corriente sólo por difusión molecular y el flujo transver
sal es tan pequeño que cuando se inyecta un tinte de color en el fluido en algún
punto, éste sigue una línea de corriente sin mezclarse de manera apreciable. Por
otro lado, en flujo turbulento el tinte se distribuirá sobre un área amplia, una dis-
tancia corta, corriente abajo, desde el punto de la inyección. El mecanismo de mez-
clado consiste en remolinos con fluctuaciones rápidas que transportan partículas del
fluido de manera irregular. Los grupos de partículas chocan unos contra otros de
forma aleatoria, establecen un flujo transversal a una escala macroscópica y mez-
clan de manera efectiva el fluido. Como el mezclado en un flujo turbulento es a
una escala macroscópica con grupos de partículas transportados en una trayectoria
zigzagueante a través del fluido, el mecanismo de intercambio es muchas veces más
efectivo que en flujo laminar. Como resultado, las tasas de transferencia de calor
y de cantidad de movimiento en flujo turbulento y la fricción asociada y los coefi-
cientes de transferencia de calor son muchas veces mayores que en flujo laminar.
Si el flujo turbulento en un punto se promedia durante un periodo prolongado
(comparado con el periodo de una sola fluctuación), las propiedades medidas con
respecto al tiempo y la velocidad del fluido son constantes si el flujo promedio per-
manece en régimen permanente. Por tanto, es posible describir cada propiedad del
fluido y la velocidad en flujo turbulento en términos de un valor medio que no varía
con el tiempo y de una componente fluctuante que es una función del tiempo. Para
simplificar el problema, considere un flujo bidimensional (figura 4.18) en el que el
valor medio de la velocidad es paralelo a la dirección x. Entonces las componentes
de la velocidad instantánea u y v se pueden expresar en la forma
u = uq + uœ (4.57)
v = vœ
donde la barra sobre un símbolo denota el valor medio con respecto al tiempo y la
comilla simple denota la desviación instantánea del valor medio. De acuerdo con el
modelo utilizado para describir el flujo,
t*
1 (4.58)
uq = u dt
t* L0
donde t* es un intervalo de tiempo prolongado comparado con el periodo de las fluc-
tuaciones. En la figura 4.19 se muestra de manera cualitativa la variación del tiempo de
u y u¿. De la ecuación (4.58) o del __ análisis de la gráfica, es evidente que el __promedio
_____
del tiempo de u¿ es cero (es decir, u¿ = 0). Un argumento similar muestra que v¿ y ( pv¿)
también son cero.
y du
l
dy
+y ′
l
u(y)
+u′
du
l
dy
x
FIGURA 4.18 Longitud de mezclado para
transferencia de calor en flujo turbulento.
u u′
u u
t
FIGURA 4.19 Variación de la velocidad instantánea
con el tiempo en flujo turbulento.
Las componentes fluctuantes de la velocidad transportan masa de manera con-

tinua y en consecuencia cantidad de movimiento, a través de un plano normal a la
dirección y. La tasa instantánea de transferencia en la dirección y de la cantidad de
movimiento x por área unitaria en cualquier punto es
-(rv)¿(ū + u¿)
donde el signo menos, como se mostrará más adelante, toma en cuenta la correlación
estadística entre u¿ y v¿.
El promedio con respecto al tiempo de la transferencia de cantidad de movimiento
x da origen a un esfuerzo cortante turbulento aparente o esfuerzo de Reynolds tt,
definido por
t*
1
tt = - (rv)œ(uq + uœ) dt (4.59)
t* L0
Descomponiendo en partes este término, el promedio con respecto al tiempo de la

primera es
t*
1
(rv)œuq dt = 0
t* L0
_
puesto que u es una constante y el promedio con respecto al tiempo de (rv)¿ es cero.
Integrando el segundo término, la ecuación (4.59) se convierte en
t*
1
tt = - (rv)œuœ dt = - (rv)œuœ (4.60)
t* L0
o, si r es constante,
tt = - r(vœuœ) (4.61)
_____
donde (v¿u¿) es el promedio con respecto al tiempo del producto de u¿ y v¿.
No es difícil visualizar que los promedios con respecto al___tiempo de los pro-
ductos mezclados de las fluctuaciones de la velocidad, como v¿u¿, difieren de cero.
De la figura 4.18, se puede observar que las partículas que viajan hacia arriba
(v¿ 7 0) llegan a una capa en el fluido en la que la velocidad media ū es mayor que
la capa de donde provienen. Suponiendo que las partículas de fluido conservan, en
promedio, su velocidad original ū durante su migración, tenderán a aminorar la
velocidad de otras partículas de fluido después de que hayan llegado a su destino

y, de esa manera, dan origen a una componente negativa u¿. Por el contrario, si v¿
____ será positivo y viceversa.
es negativa, el valor observado de u¿ en el nuevo destino
Por tanto, el promedio con respecto al tiempo de u¿v¿ no es en promedio distinto
de cero, sino una cantidad negativa. El esfuerzo cortante turbulento definido por
la ecuación (4.61) es, por tanto, positivo y tiene el mismo signo que el esfuerzo
cortante laminar correspondiente.
duq duq
tl = m = rn
dy dy
Sin embargo, se debe observar que el esfuerzo cortante laminar es un esfuerzo

verdadero, en tanto que el esfuerzo cortante turbulento aparente es simplemente un
concepto introducido para tomar en cuenta los efectos de la transferencia de canti-
dad de movimiento por fluctuaciones turbulentas. Este concepto permite expresar el
esfuerzo cortante total en flujo turbulento como
fuerza viscosa flujo de cantidad de (4.62)

t = +
área movimiento turbulento
Para relacionar el flujo de la cantidad de movimiento turbulento con el gradiente

de velocidad promedio con respecto al tiempo d u>dy, se postula que las fluctuacio-
_
nes de las partículas macroscópicas del fluido en flujo turbulento son, en promedio,
similares al movimiento de las moléculas
_
en un gas [es decir, viajan, en prome-
dio, una distancia l perpendicular a u (figura 4.18) antes de llegar al reposo en otro
plano y]. Esta distancia l se conoce como longitud de mezclado de Prandtl [16, 17]
y corresponde cualitativamente a la trayectoria libre media de una molécula de gas.
Suponiendo que las partículas de fluido retienen su identidad y propiedades físicas
durante el movimiento transversal y que las fluctuaciones turbulentas se originan
principalmente de la diferencia en las propiedades medias con respecto al tiempo
entre los planos y espaciados una distancia l, si una partícula de fluido viaja de una
capa y a una capa y + l,
duq
u¿ M l (4.63)
dy
Con este modelo, el esfuerzo cortante turbulento tt en una forma análoga al esfuerzo
cortante laminar es
duq
tt = - rvœuœ = reM (4.64)
dy
donde el símbolo eM se denomina viscosidad eddy o coeficiente de intercambio

turbulento para la cantidad de movimiento. La viscosidad eddy eM es formalmente
análoga a la viscosidad cinemática v, pero, mientras que v es una propiedad física,
eM depende de la dinámica del flujo. Al combinar las ecuaciones (4.63) y (4.64)
se demuestra que eM = -v¿l, y la ecuación (4.62) da el esfuerzo cortante total en la
forma
duq
t = r(n + eM ) (4.65)
dy
y dT
l
dy
+υ′
l +qt
T(y)
+T′
dT
l
dy
x
FIGURA 4.20 Longitud de mezclado para transferencia
de energía en flujo turbulento.
En flujo turbulento, eM es mucho mayor que v y, por tanto, el término viscoso se

puede ignorar.
La transferencia de energía como calor en un flujo turbulento se puede repre-
sentar de manera análoga. Considere una distribución de temperatura media con res-
pecto al tiempo bidimensional, como se muestra en la figura 4.20. Las componentes
fluctuantes de la velocidad transportan partículas de fluido y la energía almacenada
en ellas a través de un plano normal a la dirección y. La tasa instantánea de transfe-
rencia de energía por área unitaria en cualquier punto en la dirección y es
(rvœ)(cpT) (4.66)
__
donde T = T + T ¿. Siguiendo la misma línea de razonamiento que condujo a la ecua-
ción (4.61), el promedio con respecto al tiempo de la transferencia de energía debida
a las fluctuaciones, tasa turbulenta de transferencia de calor, qt, es
qt = ArcpvœT œ (4.67)
Utilizando el concepto de la longitud de mezclado de Prandtl, se puede relacionar la
fluctuación de la temperatura con el gradiente de temperatura medio con respecto
al tiempo mediante la ecuación
dTq
Tœ M l (4.68)
dy
Físicamente esto significa que cuando una partícula de fluido migra de una capa y
a otra capa a una distancia l arriba o abajo, la fluctuación resultante de la tempera-
tura se ocasiona principalmente por la diferencia entre las temperaturas medias con
respecto al tiempo en las capas. Suponiendo que los mecanismos de transporte de
____
la temperatura (o energía) y la velocidad son similares, las longitudes de mezclado
en las ecuaciones (4.63) y (4.68) son iguales. Sin embargo, el producto v¿T ¿, es
positivo en promedio debido a que una v¿ positiva se acompaña por una T ¿, positiva
y viceversa.
Combinando las ecuaciones (4.67) y (4.68), la tasa turbulenta de transferencia
de calor por área unitaria se convierte en
qt dTq
qtœœ = = cprvœT œ = - cprvœl (4.69)
A dy
donde el signo menos es una consecuencia de la segunda ley de la termodinámica (con-

sulte el capítulo 1). Para expresar el flujo de calor turbulento en una forma análoga a la
ecuación de conducción de Fourier, se define eH, una cantidad denominada coeficien-
te de intercambio turbulento__ para temperatura o difusividad
__ eddy de transferencia de
calor por la ecuación eH = v¿l. Sustituyendo eH por v¿l en la ecuación (4.69) da
dTq
qtœœ = - cpreH (4.70)
dy
La tasa total de transferencia de calor por área unitaria normal a la velocidad de
corriente media entonces se puede escribir como
q conducción molecular transferencia turbulenta
qœœ = = +
A área unitaria área unitaria
o en forma simbólica como
dTq
qœœ = - cpr(a + eH) (4.71)
dy
donde a = k>cpr, que es la difusividad molecular de transferencia de calor. La contri-
bución a la transferencia de calor por conducción molecular es proporcional a a, y la
contribución turbulenta es proporcional a eH. Para todos los fluidos excepto metales
líquidos, eH es mucho mayor que a en flujo turbulento. La relación de la viscosidad
cinemática molecular a la difusividad molecular del calor, v>a, se ha nombrado con
anterioridad número de Prandtl. De manera similar, la relación de la viscosidad
eddy turbulenta a la difusividad eddy, eM > eH se podría considerar como un número
de Prandtl turbulento Prt. De acuerdo con la teoría de la longitud de mezclado
___ de
Prandtl, el número de Prandtl turbulento es la unidad ya que eM= eH = v¿l.
Si bien este tratamiento del flujo turbulento está demasiado simplificado, resul-
tados experimentales indican que al menos es cualitativamente correcto. Isakoff y
Drew [18] determinaron que Prt, para el calentamiento del mercurio en flujo turbu-
lento dentro de un tubo puede variar de 1.0 a 1.6; Forstall y Shapiro [19] determi-
naron que Prt es aproximadamente de 0.7 para gases. Estos últimos investigadores
también demostraron que Prt es sustancialmente independiente del valor del número
de Prandtl laminar así como del tipo de experimento. Suponiendo que Prt es igual a
la unidad, el flujo de calor turbulento se puede relacionar con el esfuerzo cortante
turbulento combinando las ecuaciones (4.64) y (4.70):
dTq
qtœœ = - tt cp (4.72)
duq
Esta relación originalmente se dedujo en 1874 por el científico británico Osborn

Reynolds y se denomina analogía de Reynolds. Es una buena aproximación cuando
el flujo es turbulento y se puede aplicar a capas límite turbulentas así como a
flujo turbulento en tubos o conductos. Sin embargo, la analogía de Reynolds no
es válida en la subcapa viscosa [20]. Puesto que esta capa ofrece una resistencia
térmica grande al flujo de calor, la ecuación (4.72) en general, no es suficiente para
describir una solución cuantitativa. Sólo puede emplearse directamente para fluidos
que tengan un número de Prandtl de la unidad para calcular la tasa de transferencia
de calor. Ahora se considerará este caso especial.
4.11 Analogía de Reynolds para flujo turbulento sobre superficies planas 273
4.11 Analogía de Reynolds para flujo turbulento sobre superficies planas

Para deducir la relación entre la transferencia de calor y la fricción superficial en
flujo sobre una superficie plana para un número de Prandtl de la unidad, se recuerda
que el esfuerzo cortante laminar t es
du
t = m
dy
y la tasa de flujo de calor por área unitaria a través de cualquier plano perpen-
dicular a la dirección y es
dT
qœœ = - k
dy
Combinando estas ecuaciones se produce

k dT
qœœ = - t (4.73)
m du
Al examinar las ecuaciones (4.72) y (4.73) se observa que si cp = k>m (es decir,
para Pr = 1), la misma ecuación de flujo de calor se aplica en las capas laminar y
turbulenta.
Para determinar la tasa de transferencia de calor de una placa plana a un fluido
con Pr = 1 fluyendo sobre ella en flujo turbulento, se remplaza k>m por cp se sepa-
ran las variables en la ecuación (4.73). Suponiendo que qƎ y t son constantes, se
obtiene
qsœœ
du = - dT (4.74)
tscp
donde el subíndice s se utiliza para indicar que tanto qƎ como t se consideran en la
superficie de la placa. Integrando la ecuación (4.74) entre los límites u = 0 cuando
T = Ts y u = Uq cuando T = Tq da
qsœœ
U = (Ts - Tq) (4.75)
tscp q
Pero como por definición los coeficientes de transferencia de calor local y de fric-
ción son
qsœœ rU2q
hcx = y tsx = Cfx
(Ts - Tq) 2
La ecuación (4.75) se puede escribir así:
hcx Nux Cfx

= = (4.76)
cprUq RexPr 2
La ecuación (4.76) es satisfactoria para gases en los que Pr es aproximadamente

igual a la unidad. Colburn [14] demostró que la ecuación (4.76) también se puede
emplear para fluidos con números de Prandtl que varíen de 0.6 a aproximadamente
50 si se modifica de acuerdo con resultados experimentales para que sea
Nux Cf x
Pr2/3 = StxPr2/3 = (4.77)
RexPr 2
donde el subíndice x denota la distancia desde el borde de ataque de la placa.
En la práctica para aplicar la analogía entre transferencia de calor y transfe-
rencia de cantidad de movimiento, se necesita conocer el coeficiente de fricción
Cfx. Para flujo turbulento sobre una superficie plana, la ecuación empírica para el
coeficiente de fricción local
Uq x -1/5
Cf x = 0.0576a b (4.78a)
n
concuerda bien con resultados experimentales en el número de Reynolds variando

entre 5 * 105 y 107 siempre que no ocurra separación. Suponiendo que la capa límite
turbulenta empieza en el borde de ataque, el coeficiente de fricción promedio sobre
una superficie plana de longitud L se puede obtener integrando la ecuación (4.78a):
UqL -1/5
L
b
1
Cqf = Cfx dx = 0.072a (4.78b)
L L0 n
Utilizando el método integral de la sección 4.9 combinado con datos experimentales

para el coeficiente de fricción, también se puede deducir una expresión para el espe-
sor de la capa límite en flujo turbulento. De acuerdo con Schlichting [1], el espesor
de la capa límite hidrodinámica se puede aproximar mediante la relación
(d/x) = 0.37/Re0.2
x (4.79)
Al comparar las ecuaciones (4.78a) y (4.79) con los resultados para flujo laminar, las
ecuaciones (4.30) y (4.28), es evidente que el decaimiento del coeficiente de fricción
con la distancia es más gradual en flujo turbulento que en flujo laminar (x-0.2 contra
x-0.5), en tanto que el espesor de la capa límite aumenta más rápidamente en flujo
turbulento que en flujo laminar (@t r x0.8 contra @ r x0.5).
En flujo turbulento, el crecimiento de la capa límite se afecta más por fluctua-
ciones aleatorias en el fluido que por difusión molecular. De aquí que el número de
Prandtl no tiene un papel importante en el proceso y para Pr 7 0.5, la ecuación (4.79)
también es una aproximación razonable para determinar el espesor de la capa límite
térmica, es decir, @ « @t en flujo turbulento [21].
4.12 Capa límite mezclada

En realidad, una capa límite laminar precede a la capa límite turbulenta entre x = 0 y
x = xc. Como la resistencia de rozamiento local de una capa límite laminar es menor
que la resistencia de rozamiento local de una capa límite turbulenta en el mismo
número de Reynolds, el rozamiento promedio calculado con la ecuación (4.78b) sin
corregir para la parte laminar de la capa límite es demasiado grande. Sin embargo,
el rozamiento real se puede estimar bastante bien suponiendo que detrás del punto
4.12 Capa límite mezclada 275
de transición, la capa límite turbulenta se comporta como si hubiera empezado en

el borde de ataque.
Sumando la resistencia de rozamiento por fricción laminar entre x = 0 y x = xc
con el rozamiento turbulento entre x = xc y x = L da, por ancho unitario,
[0.072ReL-1/5L - 0.072Rex-1/5 xc + 1.33Rex-1/2 xc]

Cqf =
c c
L
Para un número de Reynolds crítico de 5 * 105, esto da
b
0.0464xc
Cqf = 0.072aReL-1/5 - (4.80)
L
Sustituyendo la ecuación (4.78a) para Cfx en la ecuación (4.77) se obtiene el

número de Nusselt local en cualquier valor de x mayor que xc:
Uqx 0.8
= 0.0288Pr1/3 a b
hcxx
Nux = (4.81)
k n
Se observa que el coeficiente de transferencia de calor local hcx para transferencia
de calor por convección a través de una capa límite turbulenta disminuye con la
distancia x como hcx r 1>x0.2. La ecuación (4.81) muestra que en comparación con
el flujo laminar, donde hcx r 1>x1/2, el coeficiente de transferencia de calor en flujo
turbulento disminuye rápidamente con x y que el coeficiente de transferencia de
calor turbulento es mucho mayor que el coeficiente de transferencia de calor laminar
a un valor dado del número de Reynolds.
El coeficiente de transferencia de calor promedio en flujo turbulento sobre una
superficie plana de longitud L se puede calcular hasta una primera aproximación
integrando la ecuación (4.81) entre x = 0 y x = L:
L
1
hqc = hcx dx
L L0
En forma adimensional, se obtiene
hqcL
NuL = = 0.036Pr1/3Re0.8
L (4.82)
k
En la ecuación (4.82) se ignora la existencia de la capa límite laminar y, por tanto, es
válida sólo cuando L » xc. La capa límite laminar se puede incluir en el análisis si la
ecuación (4.56) se utiliza entre x = 0 y x = xc y si la ecuación (4.81) se utiliza entre
x = xc y x = L para la integración de hcx. Esto produce, con Rec = 5 * 105,
NuL = 0.036Pr1/31Re0.8
L - 23200 2 (4.83)
EJEMPLO 4.4 El cárter de un automóvil tiene dimensiones aproximadas de 0.6 m de longitud, 0.2 m
de ancho y 0.1 m de profundidad (consulte la figura 4.21). Suponiendo que la tempe-
ratura superficial del cárter es 350 K, estime la tasa de flujo de calor del cárter al aire
0.2 m
0.6 m
Parte inferior
Aire a del cárter, 350 K
30 m/s
FIGURA 4.21 Cárter automotriz del ejemplo 4.4.
atmosférico a 276 K a una velocidad de marcha de 30 m/s. Suponga que las vibraciones
del motor y del chasís inducen la transición de flujo laminar a turbulento tan cerca del
borde de ataque que, para fines prácticos, la capa límite es turbulenta sobre toda la
superficie. Ignore la radiación y utilice para las superficies frontal y posterior el mismo
coeficiente de transferencia de calor por convección promedio que para los lados infe-
rior y laterales.
SOLUCIÓN Utilizando las propiedades del aire a 313 K, el número de Reynolds es
rUqL 1.092 * 30 * 0.6

ReL = = = 1.03 * 106
m 19.123 * 10-6
De la ecuación (4.82), el número de Nusselt promedio es
NuL = 0.036Pr1/3Re0.8
L
= 0.036(0.71)1/3(1.03 * 106)0.8
= 2075
y el coeficiente de transferencia de calor por convección promedio se convierte en
NuLk 2075 * (0.0265 W/mK)

hqc = = = 91.6 W/m2K
L (0.6 m)
El área superficial que disipa calor es 0.28 m2 y, por tanto, la tasa de pérdida de
calor del cárter es
qc = hqc A(Ts - Tq) = (91.6 W/m2K)(0.28 m2)(350 - 276)(K)

= 1898 W
4.13 Condiciones de frontera especiales y flujo a alta velocidad 277
4.13* Condiciones de frontera especiales y flujo a alta velocidad

En las secciones anteriores de este capítulo, la condición de frontera en la superficie
supuso una temperatura constante y uniforme. Sin embargo, existen otras situacio-
nes de interés práctico, por ejemplo, una superficie con flujo de calor constante y
una placa plana con una longitud inicial sin calentar seguida de una temperatura
superficial uniforme mayor que la temperatura ambiente.
Las condiciones de frontera para una placa plana con flujo sobre su superficie y
una longitud inicial sin calentar se ilustra en la figura 4.22. La capa límite hidrodiná-
mica inicia en x = 0, en tanto que la capa límite térmica comienza en x = z, donde la
temperatura superficial aumenta de Tq a Ts. No ocurre transferencia de calor entre
0 … x 6 z. Para este caso, suponiendo flujo laminar, el método integral da
Nux/Nux,z=0 = {1 - (z/x)3/4}-0.33 (4.84)
donde Nux,z=0 se determina a partir de la ecuación (4.56).

Para flujo turbulento con una longitud inicial sin calentar,
Nux/Nux,z=0 = {1 - (z/x)0.9}-0.9 (4.85)
donde Nux,z=0 se obtiene con la ecuación (4.81).

Cuando el flujo ocurre sobre una placa plana con un flujo de calor constante
impuesto (por ejemplo, por calentamiento eléctrico), el número de Nusselt en flujo
laminar está dado por
Nux = 0.453Re0.5
x Pr
0.33
(4.86)
en tanto que para flujo turbulento
Nux = 0.0308Re0.8
x Pr
0.33
(4.87)
Para una longitud inicial sin calentar con flujo de calor constante en x 7 z, las ecua-
ciones (4.86) y (4.87) se pueden modificar utilizando factores de corrección iguales
en los lados derechos de las ecuaciones (4.84) y (4.85), respectivamente.
Capa límite hidrodinámica

δ
U∞, T∞
Capa límite térmica

δt
Ts = T∞ Ts > T∞
ζ qsʹʹ
x
FIGURA 4.22 Placa paralela en flujo paralelo con longitud inicial
sin calentar.
Con flujo de calor constante, la temperatura superficial no es uniforme, pero una

vez que se conoce el coeficiente de transferencia de calor, la temperatura superficial
se puede obtener con
Ts(x) = Tq + (qs"/hx) (4.88)
en cualquier ubicación dada.
Otra condición superficial especial se tiene en el flujo a alta velocidad. La
transferencia de calor por convección en flujo a alta velocidad es importante en
aeronaves y misiles cuando la velocidad se aproxima o excede la del sonido. Para
un gas perfecto la velocidad del sonido a se puede obtener de
g uT
a = (4.89)
C
donde g = relación de calores específicos (casi 1.4 para aire)

u = constante universal de los gases
T = temperatura absoluta
= peso molecular del gas
Cuando la velocidad de un gas que fluye sobre una superficie calentada o enfriada
es del orden de la velocidad del sonido o mayor, el campo de flujo ya no se puede
describir únicamente en términos del número de Reynolds; en cambio, también
se debe considerar la relación de la velocidad del flujo de gas a la velocidad acús-
tica: el número de Mach, M = Uq>aq. Cuando la velocidad del gas en sistema de
flujo alcanza un valor de aproximadamente la mitad de la velocidad del sonido,
los efectos de la disipación viscosa en la capa límite se vuelven importantes. En
esas condiciones, la temperatura de la superficie sobre la cual fluye un gas, en
realidad puede exceder la temperatura de corriente libre. Para flujo sobre una
superficie adiabática, como una pared perfectamente aislada, la figura 4.23 mues-
tra de manera esquemática las distribuciones de velocidad y temperatura. La alta
temperatura en la superficie es el resultado combinado del calentamiento debido
a la disipación viscosa y al aumento de temperatura del fluido cuando la energía
cinética del flujo se convierte en energía interna, mientras que el flujo desacelera
a través de la capa límite. La forma real del perfil de temperatura depende de
la relación entre la tasa a la que el trabajo cortante viscoso aumenta la energía
interna del fluido y de la tasa a la que el calor se conduce hacia la corriente libre.
T∞
U∞
Perfil de velocidad
Perfil de temperatura
Superficie de la placa
Tas
FIGURA 4.23 Distribuciones de

velocidad y temperatura en flujo a
alta velocidad sobre una placa
aislada.
Aunque los procesos en una capa límite a alta velocidad no son adiabáticos, es
práctica general relacionarlos con procesos adiabáticos. La conversión de energía ciné-
tica en un gas que se desacelera adiabáticamente a una velocidad cero se describe por
U2q
i0 = iq + (4.90)
2gc
donde i0 es la entalpía de estancamiento e iq es la entalpía del gas en la corriente
libre. Para un gas ideal, la ecuación (4.90) se convierte en
U2q
T0 = Tq + (4.91)
2gccp
o, en términos del número de Mach,
T0 g - 1 2 (4.92)
= 1 + Mq
Tq 2
donde T0 es la temperatura de estancamiento y Tq es la temperatura de corriente
libre.
En una capa límite real, el fluido no se lleva al reposo reversiblemente debido a
que el proceso de cortante viscoso es termodinámicamente irreversible. Para tomar
en cuenta la irreversibilidad en un flujo de una capa límite, se define un factor de
recuperación r como
Tas - Tq
r = (4.93)
T0 - Tq
donde Tas es la temperatura superficial adiabática.
Mediante experimentos [22] se ha demostrado que en flujo laminar:
r = Pr1/2 (4.94)
en tanto que en flujo turbulento:
r = Pr1/3 (4.95)
Cuando una superficie no está aislada, la tasa de transferencia de calor por con-
vección entre un gas a alta velocidad y esa superficie está gobernada por la relación
0T
qcœœ = - k `
0y y = 0
La influencia de la transferencia de calor hacia y desde la superficie en la distribución

de la temperatura se ilustra en la figura 4.24. Se observa que en flujo a alta velocidad el calor
se puede transferir hacia la superficie incluso cuando la temperatura superficial esté arriba
de la temperatura de corriente libre. Este fenómeno es el resultado del cortante viscoso
y con frecuencia se denomina calentamiento aerodinámico. La tasa de transferencia de
calor en flujo a alta velocidad sobre una superficie plana se puede predecir [1] a partir
de la ecuación de energía de la capa límite:
0T 0T 0 2T m 0u 2
u + v = a 2 + a b
0x 0y 0y rcp 0x
donde en el último término se toma en cuenta la disipación viscosa. Sin embargo, para
la mayoría de los fines prácticos, la tasa de transferencia de calor se puede calcular
T0
T∞ Perfil de temperatura
Perfil de para flujo de calor desde
temperatura para la superficie
flujo de calor hacia
la superficie
Tas
FIGURA 4.24 Perfiles de temperatura en una
capa límite a alta velocidad para calentamiento
y enfriamiento.
con las mismas relaciones utilizadas para flujo a baja velocidad si el coeficiente de
transferencia de calor por convección promedio se redefine mediante la relación
_
qcƎ = hc(Ts - Tas) (4.96)
que producirá flujo de calor cero cuando la temperatura superficial Ts sea igual a la
temperatura superficial adiabática.
Puesto que en el flujo a alta velocidad los gradientes de temperatura en una
capa límite son grandes, las variaciones en las propiedades físicas del fluido tam-
bién serán sustanciales. Sin embargo, Eckert [23] demostró que las ecuaciones de
transferencia de calor de propiedades constantes aún se pueden utilizar si todas las
propiedades se evalúan a una temperatura de referencia T* dada por
T* = Tq + 0.5(Ts - Tq) + 0.22(Tas - Tq) (4.97)
Los valores locales del coeficiente de transferencia de calor, definidos por
qcœœ
hcx = (4.98)
Ts - Tas
se pueden obtener a partir de las ecuaciones siguientes:

Capa límite laminar (Re*x 6 105 ):
St*x = a b = 0.332(Re*x )-1/2(Pr*)-2/3

hcx *
(4.99)
cprUq
Capa límite turbulenta (105 6 Rex* 6 107 ):
St*x = a b = 0.0288(Re*)
hcx * -1/5
x (Pr*)-2/3 (4.100)
cprUq
Capa límite turbulenta (107 6 Re*x 6 109 ):
St*x = a b =
hcx * 2.46
(Pr*)-2/3 (4.101)
cprUq (ln Rex*)2.584
La ecuación (4.101) se basa en datos experimentales de los coeficientes de fricción

locales en flujo de gas a alta velocidad [23] en el intervalo del número de Reynolds
de 107 a 109 que se correlacionan mediante
4.92
Cfx = (4.102)
(ln Re*x)2.584
Si se tiene que determinar el valor promedio de los coeficientes de transferencia
de calor, las expresiones anteriores se tienen que integrar entre x = 0 y x = L. Sin
embargo, la integración puede que se tenga que efectuar numéricamente en la mayo-
ría de los casos prácticos ya que la temperatura de referencia T* no es la misma para
las partes laminar y turbulenta de la capa límite, como se demuestra mediante las
ecuaciones (4.94) y (4.95).
Cuando la velocidad de un gas es extremadamente alta, la capa límite se vuelve
tan caliente que el gas comienza a desasociarse. En esas situaciones, Eckert [23]
recomienda que el coeficiente de transferencia de calor se apoye en la diferencia de
entalpía entre la superficie y el estado adiabático, (is - ias) y se puede definir por
qcƎ = hci(is - ias) (4.103)
Si se define un factor de recuperación de la entalpía por

ias - iq
ri =
i0 - iq
se puede utilizar la misma relación empleada antes para calcular la temperatura de

referencia para calcular una entalpía de referencia:
i* = iq + 0.5(is - iq) + 0.22(ias - iq) (4.104)
Entonces el número local de Stanton se redefine como
hc,i
x, t =
St* (4.105)
r*uq
y se utiliza en las ecuaciones (4.99), (4.100) y (4.101) para calcular el coeficiente de

transferencia de calor. Se debe observar que las entalpías en las relaciones anteriores
son los valores totales, que incluyen la energía química de disociación así como la
energía interna. Como se demuestra en [23], este método de cálculo concuerda bas-
tante bien con datos experimentales.
En algunas situaciones, por ejemplo, a altitudes extremadamente altas, la densi-
dad del fluido puede ser tan pequeña que la distancia entre las moléculas de gas se
vuelve del mismo orden de magnitud que la capa límite. En esos casos, el fluido no
se puede tratar como un continuo y se necesitan subdividir los procesos de flujo en
regímenes. Estos regímenes de flujo están caracterizados por la relación de la tra-
yectoria libre molecular y una escala física significativa del sistema; esta relación se
denomina número de Knudsen, Kn. Un flujo continuo corresponde a valores pequeños
de Kn, en tanto que a valores mayores de Kn, los choques moleculares ocurren prin-
cipalmente en la superficie y en la corriente principal. Como el transporte de energía es
por el movimiento libre de moléculas entre la superficie y la corriente principal,
este régimen se denomina régimen de “moléculas libres”. Entre los regímenes de
100
Satélite de 3 m a una
altitud de 160 km
Número de Mach
10
lib jo
Transición nte
s
as lu
liza
re
ul f
1.0 es
éc de
jo d
ol n
flu Aeronave
m me
e
nd
de égi
ime comercial
R
0.1 g
Ré
Región de flujo continuo
0.01
0.001 0.01 0.1 1 10 102 103 104 105 106 107
Número de Reynolds
FIGURA 4.25 Regímenes de flujo a alta velocidad.
moléculas libres y continuo existe un intervalo de transición, denominado régimen de

“flujo deslizante” debido a que se trata suponiendo que la temperatura y la velocidad
se “deslizan” en las interfaces fluido-sólido. En la figura 4.25 se muestra un mapa
de estos regímenes. Para consultar un tratamiento de transferencia de calor y fricción
en estos sistemas de flujo especializados al lector se le refiere a [24-27].
4.14 Comentarios finales

En este capítulo, se estudiaron los principios de la transferencia de calor por con-
vección forzada. Se vio que la transferencia de calor por convección está muy rela-
cionada con la mecánica del flujo de fluido, en particular con el flujo en la vecindad
de la superficie de transferencia. También se observó que la naturaleza de la trans-
ferencia de calor, así como los fenómenos de flujo, depende en gran medida de si el
fluido muy alejado de la superficie está en flujo laminar o en turbulento.
Para familiarizarse con los principios básicos de la teoría de la capa límite y
de la transferencia de calor por convección forzada, se consideró el problema de
convección en flujo sobre una placa plana con cierto detalle. Este sistema es geomé-
tricamente simple, pero ilustra las características más importantes de la convección
forzada. En capítulos subsiguientes, se abordará la transferencia de calor por con-
vección en sistemas geométricamente más complicados. En el capítulo siguiente,
se examinan los fenómenos de convección natural. En el capítulo 6 se abordará
la transferencia de calor por convección hacia y desde fluidos circulando dentro
de tubos y conductos. En el capítulo 7 se considerará la convección forzada en flujo
sobre las superficies exteriores de cuerpos como cilindros, esferas, tubos y paquetes
de tubos. La aplicación de los principios de la transferencia de calor por convección
forzada para el diseño y selección de equipo de transferencia de calor se abordará
en el capítulo 8.
Para la conveniencia del lector, en la tabla 4.5 se presenta un resumen de las
ecuaciones utilizadas para calcular la transferencia de calor y los coeficientes de
fricción en flujo a baja velocidad de gases y líquidos sobre superficies planas o
planas sólo ligeramente curvas. Para consultar información adicional y correla-
ciones o ecuaciones para otras condiciones y geometrías, al lector se le refiere a
[3, 18, 21 y 28].
Referencias 283
TABLA 4.5 Resumen de ecuaciones empíricas útiles para calcular factores de fricción y coeficientes de trans-
ferencia de calor en flujo sobre superficies planas o ligeramente curvas con un ángulo de ataque igual a ceroa
Coeficiente Ecuación Condiciones
Flujo laminar
Coeficiente de fricción local Cfx = 0.664Re-0.5
x Rex 6 5 * 105
Número de Nusselt local a una distancia Nux = 0.332Re0.5
x Pr
0.33 Pr 7 0.5, Rex 6 5 * 105
x del borde de ataque Nux = 0.565(RexPr)0.5 Pr 6 0.1, Rex 6 5 * 105
_
Coeficiente de fricción promedio Cf = 1.33Re-0.5 ReL 6 5 * 105
___ L
Número de Nusselt promedio entre NuL = 0.664Re0.5
L Pr
0.33 Pr 7 0.5, ReL 6 5 * 105
x=0 y x=L
Flujo turbulento
Coeficiente de fricción local
Número de Nusselt local a una distancia
x del borde de ataque
Cfx = 0.0576Re-0.2
x
Nux = 0.0288Re0.8
x Pr
0.33 } Rex 7 5 * 105, Pr 7 0.5
_
Cf = 0.072[Re-0.2 - 0.0464(xc /L)]
Coeficiente de fricción promedio
Número de Nusselt promedio entre
x = 0 y x = L con transición en
___ L
NuL = 0.036Pr 0.33[Re0.8

L - 23,2000] } ReL 7 5 * 105, Pr 7 0.5
Rex,c = 5 * 105
a Aplicablea flujo a baja velocidad (número de Mach 60.5) de gases y líquidos con todas las propiedades físicas evaluadas a la temper-
atura pelicular media, Tf = (Ts + Tq)/2.
L
Cfx = ts >1rU2q/2gc2 Cqf = (1/L) Cfx dx Pr = cpm/k
L0
L
Nux = hcx/k Nu = hqcL/k hqc = (1/L) hc(x)dx
L0
Rex = rUqx/m ReL = rUqL/m
Referencias
1. H. Schlichting, Boundary Layer Theory, 6a. ed., J. 6. E. R. Van Driest, “On Dimensional Analysis and the
Kestin, transl., McGraw-Hill, Nueva York, 1968. Presentation of Data in Fluid Flow Problems”, J. Appl.
2. E. R. Van Driest, “Calculation of the Stability of the Mech., vol. 13, p. A-34, 1940.
Laminar Boundary Layer in a Compressible Fluid on 7. Handbook of Heat Transfer, W. M. Rohsenow, J. P.
a Flat Plate with Heat Transfer”, J. Aero. Sci., vol. 19, Hartnett, y Y. I. Cho, eds., McGraw-Hill, Nueva York,
pp. 801-813, 1952. 1998.
3. S. Kakaç, R. K. Shah y W. Aung, Handbook of Single- 8. S. V. Patankar y D. B. Spalding, Heat and Mass
Phase Convective Heat Transfer, Wiley, Nueva York, Transfer in Boundary Layers, 2a. ed., International
1987. Textbook Co., Londres, 1970.
4. L. D. Landau y E. M. Lifshitz, Fluid Mechanics, 9. P. Majumdar, Computational Methods for Heat and Mass
Pergamon Press, Nueva York, 1959. Transfer, Taylor & Francis, Nueva York, NY, 2005.
5. H. L. Langhaar, Dimensional Analysis and Theory of 10. M. Blasius, “Grenzschichten in Flüssigkeiten mit Kleiner
Models, Wiley, Nueva York, 1951. Reibung”, Z. Math. Phys., vol. 56, núm. 1, 1908.
11. M. Hansen, “Velocity Distribution in the Boundary 20. R. C. Martinelli, “Heat Transfer to Molten Metals”,
Layer of a Submerged Plate”, NACA TM 582, 1930. Trans. ASME, vol. 69, pp. 947-959, 1947.
12. E. Pohlhausen, “Der Wärmeaustausch zwischen festen 21. W. Kays, M. Crawford y B. Weigand, Convective Heat
Körpern und Flüssigkeiten mit kleiner Reibung und and Mass Transfer, 4a. ed., McGraw-Hill, Nueva York,
kleiner Wärmeleitung”, Z. Angew. Math. Mech., vol. NY, 2005.
1, p. 115, 1921. 22. J. Kaye, “Survey of Friction Coefficients, Recovery
13. B. Gebhart, Heat Transfer, 2a. ed., McGraw-Hill, Factors, and Heat Transfer Coefficients for Supersonic
Nueva York, 1971. Flow”, J. Aeronaut. Sci., vol. 21, núm. 2, pp. 117-229,
14. A. P. Colburn, “A Method of Correlating Forced 1954.
Convection Heat Transfer Data and a Comparison with 23. E. R. G. Eckert, “Engineering Relations for Heat
Fluid Friction”, Trans. AIChE, vol. 29, pp. 174 -210, Transfer and Friction in High-Velocity Laminar and
1993. Turbulent. Boundary Layer Flow over Surfaces with
15. E. R. G. Eckert y R. M. Drake, Heat and Mass Transfer, Constant Pressure and Temperature”, Trans. ASME, vol.
2a. ed., McGraw-Hill, Nueva York, 1959. 78, pp. 1273-1284, 1956.
16. L. Prandtl, “Bemerkungen über den Wärmeübergang 24. E. R. Van Driest, “Turbulent Boundary Layer in
irn Rohr”, Phys. Zeit., vol. 29, p. 487, 1928. Compressible Fluids”, J. Aeronaut. Sci., vol. 18, núm. 3,
17. L. Prandtl, “Eine Beziehung zwischen Wärmeaustauch pp. 145-161, 1951.
und Ströhmungswiederstand der Flüssigkeiten”, Phys. 25. A. K. Oppenheim, “Generalized Theory of Convective
Zeit., vol. 10, p. 1072, 1910. Heat Transfer in a Free-Molecule Flow”, J. Aeronaut.
18. S. E. Isakoff y T. B. Drew, “Heat and Momentum Sci., vol. 20, pp. 49-57, 1953.
Transfer in Turbulent Flow of Mercury”, Institute of 26. W. D. Hayes y R. F. Probstein, Hypersonic Flow
Mechanical Engineers and ASME, Proceedings, General Theory, Academic Press, Nueva York, 1959.
Discussion on Heat Transfer, pp. 405-409, 1951. 27. F. M. White, Viscous Fluid Flow, 2a. ed., McGraw-
19. W. Forstall, Jr., y A. H. Shapiro, “Momentum and Mass Hill, Nueva York, 1991.
Transfer in Co-axial Gas Jets”, J. Appl. Mech., vol. 17, 28. R. B. Bird, W. E. Stewart y E. N. Lightfoot, Transport
p. 399, 1950. Phenomena, 2a. ed., Wiley, Nueva York, NY, 2007.
Problemas
Los problemas de este capítulo están organizados por tema 4.2 Evalúe el número de Prandtl a partir de los datos siguien-
como se muestra a continuación. tes: cp = 0.5 Btu/lbm °F, k = 2 Btu/ft h °F, m = 0.3 lbm/ft s.
4.3 Evalúe el número de Nusselt para flujo sobre una esfera
para
_ las condiciones siguientes: D = 6 in, k = 0.2 W/m K,
Tema Problema número h c= 18 Btu/h ft2 °F.
Números adimensionales 4.1-4.6 4.4 Evalúe el número de Stanton para flujo sobre un tubo
Análisis dimensional 4.7-4.19 a partir de los datos siguientes: D = 10 cm, Uq = 4 m/s,
_ = 13 000 kg/m , m = 1 * 10 N s/m , cp = 140 J/kg K,
r 3 -3 2
Capas límites 4.20-4.28
hc= 1000 W/m K. 2
_
Flujo sobre una placa plana 4.29-4.40
4.5 Evalúe los grupos adimensionales, hcD>k, UqDr>m y
Analogía entre transferencia de calor 4.41-4.44
cpm>k para agua, alcohol n-butílico, mercurio, hidrógeno,
y de cantidad de movimiento
aire y vapor saturado a una _ temperatura de 100 °C. Sea
Disipación viscosa 4.45-4.48 D = 1 m, Uq = 1 m/sec y hc = 1 W/m2 K.
Problemas de diseño 4.49-4.59 4.6 Un fluido circula a 5 m/s sobre una placa plana ancha
Flujo a alta velocidad 4.60-4.66 de 15 cm de longitud. Para cada elemento de la lista
siguiente, calcule el número de Reynolds en el extremo
corriente abajo de la placa. Indique si el flujo en ese
4.1 Evalúe el número de Reynolds para flujo sobre un tubo punto es laminar, de transición o turbulento. Suponga que
a partir de los datos siguientes: D = 6 cm, Uq = 1.0 m/s, todos los fluidos están a 40 °C. a) aire, b) CO2, c) agua,
r = 300 kg/m3, m = 0.04 N s/m2. d) aceite para motores.
Problemas 285
4.10 El par de torsión debido a la resistencia de fricción de la

película de aceite entre un eje giratorio y su cojinete
15 cm depende de la fuerza F normal al eje, de la velocidad
de rotación N del eje, de la viscosidad dinámica m del
aceite y del diámetro del eje D. Establezca una correlación
∞ entre estas variables empleando el análisis dimensional.
4.11 Cuando una esfera cae libremente a través de un fluido
Fluido a homogéneo, alcanza una velocidad terminal a la que el
5 m/s peso de la esfera se equilibra por la fuerza de flotación y
Problema 4.6 por la resistencia de fricción del fluido. Haga un análisis
4.7 Trace de nuevo los puntos de datos de la figura 4.9b) dimensional de este problema e indique cómo se podrían
en papel logarítmico y determine una ecuación que correlacionar los datos experimentales para este problema.
aproxime la mejor línea de correlación. Compare sus Ignore los efectos de compresibilidad y la influencia de la
resultados con la figura 4.10. Luego suponga que rugosidad superficial.
vapor a 1 atm y a 100 °C fluye a través de un tubo con 4.12 Se han realizado experimentos sobre la distribu-
diámetro exterior de 5 cm a una velocidad de 1 m/s. ción de temperatura en un cilindro largo homogéneo
Utilizando los datos de la figura 4.10, estime el número (0.1 m de diámetro, conductividad térmica de 0.2 W/
de Nusselt, el coeficiente de transferencia de calor y la m K) con generación de calor interna uniforme. Mediante
tasa de transferencia de calor por metro de longitud de un análisis dimensional, determine la relación entre
tubo si éste está a 200 °C; compare estos resultados con la temperatura en régimen permanente en el centro
predicciones de su ecuación de correlación. del cilindro Tc, el diámetro, conductividad térmica y
4.8 El número de Reynolds promedio para aire que pasa tasa de generación de calor. Tome la temperatura en
con flujo turbulento sobre una placa plana de 2 m de la superficie como referencia. ¿Cuál es la ecuación para la
longitud es 2.4 * 106. En estas condiciones, el número temperatura del centro si la diferencia entre la tempera-
de Nusselt promedio se determinó igual a 4 150. Determine tura en el centro y en la superficie es 30 °C cuando la
el coeficiente de transferencia de calor promedio para un generación de calor es 3000 W/m3?
aceite que tiene propiedades térmicas similares a las que 4.13 Las ecuaciones de convección que relacionan los núme-
se encuentran en el apéndice 2, tabla 18, a 30 °C al mismo ros de Nusselt, Reynolds y Prandtl se pueden reaco-
número de Reynolds y que fluye__sobre la misma placa. modar para demostrar que para gases el coeficiente de
4.9 La relación adimensional Uq>¥ Lg , denominada número transferencia de calor hc depende
_____ de la temperatura
de Froude, es una medida de similitud entre un barco absoluta T y del grupo ¥
_____Uq >x . Esta formulación es
que navega por el océano y un modelo a escala del de la forma hc,x CT n¥ Uq >x , donde n y C son constan-
barco que se probará en un canal de agua en un labo- tes. Indique con claridad cómo se podría obtener esa
ratorio. Un barco de carga de 500 ft de longitud está relación para el caso de flujo laminar a partir de Nux =
0.332 Re0.5 0.333 para la condición 0.5 6 Pr 6 5.0. De ser
diseñado para que navegue a 20 nudos y un modelo de x Pr
5 ft geométricamente similar se remolca en un canal necesario establezca restricciones.
de agua para estudiar la resistencia de las olas. ¿Cuál debe 4.14 Mediante una serie de pruebas en las que se calentó agua
ser la velocidad de remolcado en ms-1? mientras fluía a través de un tubo de 38.6 in de longitud y
0.527 in de diámetro interior calentado eléctricamente se
20 nudos obtuvieron los datos de caída de presión experimentales
Prototipo que se muestran en la tabla siguiente.
Tempe- Tempe- Caída de

ratura ratura presión con la
Flujo en la masa superficial transferencia
500 ft másico m∙ del fluido T del tubo T de calor, ¢p
b s ht
U∞ = ¿? (lb/sec) (°F) (°F) (psi)
3.04 90 126 9.56
Modelo a escala 2.16 114 202 4.74
1.82 97 219 3.22
5 ft 3.06 99 248 8.34
2.15 107 283 4.45
Problema 4.9
Los datos de la caída de presión isotérmica para el 4.16 Los datos de pruebas que se muestran en la tabla
mismo tubo se dan a continuación en términos del siguiente se obtuvieron a partir de mediciones hechas
factor de fricción adimensional f = (¢p>rū 2)(2D>L)gc para determinar el coeficiente de transferencia de calor
y del número de Reynolds con base en el diámetro del dentro de tubos a números de Reynolds sólo ligera-
tubo, ReD = 4m∙>pDm. El símbolo ū denota la velocidad mente arriba de la transición y a números de Prandtl
promedio en el tubo. relativamente altos (como los asociados con aceites).
Las pruebas se realizaron en un intercambiador de tubo
ReD f dentro de tubo con un flujo transversal de agua para
fines de enfriamiento. El tubo utilizado para transportar
1.71 * 105 0.0189 los aceites era de 5/8 in de diámetro exterior, 18 BWG
1.05 * 105 0.0205 y 121 in de longitud. Correlacione los datos en términos
1.9 * 105 0.0185 de parámetros adimensionales apropiados.
2.41 * 105 0.0178
Prueba __
Comparando los coeficientes de fricción isotérmica con núm. Fluido hc ru cp kf mb mf
los no isotérmicos a números de Reynolds de la masa 11 Aceite 10C 87.0 1 072 000 0.471 0.0779 13.7 19.5
similares deduzca una ecuación adimensional para los 19 Aceite 10C 128.2 1 504 000 0.472 0.0779 13.3 19.1
coeficientes de fricción no isotérmicos en la forma 21 Aceite 10C 264.8 2 460 000 0.486 0.0776 9.60 14.0
f = constante * ReDn(ms>mb)m
23 Aceite 10C 143.8 1 071 000 0.495 0.0773 7.42 9.95
24 Aceite 10C 166.5 2 950 000 0.453 0.0784 23.9 27.3
25 Aceite 10C 136.3 1 037 000 0.496 0.0773 7.27 11.7
donde m s = viscosidad a temperatura superficial,
36 Piranol 140.7 1 795 000 0.260 0.0736 12.1 16.9
mb = viscosidad a temperatura de masa y n y m = cons-
1488
tantes empíricas.
39 Piranol 133.8 2 840 000 0.260 0.0740 23.0 29.2
Ts 1488
45 Piranol 181.4 1 985 000 0.260 0.0735 10.3 12.9
1488
Fluido 48 Piranol 126.4 3 835 000 0.260 0.0743 40.2 53.5
0.527 in 1488
m, Tb
49 Piranol 105.8 3 235 000 0.260 0.0743 39.7 45.7
1488
38.6 in
_
Problema 4.14 donde hc = coeficiente de transferencia de calor
superficial promedio basado en la dife-
rencia de temperatura media, Btu/h ft2°F
4.15 Los datos experimentales que se muestran en la tabla ru = velocidad de la masa, lbm/h ft2
siguientes se obtuvieron haciendo pasar alcohol n-bu- cp = calor específico, Btu/lbm °F
tílico a una temperatura de masa de 15 °C sobre una kf = conductividad térmica, Btu/h ft °F
placa plana calentada de 0.3 m de longitud, 0.9 m de (basada en la temperatura superficial)
ancho y con una temperatura superficial de 60 °C. mb = viscosidad, basada en la temperatura de
Correlacione los datos experimentales utilizando núme- masa promedio (media mezclada), lbm/h ft
ros adimensionales apropiados y compare la línea de mf = viscosidad, basada en la temperatura
mejor ajuste de los datos con la ecuación (4.38). pelicular promedio, lbm/h ft
Sugerencia: Inicie correlacionando Nu y ReD
Coeficiente de independientemente de los números de Prantl ya que
transferencia la influencia del número de Prandtl en el número de
Velocidad de calor promedio Nusselt se espera que sea relativamente pequeña. Al
(m/s) (W/m2°C) trazar Nu contra Re en papel logarítmico, se puede
suponer la naturaleza de la ecuación de correlación,
0.089 121
Nu = f1(Re). Después una gráfica de Nu/f1(Re) contra
0.305 218
Pr revelará la dependencia en Pr. Para la ecuación
0.488 282
final, también se considerará la influencia de la varia-
1.14 425
ción de la viscosidad.
Problemas 287
4.17 Un álabe de una turbina con una longitud característica de Escriba una correlación adimensional apropiada para
1 m se enfría en un túnel de viento a presión atmosférica el número de Nusselt promedio para estos datos y
por aire a 40 °C con una velocidad de 100 m/s. Para una establezca cualesquiera limitaciones de su ecuación.
temperatura superficial de 500 K, la tasa de enfriamiento 4.20 Sobre una superficie plana y de manera paralela a ella
es de 10 000 watts. Utilice estos resultados para estimar fluye aceite para motores a 100 °C a una velocidad de
la tasa de enfriamiento de otro álabe de turbina de forma 3 m/s. Calcule el espesor de la capa límite hidrodiná-
similar, pero con una longitud característica de 0.5 m y mica a una distancia de 0.3 m del borde de ataque de
que funciona con una temperatura superficial de 600 K la superficie.
en aire a 40 °C con una velocidad de 200 m/s.
Aceite
3 m/s
100 °C δ =?
0.3 m
Problema 4.20
1.0 m 0.5 m
4.21 Suponiendo una distribución lineal de velocidad y
distribución lineal de temperatura en la capa límite
Aire 600 K
Aire sobre una placa plana, deduzca una relación entre los
200 m/s,
100 m/s, 40 °C 500 K 40 °C espesores de las capas límites, térmica e hidrodinámica
y el número de Prandtl.
Problema 4.17 4.22 Entre dos placas planas paralelas separadas 5 cm
4.18 El rozamiento sobre la ala de un avión en vuelo se sabe fluye aire a 20 °C a 1 m/s. Estime la distancia desde la
que es una función de la densidad del aire (r), de la entrada hasta el punto en el que convergen las capas
viscosidad del aire (m), de la velocidad de corriente libre límites hidrodinámicas.
(Uq), de una dimensión característica del ala (s) y del
esfuerzo cortante sobre su superficie (ts).
Demuestre que el rozamiento adimensional, Aire
ts
____ 1 m/s 5 cm
rUq2
, se puede expresar como una función del número 20 °C
rUqs
de Reynolds, ____ m .
?
4.19 Suponga que la gráfica siguiente muestra valores
medidos de hc para aire en convección forzada sobre Problema 4.22
un cilindro de diámetro D, trazados en una gráfica
logarítmica de NuD como una función de ReDPr. 4.23 Un fluido a temperatura Tq circula a una velocidad Uq
sobre una placa plana que está a la misma temperatura
que el fluido a lo largo de una distancia x0 desde el
103 borde de ataque, pero a una temperatura mayor Ts más
allá de este punto.
102 Fluido
NuD
U∞ δ Ts
δt
T∞
101 Sin calentar Calentada

x0
1 101 102 103 Problema 4.23

(ReD)(Pr)
Demuestre mediante las ecuaciones integrales de la capa
Problema 4.19 límite que z, la relación del espesor de la capa límite
térmica al espesor de la capa límite hidrodinámica, sobre Siguiendo el enfoque descrito en la sección 4.9.1 para
la parte calentada de la placa es aproximadamente de condiciones laminares, sustituya las relaciones en las
ecuaciones integrales de la capa límite, cantidad de
x0 3/4 1/3
z L Pr-1/3 c 1 - a b d
movimiento, energía y deduzca ecuaciones para: a) el
x espesor de la capa límite, b) el coeficiente de fricción
local y c) el número de Nusselt local. Suponga que
si el flujo es laminar.
@ = @t y explique las limitaciones de sus resultados.
4.24 Entre dos placas paralelas separadas 1 cm fluye aire a
4.28 Para metales líquidos con números de Prandtl mucho
1 000 °C a una velocidad de entrada de 2 m/s. Estime
menores que la unidad, la capa límite hidrodinámica
la distancia desde la entrada hasta el punto donde con-
es mucho más delgada que la capa límite térmica. Como
vergen la capas límites.
resultado, se podría suponer que la velocidad en la capa
4.25 Mediciones experimentales de la distribución de
límite es uniforme [u = Uq y v = 0]. Partiendo con la
temperatura del flujo de aire a presión atmosférica
ecuación (4.7b), demuestre que la ecuación de energía y su
sobre el ala de un avión indican que la distribución
condición límite son análogas a las correspondientes para
de temperatura cerca de la superficie se puede aproxi-
una losa semiinfinita con un cambio repentino en la tem-
mar mediante una ecuación lineal:
peratura superficial [consulte la ecuación (2.93)]. Después
(T - Ts) = ay(Tq - Ts) demuestre que el número de Nusselt local está dado por
donde a = una constante = 2 m-1, Ts = temperatura Nux = 0.56 (RexPr)0.5
superficial, K, Tq = temperatura de corriente libre, K,
Compare esta ecuación con la relación apropiada en
y y = distancia perpendicular desde la superficie (mm).
la tabla 4.5.
a) Estime el coeficiente de transferencia de calor por
4.29 A lo largo de una placa plana fluye hidrógeno a 15 °C y
convección si Ts = 50 °C y Tq = -50 °C. b) Calcule el
a una presión de 1 atm, con una velocidad de 3 m/s. Si la
flujo de calor en W/m2.
placa tiene un ancho de 0.3 m y está a 71 °C, calcule las can-
4.26 Para flujo sobre una superficie isotérmica ligeramente
tidades siguientes en x = 0.3 m y a la distancia correspon-
curva, la distribución de temperatura dentro de la capa
diente al punto de transición, es decir, Rex = 5 * 105 (utilice
límite @t se puede aproximar mediante el polinomio
propiedades a 43 °C): a) el espesor de la capa límite hidro-
T(y) = a + by + cy2 + dy3(y 6 @t), donde y es la distan-
dinámica, en cm, b) el espesor de la capa límite térmica, en
cia normal a la superficie. a) Aplicando condiciones
cm, c) el coeficiente de fricción local, adimensional, d) el
límites apropiada evalúe las constantes a, b, c y d.
coeficiente de fricción promedio, adimensional, e) la fuerza
de rozamiento, en N, f) el coeficiente de transferencia de
Fluido δt y calor por convección, en W/m2°C, g) el coeficiente de trans-
Ts ferencia de calor por convección promedio, en W/m2 °C
y h) la tasa de transferencia de calor, en W.
4.30 Repita el problema 4.29, incisos d), e), g) y h) con
Problema 4.26 x = 4.0 m y Uq = 80 m/s, a) tomando en cuenta la capa
límite laminar y b) suponiendo que la capa límite turbu-
b) Después obtenga una relación adimensional para la
lenta empieza en el borde de ataque.
distribución de temperatura en la capa límite.
4.31 Determine la tasa de pérdida de calor en Btu/h de la pared
4.27 El método integral también se puede aplicar a condi-
de un edificio que resulta de un viento a 10 mph que sopla
ciones de flujo turbulento si se dispone de datos expe-
paralelo a su superficie. La pared es de 80 ft de longitud
rimentales del esfuerzo cortante en la pared. En uno de
y 20 ft de altura, su temperatura superficial es 80 °F y la
los primeros intentos para analizar el flujo turbulento
temperatura del aire ambiente es 40 °F.
sobre una placa plana, Ludwig Prandtl propuso en
4.32 Un intercambiador de calor de placa plana funcionará en
1921 las relaciones siguientes para la velocidad adi-
una atmósfera de nitrógeno a una presión de aproxima-
mensional y distribuciones de temperatura:
damente 104 N/m2 y 38 °C. El intercambiador de calor
u> Uq = (y> d)1/7 de placa plana originalmente se diseñó para funcionar
(T - Tq)> (Ts - Tq) = 1 - (y> dt)1/7 (Ts 7 T 7 Tq)

en aire a 1 atm y 38 °C en flujo turbulento. Estime la
relación del coeficiente de transferencia de calor en aire
al correspondiente en nitrógeno, suponiendo enfriamiento
A partir de datos experimentales, una relación empírica por circulación forzada de la superficie de la placa plana
que relaciona el esfuerzo cortante en la pared con el a la misma velocidad en los dos casos.
espesor de la capa límite es: 4.33 Un intercambiador de calor para calentar mercurio líquido
ts = 0.023rU2q > Re0.25
d donde Red = Uqd> n está en desarrollo. El intercambiador se puede visualizar
Problemas 289
como una placa plana de 6 in de longitud y 1 ft de ancho. La 0.1, 1 y 5 m desde el borde de ataque. Explique el efecto del
placa se mantiene a 160 °F y el mercurio fluye paralelo al gradiente de temperatura a lo largo de la cuerda.
lado corto a 60 °F con una velocidad de 1 ft/s. a) Determine
el coeficiente de fricción local en el punto medio de la placa
y la fuerza de rozamiento total sobre la placa. b) Determi-
ne la temperatura del mercurio en un punto a 4 in del borde
de ataque y a 0.05 in de la superficie de la placa. c) Calcu-
le el número de Nusselt en el extremo de la placa.
Mercurio líquido 160 °F

1 ft/s
60 °F
1 ft
6 in Radiación solar
100 W/m2
Problema 4.33
4.34 Agua a una velocidad de 2.5 m/s fluye paralela a una Aire
Aislada
placa plana y delgada de 1 m de longitud, horizontal y lisa. 150 m/s
Determine los espesores de la capa límite térmica local e 6m
hidrodinámica y el coeficiente de fricción local en el punto
medio de la placa. ¿Cuál es la tasa de transferencia de calor Problema 4.37
de un lado de la placa al agua por ancho unitario de la placa 4.38 Una aleta de enfriamiento de aluminio para un intercam-
si la temperatura superficial se mantiene uniforme a 150 °C biador de calor está ubicada paralela a una corriente de
y la temperatura de la corriente de agua principal es 15 °C? aire a presión atmosférica. La aleta, como se muestra en
4.35 Una placa plana delgada se coloca en una corriente de aire el bosquejo, tiene una altura de 0.075 m, un espesor de
a presión atmosférica que fluye paralela a ella a una veloci- 0.005 m y 0.45 m en la dirección del flujo. Su temperatura
dad de 5 m/s. La temperatura en la superficie de la placa se en la base es 88 °C y el aire está a 10 °C. La velocidad del
mantiene uniformemente a 200 °C y la correspondiente a la aire es 27 m/s. Determine la fuerza de rozamiento total y
corriente de aire principal es 30 °C. Calcule la temperatura la tasa total de transferencia de calor de la aleta al aire.
y la velocidad horizontal en un punto a 30 cm del borde
de ataque y a 4 mm arriba de la superficie de la placa.
4.36 La temperatura superficial de una placa plana delgada ubi-
cada paralela a una corriente de aire es de 90 °C. La velo- 0.005 m
cidad de corriente libre es 60 m/s y la temperatura del aire
es 0 °C. La placa tiene un ancho de 60 cm y una longitud
de 45 cm en la dirección de la corriente de aire. Ignorando 0.075 m 0.45 m
el efecto del extremo de la placa y suponiendo que el
flujo en la capa límite cambia abruptamente de laminar
a turbulento a un número de Reynolds de Retr = 4 * 105, Aire 27 m/s
determine: a) el coeficiente de transferencia de calor 10 °C
promedio en las regiones laminar y turbulenta, b) la tasa Problema 4.38
de transferencia de calor para toda la placa, considerando
los dos lados, c) el coeficiente de fricción promedio en las 4.39 Se utiliza aire a 320 K con una velocidad de corriente libre
regiones laminar y turbulenta y d) la fuerza de rozamiento de 10 m/s para enfriar dispositivos electrónicos montados
total. Además, trace el coeficiente de transferencia de sobre una tarjeta de circuitos impresos, como se mues-
calor y el de fricción local como una función de la distan- tra en el bosquejo siguiente. Cada dispositivo tiene una
cia desde el borde de ataque de la placa. plataforma cuadrada de 5 * 5 mm y disipa 60 mW. Un
4.37 El ala de un avión tiene una cubierta de aluminio pulido. dispositivo para crear una corriente turbulenta está ubicado
A una altitud de 1 500 m, absorbe 100 W/m2 por radiación. en el borde de ataque para ocasionar que la capa límite se
Suponiendo que la superficie interior de la cubierta del ala vuelva turbulenta. Suponiendo que las superficies inferio-
está bien aislada y que tiene una cuerda de 6 m de longitud res de los dispositivos electrónicos están aisladas, estime la
(es decir, L = 6 m), estime la temperatura de equilibrio del temperatura superficial en el centro del quinto dispositivo
ala a una velocidad de vuelo de 150 m/s a distancias de en la tarjeta de circuitos (consulte la página siguiente).
Aire Dispositivos 50 m/s, estime el flujo de calor a 1 m del borde de ataque

10 m/s 0.005 m electrónicos para una temperatura superficial de la placa de 300 K.
320 K 4.45 Cuando la disipación viscosa es apreciable, en la
conservación de la energía [ecuación (4.6)] se debe
tomar en cuenta la tasa a la que la energía mecánica se
Dispositivo para crear una corriente turbulenta convierte irreversiblemente en energía térmica debido
Problema 4.39 a los efectos viscosos en el fluido. Esto da origen a un
término adicional, f, la disipación viscosa, en el lado
4.40 El coeficiente de fricción promedio para flujo sobre una derecho, donde:
placa larga de 0.6 m de longitud es 0.01. ¿Cuál es el valor
0u 0v 2 0u 2 0v 2
= a b + 2c a b + a b d
de la fuerza de rozamiento para los siguientes fluidos? f
+
a) Aire a 15 °C, b) vapor a 100 °C y presión atmosfé- m 0y 0x 0x 0y
rica, c) agua a 40 °C, d) mercurio a 100 °C y e) alcohol
2 0u 0v 2
n-butílico a 100 °C. - a + b
4.41 Una placa plana, delgada y cuadrada de 6 in por lado 3 0x 0y
se someterá a una prueba de rozamiento en un túnel de Aplique la ecuación resultante a un flujo laminar
viento con aire a 100 ft/s, 14.7 psia y 60 °F que fluye en entre dos placas paralelas infinitas, en donde la placa
sentido transversal y paralelo a las superficies superior superior se mueve a una velocidad U. Suponiendo
e inferior. La fuerza de rozamiento total observada es propiedades físicas constantes (r, cp, k y m), obtenga
0.0135 lb. Utilizando la definición del coeficiente de fric- expresiones para las distribuciones de velocidad y
ción [ecuación (4.13)] y la analogía de Reynolds, calcule temperatura. Compare la solución con el término
la tasa de transferencia de calor de esta placa cuando la de disipación incluido con los resultados cuando se
temperatura superficial se mantiene a 250 °F. ignora la disipación. Encuentre la velocidad de la
4.42 Una placa plana, delgada y cuadrada de 15 cm por placa necesaria para producir un aumento de tempera-
lado está suspendida de una balanza en una corriente tura de 1 K en aire nominalmente a 40 °C relativa al
de aceite para motores que fluye uniformemente de tal caso donde se ignora la disipación.
manera que el aceite fluye paralelo a y a lo largo de las 4.46 Una chumacera se puede idealizar como una placa plana
dos superficies de la placa. El rozamiento total sobre la estacionaria y una placa plana móvil que se mueve paralela
placa se midió igual a 55.5 N. Si el aceite fluye a una a ella. El espacio entre las dos placas está lleno con un
velocidad de 15 m/s con una temperatura de 45 °C, fluido incompresible. Considere un cojinete en el que las
calcule el coeficiente de transferencia de calor utili- placas estacionaria y móvil están a 10 y 20 °C, respecti-
zando la analogía de Reynolds. vamente, la distancia entre ellas es 3 mm, la velocidad de
4.43 Para un estudio sobre el calentamiento global, se ha la placa móvil es 5 m/s y hay aceite para motores entre las
diseñado un instrumento para mapear las características placas. a) Calcule el flujo de calor hacia las placas superior
de absorción de CO2 del Océano Pacífico. El paquete de e inferior. b) Determine la temperatura máxima del aceite.
instrumentos se parece a una placa plana con un área
superficial total (superior e inferior) de 2 m2. Para ope-
Eje o muñón
rarlo con seguridad, su temperatura superficial no debe
Cojinete
exceder la temperatura del océano en más de 2 °C. Para
monitorear la temperatura del paquete de instrumentos,
que se remolca por un barco que se navega a 20 m/s, se
mide la tensión en el cable de remolque. Si la tensión
es 400 N, calcule la tasa de generación de calor máxima
permisible del paquete de instrumentos.
4.44 Para flujo de gas sobre una superficie plana que se
asperezó de manera artificial con un chorro de arena, el 20 °C
coeficiente de transferencia de calor local por convección
se puede correlacionar mediante la relación adimen-
sional Nux= 0.05Re0.9x . a) Deduzca una relación para el 3 mm Fluido
coeficiente de transferencia de calor promedio en flujo
sobre una placa plana de longitud L. b) Suponiendo que
es válida la analogía entre la transferencia de calor y
de cantidad de movimiento, deduzca una relación para
10 °C
el coeficiente de fricción local. c) Si el gas es aire a
una temperatura de 400 K fluyendo a una velocidad de Problema 4.46
Problemas 291
4.47 Una chumacera tiene una holgura de 0.5 mm. El muñón 20 °C sopla sobre el techo a una velocidad de 15 mph.
tiene un diámetro de 100 mm y gira a 3 600 rpm dentro Estime la pérdida de calor por convección del colector
del cojinete. Está lubricado con un aceite que tiene al aire cuando el colector está montado: a) en el bor-
una densidad de 800 kg/m3, una viscosidad de 0.01 kg/ms de de ataque del techo (Lc = 0) y b) a una distancia de
y una conductividad térmica de 0.14 W/m K. Si la super- 10 m del borde de ataque.
ficie del cojinete está a 60 °C, determine la distribución
de temperatura en la película de aceite, suponiendo que Aire
la superficie de la chumacera está aislada. 15 mph Colector solar
20 °C
4.48 Una chumacera tiene una holgura de 0.5 mm. El muñón Techo plano
tiene un diámetro de 100 mm y gira a 3 600 rpm dentro del
cojinete. Está lubricado con un aceite que tiene una den-
sidad de 800 kg/m3, una viscosidad de 0.01 kg/ms y una
conductividad térmica de 0.14 W/m K. Si las temperaturas
del muñón y del cojinete se mantienen a 60 °C, calcule
la tasa de transferencia de calor del cojinete y la potencia
necesaria por rotación por unidad de longitud.
4.49 Un camión refrigerado viaja a 80 mph por una carretera
en el desierto donde la temperatura del aire es 50 °C.
El cuerpo del camión se puede idealizar como una caja
rectangular de 3 m de ancho, 2.1 m de altura y 6 m de
longitud, a una temperatura superficial de 10 °C. Suponga
Problema 4.50
que: 1) la transferencia de calor de la parte frontal y poste- 4.51 Un dispositivo electrónico se enfriará por aire flu-
rior del camión se puede ignorar, 2) la corriente de aire no yendo sobre aletas de aluminio fijas a la superficie
se separa de la superficie y 3) la capa límite es turbulenta inferior como se muestra:
sobre toda la superficie. Si se requiere una capacidad de
1 ton por cada 3600 W de pérdida de calor, calcule el
U∞ = 10 m/s
tonelaje necesario de la unidad de refrigeración. T∞(aire) = 20 °C
3m Dispositivo
Aislamiento
electrónico
Aletas de
2.1 m aluminio
3 mm
Aire 6m 6 mm
80 mph
50 °C 10 mm
1 mm
3 mm
U∞
9 mm
Problema 4.51
El dispositivo disipa 5 W y la resistencia por contacto
térmico entre la superficie inferior del dispositivo y la
superficie superior del conjunto de aletas de enfria-
Problema 4.49 miento es 0.1 cm2 K/W. Si el dispositivo está a una
4.50 En el ecuador, donde el sol está aproximadamente arriba temperatura uniforme y aislado en la parte superior,
de la cabeza al medio día, una orientación casi óptima estime la temperatura en régimen permanente.
para un calentador solar de placa plana es en la posi- 4.52 Un arreglo de 16 láminas de silicio configuradas en
ción horizontal. Suponga que para calentar agua para dos filas está aislado en la parte inferior y se enfría por
uso doméstico se utiliza un colector solar de 4 * 4 m aire que fluye por convección forzada sobre la parte
montado sobre un techo horizontal, como se muestra superior. El arreglo se puede colocar con su lado largo
en el bosquejo siguiente. La temperatura superficial de o bien con el corto de frente al aire de enfriamiento. Si
la cubierta de cristal se estima que es 40 °C y aire a cada chip tiene un área superficial de 10 * 10 mm y
disipa la misma energía, calcule la tasa de disipación de 4.54 Las seis aletas idénticas de aluminio que se muestran en
energía máxima permisible para los dos arreglos si la tem- el siguiente bosquejo están fijas a un dispositivo elec-
peratura superficial máxima permisible de las láminas es trónico para su enfriamiento. El aire de enfriamiento
100 °C. ¿Cuál sería el efecto de un dispositivo para crear está disponible a una velocidad de 5 m/s de un ventila-
una corriente turbulenta en el borde de ataque para ocasio- dor a 20 °C. Si la temperatura promedio en la base de
nar la capa límite en flujo turbulento? La temperatura del una aleta no debe exceder 100 °C, estime la disipación
aire es 30 °C y su velocidad es 25 m/s. de energía máxima permisible del dispositivo.
Láminas de silicio
Aire, caso b 1 mm
6 cm
6 mm
Aislado 1 cm
Aire
Aire, caso a Problema 4.54
25 m/s
30 °C
4.55 Veinticinco láminas cuadradas de computadora, cada
Problema 4.52 una con tamaño de 10 * 10 mm y 1 mm de espesor,
4.53 Se quiere dimensionar el sistema de acondicionamien- están montadas separadas 1 mm en un sustrato plástico
to de aire en una furgoneta Chevrolet para su uso en aislante como se muestra en el siguiente bosquejo. Las
climas desérticos. El sistema debe mantener una tem- láminas se enfriarán con nitrógeno fluyendo a lo largo
peratura interior de 20 °C cuando la furgoneta viaje a de la longitud de la fila a -40 °C y presión atmosférica
100 km/h a través de aire seco a 30 °C en la noche. Si la para evitar que su temperatura exceda 30 °C. El diseño
parte superior de la furgoneta se puede idealizar como debe proporcionar una tasa de disipación de 30 mW
una placa plana de 6 m de longitud y 2 m de ancho por chip. Estime la velocidad de corriente libre mínima
y los lados como placas planas de 3 m de altura y 6 m de necesaria para proporcionar condiciones de seguridad
longitud, estime la tasa a la que el calor se debe remover para cada chip en el arreglo.
del interior para mantener las condiciones de comodidad
especificadas. Suponga que el coeficiente de transferencia Nitrógeno Ts = 30 °C
de calor en el interior de la furgoneta es 10 W/m2 K. T∞ = −40 °C 10 mm
1 mm
Plástico
2m L = 25 (10 mm) + 24 (1 mm) = 274 mm
Problema 4.55
3m
4.56 Se propone remolcar témpanos de hielo de la región
polar al Medio Este a fin de suministrar agua potable a
regiones áridas allí. Un témpano común adecuado para
6m remolcarlo deberá ser relativamente ancho y plano.
Aire
Considere un témpano de 0.25 km de espesor y 1 km
100 km/h
30 °C
por lado. El témpano se remolcará a 1 km/h una dis-
Furgoneta idealizada tancia de 6 000 km a través de agua cuya temperatura
promedio durante el viaje es 8 °C. Suponiendo que la
interacción del témpano con sus alrededores se puede
aproximar por la transferencia de calor y fricción en su
superficie inferior, calcule los parámetros siguientes:
a) la tasa a la que el hielo se fundirá en la superficie
inferior y b) la potencia requerida para remolcar el
témpano a la velocidad designada. c) Si la energía de
Problema 4.53 remolcar cuesta aproximadamente 50 ¢/kW ∙ h por
Problemas 293
hora de potencia y el costo de suministrar agua en el qs /A

destino se puede aproximar con la misma cifra, calcule
el costo de agua dulce. El calor latente de fusión del Flujo de aire
Cubierta de cristal
hielo es 334 kJ/kg y su densidad es 900 kg/m3.
4.57 En un proceso de manufactura, una tira larga de plan-
cha metálica se transporta en una banda de transporte a Espacio de aire Placa absorbente
una velocidad de 2 m/s mientras un recubrimiento en
su parte superior tiene que curarse por calentamiento Conductos
radiante. Suponga que lámparas infrarrojas montadas para fluido
arriba de la banda proporcionan un flujo radiante de
2500 W/m2 sobre el recubrimiento. El recubrimiento Aislamiento
absorbe 50% del flujo radiante incidente, tiene una emi- Problema 4.58
sividad de 0.5 e irradia hacia los alrededores a una
temperatura de 25 °C. Además, el recubrimiento pierde
calor por convección a través de un coeficiente de trans- c = 0.20 Btu/lb °F) en el extremo de una barra de
ferencia de calor por convección tanto en la superficie madera de 1 in de diámetro a una temperatura uni-
superior como en la inferior y el aire ambiente, que se forme de 212 °F repentinamente se coloca en una
puede suponer que está a la misma temperatura que el corriente de vapor a 60 °F y 100 ft/s que fluye para-
ambiente. Estime la temperatura del recubrimiento en lela al eje de la barra. Estime la temperatura prome-
condiciones de régimen permanente. dio en el eje de la barra de tránsito 8 min después
de que inicia el enfriamiento. Suponga conducción de
calor radial, pero incluya las pérdidas por radiación
basadas en una emisividad de 0.90 a alrededores
negros a la temperatura del aire.
4.60 Una placa plana de cromo altamente pulido se
Temperatura coloca en un túnel de viento de alta velocidad para
superficial = ¿?
simular el flujo sobre el fuselaje de una aeronave
Flujo radiante = 2500 W/m2
supersónica. El aire que fluye en el túnel de viento
2 m/s
está a una temperatura de 0 °C y a una presión de
3500 N/m2 y tiene una velocidad paralela a la placa
Problema 4.57 de 800 m/s. ¿Cuál es la temperatura en la pared
adiabática en la región laminar y cuál es la longitud
4.58 El siguiente esquema es de un colector solar de placa de la capa laminar?
plana de 4 m2 para calentamiento de agua para uso 4.61 Fluye aire, a una temperatura estática de 70 °F y a una
doméstico. Sobre la cubierta de cristal incide una radia- presión estática de 0.1 psia, a un ángulo de ataque cero
ción solar a una tasa de 750 W/m2, que transmite 90% sobre una placa plana calentada eléctricamente a una
del flujo incidente. Fluye agua a través de los tubos velocidad de 20 fps. Si la placa tiene una longitud de
soldados en la parte posterior de la placa absorbente, 4 in en la dirección del flujo y de 24 in en la dirección
entrando con una temperatura de 25 °C. La cubierta normal al flujo, determine la tasa de disipación de
de cristal tiene una temperatura de 27 °C en régimen calor eléctrico necesaria para mantener la placa a una
permanente e irradia calor con una emisividad de 0.92 temperatura promedio de 130 °F.
al cielo a -50 °C. Además, la cubierta de cristal pierde 4.62 El rechazo de calor de automóviles de carreras es un
calor por convección al aire a 20 °C que fluye sobre su problema debido a que los intercambiadores de calor
superficie a 20 mph. requeridos por lo general crean un rozamiento adicio-
a) Calcule la tasa a la que el calor se colecta por nal. Se ha propuesto probar la integración del rechazo
el fluido de trabajo, es decir, el agua en los tubos, de calor hacia la cubierta del vehículo para un automó-
por área unitaria de la placa absorbente. b) Calcule la vil en el lago salado Bonneville Salt Flats. Las pruebas
eficiencia del colector hc, que se define como la rela- preliminares se realizarán en un túnel de viento en una
ción de la energía útil transferida al agua en los tubos placa plana sin rechazo de calor. El aire atmosférico
a la energía solar incidente en la placa de cubierta del en el túnel está a 10 °C y fluye a 250 m/s-1 sobre la
colector. c) Calcule la temperatura de salida del agua placa plana térmicamente no conductora de 3 m de
si su tasa de flujo a través del colector es 0.02 kg/s. longitud. ¿Cuál es la temperatura de la placa a 1 m
El calor específico del agua es 4 179 J/kg K. corriente abajo desde el borde de ataque? ¿Qué tanto
4.59 Una barra de tránsito de 6 in de longitud y 1 in de difiere esta temperatura de la temperatura a 0.005 m
diámetro (k = 0.56, Btu/h ft °F, r = 100 lb/cu ft, del borde de ataque?
4.63 Fluye aire sobre una tira plana metálica de 0.1 m de

longitud a 15 °C y a una presión atmosférica de 0.01
a una velocidad de 250 m/s en la dirección del flujo.
Determine: a) la temperatura superficial de la placa en
equilibrio y b) la tasa de remoción de calor necesaria
por metro de ancho si la temperatura superficial se
debe mantener a 30 °C.
4.64 Una placa plana se coloca en un túnel de viento supersó-
nico con aire que fluye sobre ella a un número de Mach
de 2.0, a una presión de 25 000 N/m2 y a una tempera-
tura ambiente de -15 °C. Si la placa tiene una longitud
de 30 cm en la dirección del flujo, calcule la tasa de
Túnel de viento
enfriamiento por área unitaria que se necesita para mante-
ner la temperatura de la placa por debajo de 120 °C.
4.65 Un satélite reingresa en la atmósfera de la Tierra a una
velocidad de 2 700 m/s. Estime la temperatura máxima
que alcanzara el blindaje de calor si no se permite que
Placa su material se desprenda y se ignoran los efectos de
la radiación. La temperatura de la capa superior de la
atmósfera es '= -50 °C.
Aire 4.66 En un túnel de viento se prueba un modelo a escala de
250 m/s una sección del ala de un avión a un número de Mach
10 °C 1.5. La presión y temperatura en la sección de prueba son
Problema 4.62 20 000 N/m2 y -30 °C, respectivamente. Si la sección
del ala debe mantenerse a una temperatura promedio de
60 °C, determine la tasa de enfriamiento requerida. El
modelo del ala se puede aproximar como una placa plana
de 0.3 m de longitud en la dirección del flujo.
Problemas de diseño
4.1 Carga de calefacción en una fábrica (continuación 4.2 Camión refrigerado (continuación del problema de
del problema de diseño 1.3) diseño 2.10)
En el problema de diseño 1.3, usted calculó la pérdida Refine el diseño térmico del contenedor de embarque
de calor de un edificio industrial pequeño en el invierno. refrigerado del problema de diseño 2.10 calculando el
En los cálculos iniciales, estimó el coeficiente de trans- coeficiente de transferencia de calor por convección sobre
ferencia de calor por convección de la tabla 1.3. Ahora los dos lados y la parte superior a 65 mph a partir de la
repita los cálculos de la pérdida de calor, pero calcule el información presentada en este capítulo. La pérdida de calor
coeficiente de transferencia de calor externo con base de las partes frontal y posterior del camión se puede estimar
en el material presentado en este capítulo. Para estimar con información de la tabla 1.3. La pérdida de calor de la
las condiciones de convección forzada en el edificio, parte inferior del contenedor será relativamente pequeña
suponga que los vientos en Denver, Colorado, pueden debido a que queda relativamente cerca de la carretera sobre
ser tan altos como 70 mph, pero en condiciones nor- la plataforma del camión.
males no sobrepasarán 20 mph. Analice el efecto de la 4.3 Alberca con calentamiento solar
iluminación en la carga de calor y estime cuál sería el Estime la pérdida de calor de la superficie de una
efecto de las luces eléctricas si se necesitan 20 bombillos alberca con área superficial de 3 * 10 m y una profun-
eléctricos de 150 W para proporcionar una iluminación didad promedio de 1.5 m y después diseñe un sistema
adecuada para este edificio industrial. ¿Puede sugerir un de calentamiento solar para ella. Con anterioridad la
método de iluminación mejorado? alberca se calentó con un calentador eléctrico que se
Problemas de diseño 295
Reacondicionamiento del sistema de calentamiento puede utilizar como respaldo cuando no haya sufi-
solar de la alberca ciente radiación solar para mantener la temperatura
Paneles del colector de la alberca, pero que pueda apagarse cuando haya
suficiente radiación para que los colectores calienten
la alberca. Los paneles del colector están hacia el sur
con una pendiente igual a la latitud menos 10 °C, como
Alberca recomiendan los expertos en calentadores solares (con-
sulte, por ejemplo, J.F. Kreider, C.J. Hoogendoorn y F.
Kreith, Solar Design: Component, Systems, Economics,
Calentador Hemisphere Publishing, 1989). La alberca está ubicada
en San Diego, California, donde la velocidad promedio
Filtro del viento es 10 mph, se debe mantener a una tempera-
Válvula
tura de 28 °C todo el año. El sistema de calentamiento
Bomba
solar es un bucle cerrado de agua con una configura-
ción de tuberías que permite que el flujo sólo llegue
Detalle del panel solar al calentador existente, el flujo a través del calentador
Radiación eléctrico, de los paneles solares y sólo hacia los paneles
solar incidente solares, como se muestra en el bosquejo. El colector
solar tiene que ser de una extrusión de plástico negro
sin una cubierta para minimizar su costo. El coeficiente
de transferencia de calor para el agua en los pasajes de
flujo rectangular del colector se puede estimar con base
en la tabla 1.3. Sugiera maneras para reducir la pérdida
de calor de la alberca de noche y estime la efectividad
Agua
en cuanto al costo del sistema.
Problema de diseño 4.3

Frank Kreith C-4 PDF

Cargado por

Copyright:

Formatos disponibles

Frank Kreith C-4 PDF

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Frank Kreith C-4 PDF

Cargado por

Copyright:

Formatos disponibles

CAPÍTULO 4

4.2 Transferencia de calor por convección

FIGURA 4.1 Distribuciones de la velocidad y la

SOLUCIÓN De la ecuación (4.1), el coeficiente de transferencia de calor se puede expresar en la

-10.028 W/m K21- 66.7 K/mm2

FIGURA 4.2 Datos experimentales sobre la distribución

FIGURA 4.3 Distribuciones de

La situación es muy similar en convección natural, como se muestra en la

4.3 Fundamentos de la capa límite

Región laminar Transición Turbulenta

FIGURA 4.4 Perfiles de velocidad en capas de límite laminar,

donde du/dy es el gradiente de velocidad y la constante de proporcionalidad m se

fuerzas viscosas e inerciales y cuyo valor determina la transición de flujo laminar a

donde Uq = velocidad de corriente libre

El régimen de transición se extiende hasta un número de Reynolds de aproximada-

4.4 Ecuaciones de conservación de masa, cantidad de movimiento

fluido/pared para evaluar el coeficiente de transferencia de calor. Un enfoque un

Flujos de cantidad de movimiento

(Â u)u dy Â (u + ∂u dx) . ( u + ∂u dx)

La fuerza cortante viscosa en la cara inferior es t ƒ y = -m (0u/0y) dx y sobre la cara

FIGURA 4.7 Volumen de control diferencial para la con-

La ecuación de la conservación de energía se dedujo con base en la suposición de

para propiedades constantes. El efecto de la disipación viscosa puede ser signi-

4.5 Ecuaciones adimensionales de la capa límite y parámetros de similitud

y temperatura. En general, estas soluciones son muy complicadas y el lector debe

donde L es una dimensión de longitud característica como la longitud de una placa,

δ = Capa límite de velocidad

FIGURA 4.8 Desarrollo de las capas límite de velocidad y térmica

geométricamente similar, siempre que los parámetros de similitud tengan el mismo

4.5.1 Coeficiente de fricción

Definiendo el coeficiente de resistencia de rozamiento local Cf como

y sustituyendo la ecuación (4.12) para ts da

De la ecuación (4.11), es aparente que el gradiente de velocidad adimensional

A partir de esta ecuación dimensional, obtenga relaciones no dimensionales para los

SOLUCIÓN De la ecuación (4.13), el coeficiente de fricción local es

Integrando el valor_local y dividiendo entre el área por ancho unitario (L * 1) da el

y por tanto el coeficiente de fricción promedio es

4.5.2 Número de Nusselt

kf (Tq - Ts) 0T* kf 0T*

4.6 Evaluación de los coeficientes de transferencia de calor por convección

el fenómeno. Este método se prefiere cuando el fenómeno se puede describir mate-

4.7 Análisis dimensional

4.7.1 Dimensiones primarias y fórmulas dimensionales

TABLA 4.1 Cantidades físicas importantes de transferencia de calor, masa

4.7.2 Teorema p de Buckingham

En un problema que comprenda cinco cantidades físicas y tres dimensiones prima-

4.7.3 Determinación de grupos adimensionales

TABLA 4.2 Cantidades físicas pertinentes en transferencia de calor

y se sustituyen las fórmulas dimensionales

Es evidente que cualquier conjunto de valores de a, b, c, d y e que satisfaga simul-

bles originales se han combinado en tres grupos adimensionales. De acuerdo con la

y ahora los datos experimentales se pueden correlacionar en términos de tres varia-

4.7.4 Correlación de datos experimentales

1.0 Símbolo Fluido Diám. (mm) Ref

NuD = Pr0.3 * 0.82 * Re0.4

4.7.5 Principio de similitud

TABLA 4.3 Grupos adimensionales de importancia en la transferencia de calor

Factor j de Colburn j ( jH) StPr 2/3 Coeficiente de transferencia de calor

4.8* Solución analítica para el flujo laminar de capa límite sobre