Incropera 01
Incropera 01
Conducción en estado
transitorio
212 C a p ít u lo ó ■ i'.undarción en rsta d o transitorio
5.1
M é t o d o d e la resisten cia in terna d e s p r e c i a b le
/< 0
7 = T,
\
\
S
\
Liquidó \
~ i- “ f/ > 0
T* < Tl T = T( t y
cialm ente uniform e en cualquier instante durante el proceso transitorio. Esta suposición
implica que los gradientes de temperatura dentro del sólido son insignificantes.
De acuerdo con la ley de Fourier, la conducción de calor en ausencia de un gra-
diente de temperatura implica la existencia de una conductividad térmica infinita. Esta
condición es claramente imposible. Sin embargo, aunque la condición nunca se satisfa-
ce de forma exacta, se acerca mucho a ello si la resistencia a la conducción dentro del
sólido es pequeña comparada con la resistencia a la transferencia de calor entre el sóli-
do y sus alrededores. Por ahora suponga que, de hecho, éste es el caso.
Al no tomar en cuenta los gradientes de temperatura dentro del sólido, ya no es po
sible considerar el problema desde dentro del marco de la ecuación de difusión de calor.
En su lugar, la respuesta de temperatura transitoria se determina realizando un balance
global de energía en el sólido. Este balance debe relacionar la velocidad de perdida de
calor en la superficie con 1a rapidez de cambio de la energía interna. Al aplicar la ecua-
ción 1.11a al volumen de control de la figura 5.1, este requerimiento toma la forma
-E s e le = E al (5.1)
o
dT
-h A J L T - T J = p V c (5.2)
dt
Al introducir la diferencia de temperaturas
6 = T - T a (5.3)
donde
(5 4)
(5.5)
hAs e
La ecuación 5.5 sirve para determinar el tiempo que requiere el solido para alcanzar al-
guna temperatura T o, a la inversa, la ecuación 5 6 es útil para calcular la temperatura
que alcanza el sólido en algún tiempo /.
Los resultados anteriores indican que la diferencia entre las temperaturas del sóli-
do y el fluido deben decaer exponencialmente a cero conforme t se aproxima a infinito.
DEPARTAMENTO DE BIBLIOTECA
Unlvssrsiddü o.i ííJm [Link] *
14 C a p it u lo 5 ■ ( ondiiccion en estad o transitorio
xt 1 xi. 2 xt 3 xt 4
Flfcl l<A 5 . 2 Kt-spiiota ilc temprnitura tranriloria de Milido* de resistencia interna despreciable
que coi responden a diferente** constanle-s (('iium us de tiempo T,.
(? = J Q d i = hA s | O dt
' / '
(? = (pV< )0, 1 -e x p
/z
- Q = A E alm (5
5.2
Validez del m é to d o d e ln resistencia interna d e sp r e c ia b le
De los resultados precedentes es fácil ver por qué hay una íuerte preferencia por el uso
del método de la resistencia interna despreciable Es en verdad el método más sencillo
y conveniente para resolver problemas de conducción transitoria Por ello es importan
te determinar en que condiciones se puede usar con precisión razonable
Para desarrollar un cntcno adecuado considere la conducción en estado estable a
través de una pared plana de arca A (figura 5 3) Aunque estamos suponiendo condicio
nes de estado estable, este criterio se extiende fácilmente a los procesos transitorios
Una superficie se mantiene a una temperatura T s | y la otra se expone a un Huido de
temperatura T «, < Ts [. La temperatura de esta ultima superficie será algún valor inter-
medio. Ts 2, para el que Tx < Ts 2 < Ts | De aquí, en condiciones de estado estable, el
balance de energía de la superficie, ecuación 1 12, se reduce a
y - (Js., - r , -2) = h
donde k es la conductividad térmica del sólido. Al reacomodar, obtenemos
^conv
F h u ji ia 5.5
r h Efecto tlel numen de Biot en la distribución
de temperaturas ríe estado establr en una pared plana
í í con conver <ión en la superficie.
DEPARTAMENTO DE BIBLIOTECA
Universidad Simón 0 '♦|4
210 Cupitulo ■ C o n d u c c ió n en e s ta d o tra n sito rio
r. h
Fvci K.-\ 5 . 1 Distribución de temperaturas transitorias para «liferentes inum-ros de Biot en mmpar
plana enfriada simétricamente mediante convección.
= 0 1
k
el error asociado con el uso del método de la resistencia interna despreciable es peq*
no Por sencillez, se acostumbra definir la longitud característica de la ecuación5u|
como la relación entre el volumen del sólido y el área de la superficie, Lc a VA.
definición facilita el cálculo de Lc en solidos de forma complicada y reduce a lai
el espesor L para una pared plana de espesor 2L (figura 5.4). a r(J 2 para un cilindrol
go y a rQt3 para una esfera. Sin embargo, si se desea aplicar el criterio en forma|
dente, Lc debe asociarse con la escala de longitud que corresponde a la dife
máxima de temperaturas espaciales. En consecuencia, para una pared plana de esp
2L calentada (o enfriada) de forma simétrica, Lc permanecería igual a la mitad i
pesor L. Sin embargo, para un cilindro o esfera largos, Lt sería igual al radio realrj
lugar de r(J2 o rJ3.
5 .2 ■ Validez del m étodo de la resistencia interna despreciable 217
hAj
= Bi • Fo (5.11)
pVc
donde
cef
TT <5 1 2 )
se denomina número de Fourier. Ks un tiempo sin dimensión que, junto con el numero
de Biot, caracteriza los problemas de conducción transitoria. Al sustituir la ecuación
5.11 en la 5 6, obtenemos
q y y
— = — ~ = exp { - B i • Fo) (5.13)
v/| ■
*j * oc
• Ej e m pl o 5 . 1
Una unión termopar. cuya forma se aproxima a una esfera, se usará para la medición
de la temperatura en un (lujo de gas. Se sabe que el coeficiente de convección entre la
superficie de unión y el gas es h = 40 W/m- • K, y que las propiedades termofísicas de
la unión son k = 20 W/m • K. c = 400 J/kg • K, y p = 8500 kg/m . Determine el diá-
metro de la unión necesario para que el termopar tenga un tiempo constante de 1 s. Si
la unión está a 25°C y se la coloca en un flujo de gas que está a 200°C, ¿cuánto tiempo
tardará la unión en alcanzar 199°C?
So l í c ió n
E n c o n tr a r :
1. Diámetro de la unión necesario para una constante de tiempo de 1 s.
2. Tiempo que se requiere para alcanzar 199°C en un llujo de gas a 200°C.
E squem a:
Conductores
K -D -H
d epa r t a m en t o de b ib l io t ec a
Universidad Simón Bolívar - Sede del Litoral
218 C a p it u lo 5 ■ Conduci ion en estallo transitorio
S u p o s ic i o n e s :
1. La temperatura de la unión es uniforme en cualquier instante.
2. El intercambio de radiación con los alrededores es insignificante.
3. Las perdidas por conducción a través de los alambres de conducción son insignifi-
cantes.
4. Propiedades constantes.
Análisis:
I. Como se desconoce el diámetro de la unión, no es posible comenzar la solucml
determinando si se satisface el criterio para usar el método de la resistencia inlernj
despreciable, ecuación 5.10. Sin embargo, un procedimiento apropiado es usard
método para encontrar el diámetro y después determinar si se satisface el criten*
De la ecuación 5.7 y del hecho de que \ irD2 y \ = 7tDV6 para una estera, sj
sigue que
1 p irD
T, = 7 X — -— c
h i t D~ 6
Al reacomodar y sustituir los valores numéricos,
6 /j t , 6 X 400 W /m 2 • K X 1 s
D = = 7.06 X 10 - 4 m
pe 8500 kg/m 1 X 400 J/kg • K
Con L( = r (,/3 se sigue de la ecuación 5.10 que
h { r j 3) 400 W /m 2 • K X 3.53 X 10 ~4 m
Bi = = 2.35 X 10~4
k 3 X 20 W /m • K
En consecuencia, la ecuación 5.10 se satisface (para L = ra, así como par
r^/3) y el método de la resistencia interna despreciable sirve para una extcli
aproximación.
2. De la ecuación 5.5, el tiempo que se requiere para que la unión alcance!
199°C es
p { 7 r D 3l 6 ) c T , - T„ pD c T, - T ,
t = —r : — — In
h(7TD¿) T - 71 u r ' ln T -T
8500 kg/m 3 X 7.06 X 10 " 4 ni x 400 J/kg • K 25 - 200
t = ln
6 X 400 W /m 2 • K 199 - 200
t = 5.2 s ~ 5 t,
La transferencia de calor debida al intercambio de radiación <
C o m e n ta r io s :
unión y los alrededores y la conducción a través de los alambres conductores afee
tiempo de respuesta de la unión y da. de hecho, una temperatura de equilibrio i
itere de 7*.
5 .3
Análisis g e n e r a l d e l m é t o d o d e re siste n c ia in tern a d e s p r e c i a b le
Alrededores
T* “
r. i . no> = t ,
Vrad
E g. E aipn
o
r oo
K u .i h a 5 .5
do se separa de sus alrededores mediante un gas o un vacío Si las temperaturas del so-
lido y los alrededores difieren, el intercambio de radiación ocasionaría que cambie la
energía térmica interna y por ello la temperatura del sólido. Los cambios de temperatu-
ra también podrían inducirse aplicando un flujo de calor a una parte de la superficie o a
toda ella y/o iniciando la generación de energía térmica dentro del solido. 1 I calenta-
miento de la superficie se aplicaría, por ejemplo, al unir un calentador eléctrico de pe-
lícula o lamina a la superficie, mientras que la energía térmica se generaría haciendo
pasar una corriente eléctrica a través del sólido.
La figura 5.5 describe una situación en que las condiciones térmicas dentro de un
sólido estarán influidas de manera simultánea por convección, radiación, un flujo de
calor aplicado a la superficie y la generación de energía interna. Se supone que. inicial-
mente (t = 0). la temperatura del sólido (T¡) difiere de la del fluido, ÍT*. y de la de los
alrededores. Ta]m, y que se inicia tanto el calentamiento superficial como el global (q"s y
cj). El flujo de calor q" y la transferencia de calor por convección-radiación ocurren en
partes mutuamente excluycntes de la superficie, As{h) y A s{{. f), y se supone que la trans-
ferencia de convección-radiación se produce desde la superficie Además, aunque la
convección y la radiación se establecen para la misma superficie, las superficies pue-
den, de hecho, diferir (As( ^ Av ). Aplicando la conservación de la energía en cual-
quier instante t, se sigue de la ecuación 1.1 la que
dT
h Ég (¿/conv ífr a d W s (c , r) pV c ^ (5.14)
‘Se puede obtener una solución aproximada de diferencia finita disaetizando la derivativa de tiempo (sección 5.9) y hm u-nJo
axan'ar la solución en el nempo
DEPARTAMENTO d e [Link]
Universidad Simón P-M.y».' ° ^ ’ ¡torsf
220 Capítulo 5 ■ C o n d u cció n en e s ta d o tran sitorio
Al separar variables c integrar desde la condición inicial hasta cualquier tiempo /, se si-
gue que
rT dT
pVC J0 JTi
•pA __ •pA (5.1;
1 alr 1
AI evaluar ambas integrales y reacomodar, el tiempo que se requiere para alcanzar fe|
temperatura T se convierte en
pVc ^alr + T 7 alr + Tt
t = ln — ln
4eAs>ro T alr ^ a lr
— T T1 alr — T1 i
T
+ 2 tan - i — tan - i I I
(5.18|
alr ^"alr
Esta expresión no sirve para evaluar 7 de forma explícita en términos de '■ T¡, y T„„
se reduce de manera fácil al resultado límite para / a|r = 0 (radiación al espacio). Dere-i
greso a la ecuación 5.17 se muestra fácilmente que, para T ü\t ~ 0 ,
pV c 1 1
t = (5.
3 e A Sm¿ r \ T 3 7?,
Una solución exacta a la ecuación 5.15 también se obtendrá si es posible no toi
en cuenta la radiación y si h es independiente del tiempo. Si se introduce una tempei
tura reducida, 0 = T — 7’00, donde dO/dt = dT id t , la ecuación 5.15 se reduce aui^
ecuación diferencial lineal de primer orden no homogénea, de la forma
dd
+ a6 — b = 0 (5:
~dt
&= e— (5.:
a
Al reconocer que d d /d t = dO/dt, la ecuación 5.21 se sustituye en (5.20) lo que da
d6r
+ a& = 0 (52
dt
Por tanto.
T - T x b/a
= exp ( ~ a t ) + [1 — exp ( ~ a t )] (5.1
T. - 71
5 .3 ■ Análisis general del melada de resistencia interna despreciable 221
Ej e m pl o 5 .2
Son « ióín
E squem a:
2 1 = 3 mm *1
ü I I I I
Recubrimiento
— epoxico. h r . = 25°C
e = 0.8
.5a p o s ie io n e s :
1. La temperatura del panel es uniforme cn cualquier instante.
2. La resistencia térmica del recubrimiento epoxico es insignificante.
3. Propiedades constantes.
•? •> C a p ít u lo ó ■ Conducción en estado transitorio
f d r = TU) - T¡ =
T,
1
f
p cL Jo
[h( T -
C o m e n ta r io s :
1. Por lo general, la precisión de una integración num érica mejora al dism inuir A t,
pero a expensas del aum ento del tiem po de cálculo. Sin em bargo, en este caso los
resultados que se obtienen para A/ = 1 s son idénticos a los que se obtienen para
A t = 10 s, lo que indica que el intervalo de tiem po más largo es suficiente para des-
cribir con exactitud la historia de la tem peratura.
2. La duración del proceso de dos pasos se puede reducir al aum entar los coeficientes
de convección y/o dism inuir el periodo de calentam iento prolongado. La segunda
opción es posible por el hecho de que, durante una parte del periodo de enfria-
m iento. la tem peratura del panel perm anece por arriba de 150°C. De aquí, para sa-
tisfacer el requerim iento de curado no es necesario prolongar el calentam iento
tanto com o 5 m inutos a partir de t = tc. Si los coeficientes de convección aum en-
tan a ha = hc = 100 W /m 2 • K y se m antiene un periodo de calentam iento prolon-
gado de 300 s, la integración num érica da tc = 58 s y t, = 445 s. El intervalo de
tiem po correspondiente sobre el cual la tem peratura del panel excede I50°C es
Af(7> | 5o°c) = 306 s (58 s < t ^ 364 s). Si el periodo de calentam iento prolongado
se reduce a 294 s, la integración num érica da tc = 58 s, /, = 439 s y Af( /> i^()cC) =
300 s. Así, el tiem po total de proceso se reduce, m ientras que el requerim iento de
curado aún se satisface.
•tos espaciales
Con frecuencia surgen situaciones para las que el m étodo de la resistencia interna des-
preciable no es apropiado y deben usarse m étodos alternativos. Sin im portar el m étodo
particular, ahora debem os enfrentar el hecho de que los gradientes dentro del m edio ya
no son insignificantes.
En su form a más general, los problem as de conducción transitoria se describen
m ediante la ecuación de calor, ecuación 2 .13 para coordenadas rectangulares o las ecua-
ciones 2.20 y 2 23, respectivam ente, para coordenadas cilindricas y estencas La solución
a estas ecuaciones diferenciales parciales proporciona la variación de la tem peratura
con el tiem po y con las coordenadas espaciales Sin em bargo, en m uchos problem as,
com o el de la pared plana de la figura 5.4, solo se necesita una coordenada espacial
para describir la distribución interna de tem peraturas. Sin generación interna y con h
suposición de conductividad térm ica constante, la ecuación 2.13 se reduce entonces a
d 2T 1 dT
— = - — (5.26)
dx a at
Para resolver la ecuación 5.26 en cuanto a la distribución de tem peraturas T { \, t ).
es necesario especificar una condición iniciaI y dos condiciones de frontera. Para el
problem a de conducción transitoria típico de la figura 5.4, la condición inicial es
7'(a, 0) = Ti (5 27)
dT
= 0 (5.28)
dx r =0
DEPARTAMENTO DE BIBLIOTECA
U n iv e r s id a d » IhJ.i [Link] - Sedo „ (.ora1
221 C apitulo .> ■ ( ontltieeión en es lado Lrunsiiorit*
d r
-k = h v a u t) - t j (5.29;
dx x^L
I a ecuación 5 2 presupone una distribución de tem peraturas uniform e en el tierp
t = 0; la ecuación 5.28 refleja el requerim iento ele sim etría para el plano medio de
la pared; y la ecuación 5.29 describe la condición de superficie experimentada en e*
tiem po / > 0 De las ecuaciones 5 26 a 5 29. es evidente que. adem as de depender
y r, las tem peraturas en la pared tam bién dependen de un num ero de parámetros
eos. Ln particular
e * , l = T - r~ ix
o. T ¡~ T t
.v* = - (5.;
L
donde L es la m itad del espesor de la pared plana, y un tiem po adim ensional sci
com o
at
r** = —r rr
= Fo (5.:
ú
0
d2 * de*
dx*2 dF o
y las condiciones inicial y de frontera son
6*(x * 0 ) = 1 ój
de*
= 0
dx*
x *= 0
de*
-B i e*(\,t*)
5 .5 ■ Pared plana ron convección 225
donde el núm ero de Biol es Bi = hL/k. Bn form a adim ensional la dependencia funcio-
nal se expresa ahora com o
= f( \* ,F o ,B i) (5.38)
R ecuerde que una dependencia funcional sim ilar, sin la variación v*. se obtuvo para el
m étodo de la resistencia interna despreciable, com o se m uestra en la ecuación 5.13.
Al com parar las ecuaciones 5.30 y 5.38, la considerable ventaja asociada con el
cam bio del problem a a una form a adim ensional se hace evidente. La ecuación 5.38 im -
plica que para una geometría establecida , la distribución de temperaturas transitoria
es una fu m ió n universal de x*. Fo y Bi. Es decir, la solución adimensional supone una
form a establecida que no depende del valor particular de 7’,, Tx , L , k, a o h. C om o esta
generalización sim plifica m uchísim o la presentación y utilización de soluciones transi-
torias, las variables adim ensionales se usan de manera intensiva en las secciones poste-
riores.
5 .5 .1 S o lu c ió n e x a c ta
C onsidere una pared plana de espesor 2 L (figura 5.6í/). Si el espesor es pequeño en re-
lación con el ancho y la altura de la pared, es razonable suponer que la conducción
ocurre exclusivam ente en la dirección x. Si la pared al principio está a una tem peratura
r 7tv. 0) = Tj r(r. 0 1 = 7 ,
(ü j)
í í t t í t
L-
(a) ib)
DEPARTAM ENTO d e b ib l io t e c a
U n iv e r s id a d S im ó n Bol iva-* S e d e d e l L lto ra
C apitulo 5 ■ ( onflucción en rstadu transitoria
4 sen Cn
C = (5.39b
2Cn + sen ( 2 £,)
y los valores característicos (e ige m a l ores) de £„ son las raíces positivas de la ecuación
trascendente
ln ta n ln = F i (5.39c]
Las prim eras cuatro raíces de esta ecuación se dan en el apéndice B 3.
5 .5 » Solución aproximada
Se puede m ostrar (problem a 5.27) que para valores de Fo > 0.2, la solución en seriein
finita, ecuación 5.39a. se aproxim a con el pr m er térm ino de la sene. Al recurrir a esta
aproxim ación, la form a adim ensional de la distribución de tem peraturas se convierteei
tí* = (£ c o s (Cix*)
donde tí* ( T0 — TX)/(T, — 7») representa la tem peratura del plano m edio (.v* - 0
Ln m uchas situaciones es útil conocer la energía total que dism inuye en la pared<
cualquier tiem po t cn el proceso transitorio. El requerim iento de conservación dej
energía, ecuación 1 1 Ib. se aplica al intervalo de tiem po lim itado por la condición i
cial (r = 0 ) y cualquier tiem po t > 0
5 .5 ■ Pared plana con convección 227
ÜB¡ = hL'k para la pared plana y hr0!k para el cilindro infinito y la esfera. Véase la figura 5.6.
d epa r t a m en t o d e b ib l io t e c a ^
Universidad Simón Bollver Sed en
Capitule» 5 ■ Conducción en oslado transitorio
Al igualar la energía que se transfiere desde la pared Q a £ sale y con £ cnl 0 y Afalm =
E(t) — E{0 ), se sigue que
Q = - [ £ ( / ) - £(0)] (5.43a)
Q= - fx [T{r, t) - T, 1 dV (5 43b)
(5441
Q_ - \ T { x , t) - T,] ü
dVv 1 r
— = — (1 - d*) d V (5.45Í
Qo - J T -T * V v Jy
Al em plear la forma aproxim ada de la distribución de tem peraturas para la pared plana,
ecuación 5.40b, la integración que establece la ecuación 5.45 se ejecuta para obtener
Q_ = _ sen jj_
(54
Qo í.
donde 0 * se determ ina de la ecuación 5.41, con la ayuda de la tabla 5.1 para los valor
de los coeficientes C\ y
Com o el problem a m atem ático es precisam ente el m ism o, los resultados preceden!
tam bién se aplican a una pared plana de espesor L. la cual está aislada en un lado (jt*
0) y experim enta transporte convectivo en el otro ( a * = + 1 ). Esta equivalencia es u
consecuencia del hecho de que, sin im portar si se establece un requisito simétrico
adiabático en .v* = 0 , la condición de frontera es de la form a dd*/clx* = 0 .
También debe advertirse que los resultados anteriores sirven para determinar
respuesta transitoria de una pared plana a un cam bio súbito en la tem peratura de la
perficie. El proceso es equivalente a tener un coeficiente de convección infinito,
cuyo caso el num ero de Biot es infinito (Bi = o°) y la tem peratura del fluido 7*
reem plaza por la tem peratura establecida de la superficie 7\.
Finalm ente, observam os que las representaciones gráficas de las aproximación
de un térm ino ya se han desarrollado [5, 61 y se presentan en el apéndice D Aunque
gráficas asociadas proporcionan un m edio conveniente para resolver problemas de c
ducción unidim ensional transitoria para Fo > 0 .2 , se logra m ayor precisión medi
las ecuaciones 5.40 y 5.46.
5.Í» ■ Sistemas radiales ron convección 229
5.0
Sistemas radiales con convección
Para un cilindro infinito o una esfera de radio r0 (ligura 5.6 b), que está a una tem pera-
tura inicial uniform e y experim enta un cam bio en las condiciones de convección, se
producen resultados sim ilares a los de la sección 5.5. t s decir, es posible una solución
en serie exacta para la dependencia con respecto al tiem po de la distribución radial de
tem peraturas, y se aprovecha la aproxim ación de un term ino para la ma>oria de las
condiciones. El cilindro infinito es una idealización que perm ite la suposición de co n -
ducción unidim ensional en la dirección radial. Esta es una aproxim ación razonable p a-
ra cilindros con U ra 2 : 10 .
5 * 6 .1 Soluciones exactas
e* = X c „ ex p ( ~ ( „
2 F o )U Í„ r* ) (5.47a)
«= I
donde Fo = at!r2O»
JÁ U
C » = — ' . 2/- >**\ i . 2 / ** * (5 4 ? b)
nJl i L
C ) +
y los valores característicos de son las raíces positivas de la ecuación trascendental
M C n)
L 77TT = (5 .4 7 0
J0\bn)
Las cantidades J i y J0 son funciones de Bessel de prim era clase y sus valores se tabu-
lan en el apéndice B.4. S chneider [2] tabuló las raíces de la ecuación trascendental
(5.47c)
donde Fo = ctrlrlo*
4 [ s c n ( ¿ , ) - Cn c ° s ( £ , ) ]
C = — — —------------------------------------------------------- (5.48b)
2 £ , - s c n ( 2 £,)
1 - £„ cot = Bi (5.48c)
5 .6 .2 Soluciones aproximadas
Para el cilindro infinito y la esfera, las soluciones en serie anteriores se aproximan nue-
vam ente m ediante un solo térm ino para Fo > 0.2. De aquí, com o para el caso de la pa-
red plana, la dependencia respecto al tiem po de la tem peratura en cualquier lugar
dentro del sistem a radial es la m ism a que la de la línea central o el punto central.
C ilindro in fin ito La aproxim ación de un térm ino para la ecuación 5.47 es
0* = Ci ex p (~ C \F o)J0(£\r*)
Los valores de los coeficientes C | y í | ya se han determ inado y se enum eran en la tabla
5.1 para un rango de núm eros de Biot.
o
1
s ir
donde 0 * representa la tem peratura del centro y es de la forma
0 * = C l ex p {~ C 2\Fo)
Los valores de los coeficientes C\ y £i ya se han determ inado y se enum eran en la tabla
5.1 para un rango de núm eros de Biot.
E sfera
Q 3 0*
=r = 1 - [ s e n ( f ,) - Ci eo s ( f,) ]
Qo íi
Los valores de la tem peratura del centro 0* se determ inan a partir de la ecuación 5.i
o 5.50c, con los coeficientes de la tabla 5.1 para el sistem a apropiado.
5 .6 ■ Sistemas radiales coa convección 231
5 .6 .1 Consideraciones adicionales
C om o en el caso de la pared plana, los resultados anteriores son útiles para predecir la
respuesta transitoria de cilindros largos y esferas sujetos a un cam bio súbito en la tem -
peratura de la superficie. Esto es, se establece un núm ero de Biot infinito, y la tem pera-
tura del fluido Toe se reem plaza con la tem peratura constante de la superficie T .
En el apéndice D se m uestran representaciones gráficas de aproxim aciones de un
térm ino.
Ej k m pl o 5 .3
C onsidere una tubería de acero (AISI 1010) que tiene 1 m de diám etro interno y una
pared con espesor de 40 mm La tubería está fuertem ente aislada en el exterior y, antes
del inicio del flujo, las paredes de la tubería se encuentran a una tem peratura uniform e
de —20 C. Con el inicio del flujo se bom bea aceite caliente a 60°C por la tubería, con
lo que se crea una condición convectiva de superficie que corresponde a /; = 500 W /m 2 • K
en la superficie interior de la tubería.
1. 6C uáles son los núm eros de Biot y de Eourier apropiados, 8 minutos después de
iniciado el flu jo 7
2. A t = 8 min, ¿cuál es la tem peratura de la superficie externa cubierta por aislante?
3. ¿Cual es el flujo de calor q" (W /m ) a la tubería desde el aceite en / = 8 m inutos 7
4. ¿C uanta energía por m etro de longitud de tubería se ha transferido del aceite a la
tubería en t = 8 m inutos 9
S o u c ió in
E n co n tra r:
1. N úm eros de Biot y de F ourier después de 8 minutos
2. Tem peratura de la superficie externa de la tubería después de 8 minutos
3. Flujo de calor a la pared en 8 m inutos.
4. Energía transferida a la tubería por unidad de longitud después de 8 minutos.
E squem a:
Ti\. 0) =
T¡ = -20°C T {L, t)
no, t) = 60°C
h = 500 W/m2 • k
Aislante
Acero AISI 1010
l. = 40 mm
X
d e pa r t a m e n t o d e b ib l io t e c a
Unlvera dad Simón Boliv». Sode ó-* i í<*r
Capítulo 5 ■ C o n d u cció n en e s ta d o tru n sitorio
Suposiciones:
1. La pared de la tubería se aproxim a com o una pared plana, pues el espesor es mu-
cho m enor que el diám etro.
2. Propiedades constantes.
3. La superficie externa de la tubería es adiabática.
A n álisis:
1. En / = 8 m inutos, los núm eros de B iot y de Fouricr se calculan de las ecuac one^
5.10 y 5.12, respectivam ente, con Lc = L. Así,
hL 5 0 0 W /m 2 • K X 0 .04 m
Bi = = 0.313 <
k 6 3 .9 W /m • K
T —T
0* = — — - = C , ex p ( - f i F o )
* Z -
8* = 1.047 ex p [ - ( 0 .5 3 1 ra d )2 X 5.64] = 0 .2 1 4
q ’^ U 4 8 0 s ) ^ q ' [ = h[T(L , 4 8 0 s) - T J
Con el uso de la aproxim ación de un térm ino para la tem peratura de la sup
la ecuación 5.40b con v* = 1 tiene la form a
0* 6* eos (£,)
T(L , t) T„ + (T, - T„)0* eos (£,)
T{L , 8 m in) 6 0 °C + ( - 2 0 - 6 0 )°C X 0 .2 1 4 X eos (0.531 rad)
T(L, 8 m in) 4 5 .2 °C
5 .6 ■ Sistemas radiales coa convección 233
se sigue que
C o m entarios:
1. El signo de m enos que se asocia con q" y Q' im plica sim plem ente que la dirección
de la transferencia de calor es del aceite a la tubería (en la pared de la tubería).
2. Los resultados anteriores tam bién se obtienen de la aplicación de las gráficas de
H eisler y G rober del apéndice D. Por ejem plo, si se usa la figura D. 1 con Bx 1 =
3.2, se sigue que 0* ~ 0.22 y el valor correspondiente de la tem peratura del plano
m edio es T0 ~ 42°C. Para a * = 1 y B r x — 3 2, la figura D 2 da 6{L, 8 m in)/
0O(8 m in) ~ 0 8 6 , de donde se sigue que T(L, 8 mm) ~ T + 0 8 6 [Y„(8 min)
T, J - 45 °C y q'¡ = - 7 5 0 0 W /m 2. Con Bi = 0.313 y B rF o = 0 55, la figura D 3
da Q Q(>= 0.78. AI sustituir de la ecuación 5.44, se sigue que
Los resultados precedentes están de acuerdo con los que se obtuvieron directa-
m ente de las aproxim aciones con un solo term ino
E j k m pl o 5 .4
P aso 1 Enfriam iento en aire a 20°C durante un periodo de tiem po ta hasta que la tem -
peratura del centro alcanza un valor crítico, 7 \/0 . /tí) = 335°C Para esta situación,
el coeficiente de calor convectivo es ha — 10 W /m ■ K.
DEPARTAMENTO DE 8 I 8 L IO »_oA
Universidad S im ó n B o lív a r Spde d e l L itoral
C apítulo r> ■ ( omlncción en estallo transitorio
D espués de que la esfera alcanza esta tem peratura critica, se inicia el segundo paso.
F aso 2 Enfriam iento en un baño de agua muy agitado a 20°C, con un coeficiente de
transferencia de calor por convección h * = 6000 W /m 2 • K.
Las propiedades term ofísicas del material son /> = 3000 kg/m 3, k = 20 W m • K. c =
1000 J/kg • K y a = 6.66 X 10 ~6 m 2/s.
1. C alcule el tiem po ta que se requiere para que se com plete el paso 1 del proceso dej
enfriam iento.
2. C alcule el tiem po tw que se requiere durante el paso 2 del proceso para que el c i-j
tro de la esfera se enfríe de 335°C (condición al final del paso 1) a 50°C.
Sol liC lÓ N
E n c o n tra r:
1. T iem po ta que se requiere para llevar a cabo el enfriam iento que se desea en aire
2. T iem po tw que se necesita para com pletar el enfriam iento en el baño de agua.
E sq u em a :
Tx = 20°C U = 20°C
ha = 10 W/m2 -K hw = 6000 W/m? •K
Aire » Agua I
Esfera, r0 = 5 mm
p = 3000 kg/m3
T, = 400°C c = l '.L ?/ I u 7/ = 335°C
I T J 0 . T J = 335°C “ I m /S T J O, T j = 50°C
Paso 1 Paso 2
S u p o sicio n es:
1. C onducción unidim ensional en r.
2. Propiedades constantes.
A nálisis:
1. Para determ inar si es posible utilizar el m étodo de la resistencia interna despn
ble, se calcula el num ero de Biot. De la ecuación 5.10, con Lc = r j 3,
har0 10 W /m 2 • K X 0.005 m -4
Bi = = 8.33 X 10
3k ~ 3 X 20 W /m " K
2. Para determ inar si el m étodo de la resistencia interna despreciable tam bién sirve
en el segundo paso del proceso de enfriam iento, de nuevo se calcula el num ero de
Biot. En este ca so
K r0 6 0 0 0 W /m 2 • K X 0.005 m
Bi = = 0 .50
3k 3 X~20 W /m ~ K
donde rH. = Fo r 2 /ot. Con el núm ero de Biot definido ahora com o
hwr„ 6 0 0 0 W /m 2 • K X 0.005 m
Bi = = 1.50
~T~ ~ 2 0 W /m " K
1 (50 - 20)°C
Fo = 4 -! x = 0 .82
(1 .8 0 0 rad) 37 6
[ 1.3 (335 - 2 0 )°C
(0 005 m)"
= F o — = 0.8 2 = 3.1 s
a 6 .6 6 X 10 6 m 2/s
A dvierta que, con Fo = 0.82, se justifica el uso de la aproxim ación con un término.
T(ra) - 7 , 0.095
ti* ( O = sen (1 .8 0 0 rad) = 0 .0 5 1 4
7 ,- 7 , 1.800 rad
d e pa r t a m e n t o d e b ib l io t e c a
U n iv e r s id a d S im a n b o n « j r - S e d * - .. >rt>
236 C apítulo 5 ■ Conducción en estado transitorio
de la serie. Para (/ — ta) > 0.75 s. la aproxim ación con un térm ino proporciona
una convergencia satisfactoria. AI calcular y presentar en form a de gráfica las
torias de las tem peraturas para r = 0 y r = r 0, obtenem os los siguientes resultado
para 0 ^ (/ — ta) ^ 5 s:
. Las gráficas de H eislcr del apéndice D tam bién sirven para analizar el procet»j
paso 2. Con Bi 1 = 0.67 y 0* = 0.095, la figura D.7 da Fo ~ 0.8 en cuyoi
tw *== 3.0 s. De la figura D . 8 , con r* = I, dir0)l6o ~ 0.52, en cuyo caso T(/*);
20°C + 0.52(50 - 20)°C « 36°C.
5 .7
Sólido semiinfinito
O tra geom etría sim ple para la que es posible obtener soluciones analíticas es e l »
semiinfinito. C om o tal sólido se extiende hasta el infinito en todas direcciones exea
una, se caracteriza por una sola superficie identificable (figura 5.7). Si se imponj
cam bio súbito de condiciones en esta superficie, ocurrirá una conducción unidimel
nal dentro del sólido. El sólido sem iinfinito proporciona una idealización útil paral
chos problem as prácticos. Se aprovecha para determ inar la transferencia de
transitoria cerca de la superficie de la tierra o para aproxim ar la respuesta transitorii
un sólido finito, com o una losa gruesa. En cuanto a esta segunda situación la aprc
ción sería razonable para la prim era parte del transitorio, durante la cual las tcnif
ras en el interior de la losa (a bastante distancia de la superficie) no son afectada* |
cam bio en las condiciones de la superficie.
La ecuación de ealor para la conducción transitoria en un sólido semiinfinita
dada por la ecuación 5.26. La condición inicial se establece m ediante la ecuación^
y la condición de frontera interior es de la forma
T U. t)
Ts
d T _ d T d t? _ 1 dT
dx d ii dx (4 a i ) 1/2 drj
d2T d dT_ 3t7 _ 1 d 2T
djc2 d7] dx dx 4 at d if
d T _ d T d 77 _ x dT
dt dT) dt 2 t(4 a t)ul d i 7
d 2T dT
^ r " 2 íí^ ( 5 -5 4 )
T( t ) = 0 ) = TX (5.55)
y con x —>ce, así com o t= 0 , que corresponde a 17 —» oo,lacondición inicial yla con-
dición de frontera interior corresponden al único requerim iento
DEPARTAMENTO DE B IB L IO itu A
U n iv e r s id a d S im ó n R oU vs ^ « d e d e l L ito ra
238 C apítulo 5 ■ Conducción en estado transitorio
d{dT/dr¡)
= ~ 2 7] dr]
0dT/dr ) r
Al integrar, se sigue que
\n(dT/dr]) = - r f + C\
o
dT
= C, exp(—17 )
dr]
T = C, f e x p ( —u2) du + Ts
Jo
De la segunda condición de frontera, ecuación 5.56, obtenem os
T, = C, í e x p ( - w ) du + Ts
Jo
o, m ediante la evaluación de la integral definida,
2(Tt - Ts)
C, = 771/2
De aquí la distribución de tem peraturas se expresa com o
T -T s
— ~r = (2/ 7r 12) f ex p (—m2) du = erf 17 (52
T, ~ Tx jo
donde la función gaussiana de error . erf 17. es una función m atem ática estándar quei
tabula en el apéndice B El flujo de calor en la superficie se obtiene con la
de la ley de Fourier en v = 0, en cuyo caso
dT ¿/(erf 17 ) 8 7 7
<7" = -k (T , - TJ
~ k lh x= 0 dr) d.x h=o
™
q" k(Ts - T'J)(2 / 7r 1 2) e x p ( - i 72)( 4 a í) 1/2
k(Ts - Tt)
y" 15,!
( 7 r a í ) 1/2
Tam bién se pueden obtener soluciones analíticas para las condiciones superfic
del caso 2 y del caso 3, y los resultados para los tres casos se resum en como sigue.
«>.7 ■ Solida seniimjmilo 239
T(x, t ) - Ts
(5.57)
7 - 7,
(5.58)
V 7raí
2q"(athr) 1/2 v
exp erfc (5.59)
4a/ 2V a/
'
,)T
C aso 3 C onvección su p erficial: —¿ = h[Tx - 7(0, /)]
dx «=0
7 (a, /) - 7, ( * _
T ._ - T , _ e 'C V 2 V í r t ,
1.0
0.5
t». T o.i
i
v-. 8
005
Fuá ha 5 .8
H is to r ia s <fe te m p e ra tu ra s e n un
0.01 só lid o s e m iin fin ito co n c o n v e c c ió n
0.5 10 1.5
s u p e r fic ia l [2].
W laptada co n p e rm iso .
DEPARTAMENTO DE B lB L lO tto A
U n iv e r s id a d S im ó n B o lív a r - B o d e d e l L itoral
240 C apítulo 5 ■ Conducción en estado transitorio
F ic l r a 5 .9
Contacto de la interia/ enlre dos
sólidos semiinfinitos a diferentes
temperaturas iniciales.
o, al resolver para T\
T = (k p c)ü 2 + ( kpc)¿
1/2
Por tanto, la cantidad m = (kpc )' 2 es un factor de peso que determ ina si Ts aproa
más de cerca a TA ¡(mA > m H) o Tfí ,(mti > mA).
E j k m pl o 5 .5
diciones en las que el suelo, inicialm ente a una tem peratura uniform e de 20°C, se so-
mete a una tem peratura superficial constante de — 15°C durante 60 días?
S o l u c ió n
Se conoce: Tem peratura im puesta a la superficie del suelo inicialm cnte a 20°C.
E squem a:
Atmósfera T¡ = -15°C
Suelo
T, = 20°C
L_ Ux„„ 60d) - 0°C
Ned de tukeró He agua
S u p o s ic io n e s :
1. C onducción unidim ensional en ,v.
2. El suelo es un m edio semiinfinito.
3. Propiedades constantes.
A n á lis is: Las condiciones establecidas corresponden a las del caso 1 de la figura
5.7, y la respuesta transitoria de tem peratura del suelo está gobernada por la ecuación
5.57. Por tanto, en el tiem po t — 60 días después del cam bio de la tem peratura de la
superficie.
T(xm, t) - Ts / x \
T, ~ T, 6 l 2V 7 t j
o
0 — ( —15) „ t x,„
= 0.429 = erf f
20 -(-1 5 ) ' yrT a t
= 0.40
jc(m)
A medida que la penetración térm ica aum enta con el increm ento del tiempo, el
diente de tem peratura cn la superficie, dTI'óx[x =0 y, por tanto, la velocidad de extrae
ción de calor del suelo, dism inuyen.
5 .8
E fecto s m u liidim en sion ales
T h
T (r.x)
T
L L
L u
4 T ?
i
/
i.
1 d i dT\ d2T 1
r dr l dr ) dx2 ex dt
donde se utiliza x en lugar de z para designar la coordenada axial Se obtiene una solu-
ción en form a cerrada a esta ecuación m ediante el m étodo de separación de variables.
A unque no considerarem os los detalles de esta solución, es im portante advertir que el
resultado final se expresaría en la siguiente forma.
T(x%t) - r*
S(.x. t) = Sólido (5 64)
I w T semiinfinito
T(x, t) - Tu
P(x, t) = — ~ Pared (5 65)
Tt ~ T x plana
T(r, t ) - 7 ;
C (r, t ) = Cilindro (5.66)
T. - 7L infinito
La coordenada v para el solido sem iinfinito se m ide desde la superficie, m ientras que
para la pared plana se mide desde el plano medio. Al usar la figura 5.11 deben obser-
varse cuidadosam ente los orígenes coordenado*. La distribución tridim ensional transi-
S(A, t) /) C(r. t)
r
.l —
- 2L, -I /77J
(a) Sólido (¿) Pared plana U ) Cilindro infinito
semiinfinito
C(r,i)S(x. o
P( V|. t)P(\2. t)
k H
»
A’’' \->
Í T '
/ 2L\
Á -L
T
- 2 L 2—
(d) Placa (c) Barra rectangular (/) Cilindro
sem infinita infinita semi-infinito
C(r, nP(x. t)
5{a'3. t)P(.\ | , t)P(x2, t) P (ij .i)P(x2, f)P( t 3 ,r)
T (x{, x 2, s 3, t) ~
= P (x j, r) • P (x 2, r) • P(a 3, 0
T, - 7L
Las distancias .v,, .v2 y v3 se miden todas con respecto a un sistem a coordenador
guiar cuyo origen está en el centro del paralelepípedo.
La cantidad de energía Q transferida hacia o desde un sólido durante un p
de conducción transitoria m ultidim ensional tam bién se determ ina mediante la
nación de resultados unidim ensionales, com o m uestra Langston [7J.
5 .8 ■ Efectos miill¡dimensionales 245
Ej k m pl o 5 . 6
Son <:ió\
S e conoce: Tem peratura inicial y dim ensiones del cilindro, y tem peratura y condi-
ciones de convección de un baño de aceite.
E sq u em a :
r0 = 40 mm 7(0. L, t)
T(r. i. 0) = T, = 600 K
t T(ro, 0 . l )
L = 30 mm
L = 30 mm
i
IA T( 0 .0 . /)
Cilindro de
AISI 304
Baño d e ► ^ - ^
aceite ---- ► h = 5qo W/m2 • K
S u p o sicio n es:
1. C onducción bidim cnsional en r y x.
2. Propiedades constantes.
P ro p ied a d es: Tabla A .l, acero inoxidable, A ISI 304[T = (600 + 300)/2 = 450 K].
p = 7900 kg/m 3, c = 526 J/kg • K, k = 17.4 W /m • K, a = k/pc = 4.19 X 10 -6 m 2/s.
T(r, x , t) - Ta
= P(x , f)C (r, t)
Z - 7L
C apítulo 5 ■ Conducción cn estado transitorio
17 4 W /m ■ K
B rx= - = 1.16
hL 5 0 0 W /m 2 • K X 0.03 m
at 4 19 X 1 0 " 6 m 2/s X 180 s
Fo — — = 0 84
L} (0.03 m )2
0* = y = C, e x p ( - £ 2Fr>)
0o T{0, 3 m in) - Ta
= 1.109 exp [ - ( 0 .8 1 4 rad )2 X 0.84] = 0.636
T -T a P a re d
plana
k 17.4 W /m • K
B F X= = 0.87
hra 5 00 W /m 2 • K X 0 0 4 m
at 4 .1 9 X 1 0 " 6 m 2/s X 1 8 0 s
Fo = = 0.47
r2 (0 .0 4 m )2
efí
0* = y = C, ex p ( —£ ,F o )
La tem peratura en el centro de una cara circular se obtiene del requisito de que
0* 0
— = - = eo s (f,* * )
(KL) T(L , 3 m m ) - Ta
Pared = eos (0.814 rad X 1) = 0.687
0O 7(0, 3 m in) - 7 C plana
D e aquí
7 (7 , 3 m in ) - 7«
Pared = 0.687 X 0 .6 3 6 = 0.437
71 r - 7 1 oo plana
Por tanto,
7 (0 , L, 3 m m ) - 7 {
= 0.437 X 0 .5 5 0 = 0 .2 4 0
71 r - 7 1 oo
La tem peratura en la altura m edia lateral se obtiene del requerim iento de que
e~* -- J 6o - U, b, rr **\
)
0{ro) T(ra, 3 m m ) - 7
cdmdro — 70( 1.307 rad X 1) — 0.616
7 (0, 3 m in) - 7« infinito
De aquí
DEPARTAMENTO d e b ib l io t e c a
U n iv e r s id a d S im ó n Bol Iv a r - Sede d e l L ito ra '
218 Capitule» f» ■ Conducción en estada transitorio
C ornan ta rio s:
1. Verifique que la tem peratura en el borde del cilindro es T(ra, L, 3 min) = 344 K.
2. Las gráficas de H cisler del apéndice D tam bién servirán para obtener los resul
dos que se desean. Al usar estas gráficas, se obtiene BJ6,|pare<i 11 ** 0.64
^«^ílcilindro infinito 0 .5 5 , 0{L)!0C>|Pared plana 0 .6 8 , y 0 .6 1 , que eslímfo
acuerdo con los resultados que se obtienen de las aproxim aciones de un tennino.
5 .9
M étodos de diferencias finitas
5 .9 .1 Discretización de la ecuación
de calor: método explícito
i dr a 2r d2r
+ (5
a dt dx2 dy
Para obtener la form a en diferencias finitas de esta ecuación, podemos usar las
m ariones de clifeiencía central para las derivadas espaciales establecidas poi
ecuaciones 4.31 y 4 32. U na vez más los subíndices m y n sirven para designar1
siciones x y y de los puntos nodales discietos. Sin em barco, ademas de discre
espacio, el problem a debe discretizarse en el tiem po. El entero p se introduce:cc
i
propósito, donde
t = pAt
dT •TP+I _ pP
m, n 1 m, n
aT m.n At
turas asociada con los tiem pos nuevo {p + 1) y anterior ( p ). Por ello los cálculos deben
llevarse a cabo en tiem pos sucesivos separados por el intervalo A t, y com o una solu-
ción en diferencias finitas restringe la determ inación de tem peraturas a puntos discretos
en el espacio, tam bién la restringe a puntos discretos en el tiempo.
Si la ecuación 5.69 se sustituye en la ecuación 5.67, la naturaleza de la solución en
diferencias finitas dependerá del tiem po específico al que se evalúan las tem peraturas
en las aproxim aciones en diferencias finitas para las derivadas espaciales. En el método
explícito de solución, estas tem peraturas se evalúan en el tiem po anterior {p). Por esto,
la ecuación 5.69 se considera que es una aproxim ación en diferencias hacia adelante
para la derivada respecto al tiem po Al evaluar los térm inos en el lado derecho de las
ecuaciones 4.31 y 4.32 en p y sustituir en la ecuación 5 67, la form a explícita
de la ecuación en diferencias finitas para el nodo interior m, n es
i
1 1m ni _
-T P + T P
* mn
Ai
(5.70)
(5.71)
a Ai
(5.72)
Las ecuaciones 5 71 y 5 73 son explícitas pues las tem peraturas nodales descono-
cidas para el tiem po nuevo se determ inan de m anera exclusiva m ediante tem peraturas
nodales conocidas en el tiem po anterior Por ello el cálculo de las tem peraturas desco-
nocidas es directo. C om o se conoce la tem peratura de cada nodo interior en A = 0 (j> =
0 ) de las condiciones iniciales establecidas, los cálculos com ienzan en t = At(p = 1 ),
donde la ecuación 5 71 o 5 73 se aplica a cada nodo interior para determ inar su tem pe-
ratura. Con tem peraturas conocidas para t — A t. la ecuación en diferencias finitas apro-
piada se aplica entonces a cada nodo para determ inar su tem peratura en t = 2 A t(p
2). De esta form a, la distribución transitoria de tem peraturas se obtiene al avanzar en
el tiempo, con el uso de intervalos de At.
La precisión de la solución en diferencias finitas se m ejora dism inuyendo los valo-
res de Aa y Ai. Por supuesto, el núm ero de puntos nodales interiores que debe conside-
rarse aum enta al dism inuir Aa. y el núm ero de intervalos de tiem po que se requieren
para llevar la solución a un tiem po final establecido aum enta al dism inuir A t. Por ello
el tiem po de cálculo aum enta al dism inuir Aa y Ai. La elección de Aa norm alm ente se
basa en un com prom iso entre la precisión y los requerim ientos de calculo. Sin embar-
go, una vez que se hace esta selección, el valor de Ar tal vez no se elija de forma inde-
pendiente Se deteim ina, en realidad, m ediante requerim ientos de estabilidad
Una característica inconveniente del m étodo explícito es que no es incondici nal-
m ente estable. L n un problem a transitorio, la solución para las temperaturas nodales
debe aproxim arse de form a continua a los valores finales (de estado estable) al aumen
tar el tiem po Sin em bargo, con el m étodo explícito, esta solución se caracteriza por
oscilaciones num éricam ente inducidas, que son físicam ente im posibles Las oscilacio-
nes se vuelven inestables , lo que ocasiona que la solución difiera de las condiciones i
estado estable reales. Para evitar este tipo resultados erróneos, el valor estableen
de A t debe m antenerse por debajo de cierto lim ite, el cual depende de Av y otros pa
m etros del sistem a. Esta dependencia se denom ina criterio de estabilidad , y se obtien
m atem áticam ente o dem uestra partir de un argum ento term odinám ico (vease el pro
blem a 5.78). Para los problem as de interés en este texto, e / criterio se determina}
quiriendo que el coeficiente asociado con el nodo de Ínteres en el tiempo anterior v
mayor que o igual a cero En general, esto se hace reuniendo todos los términos
incluyen T& „ para obtener la form a del coeficiente. Este resultado sirve entonces pa
obtener una relación lím ite que incluya Fo, del cual se determ ina el máximo valor i
m isible de A t. Por ejem plo, con las ecuaciones 5.71 y 5 73 ya expresadas en la fon
que se desea, se sigue que el criterio de estabilidad para un nodo interior umdimenj
nal es (1 — 2 Fo) > 0 , o
(5J4
5‘
Para los valores establecidos de Aa y a . estos criterios sirven para determinar lír
superiores al valor de A t.
Las ecuaciones 5.71 y 5.73 tam bién se derivan al aplicar el método del balanceí
energía de la sección 4.4.3 a un volum en de control alrededor del nodo interior,
explicar cam bios en el alm acenam iento de energía térm ica, una forma general de',
ecuación de balance de energía se expresa com o
Con el interés de adoptar una m etodología congruente, de nuevo se supone que todo;
flujo de calor esta adentro del nodo.
Para ilustrar la aplicación de la ecuación 5.76, considere el nodo superficial)
sistem a unidim ensional que se m uestra en la figura 5.12. Para determ inar másprec
m ente las condiciones térm icas cerca de la superficie, a este nodo se le asigna i
pesor de la m itad del que tienen los nodos interiores. Al suponer transferencia:
convección desde un Huido contiguo y ninguna generación, se sigue de la ecua
5 76 que
kA Sx rg+l - re
hA(T„ r~g ) + — (
T ,.h
7o r,* t 2. 7b*
í í
I rfcalm '«*■
I
^cortv I — .li —J
i iI *?cond
2h A t 2 a At
7 T 1 = (T’oc - m - TQ +
Al reconocer que (2/í A f/pr Av) = 2(/z Av/A)(a At/Ax2) = 2 BiFo y reagrupar térm inos
que incluyen Tq. , se sigue que
h Ax
Bi = — — (5.78)
k
A l recordar el procedim iento para determ inar el criterio de estabilidad, requerim os que
el coeficiente para T q sea m ayor que o igual a cero, de aquí
1 — 2 Fo — 2BiFo ^ 0
D ado que la solución cn diferencias finitas com pleta exige la ecuación 5.73 para los
nodos interiores, así tam bién la ecuación 5.77 se requiere para el nodo superficial. La
ecuación 5.79 debe contrastarse con la ecuación 5.74 para determ inar cuál requisito es
el más riguroso. C om o Bi ^ 0, es evidente que el valor límite de Fo para la ecuación
5.79 es m enor que el de la ecuación 5.74. Por tanto, para asegurar la estabilidad en to-
dos los nodos hay que usar la ecuación 5.79 a fin de seleccionar el valor m áxim o per-
m isible de Fo , y de aquí At, para ser utilizados en los cálculos.
Las form as de la ecuación explícita cn diferencias finitas para varias geom etrías
com unes se presentan en la tabla 5.2. C ada ecuación se deriva aplicando elm étodo del
balance de energía a un volum en de control alrededor del nodo correspondiente. Con el
propósito de desarrollar confianza en su habilidad para aplicar este m étodo, intente ve-
rificar al m enos una de estas ecuaciones.
d e pa r t a m e n t o d e b ib l io t e c a
U n iv e r s id a d o iiiiu n B o lív a r - S e d e C LUoral
252 C apítulo 5 ■ Conducción en estado transitorio
—c
I
r- + . O CN
oo Cs ni On ON
lO »o — i un
4- c
un o
+ u
+I + O.E
—c he
£eE + .
c-aiiibiiiciud p¡mi unn mipcrflcie adiabática (o superficie de simetría), simplemente haga fíi igual
+
c +
P£4IC0 CN — —cI
[u + c# + -
“EL ++ h, i4
I —c
+ . — c ¡Sí —c
E y -+ c + + £ +-
P-E
t". + l «2 £
c c
■o _ GQ
o
-w GQ
CM n*o GQ
J. 1 CN TI 3
O) — c eNco
+
GQ
+
GQ
+ + + c + N-/
—
+ «.
-
w
<N 'o
+• «¡.E
k. h- £
CM
h,E 0.E
£ «2 h* «2 £
T}- rf CN 1 II rf
r
—H
Lh
-o un en un
r~ 00 oo 00
un m un un
-3 z -*N
4U-- «
y) VI VI VI VI
d>
X
w n> ■
*M ■
o I •« *2 GQ GQ
i—• .ts o +
+
/r*
a
—
‘c<u en CN
Qi
i' X
41 l* «2
C u r— eC r-
■o '2 '2
o '2 "G G
o4)
od
4> c
r- oy—
00
V
4 »
CN
OO
o
4) OO
4))
4>
s un, > m >
CC
wnI—1 e g un C
'w^ h.1 O w he W4)
cr O 4) •»N
Vi
« Q c r-<
GN c C GNQ c C
r/3
C O p.S C p.5
kV O he 4) + he 4)
c + nh. —v i-, 3r— c.
(R O /-s w1
O « £ e £ *C
- i 511 £ 4)
c
h.1 GQ
3 4)
c íkE OQ CN "Se o.S rfGQ X
4u) 0.E CN hf NICP .2 GQ .2
he 4>
-1» h- V, + + o1 C 3 + CN 1 ‘5 + 3
+ + /cT j; C c + rQ¡ tnUn c O C
k. t_ + rfk. 3 * k. 4) ° k. *5
+ rf O + C O * 1 1 rf a. Tj- 4T
4) c 1 Un Vi c 1 3Vi Vi
4) 1 4) 1
CN o.E
ti
>«C
— y
i■*
hsE C
c CN C he 'o C
■
a
U 4) 4) k. 4)
it o e*n + + C i2 + + C CN + O
3 3 •a ■a *3 •O
4> ii II II
w c O O — o
Z —
+£.5c- Z —
+[Link]. Z +p. cE- Z
•—t CN h. en he
u-
i
_ C.
1E i
1 1
+
h- -rH
r*- -H -
5 .9 ■ Métodos de direfencias finitas 253
E je m pl o 5 .7
Un elem ento com bustible de un reactor nuclear en la form a de pared plana de espesor
2L = 20 m m se enfría convectivam ente en am bas superficies, con h = 1100 W /m 2 • K
y Tos = 250°C. A potencia normal de operación, el calor se genera de modo uniform e
dentro del elem ento a una rapidez volum étrica de c¡j = 107 W /m 3. Si hay un cam bio en
la rapidez de generación, ocurrirá una desviación de las condiciones de estado estable
asociada con la operación norm al. C onsidere un cam bio súbito a tj2 = 2 X 107 W /m 3, y
use el m étodo explícito de diferencias finitas para determ inar la distribución de tem pe-
raturas del elem ento com bustible después de 1.5 s. Las propiedades térm icas del ele-
m ento com bustible son k = 30 W /m • K y a = 5 X 10 ~6 m 2/s.
S o l u c ió n
E sq u em a :
Elemento combustible
?! = l x 107 W/m3 = 250°C
¿l2 = 2 x 107 W/m3 h = 1100 W/m2 -K
a = 5 x 10' 6 m2/s Adiabática de
Fluido
* = 30 W/m • K simetría refrigerante
i i
i l
m —1 I m in i+ 1
• l i •
l l
ep k
^cood ^cond E airt íc o n d
A.v L
¿Ur=-
~1 Yo"
S u p o sicio n es:
1. C onducción unidim ensional cn x.
2. G eneración uniform e.
3. Propiedades constantes.
A nálisis: Se obtendrá una solución num érica con un increm ento espacial de A.v =
2 mm. C om o hay sim etría alrededor del plano m edio, la red nodal da seis tem peraturas
nodales desconocidas. Con el m étodo de balance de energía, ecuación 5.76, se deriva
una ecuación explícita en diferencias finitas para cualquier nodo interior m.
7-p+l _ y??
T P T V T P — T P
1m m
Ax
DEPARTAMENTO DE BIBLIOTECA
U n iv b rsid u d S im ó n liv o r S ede ' 1 kn rn '
254 C apítulo 5 ■ Conducción en estado transitorio
Esta ecuación sirve para el nodo 0, con Tpm- X = T^,+ 1, así com o para los nodos 1. 2,3
4. Al aplicar la conservación de la energía a un volumen de control alrededor del nodo:
T'¡ - Tp Ax Ax Tp5+1 -
hA(T^ - T%) + k A 7 -— + qA — = pA — c ~
Ax Ar
o
q ( lx ) 2 -,
T Y 1 = 2 Fo VI + B iT ^ +
2k
+ (1 - 2 Fo - 2B i F o m (2
C om o el criterio de estabilidad más restrictivo se asocia con la ecuación 2, selecciona-
m os Fo del requerim iento
Fo{ 1 + B i ) < \
De aquí, con
h Ax 1 100 W /m 2 • K (0.002 m)
Bi = = 0 0733
3 0 W /m • K
se sigue que
l o < 0.466
F ü (A x )2 0 .4 6 6 (2 X 10 ~ 3 m )2
At = < 0.373 s
a 5 X 10 6 m /s
5 X 10 -6 m 2/s(0 .3 s)
Fo = = 0.375
(2 X 10 - 3 m)
qL 10 7 W /m 3 X 0.01 m
7 , = 7? + 1— = 2 5 0 °C + — ............... , = 340.91 C
1100 W /m 2 ■ K
5 .9 ■ Métodos de direfencias finitas 255
se sigue que
Las tem peraturas calculadas para los puntos nodales de interés se muestran en el pri-
m er renglón de la tabla adjunta.
Con el uso de ecuaciones en diferencias finitas, las tem peraturas nodales se calcu-
lan de m anera consecutiva con un increm ento de 0.3 s hasta que se alcanza el tiempo
final deseado. Los resultados se ilustran en los renglones 2 a 6 de la tabla y se pueden
contrastar con la nueva condición de estado estable (renglón 7), que se obtuvo con las
ecuaciones 3.42 y 3.46 donde q = c¡2:
T e m p e ra tu ra s n o d ales ta b u la d a s
p / ( s) T0 Ti t 2 *3 T, Ts
0 0 357.58 356.91 354.91 351.58 346.91 340 91
1 0.3 358.08 357 41 355.41 352.08 347.41 341.41
2 0.6 358.58 357.91 355.91 352.58 .347.91 341.88
3 0.9 359.08 358.41 356.41 353.08 348.41 342.35
4 1.2 359.58 358.91 356.91 353.58 348.89 342.82
5 1.5 360.08 359.41 357.41 354.07 349.37 343 27
oc oc 465.15 463.82 459.82 453.15 443.82 431.82
C o m entarios: Es evidente que a 1 .5 s la pared está en las prim eras etapas del proce-
so transitorio y que se tendrían que hacer m uchos cálculos adicionales para alcanzar las
condiciones de estado estable con la solución en diferencias finitas. El tiem po de cálcu-
lo se reduce ligeram ente usando el increm ento de tiem po m áxim o perm isible (Ar =
0.373 s), pero con alguna pérdida de precisión. Con el interés de m axim izar la preci-
sión, debe reducirse el intervalo de tiem po hasta que los resultados calculados se hagan
independientes de reducciones posteriores de A t.
Al extender la solución en diferencias finitas, es posible determ inar el tiem po que
se requiere para alcanzar la nueva condición de estado estable, con historias de tem pe-
raturas calculadas para los nodos del plano m edio (0) y de superficie (5) que tienen las
siguientes formas:
/(s)
DEPARTAMENTO DE BIBLIOTECA
U n iv e rs id a d S im ó n B o lív a r - S e d e ci*« L itoral
256 C apítulo 5 ■ Cttndurrión en estado transitoriu
método implícito
a At (S x )2
77?+ 1 i 'TP + 1 ' ) rp p + 1
* m.n+ I * m.n—1 m,n
+
Para cualquier nodo interior de la figura 5.12, tam bién se m uestra que
E je m p lo 5 . 8
U na placa gruesa de cobre que inicialm entc está a tem peratura uniform e de 20°C
se expone de súbito a radiación en una superficie de m odo que el flujo neto de calor se
m antiene a un valor constante de 3 X 105 W /m 2 Con las técnicas en diferencias finitas
explícita e im plícita y un increm ento espacial de Aa = 75 mm, determ ine la tem peratu-
ra en la superficie irradiada y en un punto interior que esté a 150 mm de la superficie
después de transcurridos 2 m inutos. C om pare los resultados con los que se obtienen de
una solución analítica apropiada.
SüLl CIÓN
Se co n o ce: Placa gruesa de cobre, inicialm ente a una tem peratura uniform e, que se
som ete a un flujo neto constante de calor en una superficie.
E n con trar:
1. C on el m étodo explícito de diferencias finitas, determ ine las tem peraturas en la
superficie y a 150 m m de la superficie después de transcurrido un tiem po de 2
m inutos.
2. R epita los cálculos con el m étodo im plícito de diferencias finitas
3. D eterm ine las m ism as tem peraturas de form a analítica.
E squem a:
S u p o sicio n es:
1. Conducción unidim ensional en x.
2. La placa gruesa se aproxim a com o un m edio sem iinfinito con flujo de calor super-
ficial constante.
3. Propiedades constantes.
P ro p ied a d es: Tabla A .l, cobre (300 K): k = 401 W /m • K, a = 117 X 10 " 6 m 2/s.
A nálisis:
1. Una form a explícita de la ecuación en diferencias finitas para el nodo superficial
se obtiene aplicando un balance de energía a un volum en de control alrededor
del nodo.
- 7'(¡ A* ’ 1 - r„ 1
q„A + k A — ^ r = PA — -c I
O
/ q"a A* \
1%*' = 2 F o l — - — + TPA + (1 - 2Fo)Tp
La ecuación en diferencias finitas para cualquier nodo interior está dada por la
ecuación 5 73. Los nodos superficial e interior están regidos por el criterio de esta-
bilidad
Fo< \
(Ax) 1 (0.075 m )2
Ar = Fo = 24 s
a 2 117 X 10 ~ 6 m 2/s
Con
7 ? +l = 56 .1 °C + Tp y
j'p+i
1 m
_ rm+ 1 -l- 7T
P
1 m—I
p *(s) T0 Ti T2 t 3
0 0 20 20 20 20
20
1 24 76 1 20 20 20 20
2 48 76 1 48.1 20 20 20
3 72 104.2 48 1 34 1 20 20
4 96 104 2 69.1 34.1 21A 20
5 12 0 125 3 69.1 48 1 27 1 20
una en la que la tem peratura cam bia de forma continua con el tiempo. La deform a-
ción se elim ina y la precisión de los cálculos se m ejora reduciendo el valor de Fo.
Para determ inar el punto al que es posible m ejorar la precisión al reducir Fo,
rehagam os los cálculos para Fo = 4 (Ar = 12 s). Las ecuaciones en diferencias fi-
nitas son entonces de la form a
7%+' = h(56.\°C + Tp) + \T P
0
T p; x = \ a pm+, + 7 ^ , ) + ^ ,
y los resultados de los cálculos se tabulan com o sigue:
p f(s) Ti t 2 Ti t 4 Ts n T-, Ts
0 0 20 20 20 20 20 20 20 20 20
1 12 48.1 20 20 20 20 20 20 20 20
2 24 62.1 27.0 20 20 20 20 20 20 20
3 36 72.6 34.0 2 1.8 20 20 20 20 20 20
4 48 81.4 40.6 24.4 20.4 20 20 20 20 20
5 60 89 0 46.7 27.5 21.3 20 1 20 20 20 20
6 72 95.9 52.5 30.7 22.6 20.4 20.0 20 20 20
7 84 102.3 57.9 34.1 24.1 20.8 2 0 .1 20.0 20 20
8 96 108.1 63.1 37.6 25.8 21.5 20.3 20.0 20.0 20
9 108 113.7 68.0 41.0 27.6 2 2 .2 20.5 2 0 .1 20.0 20.0
10 12 0 118.9 72.6 44.4 29.6 23.2 20.8 202 20.0 20 0
C om o tratam os con un solido sem iinfinito, el núm ero de nodos es, en princi-
pio, infinito. En la practica, el núm ero está lim itado a los nodos que están afecta-
dos por el cam bio en la condición de frontera para el periodo de tiem po de interés.
De los resultados del m étodo explícito, es evidente que se pueden elegir con segu-
ridad nueve nodos correspondientes a 7’0, Tj,..., T8. Suponem os que en t = 120 s
no hay cam bio en ÜTg.
Tenem os ahora un conjunto de nueve ecuaciones que deben resolverse simul-
táneam ente para cada increm ento de tiem po. Con el m étodo de inversión de matri-
ces, expresam os las ecuaciones en la form a |AJ[TJ = [CJ, donde
2 -1 0 0 0 0 0 0 0
-1 4 -1 0 0 0 0 0 0
0 -1 4 -1 0 0 0 0 0
0 0 -1 4 -1 0 0 0 0
[A] 0 0 0 -1 4 -1 0 0 0
0 0 0 0 -1 4 -1 0 0
0 0 0 0 0 -1 4 -1 0
0 0 0 0 0 0 -1 4 -1
0 0 0 0 0 0 0 -1 4
56.1 + 7£
2 T px
2 r p2
2T%
2r j
\C] = 2T p
5
2n
2 T pn
+ TP + I
- 8 0
O bserve que los valores num éricos para los com ponentes de [CJ se determinan 4
los valores anteriores de las tem peraturas nodales. A dvierta tam bién cómo laec
ción en diferencias finitas para el nodo 8 aparece en las m atrices [A] y [C].
Se puede arm ar una tabla de tem peraturas nodales, que comience con el p¿
m er renglón {p = 0 ) correspondiente a la condición inicial que se establece,
obtener tem peraturas nodales de los tiem pos siguientes, debe encontrarse prir
la inversa de la m atriz coeficiente [A]-1. En cada tiem po p + 1, se multiplica)
tonces por el vector colum na [C], que se evalúa en p , para obtener las temperatu
T¡¡+1, T T T¡¡+ . Por ejem plo, al m ultiplicar [A] - 1 por el vector colui
correspondiente a p = 0 ,
76.1 “
40
40
40
[C]p = o = 40
40
40
40
L 60 J
p / (s ) T0 Ti Ti Ty t 4 Ts Tff t 7 T8
0 0 20 .0 20.0 20 0 20 0 20.0 20 0 20 0 20.0 20.0
1 24 52.4 28.7 22.3 20.6 20 .2 20.0 20.0 20.0 20.0
2 48 74.0 39 5 26 6 2 2 .1 20.7 20 2 2 0 .1 20.0 20.0
3 72 90.2 50.3 32.0 24.4 21 6 20.6 20 .2 20 .1 20.0
4 96 103.4 60 5 38 0 27.4 22.9 21.1 20 4 20 .2 20 1
5 12 0 114.7 70.0 44.2 30.9 24.7 21.9 20.8 20.3 2 0 .1
2 X 3 X 10 5 W /m 2
7 (0 , 1 2 0 s ) - 2 0 °C = — ■■■— — (117 X 10 - 6 m 2/ s X 120 s /t t ) 2
401 W /m • K
2 X 3 X 10 5 W /m 2
7X0.15 m , 120 s) — 2 0 °C = 4oT W /m • K
(0.15 m )2 1 3 X 10 5 W /m 2 X 0.15 m
X exp
4 X 117 X 10 -6 m 2/s X 120 s 401 W /m • K
0.15 m
x 1 — e rf = 4 5 .4 °C
2V 117 X 10 -6 m 2/s X 120 s
C o m en ta rio s:
1 . C om parando los resultados exactos con los que se obtienen de las tres soluciones
aproxim adas, es claro que el m étodo explícito con Fo = \ proporciona prediccio-
nes más precisas.
d ep a r t a m en t o DE BIBLIOTECA
U n iv e rs id a d S im ó n R nllvs- Sede del Litor»'
202 C apítulo 5 ■ Conducción en estado transitorio
N ote que, para t = 120 s, la suposición de un m edio sem iinfinito seguiría sie
válida si el esp eso r de la placa excede aproxim adam ente 500 mm.
3. O bserve que la m atn z coeficiente fA] es tridiagonul. E s decir, todos los elemu
son cero excepto los que están en la diagonal principal o en cualquier lado de<
Las m atrices tridiagonalcs están asociadas con problem as de conducción
m ensiona!.
4. U na condición de calentam iento radiativo m ás general sería aquella en la qir|
superficie se expone súbitam ente a los alrededores a una tem peratura elevada!
(problem a 5.91). La transferencia neta por radiación a la superficie se calcula^
toncos a partir de la ecuación 1.7. A l perm itir la transferencia de calor por con
ción a la superficie, la aplicación de la conservación de la energía al nodo|
superficie da una ecuación explícita en diferencias finitas de la forma
F '-F g Ax F g+I - F£
e o t F ^ - (Fg)4l + h(Tx - Fg) + k = P
Ax Ai
A plicar esta ecuación en diferencias finitas en una solución num érica es cor
do por el hecho de que es no lineal Sin em bargo, la ecuación se linealiia mi
la introducción del coeficiente de transferencia de calor hr definido por laca
1.9, y la ecuación en diferencias finita^ es
Tp _ Tp A x- F g + I - Fg
~ Fg) + h{T„ — Fg) + k
Ax = P T ’ C Ai
5.10
Resumen
B ibliografía
1. Carslaw, H. S. y J. C. Jaeger, Conduction q f Heat in 5. Hcisler, M. P , “Tem perature C harts for Induction
SolidSs 2a. cd., O xford U m versity Press, Londres, and C onstant Tem perature Heating” , Trans. ASME ,
1959. 6 9 , 227-236, 1947.
2. Schneider, P. J., Conduction Heat Transfer, Addi- 6 . G róber, H., S , Erk y U. G rigull, fundam entáis o f
son-Weslcy, Rcading. M A, 1955. Heat Transfer, M cG raw -H ill, N ueva York, 1961.
3 Kakac, S y Y. Yener, Heat Conduction , H em ispherc 7. Langston, L. S., ‘‘Heat Transfer from M ultidim en-
Publishing, Nueva York, 1985. sional O bjets Using O ne-D im ensional Solutions for
4 Poulikakos, D., Conduction Heat Transfer, Prenti- H eat L oss” , Int. J. Heat Mass Transfer, 25, 149-
ce-Hall, Englewood Cliffs, NJ, 1994. 150, 1982.
Problem as
X = L
(al Dibuje y acote, en coordenadas T -.\\ las distribu- (a) Dibuje y acote, en coordenadas T -x . las distribu-
ciones de temperaturas inicial, de estado estable ciones de temperaturas: inicial, de estado estable
y en dos tiempos intermedios. y en dos tiempos intermedios.
(b) Trace e flujo de calor en la superficie externa diante la observación de la historia de tempenum»
q (L , /) com o Iunción del tiempo de una estera fabricada con cobre puro La esfera, yJ
tiene 12 7 mm de diám etro está a 66 C antes dec ln
5.3 Un hom o de microondas opera sobre el principio de
caria en un flujo de aire que tiene una temperatura 4
que la aplicación de un campo de alta frecuencia oca
27 C Un term opar en la superficie externa de laeifc.
siona que oscilen as moléculas eléctricamente polari
ra indica 55 C 69 s después de que se inserta laesf®
zudas de los alimentos El efecto neto es la gene i u i m
en el flujo de a re Suponga y después justifique fd
casi uniforme de energía térmica dentro de los ahmen
la estera se com porta com o un objeto espacial isom
tos Considere el proceso de cocinar un trozo de carne
mico, y calcule el coeficiente de translercncia decaí#
de 2L de espesor en un hom o de microondas y compá
relo con cocinar en un hom o convencional, donde cada 5.7 Una esfera solida de acero (AISI 1010), d 300 mnu
lad se calienta por radía íón. En cada caso la carne diámetro, se recubre con una capa de material diel
se calentara de 0 C a una temperatura n mima de 90°C. trico de 2 mm de espesor y una conducto dad térqti
Base su comparación en una gra ca de la distribución de 0.04 W /n • K La esfera recubierta está miciali
de temperaturas en tiempos seleccionados para cada te a una temperatura un forme de 500 C y de proni ,
uno de los procesos de cocn ado En particular, consi templa en un baño de aceite para el que Tx = 100°Ci
dere el tiempo q al que se inicia el calentamiento, un h = 3300 W m • K Estime el tiempo que se rec
tiempo /, durante el proceso de calentamiento, el tiem p a n que 1 1 esfera recubierta alcance 140 C. Stigtid
po t que corresponde a la conc usión del calcntamicn cía: Deje de lado el efecto de almacenainiei
to y un tiempo /3 completamente dcntio del subsi ener ía en el m atcnal dieléctrico, puesto que su
guíente proceso de enfriamiento. citancia térmica (pcV) es pequeña com arada a»
de la esfera de acero
5.4 Una placa de espesor 2L, área superficial A s, masa V/
y calor especí ico i p, inicialmente a una tem peratura 5.8 Una bala esférica de plomo de 6 mm de diámetn*i
uniforme T , se calienta de pronto en am bas superfi mueve aproxim adam ente a Mach 3 La onda de i
cíes med ante un proceso de convección (T , h) du que rem ítante cal cnta el aire alrededor de la
rante un periodo t0, después del cual la placa se ai la. 700 K, y el coeficiente de convección promedio i
Supon 7a que a tem peran ra del pk no medio no al la trans crcncia de calor ntre el aire y la bala es3
canza 7 dentro de este periodo W n r • K Si la ba a sale del barril a 300 K y el |
a) Suponiendo B > 1 para el proceso de calenta po de vuelo es 0.4 s, ¿ cuál es la temperatura en la
miento, dibuje y acote, en coordenadas T - \. las si perlicie en el momento del impacto9
guientcs distribuciones de tem peraturas inicial, 5.9 Unos ejes de maqu nana de acero al carbón (
de estado estable (t —> <*), T{\, tfí) y en dos tiem 1010) de 0.1 m de diámetro se tratan con calor<
pos intermedios entre t = t y t —> *>. hom o calentado por gas cuyos gases están a 12001
(b) Trace y acote, ei coordenadas T /, las distribu y proporc onan un coeficiente de convección de|
ciones de ten peratura del plano medio y de la su- W m • K Si los ejes ent an en el homo a
perficie expuesta ¿cuánto tiempo deben permanecer en el homo i
(c) Repita las partes (a) y (b). suponga B, < 1 para la canzar una temperatura en la linea central de 8001
placa.
(d) Derive una expresión para la temperatura de estado 5.10 Una unidad de almacci am cnto de eneigía t
estable T a , oo) = T , y deje los resultados en térmi consi te en un canal rectangular largo, que está|
nos de los parámetros de la placa (M, c , condicio ai lado en la superficie extem a y encierra capas i
nes térmicas (7~, Tx , h), la temperatura de la super nadas del material de almacenamiento y rejilla»,!
licie 7 (L, t) y el tiempo de calentamiento t , el flujo
Cada capa del material de alm acenam iento es una tuerza electrostática y el flujo de gas helio. Después
plancha de aluminio de ancho W = 0.05 m, que está a de 15 s, se determ ina que la temperatura de la chapa
una temperatura inicial de 25°C Considere condicio- es 33°C ¿Cuál es la resistencia de contacto térmico R'¡
nes en las que la unidad de alm acenam iento se carga c (m ■ K/W) entre la chapa y el portaherramientas?
con el paso de un gas caliente a través de las rejillas, ¿El valor de R'¡ c aumentará, dism inuirá o perm anece-
suponiendo que la tem peratura del gas y el coeficiente rá igual si se usa aire, en lugar de helio, como gas de
de convección tienen valores constantes de T« = purga?
600°C y h — 100 W /m 2 • K a lo largo del canal
5.13 Del ultimo grupo de rodillos cn un tren de laminación
Cuanto tiempo se tardará cn alcanzar 7 5 $ del alm a-
en callente emergen tiras de acero y se enfrían (por
cenamiento máximo posible de energía? ¿Cuál es la
ambas superficies) mediante transferencia de calor
temperatura del alum inio cn ese m om ento9
por convección y radiación al aire ambiente y los al-
5.11 Un resorte de hojas cuyas dimensiones son 32 mm por rededores, respectivamente, donde T« = Talr = 300 K
10 mm por 1 1 m se rocía con un recubrimiento antico-
rrosivo delgado, que se trata con calor suspendiendo el
resolte de forma vertical en la dirección de su longitud r * = 300K
Aire
y pasándolo a través de un homo transportador que
= 300 K
mantiene el aire a una temperatura de 175°C Se han
obtenido recubrimientos satisfactorios cn resortes, mi-
cialmente a 25°C, con un tiempo de permanencia en el
homo de 35 min. El proveedor del recubrimiento espe- Aire
cifica que el recubrimiento debe tratarse durante 10 min r . = 300 K
por amba de una temperatura de 140°C. ¿Cuánto tiem- J
po permanecerá en el homo un resorte con dimensiones
de 76 mm por 35 mm por 1.6 m a fin de tratar térmica- El espesor de las tiras es 8 = 5 mm, y su densidad,
mente el recubrimiento de manera apropiada9 Las pro- calor específico, conductividad térm ica y emisividad
piedades termofísicas del material de resorte son p = son p = 7900 kg/m 3, cp = 640 J/kg ■ K, k = 30 W/m •
8131 kg/m3. cp = 473 J/kg • K, y k = 42 W/m • K. K y e = 0.7, respectivamente.
5.12 Una herramienta que se utiliza para fabricar dispositi- (a) Suponiendo un coeficiente de convección unifor-
vos semiconductores consiste en un portaherram icnta me de h = 25 W /m 2 • K, determine el tiempo que
(disco cilindrico metálico grueso) en el que un brazo se requiere para enfriar la tira de una temperatura
robótico coloca una chapa de silicio muy delgada (p inicial de 940°C a 540°C, punto en el que se pue-
= 2700 kg/m3. c = 875 J/kg • K, k = 177 W /m • K) de enrollar para su embarque. Si la tira se mueve
Una vez en su posición, se energiza un cam po eléctri- a 10 m/s, ¿cuán larga debe ser la sección de en-
co en el portaherramientas, lo que crea una fuerza friamiento? ¿Qué concluye usted acerca de la
electrostática que mantiene la chapa firmemente fija efectividad de este método de enfriamiento?
al portaherramientas. Para asegurar una resistencia de
(b) Determ ine el tiempo que se requiere para enfriar
contacto reproducible entre el portaherram ientas y en-
la tira, primero suponiendo una transferencia de
tre ciclo y ciclo, se introduce gas helio presurizado en
calor sólo por radiación, y después una sólo por
el centro del portaherramientas y fluye (muy lenta-
convección. Para cada uno de los tres casos (con-
mente) de forma radial hacia afuera entre las imper-
vección y radiación, sólo radiación y sólo
fecciones de la región de interfaz
Chapa Tjt). convección), elabore una gráfica de la temperatura de
TJ0) = Tw, = 100X la tira como función del tiempo cn el rango 54f)°C
w = 0.758 mm < T ^ 940°C. En este rango, también trace el
coeficiente de transferencia de calor por radia-
t ción, h r%como función del tiempo.
Gas helio 5.14 La pared plana de un hom o se fabrica de acero al
de purga carbón simple (A = 60 W /m • K, p = 7850 kg/m , c
Región de la interfaz Porta
= 430 J/kg • K) y tiene un espesor L = 10 mm. Para
muy exagerada herramientas. Xt = 23°C
protegerla de los efectos corrosivos de los gases de
Se rea iza un experimento en condiciones en las que combustión del horno, una superficie de la pared se
U chapa, micialmentc a una tem peratura uniforme Tw cuore con una película delgada de cerámica que. para
- 100°C. se coloca de pronto en el portaherram ien- un area superficial unitaria, tiene una resistencia térmica
tas. que está a una temperatura uniforme y constante de R'¡ f = 0.01 m • K/W. La superficie opuesta esta
T - 23°C. Con la chapa en su lugar, se aplican la bien aislada de los alrededores.
DEPARTAMENTO DE BIBLIOTECA
U n lv e rs d a d u u n « jr - S o d ii. jf
266 C apítulo 5 ■ Conducción en estado transitorio
Película las tem peraturas, 7\/), para los dos casos y sel
i
cerámica, Gases del horno, cionc el sistema de respaldo más adecuado.
p. c, k. T¡
5.16 Una tira de acero de espesor 5 = 1 2 mm se recu
Acero al haciéndola pasar a través de un horno grande cu
carbón paredes se mantienen a una temperatura TH
T„,h
corresponde a la de los gases de combustión que fluv
a través del horno (7* = 7V„). La tira, cuya densi
T,
calor específico, conductividad térmica y emisivr
son p = 7900 kg/m , cp = 640 J/kg • K. k = 30 Y'
U ' K y e = 0.7, respectivamente, se calentará de 3()05C
600°C.
Al poner en funcionamiento el horno, la pared está a
1
una tem peratura inicial de /', = 300 K, y los gases de
combustión entran en el hom o a Tx = 1300 K. con lo Paredes del horno. Tw
que proporcionan un coeficiente de convección h =
25 W /m 2 • K en la película cerámica. Suponiendo que ,***Gases de *\
la película tiene una resistencia térmica interna insig-
nificante, ¿cuánto tiempo tardará la superficie interior
del acero en alcanzar una tem peratura de Ts ¡ = 1200
K? ¿Cuál es la temperatura Ts t> de la superficie ex-
puesta de la película cerám ica en esc momento? Gases de
5.15 hn un proceso industrial que requiere altas corrientes combustión
5.18 Considere el sistema del problema 5 1 donde la tem 5.21 Una esfera de metal de diámetro D, que está a tem pe-
peratura de la placa es isotérm ica espacial durante el ratura uniforme T¡. se quita súbitamente de un horno y
proceso transitorio. se cuelga de un alambre fino en un cuarto amplio con
aire a una tem peratura uniforme 7 x y las paredes que
(a) Obtenga una expresión para la tem peratura de la
lo rodean a una temperatura T a l r -
placa como función del tiem po T(t) en términos
de q"), 7*, /í, L, y las propiedades p y i de la placa. (a) Sin tomar en cuenta la transferencia de calor por
radiación, obtenga una expresión del tiempo que se
(b) Determine la constante térmica de tiempo y la
requiere para enfriar la esfera a alguna temperatura T
temperatura de estado estable para una placa de
12 mm de espesor de cobre puro cuando 7 X = (b) Sin tomar en cuenta la transferencia de calor por
27°C. h = 50 VV/m2 • K, y q"K = 5000 W /m 2. Es- convección, obtenga una expresión del tiempo que
time el tiempo que se requiere para lograr condi- se requiere para enfriar la esfera a la temperatura T
ciones de estado estable.
(c) ¿Cómo determinaría el tiempo que se requiere para
(c) Para las condiciones de la parte (b). así com o pa- que la esfera se enfríe a la temperatura T si la con-
ra h = KM) y 200 W/m • K. calcule y trace las vección y radiación son del mismo orden de mag-
historias de tem peratura correspondientes de la nitud?
placa para 0 < r < 2500 s.
(d) Considere una esfera de aluminio anodi/ado (e =
j.19 Un dispositivo electrónico, com o un transistor de po- 0.75) de 50 mm de diám etro, que está a una tem -
tencia montado sobre un disipador de calor con aletas, peratura inicial de 7, = 800 K. Tanto el aire como
se modela como un objeto cspacialm cnte isotérmico los alrededores están a 300 K, y el coeficiente de
con generación interna de calor y una resistencia de convección es 10 W /m 2 • K. Para las condiciones
convección externa. de las partes (a), (b) y (c). determine el tiempo
que se requiere para que la esfera se enfríe a 400
(a) Considere un sistema de masa M, calor especílico
K Elabore una gráfica de las series históricas de
c y área superficial A s, que inicialmente está en
temperatura correspondientes. Repita los cálculos
equilibrio con el medio a Tx. De súbito se energi-
para una esfera de aluminio pulido (e = 0.1)
/.a el dispositivo electrónico, de modo que ocurre
una generación de calor constante t (W). M ues- 5.22 A medida que las estaciones espaciales permanentes
tre que la respuesta de tem peratura del dispositi- aumentan de tamaño, hay un incremento concom itan-
vo es te en la cantidad de potencia eléctrica que disipan Pa-
ra prevenir las temperaturas de los compartimientos
6 l t \
de la estación de modo que no excedan los límites es-
tablecidos, es necesario transferir el calor disipado al
espacio. Un nuevo esquem a de rechazo de calor que
donde 0 = T - y r ( x ) es la tem peratura de se propone para este propósito se denomina Radiador
estado estable que corresponde a t ^ í), = 7, de gotitas líquidas (LDR. Liquid Droplet Radiator).
- 7(°°): T = tem peratura inicial del dispositivo; Prim ero se transfiere el calor a un aceite de alto vacío,
R = resistencia térmica \/h A s: y C = resistencia que después se inyecta al espacio exterior como un
térmica interna M c. flujo de pequeñas gotas. Se perm ite que el flujo atra-
(b) Un dispositivo electrónico, que genera 60 W de viese una distancia L. en la que se enfría por radiación
calor, se monta en un disipador de calor de alum i- de energía al espacio exterior a temperatura del cero
nio que pesa 0 31 kg y alcanza una temperatura de absoluto. Las gotitas entonces se reúnen y se devuel-
100°C en aire ambiente a 20°C en condiciones ven a la estación espacial.
de estado estable. Si el dispositivo está inicialmen-
te a 20°C, ¿qué temperatura alcanzará 5 min des- i
r
Espacio exterior
pués de que se conecta la potencia? Inyector Colector
Tíllr= OK
de gotitas de gotitas
520 Antes de ser inyectado en un horno, se precalienta car-
bón pulverizado haciéndolo pasar a través de un tubo %
cilindrico cuja superficie se mantiene a / aIr =
100()°C Si los granos se aproximan como esferas de 1
mnt de diámetro y se puede suponer que se calientan
por transferencia de radiación de la superficie del tubo,
¿cuán largo debe ser el tubo para calentar el carbón D
que entra a 25°C a una temperatura de 600°C? ¿Se
justifica el uso de la resistencia interna despreciable? Retorno de aceite trio
DEPARTAMENTO DE BIBLIOTECA
ünivó.8ik.üü óiiiioii bolívar - Sedo «. .urafl
268 C apítulo 5 ■ Conducción en estado transitorio
Considere condiciones en que las gotitas con em isivi- (a) Sin tom ar en cuenta la radiación, obtenga unaet*
dad e = 0.95 y diámetro D = 0.5 mm se inyectan a presión del tiempo de vuelo, r,_y que se requien
una temperatura de T¡ = 500 K y una velocidad V = para calentar una partícula desde su tempraaturai
0.1 m/s. Las propiedades del aceite son p = 885 inicial T a su punto de fusión 7 ^ y, una vez enel
kg/m 3, c = 1900 J/kg • K. y k = 0.145 W /m • K. Su- punto de fusión, para que la partícula expcrime¿-
poniendo que cada gota radia a la profundidad del es- tc la fusión completa. Evalúe t, para T¡ = w |
pacio a 7a|r = 0 K. determine la distancia L que se y las condiciones de calentamiento estab e ¡(fe
requiere para que las gotitas impacten al colector a
(b) Suponiendo que la alúmina tiene una emisivid^
una temperatura final Tj = 300 K. ¿Cuál es la cantidad
de £p = 0.4 y que las partículas intercambian
de energía térm ica rechazada por cada gotita?
diación con los alrededores a 7^,1,. = 300 K.
5.23 A menudo se utilizan procesos de recubrimiento por lúe la validez de dejar de lado la radiación.
pulverizado de plasma para proporcionar protección su- 5.24 Unas varillas metálicas largas de sección trunsvaji
perficial a materiales expuestos a medios hostiles que circular se tratan por calentamiento haciendo pasar ug
inducen degradación a través de factores como uso. co- corriente eléctrica por ellas para proporcionar g
rrosión o falla térmica absoluta. Los recubrimientos ce- ción volumétrica uniforme a una rapidez q (W/m’).
rámicos se usan normalmente con este propósito. varillas son de diámetro D y se colocan en un c
Mediante la inyección de polvo cerámico a través de la grande cuyas paredes se mantienen a la misma te
boquilla (ánodo) de un soplete de plasma, las partículas ratura 7* que el aire encerrado. La convección de la
se alinean por el flujo de plasma, dentro del cual se ace- perficie de las varillas al aire se caracteriza
leran y calientan. coeficiente h.
Recubrimiento
cerámico
Sustrato
.26 Considere las condiciones del problem a 5.25. Además ¿Cuánto tiempo tardará la segunda pared en alcanzar
*J\
de tratar la transferencia de calor por convección di- 28.5°C en la posición x = L 2 ? Use como base del
rectamente del chip al fluido refrigerante, un análisis análisis la dependencia funcional adimcnsional para
más realista explicaría la transferencia indirecta del la distr bución de temperaturas transitoria que se ex-
chip al sustrato y luego del sustrato al fluido refrige- presa en la ecuación 5.38.
rante. La resistencia térm ica total asociada con esta
5.29 Con referencia a la herram ienta para procesar sem i-
ruta indirecta incluye contribuciones debidas a la in-
conductores del problem a 5.12, se desea en algún
terfaz chip-sustrato (una resistencia de contacto), con-
m omento del ciclo de fabricación enfriar el portahe-
ducción unidimensional en el sustrato y convección
rram ienta, que está fabricado con aleación de alum i-
de la superficie del sustrato al fluido refrigerante. Si
nio 2024. El esquem a de enfriam iento que se
esta resistencia térm ica total es R, = 200 K/W, ¿cuál
propone pasa aire a 15°C entre la cabeza del sum inis-
es la temperatura de estado estable del chip Tj'l Des-
tro de aire y la superficie del portaherramientas.
pués de la activación del chip, ¿cuánto tiem po lc toma
llegar a 1°C de esta tem peratura? I Suministro de
» aire , 2Q°C
l a k T¡ Tx h
Casi» (m) (n rte ) (YV/m • K ) (°C ) (°C ) (YV/in2 • K)
5.31 Una pared plana unidimensional con espesor de 0.1 m rizados por una tem peratura de 2300 K y un c
inicialmcnte a una tem peratura uniforme de 250°C se cíente de convección de 5000 W /m2 • K Para au
sumerge de pronto en un baño de aceite a 30°C. Su- tar la duración de la operación del motor, se pro
poniendo que el coeficiente de transferencia de calor que se aplique un recubrimiento de barrera térmica
por convección para la pared cn el baño es 500 W /m 2 = 10 W/m • K, a = 6 X 10~6 m 2, s) a la superficie
• K, calcule la temperatura de la superficie de la pared terior de la tobera.
9 nuil después de la inmersión Las propiedades de la (a) Si el recubrimiento cerámico tiene 10 mmde
pared son k — 50 W/m • K, p = 7835 kg/m3 y c = 465 pesor y una tem peratura inicial de 300 K. oble
J/kg • K. una estim ación conservadora de la máxima d
5.32 Considere la unidad de alm acenam iento de energía cion permisible de la operación del motor El
térmica del problema 5 10. pero con un material de dio de la tobera es mucho mayor que el es-
manipostería de p = 1900 kg/m 3, c = 800 J/kg • k y com binado de la pared y el recubrimiento
k = 0.70 W/m ■ K usado en lugar del aluminio. (b) Calcule y trace las temperaturas de la supe
¿Cuánto tiempo tomará alcanzar 75% del máximo a l- interna y externa del recubrimiento como fir
macenamiento posible de energía? ¿Cuáles son las dcl tiem po para 0 < / < 150 s. Repita los c
temperaturas máxima y mínima de la m anipostería en los para un espesor de 40 mm de recubrimie
este tiempo?
5.35 En un proceso de templado, una placa de vidri
5.33 La pared de la tobera de un cohete tiene un espesor L inicialmcnte está a una temperatura uniforme T
= 25 mm y esta fabricada de una alta aleación de ace- enfría mediante la reducción súbita de la tempe
ro para la que p = 8000 kg/m , i = 500 J/kg • K, y k de ambas superficies a Ty La placa tiene 20 nrn
= 25 W /m • K. Durante una prueba de encendido, la espesor, y el vidrio tiene una difusividad térmica dr
pared esta inicialmcnte a T, = 25°C y la superficie in- X 10~" m2/s.
terna se expone a los gases calientes de combustión (a) ¿Cuánto tiem po pasará para que la tempera
para los que h = 500 W /m 2 • K y T x = 1750°C. La del plano medio alcance 50% de su red
superficie externa esta bien aislada máxima posible de temperatura?
(b) Si (77 — 77) = 300°C, ¿cuál es el gradiente
tem peratura máximo en el vidrio en el tiempo
Pared de la tobera terior?
5.36 La resistencia y estabilidad de neumáticos se au
calentando ambos lados del hule (k — 0.14 Wm
a — 6.35 X 10~sm2/s) en una cámara de vaporp
que 77c = 200°C. En el proceso de calentamlento,
pared de hule de 20 mm de espesor (que se s
deshebrado) se lleva de una temperatura ¡ni
25°C a una tem peratura del plano medio de 150;C,
(a) Si la pared debe mantenerse al menos a 100°C por (a) Si el flujo de vapor sobre las superficies Je
debajo de su punto de fusión r ,)t = I600°C, ¿cuál neumáticos mantiene un coeficiente dt ::
es el tiempo de encendido máximo permisible tf? El ción h = 200 W/m • K, ¿cuánto tiempo i
diámetro de la tobera es mucho mayor que su espe- en alcanzar la temperatura del plano medio
sor. se desea?
(b) Para aum entar /y. se considera cam biar el espesor (b) A fin de acelerar el proceso de calentanu
de la pared L ¿Se debe aum entar o dism inuir L? recom ienda que el flujo de vapor sea sufi
¿Por qué? Para espesores de la pared de L = 10, mente vigoroso com o para mantener las si
25 y 50 mm. calcule y dibuje las historias de tem - cies de los neumáticos a 200°C a trav.
peratura de las superficies interna y externa cn el proceso Calcule y trace las temperaturas del
periodo 0 < r < 600 s. El valor de tf se aum enta- no medio y de la superficie para este caso,
ría seleccionando un material con diferentes pro- mo para las condiciones de la parte (a).
piedades termofísicas. ¿Hay que elegir materiales 5.37 Unas tarjetas de circuitos de fibra de vidrio r
de valores de p. c y k mayores o m enores9 con epóxico y recubiertas de cobre se tratan nr
5.34 Durante la operación transitoria, la tobera de acero el calentam iento de una pila de ellas a alta
del motor de un cohete no debe exceder una tem pera- como se muestra en el dibujo. El propósito de la
tura de operación máxima permisible de 1500 K ración de calentam iento-prensado es curar el
cuando se expone a los gases de combustión caracte- que une las hojas de fibra de vidrio, impartiendo
■ Problemas 271
de/ a las tarjetas. La pila, denom inada libro, se com - se comporta como una capa aislante en la superficie
pone de 10 tarjetas y 11 placas ele prensado, que evi- externa, ( quc coeficiente de convección interior per-
tan que el epóxico (luya entre las tarjetas e imparten mitiría a la superficie exterior alcanzar 0°C en 60 s?
un acabado suave a las tarjetas curadas. A fin de llevar Las propiedades termofísicas del parabrisas son p =
a cabo un análisis térmico simplificado, es razonable 2200 kg/m3, cp = 830 J/kg • K. y k = 12 W/m • K
iprox nar el libro como si tuviera una conductividad
ttimica electiva (k) y una capacitancia térmica efecti-
va pe,,). Calcule las propiedades efectivas si cada una C o iu liie e ió ti iiiii<liineii>ional:
de las tarjetas y placas tiene un espesor de 2.36 mm y c ilin d r o la r g o
las siguientes propiedades termofísicas tarjeta (b) ph =
1000 kg/m . cph = 1500 J/kg • K , k h = 0 30 W/m • K, 5.40 Unas varillas cilindricas de acero (AISI 1010), de 50
placa (p) p,, = 8000 kg/m . cp p = 480 J/kg • K. kp = mm de diámetro, se tratan por calentam iento hacién-
12 W/m * K dolas pasar a través de un horno de 5 m de longitud
en el que el aire se mantiene a 750°C I-as varillas en-
Fuerza aplicac Placas de compresión tran a 50°C y alcanzan una tem peratura en la línea
con fluido circulante central de 600°C antes de salir. Para un coeficiente de
1 1 1
Placa metálica convección de 125 W/m • K. estime la velocidad a la
de presión que deben hacerse pasar las varillas a través del homo
t_J 5.41 Estime el tiempo que se requiere para cocinar un hut
-5^mro dog o salchicha en agua hirviendo Suponga que el
Tarjeta de hot dog está inicialmente a 6°C, que el coeficiente de
circuitos transferencia de calor por convección es 100 W /m2 •
t Placa de compresión
K y que la tem peratura final es 80°C en la línea cen-
tral. Trate el hot dog com o un cilindro largo de 20
5.38 Unas tarjetas de circuitos se tratan mediante el calen- mm de diám etro que tiene las propiedades: p — 880
tamiento de una pila de ellas bajo alta presión como kg/m 3, c = 3350 J/kg • K y k = 0 52 W/m • K.
se ilustra en el problema 5.37 Las placas de com pre-
sión en la parte superior e inferior de la pila se m an- 5.42 Una varilla larga de 60 mm de diámetro y propieda-
tienen a una temperatura uniforme med ante un fluido des termofísicas p = 8000 kg/m , c = 500 J/kg • K y
circulante El propósito de la operación de prensado- k = 50 W /m • K. está inicialmente a una temperatura
calentamiento es curar el epóxico. que une las hojas uniforme y se calienta en un horno de convección for-
de fibra de vidrio, e impartir rig dez a las tarjetas. I.a zada que se mantiene a 750 K. Se estima que el coefi-
condición de curado se logra cuando el epóxico se ciente de convección es 1000 W /m 2 • K
mantiene en o por arriba de 170°C durante al menos 5 (a) ¿Cuál es la temperatura de la línea central de la
min. [Link] propiedades term ofísicas efectivas de la pila varilla cuando la temperatura de la superficie es
o libro (tarjetas y placas metálicas de presión) son k 550 K?
= 0.613 W/m • K y pcp = 2.73 X 106 J/m 3 • K. (b )| En un proceso de tratam iento con calor, la tempe-
(a) Si el libro está inicialmente a 15°C y, después de ratura de la linea central de la varilla debe
la aplicación de presión, las placas de compresión aum entar de T, = 300 K a T = 500 K. Calcule y
se llevan de manera súbita a una tem peratura u n i- trace las series históricas de temperaturas de la li-
forme de 190°C, calcule el tiempo te que transcu- nea central para h = 100. 500 y 1000 W /m 2 • K.
rre para que el plano medio del libro alcance la En cada caso el calculo termina cuando T = 500 K.
temperatura de curado de 170°C. 5.43 Un cilindro largo de 30 mm de diámetro, inicialmente
(b) Si en este instante, t = te, la temperatura de las a una tem peratura uniforme de 1000 K, se templa de
placas de compresión se reduce súbitamente a pronto en un gran baño de aceite de temperatura cons-
15°C, ¿cuanta energía tendría que elim inar del li- tante a 350 K. Las propiedades del cilindro son k =
bro el Huido refrigerante en circulación en las pla- 1.7 W /m ■ K, t = 1600 J/kg • K y p = 400 kg/nr3,
cas. a fin de regresar la pila a su temperatura mientras el coeficiente de convección es 50 W/m • K
inicial uniforme? (a) Calcule el tiempo que se requicic para que la su-
perficie del cilindro alcance 500 K.
5.39 Se fonna una capa de hielo en el parabrisas de 5 mm
de espesor de un auto mientras se encuentra estaciona- (b) Calcule y elabore una gráfica de la serie histórica
do durante una noche fría en la que la temperatura am- de temperaturas de la superficie para 0 ^ t ^ 300
biental es -2 0 °C . Al arrancar, con un nuevo sistema s Si se agitara el aceite, y proporcionara un coefi-
desempañante, la superficie interior se expone súbita- ciente de convección de 250 W /m 2 • K, ¿cómo
mente a un flujo de aire a 30°C. Suponga que el hielo cam biaría la historia de tem peraturas}
DEPARTAMENTO DE B IB L IG itü A
Universidad S im ó n B o lív a r - Sede d e l L ito ra
272 C apítulo 5 ■ Conducción en estado transitorio
5.44 Una varilla larga de piroccrám ica de 20 mm de diá- durante 3 min mientras se expone a enfriamiento por
metro se reviste con un tubo metálico muy delgado convección con aire ambiente a 25°C y con un coefi-
para protección mecánica. La unión entre la varilla y ciente de convección de 8 W /m 2 • K? Una condición
el tubo tiene una resistencia de contacto térmico R'¡ c adicional para obtener buenos resultados es una dife
= 0.12 m - K/W. re n d a de tem peraturas máxima-mínima de menos de
Tubo metálico sólido 10°C. ¿Se satisface esta condición? Y, si no, ¿qué ha
ccr para satisfacerla?
Varilla de cerámica
Interfaz de unión
Conducción unidimensional: esfera
.51 Unos cojinetes de bolas de acero inoxidable (AISI ciente de convección a h = 200 W /m 2 • K sobre
304). que se calientan de manera uniforme a 850°C, las historias de las temperaturas anteriores? Para
se endurecen al templarlos en un baño de aceite que se h = 75 y 200 W /m2 • K, calcule y elabore una
mantiene a 40°C. El diám etro de la bola es 20 mm. gráfica del flujo de calor en la superficie como
y el coeficiente de convección asociado con el baño función del tiempo para 0 ^ t ^ 180 s.
de aceite es 1000 W /m 2 • K.
5.54 Las esferas A y B están inicialmente a 800 K. y se
(a) Si el templado no va a ocurrir sino hasta que la templan de manera simultánea en baños de tem pera-
temperatura de la superficie de las bolas alcance tura constante, cada una con temperatura de 320 K.
100°C, ¿cuánto tiempo deben perm anecer éstas Los siguientes parámetros es ín asociados con cada
en el aceite? ¿Cuál es la tem peratura del centro al una de las esferas y sus procesos de enfriamiento.
final del periodo de enfriam iento?
(b) Si se templan 10,000 bolas por hora, ¿cuál es la E sfera A E sfera B
velocidad a la que el sistem a de enfriam iento de
baño de aceite debe quitar energía a fin de mante- Diámetro (mm) 300 30
ner su temperatura a 40°C? Densidad (kg/m 3) 1600 400
5.52 Un granizo esférico de 5 mm de diám etro se forma en Calor específico (kJ/kg • K) 0.400 1.60
una nube de gran altitud a —30°C. Si el granizo co- Conductividad térm ica (W /m • K) 170 1.70
mienza a caer a través de aire mas caliente que está a Coeficiente de convección (W /m2 ■ K) 5 50
5°C, ¿cuánto tiempo transcurrirá antes de que la super-
ficie externa comience a derretirse9 ¿Cuál es la tem pe-
ratura del centro del granizo en este tiempo, y cuánta (a) Muestre de manera cualitativa , en coordenadas T
energía (J) se transfiere al granizo? Se supone un coe- contra t, la temperatura en el centro y en la superfi-
ficiente de transferencia de calor por convección de cie para cada esfera como función del tiempo. Ex-
250 W/m2 • K. y las propiedades del granizo se toman plique de forma breve el razonamiento por el que
de las del hielo. determina las posiciones relativas de las curvas.
(b) Calcule el tiempo que se requiere para que la su-
5.53 Una esfera de 30 mm de diám etro inicialmente a 800
perficie de cada esfera alcance 415 K.
K se templa en un baño que tiene una tem peratura
constante de 320 K con un coeficiente de transferencia (c) Determ ine la energía ganada por cada uno de los
de calor por convección de 75 W /m2 • K. Las propie- baños durante el proceso de enfriamiento de las
dades termofísicas del material de la esfera son: p = esferas a 415 K
400 kg/m. c = 1600 J/kg • K y k = 1.7 W /m • K. 5.55 El coeficiente de convección para el flujo sobre una
(a) Muestre, de forma cualitativa en coordenadas esfera sólida, se determina al sumergiendo la esfera,
T -t, las temperaturas del centro y en la superficie que inicialmente está a 25°C, en el flujo, que está a
de la esfera como función del tiempo 75°C, y m idiendo su temperatura superficial en algún
m omento durante el proceso de calentam iento transi-
(b) Calcule el tiempo que se requiere para que la su-
torio.
perficie de la esfera alcance 415 K.
(a) Si la esfera tiene un diámetro de 0 .1 m, una con-
(c) Determine el flujo de calor (W /m ) en la superfi-
ductividad térmica de 15 W/m • K y una difusivi-
cie externa de la esfera en el tiempo determinado dad térmica de 10~5 m2/s, ¿en qué tiempo se
en la parte (b). registrará una tem peratura superficial de 60°C si
(d) Determine la energía (J) que pierde la esfera du- el coeficiente de convección es 300 W /n r • K?
rante el proceso de enfriado a la tem peratura de la (b )| Evalúe el efecto de la difusividad térmica sobre
superficie de 415 K. la respuesta térmica del material mediante el
(e) Al tiempo determinado por la parte (b), la esfera cálculo de las historias de temperatura en el centro y
se quita rápidamente del baño y se cubre con un la superficie para a = 10 6, 10 5 y 10-4 m2/s.
aislante perfecto, de modo que no hay pérdida de Elabore una gráfica con los resultados para el pe-
calor desde la superficie de la estera. ¿Cuál será riodo 0 ^ t ^ 300 s. De manera similar, evalúe el
la temperatura de la estera después de que trans- efecto de la conductividad térmica mediante la
curre un largo tiempo? consideración de valores de k = 1.5, 15. y 150
W m ■ K.
|(f)~| Calcule y trace las historias de las temperaturas
del centro y de la superficie en un periodo 0 < t ^ 5.56 En un proceso para fabricar cuentas de vidrio (k =
150 s. ¿Qué efecto tendrá un aumento en el coefi- 1.4 W m • K, p = 2200 kg/m3 y cp = 800 J/kg • K) de
miintinito que tiene las mismas propiedades termofísi- lesiones inducidas cuando una parte del cuerpo de un
cas que el medidor. trabajador llega a hacer contacto con maquinaria que
esta a temperaturas elevadas en el rango de 50 a
0 d
—i- 100°C El asesor médico les informa que ocurrirá una
lesión térmica irreversible (muerte de la célula) en
cualquier tejido vivo que se mantenga a T s 48°C du-
rante Ar > 10 s Quieren información con respecto al
_Mismo material grado de daño irreversible del tejido (medido por la
que el medidor distancia desde la superficie de la piel) como función
de la temperatura de la maquinaria y el tiempo duran-
te el cual se tiene contacto entre la piel y la m aquina-
ria. Suponga que el tejido vivo tiene una temperatura
normal de 37°C. es isotrópico y tiene propiedades
Desarrolle relaciones para predecir la constante de constantes equivalentes a las del agua líquida
tiempo del medidor en los dos casos y compárelas
con la relación del fabricante ¿Que conclusión extrae (a) Para evaluar la seriedad del problema, calcule lu-
de este análisis con respecto a la respuesta transitoria de gares en el tejido en los que la temperatura alcan-
los medidores para diferentes aplicaciones? zará 48°C después de 10 s de exposición a la
maquinaria a 50°C y 100°C
.63 Un procedimiento simple para m edir coeficientes de
\r.
transferencia superficial de calor por convección im (b) Para una temperatura de la maquinaria de 100°C
plica cubrir la superficie con una capa delgada de m a- y 0 ^ t ^ 30 s, calcule y trace las historias de la
terial que tenga una tem peratura precisa del punto de temperatura cn lugares de tejido a 0.5, I, 2 y 5
fusión. Después se calienta la superficie y. mediante mm de la piel
la determinación del tiem po que se requiere para que 5.65 Un procedimiento para determ inar la conductividad
ocurra la fusión, se determina el coeficiente de con- térmica de un material sólido implica incrustar un ter-
vección. El sigu ente arreglo experimental utiliza el m opar en una placa gruesa del solido y medir la res-
procedimiento para determ inar el coeficiente de con- puesta a un cam bio establecido cn la tem peratura en
vección para un llujo de gas normal a la superficie. una superficie. Considere un arreglo en que el termo-
De manera específica, una varilla larga de cobre se par se incrusta 10 mm desde una superficie que de sú-
cubre con un superaislante de conductividad térmica bito se lleva a una tem peratura de 100°C mediante la
muy baja, y se aplica una capa muy delgada a la su- exposición a agua en ebullición. Si la temperatura ini-
perficie expuesta. cial de la placa fue 30°C y el term opar mide una tem -
peratura de 65°C, 2 min después de que la superficie
se lleva a I0()°C, ¿cuál es su conductividad térmica?
Rujo de ga*
k L ík ? Se sabe que la densidad y el calor específico del sóli-
do son 2200 kg/m3 y 700 J/kg • K.
5.66 Un calentador eléctrico en forma de lámina se coloca
Superficie recubierta en contacto firme con la superficie de una placa grue-
sa de baquelita que tiene una temperatura uniforme de
Varilla de cobre.
k = 400 W/m • K. « = 10< m2/s 300 K Determine la temperatura de la placa cn la su-
perficie y a una profundidad de 25 mm, 10 min des-
Superaislante pués de que se energiza el calentador y proporciona un
llujo de calor constante a la superficie de 2500 W /m2.
5.67 Una placa muy gruesa con difusividad térmica 5.6 X
Si la varilla esta inicialmente a 25°C y se pasa un flu- 10~6 m2/s y conductividad térmica 20 W/m • K está
jo de gas para el que h = 200 W m- • K y f » = inicialmcnte a una temperatura uniforme de 325°C.
300°C, ¿cuál es la tem peratura del punto de fusión del De pronto, la superficie se expone a un fiuido refrige-
recubrimiento si se observa que la fusión ocurre en rante a 15°C cuyo coeficiente de transferencia de ca-
l = 400 s? lor por convección es 100 W/m • K
564] Una compañía de seguros lo contrata a usted como (a) Determine las temperaturas en la superficie y a
asesor para entender y saber más sobre las lesiones una profundidad de 45 mm después de transcurri-
por quemaduras Están interesados en especial en las dos 3 min.
DEPARTAMENTO DE BIBLIOTECA
U n iv e rs id a d S im ó n B o lív a r • S e d e ^ 1 * 5
276 C apítulo 5 ■ Conducción en estado transitorio
(b) Calcule y elabore una gráfica de las historias de sustrato (k¡, a f), la densidad y el calor latente de
tem peratura (0 < t < 300 s) en i = 0 y x = 45 fusión del depósito (p, /ivy), el espesor del áep
mm para las siguientes variaciones de los para to 8 y las temperaturas relevantes ( 7 . 7 )
metros (i) a = 5.6 X 10-7, 5 6 X 10-6, y 5.6 X
(b) El proceso de deposición de plasma pulverizado
I0 “5 m 2/s, y (2) k = 2, 20; y 200 W /m • K
del problema 5 23 se usa para aplicar un recubi
5.68 Una pared gruesa de roble inicialmente a 25°C, se ex miento de alúmina delgado 6 = 2 mm) sobra i
pone de pronto a productos de com bustión para los sustrato grueso de tungsteno. El sustrato tiene
que 7 , = 800°C y h = 20 W /m • K. una tem peratura mcial uniforme 7 = 300 K, v
(a) Determine el tiempo de exposición que se requie su conductividad y ditusividad térmicas se ap
re para que la superficie alcance la temperatura ximan como k = 1 2 0 W/m • K y a - 4.0 X
de ignición de 400°C. 10 5 n r s, respectivamente. La densidad y ele»-
(b) Elabore una gráfica de la distribución de temperatu- lor latente de la alumina son p = 3970 kgm!
ras 7 ( 0 en el medio a / = 325 s La distribución de- hsf = 3577 kJ/kg. respectivamente, y la alumi
be extenderse a una posición en la que t ~ 25°C. se solidifica a su tem peratura de fusión (T
2318 K) Suponiendo que la capa fundida se dt
5.69 Es bien sabido que, aunque dos materiales esten a la
posita instantáneamente sobre el sustrato esf
misma tem peratura, uno se siente más trio al tacto
el tiempo que se requiere para que el depos tu
que el otro Considere placas gruesas de cobre y vi
solidifique.
drio, cada cual a una tem peratura inicial de 300 K
Suponiendo que el dedo está a una temperatura inicial 5.72 emando un metal tundido se vacía en un molde qu-
de 310 K y que tiene las propiedades term ofísicas p = un conductor pobre, la resistencia al flujo de caloré
1000 kg/m 3, r = 4180 J/kg • K y k = 0.625 W/m • K, minante esta dentro de la pared del molde. Cons
determine si el cobre o el vidrio se sentirá más frío al condiciones en las que un metal líquido se solí
tacto. en un molde de paredes gruesas de conductividad
5.70 Dos placas de acero inoxidable (p = 8000 kg m . c = mica kw y difusividad térmica La dens dad )
500 J/kg • K, k = 15 W /m • K). cada una de 20 mm calor latente de tusion del metal se designan CUfl.
de espesor y aisladas en una superficie, están inicial- y respectivamente y en ambos estados, i'undii
mente a 400 y 300 K cuando se presionan una a otra sólido, la conductividad térmica del metal es mu
por sus superficies no aisladas. ¿Cuál es la temperatu- mas grande que la del molde.
ra de la superficie aislada de la placa caliente después
de transcurrido 1 min?
5.71 Los recubrimientos especiales a menudo se forman
depositando capas delgadas de un material fundido
sobre un sustrato solido La solidificación comienza
en la superficie del sustrato y continua hasta que el
espesor S de la capa solida se hace igual al espesor 8 Pared del molde
del depósito M*W
Líquido
resultado en térm inos de los parám etros apropia- coeficiente de transferencia de calor por convección
dos del sistema. de 15 W /n r • k Trate la carne com o un cilindro con
propiedades de agua líquida, que tiene un diámetro
igual a su longitud.
o i k I ik m m Ó i i niiillidiitiriisioiiiil
5.77 L na varilla larga de 20 mm de diám etro fabricada de
5.73 Ln lingote largo de acero (carbón hom ogéneo) de
alum ina (óxido de alum inio policristalino) está ini-
sección transversal cuadrada de 0.3 m por 0 3 m, ini-
cialm ente a una tem peratura uniform e de 850 k La
cialm ente a una tem peratura uniform e de 30°C, se co-
varilla se expone de súbito a un fluido a 350 K con h =
loca en un horno de im pregnación térm ica que tiene
500 W /n r • K Estim e la tem peratura de la línea cen-
una tem peratura de 750°C Si el coeficiente de trans-
tral de la varilla después de 30 cn un extrem o e x -
ferencia de calor por convección para el proceso de
puesto y a una distancia axial de 6 mm del extrem o
calentam iento es 100 W m • k . ¿cuánto tiem po debe
perm anecer el lingote en el horno antes de que la tem -
peratura de su centro alcance 600°C ? cn d ife r e n c ia * finitas: d e r iv a c ió n
5.74 Un ladrillo refractario de dimensiones 0 06 m X 0 09 m 5.78 El criterio de estabilidad para el m étodo explícito re-
X 0.20 m se quita de un horno de calcinación a 1600 quiere que el coeficiente del térm ino Tf„ de la ecua-
K y se enfría en aire a 4tUC con h = 50 VV m ■ K ción unidim ensional en diferencias finitas sea cero o
¿Cuál es la tem peratura cn el centro y en las esquinas positivo. Considere la situación cn que las temperara
del ladrillo después de 50 min de enfriam iento? ras en los dos nodos vecinos (7¿’_ |. / ’£+ i) s*>n 100°C
5.75 Lna punta cilindrica de cobre de 100 mm de longitud y m ientras que el nodo central (Tj¡,) está a 50°C M ues-
50 mm de diám etro está inicialm cnte a una tem peratura tre que para valores de f o > \ . la ecuación en dife-
unilormc de 20°C. Las caras de los extrem os se somc rencias finitas predecirá un valor de T £ ' 1 que viola la
ten de pronto a una intensa rapidez de calentam iento segunda le> de la term odinám ica
que las e le \a a una tem peratura de 500 C. Al mismo 5.79 Una varilla delgada de diám etro D esta inicialmentc
tiempo, la su|>erlicie cilindrica se som ete a calenta- cn equilibrio con sus alrededores, un recinto grande al
miento por un flujo de gas con una tem peratura de vacío a tem peratura T ^ . De súbito se hace pasar una
5LK)°C y un coeficiente de transferencia de calor de 100 co m en te eléctrica / (A) por la varilla que tiene una re-
W /nr • k sistividad eléctrica pt y em isividad f . En la hgura se
o de identifican otras propiedades term ofísicas pertinentes
Derive la ecuación en diferencias finitas transitoria
para el nodo m.
i i 1 1 1
A
7 *-
Escriba la ecuación en diferencias imitas que exprese (a) Con el método explícito con Fo = 2 derive las
la conservación de la energía para el nodo 0 que se lo- ecuaciones en diferencias finitas para los nodos J.
caliza en la superficie externa en x = —L. Rcacomode 2. 3 y 4.
su ecuación e identifique cualesquiera coeficientes (b) Construya una tabla que tenga como encabezados
adim ensionales importantes. /?, r, y las temperaturas nodales T0 a T4 Detenn ie
5.81 Una paicd plana (p = 4000 kg/m , cp = 500 J/kg * K. k = la temperatura de la superfic c T0 cuando f4 =
10 W /m • K) de espesor L = 20 mm tiene inicialmente 35°C
una distribución lineal de tem peraturas en estado esta-
ble con límites que se mantienen a = 0°C y T2 = Solucionéis cn diferencia* finita*:
100 C. De sub to. se hace pasar una corriente eléctrica sistema* unidimensionales
a través de la pared que ocasiona una generación uni-
forme de energía a una razón q = 2 X 10 W /m3. Las 5.83 Una pared de 0 12 m de espesor que tiene una difusj.
condiciones de frontera T ] y T2 perm anecen fi jas. vidad térm ica de 1.5 X 10~6 n r/s está inicialmcnte;
una tem peratura uniforme de 85 C. De pronto una ca-
A = 0°C- ■72= 100°C ra se baja a una tem peratura de 20 C mientras la oís
cara queda perfectamente aislada.
1 m
- t i rO.c/ = 2 x 107W'm3 (a) Con la técnica de diferencias finitas explícita'
con incrementos de espacio y tiempo de 30
y 300 s, respectivamente, determine la distrfi
u L = 20 mm
ción de temperaturas a t = 45 mm.
(a) En coordenadas T -x . dibuje las distribuciones de
(b) Con Av = 30 mm y At = 300 s calcule T[ ¡
tem peraturas para los siguientes casos: (i) condi-
para 0 ^ t ^ tss, donde /ss es el tiempo que se
ción inicial (t ^ 0); (u) condiciones de estado esta
quiere para que la temperatura en cada punto 1
ble (/ —» oc). suponiendo que la temperatura
dal alcance un valor que esté a 1°C de
máxima en la pared excede 7 \; y (iii) para dos
tem peratura de estado estable. Repita los cálen
tiempos intermedios. Acote todas las característi-
precedentes para At = 75 s Rara cada valor
cas importantes de las distribuciones.
Ar, trace las historias de la temperatura para 1
(b) Para el sistema de tres puntos nodales que se mues- cara y para el plano medio.
tra de forma esquemática (1. m, 2), defina un volu-
men de control apropiado para el nodo m y, con la 5.84 La pared plana del problema 2.43 (k = 1.5 W/m • |J
identificación de todos los procesos relevantes de- a = 7 5 X 10~6 m2/s) tiene un espesor de L - 501
rive la ecuación en diferencias finitas correspon- y una tem peratura inicial uniforme de 25°C. De si
diente siguiendo el método expl cito o el implícito. to, el límite en v = L experim enta calentamiento o
b do a un fluido para el que T x = 50°C y h =
(c) Con un incremento temporal de Ar = 5 s. use el mé-
W /m 2 • K. mientras que el otro límite en x = 0e$
todo de diferencias finitas para obtener valores de
rim enta un flujo de calor aplicado de </¿ =
T„, de los primeros 45 s de tiempo transcurrido De-
W m 2.
termine los flujos de calor correspondientes en las
fronteras, es decir, c¡\ (0 ,4 5 s) y í/' (20 mm, 45 s). (a) Con Av = 5 mm y Ar = 20 s. calcule y dibuje (
distribuciones de temperatura^ en la pared
(d) Para determ inar el efecto del tamaño de la malla
i) la condición nicial, (ii) la condición de 1
repita su análisis con rejillas de 5 y 11 puntos no-
estable y ( 111 ) dos tiempos intermedios.
dales Av = 5.0 y 2.0 mm. respectivamente)
(b) bn coordenadas q'[ —x. trace las distribucioi
5.82 Un cilindro circular solido de material plástico (a = 6 X
del flujo de calor que corresponden a las
10“7 m2/s) está rtiicialmente a una temperatura uniforme
distribuciones de temperaturas que se represe
de 20°C y esta bien aislado a lo largo de su superficie la-
en la parte (a).
teral y en un extremo AI tiempo t = 0, se le aplica calor
(c) En coordenadas q"x —/. elabore una gráfica 1
al extremo izquierdo lo que ocasiona que T0 aumente li-
jo de calor en v = 0 y x = L
nealmente con el tiempo a una razón de I°C/s.
.. . «I*.", .v 5.85 La pared plana del problema 2.44 (k = 50 W/m ■y
= 1.5 X 10 6 m2/s) tiene un espesor de L = 401
una tem peratura inicial uniforme de T0 = 25°C.
tamente, la frontera en x = L experimenta ca
miento por un fluido para el que T» = 50°C \l
1000 W /m2 ■ K, mientras se genera calor de
uniforme dentro de la pared a q = I X 107 Wm
L = 24 mm frontera en x — 0 permanece a Ta.
Problemas 279
(a) Con Av = 4 mm y Ai = 1 s. trace las distribucio- vana y la velocidad del ventilador cambia para aumen-
nes de tem peraturas en la pared para (i) la condi- tar sustancialmente la circulación del aire dentro del
ción inicial, (ii) la condición de estado estable y horno, la superficie interna del horno experimenta
(iii) dos tiempos intermedios. un cam bio de tem peratura súbito a 100°C. El aislan-
te tiene una conductividad térm ica de 0.03 W /m • K
(b) En coordenadas q"x—t, dibuje el (lujo de calor en x
y una difusividad térm ica de 7.5 X 10-7 n r/s . Para
■ 0 y v = L. ¿En que tiem po transcurrido hay
la solución en diferencias finitas, use un incremento
flu jo de calor cero en x = L?
espacial de 6 mm. Suponga que el efecto de las ho-
5.86 Considere el elemento de combustible del ejemplo jas de acero inoxidable es insignificante y que el
5.7. Inicialmente, el elemento está a una temperatura coeficiente de transferencia de calor por convección
uniforme de 250 C sin generación de calor. De súbito, exterior ha perm anece sin cam bio. Estime el tiempo
el elemento se inserta en el núcleo del reactor, lo que que se requiere para que la pared del horno se apro-
ocasiona una rapidez de generación volumétrica de ca- xime a condiciones de estado estable después de que
lor de q = 108 W /m 3. Las superficies se enfrían de la tem peratura de la pared interior cam bia a 100°C.
forma convectiva con Tx = 250°C y Ii = 1100 W/m2 • K
Con el método explícito y un incremento espacial de Ventilador circular
2 mm. determine la distribución de tem peraturas 1,5 s "l |~ Montaje del calentador
—HI mm h—
I
después de que el elem ento se inserta en el núcleo.
DEPARTAMENTO DE BIBLIOTECA
U n iv e r s id a d S im ó n B o lív a r • B edti u
280 C apítulo 5 ■ Conducción en estado transitorio
5.92 Considere la operación de unión que se describe en el 5.94 Una varilla de tantalio de 3 mm de diámetro y 120
problema 3.97, que se analizó en condiciones de esta- mm de longitud se sostiene con dos electrodos dentro
do estable. En este caso, el láser se utilizará para ca- de un recinto grande al vacio Inicialmente la varilla
lentar la película en un tiempo establecido, con lo que esta en equilibrio con los electrodos y sus alrededo-
se crea la situación de calentam iento transitorio que se res, que se mantienen a 300 K De súbito, una corriente
muestra en el dibujo. eléctrica, I = 80 A. pasa a través de la var lia vSupon-
ga que la emisividad de la varilla es 0 1 y que la resis-
tividad eléctrica 95 X !0 -8 i l • ni. Use la tabla A.l
para obtener las otras propiedades termofísicas que se
requieren para la solución. Utilice un método en dite
nuil
Fuente láser, q''
rencias finitas con un incremento espacial de 10 mm,
Película plástica
Tira metálicT
( Electrodo,
300 K
Vanlla
La tira está inicialmente a 25°C y el láser proporcio-
na un flujo uniforme de 85,000 W /m en un tiem po
r
A/on = 10 s Las dim ensiones del sistem a y propieda- Electrodo,
des term ofísicas perm anecen iguales, pero el coefi- 300 K
ciente de convección al aire am biente a 25°C es Alrededores, T&
ahora 100 W /m 2 ■K.
(a) Con el método implícito en diferencias finitas para (a) Estime el tiempo que se requiere para que la i
Aa = 4 mm y A/ = 1 s, obtenga las series históri- tad de la longitud de la vanlla alcance 1000 K
cas de temperatura para 0 < t < 30 s en el centro y (b) Determine la distribución de temperaturas dei
el borde de la película. 7(0. t) y T(w {/2, r), respec- tado estable y calcule de manera aproximad!
tivamente. para determinar si el adhesivo se cura cuánto tiem po tardará la varilla en alcanzares
de manera satisfactoria por arriba de 90°C para 10 s condición.
y si excede su temperatura de degradación de
5.95 Una vanlla de apoyo (k = 15 W /m • K, a = 4.0
200°C.
10-6 m 2/s) de diám etro D = 15 mm y longitud/
(b) Valide el código de su program a m ediante la 100 mm atraviesa un canal cuyas paredes se manti
comparación con los resultados de estado estable nen a una tem peratura Th = 300 K. Súbitamente I
del problema 3.97. ¿Qué tipo de solución analíti- varilla se expone a un flujo cruzado de gases calien
ca buscaría a fin de probar el com portam iento para los que Tx = 600 K y h = 75 W/m • K Las i
transitorio apropiado de su código? redes del canal se enfrían y permanecen a 300 K
5.93 Un extremo de una varilla de acero inoxidable (AISI
316) de 10 mm de diám etro y 0.16 m de longitud se
inserta en una montura que se mantiene a 200°C. La va-
rilla, cubierta con una manga aislante, alcanza una
7>,= 3iX'i
tem peratura uniforme en toda su longitud Cuando se / / /
quita la manga la varilla queda expuesta al aire am -
biente a 25°C de modo que el coeficiente de transfe-
rencia de calor por convección es 30 W /m2 • K.
— Varilla, D = 15 mm, L = 100 mm
(a) Con la técnica de diferencias finitas explícita y un
incremento espacial Aa = 0 016 m, estim e el (a) Con una técnica numérica apropiada detc n
tiempo que se requiere para que la parte media de respuesta ténnica de la varilla al calent mx
la longitud de la varilla alcance 100°C. convectivo. Dibuje la temperatura de la mitadc
(b) Con Aa = 0.016 m y Ai = 10 s. calcule 7"(a , t) mo función del tiempo transcurrido. Con uní
para 0 < t < t\, donde t l es el tiempo que se re- lo analítico apropiado de la vanlla, determiij
quiere para que la parte media de la longitud de distribución de temperaturas de estado estafci
la varilla alcance 50°C. Dibuje la distribución compare el resultado con el que se obtiene ni
de temperaturas para t = 0, 200 s, 400 s y r,. camente para tiempos transcumdos muy
■ Problemas 281
(b) Después de que la varilla alcanza las condiciones (a) Usando un incremento de tiempo At = 60 s y el
de estado estable, el flujo de gases calientes se sus- método implícito, determine la historia de la tem -
pende súbitamente y la varilla se enfria por con- peratura del plano medio del libro y vea si ocurri-
vección libre al aire ambiental a i » = 300 K y por rá el curado (170°C por 5 min).
intercambio de radiación con los alrededores a Talr —
300 K. El coeficiente de convección libre se expre- ib) Siguiendo la reducción de las temperaturas de los
rodillos a 15°C (i = 50 min). ¿cuánto tiempo tar-
sa como /i(W /m2 • K ) = C A T". donde C = 4.4
dará el plano medio del libro en alcanzar 37°C,
W/m2 • K l 188 y n = 0.188. La emisividad de la va-
una temperatura segura a la que el operador pue-
rilla es 0.5. Determine la respuesta térmica poste-
de com enzar a descargar la prensa?
rior de la varilla. Trace la tem peratura del medio
como función del tiempo de enfriamiento, y deter- (c) Valide su código de programa con el calendario
mine el tiem po que se requiere para que la varilla de calentam iento de un cambio súbito de la tem-
alcance una temperatura segura al tacto de 315 K peratura del rodillo de 15 a 190°C y compare sus
Considere la hoja rejilla de aceleración (k = 40 W/m resultados con los de una solución analítica apro-
• K a = 3 X 10-5 m2/s, e = 0.45) del problem a 4.70. piada (véase el problema 5.38).
Desarrolle un modelo en diferencias finitas implícito
de la hoja que sirve para los siguientes propósitos.
(a) Suponiendo que la hoja está a una temperatura K ru a rio iie s en clifereneias finita*:
uniforme de 300 K cuando se activa la fuente del coordenadas cilindricas
haz de iones, obtenga una gráfica de histórica pos-
5.98 Un disco circular delgado está sujeto a calentamiento
terior de tem peratura-tiem po del tramo medio. ¿A
por inducción de una bobina, el efecto de la cual es
qué tiempo transcurrido alcanza este punto de la
proporcionar una generación de calor uniforme dentro
hoja una tem peratura a 1 K del valor de estado
de una sección anular, como se muestra. La convec-
estable?
ción ocurre en la superficie superior, mientras la su-
(b) La hoja se opera en condiciones de estado estable perficie inferior está bien aislada.
cuando, de pronto, se desactiva el haz de iones.
Obtenga una gráfica de la historia temperaturas-
tiempo del medio. ¿Cuánto tiempo transcurre pa-
ra que el punto más caliente de la hoja se enfríe a
315 K. que es un estado seguro al tacto?
5.97 l'nas tarjetas de circuitos se tratan m ediante el calen-
tamiento de una pila de ellas bajo alta presión como
se ilustra en el problema 5.37 y se describe, además,
en el problema 5.38. Se busca un método de solución en
diferencias finitas con dos consideraciones adiciona-
les. Primera, el libro se tratará como si tuviese caracte-
rísticas distribuidas, en lugar de concentradas, con un
espaciado de rejilla A x = 2.36 mm con nodos en el
centro de la tarjeta o placa de circuitos individual. Se-
gunda. en lugar de elevar la temperatura de las placas
a 190°C mediante un cam bio súbito, se utilizará el
programa de calentamiento que se muestra a continua-
ción a fin de minimizar las tensiones térmicas excesi-
vas inducidas por los gradientes térmicos rápidamente (a) Derive la ecuación en diferencias finitas, transito-
cambiantes en la vecindad de los rodillos. ria, para el nodo m. que está dentro de la región
sujeta al calentam iento por inducción.
v. n
■. 2, 1 k.------ 3-’t 111
x, m
A* = Av = 20 mm
CAPÍTULO
Introducción
a la convección
284 (Capítulo 6 ■ Introducción a la convección
q = (Ts - TV)
(a) (/>)
Fl<;i HA 6 . 1 Efe( tos d r la Irunsferenc id local > total do calor por ronvereirin.
(«) Superficie de forma arbitraria. (b) Flaca plana.
o, de la ecuación 6 .1
Definiendo un coeficiente de convección prom edio h para lodu la superficie, el calor to-
tal transferido se expresa como
Al igualar las ecuaciones 6.3 y 6.4, se sigue que los coeficientes de convección prome-
dio y local están relacionados por una expresión de la forma
h — — | h dA s (6.5)
a*
Advierta que para el caso especial de flujo sobre una placa plana (figura 6.16), h varía
con la distancia x desde la primera orilla y la ecuación se reduce a
h = —1 IrLh d x (6.6)
L Jo
Es posible obtener resultados similares para la transferencia de masa por convec-
ción. Si un fluido de concentración molar de especies CA<« fluye sobre una superficie
en la que la concentración de especies se mantiene en algún valor uniforme CA 5 #
C a . « (figura 6.2a), ocurrirá una transferencia de especies por convección. La especie A
es normalmente un vapor que se transfiere en un flujo de gas debido a la evaporación o
sublimación de una superficie líquida o sólida, respectivamente, y estamos interesados
en determinar la velocidad a la que ocurre esta transferencia. Como en el caso de la
transferencia de calor, este cálculo se basa en el uso de un coeficiente de convección
[Link] -
(a) (h)
F l C I RA 6 . 2 E f e c l o » d e la t r a n s f e r e n c ia Itx a l y lo la l d e e sp r< íe s p o r
e o n v e r riilii. (a) S u p e r f i c i e d e fo rm o a r b it r a r ia . ( ¿ ) P iu c a p la n a .
286 C apítulo 6 ■ Introducción a la convección
N * = h „ (C (6
6
( .
A s3 JA
*
Para la placa plana de la figura 6.2 b, se sigue que
11 ff L
hm = 7 1
h.„ d x (6.1
La transferencia de especies también se expresa como un llujo de masa.
n r), o como una transferencia de masa. nA (kg/s), multiplicando ambos lados de!
ecuaciones 6.7 y 6.8, respectivamente, por el peso molecular J t A (kg/kmol) de lae
e'ie A. En consecuencia.
^A h/tf p A s Pa. dc) 6.
y
n A = hmA s(pA s - pA, J
donde p A (kg/m3) es la densidad de masa de la especie A.1
Para llevar a cabo un cálculo de transferencia de masa por convección, es ne
rio determinar el valor de CA ^ o p A s. Tal determinación se hace fácilmente n
que existe equilibrio termodinámico en la interfaz entre el gas y el líquido o fase
da. Una consecuencia de esta condición de equilibrio es que la temperatura del
en la interfaz es igual a la temperatura superficial 7's. Una segunda implicación e*
el vapor se encuentra en un estado saturado , en cuyo caso las tablas termodinám!
como la tabla A.6 para el agua, sirven para obtener su densidad a partir del cq
miento dc T y Con una buena aproximación, la concentración molar del vapórenla
períicie también se determina de la presión de vapor a través de la aplicación
ecuación de estado para un gas ideal. Es decir.
P sJT J
Ca .s
9hT..
'Aunque la nomenclatura anterior es bastante adecuada para caracterizar los procesos de transferencia de maiat
este texto, no existe nomenclatura estándar alguna, por lo que a menudo es difícil reconciliar los resultados de
blicaciones. Webb 11 ] proporciona una revisión de las diferentes formas en las que es posible formular poiencute
tores. flujos y coeficientes de convección.
6 .1 ■ El problema de la transferencia de calor por convección 287
Kjivviri.o 6 .1
SOLI CIÓN
E n c o n tr a r :
1. La razón del coeficiente promedio de transferencia de calor h(x) al valor local h x(x).
2. Dibuje la variación de h x y hx con v.
E squem a:
Capa límite
hx - as 01
Análisis:
1. De la ecuación 6.6 el valor promedio del coeficiente de transferencia de calor por
convección sobre la región de 0 a v es
- 1 (x
h x = hx(x) = — h fx ) dx
X Jo
DEPARTAMENTO DE BIBLIOTECA
Universidad Sunco boi*[Link] - Sed*- ít cwp'
C apítulo 6 ■ Introducción a la convección
e integrar, obtenemos
1 rx a cx ,+09
hx = — I a x 01 d x = - x~° 1 d x = - = 1.11 a x 01
x •'O X Jo X 0 .9
o
iix = l . l l / i v
0
0 1 2 3 4
a (m)
C om entarios: El desarrollo de la capa limite ocasiona que los coeficientes local y p
medio disminuyan al aumentar la distancia desde la orilla, por tanto, el coeiici
promedio hasta \ debe excedei el valor local en x.
Ej e m pl o 6 .2
Son <;ioi\
S a p o sic io n e s:
a
Capas límite de convección
Para introducir el concepto de una capa límite, considere el flujo sobre la placa plana
de la figura 6.3. Cuando las partículas del fluido hacen contacto con la superficie, ad-
quieren una velocidad cero. Estas partículas actúan entonces para retardar el movi-
miento de partículas en la capa contigua del fluido, que a su vez actúa para retardar el
movimiento de las partículas en la siguiente capa, y así sucesivamente hasta que, a una
distancia y = 8 de la superficie, el efecto se hace insignificante. Este retardo o desace-
leración del movimiento del fluido se asocia con los esfuerzos cortantes r que actúan
en planos que son paralelos a la velocidad del fluido (figura 6.3). Al aumentar la dis-
tancia y desde la superficie, el componente .v de la velocidad del fluido, u, debe enton-
ces aumentar hasta que se aproxima al valor del flujo libre ux . Se usa el subíndice °c
para designar las condiciones en el flu jo libre fuera de la capa limite
La cantidad 8 se denomina espesor de la capa límite y normalmente se define co-
mo el valor de y para el que u — 0 .99u<x>- El perfil de velocidad de la capa lim ite se re-
fiere a la forma cn la que u varía con y a través de la capa limite. En consecuencia, el
flujo del fluido se caracteriza por dos regiones distintas, una capa fluida delgada (capa
límite) cn la que los gradientes de velocidades y los estuerzos cortantes son grandes y
Flujo libre
8 (A)
«oo
,iú f Capa límite de
- * J_ _ J velocidad o F ic ü k a 6 .3
—►
—► i y ~ hidrodinámica D e s a r r o llo d e la c a p a lím ite d e
_ _ v e l o c id a d o h id r o d in á m ic a s o h r e u n a
p la n a .
DEPARTAMENTO DE BIBLIOTECA
UnlvarsiucKi o ih iu i L>oii«jr - Serie d ‘ 'oral
290 C apítulo 6 ■ introducción a la convección
una región fuera de la capa limite en la que los gradientes de velocidad y los esfuerce*
cortantes son insignificantes. Con el aumento de la distancia desde la primera orilla, los
efectos de la viscosidad penetran más en el flujo libre y la capa limite crece (fi aumen
la con x).
Como se relaciona con la velocidad del fluido, la capa límite anterior se denomina,
de manera mas específica, capa lim ite de velocidad o hidrodinám ica. Se produce siem-í
pre que hay un Mujo de fluido sobre una superficie y es de fundamental importancia pa-
ra problemas que incluyen transporte por convección. En la mccanica de fluidos, su]
significado para el ingeniero proviene de su relación con el esfuerzo cortante supertjJ
cial r y, en consecuencia, con los efectos de fricción de la superficie En cuanto a lid
jos externos, proporciona la base para determinar el coeficiente clefricci n local
Cf = — ;—
puJ2
du'
(6.
v= D
'
donde /x es una propiedad del fluido que se conoce como vise isid a d dinámica.
F ig l r \ 6 .1
l’r t i i i ( t* a ( apa límite
soh re una { a a plana ¡sotrnuJ
6 .2 ■ Capas limite de convección 291
(6.16)
Por ello las condiciones en la capa límite térmica, que influyen fuertemente en el gra
diente de temperatura de la pared r)7Yr)v|v= ü, deteiminan la transferencia de calor a ira
ves de la capa límite Como (T s — T f) es una constante, independiente de v, mientras 8,
se incrementa al aumentar los gradientes de temperatura en la capa limite deben dis
minuir al aumentar x. En consecuencia, la magnitud de c)T/dv\v =0 disminuye al aumen-
tar v, y se sigue que q" y h disminuyen al aumentar .v.
Así como las capas límites hidrodinámica y térmica determinan la fricción de la pared
y la transferencia de calor por convección, la capa ¡imite de concentración determina
la transferencia de masa por convección. Si una mezcla binaria de las especies quimi
cas A y B fluye sobre una superficie y la concentración de la especie A en la superficie,
C A 5, difiere de la concentración en el flujo libre, CA «, (figura 6.i), se producirá una
capa limite de concentración. Esta es la región del fluido en la que existen gradientes
de concentración, y su espesor d se define normalmente como el valor de y para el que
[(CA. S — CA)/(CA 5 — CA ce)] ~ 0 99 La transferencia de especies por convección en-
tre la superficie y el flujo libre de fluido está determinada por las condiciones en esta
capa limite
La ielación entre transferencia de especies por convección y la capa limite de con-
centración se demuestra al reconocer primero que el flujo molar asociado con la trans-
ferencia por difusión se determina mediante una expresión analoga a la ley de Fouricr.
Para las condiciones de ínteres en este texto, la expresión, que se denomina ley de
F ick , tiene la forma
dCA
(6 18)2
¿y
:Estn expresión resulta de una forma más general de la ley de difusión de Fick (sección 14.1.2) cuando la concentración mo-
lar total de la mezcla. C = CA + CB. es una constante.
Flujo libre
Mezcla —
de A B .. C.a, «>
8 c(x)
F it a ha 6 .5
E v o l u c i ó n tic la c a p a lim ito d e
c o n c e n tra c ió n d e e s p e c ie s so b re un a
p l a c a p la n a .
292 C apítulo 6 ■ Introducción a la convección
donde Dab es una propiedad de la mezcla binaria que se conoce como coeficiente i
difusión binario. En cualquier punto que corresponda a y > 0 en la capa límite de i
centración de la figura 6.5, la transferencia de especies se debe al movimiento gl<
de fluido y a la difusión. Sin embargo, en y = 0 no hay movimiento de fluido yB
transferencia de especies es sólo por difusión. Al aplicar la ley de Fick en y = 0. el flu-
jo de especies a cualquier distancia desde el inicio de la superficie es entonces
bj" — _ n ¿_
_CAA
A — AH í(
¿y v=0
— ^ A B ¿ ^ A !d y \ v= o
C a . s ~ C a . oo
Por tanto, las condiciones en la capa límite de concentración, que influyen grander
en el gradiente de concentración de la superficie dCA/dy|v= 0, influirán en el coefic
te de transferencia de masa por convección y, por ello, en la transferencia de espeo
la capa límite.
Los resultados anteriores también se expresan en una base de masa, en lugar!
molar. Al multiplicar ambos lados de la ecuación 6.18 por el peso molecular de last
pecies M a , el flujo de masa de especies debido a la difusión es
¿P a
n 'k — £*ab
¿y
Con la aplicación de esta ecuación en y = 0 y al combinar con la ecuación 6.11,<
nemos
‘^ a b ¿ P a / ¿ > ' U o
P a , s ~ P a ,«o
Ej empl o 6 .3
So l u c ió n
E sq u em a :
S u p o s ic io n e s :
1. El vapor de agua se aproxima como un gas ideal.
2. Las condiciones son isotérmicas.
Tabla A.6, vapor saturado (0.1 atm = 0.101 bar): Ts = 319 K. Tabla
P r o p ie d a d e s :
A.8. vapor de agua-aire (319 K): D AB(319 K) = DAB(298 K) X (319 K/298 K )m =
0.288 X 10 4 m2/s.
_ ~ ¿ > A B d / ? A % ly = 0
™m%x_____________
Pa,s
_
P a. 20
P a = Pa R T
con 1 constante (condiciones isotérmicas),
^ a b d p A/dy\y=o
km x P A. 5 P A. o»
De la distribución de presión de vapor medida
3p j (0 — 0.1) atm
= —45.5 atm/m
dy ,0 = (0.0022 - 0) m
De aquí
- 0 .2 8 8 X 10 4 m /s ( - 4 5 .5 atm/m)
h =■ 0.0164 m/s
(0.1 - 0.02) atm
\
F lujo lam iruir y tu rb u len to
Linea de flujo
turbulenta
Capafe
} amortigm
; ubeapsi
Transición
F l G l HA 6 . 6 D e s a r r o llo d e la c a p a lím ite h id r o d in á m ic a s o b r e u n a p l a c a p la n a .
6 .3 ■ Flujo laminar y turbulento 295
Transición c a lo r h p a r a e l flu jo s o b r e u n a p l a c a p la n a is o té r m ic a
DEPARTAMENTO d e b i b l i o t e c a
Universidad Simón Bolívar - Sede ^ .oral
296 C apítulo 6 ■ Introducción a la convección
para cálculos de la capa límite y, a menos que se señale otra cosa, se utiliza para
cálculos de este texto.
6 .4
Ecuaciones pa ra la transferencia p o r convección
p v + — -(pv) dy
í
v i" j (b
i
! *■ pu + Jj-( p u) d x
I
. ————i F ig u r a 6 .9 V o lu m e n ele c o n tro l
4 d if e r e n c ia l Ulx • (h 1) p a r a
I c o n t e n í'a c ió n d e la tn a s a e n la r a p a
pv lím it e h id r o d in á m u a U id im e n s io n a l.
miento a menudo se denomina culveccwn. Si una esquina del volumen de control se lo-
caliza en ( a , y ) , la velocidad a la que entra la masa al volumen de control a través de la
superficie perpendicular a x se expresa como (pii ) dy, donde p es la densidad total de
masa (p = pA + pb) y u es el componente \ de la velocidad de m asa prom edio. El vo-
lumen de control tiene profundidad unitaria en la dirección z. Como p y u varían con v,
la velocidad a la que la masa sale de la superficie en \ 4- cLx se expresa mediante una
expansión en serie de Taylor de la forma
3 (pu)
(pu) + dx dy
3 (pu) 3 (p v)
{pu) d y + (p v ) d x - pu + dx dy - pv + ^ d y dx = 0
dx
3 (pu) 3 (p v)
" Td x“ + " Idv
T" = 0 (625)
- ( C T y y )d y
V I
I Tx y +
l
(T, I
I
I F k .I 6.10 Esfucrafts
Txy i l ha
d\ n o r m a l y c o r ta n te p a r a un voluni
(v. >•) dx
d e c o n tr o l {dx • d y • 1) c n l a c a -
lím ite h id r o d in á m ic a
°Vy 1>H u n c í sio iid l.
la capa límite, el esfuerzo viscoso (una fuerza por unidad de área) se descompone
dos componentes perpendiculares, que incluyen un esfuerzo norm al cr y un esfue
cortante r (figura 6.10).
Se utiliza una notación de doble subíndice para especificar los componemos
esfuerzo. L1 primer subíndice indica la orientación de la superficie al proporcionar
dirección de su normal hacia afuera, y el segundo señala la dirección del compon
de la fuerza. En consecuencia, para la superficie x de la íigura 6 10, el esfuerzo no
er„ corresponde a un componente de la fuerza normal a la superficie, y el esfuerzo
tante t vv, a una fuerza en la dirección y a lo largo de la superficie. Todos los compo
tes de esfuerzo que se muestran son positivos en el sentido de que tanto la nonnala
superficie com o el componente de la fuerza están en la misma dirección Es decir,
bos están en la dirección coordenada positiva o en la dirección coordenada negaf
Mediante esta convención los esfuerzos viscosos normales son esfuerzos de trac
En cambio, la presión estática se origina a partir de una fuerza externa que actúa so
el fluido en el volumen de control y es. por tanto, un esfuerzo de compresión.
Hay que señalar vanas características del esfuerzo viscoso La fuerza asociada
da entre los elementos contiguos de fluido y es una consecuencia natural del mo
miento del fluido y la viscosidad. Por tanto, se supone que las fuerzas superficiales
la figura 6 10 actúan sobre el fluido dentro del volumen de control y se atribuyen a
interacción con el fluido de los alrededores. Estos esfuerzos desaparecerían si la vel
dad del fluido, o el gradiente de velocidad, se hicieran cero En lo que respecta a los
fuerzos viscosos normales (íjcv y cr y) no deben confundirse con la presión estática,
no desaparece cuando la velocidad es cero
Cada uno de los esfuerzos cambia de forma continua en cada una de las direc
nes coordenadas. Con una expansión en serie de Taylor para los esfuerzos, la fu
superficial neta para cada una de las dos direcciones se expresa como
d a xx , d xy A dx dy
5,X (6.
\ dx dx dv
' d i„ | d a ^ dp'
S. V
dx dy 6.
dx dy dy j
Para aplicar la segunda ley de Newton, también deben evaluarse los flujos de mom
del fluido a través del volumen de control. Si nos concentramos en la dirección.!
flujos relevantes son como se muestra en la figura 6.11. El flujo de masa en cadau-
las dos direcciones hace una contribución al flujo de momento .v total Por ejemplo,
flujo de masa a través de la superficie x (en el plano y-z) es (pu), y el flujo de mo
6 .4 ■ Ecuaciones para la transferencia por convección 299
( pv)it + [ ( p v )u ]d \
i . f,y
I dy
i
j * (p u)tt + - — ■[( p v )u ]d \
I d>'
1 - f T . J
U.y) ( p V)u ^
F lC .I K \ 6 . 1 1 F lu jo s d r m o m en to p a ra u n c o n tro l d e
v o lu m e n d if e r e n c ia l ( dx • d \ • 1) en la c a p a lím ite
lu d ro d in á n iK a b id im e n s io n a l.
3 [(p « )« J 3 Kpv)u]
— r ------ d x(dx) + d x(d x)
3a " dy
Al igualar la velocidad de cambio del momento v del fluido con la suma de las
fuerzas en la dirección .i, obtenemos
3 l(p w )« ] 9[(pü)u] ¿p , . v ,,
------- - f -------------------- = - r ----------~r b — + X ( 6 .2 8 )
dx dy dx dx dy
3u 3u\ 3 . 3rv
p [u ^ + v ^ ) = Y x {(T" ^ p ) + ~ s i + x (629)
/ dv dv\ 3 r^, 3
(630)
No hay que perder de vista la física representada por las ecuaciones 6.29 y 6.30.
Los dos términos en el lado izquierdo de cada ecuación representan la velocidad neta
de flujo de momento del volumen de control. Los términos del lado derecho explican
la fuerza neta de viscosidad debida a la presión, así como la fuerza de cuerpo. Estas
ecuaciones deben satisfacerse en cada punto de la capa límite, y con la ecuación 6.25
se resuelven para el campo de velocidades.
Antes de obtener una solución para las ecuaciones precedentes, es necesario rela-
cionar los esfuerzos viscosos con las otras variables de flujo. Estos esfuerzos están aso-
ciados con la deformación del fluido y son una función de la viscosidad del fluido y de
los gradientes de velocidad. De la figura 6.12 se deduce que un esfu erzo n o rm a l debe
producir una d efo rm a ció n lin e a l del fluido, mientras que un esfu erzo co rta n te produce
una d efo rm a ció n a n gu lar. Además, la magnitud de un esfuerzo es proporcional a la v e -
lo c id a d a la que ocurre la deformación. La velocidad de deformación está relacionada,
a su vez, con la viscosidad del fluido y con los gradientes de velocidad en el flujo. Para
D€PARTAMENTO DE BI&LIOitLA
Universidad Simón Bolívar «tede del Litorf
3l>0 C apitu lo 6 ■ Introducción a la convección
«
V t
/ / Tr,
(7,
_l
I ■<r.
f 3r*Í
(ü ) (h )
3 u 2 / 3u dv\
crX X (6JI
= 2'i ^ ' I M l a J + a7)
do 2 / 3 u dv
°yy = 2 fi 3y 3 \ 3* 3y
/ óu
3u dv \
dv
’xy
= T” = M l a T + ^ j
Al sustituir las ecuaciones 6.31 a 6.33 en las ecuaciones 6.29 y 6.30, las eci
nes de los momentos y y se convierten en
/ 3 u 3 u\ 3p d_ 3 u 2 / 3 u dv
Y & + D ¥ r
dx dx 3 jc 3 y 3 .r 3y
( du dv \
+ + X
dy 37 + 37 j
dv dv 3p d_ ' ■ ^ 2 f du + d v MI
p /u
dy dy dy 3 !Kd x b ) ..
(du 3d \ 1
+ + Y
dx K a 7 + 3x (í
Va V2
£ a d v .* - E adv = pw(e + ~ ^ \ dy ~ ■pul e +
a ■"I >
i v2\
d x dy
d x [p u ( e + 2 )
a ( \l
p u [e + T j dx dy (6.36)
a
r _ r dT a \
^cond, x *^cond, x+ d x ~k k | d x dy
dx dx \ dx J
d ( dT\
(6.37)
£«Jv, >■+ dy
1cond1v + í / v
« V *
W t 1
r dy
i
Eaxti. x * » ^cond. x+dx
^ i
^adv. t ^ adv, x + dx
i
4.
y dx
F ig u r a 6 . 1 3 V o lu m e n d e
DEPARTAMENTO DE BIBLIOTECA
OímAn C.-wU J a I I S+j
atpilulo 6 ■ In tr o d u c c ió n a la co n rvvciém
3 3
WI>C1A = (Xu) d x d y + — [(</„ - p )u ] d x d y + ^ ( r fxu) d x d y (6.,
d_ v2
pv¡e + —
dv
3 / 37*\ d / dT\ 3 d
T -(P V )
dy
d d
+ Td x* (<7~ M + + ay + + 4 = 0
donde t) es la rapidez a la que se genera energía por unidad de volumen. Esta expn
proporciona una forma general del requerimiento de conservación de la energía
capa límite térmica.
Como la ecuación 6.39 representa la conservación de las energías cinética \
ca interna, rara vez se usa en la solución de problemas de transferencia de calos
lugar, se obtiene una forma más conveniente, denominada ecuación de energía t t
c a , multiplicando las ecuaciones 6.29 y 6.30 por ti y v , respectivamente, y restan
resultados de la ecuación 6.39. Después de muchas manipulaciones, se sigue que^
de de ó / dT\ d / dT\ ,d u dv ,
pU d.x + p l d y dx \ d x ) + 3y \ P \d x dy
du dv \2 r/0n\2 / d v \ 2' 2 / du d v \ 2)
3 \ 3.x 3y í
37 + 37 3.x j \ dy
du dv
— + — = 0 (6.45)
áx dy
+ = + W*
+ p$> + q (6.46)
■
‘^A. adv. y +dy
| | MKdit. y +dy
r“ -j dy
I
^ J I * *^A. adv. x +■av
ü
[Link].-r ** |I ! ^ *^A, [Link] +■d.\
I
l J F l C l KA 6 . 1 4 V o lu m e n d e c o n tro l
x. y d i f e r e n c ia l [dx • dy * 1 p a r a la c o n s e r v a c ió n d e
í t* e s [ e c i e s e n la c a | a lim ite d e c o n c e n t r a c i ó n
[Link]. V
MA,adv. v b id im e n s in n a l.
DEPARTAMENTO DE BIBLIOTECA
Universidad Simón Bolívar - Sede u ...urai
301 C apítulo 6 ■ Introducción a la convección
3 ( p Atf)
dx dy
dx
De manera similar, al suponer un fluido incompresible (p constante) y usar la ley
Fick (ecuación 6.21) para evaluar el flujo de difusión, la velocidad neta a la que lai
pecie A ingresa en el volumen de control debido a la difusión en la dirección x es
dpA 3pA
A*a . dif, * - M A. dif, x + d x = —D AB dy -D
dx AB dx ,
3 / ^ ^P a . j d ( dwpA
p \
dy = dx dy (M
dx
Al sustituir de las ecuaciones 6.47 y 6.48, así como de las formas similares parala
rección >, se sigue que
Una forma más útil de esta ecuación se obtiene al expandir los términos del lado
quierdo y sustituir de la ecuación global de continuidad (6.25). Si la densidad de
total p se supone constante, la ecuación 6.50 se reduce a
dCA 3C a ^ 3 / 3C A 3 / 3C a
U— + v (6.
C/X ~ 3 * \ AB d x + ^ K - a T 1 + WA
E jE M IT O 6 . 4
Una de las pocas situaciones en la que es posible obtener soluciones exactas pan,
ecuaciones de transferencia por convección incluye lo que se denomina jlujo par,
En este caso el movimiento del fluido es sólo en una dirección. Considere un c
pccial de flujo paralelo que incluye una placa estacionaria y una móvil de exten
6 . t ■ Ecuaciones para la transferencia por convección 305
finita separadas por una distancia L, con el espacio de en medio lleno con un fluido in-
compresible. bsta situación se denomina flujo de Couctte y ocurre, por ejemplo, en una
chumacera.
1. ¿Cuál es la forma apropiada de la ecuación de continuidad (ecuación 6.25)?
2. Comenzando con la ecuación de momento (ecuación 6.34), determine la distribu-
ción de velocidades entre las placas.
3. Comenzando con la ecuación de energía (ecuación 6.46), determine la distribución
de temperaturas entre las placas.
4. Considere condiciones en las que el fluido es aceite de motor con L = 3 mm. La
velocidad de la placa móvil es U = 10 m/s, y las temperaturas de las placas esta-
cionaria y móvil son T0 = 10°C y TL — 30°C, respectivamente. Calcule el flujo de
calor para cada una de las placas y determine la temperatura máxima en el aceite.
So l u c ió n
E n co n tra r:
1. La forma de la ecuación de continuidad.
2. Distribución de la velocidad.
3. Distribución de temperaturas.
4. Flujos de calor superficiales y temperatura máxima para las condiciones estable-
cidas
E squem a:
Placa
Placa
estacionaria
S u posicion es:
1. Condiciones de estado estable.
2. Flujo bidimensional (sin variaciones en z).
3. Fluido incompresible con propiedades constantes.
4. No hay fuerzas de cuerpo.
5. No hay generación interna de energía
P ro p ied a d es: Tabla A.8, aceite de motor (20°C): p — 888.2 kg/m’, k — 0.145 W/m ♦
K, u = 900 X 10 '6 m2/s, p = v p = 0.799 N • s/m2.
A nálisis:
1. Para un fluido incompresible (p constante) y flujo paralelo (v = 0), la ecuación
6.25 se reduce a
3u
DEPARTAMENTO DE BIBLIOTECA
Universidad Simón Bolívar - Sude leí Litoral
C apítulo 6 ■ Introducción a la convección
»-!4
du
P
d\ j
u(y) = C,v + C2 |
z/(0) = 0 u(L) — U
«('■ )=—t/
L
3. La ecuación de energía (6.46) se simplifica para las condiciones establecí
particular para condiciones de estado estable bidimensionales con v = 0, (elui
y c¡ = 0, se sigue que
dT d
k— +— k — ---
p c r u T x = l h < dx , dy , d y )
Sin embargo, como las placas superior e inferior están a temperatura unifef
campo de temperaturas también debe estar desarrollado por completo, ene-
so ( dTIdx) = 0. Para conductividad térmica constante la forma apropiada
ecuación de energía es, entonces,
d 2r o ,e 2
0 = k ~—
~ +p
dy dy)
La distribución dc temperaturas deseada se obtiene resolviendo esta ecua
reacomodar y sustituir para la distribución de velocidades.
d 2T (d u \2 (U ^2
k i / = ~ p dv U
Al integrar dos veces, obtenemos
6 .4 ■ Ecuaciones para la transferencia por convección 307
C4 - T0 Cl = Z W o + ^ -
y
2k L
y
y Í2 ' v
T (y ) = T0 + - ^ - U z H T l -T 0) <3
2k Z
4. Conociendo la distribución de temperaturas, los flujos de calor de la superficie se
obtienen al aplicar la ley de Fouricr. De aquí
q =
■
—i ,tT
= i
ll\ 2k y
En las superficies inferior y superior, respectivamente, se sigue que
LiU2 k liü 2 k
q ---------- (7/ -7 o )
' 2L L
T, - T0
! £ = J L U 2 (] l _ W + — —= 0
dx 2k 1 L L2
.Vmáx k 2 ( TL - T , ) + X
-
HU
0.145 W/m -K
^'rriáx — (30-10)°C + - L = 0.536L
0.799 Ns / m~ x 100m 2/s2
DEPARTAMENTO DE BIBLIOTtUA
llm f^ m ll' Universidad Simón Bolívar Sede del Litoral
308 C apítulo 6 ■ Introducción a la convección
6*5
Aproximaciones y condiciones especiales
u> v
du du dv dv Capa límite de velocidad
o hidrodinámica
dy d x ' d y ' dx
dT dT Capa límite
dy dx térmica
d e ,a 3C / Capa límite de
dy dx concentración
6 .5 ■ Aproximación*'* > condiciones nsperiales 300
r iV = T =// (6 53)
(hl Bu , , _ ..
— + —=0 (6 5 4 )
<7i dy
Bu Bu 1 dp B2u
" * +t^ =- p * +uó F (6 5 5 )
Adem ás, a partir de un análisis del orden de magnitud que lisa las aproxim aciones
de la capa límite de velocidad [2], se m uestra que la ecuación del momento v (6.35)
se reduce a
~ =0 (6.56)
B\
*
d epa r t a m en t o d e b ib l io t e c a
C apítu lo (» ■ Introducción o la convección
dC a dC ^ _ d 2C f
ti-— + = D AB (6.5
dx ds dy2
Advierta que el último término del lado derecho de la ecuación 6.57 es lo que queda
la disipación viscosa, ecuación 6.41. En la mayor parte de las situaciones, estelen:
se deja de lado en relación con los que explican la adveccion (lado izquierdo de
ecuación) y la conducción ( primer termino del lado derecho) De hecho, es sólo
llujos sónicos o para el movimiento de alta velocidad de aceites lubricantes que U
pación viscosa no se puede dejar de lado.
Las ecuaciones 6 54, 6.55 y 6.58 se resuelven para determinar las variaciones
pacíales de u. v. T y CA en las diferentes capas límite Para un flujo incompresible
propiedades constantes, las ecuaciones 6.54 y 6.55 están desacoplada ,v de 6.57 j 6
Es decir, las ecuaciones 6.54 y 6.55 se resuelven para el cam po de velot ¡dad. u(x, 4
i/(.v, y), para excluir las ecuaciones 6.57 y 6.58. Del conocimiento de z/(.v. y), seev
el gradiente de velocidad (du/<)y)y - 0, y el esfuer/o cortante en la pared se obtiene de
ecuación 6.15 En cambio, a través de la aparición de // y u e n las ecuaciones 6.5’
6.58, la temperatura y la concentración de especies están acoplada\ con el campo
velocidad Por ello deben conocerse las ecuaciones 6.57 v 6.58 antes de que u(x,í
v (x , y ) se resuelvan para T(x. y) y CA(a. y). Una vez que se obtienen T(x%y) y Q
de tales soluciones, se determinan los coeficientes de transferencia de calor
convección y de transferencia de masa a partir de las ecuaciones 6.17 y 6.20. re
vamente. Se sigue que estos coeficientes dependen en gran medida del campo
velocidad.
Como las soluciones de la capa límite por lo general implican matemáticas
allá del alcance de este libro, nuestro tratamiento al respecto se restringirá al añil
del flujo paralelo M^bre una placa plana isotérmica (sección 7.2 y apéndice E) Sin
bargo, en textos avanzados de convección [5-7J se discuten otras soluciones
y se obtienen soluciones detalladas de la capa límite mediante el uso de técnicas
ricas (diferencias finitas o elemento finito) [8].
No sólo desarrollamos las ecuaciones de la capa límite con el propósito de
soluciones. De hecho, nos motivaron principalmente otras dos consideraciones,
ellas es cultivar una apreciación de los diferentes procesos físicos que ocurren
capas limite. Estos procesos, por supuesto, afectarán Ja fricción de la pared, asi
la transferencia de energía y especies en las capas limite Una segunda motiv
surge del hecho de que las ecuaciones sirven para identificar los parámetros el
sim ilitu d de la capa lím ite , asi como a n a lo g ía s importantes entre momento,
transferencia de m asa
6 .6 ■ Similitud de rapas límite 311
6.6
Similitud d e c a p a s lím ite : e c u a c io n e s d e tr a n s fe r e n c ia
por co n vecció n n o r m a liz a d a s
Si examinamos las ecuaciones 6.55, 6.57 y 6 58 con más cuidado, reparamos cn una
fuerte similitud. De hecho, si el gradiente de presión que aparece en la ecuación 6.55 y
el termino de disipación viscosa de la ecuación 6.57 son insignificantes, las tres ecua-
ciones son de la misma forma. Cada ecuación se caí ac te riza por térm inos de advec-
c ion sobre el lado izquierdo y un térm ino de difusión en el lado derecho Esta situación
describe flu jo s de convección fo rza d a de baja velocidad , que se encuentran en muchas
aplicaciones de ingeniería y ocuparán gran parte de nuestra atención en este texto. Es
posible desarrollar las implicaciones de esta similitud de manera racional haciendo pri-
mero adim ensionales las ecuaciones gobernantes.
a * = — y y* = J ( 6 .5 9 )
u v
y U* = — (6 6 0 )
CA - CA s
c a 38 r (6 6 2 )
Las ecuaciones 6.59 a 6.62 se sustituyen en las ecuaciones 6.55. 6.57 y 6.58 para obte-
ner las formas adimensionales de las ecuaciones de conservación que se muestran en la
tabla 6.1. Advierta que no se toma en cuenta la disipación viscosa y que p* = (p /p V 2)
es una presión adimensional. Las condiciones de frontera que se requieren para resol-
ver las ecuaciones también se muestran en la tabla.
De la forma de las ecuaciones 6.63 a 6.65, se infieren tres parám etros de similitud.
Los parámetros de similitud son importantes pues nos permiten aplicar los resultados
obtenidos para una superficie que experimenta un conjunto de condiciones a superfi-
cies geom étricam ente sim ilares que experimentan condiciones por completo diferen-
tes. Estas condiciones varían, por ejemplo, con la naturaleza del fluido, la velocidad del
fluido y/o con el tamaño de la superficie (determinada por L).
d epa r t a m en t o de b ib l io t ec a
Universidad Simón Bol ivar c.’ ríp del Llton
312 C apítulo 6 ■ introducción a la convección
du* du *
Velocidad u + v*
dx* 3yJ
dp* v d2u* u*(x*, 0) = 0 UrxjX*)
+ u*(x*, =c) = Re,
dx* VL dy*2 v*(x*, 0) = 0 V
(6.63) (6.66)
dr* dT* a d2T*
Térmica u*~ -----1- v* ----- = ------------ T*(x*, 0) = 0 T*(x*, oo) = i ReL,¡
3.1 3v* VL dy*
(6.64) (6.67)
Concentración
dc% . 3c*
+ v*
a!c* C*(jc*, 0) = 0 C*(x *. ^) = 1 Re[, Se
3jc* dy* VL dy*2
(6.65) (6 .68 )
VL
Re¡ = ( 6.6
v
Pr = — ( 6.1
a
Sc = (6.1
D AH
Con las ecuaciones 6.69 a 6.71 y las ecuaciones de capa límite, ecuaciones 6 ai
y la inclusión de la forma adimensional de la ecuación de continuidad (6.54). elt
junto completo de ecuaciones de capa límite viene a ser
dit * d\ * „
—-----h ~ — —0
dx * dv *
6 .6 ■ Similitud de copas límite 313
(6.73)
dr*
— + V*~z
¿rr*
=
i d2r*
r-
(6.74)
cbc* Dy* R eLP r <7v*2
(6.75)
Las ecuaciones anteriores son muy útiles desde el punto de vista de que indican cómo
se simplifican y generalizan los resultados importantes de capa limite. La ecuación de
momento (6.73) indica que, aunque las condiciones en la capa límite hidrodinámica
dependen de las propiedades del Huido p y ¡x. la velocidad V y la escala de longitud L,
es posible simplificar esta dependencia agrupando estas variables en la forma del nu-
mero de Reynolds. Por tanto, anticipamos que la solución a la ecuación 6.73 sera de la
forma funcional
(6.76)
2 du*
(6.77)
2
Cf = — f 2{x** R eD (6.78)
KeL
La importancia de este resultado no debe pasarse por alto. La ecuación 6.78 afirma
que el coeficiente de fricción, parámetro adimensional de importancia considerable pa
ra el ingeniero, se expresa exclusivamente en términos de una coordenada espacial adi-
314 C apítulo 6 ■ Introducción a la convección
mensional y del numero de Reynolds Por consiguiente, para una geometría establecida
esperamos que la función que relaciona C j con x* y R eL se aplique universalmente h
decir, esperamos que se aplique a diferentes fluidos y sobre un amplio intervalo de va-
lores de V7 y L.
Resultados similares se obtienen para los coeficientes de convección dc calor yf
transferencia de masa. De manera intuitiva, es posible anticipar que /? depende de i
propiedades del fluido (k, cp, ¡jl , y p), la velocidad del fluido V, la escala de longitud!
y la geometría de la superficie. Sin embargo. la ecuación 6.74 sugiere la manera en|
que se simplifica esta dependencia. I n particular, la solución a esta ecuación se exp
sa en la forma
dp*
T * = f 3 [ x * i y * ,R e L, P r , (67!
dx*
k f (T« - T.) 9 kf dT *
h =
L (T . - T J dy* L dy* v*=0
N ú m e ro tle N usselt:
hL dr *
Nu ------- +
ch *
Nu = = f<i(R e L, P t )
kf
De manera similar, se argumenta que, para la transferencia de masa en uní
gas sobre un liquido que se evapora o un sólido que se sublima, el coeficientedet
6 .6 ■ Similitud de capas límite 315
dp*
CX = fjx*.y*, R eL, (6.83)
dx*
S ¡ ,s '2 s L = 3 C \ (6.84)
D AB 3\ *
v*=0
Sh = = fx (R e L ,5 0 (6 86)
^AB
Del desarrollo anterior obtuvimos los parámetros adimensionales relevantes para
capas límite de convección forzada y de baja velocidad. Lo realizamos al expresar en
forma adimensional las ecuaciones diferenciales que describen los procesos físicos
dentro de las capas limite. Un enfoque alternativo incluiría el uso del análisis dimen-
sional en la forma del teorema pi de Buckingham [9|. Sin embargo, el éxito de este
método depende de la habilidad para seleccionar, principalmente de la intuición, los di-
versos parámetros que influyen en un problema. Por ejemplo, al conocer de antemano
que h = j{ k , cp, p, p., V , L), se utiliza el teorema pi de Buckingham para obtener la
ecuación 6.82. Sin embargo, al comenzar con la forma diferencial de las ecuaciones de
conservación, eliminamos el trabajo de adivinar y establecimos los parámetros de simi-
litud de forma rigurosa.
El valor de una expresión como la ecuación 6.82 debe apreciarse por completo
Establece que los resultados de la transferencia de calor por convección, obtenidos teó-
rica o experimentalmente, se representan en términos de tres grupos adimensionales,
DEPARTAMENTO DE BIBLIOTECA
Un v/„is i Uuü oriiion oonvur - Sede <
C apítulo 6 ■ Introducción a la convección
Ej e m pl o 6 .5
Las pruebas experimentales sobre una parte del álabe de turbina que se muestra indic*^
un flujo de calor hacia la hoja de q" = 95,000 W/m2. Para mantener una tempera
superficial en estado estable de 800 C, se elimina el calor que se transfiere al álabe
ciendo circular un fluido refrigerante dentro del mismo.
q" = 95 kW/m2
Condiciones
originales
Son p ió n
E n c o n tr a r :
1. Flujo de calor hacia el alabe cuando se reduce la temperatura de la superficie.
2. Flujo de calor hacia un álabe de turbina mas larga de la misma forma [Link]
dad de aire reducida.
E squem a:
— L 40 mm
Caso 1 Caso 2 V
6 .6 ■ Similitud de capas límite 317
S u p o s ic io n e s :
1. Condiciones de estado estable
2. Propiedades del aire constantes.
A n á lisis:
1. De la ecuación 6.81 se sigue que, para la geometría establecida,
hL
N u= — = f 4(x*, R eL, P r )
k
donde
_ q"
H'~ ( T „
De aquí
2. Para determinar el flujo de calor asociado con el álabe más largo y el flujo de aire
reducido (caso 2), advertimos primero que, aunque L aumenta por un factor de 2,
la velocidad disminuye por el mismo factor y el número de Reynolds no cambia.
Es decir,
v 2l 2 vl
R e 1 ^2 ~ R ^l
N u 2 —N u
h 2L 2 hL . L q" L
-L -L = o h2 = h— = — ------— —
k k L2 (Too T9) L2
0.04 m .
A = 95,000 W/m2 X — — = 47,500 W/nr <1
*1 0.08 m
DEPARTAMENTO DE BlBLlOítOA
Universidad Simón Bol iva- Bode dol Utora'
318 C a p ítu lo 6 ■ introducción a la convección
6 .7
Sifpiijti'Ctdo Jísiro de los parám etros adimensionales
F ,a p y ^ L = pyL =
Fx pV / L ¡j
Esperamos entonces que las fuerzas de inercia dominen para valores grandes de
que las fuerzas viscosas dominen para R e pequeños
Hay varias consecuencias importantes de este resultado. Recuerde que el nú
de Reynolds determina la existencia de flujo laminar o turbulento. En cualquier
existen pequeñas perturbaciones que se pueden amplificar para producir condic
turbulentas. Sin embargo, para R e pequeños, las fuerzas viscosas son suficiente
grandes con relación a las fuerzas de inercia para evitar esta amplificación. Por ello
mantiene el flujo laminar. Pero, al aumentar Re, los efectos viscosos se hacen
importantes de manera progresiva en relación con los efectos de inercia, y las peq
perturbaciones se amplifican a un punto en el que ocurre la transición. Debemos
rar también que la magnitud del número de Reynolds influya en el espesor 8 de la
límite hidrodinámica. Al aumentar Re en una posición fija sobre una superficie e
mos que las fuerzas viscosas se vuelvan menos influyentes en relación con las fue
de inercia. Por ello, los efectos de la viscosidad no penetran tan lejos en el flujolí
el valor de 8 disminuye.
La interpretación física del número de Prandtl se sigue de su definición como
razón de la difusividad del momento u a la difusividad térmica a. El numero de
proporciona una m edida de la efecto idad relativa del transporte de momento ve
pot difusión en las capas lím ite hidrodinám ica y térm ica . respectivamente Déla
A.4 vemos que el número de Prandtl de los gases es cercano a la unidad, cn cuyo
la transferencia de energía y momento por difusión son comparables. En un mcüi
quido (tabla A.7). P r < 1 y la velocidad de difusión de energía excede granéeme
velocidad de difusión de momento. Lo opuesto es cierto para aceites (tabla A.5l
6 .7 ■ Significado físico de los parámetros adimensionales 319
los que P r > 1 . De esta interpretación se sigue que el valor de P r influye fuertemente
en el crecimiento relativo de las capas límite hidrodinámica y térmica. De hecho, para
capas límite laminares (en las que el transporte por difusión no se oscurece por la mez-
cla turbulenta), es razonable esperar que
8
— « P rn (6.87)
O/
8
— « S cn (6.88)
Otro parámetro, que está relacionado con P r y 5c. es el núm ero de Lew is (Le). Se defi-
ne como
a Se
Le" B = Tr <6 í i 9 >
5,
— « L en (6.90)
uc
Así. el número de Lewis es una medida de los espesores relativos de las capas límite
térmica y de concentración. Para la mayor parte de las aplicaciones es razonable supo-
ner un valor de n = 1/3 en las ecuaciones 6.87, 6.88 y 6.90.
La tabla 6.2 enumera los grupos adimensionales que aparecen con frecuencia en
los textos sobre transferencia de calor y masa. La lista incluye grupos ya considerados,
así como los que se introducirán para condiciones especiales. Conforme se enfrente a
un grupo nuevo, apréndase de memoria su definición e interpretación. Advierta que el
num ero de G ra sh o f proporciona una medida de la razón de las fuerzas de empuje a las
fuerzas viscosas cn la capa límite hidrodinámica. Su papel en la convección libre (capí-
tulo 9) es, con mucho, el mismo que tiene el número de Reynolds en la convección for-
zada. LI núm ero de Eckert proporciona una medida de la energía cinética del flujo en
relación con la diferencia de entalpias a través de la capa límite térmica. Juega un pa-
pel importante en flujos de alta velocidad para los que la disipación viscosa es signifi-
cativa. Tenga en cuenta también que. aunque similares en forma, los números de
Nusselt y Biot difieren en definición e interpretación. Mientras que el número de Nus-
selt se define cn términos de la conductividad térmica del fluido, el número de Biot se
basa en la conductividad térmica del sólido, ecuación 5.9.
320 C a p ítu lo 6 ■ Introducción a ln convección
Número de Schmidt
V Razón de las difusividades de momento y
(Se) Da b de masa.
Número de Shcrwotxi KL Gradiente de concentración adimensional en
(Sin) D AB la superficie.
Número de Stanton h Nul Numero de Nusselt modificado.
(St) pVcp R eL P r
6 .8 ■ A n a lo g ía s d e la c a p a lím ite 321
T a b la 6 .2 C o n tin u a c ió n
<U
\nalogias d e la c o p a lím ite
DEPARTAMENTO DE BIBLIOTECA
Universidad Simón Solivar - Sede o. a 1
322 Capítulo 6 ■ In tro d u c c ió n a la c o n ve c c ió n
TABLA 6 . 3 R elacion es funcionales relacionadas con las analogías de las capas lím ite
dp* dp * \ / dp*\
u* = / , **, v*, R e ,, T* = f J x*, y * , Re, , Pr, C% = / 6Íx*, y*, R e „ Se, — J
dx* dx* )
(6.76) (6.79) (6.83)
2 du* hl dT*
\r. __ — i. 5 / j = hmL = + ^C *
Ct - N k d y* DAB ay [
Ret dy* v* = 0 y * ” 0
La analogía también es útil para relacionar de forma directa los dos coeficientes
convección. En los capítulos siguientes encontraremos que N u y Sh por lo general s
proporcionales a Pr11 y Se", respectivamente, donde n es un exponente positivo menor
que 1. Para anticipar esta dependencia, usamos las ecuaciones 6.81 y 6.85 para obtener
Nu = f 4 ( a * . ReL)Prn Sh = f 1 (x * .R e L)Scn
en cuyo caso
Nu * ' Sh
7T = / 4(-* ,R e L ) = f 1( x =— (6.911
Pr
h L /k hmL /D AB
p rn Scn
o, dc la ecuación 6.89.
h \-n
p e ¡} L e (6.92)
0 M Len
Este resultado a menudo sirve para determinar un coeficiente de convección, por ejem-
plo, hm, a partir del conocimiento del otro coeficiente. La misma relación se aplical
los coeficientes promedio h y h m, y sirve en el flujo turbulento, así como cn el lamin*
Para la mayor parte de las aplicaciones es razonable suponer un valor de n =
Ej e m pl o 6 .6
Un sólido de forma arbitraria se suspende cn aire atmosférico que tiene una temperan
ra de flujo libre y velocidad de 20°C y 100 m/s, respectivamente. El sólido tiene j
longitud característica de 1 m, y la superficie se mantiene a 80°C. En estas concL
nes. las mediciones del flujo de calor en un punto particular Cv*) sobre la superi*¡
y de la tem peratura en la capa límite sobre este punto ( a * , y * ) revelan valores dt
6 .8 ■ Analogías de la capa límite 323
So l u c ió n
E squem a:
\ / ~ C h {x*,y*)
S u p o s ic io n e s :
1. Comportamiento de capa límite incompresible bidimensional de estado estable;
propiedades constantes.
2. Las aproximaciones de capa límite son válidas.
3. Disipación viscosa insignificante.
4. La fracción molar del vapor de agua en la capa limite de concentración es mucho
menor que la unidad.
T -T s ( dp*
DEPARTAMENTO DE BIBLIOTECA
Universidad Simón Bolívar - Sede del Litoral
C a p itu la 6 ■ introducción a la conrección
Ca Ca .s
C*A~
A
^
Ct z ^ r hix*'y*-ReL'Sc-^)
Sin embargo, para el caso 1
„ „ P a .sat 0 .0 8 2 kg/m 3 3
CA ,j= [Link],(S° C) = — ----- = — — = 0.0046 kmol/m
11A 18 kg/kmol
se sigue que
CA(**, >’*) = 0.67(0.0046 kmol/m3) = 0.0031 kmol/m3
El flujo molar se obtiene de la ecuación 6.7
o
N'Á(x*) = 3.54 X 10 4 kmol/s • m2
6 .8 ■ Analogías de la capa límite 325
donde <y"vap se aproxima como el producto del flujo de masa evaporativo y el calor la-
tente de vaporización
h ( 7'os - Ts) = V U P a . s a t( E ) - P a .o c I (6 .9 5 )
donde la densidad de vapor en la superficie es la que se asocia con las condiciones sa-
turadas en T . Por ello, la magnitud del efecto de enfriamiento se expresa como
(h
Ts hfg J [ P a , sai (^Ty ) P a .<» (6 .9 6 )
Flujo de gas
(especie B)
Q onv <7evap Interfaz
gas-liquido
/ /
4
15^.3 = = -íCapa líquida M "tfagr F ig u r a 6 .1 5
^ (especia A) rJ| Jl
=I I n t e r c a m b io d e c a l o r la t e n te y
J
s e n s i b l e e n u n a in t e r fa z g a s - l íq u i d o .
DEPARTAMENTO DE BIBLIOTECA
Un vwrsidad Simón Bolívar - Sede oral
C apítulo 6 ■ Introducción a ln convección
P a . s«i(7í) Pa
(T x ~ Ts) 2/3 (6.
cP Le
Ej e m pl o 6 .7
Un recipiente, que se envuelve en una tela humedecida de forma continua con uní
do altamente volátil, se utiliza para conservar bebidas frías en regiones ár das calien
Suponga que el recipiente se coloca en aire ambiental seco a 4ܰC, y que la transfe
cia de calor y masa entre el agente humedecedor y el aire ocurre por convección fr
da Se sabe que el agente humedecedor tiene un peso molecular de 200 kg/mol \
calor latente de vaporización de 100 kJ/kg Su presión de vapor saturado para las
diciones que se establecen es aproximadamente 5000 N /irr, y el coeficiente de difi
del vapor en aire es 0.2 X 10-4 m2/s. ¿Cuál es la temperatura de estado estable de
bebida?
S o l u c ió n
E squem a:
S u p o s ic io n e s :
1. La analogía de transferencia de calor y masa es aplicable.
2. El vapor muestra un comportamiento de gas ideal.
3. Los efectos de radiación son insignificantes.
4. Las propiedades del aire se evalúan en una temperatura media de la capa limite
que se supone a 300 K
h fg P a . sat(^í) PA
(7U - Tf ) = 2/3
t^lpc L e
t 2
s - TccT; + b = o
donde el coeficiente B es
^A^fgPA. sat
B =
t^lpCp L e 213
o
B = [200 kg/kmol X 100 kJ/kg X 5000 N/m 2 X 1 0 '3 kJ/N • m]
22.5 x 10 m /s\2/3
X = 9514 K 2
20 X 1(T6 m2/s
De aquí
r ± V r i - 4B 313 K ± V (3 1 3 )2 - 4(9514) K
T =
Al rechazar el signo de menos sobre bases físicas ( Ts debe ser igual a T «, si no hay eva-
poración, en cuyo caso P A sa, = 0 y B = 0), se sigue que
6 .8 .3 Analogía de Reynolds
Es posible obtener una segunda analogía de capa límite al observar en la tabla 6 1 que,
para d p * /d \* = 0 y P r = Se = 1, las ecuaciones de conservación, ecuaciones 6.63 a
6 65, son precisamente de la misma forma. Ademas, como «oc = V si dp*/cLx* = 0, las
condiciones de frontera, ecuaciones 6.66 a 6.68, también tienen igual forma Por consi-
d epa r t a m en t o de b ib l io t ec a
328 C apítulo 6 ■ Introducción a la convección
guíente, las soluciones para //*, T* y deben ser equivalentes. Es decir, de las ecu®
ciones 6.76, 6.79 y 6.83 de la tabla 6.3, f \ = / 3 = /(,. Además, el coeficiente de fricción,
el numero de Nusselt y el numero de Sherwood están relacionados por el requisito que
f 2 — / a = f , y, de las ecuaciones 6.78, 6.81 y 6.85, concluimos que
Rp
C f — - = N u = Sh (6.99}
f l
Nu
St = 6.101
p V cp Re Pr
h.m Sh
S tm = 6 . 18!
V R e Se
La ecuación 6.99 también se expresa en la forma
C\
'/ St = St m
C
—= S t P r 1 3 =
~ Jh 0.6 < P r < 60 ( 6.1
C,
= S tm S e 2 3 = Jm 0.6 < Se < 3000 (6.
donde j H y j m son los fa cto res j de Colburn para transferencia de calor y de mnsíjj
pectivamente Para el flujo laminar las ecuaciones 6.103 y 6.104 sólo son apra
cuando dp*ldx* ~ 0, pero en el flujo turbulento las condiciones son menos sensihl
efecto de los gradientes de presión y estas ecuaciones siguen siendo aproximad
válidas. Si la analogía es útil para cualquier punto sobre una superficie, se aplican
coeficientes promedio de la superficie.
6 .»
Efectos de la turbulencia
Tiempo, i
F l í . l HA C». 1 6 V a r ia c ió n d e u n a p ro | e d a d c o n e l tie m p u e n
a lq u il p u n to e n u n a c a p a im ite t u r b u le n ta .
ñas perturbaciones asociadas con distorsiones cn las líneas de fluido de un flujo laminar
finalmente conducen a condiciones turbulentas Estas perturbaciones se pueden originar
desde el flujo libre, o se inducen por la aspereza de la superficie. El comienzo de la tur-
bulencia depende de si estas perturbaciones se amplían o se atenúan en la dirección de
flujo del fluido, que a su vez depende de la razón de la fuerza de inercia a la viscosa (nú-
mero de Reynolds). Recuerde que si el número de Reynolds es pequeño, las fuerzas
de inercia son pequeñas en relación con las fuerzas viscosas. Las perturbaciones que
ocurren de manera natural se disipan entonces, y el flujo permanece como laminar. Sin
embargo, para un número de Reynolds grande, las fuerzas de inercia son suficientemen-
te grandes para amplificar las perturbaciones, y ocurre una transición a la turbulencia.
Fn la sección 6.3 observamos que el numero de Reynolds critico. Re , que se requiere
para la transición es aproximadamente 5 X 10‘ para un flujo sobre una placa plana.
La turbulencia se asocia con la existencia de fluctuaciones aleatorias cn el fluido
y, al menos en pequeña escala, el flujo es inherentemente inestable Este comporta-
miento se muestra en la figuia 6.16, donde la variación sobre una propiedad del flujo
arbitraria P se traza como función del tiempo en alguna posición en una capa límite
turbulenta La propiedad P es un componente de la velocidad, la temperatura del fluido
o una concentración de especies, y en cualquier instante se representa como la suma
de un valor m edio respecto al tiem po P y un componente de fluctuación P '. El prome-
dio se toma sobre un tiempo grande comparado con el periodo de una fluctuación típi-
ca, y si P es independiente del tiempo, se dice que el flujo medio respecto al tiempo es
estable.
La existencia de flujo turbulento será ventajosa en el sentido de que proporciona
transferencia de calor y de masa aumentadas. Sin embargo, el movimiento es en extre-
mo complicado y difícil de describir de forma teórica Aunque las ecuaciones de capa
límite que se desarrollaron en las secciones anteriores son aun aplicables, las variables
dependientes (u, v, T, CA) deben interpretarse como valores instantáneos , y es imposi
ble predecir su variación exacta en el tiempo. En sentido práctico, sin embargo, esa in-
capacidad para determina! la variación de las propiedades instantáneas P con el tiempo
no es una restricción seria, pues el ingeniero por lo general se ocupa sólo de las propie-
dades medias respecto al tiempo P . Las ecuaciones de la forma P = P + P ' se sustitu-
yen por cada una de las variables de flujo en las ecuaciones de capa límite Para un
flujo de propiedades constantes, incompresible y estable, con los procedimientos de
promedio temporal ya establecidos [2, 12], se obtienen las siguientes formas de las
ecuaciones del momento v, energía y conservación de especies:
/ du du\ f du ~ r~ Á
(6.105)
r a í + Vd y ) ~ dx dyi » d y - p U V )
DEPARTAMENTO DE BIBLIOTECA
Universidad Simón Bolívar Sede do dtora
330 Capitulo 6 ■ i nlrotlncción n la convección
st ar\ a
" 3 7 + e a7 (6.106
dCA ac,
U + l>- D AB - u 'C ' (6.l07i
a.x c)v 3v av
Las ecuaciones son como las de la capa límite laminar, excepto por la presencia de té
minos adicionales de la forma a 'b '. Estos términos explican el efecto de las fluctuacio*
nes de turbulencia sobre el transporte de momento, energía y especies.
Sobre la base de los resultados precedentes, es normal hablar de un esfuerzo er
tante total y de flujos totales, que se definen como
dü
Ttot IM" pu ' v ' (6.1
dy
3f
q’L = ~ pcpV'T' (61
ac,
N"
i v A tot
- - n
i y AB - v 'C (6.11
dy
y consisten en contribuciones debidas a difusión molecular y mezcla turbulenta. A
tir de la forma dc estas ecuaciones vemos cómo las transferencias de momento, en^
y especies aumentan por la existencia de la turbulencia. El término p u 'v , que ap-
en la ecuación 6.108, representa el llujo de momento debido a las fluctuaciones tur
lentas y a menudo se le denomina esfuerzo Je Reynolds.
Un modelo conceptual sencillo atribuye el transporte de momento, calor y masa
una capa limite turbulenta al movimiento de rem olinos, pequeñas porciones de fi-
en la capa límite que se mueven por un tiempo corto antes de perder su identidad,
bido a este movimiento, el transporte de momento, energía y especies aumenta mu
La noción de transporte por remolinos implica la introducción de un coeficiente
transporte que se define como la difusividad parásita para la transferencia de mo-
to e,w, que tiene la forma
du
pu ' v ' (61
PEm 3 7 ^
Por eso el esfuerzo cortante total se expresa como
du
= p ( v + e M) (611
dy
De manera similar, se definen las difuso idades parásitas para la transferencia ¿
e h y e m mediante las relaciones
lor y de m asa
dT ______
e „ — = —v T ' ( 6.1
H dy
dC_A
- - v 'C 'A (61
dy
en cuyo caso
dT
tfíó. = ~ p c p( a + e„) ( 6.1
dy
6 .1 0 ■ Coeficientes de convección
ÓU <
dV >■= 0, lam dv \ = O turb
=0
dCA
N a toi = ~ ( d a b + e rn) “3 — (6 .1 1 6 )
d epa r t a m en t o de b ib l io t ec a
Universidad Simón Bolívar <• iv>r'
332 C apítulo 6 ■ Introducción a ln convección
situaciones de flujo simple que tales soluciones se llevan a cabo fácilmente. 1:1 méttxl
más práctico a menudo implica el cálculo de h y hm a partir de relaciones empíricas dr
la forma dada por las ecuaciones 6 81 y 6 85. La forma particular de estas ecuador
se obtiene c orrelac tonando resultados de mediciones de transferencia de calor i
masa por convección en términos de grupos adimensionales apropiados. Este es ele
foque en el que se hace énfasis en los capítulos siguientes
B ib lio g r a fía
1 Webb, R L., Int. Cornrn. H eat M ass T r a n s 17, 7. Kaviani, M , P rincipies o f Convectiva
529. 1990. Transfer, Springer-Verlag, Nueva York. 1
2. Schlichting, H., B oundary Lciyer T hcory , 7a ed., 8 Patankar, S. V., N um érica! Heat Tra
McGraw Hill, Nueva York, 1979. F luid Flcm\ Hemisphere Publishing, Nu
3. Bird, R. B., W. E. Stewart y E. N Lightfoot, 1980.
T ra n sp o n P heiw m ena, caps. 10 y 18, Wiley, Nue- 9 Fox, R W. y A T. McDonald. Introd
va York, 1966. F luid M echan i es. Wiley, Nueva York. 1
4. Hartnett, J. P., “Mass Transfer Cooling”, en W. M. 10. Colburn, A. P , Trans. Am. Inst Chent.
Rohsenow y J. P. Hartnett. Editores, H andhook o f 174, 1933.
H eat Transfer , McGravv-Hill. Nueva York, 1973. 11 Chilton, T H y A. P Colburn, huí. E
5. Kays, W M. y M E. Cravvtord, C onvective H eat 26,1183,1934
a nd M ass Transfer , McGraw-Hill, Nueva York, 12. Hinze, J. O., Turhulence. 2a ed.. Me
1980. Nueva York, 1975.
6. Burmeister, L. C., Convectivo H eat Transfer , Wi-
ley- Interscience, Nueva York, 1983.
■ Problemas 333
Problemas
(.««‘íic irn te - «l«* IraiisiVriMK'ia «le c a lo r están correlacionadas por una expresión de la forma
T(°C) = 20 + 70 exp( —600xx). donde v y y están cn
6.1 Se sabe que el coeficiente local de transferencia de ca-
metros Determine y elabore una gráfica de la forma
lor h,. para el flujo laminar sobre una placa plana varía
en la que varia el coeficiente de convección local h con
av donde x es la distancia medida desde el inicio (x =
x Evalué el coeficiente de convección promedio h para
0 ) de la placa ¿Cuál es la razón del coeiiciente prome
la placa.
d o entre el inicio y algún lugar x sobre la placa al coe-
ficiente local en .v? 6 .6 La transferencia de calor por unidad de anchura (nor-
mal a la página) desde una sección longitudinal, xs — *i.
6.: Para la convección laminar libre de una superficie ver- se expresa com o q \ 2 = /»i2(x2 — .V |)(T , — T » ) . donde
tical caliente, el coeficiente de convección local se ex
/i 12 es el coeficiente promedio para la sección de lon-
presa como J¡x = C v~‘ , donde h x es el coeficiente en la
gitud (x2 — Vi) Con idere un flujo laminar sobre una
distancia x desde el inicio de la superficie y la cantidad
placa plana con una temperatura uniforme / , . La va-
( que depende de las propiedades del fluido, es inde-
riación espacial del coeficiente de convección local
pendiente de .v. Obtenga una expresión para la razón
es de la forma h x = C x 12. donde C es una constante.
/ //q. donde h es el coeficiente promedio entre el micio
i 0) y la posición x Dibuje la variación de h x y h x
con x.
DEPARTAMENTO DE BIBLIOTECA
Universidad «Jil itJll 4* - • «JH
331 C apítulo 6 ■ Intnulucción a la convección
6.8 Un procedimiento experimental para validar los resulta- donde y es la distancia normal a la superficie} el
dos del problema anterior implica precalentar un diseo mero de Prandtl, P r = cr pJk = 0.7. es una pro
de cobre a una temperatura inicial T elevada y registrar adimensional del fluido. Si Tx = 400K. 7, = 3{X)
su historia de temperaturas T(t) conforme se enfría me- u j v = 5000 m" . ¿cuál es el flujo de calor poru
diante el flujo que choca a una temperatura final 7 La de arca en la superficie ?
disminución de temperatura medida se compara enton-
ces con las predicciones basadas en la correlación para T ra n sició n d e la c a p a lím ite
NuD. Suponga que valores de a = 0.30 y n = 2 se aso-
6.12 Considere un flujo de aire sobre una placa plana
cian con la correlación. longitud L = 1 m en condiciones para las queoc
Cons dere las condiciones experimentales para las transición en x = 0.5 m con base en el nur
que un disco de diámetro D = 50 mm y longitud L = Reynolds crítico. R e v c = 5 X I05. En las regicr
25 mm se precalienta a 7, = 1000 K y se enfría a Tf - minar y turbulenta, los coeficientes de convección
400 K mediante un flujo de aire que choca a 7* = 300 K. son. respectivamente.
La superficie enfriada del disco tiene una cm sividad
-0..S
de e = 0 . 8 y se expone a alrededores isotérmicos leja- ^lamCO ^turb Cur^V
nos para los que 7j)ll — T*. Las otras superficies del
donde CIam = 8.845 W/m2 • K ° \ CHlrb = 49.75 Wf
disco están bien aisladas, y la trasferencia de calor a
K0 8. y a tiene unidades de m.
través de la varilla de apoyo se considera insignifican-
(a) Mediante la evaluación de las propiedades te
te. Con los resultados del problema anterior, calcule y
sicas del aire a 350 K. determine la veloci
trace las historias de temperatura que corresponden a
flujo de aire.
velocidades del aire de V7 = 4, 20 y 50 m/s. Se suponen
propiedades constantes para el cobre (p = 8933 k g/m \ (b) Desarrolle una expresión para el coeficiente de
cp = 425 J/kg • K, k = 386 W/m • K) y para el aire vccción promedio /íjam(.v), como función de
(u = 38 8 X 10 6 nr/s, k = 0.0407 W m • K P r = tancia desde el inicio de la placa v. para la
0.684). laminar. 0 < x ^ x L.
(c) Desarrolle una expresión para el coeficiente
P e rfile s de la r a p a lim ite vccción promedio, h tuIb(.\). como función de
tanda desde el inicio de la placa, v, para a
6.9 Ln un flujo sobre una superficie, los perfiles de veloci- turbulenta. xc < x ^ L.
dad y temperatura son de la forma (d) En las mismas coordenadas, trace los c
u(y) = Ay + B y2 - Cv3 y de convección local y promedio. hÁ. ) hr
T(y) = D + Ey + F y2 - G y 3 vamente, como función de x para 0 < x < ¿
donde los coeficientes A a G son constantes. Obtenga 6.13 Un ventilador que proporciona velocidades de1
expresiones para el coeficiente de fricción C f y el coefi- ta de 50 m/s se utilizará en un túnel de viento
ciente de convección /? en términos de mx, 7* v los velocidad con aire atmosférico a 25°C. Si sed
Problemas 333
el túnel de viento para estudiar el comportamiento de que las placas están separadas por agua y otro en que
capa límite de una placa plana hasta números de Rey- las placas están separadas por aire.
nolds R cx = 10 , ¿cual es la longitud de placa mínima (a) Para cada uno de los dos fluidos, ¿cuál es la fuerza
que debe utilizarse? ¿A que distancia desde el inicio de por unidad de arca superficial que se requiere para
la placa ocurriría la transición si el numero de Reynolds mantener la condición anterior? ¿Cuál es el reque-
crítico fuera R e ( = 5 X 10S? rimiento de potencia correspondiente?
6.14 Suponiendo un número de Reynolds de transición de (b) ¿C ual es la disipación viscosa asociada con cada
5 X 10"', determine la distancia desde el inicio de una uno de los dos fluidos?
placa plana a la que ocurrirá la transición para cada (c) ¿Cuál es la temperatura máxima en cada uno de los
uno de los siguientes fluidos cuando ux = 1 m/s: aire dos fluidos?
atmosférico, aceite de motor y mercurio En cada caso
6.19 Se hace un juicio con respecto a la influencia de la disi-
la temperatura del fluido es 27°C.
pación viscosa en la transferencia de calor por convec-
ción forzada mediante el cálculo de la cantidad P r Ec.
Ecuaciones d e c o n s e r v a c ió n y s o lu c io n e s donde el número de Prandtl Pi = c mIk y el numero de
6.15 Considere el volumen de control que se muestra para el Eckert Ec = LPlcyDT son grupos adim cnsionales. La
caso especial de condiciones de estado estable con v = velocidad característica y la diferencia de temperaturas
0. T - T(y) y p es constante. del problema se designan como U y A7. respectiva-
mente. Si P r E c < I. no se toman en cuenta los efectos
de disipación. Considere el flujo de Couette para el que
t + dy una placa se mueve a 1 0 m/s y se mant ene una diferen-
dy • cia de temperaturas de 25°C entre las placas Mediante
I(l\ ■*— f>p ,
la evaluación de las propiedades a 27°C, determine el
.J l- l
X, U valor de P r Ec para aire agua, y aceite de motor. ¿Cuál
es el valor de Pi be para aire si la placa se mueve a la
velocidad del sonido?
(a) Pruebe que u = u(y) si v = 0 en cualquier lugar.
6.20 Considere el flujo de Couette para el que la placa móvil
(b) Derive la ecuación del momento en v y simplifí- se mantiene a una temperatura uniforme y la placa es-
quela tanto como sea posible. tacionaria esta aislada. Determine la temperatura de la
(c) Derive la ecuación de energía y simplifíqucla tanto placa aislada, exprese su resultado en términos de las
como sea posible. propiedades del fluido y la temperatura y velocidad de
la placa móvil. Obtenga una expresión para el flujo
6.16 Considere una chumacera ligeramente cargada que usa
de calor por unidad de área en la placa móvil
aceite y tiene las propiedades constantes m = 1 0 - 2 kg/s ■
m y k = 0.15 W m • K. Si el cojinete y la chumacera se 6.21 Considere el flujo de Couette con transferencia de calor
mantienen cada uno a una temperatura de 40°C, ¿cuál para el cual la placa inferior (placa móvil) se mueve
es la temperatura máxima en el aceite cuando el cojine- con una velocidad de U = 5 m/s y está perfectamente
te gira a 1 0 m/s? aislada. La placa superior (placa estacionaria) es esta-
cionaria y esta construida de un material con conducti-
6.17 Considere una chumacera ligeramente cargada que usa
vidad térmica = 1.5 W/m • K y espesor I.]K = 3 mm.
aceite y tiene las propiedades constantes p = 800
La superficie externa se mantiene a T e = 40°C Las
kg/m\ v = 10 5 m2/s, y k = 0.13 Wm • K. El diámetro
placas están separadas por una distancia L„ = 5 mm,
de la chumacera es 5 mm, el espacio es 0.25 mm y el
que se llena con un aceite de motor de viscosidad m =
cojinete opera a 3600 rpm
0.799 N • s/m 2 y conductividad térmica kG — 0 145
ta) Determine la distribución de temperaturas en la pe
W/m • K.
lícula de aceite, suponiendo que no hay transferen-
cia de calor en la chumacera y que la superficie del -pe
d epa r t a m en t o de b ib l io t ec a
Universidad Simón Bolívar Seded*' oral
336 C apítulo 6 ■ Introducción a la convección
(a) En coordenadas Y (\)—y, dibuje la distribución de tem- (d) Genere una gráfica de 6 contra 17 para 0 ^ 17 < |
peraturas en la película de aceite y la placa móvil valores de P r E c = 0, 1, 2, 4, 10 Explique lasca
(b) Obtenga una expresión para la temperatura en la ractensticas clave de las distribuciones de tempe
superficie inferior de la película de aceite, 7(0) = ratura
Tth en términos de la velocidad de la placa U , los
6.24 Considere el problema de un flujo laminar incompresi
parámetros de la placa estacionaria (7 pc, k^., L ^ ) y
ble estable entre dos placas paralelas infinitas estado
los parámetros del aceite (/x . k , l. ). Calcule esta
nanas que se mantienen a diferentes temperaturas
temperatura para las condiciones establecidas.
r 71 Placa
6.22 Un eje con un diámetro de 100 mm gira a 9000 rpm en .
paralelas
-=---- »»r"*!raJ—■
una chumacera de 70 mm de longitud. Un hueco uni infinita
torme con lubricante separa el eje y la chumacera Las
propiedades del lubricante son ¡x = 0 03 N • s m' y £ j x <G
k — 0.15 W/m • K, mientras que el material del cojine rr-
te tiene una conductividad térmica de kc = 45 Wm • K.
m
Cojinete, k
T*C Denominado flu jo d e P oiseuille con transferencia tfcj
calor, este caso especial de flujo paralelo es uno para<
Cojinete, kt Lubricante que el componente .1 de la velocidad es finita, perol
Te componentes y y 2 (i/y w) son cero.
Eje (a) ¿Cuál es la forma de la ecuación de continui
Eje 100 Lubricante para este caso? ¿En qué dirección el flujo se
de cuentra com p leta m en te d e sarrollad o ?
diámetro
Superficie enfriada con agua, (b) ¿Que formas toman las ecuaciones de momento.
T , = 30 C \ y y* 6Cuál es la forma del perfil de velocidad!;
Observe que, a diferencia del flujo de Couena.
movimiento del fluido entre las placas ahora se
tiene mediante un gradiente de presión finito. L
(a) Dcteim ne la disipación viscosa. /x<t>( W/m ). en el es este gradiente de presión en relación con la v
lubricante dad máxima del fluido?
(b) Determine la transferencia de calor (W ) del lubri- (c) Suponiendo que la disipación viscosa es Mgn
cante. suponiendo que no se pierde calor a través tiva y reconociendo que las condiciones deben
del eje tar térmicamente desarrolladas por completo,
(c) Si la cubierta del cojinete se enfría con agua, de es la forma apropiada de la ecuación de ene
modo que la superficie externa del cojinete se man- Resuelva esta ecuación para la distribución deti
tiene a 30°C, determine las temperaturas del cojine- peraturas. ¿Cuál es el flujo de calor en la supe
superior (y — £.)?
te y del eje, 7 y 7
6.23 Considere el flujo de Couctte con transferencia de calor 6.25 Considere las ecuaciones de conservai ion (6.25.6.,
como se describe en el ejemplo 6.4 6 43)
a) Reacomode la distribución de temperaturas para ob- (a) Identifique cada ecuación y describa con bre
tener la forma adimensional el significado físico de cada término.
(b) Identifique las aproximaciones y cond cion
Qiv) = v U + í P r l - y)\ cialcs que se hacen para reducir estas expr
a las ecuaciones de capa limite (6.54 6.55 y
donde 6 = [7(y) — T{)]/[T¡ — 7 ) y r¡ = \ L Los
(c) Compare las ecuaciones 6.55 y 6.57. e iden
r ipos a limen tonales son el numero de Prandtl
Pt = p a p!k y el número de bckert Ec = U t p
las condiciones para las que las ecuaciones
la misma forma Comente la existencia de
(7, - 7( )
logia de transferencia de calor y de momento.
(b) Derive una expresión que establezca las condicio
nes en las que no habrá transferencia de calor a la
placa superior Similitud y parámetros adimensionales
(c) Derive una expresión para la transferencia de calor 6.26 En una buena aproximación, la viscosidad diná
a la p h ta inferior según las condiciones identifica- la conductividad térmica k y el calor especifico e
cas en la parte (b). independientes de la presión. ¿De qué forma la
■ Problemas 337
dad cinemática n y la difusividad térmica a varían con 6.30 Se encontró que lo> resultados experimentales para la
la presión para un líquido incompresible y para un gas transferencia de calor sobre una placa plana con una
ideal? Determine v y a del aire a 350 K para presiones superficie en extremo áspera están correlacionados por
de 1 y 1 0 atm. una expresión de la forma
DEPARTAMENTO DE BIBLIOTECA
Universidad Simón Bolívar S«de del Litoral
338 C apítulo 6 ■ Introducción a la convección
6.36 Considere el chip sobie la tarjeta de circuitos del pro- (c) Est me la razón de pérdida de masa por eva
blema anterior. Para asegurar una operación confiable ción cuando las temperaturas del agua y del i
en periodos extensos, la temperatura del chip no debe ambiente son 47°C, suponiendo que el coeficwj
exceder 85°C Suponiendo la disponibil dad de aire de transferencia de masa por convección perm
forzado a í « = 25°C y la aplieabilidad de la correla- ce sin cambio y el aire ambiental está seco
ción de transferencia de calor establecida, calcule y di-
6.42 La razón a la que se pierde agua debido a la cvaponl
buje la disipación de potencia máxima permisible del
cion de la superficie de un cuerpo de agua se detemí
chip Pc como función de la velocidad del aire para 1 ^
al medir la velocidad de retroceso de la superti
V ^ 25 m/s. Si la superficie del chip tiene una ennsivi-
Considere un día de verano para el que la tempe
dad de 0 80 y la tarjeta esta montada en un recinto
del agua y del aire ambiente es 305 K y la humedad
grande cuyas paredes están a 25°C, ¿cuál es el efecto
lativa del aire es 40%. ¿Si se sabe que la velocidad
de la radiación sobre la gráfica P, —V?
retroceso de la superficie es 0 I nun/h, cual es I; r
a la que se picide masa debido a la evaporación
A n alo g ía dt K cvu
tt o lil»
área de superficie unitaria? ¿Cuál es el coeficiente
6.37 Una placa plana delgada de 0.2 m por 0.2 m de lado se transferencia de masa por convección?
orienta de forma paralela a un flujo de aire atmosférico
con velocidad 40 m/s El aire esta a una temperatura 6.43 La fotosíntesis, como ocurre cn las hojas de una p
T x = 20°C, mientras la placa se mantiene a T s = verde, implica el transporte de dióxido de c
120°C. El aire fluye sobre las superficies superior e in- (C 0 2) de la atmósfera a los cloroplastos de las hi
ferior de la placa, y la medición de la fuerza de arrastre la rapidez de fotosíntesis se cuantifica en términos £
revela un valor de 0 057 N. ¿Cual es la transferencia de rapidez de asimilación de C ü 2 por los cloroplasto
calor desde ambos lados de la placa al aire? ta asimilación está fuertemente influenciada p<v
transferencia de C 0 2 a través de la capa límite que
6.38 El aire atmosférico esta en flujo paralelo (//<* = 15 m rs, produce cn la superficie de la hoja Bajo condi
Toe = 15 C) sobre una superficie plana de calentamien- para las que la densidad del C 0 2 es 6 X 10-4 kg m
to que se mantendrá a una temperatura de I40°C. El el aire y 5 X I0- 4 kg/nL en la superficie de la hoja
área de la superficie de calentamiento es 0.25 m2, y se coeficiente de transferencia de masa por cornees
sabe que el flujo de aire induce una fuerza de arrastre 1CT2 m/s, ( cuál es la rapidez de fotosíntesis en
de 0 25 N sobre el calentador. ¿Cuál es la potencia eléc- nos de k logramos de CO? asimilados por uní
trica necesaria para mantener la temperatura superficial tiempo y área de la superficie de la hoja?
establecida?
6.44 La especie A se evapora desde una superficie p
6.39 Para el flujo sobre una placa plana con una superficie la especie B Suponga que el perfil de concern
extremadamente áspera, se sabe que los efectos de la para la especie A en la capa limite de concentra
transferencia de calor por convección están correlacio- de la forma CA(y) = D y 2 + E y + F. donde D
nados mediante la expresión del problema 6.30. Para son constantes en cualquier posición v y y se
un flujo de aire a 50 m/s, ¿cuál es el esfuerzo cortante lo largo de una normal desde la superficie. D
de la superficie en v = 1 m desde el inicio de la placa? una expresión para el coeficiente de transiere
Suponga que el aire está a una temperatura de 300 K. masa por convección h„, en términos de estas c
tes. la concentración de A cn el flujo libre C\ ,
C o e fir ie n te s d e tr a n s fe r e n c ia d e m a sa fusividad de masa D AB. Escriba una expresión
6.40 En un día de verano la temperatura del aire es 27°C y flujo molar de transferencia de masa por c<r
la humedad relativa es 30%. El agua de la superficie de para la especie A
un lago se evapora a razón de 0 . 1 0 kg/h por metro cua-
drado de la superficie del área. La temperatura del agua Ecuación y solución «le conservación <lc <>|
también es 27°C. Determine el valor del coeficiente de 6.45 Cons dere el problema 6.24. cuando el fluido
transferencia de masa por convección mezcla binaria con diferentes concentraciones
6.41 Se observa que un contenedor de agua de 230 mm de diá- CA. y Ca. 2 en Ias superficies superior e inte
i
metro a 23°C tiene una razón de pérdida de masa de 1.5 X pectivamente. Para la región entre las placas,
10 kg/s cuando el aire ambiental esta seco y a 23°C la forma apropiada de la ecuación de conlinui1
especie A? Obtenga expresiones para la distrib'
(a) Determine el coeficiente de transferencia de masa
concentración de especies y el flujo de especia ^
por convección para esta situación.
superficie superior.
(b) Estime la razón de pérdida de masa por evapora-
ción cuando el aire ambiental tiene una humedad 6.46 Un esquema simple para desalineación implica
relativa de 50 por ciento 11er una película delgada de agua salada en las
■ ProblettMts 339
Condensado
Flujo de aire L
-7 0
l ; Película
delgada de
El flujo se mantiene por gravedad, y la especie gaseosa
agua salada
A fuera de la película es absorbida en la intcifaz líqui-
do-gas. La película está en un flujo laminar completa
mente desarrollado sobre toda la placa, de modo que
Hxiste un flujo de aire laminar incompresiblemente len-
sus componentes de la velocidad en las direcciones \ y
to entre las placas, de modo que el componente x de la
r son cero. La densidad de masa de A en y = 0 en el li
velocidad es finito mientras que los componentes y y r
quido es una constante pA independiente de .v
son cero Ocurre evaporación de la película líquida so-
bre la superficie inferior, que se mantiene a una tempe- (a) Escriba la forma apropiada de la ecuación de mo-
ratura elevada T0. mientras sucede la condensación en mento en v para la película. Resuelva esta ecuación
la superficie superior, que se mantiene a una temperatu- para la distribución del componente de la veloci
ra reducida Tl . Las concentraciones molares correspon- dad en x, u{y). en la película. Exprese el resultado
dientes de vapor de agua en las superficies inferior y en términos de 5, g. é , y las propiedades del liqui-
superior se designan como CA 0 y CA / , respectivamen- do ix y p. Escriba una expresión para la velocidad
te. La concentración y temperatura de especies se supo- máxima ur ‘ max*
5,47 Considere las ecuaciones de conservación (6.43) y (c) Si un coeficiente de transferencia de masa por con-
(6.52). vección se define como
DEPARTAMENTO d e b i b l i o t e c a
U n iv e rsid a d Simón Bouvar - Sede w
340 C apitulo 6 ■ Introducción a la convección
flujo de aire que contiene amoniaco (NH3) se mue- terior como caso 1 . considere los siguientes casos y de-
ve por el tubo, de modo que la densidad de masa de termine si las condiciones son análogas a las del caso
NH3 en la interfaz gas-líquido (no en el líquido) es 1 . Cada caso incluye un objeto de la misma forma, que
25 kg/m3. Se forma una solución diluida de amo- se suspende en un flujo de aire de la misma manera.
niaco en agua, y el cocfic ente de difusión es 2 X Donde exista un comportan! ento análogo, determine el
10 9 m2 /s. ¿Cuál es la rapidez de eliminación de valor correspondiente del coeficiente promedio de con-
masa de NH3 por absorción? vección.
(a) Los valores de 7 k, 7* y p permanecen igual. pero¿ =
S im ilit u d y a n a lo g ía d e t r a n s f e r e n c ia d e c a lo r 2 m y V' = 50 m/s.
\•f d e m a sa (b) Los valores de Ts y 7* permanecen igual, peroí.
2 m. V = 50 m/s. y p = 0.2 atm.
6.49 Considere el flujo cruzado del gas X sobre un objeto
que tiene una longitud característica L = 0.1 m Para (c) La superficie se cubre con una película líquida que
un número de Reynolds de 1 X 104, el coeficiente se evapora en el aire. Todo el sistema está a 300K,
promedio de transferencia de calor es 25 W/m • K. y el coeficiente de difusión para la mezcla aire \
t i mismo objeto después se impregna con el liquido por es Dab = 1.12 X 10- 4 m s. También. L - 2
y y se sujeta a las mismas condiciones de flujo D a- V = 50 m/s y p = 1 atm.
das las siguientes propiedades termofísicas, ¿cual es (d) La superficie se cubre con otra película liquida
el coeficiente promedio de transferencia de masa por ra la que 2>AB = 1.12 X IO- 4 m2/s. y el sistema
convección? tá a 300 K En este caso L = 2 ni. \ = 250 m.'i
p = 0 2 atm
v (mJ/s) k (W/m • K) a (m2/s)
6.53 Fn un día frío de abril se sabe que un corredor ve ■
Gas X 21 X I O 6 0.030 29 X I O 6
ligeramente pierde calor a razón de 500 W debido a
Líquido Y 3.75 X IO” 7 0.665 1.65 X IO- 7 convección al aire de los alrededores a 7 X = |()0C.
Vapor Y 4.25 X 10~ 5 0 023 4 55 X IO" 5 piel del corredor permanece seca y a una tempe
Mezcla de gas X vapor Y Se = 0 72 7 = 3ü°C Tres meses después, el corredor se mué
la misma velocidad, pero el día es caluroso y hu~
con una temperatura Tx = 30°C y una humedad re
6.50 Considere condiciones para las que un fluido con una
velocidad de corriente libre 1 = 1 m/s fluye sobre va de d>x = 60%. El corredor está ahora empapado
sudor y tiene una temperatura superficial uniforme
una superficie que se evapora o se sublima con una lon-
35 C. Bajo ambas condiciones las propiedades del
gitud característica L 1 m. lo que proporciona un
se suponen constantes con v = 1.6 X 10 5 m s.
coeficiente promedio de transferencia de masa por con-
vección lt„, — 10 2 m/s. Calcule los parámetros adimen- 0 026 W/m • K. Pr = 0.70. y D AB (vapor de agua-ai
2.3 X 10“ 5 m2/s.
sionales ShL, ReL, Se y j m para las siguientes
combinaciones: flujo de aire sobre agua, flujo de aire so- (a) ¿Cuál es la razón de pérdida de agua debida a
bre naftalina y glicerol caliente sobre hielo. Suponga una evaporación en el día de verano?
temperatura del fluido de 300 K y una presión de I atm. (b) ¿Cuál es la pérdida total de calor por conv
en el día de verano9
6.51 Dibuje la variación de los espesores de la capas límite
de velocidad y concentración con la distancia desde el
6.54 Ln objeto de forma irregular de 1 m de longitud
inicio de una placa plana para las siguientes condicio-
mantiene a una temperatura constante de IOO°C
nes de flujo laminar: flujo de aire sobre una película de
pende en un flujo de aire que tiene una temper
agua, flujo de aire sobre una capa de hielo seco, flujo
flujo libre de 0 C, presión atmosférica de I atm
de aire sobre una placa de naftalina y flujo de glicerol
velocidad de 120 m/s. La temperatura del aire
caliente sobre una capa de hielo, que se funde y disuel-
en un punto cercano al objeto en el flujo de airees
ve en el glicerol Para cada caso suponga una tempera-
Un segundo objeto de la misma forma tiene 2
tura media del fluido de 300 K.
longitud y se suspende en un flujo de aire de I:
6.52 Un objeto de lonna irregular tiene una longitud carac- forma. La velocidad de flujo libre del aire es 60 «
terística L = 1 m y se mantiene a una temperatura su- aire y el objeto están a 50°C. y la presión total«
perficial uniforme 7 V= 325 K. Este se suspende en un Un recubrimiento plástico sobre la superficie del
lujo de aire que está a presión atmosférica (p = 1 atm) se seca mediante este proceso El peso molecular
y tiene una velocidad V = 100 m/s y una temperatura por es 82 y la presión de saturación a 50°C para
Too = 275 K. El flujo promedio de calor de la superficie terial plástico es 0.0323 atm La difusividad
a\ aire es 12.000 W/m2 Refiriéndose a la situación an- para el vapor en aire a 50°C es 2.60 X JO' 5
Problemas 34.1
(a) Para el segundo objeto, cn una posición que corres- (b) Determine el coeficiente promedio de transferencia
ponde al punto de medición del primer objeto, deter- de calor por convección h para el perfil de tamaño
mine la concentración de vapor y la presión parcial. completo. Ln = 60 mm, cuando se expone a un
(b) Si el flujo de calor promedio q" es 2(KX) W/m2 para flujo de aiic libre con V = 60 m/s, = 184°C, y
el primer objeto, determine el flujo de masa pro- p 00 = 1 atm cuando Ts = 70°C.
medio n"A (kg/s • m2) para el segundo objeto. (c) El área de la superficie del perfil se expresa como
6.55 Un proceso industrial implica la evaporación de agua A s = 2.2 L¡¡ • /, donde / es la longitud normal a la
de una película líquida que se forma sobre el perímetro de página. Para las condiciones de la parte (b). ¿cuál
una superficie. Se hace pasar aire seco sobre la superfi- es el cambio en la transferencia de calor al perfil si
cie, y de mediciones de laboratorio la correlación de la longitud característica LH se duplica?
transferencia de calor por convección es de la forma 6.57 Glicerol caliente fluye sobre una capa de hielo cuya
forma es tal que los efectos de transferencia de calor
ÑuL = 0 A 3 R elx Pr°-A
por convección se correlacionan mediante una ecua-
(a) Para una temperatura del aire y velocidad de 27°C y ción de la forma N uL = 0.25 Re¡01 P r m . El hielo se
1 0 m/s, respectivamente, ¿cuál es la rapidez de eva- funde y disuelve en el glicerol. Bajo condiciones para
poración de una superficie de I m2 de área y longi- las que el coeficiente promedio de transferencia de ca-
tud característica L = I m? Aproxime la densidad lor por convección es /// = 100 W/m2 • K, ¿cuál es el
de vapor saturado como pA sa, = 0.0077 kg/m3. coeficiente promedio de transferencia de njti a por con-
(b) ¿Cuál es la temperatura de estado estable de la pe- vección?
lícula líquida? 6.58 Considere las condiciones del problema 6.12, peí o con
656 Un perfil aerodinámico que sostiene la envoltura de un una película delgada de agua sobre la superficie. Si el
cojinete se expone a un flujo de aire caliente del escape aire es seco y el número de Schmidt Se es 0.6, ¿cuál es
de un motor. Es necesario llevar a cabo experimentos el flujo de masa evaporativo? ¿Hay transferencia neta
pura determinar el coeficiente promedio de transferen- de energía hacia o desde el agua?
cia de calor por convección h del aire al perfil a fin de 6.59 Considere las condiciones del problema 6.4. para el
poder enfriarlo a la temperatura superficial que se de- que un experimento de transferencia de calor dio la
sea Ts. Se decide realizar experimentos de transferencia distribución que se establece del coeficiente local de
de masa sobre un objeto de la misma forma y para ob- convección. h x(x). El experimento se realizó para tem-
tener los resultados de transferencia de calor que se de- peraturas de superficie y flujo libre de 310 y 290 K,
sean mediante el uso de la analogía de transferencia de respectivamente. Ahora considere repetir el experimen-
calor y de masa. to bajo condiciones cn que la superficie se cubre con
una capa delgada de naftalina y tanto la superficie co-
= 60 mm
V mo el aire están a 300 K. ¿Cual es el valor correspon-
diente del coeficiente promedio de transfeiencia de
r. masa por convección, hm ¡1
P ee
6.60 Con el uso de la técnica de sublimación de naftalina, la
distribución radial del coeficiente local de transferencia
os experimentos de transferencia de masa se llevaron de masa por convección para un flujo uniforme normal
acabo con un modelo del perfil de la mitad del tamaño a un disco circular está correlacionada por una expre-
construido de naftalina expuesta a un flujo de aire a sión de la forma
27°C. Las mediciones de transferencia de masa dan es-
tos resultados: h J r )D T / r \n
ShL Re¡
282 60,0(H)
El número de Sbtrvvood (Shv) del punto de estanca-
491 120 .0 0 0 miento depende de los números de Reynolds (R eD =
568 144.000 VD/v) y de Schmidt (SC = tVDAn), y los datos se corre-
989 288.000 lacionan mediante la siguiente expresión:
d epa r t a m en t o de b ib l io t ec a
Universidad «Jiiujit uOn«ar - Sede «
342 C apitulo 6 ■ Introducción a ¡a convección
ción significativa en el contenido de agua del plumaje y (a) Determine la velocidad a la que el agua seev»
los polluelos sucumbirán por deshidratación. de la toalla, nA (kg/s).
(b) Lleve a cabo un balance de energía sobre la ft
V para determinar la transferencia neta de radir
<7itk1(W). a la toalla. Determine la irradiación G{\V
1 OQ
T
de- 6.63 Ln primavera, las superficies como aceras y ca_
asfálticas están algunas veces húmedas en la ma
aun cuando no llueva durante la noche. Lasco
nes nocturnas típicas se muestran en el dibujo.
para comprender mejor la transferencia convectiva du-
rante el vuelo, se llevaron a cabo estudios en un túnel
de viento con el uso de modelos del urogallo. Mediante
el calentamiento de la parte del modelo que correspon-
de al plumaje donde se recoge el agua, se determinó el
/ I _____
coeficiente de transferencia de calor por convección. Tx = 290 K, 0 x = 0.7 t- = 240 K
Los resultados para diferentes velocidades de aire y ta- h = 53 W/m2 ■K
maños de modelo se usaron para desarrollar una corre
laeión empírica de la forma B risa ^ T< : ™
ÑiiL = 0.0 M R e f 5 P r m
f' f '€'€ : f f ü l C f f r i ' f Capa
El arca efectiva de la superficie de la parte que recoge de agua
agua del plumaje se designa como Av y la longitud ca- Concreto |¡[Link]
racterística se define como L = (A J12.
Problemas 343
(a) Determine los flujos dc calor asociados con la con- 6.66 Se sabe que en noches claras la temperatura del aire no
vección. í/”onv’ evaporación. <7eVap. e intercambio dc necesita caer por debajo dc 0 C para que una capa del-
radiación con los al rededo re s.í/'^d gada de agua sobre la tierra se congele. Considere tal
(b) ¿Sus cálculos explican por qué el concreto está capa de agua en una noche clara para la que la tem pe-
mojado en lugar de seco? Explique con brevedad ratura efectiva del cielo es —30 C y el coeficiente dc
transferencia de calor por convección debido al m ovi-
te) Fluye calor de la capa liquida al c o n c r e t o <O del
concreto a la capa líquida? Determine el (lujo dc miento del viento es h = 25 W /m 2 • K. Se supone que
calor por conducción hacia adentro o hacia afuera el agua tiene una emisividad de 1.0 y que está aislada
del concreto. de la tierra en lo que a la conducción respecta
(a) Sin tom ar en cuenta la evaporación, determine la
6.64 Aire seco a 32°C fluye sobre una placa mojada de lon- tem peratura más baja que tendrá el aire sin que se
gitud 200 mm y ancho I ni (caso A). Un calentador congele el agua.
eléctrico empotrado sum inistra 432 W y la temperatura
(b) Para las condiciones dadas, estime el coeficiente
de la superficie es 27 C.
de transferencia de masa para evaporación del agua
(m/s).
l r = 32°C 7 , = 32°C
(c) Tome en cuenta ahora el efecto de la evaporación,
Película de agua ¿cual es la tem peratura más baja que tendrá el aire
Aire
r 27°C
L sin congelamiento del agua? Suponga que el aire esta
seco
u ,
6.67 Una expresión para la presión parcial real dc vapor de
Caso A Caso B agua en términos de las temperaturas de bulbo húmedo
y bulbo seco, denominada ecuación de Carrier. está da-
(a) ¿Cual es la velocidad de evaporación del agua des- da como
de la placa (kg/h)?
(b) Después de un largo periodo de operación, toda el
agua se evapora de la placa y la superficie está se- Pv = P g\v
1810 —7:wb
ca (caso B). Para las mismas condiciones de flujo
libre y la misma potencia de calentam iento que en donde p „ p ^ , y p son la presión parcial real, la presión
el caso A. estime la tem peratura de la placa, Ts. de saturación a la temperatura de bulbo húmedo y la
presión total (todos en bar), mientras que Ths y Tbh son
las temperaturas dc bulbo seco y húmedo en kclvin.
Enfriamiento e v a p o r a t i v o
Considere aire a 1 atm y 37.8°C que fluye sobre un ter-
6.65 Un exitoso ingeniero de California instaló una piscina m ómetro de bulbo húmedo que indica 21.1°C.
caliente circular en su patio y encontró que, para las (a) Con la ecuación de Carrier. calcule la presión par-
condiciones típicas de operación que se muestran aba- cial del vapor de agua en el flujo libre. ¿Cuál es la
jo. debe agregar agua a ra/ó n dc 0.001 kg/s a fin dc humedad relativa?
mantener un nivel fijo de agua en la piscina.
(b) Refiérase a una carta psicrométrica y obtenga la
T* = 290 K humedad relativa de forma directa para las condi-
Aire ciones que se indican, y compare el resultado con
<6X = 0.30
la parte (a).
(c) Use el valor dc la presión de vapor que se obtiene
en la parte (a) con la relación de enfriamiento eva-
porativo, ecuación 6.98, para calcular una diferen-
cia de temperaturas entre el aire global y el bulbo
húmedo Compare esta diferencia de temperatura
con la diferencia dc temperaturas real. ¿Cuál es el
porcentaje de diferencia entre los dos valores?
6.68 Un termómetro dc bulbo húmedo consiste en un termó-
Si a piscina esta bien aislada cn sus lados y en la parte metro de mercurio en vidrio cubierto con una tela hú-
inferior y la temperatura del agua dc la toma es igual a meda (agua). Cuando se suspende en un flujo de aire,
tedel agua de la piscina, ¿a qué rapidez deben sum inis- la lectura de estado estable del term ómetro indica la
trar energía los calentadores eléctricos para m antener el tem peratura de bulbo húmedo obtenga una expre-
agua a 310 K sión para determ inar la humedad relativa del aire a par-
UCPARTAMENTO DE BIBLIOTECA
Universidad Simón bolívar - Sede u . Litoral
341 C apítulo 6 ■ Introducción a la convección
tir del conocim iento de la tem peratura del aire (7®), la 6.71 Se conduce un experim ento para determinar el c
temperatura de bulbo húmedo, y las propiedades de cíente de transferencia de masa por convección de
aire y vapor de agua apropiadas. Si 7® = 45°C y 7 bh = pequeña gota, con un calentador controlado para o
25°C, ¿cuál es la humedad relativa del flujo de aire? a una tem peratura constante. La historia de la potcnc*
requerida para evaporar por completo la gota a una
6.f>9 Un proceso industrial implica la evaporación de una pe-
tem peratura de 37°C, se m uestra en el dibujo. Se ob
lícula delgada de agua de una superficie pennietrica al
va que, a medida que la gota se seca, el diámetro
calentarla desde abajo y al forzar aire a través de ésta.
medo sobre la superficie del calentador permanece c%i
Mediciones de laboratorio para esta superficie propor-
constante cn un valor de 4 mm.
cionan la siguiente correlación de transferencia de calor.