Problemas transitorio

PROBLEMAS SOBRE
REGIMEN TRANSITORIO
pfernandezdiez.es
pfernandezdiez.es Régimen transitorio.III.-37

III.1.- Se tiene una placa de 10 cm de espesor, a una temperatura uniforme de 20ºC, que se introduce en un medio
a 100ºC, adquiriendo instantáneamente esta temperatura.
Determinar mediante técnicas numéricas, el tiempo necesario para que el plano medio de la pieza alcance una
temperatura de 60ºC.
α = 6.10-6
(m2
/seg)
_________________________________________________________________________________________
RESOLUCIÓN
a) Tiempo necesario para que el plano medio de la pieza alcance una temperatura de T = 60ºC.
Ecuación diferencial unidimensional
∂T
∂t
= α
∂2T
∂x2
= α T1 - 2 T + T2
Δx2
= T1 = T2 = α 2 T1 - 2 T
Δx2
= 2 α T1 - T
Δx2
€
Integrándola:
dT
T1 - T
=
2 α
Δx2
dt ;
293
333
∫
dT
T1 - T
=
2 α
Δx2 0
t
∫ dt
- ln (T1 - T)293
333
=
2 α t
Δx2
; - ln
373 - 333
373 - 293
=
2 x 6.10-6
t
0,052
⇒ t = 144,4 seg
*****************************************************************************************
III.2.- Un calentador plano y grande, con una temperatura superficial de 100°C, se coloca en contacto directo con
la tierra k = 2,0 W/m°K,α = 5.10-7 m2/seg. Si el suelo estaba originariamente a una temperatura uniforme de
15°C, determinar:
a) La temperatura del plano situado a 15 cm por debajo del calentador, dos horas después de que se aplique
calor por contacto
b) La cantidad de calor total conducido hacia el interior de la tierra por unidad de área en las dos primeras
horas
_________________________________________________________________________________________
RESOLUCIÓN
a) Temperatura del plano situado a 15 cm por debajo del ca-
lentador, dos horas después de que se aplique calor
Dado que se puede considerar a la Tierra como un sólido se-
miinfinito, cuando el calentador entra en contacto con ella la
condición de contorno se puede considerar isotérmica,por cuan-
to hay un cambio brusco de la temperatura de la superficie; la
distribución de temperaturas es:
T( x,t) - Ts
T0 - Ts
= G ( u) = G (
x
2 α t
) = G (
0,15
2 5.10 −7
x 2 x 3600
) = G (1,25) = 0,92287
T(x,t) = Ts + G(u) (T0 - Ts) = 100 + 0,92287 x (15 - 100) = 21,55°C
b) Cantidad de calor total conducido hacia el interior de la tierra por unidad de área en las dos primeras horas:
Q(t) = 2 k (Ts - T0 ) t
α π
= 2 x 2 W
mºC
x (100 - 15)ºC
7200 seg
5.10-7
π (m2
/seg )
= 6,394 kW
m2
*****************************************************************************************
III.3.-Una placa gruesa de hormigón (α = 0,049.10-5 m2/seg) que esta inicialmente a 400ºK se rocía con una gran
cantidad de agua a 300ºK. Determinar el tiempo necesario para que un punto situado a 5 cm
de la superficie se enfríe hasta 320ºK.
_____________________________________________________________________________
RESOLUCIÓN
Hay que suponer que la velocidad del rociado permite mantener la superficie a 300ºK; la placa se
puede tratar como un sólido semi ∞ con condición de contorno isotérmica, dados los posibles va-
lores de hcF y k.
T( x,t) - Ts
T0 - Ts
= G ( u) = 320 - 300
400 - 300
= 0,2 ⇒ u = 0,18 = x
2 α t
=
0,05 m
2 0,049.10-5
(m2
/seg) t
⇒ t = 39370 seg

La profundidad xt de penetración de la temperatura es la posición en la que la tangente al perfil de temperatura en x =
0 corta a la recta T0 = 400ºK.
x t = π α t = π x 0,049.10-5
x 39370 = 0,246 m
*****************************************************************************************
III.4.- Una pared de hormigón α = 7.10-7 m2/seg está inicialmente a una temperatura uniforme de T0 = 25°C. Re-
pentinamente una de sus caras aumenta su temperatura y pasa a 125°C, manteniéndose esta temperatura indefi-
nidamente. Calcular la temperatura a 5, 10 y 15 cm de la superficie caliente, 30 minutos después del ascenso
brusco de la temperatura superficial de la pared.
_________________________________________________________________________________________
RESOLUCIÓN
Al haber un cambio repentino en la temperatura de la superficie, la condición de contorno es isotérmica
Para 30 minutos = 1800 segundos, el valor de “u” es:
u = x
2 α t
= x
2 7.10-7
(m2
/seg) x 1800 seg
= 14,08 x
€
a) Para: x = 0,05 m ; u =14,08 x 0,05 = 0,7042 ⇒ G(u) = 0,6806
T(x,t) = Ts + G(u)(T0 - Ts) = 125 + 0,68 x (25 - 125) = 56,9°C
€
b) Para: x = 0,1 m ; u =14,08 x 0,1 =1,408 ⇒ G(u) = 0,95358
T(x,t) = Ts + G(u)(T0 - Ts) = 125 + 0,95358 x (25 - 125) = 29,64°C
€
c) Para: x = 0,15 m ; u =14,08 x 0,15 = 2,1126 ⇒ G(u) = 0,9971
T(x,t) = Ts + G(u)(T0 - Ts) = 125 + 0,9971 x (25 - 125) = 25,29°C
Este último resultado está muy próximo a la T0 inicial, por lo que el calor no ha penetrado
en los 30 minutos a esa profundidad.
*****************************************************************************************
III.5.- Una placa de acero de 1 cm de espesor se saca de un horno a 600ºC y se sumerge en un baño de aceite a
30ºC. Sí se calcula que el coeficiente de transferencia de calor es de 400 W/m2ºK, ¿cuánto tiempo tarda la placa
en enfriarse hasta 100ºC?
Tomar para el acero: k= 50 W/mºK, ρ = 7800 kg/m3 y cp = 450 J/kgºK.
_________________________________________________________________________________________
RESOLUCIÓN
Condición de contorno a utilizar:
hC /k = 400/50 = 8 ⇒ c.c. convección
Bi = hC (e/2)/k = 400 x 0,005/50 = 0,04 < 0,1 ⇒ c.c. R.T.I. despreciable
⎧
⎨
⎩
por lo que admite, al menos, dos métodos de resolución; elegiremos la condición de resistencia térmica interna des-
preciable, en la que para placa plana el nº de Bi es el mismo.
Fo = α t
(V/A)2
= α t
(e/2)2
= α = k
ρ cp
=
50 (W/mK)
7800 (kg/m3
) x 450 (J/kgK)
= 1,4245.10-5 m2
seg
=
1,4245.10-5
( m2
seg
) t
(0,01/2)2
= 0,5698 t
Distribución de temperaturas:
T - TF
T0 - TF
= e- Bi Fo
⇒
100 - 30
600 - 30
= 0,1228 = e-0,04 x 0,5698 t
⇒ t = 92 seg
*****************************************************************************************
III.6.- Un experimento sobre procesado de materiales en condiciones de microgravedad en un transbordador es-
pacial exige el enfriamiento de una pieza en un flujo forzado de algún gas inerte. Una esfera metálica de 1 cm de
diámetro sale de un horno a 800ºC y se enfria hasta 500ºC mediante un flujo de gas nitrógeno a 25ºC. El coefi-
ciente de transferencia de calor por convección para el enfriamiento es de 25 W/m2ºK.
Las propiedades del material de la esfera son: ρ = 14000 kg/m3 ; cp = 140 J/kgºK ; emisividad,ε = 0,1.
Las superficies del entorno pueden considerarse como casi negras a 25ºC.
Determinar el tiempo necesario del experimento
_________________________________________________________________________________________
RESOLUCIÓN
Al tratarse de un enfriamiento en un gas, el modelo de enfriamiento con resistencia térmica interna despreciable
(r.T.I.D.) se puede suponer válido, aunque a posteriori habrá que comprobar.

El coeficiente de transferencia de calor lo será por convección y por radiación h = hc + hr
Distribución de temperaturas:
T(t) - TF
T0 - TF
= e-(Bi)(Fo)
h r = ε σ (T0
2
+ TF
2
) (T0 + TF ) = 0,1 x 5,67.10-8 W
m2
ºK 4
(10732
+ 2982
) (1073 + 298)ºK3
= 9,64
W
m2
ºK
Bi Fo =
h (V/A)
k
α t
(V/A)2
=
h c + h r
k
(k/ρ cp ) t
(4 /3) π R3
4 π R 2
=
3 (h c + h r )
ρ cp R
=
3 (25 + 9,64) t
14000 x 140 x 0,005
= 0,010604 t
500 - 25
800 - 25
= 0,613 = e- Bi Fo
= e- 0 ,010604 t
⇒ t = 46,15 seg
Comprobación: Bi =
h (R/3)
k
=
34,6 x (0,005/3)
k
=
0,05767
k
0,1 (ya que para un metal k 10 W/mºC)
Fo =
9 k t
ρ cp R2
=
9 k x 46,15 seg
14000 (kg/m3
) x 140 (J/kgºC) x 0,0052
m2
= 8,47 k
*****************************************************************************************
III.7.- Una conducción de 10 cm de diámetro exterior está enterrada en el suelo, de forma que éste se encuentra a
una temperatura uniforme de 10ºC. En el invierno se produce una caída de temperatura a nivel del suelo del or-
den de -15ºC y se mantiene así durante 50 días. Datos: α = 0,2 .10-6 m2/seg ; k = 1,5 W/mºK
Determinar
a) La profundidad mínima a colocar la tubería para que el agua que circula por su interior no congele
b) Calor que se disipa al exterior cada día
_________________________________________________________________________________________
RESOLUCIÓN
a) Profundidad mínima a la que habrá que colocar la tubería para que el agua que circula por su interior no con-
gele.- Se puede interpretar la Tierra como un sólido semiinfinito, de forma que el enfriamiento provoque una situa-
ción de contorno isotérmica, ya que su efecto a mucha profundidad sería inapreciable
€
u =
x
2 α t
=
α = 0,2.10-6 m2
seg
t = 50 x 24 x 3600 seg = 4,32.106
seg
=
x
2 0,2.10-6
x 4,32.106
= 0,5379 x
G(u) =
T(x,t) - Ts
T0 - Ts
=
0 - (-15)
10 - (-15)
= 15
25
= 0,6 ⇒
G(u) = 0,60386 ; u = 0,60
G(u) = 0,58792 ; u = 0,58
⇒ u = 0,59515
luego la profundidad mínima para que el agua no congele es: u = 0,59515 = 0,5379 x ; x = 1,1064 m
b) Calor que se disipa al exterior cada día
q( t) =
- k (T0 - Ts )
π α t
=
- 1,5 {10 - (-15)} (W/m)
π x 0,2.10 -6
t m
=
- 47308
t
W
m2
€
Q(t) =
0
t
∫ q(t) dt =
0
86400
∫
- 47308
t
dt = 27,81.106 W.seg
m2
= 7,725
kW
m2
Superficie A de evacuación por metro de longitud de tubería:A = 2π R L = π x 0,1 x 1 = 0,31416 m2
Q(t) = 7,725 kW
m2
x 0,31416 m = 2,427 kW
metro de longitud de tubería
*****************************************************************************************
III.8.- Una pared de hormigón a prueba de incendios, de 15 cm de espesor, está reforzada por una cubierta de
acero de 4 mm de espesor. El conjunto queda repentinamente expuesto a la acción de una fuente de calor radian-
te, que se puede suponer como un cuerpo negro, a 1000ºK.
Determinar el tiempo necesario para que la superficie alcance 500ºK si la temperatura
inicial de la pared es de 300ºK.
khormigón = 1,4W/mºK ; αhormigon = 0,75.10-6
m2
/seg
k acero= 59 W/mºK ; absortancia de la superficie muy oxidada del acero = 0,9

_________________________________________________________________________________________
RESOLUCIÓN
Se pueden hacer las siguientes suposiciones:
a) La pared de hormigón se puede considerar como un sólido semiinfinito que inicialmente está a una temperatura T0
b) El valor q0 es el calor que absorbe la plancha de acero y que va a transmitir al hormigón en el intervalo de tempe-
raturas considerado.
q 0 = αabsortividad σ (Trad
4
- T0
4
) = 0,9 x 5,67.10-8
x (1000 4
- 300 4
) = 50616 (W/m 2
)
c) La variación de temperaturas ΔT a través de la cubierta de acero de espesor L = 0,004 m, en régimen cuasi-estacio-
nario es:
ΔT =
q 0
k/L
=
47840,6
59/0,004
= 3,24ºK (que es despreciable frente a 1000ºK)
que dice que la caida de temperatura a través de la cubierta de acero es insignificante; el flujo de calor q0 se transmi-
te al hormigón a través de su superficie (x = 0) no influyendo en ello la cubierta de acero.
La distribución de temperaturas en el hormigón, para esta situación, es:
T(x,t) = T0 +
2 q0
k
α t [ e-u 2
π
- u {1 - G(u)}] , con: u = x
2 α t
que evaluada para x = 0 resulta:
T(0,t) = Ts = T0 +
2 q 0
k
α t
π
⇒ t = π
α
(
k (Ts - T0 )
2 q 0
)2
= π
0,75.10 -6
(
1,4 x (500 - 300)
2 x 50616
)2
= 32 seg
Para comprobar la suposición de sólido semi ∞, se calcula la profundidad de penetra-
ción xt, y se tiene:
xt =
Ts - T0
- (
∂T
∂x
)x=0
= k
Ts - T0
qs
= 1,4
500 - 300
50616
= 5,53.10-3 m = 0,553 cm
que es menor que los 15 cm de espesor de la pared
*****************************************************************************************
III.9.- Un habitáculo refrigerado a 5ºC ha estado funcionando durante varios años y se le pone fuera de servicio,
de forma que a partir de este momento se permite que el aire del ambiente a 27ºC circule libremente por él. Deter-
minar el tiempo que tardará la tierra situada a 1 m de profundidad en alcanzar una temperatura de 15ºC. Se su-
pondrá que el coeficiente de transferencia de calor por convección es de 3 W/m2ºK, siendo las propiedades térmi-
cas para la tierra mojada (k = 2,6 W/mºK, α = 0,45.10-6 m2/seg).
_________________________________________________________________________________________
RESOLUCIÓN
A la tierra se la puede suponer, en este caso, como un sólido semi∞ que se expone al aire a 27ºC, con una temperatu-
ra inicial de 5ºC a una profundidad apreciable. Como hC/k = 3/2,6 = 1,15 se puede suponer que la condición de con-
torno es por convección. La distribución de temperaturas es:
T - T0
TF - T0
= 1 - G (u ) - {1 - G (u + η )} e- Bix + η
=
15 - 5
27 - 5
= 0,4545 =
€
=
η = Bi2
Fo =
hC
2
α t
k2
=
32
x 0,45.10-6
t
2,62
= 6.10-7
t
u =
x
2 α t
=
1 m
2 0,45.10-6
t
=
745,35
t
Bix =
hC x
k
=
3 (W /mºK) x 1 m
2,6
= 1,154
= 1 - G (
745,35
t
) - {1 - G (
745,35
t
+ 6.10-7
t} e- 1,154 + 6.10-7
t
= 0,4545
Su resolución es iterativa. Para tener una referencia inicial hay que partir de un valor de t razonable, que puede ser el
correspondiente a la condición isotérmica, en la forma:
0,4545 = 1 - G( u) ⇒ G( u) = 1 - 0,4545 = 0,5455 ⇒ 0 = 0,53 =
745,35
t
⇒ t = 1,97.106
seg
El tiempo real será superior a éste.
En la iteracción se van dando valores a t obteniéndose los correspondientes a
T - T0
TF - T0
:
4 5 6 6,05 6,06
0,368 0,415 0,452 0,454 0,4545
10 −6
t
(T - T0 )/(TF - T0 )
luego: t = 6,06.106 segundos = 70,13 días

*****************************************************************************************
III.10.- Para medir las fluctuaciones de la temperatura en la pared del cilindro de un motor diesel de 5 mm de es-
pesor, se instala un termopar a 1 mm de la superficie interior de la pared. En un determinado ensayo el motor
funciona a 1000 rpm, y se observa que la media de la lectura del termopar es de 322ºC con una amplitud de
0,79ºC.
Si se supone que la variación de la temperatura es aproximadamente senoidal, calcular la amplitud y diferencia
de fase de la variación de temperatura en la superficie interior.
Datos de la pared: α = 12.10-6 m2/seg ; k = 40 W/mºK
_________________________________________________________________________________________
RESOLUCIÓN
La impresión inicial es que al tener la pared del motor un espesor pequeño y dado además que la superficie exterior
está refrigerada, existirá un gradiente de temperatura a través de la pared, lo que podría eliminar la condición de sóli-
do semi∞. Sin embargo, si la onda térmica se amortigua hasta anularse a una distancia muy corta de la superficie, sí
se puede considerar a la pared del cilindro como sólido semi∞ sometido a una variación periódica de su temperatura
superficial (x = 0).
La variación periódica de la temperatura superficial es:
Tsup - Tm
Tmáx - Tm
= sen w t, siendo la temperatura media: Tm =
Tmáx + Tmín
2
Distribución de temperaturas en la pared, (es la respuesta de la temperatura resultante en cualquier punto):
T( x,t) - Tm
Tmáx - Tm
= e
- x w
2 α sen (w t - x w
2 α
), en la que :
Función de amortiguación: e
- x w
2 α
Función periódica: sen (x w
2 α
- w t)
⎧
⎨
⎪
⎩
⎪
La amplitud de la variación de la temperatura disminuye exponencialmente a medida que penetra en la pared, de-
sarrollándose con un desfase igual a: x w
2 α
, por lo que la amplitud de la variación de temperatura (Tmáx- Tm) en
la superficie interior (x = 0) es:
T( x,t) - Tm
Tmáx - Tm
=
T* - Tm
Tmáx - Tm
= T* - Tm = 0,79ºC =
0,79
Tmáx - Tm
= e
-0 ,001 w
2α sen ( w t - 0,001 w
2 α
) =
=
Período: T* = 2 π
w
= 2 π
πn/30
= 60
n
= 60
1000
= 0,06 seg ⇒ w = 104,7 rad
seg
w
2 α
=
104,7
2 x 12.10 -6
= 2089 m-1
= e
-0 ,001 w
2 α = e- 0,001x 2089
= 0,1238
La amplitud de la variación de temperatura en la superficie interior es:
0,79
Tmáx - Tm
= e
-0 ,001 w
2α = 0,1238 ⇒ Tmáx - Tm = 6,38ºC
La diferencia de fase en la superficie (x=0) interior es: x
w
2 α
= 0,001 x 2089 = 2,09 rad = 119,75º
*****************************************************************************************
III.11.- La temperatura de un chorro de gas se mide con un termopar; la unión del mismo se puede considerar
aproximadamente como una esfera de diámetro d = 1 mm, k = 25 W/m°C, ρ = 8400 kg/m3, cp = 0,4 kJ/kg°C. El
coeficiente de película entre el chorro de gas y la unión del termopar es hC = 560 W/m2°C.
Determinar:
a) El tiempo que tardará el termopar en registrar el 99% de la temperatura del chorro
b) El tiempo que tardará el termopar en registrar el 99,99% de la temperatura del chorro. ¿Podría medir la

temperatura exacta?¿Por qué?
_________________________________________________________________________________________
RESOLUCIÓN
Dado que el fluido exterior es un gas, y el sólido un metal, es muy probable que la condición de contorno sea de re-
sistencia térmica interna despreciable; en efecto:
hC /k = 560/25 = 22,4 ⇒ c.c. convección
Bi = hC R/k = 560 x 0,5.10-3
/25 = 0,0112 0,1 ⇒ c.c. R.T.I. despreciable
⎧
⎨
⎩
luego para la condición de contorno con resistencia térmica interna despreciable :
Bi =
hC (V/A)
k
=
h C
(4/3) π R3
4 π R 2
k
=
h C
R
3
k
=
560 x
1.10-3
6
25
= 3,73.10-3
T - TF
T0 - TF
= e- Bi Fo
= e
-
h C A t
ρ V cp
= e
-
hC t
ρ
R
3
cp
= e
-
560 ( Wm2 ºC /)t
8400 ( kg/m3 ) x
0,001
6
m x 400 ( J/kg ºC )
= e- t
, con t en segundos
a) Tiempo que tardará el termopar en registrar el 99% de la temperatura del chorro
Para que el termopar registre el 99% de la temperatura del gas, el valor de (T - TF) tiene que ser igual a la fracción de
temperatura que queda sin precisar: % (T0 - TF), es decir:
T - TF
T0 - TF
= error que se comete
T - TF = 1% (T0 - TF ) =
1
100
( T0 - TF ) ⇒
T - TF
T0 - TF
=
1
100
= e-t
⇒ t = 4,6 seg
b) Tiempo que tardará el termopar en registrar el 99,99% de la temperatura del chorro. ¿Podría medir la tempera-
tura exacta?¿Por qué?
Para un 99,99% de precisión:
T - TF
T0 - TF
= 1
10000
= e-t ; t = 9,21 segundos
Para un 100% de precisión, la temperatura medida sería la exacta:
T - TF
T0 - TF
= 1
∞
= e-t ; t = ∞ segundos
*****************************************************************************************
III.12.- Una placa de Al, k = 160 W/m°C, ρ = 2790 kg/m3, cp = 0,88 kJ/kg°C, de espesor (e = 3 cm) y con una tem-
peratura uniforme T0 = 225°C se sumerge, agitándola, en un fluido que se encuentra a una temperatura constan-
te TF = 25°C. El coeficiente de transmisión de calor entre la placa y el fluido es hC = 320 W/m2°C.
Determinar
a) El tiempo requerido para que el centro de la lámina alcance 50°C
b) El calor intercambiado (en Kcal), si la superficie de la lámina es de 1 m2.
_________________________________________________________________________________________
RESOLUCIÓN
a) Tiempo requerido para que el centro de la lámina
alcance 50°C.- La convección se realiza por las dos
caras de superficie A, luego:
Bi =
hC (e/2)
k
=
320 x (0,03 /2)
160
= 0,03 0,1
por lo que el sólido se comporta con resistencia tér-
mica interna despreciable.
T - TF
T0 - TF
= 50 - 25
225 - 25
= 0,125 = e- Bi Fo
=
α =
k
ρ cp
=
160 ( W/mºC)
2790 ( kg/m3
) x 0,88 (kJ/kgºC)
= 0,06517.10-3 m2
seg
Fo =
α t
(L/2 )2
=
0,06517.10-3
(m2
/seg) t
0,0152
m2
= 0,2896 t
=
= e- 0 ,03x 0 ,2896 t = 0,125 ⇒ t = 240 seg
que es el tiempo requerido para que tanto el centro de la placa, como cualquier otro punto de ella, incluida la superfi-
cie, alcance 50°C.
e

b) Calor intercambiado, si la superficie de la lámina es de 1 m2.
Calor intercambiado en Kcal si: A = 2 x 1 m2 ; t = 240 segundos; Fo = 0,2896 x 240 = 69,50
Q(t) = hC A (T0 - TF) t 1 - e-Bi.Fo
Bi.Fo
= 320 (W/m2°C)x 2 (m2) (225 - 25)°Cx 240 (seg)
1 - exp(-3 x 10-2 x 69,50)
3 x 10-2 x 69,50
=
= 121902.295 (W.seg) = 121902.295 (Joules) =
12.902,295 (kJoul)
4,186 kJoul
Kcal
= 3.082 Kcal
*****************************************************************************************
III.13.- Una plancha eléctrica tiene una base de acero k = 70 W/m°C, ρ = 7840 kg/m3, cp = 0,45 kJ/kg°C y pesa 1
kg. Dicha base tiene una superficie de calentamiento A = 0,025 m2 y se calienta con el calor procedente de una re-
sistencia eléctrica que proporciona Q = 250 W. Inicialmente la temperatura de la placa base es T0 = 20°C; a medi-
da que la placa se va calentando disipa calor por convección a un fluido exterior que se encuentra a TF = 20°C,
con un coeficiente hcF = 50 W/m2°C.
Calcular la temperatura de la placa al cabo de 5 minutos después de iniciado el calentamiento de la misma.
_________________________________________________________________________________________
RESOLUCIÓN
Flujo de calor: q =
Q
A
= 250 W
0,025 m2
= 10.000 W
m2
Cálculo del espesor e: Masa = V ρ = (Espesor x A) ρ = e A ρ
e =
Masa
A ρ
=
1 kg
0,025 m2
x 7840 (kg/m3
)
= 0,51.10-2
m
Φ =
X
m
+ (Φ 0 -
X
m
) e- mt
=
X
m
=
q
h C
= 10000
50
= 200
m =
h C
ρ e cp
= 50
7840 x 0,51.10 -2 x 450
= 0,002778
=
= 200 + (0 - 200) e-0,002778 x 5 x 60
= 113,09 = T - TF = T - 20 ⇒ T = 133,1ºC
Temperatura de equilibrio: Φ∞ = X
m
= 200 ; T(t) - TF = T( t) - 20 = 200 ; T( t) = 220 ºC
*****************************************************************************************
III.14.- Una placa grande tiene una difusividad térmica α = 5.10-5 m2/seg , k = 20 W/mºC y un espesor de 0,30 m.
La placa se encuentra inicialmente a la temperatura de 0°C. Su superficie se pone en contacto térmico con un
cuerpo sólido a 500°C, de forma que casi instantáneamente, la superficie del cuerpo frío adquiere la temperatura
de 500ºC. Determinar:
a) La distribución de temperaturas
b) Su historia, a lo largo de los 100 primeros segundos de contacto con los cuerpos calientes.
c) Tiempo que debe transcurrir para que la temperatura en el plano central sea de 300ºC
_________________________________________________________________________________________
RESOLUCIÓN
a) Distribución de temperaturas
Condición de contorno isotérmica (contacto de 2 cuerpos sólidos en donde la superficie del cuerpo frío se pone casi
instantáneamente a la temperatura del otro, 500ºC). Longitud característica = espesor
Fo =
α t
L2 =
5.10-5 t
0,32 = 5,55.10-4 t
T - Ts
T0 - Ts
=
4
π ∑ e- λn
2 αt sen (λn x)
n
= λn =
π n
L
; ξ =
x
L
=
4
π ∑ e
- π n
L2
α t sen (π n ξ)
n
T - 500
0 - 500
=
4
π
n=1,3,5
∞
∑ e
-
(π n )2 x 5.10-5 t (seg)
0,32
sen
π n x
0,3
n
=
4
π
n=1,3,5
∞
∑ e-0,005483 t n2 sen (10,47n x)
n

T = 500 -
2000
π n= 1,3,5
∞
∑ e- 0 ,005483 t n 2 sen (10,47 n x)
n
b) Su historia, a lo largo de los 100 primeros segundos de contacto de los cuerpos para diversas profundidades.
Para: t = 100 seg ⇒ Fo =
α t
L2
=
5.10-5 x 100
0,32
= 0,0555 0,2 (varios sumandos)
ξ= 0,1 ⇒ Φ
Φ0
= 0,24 ; ξ= 0,3 ⇒ Φ
Φ0
= 0,6 ; ξ= 0,4 ⇒ Φ
Φ0
= 0,69 ; ξ= 0,5 ⇒ Φ
Φ0
= 0,72
c) Tiempo que debe transcurrir para que la temperatura en el plano central sea de 300ºC:
TC = 500 -
2000
π
e-0 ,005483 t
sen (10,47 x 0,15) + ... = 300 ⇒ t = 212,5 seg
*****************************************************************************************
III.15.- Una esfera de aluminio de 6 cm de diámetro se calienta en un horno hasta 200°C y mediante un dispositi-
vo adecuado se introduce rápidamente en un conducto por donde circula aire a 50°C. Haciendo medidas de tem-
peraturas en la esfera se obtienen los siguientes resultados:
- Un minuto después de colocar la esfera en el conducto, T1=182°C
- Cinco minutos después de colocar la esfera en el conducto, T2=129,5°C
- Diez minutos después de colocar la esfera en el conducto, T3=92,5°C
Los datos del aluminio son constantes en el intervalo 50ºC÷200°C: ρ = 2690 kg/m3 ; cp = 0,248 Kcal/kg°C ;
k=235 W/m°C
Determinar:
a) La distribución de temperaturas
b) El coeficiente de transmisión de calor medio, esfera-aire en Kcal/h.m2.°C
c) El flujo térmico instantáneo q/A en la superficie de la esfera en los tiempos 1 minuto, 5 minutos y 10 minu-
tos y la relación Q/Q0 en los tiempos citados.
_________________________________________________________________________________________
RESOLUCIÓN
Dado que el medio de enfriamiento es el aire, se puede suponer que la condición de contorno es con resistencia tér-
mica interna despreciable, que se tendrá que confirmar con Bi 0,1
Bi =
h cF (R/3)
k
=
hcF x 0,01
235
= 0,00004255 h cF
a) Distribución de temperaturas:
Φ
Φ 0
= e- Bi Fo
= e
-
h cF A
ρ V c p
= e
-
3 h cFt
ρ R c p
T - 50
200 - 50
= e- Bi Fo
⇒ T = 50 + 150 e
-
3 h cFt
ρ R c p
= 50 + 150 e
-
3 h cFt
2690 x 0 ,03 x 0,248
= 50 + 150 e- 0,15 h cFt
que predice la historia de la relación entre el tiempo y la temperatura siempre que Bi0,1
b) El coeficiente de transmisión de calor medio, esfera-aire en Kcal/hm2°C
Para: t = 60 seg ; T1 = 182 = 50 + 150 e- 0,15 h cF1 60
⇒ h cF1
= 0,0142
Kcal
seg.m2 ºC
= 51,13
Kcal
h.m2 ºC
Para: t = 5 min ; T2 = 129,5 = 50 + 150 e- 0 ,15 h cF2 300
⇒ h cF2
= 0,0141
Kcal
seg.m2 ºC
= 50,79
Kcal
h.m2 ºC
Para: t = 10 min ; T3 = 92,5 = 50 + 150 e- 0,15 h cF3 600
⇒ hcF3
= 0,01401
Kcal
seg.m2 ºC
= 50,436
Kcal
h.m2 ºC
Valor medio de hcF =
51,13 + 50,79 + 50,436
3
= 50,78 Kcal
h.m2.ºC
= 59,04 W
m2.ºC
Bi = 0,00004255 hcF = 0,00004255 x 59,04 W
m2.ºC
= 0,00251 0,1
luego el sólido sí tiene resistencia térmica interna despreciable.
c) Flujo térmico instantáneo en la superficie de la esfera:
Q
A
= hcF (T(t) - TF) = hcF (T0 - TF) exp (- Bi Fo) =

= hcF ( T0 - TF ) exp (
h cF t
ρ cp (R/3)
) = 50,78 Kcal (200 - 50) exp (
50,78 t
2690 x 0,248 (0,03/3)
) = 7618,5 exp (- 7,61 t)
Para: t = 1 minuto =
1
60
hora ⇒ q( t)/A = 7618,5 e
-7,61 x
1
60 = 6710
Kcal
h m2
Para: t = 5 minutos =
1
12
hora ⇒ q(t)/A = 7618,5 e
-7,61 x
1
12 = 4039,5
Kcal
h m2
1
6
hora ⇒ q(t)/A = 7618,5 e
-7,61 x
1
6 = 2142
Kcal
h m2
d) Hallar la relación Q/Q0 en los tiempos citados.
Calor disipado en estos tiempos: Q(t) = hcF A (T0 - TF) t
1 - exp (- Bi Fo)
Bi Fo
=
= 50,78 Kcal
h.m2.ºC
x (4 π x 0,0 32) x 150 t x
1 - exp (- 7,61 t)
7,61 t
= 11,32 [ 1 - e- 7,61 t]
Q0 = m cp (T0 - TF ) = V ρ cp (T0 - TF ) =
4 π R 3
3
ρ cp (T0 - TF ) =
=
4 π 0,033 m3
3
x 2690
kg
m3
x 0,248
Kcal
kgºC
(200 - 50 )ºC = 11,32
Kcal
hora
€
Para: t = 1 minuto =
1
60
hora ⇒ Q(t) = 11,32 (1 - e
-7,61 x
1
60 ) = 1,34
Kcal
hora
⇒
Q
Q0
= 0,1183
€
1
12
hora ⇒ Q(t) = 11,32 (1 - e
-7,61 x
1
12 ) = 5,31
Kcal
hora
⇒
Q
Q0
= 0,4696
€
1
6
hora ⇒ Q(t) = 11,32 (1 - e
-7,61 x
1
6 ) = 8,13
Kcal
hora
⇒
Q
Q0
= 0,7186
*****************************************************************************************
III.16.- Un cilindro de fundición de gran longitud, k= 70 W/m°K, α= 2.10-5 m2/seg, con un diámetro de 20 cm,
está originariamente a una temperatura uniforme de 400°C. Se enfría la superficie exterior del cilindro con un
fluido a la temperatura de 50ºC y coeficiente de convección hC = 420 W/m2°K. Si el fluido se dirige hacia la super-
ficie del cilindro durante un tiempo de 20 minutos, determinar:
a) La temperatura en el eje y en la superficie del cilindro en ese instante
b) El calor transferido por unidad de longitud desde el cilindro, durante el periodo de los 20 minutos
_________________________________________________________________________________________
RESOLUCIÓN
Bi = hC R/k = 420 x 0,1/70 = 0,6 0,1 ⇒ c.c. convección
⎧
⎨
⎩
Fo = α t
R 2
=
2.10 -5
(20 x 60)
0,12
= 2,4 ; Φ 0 = T0 - TF = 400 - 50 = 350ºC
a) Temperatura en el eje del cilindro en ese instante
Φ C
Φ0
=
TC - TF
T0 - TF
=
TC - 50
400 - 50
=
Fo = 2,40
Bi = 0,6 ; Bi-1 = 1,666
= 0,095 ⇒ TC = 83,25ºC
Temperatura de la superficie del cilindro en ese instante
Φ(R,t)
Φ(0,t)
=
Tsup - TF
TC - TF
=
Tsup - 50
83,25 - 50
=
r
R
= 1 = 0,74 ⇒ Tsup = 74,6ºC
ó también:
Φ(R,t)
Φ 0
=
Φsup
ΦC
ΦC
Φ 0
= 0,74 x 0,095 = 0,0703 ⇒ Tsup = 50 + (0,0703 x 350) = 74,6ºC
b) Cantidad de calor transferido por unidad de longitud desde el cilindro, durante 20 minutos
€
Bi2
Fo = 0,62
x 2,4 = 0,864

Q(t)/Q0 = 0,88 (Gráfica de Gröber)
Q0 = ρ cpF π R 2
(T0 - TF ) =
k
α
π R 2
( T0 - TF ) =
70
2.10 −5
π x 0,12
x 350 = 3,848.10 7 W.seg
m
Q(t) = 0,88 Q0 = 0,88 x 3,848.107 W seg
m
= 3,388.107 W seg
m
=
3,388.10 7
3600
W hora
m
= 9,41 kW
m
*****************************************************************************************
III.17.- Una lámina de un cierto material de propiedades térmicas, k = 2 Kcal/hmºC, α = 0,0002 m2/hora, tiene un
espesor de 3 cm, y se encuentra a una temperatura de 24ºC. En estas condiciones se introduce en un horno a
400ºC. El coeficiente de convección es hcF = 60 Kcal/hm2ºC.
Determinar
a) El tiempo que tardará el centro geométrico en alcanzar la temperatura de 300ºC y la temperatura que se al-
canza en ese instante en un plano situado a 1 cm del plano central
b) Si el calentamiento continúa, el tiempo necesario para que en el plano situado a 1 cm del plano central se
alcancen 350ºC
c) Si existiese un aislamiento perfecto en una de las caras de la lámina, el tiempo que deberá transcurrir para
que en el plano central se alcancen 200ºC.
_________________________________________________________________________________________
RESOLUCIÓN
a) Tiempo que tardará el centro en alcanzar la temperatura de 300ºC. En la condición de contorno de convección,
la longitud característica es la mitad del espesor
Bi = hC (e/2)/k = 60 x 0,015/2 = 0,45 0,1 ⇒ c.c. convección
⎧
⎨
⎩
Mediante la gráfica de Heysler se obtiene:
Bi = 0,45 ; Bi-1 = 2,22
Φ C
Φ0
=
TC - TF
T0 - TF
=
300 - 400
24 - 400
= 0,27
⎫
⎬
⎪
⎭
⎪
⇒ Fo = 3,56 =
α t
L2
=
2.10-4 m2
hora
t
0,0152 m2
⇒ t = 4 horas
Temperatura que se alcanza en ese instante en un plano situado a 1 cm del plano central
Φ
Φ C
=
T1cm - TF
T(0,t) - TF
=
T1cm - 400
300 - 400
= Bi-1
= 2,22 ;
x
L
=
1
1,5
= 0,66 = 0,92 ⇒ T1cm = 400 - 0,92 x 100 = 308ºC
b) Si el calentamiento continúa, el tiempo necesario para que en el plano situado a 1 cm del plano central se al-
cancen 350ºC
Temperatura en el plano central en estas condiciones (1ª gráfica de Heysler)
Φ
Φ C
=
T1cm - TF
T(0,t) - TF
=
350 - 400
TC - 400
= Bi = 0,45 ;
x
L
=
1
1,5
= 0,66 = 0,92 ⇒ TC = 345,6ºC
por lo que entrando en la 2ª gráfica de Heysler se obtiene:
Bi-1
= 2,22
Φ C
Φ0
=
T(0,t) - TF
T0 - TF
=
345,6 - 400
24 - 400
= 0,145
⎫
⎬
⎭
⇒ Fo = 5,1 =
α t
L2
=
2.10-4 m2
hora
t
0,0152
m2
; t = 5,7 horas
c) Si existe un aislamiento perfecto en una de las caras de la lámina, el tiempo que deberá transcurrir para que en
el plano central se alcancen 200ºC se calcula considerando el problema en forma similar a los casos anteriores pero
con un espesor doble = 6 cm, por cuanto el plano aislado se considera como plano adiabático.
La longitud característica es: L* = 3 cm.
La temperatura en la cara aislada, que ahora es el plano central de la nueva configuración, se calcula a partir de la
primera gráfica de Heysler:
Φ
Φ C
=
T(x,t) - TF
T(0,t) - TF
=
200 - 400
T(0,t) - 400
=
Bi =
hC L*
k
=
60 x 0,03
2
= 0,9
Bi-1
= 1,11 ; Plano central: x
L*
= 0,5
= 0,915 ⇒ T(0,t) = 181,4ºC

Con esta temperatura y la segunda gráfica de Heysler se obtiene el nº de Fo y el tiempo correspondiente:
Bi-1
= 1,11
Φ C
Φ0
=
T(0,t) - TF
T0 - TF
=
181,4 - 400
24 - 400
= 0,58
⎫
⎬
⎭
⇒ Fo = 3,83 = α t
L2
=
2.10-4
(m2
/hora) t
0,032
m2
⇒ t = 3,735 horas
*****************************************************************************************
III.18.- Una placa grande de acero aleado k= 30 W/m°K y α= 1,5.10-5 m2/seg sale de un tren de laminación a una
temperatura constante de 800°C. La placa tiene un espesor de 30 cm. Para conseguir enfriar la placa se dirigen
sendos chorros de aire a gran velocidad sobre ambas caras de la placa. El aire tiene una temperatura de 30°C y el
coeficiente de transferencia de calor entre el aire y la superficie de la placa es de 500 W/m°K. A la superficie de la
placa ha de aplicarse una capa de material aislante, pero la temperatura de la misma tiene que estar por debajo de
los 200°C, antes de que pueda aplicarse la capa.
Determinar el tiempo mínimo necesario durante el cual deberá estar enfriándose la placa antes de que pueda apli-
carse la capa aislante.
_________________________________________________________________________________________
RESOLUCIÓN
⎧
⎨
⎩
(En este ejercicio al aire se le impone un coeficiente de convección muy elevado)
Φ
Φ C
=
T(x,t) - TF
TC - TF
=
200 - 30
TC - 30
=
Bi =
hC L
k
=
500 x 0,15
30
= 2,5
Bi-1 = 0,4 ; x/L = 1
= 0,415 ⇒ TC = T(0,t) = 439,6ºC
Tiempo mínimo necesario durante el cual deberá estar enfriándose la placa antes de que pueda aplicarse la capa
aislante.
Bi-1
= 0,4
Φ C
Φ0
=
T(0,t) - TF
T0 - TF
=
439,6 - 30
800 - 30
= 0,532
⎫
⎬
⎭
⇒ Fo = 0,59 =
α t
L2
=
1,5.10 -5
(m2
/seg)
0,152
m2
; t = 885 seg
*****************************************************************************************
III.19.- Las paredes de un recuperador de calor están formadas por ladrillos de sílice de 40 cm de espesor. Una de
las caras se encuentra aislada térmicamente, mientras que la otra se calienta por gases que están a 900°C con un
coeficiente hcF = 100 Kcal/hm2°C. Si la T0 es de 700°C, se desea conocer:
a) La temperatura en un plano situado a 20 cm de la superficie caliente al cabo de 1 hora y 10 horas
b) Los flujos térmicos para el apartado (a)
c) La energía almacenada en la pared
_________________________________________________________________________________________
RESOLUCIÓN
Datos de la sílice:
cp = 0,24 Kcal
Kg°C
; k = 1 Kcal
h m °C
; ρ = 2000
Kg
m3
; α = k
ρ cp
= 1
2000 x 0,24
= 0,0020 m2
hora
a) Temperatura en un plano situado a 20 cm de la superficie caliente al cabo de 1 hora.- Al tratarse de un material
refractario, con un bajo coeficiente k, la condición de contorno puede ser isotérmica ya que:
hC /k = 100/1 = 100 ⇒ c.c. isotérmica
Fo =
α t
L2
=
0,002 (m2 /hora) x 1 hora
0,4 2
m2
= 0,0125 0,5 ⇒ c.c. isotérmica
⎧
⎨
⎪
⎩
⎪
Fo = 0,0125 confirma que se puede tratar el problema como sólido semiinfinito, siempre que lejos de la superficie se
mantenga la temperatura aproximadamente inicial a la inicial.
A) Tratamiento como sólido semiinfinito con condición de contorno isotérmica
Para t = 1 hora, en un plano situado a la distancia x = 0,2 m de la superficie, se tiene:

u =
x
2 α t
=
0,20 m
2 0,002
m2
hora x 1 hora
= 2,236 ⇒ G(2,236) = 0,9982
T0,2 m; 1h = TF + G(u) (T0 - TF ) = 900 + 0,9982 (700 - 900) = 700,36ºC
Para t = 10 horas, en un plano situado a la distancia x = 0,2 m de la superficie, se tiene:
€
Fo =
α t
L2 =
0,002 (m2
/hora) x 10 h
0,42
m2 = 0,125 ; u =
x
2 α t
=
0,20
2 0,002 x10
= 0,707 ⇒ G(0,707) = 0,6825
T0,2 m; 10h = TF + G(u) (T0 - TF ) = 900 + 0,6825 (700 - 900) = 763,5ºC
Para ver si éste resultado es válido, hay que comprobar que lejos de la superficie, (en este caso el espesor correspon-
diente a la cara aislada x = 0,4 m) la temperatura prácticamente no se ha modificado:
u =
x
2 α t
=
0,40
2 0,002 x 10
= 1,414 ⇒ G(1,414) = 0,95445
T0,4 m; 10h = TF + G(u) ( T0 - TF ) = 900 + 0,95445 (700 - 900) = 709,1ºC
valor próximo a 700°C, por lo que puede considerarse como una buena aproximación.
También se puede comprobar utilizando la expresión de la profundidad de penetración: xt = π α t
¿Hasta cuando es sólido semi-∞?; cuando: xt = 0,4 = π α t ⇒ t ≅ 25 horas
b) Los flujos térmicos para el apartado (a) , sólido semiinfinito:
q( t) =
- k ( T0 - TF )
π α t
⇒
q1 hora =
- 1 (Kcal/h mºC) (700 - 900)ºC
π x 0,002 (m2
/hora) x 1 hora
= 2523,1 Kcal
h m2
q10 horas =
- 1 ( Kcal/h mºC) (700 - 900)ºC
π x 0,002 (m2
/hora) x 10 hora
= 797,88
Kcal
h m2
⎧
⎨
⎪
⎩
⎪
c) Energía almacenada en la pared
Q1hora = 2
π
( TF - T0 ) ρ cp k t = 2
π
(900 - 700)ºC 2000
kg
m3
x 0,24 Kcal
kgºC
x 1 Kcal
h mºC
x 1 hora = 4944 Kcal
m2
B) Tratamiento como sólido infinito con condición de contorno isotérmica.- En estas condiciones, el plano aislado
se considera como plano de simetría del sólido, y éste se trata como placa plana de espesor L = 0,80 m.
T( x,t) - TF
T0 - TF
= 4
π n=1,3,5..
∞
∑ e- λ n
2
α t sen (λ n x)
n
= λ n = π n
L
= 4
π n=1,3,5..
∞
∑ e- (π n )2
Fo sen (π n ξ)
n
Para: x = 0,20 m ; ξ = 0,2/0,8 = 1/4
Para t = 1 hora ⇒ Fo = α t
L2 =
0,002 (m2 /hora) 1 h
0,82 m2 = 0,003125 0,2 (Se necesitan varios sumandos)
n = 1 ⇒ exp{-(1 π)2
x 0,003125}
sen (1 π x 0,25)
1
= 0,6856
n = 3 ⇒ exp{-(3 π)2
x 0,003125}
sen (3 π x 0,25)
3
= 0,17857
n = 5 ⇒ exp{-(5 π)2
x 0,003125}
sen (5 π x 0,25)
5
= - 0,0654
n = 7 ⇒ exp{-(7 π)2
x 0,003125}
sen (7 π x 0,25)
7
= - 0,02228
n = 9 ⇒ exp{-(9 π)2
x 0,003125}
sen (9 π x 0,25)
9
= 0,0083
La suma del desarrollo es: 0,7847, y la temperatura a la distancia x = 0,2 m de la superficie, para t = 1 hora es:
T0,2 m; 1 h = TF +
4 (T0 - TF )
π x 0,7847 = 900 +
4 (700 - 900)
π x 0,7847 = 700,16ºC
que coincide con la del apartado (a).
Para: t = 10 horas ⇒ Fo = 0,03125
n = 1 ; e- (π)2 x 0,03125 sen (π x 0,25)
1
= 0,5194
n = 3 ; e- (3π)2 x 0,03125 sen (3π x 0,25)
3
= 0,01468
n = 5 ; e- (5π)2 x 0,03125 sen (5π x 0,25)
5
= - 0,00006338
⎧
⎨
⎪
⎩
⎪
x

La suma de estos valores es: 0,5340166, por lo que: T(x,t) = 900 +
0,5340166 x 4
π
(700 - 900) = 764,01°C
Flujo de calor en la mitad del sólido:
€
q(t) =
4 Φ0 k
L n=1,3,...
∑ e-(π n)2
Fo
=
n = 1 ; e-(π)2
0,03125
= 0,7346
n = 3 ; e-(3π)2
0,03125
= 0,06229
n = 5 ; e-(5π)2
0,03125
= 0,000448
=
4 x (700 - 900) x 1
0,8
0,7973 = 797,3
Kcal
h m2
*****************************************************************************************
III.20.- Una placa de acero está inicialmente a una temperatura uniforme de 300°C. La placa se enfría enviando
sobre ella una corriente de aire a 50°C, de modo que incida sobre una de las superficies de la placa. Esta tiene un
espesor de 10 cm, difusividad térmica α = 10-5 m2/seg, conductividad térmica k = 40 W/m°K, siendo el coeficiente
de película aire- placa hC = 400 W/m2°K.
Determinar
a) El tiempo necesario para que la superficie de la placa sobre la que incide el aire se enfríe hasta 200°C
b) Cuando la superficie de la placa alcance los 200°C, las temperaturas en los planos situados a 1 cm. y a 10
cm. de la misma.
_________________________________________________________________________________________
RESOLUCIÓN.- Si Fo1 se puede tratar como sólido semiinfinito siempre que a una gran distancia de la superficie
la temperatura sea igual a la temperatura inicial (dato que se pide en el apartado (b) para la temperatura en el plano
situado a 10 cm de la superficie)
Bi =
hC e
k
=
400 x 0,1
40
= 1 ; Fo =
α t
e2
=
10-5
t
0,12
= 10 -3
t
por lo que durante 10-3 t 1, (1000 seg), el sólido se comporta como sólido semiinfinito y la temperatura T0 se man-
tiene prácticamente constante lejos de la superficie; por encima de 1000 seg se comportaría como infinito.
a) Tiempo necesario para que la superficie de la placa sobre la que incide el aire se enfríe hasta 200°C.- Del dia-
grama correspondiente se obtiene:
Para: x = 0 ⇒
x
2 αt
= 0 ;
T( 0,t) - T0
TF - T0
=
200 - 300
50 - 300
= 0,40 Gráfica
⎯ →⎯⎯
Bi Fo = 0,5
Fo = 0,25 = 10-3
t ⇒ t = 250 seg
⎧
⎨
⎩
Distribución de temperaturas en un sólido semiinfinito sometido a convección
b) Cuando la superficie de la placa alcance los 200°C, las temperaturas en los planos situados a 1 cm. y a 10 cm.
de la misma son:
€
Temperatura para x = 0,1 m:
x
2 α t
=
0,1
2 10-5
x 250
= 1
hC
k
α t = 0,5
⎧
⎨
⎪
⎩
⎪
⎫
⎬
⎪
⎭
⎪
⇒
Φ
Φ0
=
T(x,t) - T0
TF - T0
= 0,035
Tx=0,1 = 300 + 0,035 (50 - 300) = 291,25°C 300°C
aproximación suficiente para considerar al sólido como semi∞
€
Temperatura para x = 0,01 m:
x
2 α t
=
0,01
2 10-5
x 250
= 0,1
hC
k
α t = 0,5
⎧
⎨
⎪
⎩
⎪
⎫
⎬
⎪
⎭
⎪
⇒
Φ
Φ0
=
T(x,t) - T0
TF - T0
= 0,325
Tx=0,01 = 300 + 0,325 (50 - 300) = 218,75°C

*****************************************************************************************
III.21.-- Unas bolas para cojinetes de acero al Cr k=50 W/m°K, α=1,3.10-5 m2/seg han de tratarse térmicamente.
Para ello se calientan hasta una temperatura de 650°C y luego se enfrían en un baño de aceite que tiene una tem-
peratura de 55°C. Las bolas de los cojinetes tienen un diámetro de 4,0 cm. El coeficiente de transferencia de calor
por convección entre las bolas y el aceite del baño es de 300 W/m2°C.
Determinar
a) El tiempo durante el cual deben permanecer las bolas en el aceite antes de que su temperatura descienda
hasta los 200°C
b) Calor eliminado en cada bola en este intervalo de tiempo
_________________________________________________________________________________________
RESOLUCIÓN
Bi = hC R/k = 300 x 0,02/50 = 0,12 0,1 ⇒ c.c. convección
⎧
⎨
⎩
mediante el método gráfico se puede suponer una temperatura periférica límite de 200ºC, ya que el núcleo estará a
mayor temperatura
a) Tiempo durante el cual deben permanecer las bolas en el aceite antes de que su temperatura descienda hasta
los 200°C .- Con la primera gráfica de Heysler, para Bi-1 = 8,33 y ξ = 1 se obtiene que Φ/ΦC = 0,95, por lo que:
Φ periferia
Φ C
=
Tperiferia - TF
TC - TF
=
200 - 55
TC - 55
= 0,95 ⇒ TC = 207,6ºC
Con la 2ª gráfica de Heysler se obtiene el nº de Fo y el tiempo.
Φ C
Φ0
=
TC - TF
T0 - TF
=
207,6 - 55
650 - 55
= 0,26 ⇒ Fo = 4,6 = α t
R2
=
1,3.10 -5
(m2
/seg) t (seg)
0,022
m2
⇒ t = 141,5ºC
Energía almacenada inicialmente en la bola:
Q0 = 4
3
π R3 ρ cpF (T0 - TF) = ρ = k
α cpF
= 4
3
π (0,02) 3 50
1,3.1 0-5 cpF
cpF (650 - 55) = 76687 J
Transferencia de calor para: t = 141,5 seg
Bi2
Fo = 0,122
x 4,6 = 0,06624
Bi = 0,12
⎫
⎬
⎭
Gröber
⎯ →⎯⎯⎯
Q
Q0
= 0,87 ⇒ Q = 0,87 Q0 = 0,87 x 76687= 66717 J
DE OTRA FORMA.- Como estamos en la frontera de Bi 0,1 también se podría considerar la condición de resisten-
cia térmica interna despreciable, cometiendo un error del orden del 5%. Si este error se admite, se tiene:

Bi =
h cF (R/3)
k
=
300 x (0,02/3)
50
= 0,04
Fo = α t
(R/3)2
=
9 x 1,3.10-5 t
0,022
= 0,2925 t
a) Tiempo durante el cual deben permanecer las bolas en el aceite antes de que su temperatura descienda hasta
los 200°C
Φ(t)
Φ 0
= e- Bi Fo
=
T(x,t) - TF
T0 - TF
⇒ e- 0,04 x 0,2925 t
=
200 - 55
650 - 55
= 0,2437 ⇒ t = 120,6 seg
Q(t) = hcF A (T0 - TF ) t
1 - e- Bi Fo
Bi Fo
= Fo = 0,2925 t = 0,2925 x 120,6 = 35,28 =
= 300 x 4 π x 0,022 (650 - 55) x 120,5 x
1 - e- (0,4 x 35,28)
0,04 x 35,28 = 57930 J
Transferencia de calor instantánea, procedentes de las bolas cuando se acaban de sumergir en el baño de aceite y
cuando alcanzan los 200°C.
En t = 0 ⇒ Fo = 0 ; q(t) = h cF A (T0 - TF ) e-Bi Fo
= 300 x 4 π x 0,02 2
(650 - 55) = 897,2 W
En t = 120,6 seg ⇒ Fo = 35,28 ; q(t) = 300 x 4 π x 0,02 2
(650 - 55) e- 80,04 x 35,28
= 218,8 W
*****************************************************************************************
III.22.- Una plancha metálica de espesor 3 cm se encuentra a una temperatura de 20ºC y en estas condiciones se
introduce en un horno a 1000ºC.
Si la plancha se considera de grandes proporciones, determinar:
a) El tiempo que debe transcurrir para que el centro alcance 500ºC
b) La temperatura que en ese instante adquiere el plano situado a 1 cm del plano central
c) El calor absorbido en ese intervalo de tiempo
d) El tiempo que debe transcurrir para que el plano situado a 1 cm del plano central alcance la temperatura
de 750ºC
e) Si se considera que una de las caras de la placa tiene un aislamiento térmico perfecto, ¿qué tiempo deberá
transcurrir para que en el plano central se alcancen 500ºC? ¿cuál sería en ese instante la temperatura en la cara
no aislada? Datos: k = 8 W/m.ºC ; α =5.10-3 m2/hora; hC = 93 W/m2 ºC
_________________________________________________________________________________________
RESOLUCIÓN
⎧
⎨
⎩
a) Tiempo que debe transcurrir para que el centro alcance 500ºC
Temperatura adimensional en la zona central:
Φ C
Φ0
=
T(0,t) - TF
T0 - TF
=
500 - 1000
20 - 1000
= 0,51
Bi-1
= 5,734
⎫
⎬
⎪
⎭⎪
⇒ Fo = 4,8 = α t
L2
= 5.10 3
t
0,0152
⇒ t = 0,216 h = 13 min
b) Temperatura que en ese instante adquiere el plano situado a 1 cm del plano central
Fo = 4,8 ; x
L
= 1 cm
1,5 cm
= 0,66
Bi-1
= 5,734
⎫
⎬
⎭
⇒
Φ
Φ C
=
T1 cm - 1000
500 - 1000
= 0,96 ; T1 cm = 520 ºC
c) Calor absorbido en ese intervalo de tiempo
Q0 = V ρ cpF (T0 - TF ) =
α =
k
ρ cp
⇒ ρ cp =
k
α
=
8 (W/mºC)
5.10-3 (m2/h)
V = a e = 0,03 A (m3 )
=
= 0,03 A(m3 )
8 (W/mºC)
5.10-3 (m2/h)
(20 - 1000)ºC = 47,04 A (
kW
m2 )

Bi2Fo = 0,17442 x 4,8 = 0,146 ⇒
Q
Q0
= 0,53
Q0 = 47,04 A ( kW
m2
)
⎫
⎬
⎪
⎭
⎪
⇒ Q = 0,53 Q0 = 0,53 x 47,04
kW
m2 = 24,93
kW
m2
d) Tiempo que debe transcurrir para que el plano situado a 1 cm del plano central alcance la temperatura de
750ºC
Φ
Φ C
=
T - TF
TC - TF
=
750 - 1000
TC - 1000
= 0,96 ⇒ TC = 1000 +
750 - 1000
0,96
= 739,6ºC
luego:
T(0,T) - 1000
20 - 1000
=
739,6 - 1000
20 - 1000
= 0,2657
Bi-1
= 5,734
⎧
⎨
⎩
⇒ Fo = 8,8 = α t
L2
=
5.10 -3 m2
hora
t
0,0152
⇒ t = 0,396 h = 23,45 min
e) Si se considera que una de las caras de la placa tiene un aislamiento térmico perfecto, el tiempo que deberá
transcurrir para que en el plano central se alcancen 500ºC
En caso de aislamiento perfecto en una de las caras, el problema es similar al caso ante-
rior pero con L*= e*/2 = 3 cm
Bi =
h cF L*
k
=
93 ( W/m2
ºC) 0,03 m
8 (W/mºC)
= 0,34875 ; Bi-1
= 2,867
Temperatura adimensional en la zona central del sólido duplicado:
T - TF
T(0,t) - TF
= 500 - 1000
T(0,t) - 1000
= 0,96 ; T(0,t) = 479,17ºC (superficie aislada)
Con esta temperatura y los gráficos de Heysler se obtiene el nº de Fo:
luego:
T - 1000
20 - 1000
=
479,17 - 1000
20 - 1000
= 0,5314
Bi-1
= 2,867
⎧
⎨
⎩
⇒ Fo = 2,1 =
α t
L* 2
=
5.10-3
t
0,032
⇒ t = 0,378 horas
Temperatura en la cara no aislada:
x/L = 1
Bi-1
= 2,867
⎧
⎨
⎩
⇒
T - TF
TC − TF
= 0,84
T(x,t) = TF + 0,84 {T(0,t) - TF ) = 1000 + 0,84 (479,17 - 1000) = 562,5ºC
*****************************************************************************************
III.23.- Un cilindro de 10 cm de diámetro y una longitud de 16 cm, k=0,5 W/m°K, α=5.10-7 m2/seg está inicial-
mente a la temperatura uniforme de 20°C. Se coloca el cilindro en un horno en el que la temperatura ambiental
es de 500°C con un hC=30 W/m2°K.
Determinar
a) Las temperaturas máxima y mínima del cilindro, treinta minutos después de su introducción en el horno
b) El calor absorbido en ese tiempo
_________________________________________________________________________________________
RESOLUCIÓN
h cF
k
=
30
0,5
= 60 ⇒ C.c. convección
Distribución de temperaturas, intersección de una placa plana ∞ de espesor 16 cm con un cilindro ∞ de diámetro 10
cm: P(x) C(r) =
Φ(x, r, t)
Φ0
PROCESO TÉRMICO:.- En el calentamiento, en el tiempo t, la temperatura mínima se localiza en el centro geomé-
trico del cilindro y la temperatura máxima se corresponde con los puntos de las circunferencias de sus bases.
Temperatura mínima en: x = 0, r = 0
Temperatura máxima en: x = L, r = R
Placa infinita ⇒ Bi =
h cF L
k =
30 x 0,08
0,5 = 4,8 ; Bi-1 = 0,208 ; Fo =
α t
L2 =
5.10-7 x 1800
0,082 = 0,14
P(0) =
Φ C
Φ 0
=
Φ(0,t)
Φ 0
=
TC - TF
T0 - TF
= 0,9

P(L) =
Φ L
Φ0
=
ΦL
ΦC
ΦC
Φ0
=
ΦL
ΦC
=
TL - TF
TC - TF
=
x
L
= 1{ }= 0,27 = 0,90 x 0,27 = 0,249
Cilindro infinito ⇒ Bi =
hcF R
k
=
30x 0,05
0,5
= 3 ; Bi-1
= 0,333 ; Fo =
α t
R2
=
5.10-7
x 1800
0,052
= 0,36
C(0) =
Φ C
Φ 0
=
Φ(0,t)
Φ 0
=
TC - TF
T0 - TF
= 0,47 (eje)
C(R ) =
Φ(R,t)
Φ0
=
Φ R
Φ C
Φ C
Φ0
=
Φ R
Φ C
= 0,36 = 0,47 x 0,36 = 0,1692
Temperatura mínima del cilindro (en el centro geométrico)
Φ mín
Φ 0
= P(0) C(0) = 0,90 x 0,47 = 0,423 =
TC - TF
T0 - TF
; TC = Tmín = 0,423 (20 - 500 ) + 500 = 297ºC
Temperatura máxima del cilindro (en la circunferencia de las bases)
Φ máx
Φ0
= P(L) C(R) = 0,249 x 0,1692 = 0,042 =
Tmáx - TF
T0 - TF
; Tmáx = 0,042 (20 - 500) + 500 = 480ºC
b) Calor absorbido en ese tiempo: Q = θ ρ cp V(T0 − TF )
Placa: Bi2 Fo = 4,82 x 0,14 = 3,2256 ⇒ Φplaca =
Q
Q0
= 0,32
Cilindro: Bi2 Fo = 32 x 0,36 = 3,24 ⇒ Φcilindro =
Q
Q0
=0,7
θ = Φplaca + Φcilindro (1 - Φ placa ) = 0,32 + 0,7 (1 - 0,32 ) = 0,796
Q = θ ρ cp V(T0 − TF ) =
α =
k
ρ cp
; ρ cp =
k
α
D = 0,1 m ; h = 0,16 m
= θ k
α
V(T0 − TF ) =
= 0,796
0,5 (W/mºC)
5.10-7 (m2 /seg)
π x 0,052 m2 x 0,16 m3 (20 - 500) ºC = 480135 W seg =
480135 Kcal
4,186 x1000
= 114,7 Kcal
*****************************************************************************************
III.24.- Una pared de ladrillo de 25 cm de espesor, que se encontraba inicialmente a una temperatura T0 = 90°C,
reduce bruscamente la temperatura de una de sus caras hasta alcanzar una temperatura Ts = 15°C.
Determinar:
a) La temperatura del plano que se encuentra a 10 cm de la superficie, cuando hayan transcurrido 2 horas
b) El calor disipado a través de la pared en ese tiempo
c) Suponiendo ahora que las paredes del ladrillo se ponen en contacto con un fluido exterior a 15°C, siendo el
coeficiente de película hcF = 10 Kcal/h.m2°C
c-1) La temperatura en un plano situado a 10 cm del plano central, al cabo de 2 horas
c-2) La temperatura en un plano situado a 10 cm de la superficie, al cabo de 2 horas
c-3) El calor disipado a través de la pared en dicho tiempo
Datos del ladrillo: cpF = 0,2 Kcal/kgºK ; α = 0,00185 m2
/hora ; ρ = 1602 kg/m3
; k = 0,6 Kcal/h.mºK
_________________________________________________________________________________________
RESOLUCIÓN
a) Temperatura del plano que se encuentra a 10 cm de la superficie, cuando hayan transcurrido 2 horas.- El va-
riar bruscamente la temperatura de una de las caras, implica una condición de contorno isotérmi-
ca; hay que comprobar que el nº de Fo 1
Fo = α t
L2
=
0,00185 (m2
/hora)x 2 horas
0,252
m2
= 0,0592
Con este valor de Fo se puede tratar también como sólido semiinfinito, ya que la condición Fo1
implica que lejos de la superficie la temperatura de la cara opuesta permanece, prácticamente, a
la temperatura inicial de 90ºC. También se puede comprobar con: xt = π α t
Para: x = 0,25 m ; u = x
2 α t
=
0,25
2 0,00185 x 2
= 2,054 ⇒ G(2,054) = 0,99615
T(0,25 m ; 2 horas) = 15 + 0,99615 x (90 - 15) = 89,71ºC ≅ 90ºC (SI sólido semiinfinito)

Para: x = 0,10 m ; u =
x
2 α t
=
0,1
2 0,00185 x 2
= 0,822 ⇒ G(0,822) = 0,754
T(0,1 m ; 2 horas) = 15 + 0,754 x (90 - 15) = 71,55ºC
b) Calor disipado a través de la pared en ese tiempo
Q(2 horas) = 2
π
(Ts - T0) ρ cpF k t = 2
π
(15 - 90) 1602 x 0,2 x 0,6 x 2 = 1659,4 Kcal
m2
¿PUEDE TRATARSE COMO SOLIDO INFINITO CON CONDICIÓN DE CONTORNO ISOTÉRMICA? Si se in-
terpreta que la temperatura de la cara opuesta prácticamente no se modifica en un cierto tiempo, esa cara se puede
considerar como adiabática en ese tiempo, por lo que se puede duplicar el espesor a 2 L = 50 cm; las condiciones de
contorno implican que la temperatura en x=0 y x=2 L tienen que ser la misma.
La x se toma a partir del origen de coordenadas M, y por la simetría de la distribución de temperaturas, su valor está
entre 0 y 0,25, mientras que L* =0,5 m
Fo =
α t
L*2
=
0,00185 (m2
/h) x 2 horas
0,52
m2
= 0,0148
ξ = x
L*
=
0,1
0,5
= 0,2
⎫
⎬
⎪
⎭
⎪
⇒
T10 cm - TF
T0 - TF
= 0,735
T10 cm = 15 + 0,735 (90 - 15) = 70,13ºC
⎧
⎨
⎪
⎩
⎪
b) Calor disipado a través de la pared en ese tiempo:
Qx=0
A
= 0,265
Q0
A
=
Q 0
A
= L ρ cpF (T0 - TF ) = 0,25 x 1602 x 0,2 x (90 - 15) = 6008
Kcal
m2
= 1592
Kcal
m2
c) Las paredes del ladrillo se ponen en contacto con un fluido exterior a 15°C, siendo el coeficiente de película
hcF = 10 Kcal/h.m2°C
c-1) Temperatura en el plano situado a 10 cm del plano central, al cabo de 2 horas
Sólido infinito con condición de contorno de convección
Bi =
hcFL
k =
10 x 0,125
0,6 = 2,083 ; Bi-1 = 0,48
Fo =
α t
L2 =
0,00185 x 2
0,1252 = 0,2368 ; ξ =
x
L
=
0,1
0,125
= 0,8
Bi2Fo = 2,0832 x 0,2368 = 1,0274
Temperatura en el plano medio:
TC (2 h ) - TF
T0 - TF
=
TC( 2 h ) - 15
90 - 15
=
Bi-1
= 0,48
Fo = 0,2368
= 0,71 ⇒ TC(2 h ) = 15 + 0,71 (90 - 15) = 68,25ºC
Temperatura en un plano situado a 10 cm del plano central, al cabo de 2 horas
€
T10cm; 2h - TF
TC - TF
=
Bi-1
= 0,48
x
L
=
0,1
0,125
= 0,8
= 0,64 ⇒ T10cm;2h = 15 + 0,64 (68,25 - 15) = 49,08ºC
Temperatura en un plano situado a 10 cm de la superficie, al cabo de 2 horas
€
T25cm; 2h - TF
TC - TF
=
Bi-1
= 0,48
x
L
=
0,025
0,125
= 0,2
= 0,98 ⇒ T25cm;2h = 15 + 0,98 (68,25 - 15) = 67,18ºC
c-3) Calor disipado a través de la pared en dicho tiempo
Q
Q0
=
Bi2Fo = 1,0274
Bi = 2,083
= 0,25 ;
Q
A
= 0,25 x 6007,5 = 1501,7 Kcal
m2

III.25.- Una barra cilíndrica de acero inoxidable 18-8, de 20 cm de diámetro, se calienta a 1000°C y a continua-
ción se enfría en un baño de aceite a 50°C, en el que el coeficiente de película es hC= 500 Kcal/hm2°C.
Determinar: a) El tiempo que transcurrirá hasta que el eje del cilindro alcance una temperatura de 250°C.
b) La temperatura que se alcanzará en r = 0,05 m al cabo de ese tiempo.
Datos del acero: ρ= 7.817 kg/m3 ; cp= 0,11 Kcal/kg°C ; k= 22,5 Kcal/m.hº.C ; α = 0,01598 m2/hora
_________________________________________________________________________________________
RESOLUCIÓN
hC /k = 500/22,5 = 22,2 ⇒ c.c. convección
Bi = hCR/k = 500 x 0,1/22,5 = 2,22 0,1 ⇒ c.c. convección
⎧
⎨
⎩
Bi-1
= 0,45 ;
T( 0,t) - TF
T0 - TF
=
250 - 50
1000 - 50
= 0,21
a) Tiempo que transcurrirá hasta que el eje del cilindro alcance una temperatura de 250°C:
Fo = 0,60 = α t
R2 ⇒ t =
0,6 R2
α
=
0,6 x 0,12 m2
0,01598 (m2 /h)
= 0,375 h = 22 min 30 seg
b) Temperatura que se alcanzará en r = 0,05 m al cabo de ese tiempo:
T(0,5; t) - TF
T(0,t) - TF
=
r
R
=
0,05
0,1
= 0,5 = 0,85 ⇒ T(0,5; t) = 50 + 0,85 (250 - 50) = 220ºC
*****************************************************************************************
III.26.- Una barra de acero de (5 x 4) cm2 de sección transversal y 4 metros de longitud, k= 60 W/m°C, α= 1,6.10-5
m2/seg está inicialmente a una temperatura uniforme de 225°C. A continuación se la introduce en un medio con-
vectivo con un coeficiente hCF= 500 W/m2°C, estando el ambiente a 25°C.
Calcular la temperatura en el centro para un tiempo de 2 minutos después de haberla sometido a la acción del
medio convector exterior
_________________________________________________________________________________________
RESOLUCIÓN
Consideraremos al sólido como intersección de 3 placas infinitas, de las que una de ellas es muy larga respecto a las
otras dos (4 metros), por lo que se puede considerar la transmisión de calor en dos direcciones correspondientes a la
sección transversal, definiendo una temperatura adimensional de la forma:
Φ(x,t)
Φ0
=
T(x,y,t) - TF
T0 - TF
= P(x) P(y)
Llamamos
ΦC1
Φ0
a la temperatura adimensional en el plano central para una placa de espesor 2 L1 = 5 cm
Fo = α t
L1
2
=
(1,6.10-5 ) (2 x 60)
2,52.10−4
= 3,1
Bi-1 =
k
hcF L1
=
60
500 x 2,5.10-2
= 4,8
⎫
⎬
⎪
⎭⎪
⇒
ΦC1
Φ0
= 0,58
Llamamos
Φ C2
Φ0
a la temperatura adimensional en el plano central para una placa de espesor 2 L2 = 4 cm
Fo =
α t
L2
2 =
(1,6.10-5 ) (2 x 60)
22.10−4 = 4,8
Bi-1 =
k
hcF L2
=
60
500 x 2.10-2
= 6
⎫
⎬
⎪
⎭⎪
⇒
ΦC2
Φ0
= 0,45
Para las dos dimensiones se tiene:
Φ(x,r,t)
Φ 0
= P(x) C(r) = 0,58 x 0,45 = 0,26
Temperatura en el centro para el problema bidimensional:

0,26 =
T(0,t) - TF
T0 - TF
=
T(0,t) - 25
225 - 25
⇒ T(0,t) = 25 + 0,26 x (225 - 25) = 77ºC
***********************************************************************************************************
III.27.-Una lata cilíndrica de 5 cm de longitud y 30 cm de diámetro contiene un determinado producto a una tem-
peratura uniforme de 15ºC, y se calienta mediante una corriente de aire a 160ºC.
Sabiendo que: α = 0,04 m2/hora , k = 4 Kcal/h.m.ºC ; hC(aire-lata) = 20 Kcal/h.m2.ºC
Determinar
a) El tiempo necesario para que la temperatura de cualquier punto del producto sea, por lo menos, de 138ºC
b) La temperatura en el centro de la base de la lata en ese instante
c) El calor absorbido por la lata durante los 36 primeros minutos del calentamiento
_________________________________________________________________________________________
RESOLUCIÓN
a) Tiempo necesario para que la temperatura del centro sea de 138ºC, ya que la de los demás puntos será superior
El diámetro del cilindro es bastante mayor que la altura, por lo que la transmisión de calor se realizará, principalmen-
te, en sentido axial, y sería equivalente al de transmisión de calor por conducción en régimen transitorio en una placa
de 5 cm de espesor y longitud infinita.
h C /k = 20/4 = 5 ⇒ c.c. convección
Bi = hC L/k = 20 x 0,025/ 4 = 0,125 0,1 ⇒ c.c. convección
⎧
⎨
⎩
Bi-1
= 8 ; Fo = 16
En el plano medio de la placa hay que conseguir una temperatura, T = 138ºC
De los gráficos de Heysler se obtiene la temperatura adimensional en el plano medio de la placa:
TC - TF
T0 - TF
=
138 - 160
15 - 160
= 0,1517
Fo = 16 = α t
L2
; t =
16 x 0,0252 m2
0,04 (m2
/h)
= 0,25 h = 15 min
que es un valor bastante aproximado.
b) Temperatura en el centro de la base de la lata en ese instante
b.1.- Conducción bidimensional
Placa infinita ⇒ Bi =
h cF L
k
=
20 x 0,025
4
= 0,125 ⇒ Bi-1
= 8 ;
h cF
k
= 5
1
m
; Fo = 16
P(0) =
Φ( 0,t)
Φ 0
=
TC - TF
T0 - TF
= 0,1517
P(L) =
Φ(L,t)
Φ0
= Φ
ΦC
Φ C
Φ 0
= Φ
Φ C
=
T( L) - TF
TC - TF
= { x
L
= 1} = 0,93 = 0,93 x 0,1517 = 0,1411
Cilindro infinito ⇒ Bi =
hcF R
k
=
20 x 0,15
4
= 0,75 ;Bi-1
= 1,333 ; Fo = α t
R 2
=
0,04 m2 /h x 0,25 h
0,152
m2
= 0,444
C(0) =
Φ(0,t)
Φ0
=
TC - TF
T0 - TF
= 0,58 (eje)
C(R ) =
Φ(R,t)
Φ 0
=
Φ( R,t)
Φeje
Φ eje
Φ 0
=
Φ( R,t)
Φ eje
= 0,62 = 0,62 x 0,58 = 0,3596
Temperatura en el centro de la base: P(L) C(0) = 0,1411x 0,58 = 0,0818
Φx
Φ0
=
Tx - TF
T0 - TF
= 0,0818 ; Tx = TF + 0,0818 (T0 - TF) = 160 + 0,0818 x (15 - 160) = 148,14ºC
b.2.- Conducción unidimensional
€
Φ(L,t)
Φ(0,t)
=
T(L) - TF
TC - TF
=
T(L) -160
138 -160
= 0,93 ⇒ T(L) = 139,5°C

c) Calor absorbido por la lata durante los 36 primeros minutos del calentamiento. Conducción bidimensional
Q0 = V ρ cpF ΔT = V
k
α
ΔT = (π R 2
L)
k
α
ΔT = (π x 0,152
x 0,05) m3 4 ( Kcal/h mºC)
0,04 (m2
/h)
(160 - 15)ºC = 51,24 Kcal
Q = θ ρ cp V(T0 − TF )
Placa:
Fo =
α t
L2
=
0,04 (m2
/h) (36/60)
0,0252
= 38,4
Bi2 Fo = 0,1252 x 38,4 = 0,6 ; Bi = 0,125
⎧
⎨
⎪
⎩
⎪
⎫
⎬
⎪
⎭
⎪
Gröber
⎯ →⎯⎯ Φplaca =
Q placa
Q
= 0,97
Cilindro:
Fo = α t
R 2
=
0,04 (m2
/h) (36 /60)
0,152
= 1,07 ; Bi =
h CR
k
=
20 x 0,15
4
= 0,75
Bi2 Fo = 0,752 x 1,07 = 0,6
⎧
⎨
⎪
⎩⎪
⎫
⎬
⎪
⎭⎪
Gröber
⎯ →⎯⎯ Φ cilindro =
Qcilindro
Q
= 0,74
θ = Φplaca + Φcilindro (1 - Φ placa ) = 0,97 + 0,74 (1 - 0,97) = 0,992
Q = 0,992 x Q0 = 50,83 Kcal
*****************************************************************************************
III.28.- Un cilindro de acero, k = 25 W/m°C, α = 1,11 x 10-5 m2/seg, tiene 10 cm de diámetro y 20 cm de longitud, y
se encuentra inicialmente a una temperatura de 850°C.
a) Tiempo que debe transcurrir para que el cilindro alcance una temperatura de 350°C en su centro de gravedad.
Se le sumerge en agua a 20°C y sabiendo que el coeficiente de transmisión de calor por convección en cualquier
punto del cilindro es de 1725 W/m2°C, determinar:
a) El tiempo que debe transcurrir para que el cilindro alcance una temperatura de 350°C en su centro de gra-
vedad.
b) La temperatura que en ese instante se alcanzará en la circunferencia de sus bases.
c) La temperatura en el centro de la base
d) La temperatura en la mitad de la pared cilíndrica
_________________________________________________________________________________________
RESOLUCIÓN
Distribución detemperaturas:
Φ(x,r,t)
Φ 0
= P(x) C(r)
€
Placa plana ∞ de espesor 20 cm: Bi =
hcF L
k
=
1725 x 0,1
25
= 6,9 ; Bi-1
= 0,145
Temperatura adimensional en el plano medio:
P(0) =
Φ(0,t)
Φ0
=
T(0,t) - TF
T0 - TF
= 350 - 20
850 - 20
= 0,3975
Entrando en la gráfica de Heysler para placa plana ∞ con P(0) = 0,3975 y
Bi = 6,9 se obtiene un nº de Fo = 0,60, por lo que:
Fo =
α t
L2
=
1,11.10-5 t
0,12
= 1,11.10-3 t = 0,60 ⇒ t = 552 seg
Cilindro infinito de radio 5 cm ⇒ Bi =
h cF R
k
=
1725x 0,05
25
= 3,45 ; Bi-1
= 0,29
La temperatura adimensional en el eje del cilindro es la misma que en el centro de la placa = C(0) = 0,3975
Entrando en la segunda gráfica de Heysler para cilindro infinito con C(0) = 0,3975 y Bi= 3,45 se obtiene un nº de Fo
= 0,38, por lo que:
Fo =
α t
R2
=
1,11.10-5 t
0,052
= 0,38 ⇒ t = 85,5 seg
que es el mínimo de los dos tiempos, y que nos da una referencia para iterar
En lo que sigue, mantendremos los resultados obtenidos en las Tablas
∞
∞
∞
∞

PLACA CILINDRO
€
Para t = 80 segundos ⇒
Placa ∞
Bi-1
= 0,145
Fo =
α t
L2 =
1,11.10-5
x 80
0,12 = 0,088
⎧
⎨
⎪
⎩
⎪
⇒
TC - TF
T0 - TF
= P(0) = 0,98
Cilindro ∞
Bi-1
= 0,29
Fo =
α t
R2 =
1,11.10-5
x 80
0,052 = 0,35
⎧
⎨
⎪
⎩
⎪
⇒
TC - TF
T0 - TF
= C(0) = 0,4336
⎧
⎨
⎪
⎪
⎪
⎩
⎪
⎪
⎪
TC - TF
T0 - TF
= 0,98 x 0,4336 = 0,425 ; TC = 20 + (850 - 20) x 0,425 = 372,75ºC
€
Para t = 90 segundos ⇒
Placa ∞
Bi-1
= 0,145
Fo =
α t
L2 =
1,11.10-5
x 90
0,12 = 0,099
⎧
⎨
⎪
⎩
⎪
⇒
TC - TF
T0 - TF
= P(0) = 0,974
Cilindro ∞
Bi-1
= 0,29
Fo =
α t
R2
=
1,11.10-5
x 90
0,052
= 0,4
⎧
⎨
⎪
⎩
⎪
⇒
TC - TF
T0 - TF
= C(0) = 0,373
⎧
⎨
⎪
⎪
⎪
⎩
⎪
⎪
⎪
€
TC − TF
T0 − TF
= 0,974 × 0,373 = 0,364 ; TC = 20 + (850 − 20)× 0,364 = 322,12°C
Para alcanzar la temperatura de 350ºC en el (c.d.g.) del cilindro, son necesa-
rios, aproximadamente, 84,5º segundos.
b) Temperatura que en ese instante se alcanza en la circunferencia de sus
bases.
P(L) =
Φ
Φ0
=
Φ
ΦC
ΦC
Φ 0
=
ΦC
Φ0
= P(0) = 0,3975
Φ
ΦC
= ξ = 1( ) = 0,22
= 0,22 x 0,3975 = 0,087
C(R) =
Φ
Φ 0
=
Φ
ΦC
ΦC
Φ0
=
ΦC
Φ0
= C(0) = 0,3975
Φ
ΦC
= ξ = 1( ) = 0,34
= 0,34 x 0,3975 = 0,1351
Φ circ. bases
Φ 0
=
Tcirc. bases - TF
T0 - TF
=
Tcirc. bases - 20
850 - 20
= P(L) C(R) = 0,087 x 0,1312 = 0,01147 ⇒ Tcirc . bases = 29,5ºC
c) Temperatura en el centro de la base
Φcentro base
Φ0
=
Tcentro base - TF
T0 - TF
=
Tcentro base- 20
850 - 20
= P(L) C(0) = 0,087 x 0,3975 = 0,03458 ⇒ Tcentro base = 48,7ºC
d) Temperatura en la mitad de la pared cilíndrica
Φ(0,R; t)
Φ 0
=
T(0,R; t) - TF
T0 - TF
=
T(0,R;t) - 20
850 - 20
= P(0) C(R) = 0,3975 x 0,1351 = 0,0537 ⇒ T(0,R; t) = 64,6ºC
*****************************************************************************************

III.29.-Se desea hacer un tratamiento térmico de unas bolas de acero al niquel-cromo 15%Cr, 10%Ni, y para ello
se calientan a una temperatura de 875°C. A continuación se templan en aceite que se encuentra a la temperatura
de 120°C, sometiéndolas a un movimiento de agitación que se puede considerar dotado de una velocidad media
u=1 m/seg. Las bolas tienen un diámetro d=3 cm.Con estos datos, determinar:
a) El coeficiente de transferencia de calor por convección entre las bolas y el aceite cuando TpF = 200ºC
b) El tiempo necesario para que las bolas alcancen en su superficie la temperatura de 200 °C
c) La cantidad de calor eliminada, en cada bola, en el intervalo de tiempo del apartado anterior.
RESOLUCIÓN
________________________________________________________________________________________
a) Coeficiente de transferencia de calor por convección entre las bolas y el aceite
Flujo en convección forzada sobre esferas
Nu = 2 + (0,4 Re + 0,06 Re2/ 3 ) Pr0,4 ηF
ηpF
4 en el intervalo:
3,5 Re 80000 ; 0,7 Pr 380
1
ηF
ηpF
3,2 (Parámetros a TF )
⎧
⎨
⎩
Propiedades del aceite de motor sin usar a 120ºC:
cp = 2307
J
kgºC
; k = 0,135
W
mºK
; α = 7,1.10 -8 m2
seg
; ν = 12,4.10 -6 m2
seg
; η = 10,3.10-3 Nseg
m2
; Pr = 175
La temperatura de referencia, para hallar ηpF será la media de la pared durante el tratamientor: 875 + 200
2
= 537ºC
Tomaremos ηpF a la máxima temperatura disponible: 4,5.10-3
(N.seg/m2
)
Acero a 20ºC y composición (15% Cr ; 10% Ni):
k = 19 w/mºC ; ρ = 7865 kg/m3
cp = 460 J/kgºC ; α = 0,526.10-6 m2 /seg
⎧
⎨
⎩
Re =
d u
ν
=
0,03 m x 1 m/seg
12,4.10-6 m2 /seg
= 2500
Nu = 2 + (0,4 Re + 0,06 Re2/ 3 ) Pr0,4 ηF
ηpF
4 = 2 + (0,4 2500 + 0,006 x 25000,66 ) 1750,4 10,3
4,5
4 = 298
h cF =
k Nu
d
=
0,135 x 298
0,03
= 1340
W
m2 ºC
b) Tiempo necesario para que las bolas alcancen en su superficie la temperatura de 200 °C
Cuando las bolas alcanzan en la superficie la temperatura de 200ºC, la temperatura en el centro es:
T - TF
TC - TF
=
200 - 120
TC - 120
=
Bi =
h cF R
k
=
1340 (0,03/2)
19
= 1,0579
Bi-1 = 0,945 ; ( r/R) = 1
= 0,63
Heysler
⎯ →⎯ TC = 247 C
y el tiempo necesario para lograrlo:
Fo =
α t
R2 =
0,526.10-5 t
0,0152 = 0,0233 t =
TC - TF
T0 - TF
=
247 - 120
875 - 120
= 0,1682
Bi-1 = 0,945
⎧
⎨
⎩
⎫
⎬
⎭
= 0,8 ⇒ t =
0,8
0,0233
= 34,4 seg
c) Calor eliminado en cada bola en el intervalo de tiempo del apartado anterior.
Q0 =
4
3
π R 3
ρ cp (T0 - TF ) =
4 π
3
(
0,03
2
)3
m3
x 7865
kg
m3
x 0,46
kJ
kgºC
(875 - 120)ºC = 38,61 kJ
Q
Q0
=
Bi2
Fo = 0,8954
Bi = 1,058
= 0,90 ⇒ Q = 0,90 Q0 = 0,90 x 38,61 kJ = 34,75 kJ
*****************************************************************************************
III.30.- Un cilindro de acero de k = 25 Kcal/h.m.ºC, y difusividad térmica α = 0,05 m2/hora, tiene 50 mm de diá-
metro y 65 mm de longitud, y está inicialmente a una temperatura de 650ºC. A continuación se enfría en un fluido
que se encuentra a 90ºC con un coeficiente de película hcF = 150 Kcal/hm2ªC.
Determinar
a) La temperatura en el centro del cilindro y en la circunferencia de las bases al cabo de 2 minutos
b) El tiempo necesario para que en el centro del cilindro se alcancen 100ºC.

_________________________________________________________________________________________
RESOLUCIÓN
El valor de hcF/k = 150/25 = 6 ⇒ condición de contorno de convección
a) Temperatura en el centro del cilindro al cabo de 2 minutos
Cilindro ∞ :
Bi =
hcF R
k
=
150 x 0,025
25
= 0,15 ; Bi-1 = 6,67
Fo = α t
R2 =
0,05 (m2 /h) x 120 seg x (1 h/3600 seg)
0,0252 m2 = 2,67
⎧
⎨
⎪
⎩
⎪
⇒ C(0) =
Φ(0,t)
Φ0
=
TC - TF
T0 - TF
= 0,43
Placa ∞ :
Bi =
h cFL
k
=
150 (0,065/2)
25
= 0,195 ; Bi-1 = 5,128
Fo = α t
L2 =
0,05 (m2 /h) x 120 seg x (1 h/3600 seg)
(0,065/2) m2 = 1,58
⎧
⎨
⎪
⎩
⎪
⇒ P(0) =
Φ(0,t)
Φ0
=
TC - TF
T0 - TF
= 0,75
TCentro cilindro - TF
T0 - TF
= 0,43 x 0,75 = 0,3225 ; TC = 90 + (650 - 90) x 0,3225 = 270,6ºC
Temperatura en la circunferencia de las bases al cabo de 2 minutos
Cilindro ∞
Bi-1
= 6,67
r/R = 1
⎧
⎨
⎩
⇒
Φ(R,t)
Φ(0,t)
=
T(0,t) - TF
T0 - TF
= 0,94
C(R) =
Φ(R,t)
Φ0
=
Φ(R,t)
Φ(0,t)
Φ(0,t)
Φ0
= 0,94 x 0,43 = 0,4042
Placa ∞
Bi-1
= 5,128
x/L = 1
⎧
⎨
⎩
⇒
Φ(L,t)
Φ(0,t)
=
TL - TF
T0 - TF
= 0,90
P(L) =
Φ(L,t)
Φ0
=
Φ(L,t)
Φ(0,t)
Φ(0,t)
Φ0
= 0,90 x 0,75 = 0,675
Φmín
Φ0
=
Tcircunf. bases - TF
T0 - TF
= P(L) C(R) = 0,675 x 0,4042 = 0,2728
Tcircunf. bases = 90 + 0,2728 x (650 - 90) = 242,77ºC
b) Tiempo necesario para que en el centro del cilindro se alcancen 100ºC. Hay que iterar
1ª iteración: Supongamos que esta temperatura se alcanza al cabo de 6 minutos = 360 seg
Cilindro ∞
Bi-1
= 6,67
Fo = α t
R 2
=
0,05(m2
/h) x 360 seg x (1 h/3600 seg)
0,0252
m2
= 8
⎧
⎨
⎪
⎩
⎪
⇒ C(0) =
Φ (0,t)
Φ 0
=
TC - TF
T0 - TF
= 0,09
Placa ∞
Bi-1
= 5,128
Fo =
α t
L2
=
0,05(m2
/h ) x 360 seg x (1 h/3600 seg)
(0,065/2)2
m2
= 4,73
⎧
⎨
⎪
⎩
⎪
⇒ P(0) =
Φ(0,t)
Φ 0
=
TC - TF
T0 - TF
= 0,44
TC2
- TF
T0 - TF
= 0,09 x 0,44 = 0,0396 ; TC2
= 90 + (650 - 90) x 0,0396 = 112,18ºC
2ª iteración: Supongamos que esta temperatura se alcanza al cabo de 8 minutos = 480 seg
Cilindro ∞
Bi-1
= 6,67
Fo =
α t
R 2
=
0,05(m2
/h) x 480 seg x (1 h/3600 seg)
0,0252
m2
= 10,67
⎧
⎨
⎪
⎩
⎪
⇒ C(0) =
Φ (0,t)
Φ 0
=
TC - TF
T0 - TF
= 0,042
Placa ∞
Bi-1
= 5,128
Fo =
α t
L2
=
0,05(m2
/h ) x 480 seg x (1 h/3600 seg)
(0,065/2 )2
m2
= 6,31
⎧
⎨
⎪
⎩
⎪
⇒ P(0) =
Φ(0,t)
Φ 0
=
TC - TF
T0 - TF
= 0,32
TC3
- TF
T0 - TF
= 0,042 x 0,32 = 0,01344 ; TC3
= 90 + (650 - 90) x 0,01344 = 97,52ºC

Representando estos valores, se obtiene un tiempo próximo a los
7,6 minutos
Mediante programa informático sale:
t = 463,7 seg = 7,7 minutos
Focilindro = 10,31 ; Foplaca = 6,1
*****************************************************************************************
III.31.- Una lata cilíndrica de 5 cm de longitud y 30 cm de diámetro contiene un determinado producto a una
temperatura uniforme de 15ºC, y se calienta mediante una corriente de aire a 160ºC.
€
Sabiendo que: α = 0,04 (m2
/ hora) ; k = 4 (Kcal / hmºC) ; hC(aire-lata) = 20 (Kcal / hm2
º C) , determinar:
a) El tiempo necesario para que la temperatura de cualquier punto del producto sea como mínimo de 138ºC
b) La temperatura en el centro de la base de la lata en ese instante
c) El calor absorbido por la lata durante los 36 primeros minutos del calentamiento
_________________________________________________________________________________________
RESOLUCIÓN
El que la temperatura mínima sea de 138ºC implica que sea el c.d.g. quien la tiene que alcanzar, por cuanto los de-
más puntos del cilindro estarán a una temperatura mayor.
a) Tiempo necesario para que la temperatura de cualquier punto del producto sea como mínimo de 138ºC
Por ser el diámetro del cilindro muy superior a la altura, la transmisión de calor se puede interpretar se realiza princi-
palmente en sentido axial, por lo que sería equivalente al de transmisión de calor por conducción en régimen transito-
rio en una placa infinita de 5 cm de espesor.
FLUJO UNIDIMENSIONAL
Bi = hC L/k = 20 x 0,025/4 = 0,125 0,1 ⇒ c.c. convección
⎧
⎨
⎩
La temperatura de 138ºC se consigue en el plano medio de la placa.
PLACA.- Haciendo uso de la 2ª gráfica de Heysler se obtiene:
Temp. adimensional en el plano medio de la placa:
TC- TF
T0 - TF
=
138 - 160
15 - 160
= 0,1517
Bi-1 = 8
⎫
⎬
⎭
⇒ Fo = 16 =
α t
L2
Tiempo que se tarda en alcanzar la temperatura de 138ºC: t = 16 L2
α
=
16 x 0,0252 m2
0,04 m2 /hora
= 0,25 horas = 15 minutos
b) Temperatura en el centro de la base de la lata en ese instante
T(x,t) - TF
TC - TF
=
T(x,t) - 160
138 - 160
= 0,93 ; T(x,t) = 160 + 0,93 (138 - 160) = 139,5ºC
c) Calor absorbido por la lata durante los 36 primeros minutos del calentamiento
Q0 = m cp ΔT = V ρ cp ΔT = ρ cp = k
α
= V k
α
ΔT = (π R2 L) k
α
ΔT =
= (π x 0,152 x 0,05) m3 4 (Kcal/h mºC)
0,04 ( m2 /hora)
(160 - 15) ºC = 51,24 Kcal
Fo =
α t
L2 =
0,04 (m2 /h) x (36/60) h
0,0252 m2 = 38,4 ; Bi2 Fo = 0,1252 x 38,4 = 0,6
y entrando en la gráfica de Gröber se obtiene:
Q(t)
Q0
= 0,98 ; Q = 0,98 x 51,24 Kcal = 50,21 Kcal
FLUJO BIDIMENSIONAL: En la intersección de un cilindro infinito y una placa infinita; la temperatura en el cen-
tro de la base de la lata se obtiene tomando como referencia el tiempo obtenido en el apartado (a), t = 0,25 horas:
Cilindro(0,25h):
Bi = 0,75 ; Bi-1 = 1,333
Fo = α t
R 2
=
0,04 (m2 /h) x 0,25 h
0,152 m2
= 0,447
⇒ C(0) =
Φ(0,t)
Φ0
=
TC - TF
T0 - TF
= 0,70
⎧
⎨
⎪
⎩
⎪
C(R) =
ΦR
Φ0
=
ΦR
ΦC
ΦC
Φ0
=
ΦR
ΦC
= (ξ = 1) = 0,71
⎧
⎨
⎩
⎫
⎬
⎭
= 0,71 x 0,7 = 0,497
⎧
⎨
⎪
⎪⎪
⎩
⎪
⎪
⎪

Placa (0,25h):
Bi = 0,125 ; Bi-1
= 8
Fo =
α t
L2
=
0,04 (m2
/h) x 0,25 h
0,0252
m2
= 16
⇒ P(0) =
Φ(0,t)
Φ 0
=
TC - TF
T0 - TF
= 0,15
⎧
⎨
⎪
⎩
⎪
P( L) =
Φ L
Φ 0
=
Φ L
Φ C
ΦC
Φ 0
=
Φ L
ΦC
= (ξ = 1) = 0,94
⎧
⎨
⎩
⎫
⎬
⎭
= 0,94 x 0,15 = 0,141
⎧
⎨
⎪
⎪
⎩
⎪
⎪
Temperatura en el c.d.g. al cabo de 0,25 h :
Tc.d.g = 160 + (15 - 160) P(0) C(0) = 160 - (145 x 0,15 x 0,70) = 144,78ºC
En una segunda aproximación lo haremos con un tiempo t = 0,2 horas
Cilindro(0,2 h):
Bi = 0,75 ; Bi-1
= 1,333
Fo =
α t
R2
=
0,04 (m2
/h) x 0,2 h
0,152
m2
= 0,355
⇒ C(0 ) =
Φ(0,t)
Φ0
=
TC - TF
T0 - TF
= 0,79
⎧
⎨
⎪
⎩
⎪
C( R) =
Φ R
Φ 0
=
Φ R
Φ C
Φ C
Φ0
=
Φ R
ΦC
= (ξ = 1) = 0,71
⎧
⎨
⎩
⎫
⎬
⎭
= 0,71 x 0,79 = 0,56
⎧
⎨
⎪
⎪
⎩
⎪
⎪
Placa (0,2h ):
Bi = 0,125 ; Bi-1
= 8
Fo =
α t
L2
=
0,04 ( m2
/h) x 0,2 h
0,0252
m2
= 12,8
⇒ P(0) =
Φ(0,t)
Φ0
=
TC - TF
T0 - TF
= 0,20
⎧
⎨
⎪
⎩
⎪
P(L) =
Φ L
Φ0
=
Φ L
Φ C
Φ C
Φ0
=
Φ L
Φ C
= (ξ = 1) = 0,94
⎧
⎨
⎩
⎫
⎬
⎭
= 0,94 x 0,2 = 0,188
⎧
⎨
⎪
⎪
⎩
⎪
⎪
Temperatura en el c.d.g. al cabo de 0,2 h :
Tc.d. g = 160 + (15 - 160) P(0) C(0) = 160 - (145 x 0,2 x 0,79) = 137,09ºC
ligeramente por debajo de la que pide el problema (138ºC). La solución es
del orden de 0,207 h
Temperatura en el centro de la base:
Tcentro base = 160 + (15 - 160) P(L) C(0) = 160 - (145 x 0,188 x 0,79) = 138,46ºC
que difiere de los 139,5ºC obtenidos como problema unidimensional en 1ºC aproximadamente.
La aproximación unidimensional es bastante buena
*****************************************************************************************
III.32.- Se dispone de una serie de botes de conserva de diámetro d= 12 cm y longitud L= 10 cm, que inicialmente
están a una temperatura de 20°C; en un momento determinado se sumerge uno de ellos en un líquido caliente a
100°C. El coeficiente de película es de 200 W/m2°C. La difusividad térmica del bote y su contenido es de 0,03
m2/hora y su conductividad térmica de 9 W/m°C.
Determinar
a) El tiempo en minutos que deberá tenerse el bote en el líquido caliente con el fin de garantizar que toda su
masa se encuentre por encima de 76°C
b) La temperatura mínima que se registrará en el bote al mantenerlo sumergido 10 minutos. ¿En donde se
tendrá la temperatura máxima y valor de ésta?
c) Si ahora se introducen los botes en el líquido caliente en grupos de 20, perfectamente alineados, confor-
mando un cilindro de 12 cm de diámetro y 2 metros de longitud, hallar los supuestos pedidos en el apartado (b).
_________________________________________________________________________________________
RESOLUCIÓN.- Vamos a determinar el tiempo en minutos que se deberá tener el bote en el líquido caliente con el
fin de garantizar que toda su masa se encuentre por encima de 80°C; para ello calcularemos el tiempo que se tardará
en alcanzar el centro del cilindro, como cilindro infinito y como placa infinita; el menor de los resultados nos orienta-
rá sobre la magnitud de este tiempo.
a) Tiempo en minutos que deberá tenerse el bote en el líquido caliente con el fin de garantizar que toda su masa
se encuentre por encima de 80°C.- No se conoce el líquido, pero la relación hcF/k = 200/9 = 22, nos indica que la
condición de contorno es de convección

€
Cilindro ∞
Bi =
hcFR
k
=
200 (W /m2
°C) 0,06 m
9 (W /m°C)
=1,333 ; Bi-1
= 0,75
Fo =
α t
R2 =
0,03 (m2
/h) t (h)
0,062
m2 = 8,333 t (horas) = 0,1389 t (min)
⎧
⎨
⎪⎪
⎩
⎪
⎪
T( 0,t) - TF
T0 - TF
= TF = 100ºC ; T0 = 20ºC ; T(0,t) = 76ºC =
76 - 100
20 - 100
= 0,3
La temperatura en el centro es T(0,t), y como toda la masa tiene que estar por encima de esta temperatura, el centro
tendrá la mínima temperatura = 76ºC
Con estos resultados entramos en el diagrama de Heysler y encontramos Fo = 0,6, por lo que:
Fo = 0,6 = 0,1389 t ⇒ t = 4,32 minutos
Placa ∞
Bi =
h cFL
k
=
200 (W/m2 ºC) 0,05 m
9 (W/mºC)
= 1,111 ; Bi-1
= 0,9
Fo = α t
L2
=
0,03 (m2
/h) t ( h)
0,052
m2
= 12 t (horas) = 0,2 t (min)
⎧
⎨
⎪
⎩
⎪
en la que se ha considerado un espesor de placa igual a la altura del cilindro = 2 L = 10 cm ⇒ L = 0,05 m
Con estos resultados entramos en el diagrama de Heysler y encontramos que Fo = 1,55, por lo que:
Fo = 1,55 = 0,2 t ⇒ t = 7,75 minutos
y la referencia es (t = 4,32 minutos = 272 seg).
A partir de esta referencia vamos a determinar la temperatura en el centro para diversos tiempos y a partir de su
representación gráfica, y para una temperatura en el centro de 76ºC, hallar la solución:
Para t = 200 seg:
Placa ∞ ⇒ Bi-1 = 0,9 ; Fo =
α t
L2
=
0,03 (m2 /h) (200/3600) (h/seg)
0,052 m2
= 0,67 ⇒
T(0,t) - TF
T0 - TF
= 0,66
Cilindro ∞ ⇒ Bi-1 = 0,75 ; Fo =
α t
R 2 =
0,03 (m2 /h) (200/3600) (h/seg)
0,062 m2 = 0,463 ⇒
T(0,t) - TF
T0 - TF
= 0,51
⎧
⎨
⎪
⎩
⎪
TC - TF
T0 - TF
= 0,66 x 0,51 = 0,3366 ⇒ TC = 100 + 0,3366 (20 - 100) = 73,07ºC
Para t = 225 seg:
Placa ∞ ⇒ Bi-1 = 0,9 ; Fo =
α t
L2
=
0,03 (m2 /h) (225/3600) (h/seg)
0,052 m2
= 0,75 ⇒
T(0,t) - TF
T0 - TF
= 0,62
Cilindro ∞ ⇒ Bi-1 = 0,75 ; Fo =
α t
R2 =
0,03 (m2 /h) (225/3600) (h/seg)
0,062 m2 = 0,52 ⇒
T(0,t) - TF
T0 - TF
= 0,46
⎧
⎨
⎪
⎩
⎪
TC - TF
T0 - TF
= 0,62 x 0,46 = 0,285 ⇒ TC = 100 + 0,285 (20 - 100) = 77,18ºC
Para t = 215 seg:
Placa ∞ ⇒ Bi-1 = 0,9 ; Fo = α t
L2
=
0,03 (m2 /h) (215/3600) (h/seg)
0,052 m2
= 0,7167 ⇒
T(0,t) - TF
T0 - TF
= 0,64
Cilindro ∞ ⇒ Bi-1 = 0,75 ; Fo = α t
R2 =
0,03 (m2 /h) (215/3600) (h/seg)
0,062 m2 = 0,498 ⇒
T(0,t) - TF
T0 - TF
= 0,47
⎧
⎨
⎪
⎩
⎪
TC - TF
T0 - TF
= 0,64 x 0,47 = 0,3008 ⇒ TC = 100 + 0,3008 (20 - 100) = 75,93ºC
Haciendo la representación gráfica de los resultados obtenidos se en-
cuentra que el tiempo que debe tenerse el bote en agua para que toda
su masa se encuentre por encima de 76ºC es de 216 segundos
b) Temperatura mínima que se registrará en el bote al mantenerlo
sumergido 2 y 10 minutos, respectivamente. ¿En donde se tendrá
la temperatura máxima y valor de ésta?
La temperatura mínima se presenta en el centro del cilindro, mien-
tras que la máxima se presenta en los bordes
€
Placa ∞
Bi-1
= 0,9
Fo = 0,2 t = 0,2 x 2 = 0,4
⎧
⎨
⎩
⇒
P(0) =
Φ(0,t)
Φ0
=
T(0,t) - TF
T0 - TF
= 0,82
P(L) =
Φ(L,t)
Φ0
=
x
L
= 1 =
Φ(L,t)
Φ(0,t)
Φ(0,t)
Φ0
= 0,515
⎧
⎨
⎪⎪
⎩
⎪
⎪
Para t = 2 minutos

€
Cilindro ∞
Bi-1
= 0,75
Fo = 0,1389 x 2 = 0,2778
⎧
⎨
⎩
⇒
C(0) =
Φ(0,t)
Φ0
=
T(0,t) - TF
T0 - TF
= 0,73
C(R) =
Φ(R,t)
Φ0
=
r
R
= 1 =
Φ(R,t)
Φ(0,t)
Φ(0,t)
Φ0
= 0,4162
⎧
⎨
⎪⎪
⎩
⎪
⎪
Temperatura mínima y máxima del cilindro al cabo de 2 minutos:
EJE:
Φmínima
Φ0
= P(0) C(0) = 0,82 x 0,73 = 0,5986 ⇒ Tmín = 100 + 0,5986 (20 - 100) = 52,11ºC
PERIFERIA:
Φmáxima
Φ0
= P(L) C(R) = 0,515 x 0,4162 = 0,2143 ⇒ Tmáx = 100 + 0,2143 (20 - 100) = 82,85ºC
Para t = 10 minutos
Placa ∞
Bi-1 = 0,9
Fo = 0,2 t = 0,2 x 10 = 2
⎧
⎨
⎩
⇒
P(0) =
Φ(0,t)
Φ0
=
T(0,t) - TF
T0 - TF
= 0,23
P(L) =
Φ(L,t)
Φ0
= x
L
= 1 =
Φ(L,t)
Φ(0,t)
Φ(0,t)
Φ0
= 0,64 x 0,23 = 0,1432
⎧
⎨
⎪
⎩
⎪
Cilindro ∞
Bi-1 = 0,75
Fo = 0,1389 x 10 = 1,389
⎧
⎨
⎩
⇒
C(0) =
Φ(0,t)
Φ0
=
T(0,t) - TF
T0 - TF
= 0,084
C(R) =
Φ(R,t)
Φ 0
= r
R
= 1 =
Φ(R,t)
Φ(0,t)
Φ(0,t)
Φ0
= 0,047
⎧
⎨
⎪
⎩
⎪
Temperaturas mínima y máxima del cilindro al cabo de 10 minutos:
Φ mínima
Φ0
= P(0) C(0) = 0,23 x 0,084 = 0,01932 ⇒ Tmín = 100 + 0,1932 (20 - 100) = 98,45ºC
Φ máxima
Φ 0
= P(L) C(R) = 0,1432 x 0,047 = 0,006739 ⇒ Tmáx = 100 + 0,006739 (20 - 100) = 99,46ºC
c) Los botes se introducen en el líquido caliente en grupos de 20, perfectamente alineados, conformando un cilin-
dro de 12 cm de diámetro y 2 metros de longitud; los supuestos pedidos en el apartado (b) se calculan como sigue:
El conjunto de los botes se puede considerar como un cilindro infinito al ser la longitud mucho mayor que el diáme-
tro. La temperatura mínima se presenta en el eje del cilindro, mientras que la máxima se presenta en la periferia. Los
datos de partida se pueden tomar del apartado (b) para el cilindro infinito.
Para t = 2 minutos
Cilindro ∞
Bi-1 = 0,75
Fo = 0,1389 x 2 = 0,2778
⎧
⎨
⎩
⇒
C(0) =
Φ(0,t)
Φ0
=
T(0,t) - TF
T0 - TF
= 0,73
C(R) =
Φ(R,t)
Φ 0
= r
R
= 1 =
Φ(R,t)
Φ(0,t)
Φ(0,t)
Φ 0
= 0,4162
⎧
⎨
⎪
⎩
⎪
Temperatura mínima (en el eje) y máxima (en la periferia) del cilindro infinito al cabo de 2 minutos:
Φmínima
Φ0
=
Teje - TF
T0 - TF
= 0,73 ⇒ Teje = 100 + 0,73 (20 - 100) = 41,6ºC
Φmáxima
Φ 0
=
TR - TF
T0 - TF
= 0,4162 ⇒ Tmáx = 100 + 0,4162 (20 - 100) = 66,7ºC
Para t = 10 minutos
Cilindro ∞
Bi-1 = 0,75
Fo = 0,1389 x 10 = 1,389
⎧
⎨
⎩
⇒
C(0) =
Φ(0,t)
Φ0
=
T(0,t) - TF
T0 - TF
= 0,084
C(R) =
Φ(R,t)
Φ 0
= r
R
= 1 =
Φ(R,t)
Φ(0,t)
Φ(0,t)
Φ0
= 0,047
⎧
⎨
⎪
⎩
⎪
Temperatura mínima (eje) y máxima (periferia) del cilindro infinito al cabo de 10 minutos:
Φmínima
Φ0
=
Teje - TF
T0 - TF
= 0,084 ⇒ Teje = 100 + 0,084 (20 - 100) = 93,3ºC
Φmáxima
Φ 0
=
TR - TF
T0 - TF
= 0,047 ⇒ Tmáx = 100 + 0,047 (20 - 100) = 96,25ºC
*****************************************************************************************
III.33.- Una esfera de hierro, maciza, de 10 cm de diámetro, que inicialmente se encuentra a 780ºC, se enfría en
aire en reposo a 20ºC, hasta que en el centro de la esfera se consigue una temperatura de 300ºC; en el instante en
que se alcanza esa temperatura en el centro de la esfera, ésta se introduce en un baño de agua en movimiento, a

la temperatura de 20ºC, de forma que el coeficiente de película es, hcF = 1950 W/m2.ºK. Propiedades físicas de la
esfera, k = 40 W/mºC ; cp = 460 Joules/kgºC ; ρ = 7850 kg/m3
Determinar
a) El tiempo (en minutos) que debe transcurrir para que en el enfriamiento en el aire seco a 20ºC se alcancen
en el centro de la esfera 300ºC
b) La temperatura en el centro de la esfera, a una distancia radial de 2,5 cm y en la periferia, al cabo de 2 mi-
nutos de ser introducida la esfera en el agua. ¿Qué temperatura en el centro de la esfera se tendrá a los 5 minu-
tos?
c) El calor evacuado al medio exterior en todo el proceso.
_________________________________________________________________________________________
RESOLUCIÓN
a) El tiempo que debe transcurrir para que en el enfriamiento en el aire en reposo a 200ºC se alcancen en el cen-
tro de la esfera 300ºC
Convección natural en esferas:Nud = 2 + 0,45 (Grd Pr)0,25
Propiedades del aire a la temperatura media de película salvo el de β que lo es a la del aire:
T =
TpF + TF
2
= 780 + 20
2
= 400 C = 673 K ⇒ k = 0,051 W
m C
; Pr = 0,68 ;
g β
ν2
=
9,8 (1/293)
(62,38.10-6
)2
= 8,595.106
Grd =
g β
ν2
ΔT d3
=
ΔT = 780 - 20 = 760 ºC
d = 0,1 m
= 8,595.106
x 760 x 0,13
= 6,53.10 6
(laminar)
Grd Pr = 6,53.106
x 0,68 = 4,44.10 6
Nu d = 2 + 0,45 Grd Pr4 = 2 + 0,45 4,44.1064
= 22,65
h cF(aire ) =
22,65 x 0,051(W/mºC)
0,1 m
= 11,56
W
m2
ºC
Al ser un enfriamiento en aire, es previsible una condición de contorno con resistencia térmica interna despreciable,
por lo que Bi 0,1. En efecto:
Bi =
h cF (R/3)
k
=
11,56 (W/m2
ºC) (0,05/3) m
40 (W/mºC)
= 0,00481 0,1
Fo =
α t
(R/3)2
= α =
k
ρ cpF
=
40 ( W/mºC)
7850 (kg/m3
) 460 (J/kg ºC)
= 1,106.10 -5 m2
seg
=
1,106.10-5
(m2
/seg) t
(0,05/3)2
= 0,03987 t
T - TF
T0 - TF
= e-Bi Fo
⇒
300 - 20
780 - 20
= 0,3684 = e- 4,81.10-3
x 0,03987 t
= e- 1,92.10-4
t
⇒ t = 9028 segundos
b) La temperatura en el centro de la esfera, a una distancia radial de 2,5 cm y en la periferia, al cabo de 2 minutos
de ser introducida la esfera en el agua.
Como el proceso anterior ha sido un enfriamiento en aire, con resistencia térmica interna despreciable, toda la esfera
se encuentra a T0 = 300ºC, que será la temperatura de partida para el enfriamiento en el líquido.
Temperatura en el centro de la esfera al cabo de 2 minutos:
h cF
k acero
=
1950 (W/m 2
ºC)
40 (W/mºC)
= 48,7 (Parece condición de contorno de convección, pero conviene hallar Fo)
Bi =
h cFR
k
=
1950 (W/m2
ºC) 0,05 m
40 ( W/mºC)
= 2,435 ; Bi-1 = 0,41
Fo = α t
R 2
=
1,106.10-5
(m2
/seg) 120 seg
0,052
m2
= 0,5308
⎧
⎨
⎪
⎩
⎪
⇒
T(0,t) - TF
T0 - TF
= 0,11
T(0,t) - 20
300 - 20
= 0,11
⎧
⎨
⎪
⎩
⎪
⇒ T(0,t) = 50,8ºC
Al ser Fo 1 puede ser también condición de contorno isotérmica; sabemos que esta condición realmente empieza a
ser significativa para Fo 0,5, por lo que también se puede intentar por este método.
Empleando la condición de contorno de convección, gráficas de Heysler:
Temperatura en r = 0,025 m al cabo de 2 minutos:
r
R
= 0,5 ⇒
T0,025 m - 20
50,8 - 20
= 0,82 ⇒ T0,025 m = 20 + 0,82 (50,8 - 20) = 45,25ºC

Temperatura en r = 0,05 m al cabo de 2 minutos:
r
R
= 1 ⇒
T0,025 m - 20
50,8 - 20
= 0,38 ⇒ T0,025 m = 20 + 0,38 (50,8 - 20) = 31,7ºC
Temperatura en el centro de la esfera al cabo de 5 minutos:
Bi-1
= 0,41
Fo =
α t
R 2
=
1,106.10-5
(m2
/seg) 300 seg
0,052
m2
= 1,3272
⎧
⎨
⎪
⎩
⎪
⇒
T(0,t) - TF
T0 - TF
= 0,0035
T(0,t) - 20
300 - 20
= 0,0035
⎧
⎨
⎪
⎩
⎪
⇒ T( 0,t) = 21ºC
c) Calor evacuado al medio exterior en todo el proceso.
Q(t) = Q0 (1 - e-Bi Fo
) =
Q 0 = ρ V cpF (T0 - TF ) = ρ
4 π R 3
3
cpF ( T0 - TF ) =
= 7850 (kg/m3 )
4 π x 0,053
3
m3x 460 (J/kg ºC) x (780 - 20)ºC = 1,437.106 J = 343,27 Kcal
=
= 343,27 Kcal (1 - 0,3684) = 216,8 Kcal
Calor evacuado al agua al cabo de 5 minutos:
En la esfera quedan almacenados: Q0* =343,27 - 216,8 = 126,45 Kcal
Bi2
Fo = 2,435 x 1,3272 = 7,87
Bi = 2,435
⎫
⎬
⎭
Gröber
⎯ →⎯⎯⎯
Q
Q0
*
= 1
luego al agua iría prácticamente toda la energía almacenada.
*****************************************************************************************
III.34.- Una esfera de hierro maciza de 10 cm de diámetro se encuentra inicialmente a una temperatura de 780ºC
y se enfría en dos etapas. En la primera etapa se utiliza un chorro de aire a 5 m/seg y 20ºC, hasta que en el centro
adquiere una temperatura de 300ºC, y a continuación, en la segunda etapa, se introduce en un baño de mercurio
a 20ºC a una velocidad de 0,1 m/seg.
Determinar
a) El tiempo que debe transcurrir para el enfriamiento en el aire
b) La temperatura en el centro de la esfera, y en la periferia al cabo de los 10 primeros segundos y al cabo de
1 minuto una vez introducida la esfera en el mercurio
Datos: Propiedades térmicas de la esfera: k = 40 W/mºC ; cp = 460 J/kgºC ; ρ = 7850 kg/m3
_________________________________________________________________________________________
RESOLUCIÓN
a) Enfriamiento en aire en convección forzada, u = 5 m/seg
Propiedades del aire a T = Tinicial + Tfinal
2
=
Tinicial =
Ti + TF
2
= 780 + 20
2
= 400
Tfinal =
Tf + TF
2
= 300 + 20
2
= 160
= 400 + 160
2
= 280ºC
o también:
Temperatura media de la pared:
780 + 300
2
= 540ºC
Temperatura del fluido (aire) = 20ºC
⎫
⎬
⎪
⎭⎪
⇒ ˆTaire =
540 + 20
2
= 280ºC
a esta temperatura se tiene:
ρ = 0,6423
kg
m3
; cp = 1,04 kJ
kg C
; ν = 44,34.10-6 m2
seg
; k = 0,0436 W
m C
; α = 0,6532.10 -4 m2
seg
; Pr = 0,68
Cálculo del coeficiente de convección para el flujo del aire sobre la esfera:
Nu = 0,37 Re0 ,6 =
17 Re 70000
Re =
d u F
ν
=
0,1 x 5
44,34.10 -6
= 11276
= 0,37 x 112760,6 = 99,88
hC aire = Nu k
d
=
99,88 x 0,0436
0,1
= 43,55 W
m2.ºC

Con este dato se puede determinar la condición de contorno correspondiente; como se trata del enfriamiento de un
metal en aire, parece lógico que la condición de contorno sea de resistencia térmica interna despreciable, para lo cual
Bi 0,1.
Bi =
hcF R
k
=
43,55 x 0,05
40
= 0,054 0,1 (Condición de contorno con R.T.I.D), por lo que:
Bi =
hcF (R/3)
k
=
43,55 (0,1/6)
40
= 0,018
Fo =
α t
(R/3)2 = α =
k
ρ cp
=
40
7850 x 460
1,106.10-5 m2
seg
=
1,106.10-5 t
(0,1/3)2 = 0,03987 t
T(x,t) - TF
T0 - TF
= e-Bi Fo ⇒
300 - 20
780 - 20
= e-(0,018 x 0,03987 t )= e-7,235.10-4 t ⇒ t = 1391 seg
b) Enfriamiento en mercurio en convección forzada, u = 0,1 m/seg.- Por tratarse de mercurio, en este enfriamiento
la condición de contorno es isotérmica. La ecuación particular a utilizar para calcular Nu dice que la temperatura para
determinar las propiedades del fluido será la media entre la pared y el fluido y en condiciones isotérmicas esa tempe-
ratura es TF
Propiedades del mercurio a 20ºC
ρ = 13579
kg
m3
; cp = 139,4
kJ
kg C
; ν = 0,114.10 -6 m2
seg
; k = 8,69
W
m C
; α = 46,06.10 -7 m2
seg
; Pr = 0,0249
Nu = 2 + 0,386 Re Pr =
30000 Re 150000
Re =
d uF
ν
=
0,1 x 0,1
0,114.10-6
= 87720 = 2 + 0,386 87720 x 0,0249 = 20,04
h C mercurio = Nu k
d
=
20,04 x 8,69
0,1
= 1741,5 W
m2 ºC
Temperatura en el centro de la esfera al cabo de los 10 primeros segundos una vez introducida la esfera en el
mercurio
Fo =
α t
R2 = α =
k
ρ cp
=
40(W/mºC)
7850(kg/m3 ) x 460 (J/kgºC)
= 1,106.10-5 m2
seg
=
1,106.10-5 t
0,052 = 0,00442 t
Al ser Fo = 0,0442 0,2 , hay que tomar varios sumandos del desarrollo:
TC - TF
T0 - TF
= -2
n=1
∞
∑ (−1)n e-λn
2 αt = λn =
π n
R ⇒
n= 1 ; λ1= 1 π/0,05 = 62,83
n= 2 ; λ2 = 2 π/0,05 = 125,61
n= 3 ; λ3 = 3 π/0,05 = 188,5
n= 4 ; λ4 = 4 π/0,05 = 251,32
⎧
⎨
⎪
⎩
⎪
=
€
= - 2 (- e- 62,832
x 1,107.10-5
x 10
+ e- 125,612
x 1,107.10-5
x 10
- e- 188,52
x 1,107.10-5
x 10
+ ...) = 0,98237
TC = TF + (T0 - TF) x 0,98237 = 20 + (300 - 20) x 0,98237 = 295,06ºC
Temperatura en el centro al cabo de 1 minuto una vez introducida la esfera en el mercurio
TC - TF
T0 - TF
= - 2 (- e- 62,832 x 1,107.10-5 x 60 + e- 125,612 x 1,107.10-5 x 60 - e- 188,52 x 1,107.10-5 x 60 + ...) = 0,1452
TC = TF + (T0 - TF) x 0,1452 = 20 + (300 - 20) x 0,1452 = 60,67ºC
La temperatura de la periferia de la esfera, en cualquier instante es la temperatura del fluido = 20ºC.
*****************************************************************************************
III.35.- Una esfera maciza de acero de diámetro 10 cm se encuentra a una temperatura de 890ºC y se introduce en
una corriente de aire a 10 m/seg y 10ºC. Las propiedades del acero son: k = 25 Kcal/m.h.ºC, ρ = 7800 kg/m3 ; cp =
0,11 Kcal/kgºC. Con estos datos determinar:
a) El coeficiente de convección
b) La temperatura en el centro de la esfera al cabo de 1 segundo y al cabo de 15 segundos.
c) El tiempo que tardará en coger en el centro 300ºC
d) La temperatura de un punto situado a una distancia radial de 2,5 cm al cabo de 15 segundos y el calor eva-
cuado en ese tiempo.
_________________________________________________________________________________________
RESOLUCIÓN

Propiedades del aire a: T =
TpF + TF
2
= 890 + 10
2
+ 273 = 723ºK
ρ = 0,4882
kg
m3
; cpF = 1,08
kJ
kgºC
; ν = 6,964.10 -5 m2
seg
; k = 0,05358
W
mºC
; Pr = 0,6854
Nu = 0,37 Re0,6 =
17 Re 70000
Re =
d uF
ν
=
0,1 m x 10 (m/seg)
6,964.10-5 (m2 /seg)
= 14357 = 0,37x143570,6 = 115,47
a) Coeficiente de convección
h cF =
115,47 x 0,05348
0,1
= 61,86
W
m2
ºC
b.1) Temperatura en el centro de la esfera al cabo de 1 segundo.- Al ser un enfriamiento en aire, lo más probable es
que sea un caso con resistencia térmica interna despreciable, por lo que si Bi 0,1 se confirma.
Bi =
h cF (R/3)
k
=
61,86 (W/m2
ºC) (0,1/6) m
25 x 1,163 (W/mºC)
= 0,035 0,1 ⇒ (Se confirma la condición de contorno)
Fo1seg = α t
( R/3)2
= α = k
ρ cpF
=
25(Kcal/h m2
ºK)
7800 (kg/m3
) x 0,11 (Kcal/kgºK)
= 0,02914 m2
hora
= 8,0937.10 -6 m2
seg
= 0,0291
T- TF
T0 - TF
= e-Bi Fo
⇒ T1 seg = 10 + (890 - 10) e- 0 ,035 x 0,0291
= 889 C
b.2) La temperatura en el centro de la esfera al cabo de 15 segundos.
T15 seg = Fo = 0,0291 x 15 = 0,437 = 10 + (890 - 10 )e- 0,035 x 0,437
= 876,6ºC
c) Tiempo que tardará en coger en el centro 300ºC (y por la condición de contorno, en cualquier punto de la esfera)
Fo = α t
(R/3)2
=
8,0937.10
-6
t
(0,05/3) 2
= 0,0291 t
T - TF
T0 - TF
= e- Bi Fo ;
300 - 10
890 - 10
= e- 0 ,035 x 0,0291 t ⇒ t = 1088 seg
d) La temperatura de un punto situado a una distancia radial de 2,5 cm al cabo de 15 segundos
€
T15 seg; 2,5 cm = 876,6C
El calor evacuado en ese tiempo es:
Q(t) = h cF A (T0 - TF ) t
1 - e- Bi Fo
Bi Fo
= 61,86 x 4 π x 0,052
x (890 - 10) x 15
1 - e-0 ,035 x 0 ,4371
0,035 x 0,4371
= 25,46 kW
*****************************************************************************************
III.36.- Una esfera de acero de diámetro 2 cm se encuentra a una temperatura de 800ºC, y se introduce con una
velocidad de 0,2 m/seg, en un recipiente que contiene mercurio a 20ºC. Las propiedades del acero son: k = 25
Kcal/h.m.ºC; ρ = 7800 kg/m3 ; cp = 0,11 Kcal/kgºC
Determinar
a) La temperatura a la que hay que tomar las propiedades del mercurio, y calcular el coeficiente de convec-
ción correspondiente.
b) La ecuación de la distribución de temperaturas en el centro de la esfera, en función del tiempo.
c) Haciendo uso del apartado (b) la temperatura en el centro de la esfera al cabo de 1, 2, 10 y 25 segundos.
d) La temperatura de un punto situado a una distancia radial de 0,5 cm al cabo de 2 segundos
e) El calor evacuado al exterior al cabo de 10 segundos
______________________________________________________________________________________
RESOLUCIÓN
a) Temperatura a la que hay que tomar las propiedades del mercurio.- Al ser el medio de enfriamiento un metal lí-
quido, con toda seguridad se va a tener una condición de contorno ISOTÉRMICA, por lo que la temperatura a la que
hay que considerar las propiedades del fluido son a TF, en este caso a 20ºC, por lo que:
T =
TpF + TF
2
= 20 + 20
2
= 20ºC

es decir, prácticamente en el mismo instante de introducir la esfera en el metal líquido, su superficie adquiere la tem-
peratura del medio de enfriamiento (metal líquido).
Datos del mercurio a 20ºC: ν = 0,114.10-6 m2
seg
; k = 8,69
W
mºC
; Pr = 0,0249
Re =
u d
ν
=
0,2 x 0,02
0,114.10-6
= 35088
Coeficiente de convección: Nu = 2 + 0,386 (Re Pr)0,5 = 2 + 0,386 (35088 x 0,0249 )0,5 = 13,41
hcF = Nu k
d
=
13,4 x 8,69
0,02
= 5826 W
m2.ºK
=
5,826
1,163
Kcal
h.m2.ºC
= 5010 Kcal
h.m2.ºC
b) Ecuación de la distribución de temperaturas en el centro de la esfera, en función del tiempo.
Fo = α t
R 2
= α = k
ρ cp
=
25 (Kcal/hm2
ºC)
7800 ( kg/m3
) 0,11 (Kcal/kgºC)
= 0,0291 m2
hora
=
0,0291 ( m2
/h ) t
0,012
m2
= 0,08 t (t en seg)
TC - TF
T0 - TF
= - 2
n=1
∞
∑ (-1)n
e- λ n
2
α t
, con:
λ n = π n
R
= π n
0,01
= 314,16 n
α =
k
ρ cp
=
25 (Kcal/h mºC)
7800 (kg/m3
) 0,11 (Kcal/kgºC)
= 0,029
m2
hora
⎧
⎨
⎪
⎩
⎪
TC - 20
800 - 20
= -2
n=1
∞
∑ (- 1)n
e- (314 ,16 n )2
0,029 t
= -2
n=1
∞
∑ (- 1)n
e- 2862 n 2
t
TC = 20 - 2 x 780
n=1
∞
∑ (- 1) n
e- (314,16 n )2
0,029 t
= 20 - 1560
n=1
∞
∑ (- 1)n
e- 2862 n 2
t
c) Haciendo uso del apartado (b) la temperatura en el centro de la esfera al cabo de 1 segundo.- (Para tiempos de 1
y 2 seg harán falta varios sumandos por ser Fo 0,2)
TC (1 seg) = 20 - 2 x 780
n=1
∞
∑ (- 1)n e- 2862 n2t = t = 1 seg =
1
3600
h = 20 - 1560
n=1
∞
∑ (- 1)n e
- 2862 n2 1
3600 =
= 20 - 1560
n=1
∞
∑ (- 1)n e- 0,795 n2
= 20 - 1560 (- e- 0,795 + e- 0795 x 4 - e- 0795 x 9 +...) =
= 20 - 1560 (-0,451 + 0,04158 - 0,00078+...) = 20 + 638 = 658ºC
En la misma forma se obtendría para 2, 10 y 25 segundos,
TC(2 seg) = 336ºC
TC(20 seg) = 20,55ºC
TC(25 seg) = 20,000003ºC
⎧
⎨
⎪
⎩
⎪
d) Temperatura en un punto situado a una distancia radial de 0,5 cm, al cabo de 2 segundos
T - TF
T0 - TF
= -
2 R
π
n=1
∞
∑
(-1)n exp (-λn
2 α t)
n r
sen
π n r
R
T - 20
800 - 20
= -
2 x 0,01
π
n=1
∞
∑
(-1)n exp (-2862 n2 t)
0,005 n
sen
π n 0,005
R
= t =
1
3600
h =
= - 6,366.10-3
n=1
∞
∑
(- 1)n exp(- 1,59n2 )
0,005 n sen (
π n
2 ) = - 6,366.10-3{
exp(- 1,59n2 )
0,005 +
exp(- 1,59x 32 )
0,005x3 -
exp(- 1,59x 52 )
0,005x5 + ... } =
= - 6,366.10-3
{- 40,78 + 0,0000406 - ...) = 0,2596 ⇒ Tr=0 ,5cm = 20 + (780 x 0,2596) = 222,5ºC
e) Calor evacuado al exterior al cabo de 10 segundos
Q
Q0
=
6
π2
n=1
∞
∑
1 - exp (-λn
2 α t)
n2 =
6
π2
n=1
∞
∑
1 - exp (-
2862
360
n2 )
n2 =
=
6
π2 {
1 - exp(-
2862
360
12 )
12 +
1 - exp(-
2862
360
22 )
22 +
1 - exp(-
2862
360
32 )
32 + ...} =
=
6
π2 (0,9996 + 0,25 + 0,111 + 0,0625 + 0,04 + ...) = 0,90642

Q = 0,90642 Q0 =
Q0 = m cp ΔT = V ρ cp ΔT =
4
3
π R3ρ cp ΔT =
=
4
3
π 0,013 m3x 7800
kg
m3 x 0,11
Kcal
kg C
x780 C = 2,8 Kcal
= 0,9 x 2,8 = 2,52 Kcal
*****************************************************************************************
III.37.- Se desea hacer un tratamiento térmico de unas bolas de diámetro d = 3 cm construidas con acero al Cr-Ni
(20/10), α = 0,415 x10-5 m2/seg, k = 15,1 W/m.ºK, y para ello se calientan a una temperatura de 975ºC. A conti-
nuación se enfrían en un baño de mercurio que se encuentra a la temperatura de 20ºC, sometiéndolas a un movi-
miento de agitación dotado de una velocidad media de 0,5 m/seg.
Determinar: a) El coeficiente de transferencia de calor
b) El tiempo necesario para que las bolas alcancen la temperatura de 200ºC en el centro
c) El calor disipado en ese tiempo
_________________________________________________________________________________________
RESOLUCIÓN
a) Coeficiente de transferencia de calor.- Al ser el medio de enfriamiento un metal líquido, casi con toda seguridad
se va a tener una condición de contorno ISOTÉRMICA
Datos de Mercurio a 20ºC:
€
ρ = 13579
kg
m3 ; cp = 139,4
kJ
kg°C
; ν = 0,114.10-6 m2
seg
; k = 8,69
W
m°C
; α = 46,06.10-7 m2
seg
; Pr = 0,0249
Re =
u F d
ν
=
0,5 (m/seg) x 0,03 m
0,114.10 -6
(m2
/seg)
= 131.578 ; Fo =
α t
R2
=
46,06.10
-7
(m
2
/seg) t
0,0152
m2
= 0,02 t
Nu = 2 + 0,386 (Re Pr)0,5 ; campo de aplicación : (3 x 104 Re 1,5 x 105) (válida)
Nu = 2 + 0,386 (Re Pr)0,5 = 2 + 0,386 (131579 x 0,0249 )0,5 = 24,1
hcF = Nu k
d
=
24,1 x 8,69
0,03
= 6981 W
m2.ºK
b) Tiempo necesario para que las bolas alcancen la temperatura de 200ºC en el centro
TC - TF
T0 - TF
= - 2
n=1
∞
∑ (-1)n
e- λn
2
α t
= λn =
π n
R
; α = 0,415.10 -5 m2
hora
=
200 - 20
975 - 20
= 0,1884
0,094 = -
n=1
∞
∑ (-1)n
e
-λn
2
x 0,415.10-5 m 2
hora
t
= ... =
λ1 = π/R = π/0,015 = 209,44
λ2 = 2π/R = 2π/0,015 = 418,88
λ3 = 3π/R = 3π/0,015 = 628,34
=
= {(-1) e- 209,442
x 0,415.10-5
t
+ (1) e- 418,882
x 0,415.10-5
t
+ (-1) e- 628,32
x 0,415.10-5
t
+ ...} =
= {(-1) e- 209,442 x 0,415.10-5 t + (1) e- 418,882 x 0,415.10-5 t + (-1) e- 628,32 x 0,415.10-5 t + ...} = e-0,182 t - e- 4 x 0,182 t + e- 9 x 0,182 t + ...
observándose que el desarrollo en serie converge muy rápidamente.
Como: Fo = 0,02 x12,37 = 0,253 0,2, resolviendo con sólo el primer sumando se obtiene:
0,1052 = e-0 ,182 t
t = 12,37 seg
⎧
⎨
⎩
c) Calor disipado en ese tiempo
Q
Q0
= 6
π2
n=1
∑ 1 - e- λn
2
α t
n 2
= 6
π2
( 1 - e- λ1
2
α t
12
+ 1 - e- λ2
2
α t
2 2
+ 1 - e- λ3
2
α t
32
+ ...) =
= 6
π2
{ 1
12
+ 1
22
+ 1
32
+ ... - ( e- λ1
2
α t
12
+ e- λ2
2
α t
22
+ e- λ 3
2
α t
32
+ ...)} = 6
π2
( π
2
-
n=1...
∑ e- λn
2
α t
n 2
) =
=
6
π2
{
π
2
- (
e-0 ,182 x 12,37
12
+
e- 4 x 0 ,182 x 12,37
2 2
+
e-9 x 0 ,182 x 12,37
32
+ ...)} =
6
π2
(
π
2
- 0,10523 - 0,000031) = 0,89
Q = 0,89 Q0 =
Q0 = m cp Δt = V ρ cp Δt = ρ cp =
kesf
αesf
=
4 π R 3
3
k ΔT
α
=
=
4 π 0,0153
3
m3
15,1 (W/mºK) x 955ºK
0,415.10 -5
(m2
/seg)
= 49,124 kJ
⎧
⎨
⎪
⎩
⎪
⎫
⎬
⎪
⎭
⎪
= 43,72 kJ

*****************************************************************************************
III.38.- Una esfera de hierro macizo de diámetro 4 cm, se encuentra a la temperatura de 390ºC y tiene las siguien-
tes constantes físicas: k = 65 W/m.ºC ; cp = 460 Joules/kgºC ; ρ =7800 kg/m3
Se la somete a un enfriamiento en un medio exterior que se encuentra a 20ºC y se desea saber:
a) En el supuesto de que el medio de enfriamiento sea aire en reposo, hallar al cabo de 1 minuto la temperatu-
ra en el centro y en la periferia. ¿Qué tiempo tardaría en alcanzar en el centro la temperatura de 100ºC?
b) Si se supone que el medio de enfriamiento es una corriente de agua que circula a 1 m/seg, determinar la
temperatura en el centro y en la periferia al cabo de 10 segundos y de 1 minuto
_________________________________________________________________________________________
RESOLUCIÓN
a) Temperatura en el centro y en la periferia, para enfriamiento en aire en reposo, al cabo de 1 minuto
Coeficiente de convección para el aire
€
Las propiedades del aire se calculan a: T =
TpF + TF
2
=
390 + 20
2
+ 273 = 478ºK
ρ = 0,744
kg
m3
; cpF = 1,025 kJ
kgºC
; k = 0,03872 W
mºC
; ν = 34,8.10-6 m2
seg
; Pr = 0,68
g β
ν2
=
9,8 (1/293ºK)
(34,8.10-6
)2
= 2,76.10 7
Gr =
g β
ν2
ΔT d3
=
ΔT = 390 - 20 = 370 ºC
d = 0,04 m
= 2,76.10 7
x 370 x 0,043
= 654000
Ra = Gr Pr= 654000x 0,68= 444725
Nu d = 2 + 0,45 Rad
0,25
= 2 + 0,45 x 4447250,25
= 13,62 ⇒ hcF =
Nud k
d
=
13,62 x 0,03872
0,04
= 13,185
W
m2
ºC
Al ser aire el medio de enfriamiento es de prever una condición de contorno con resistencia térmica interna despre-
ciable, por lo que:
Bi =
h cF (R/3)
k hierro
=
13,185 x (0,02/3)
65
= 0,00135
Para t = 1 min ⇒ Fo1 min =
α t
(R/3)2
= α =
k
ρ cpF
=
65
7800 x 460
= 1,81.10 -5
(m2
/seg) =
1,81.10 -5
(m2
/seg) 60 seg
(0,02/3)2
= 24,43
T( x,t) - TF
T0 - TF
= e- Bi Fo
;
T(x,t) - 20
390 - 20
= e- 1,219.10-3
x 24,43
= 0,9706 ⇒ T(x,t) = 379,14ºC
Tiempo que se emplea en alcanzar en el centro la temperatura de 100ºC
100 - 20
390 - 20
= 0,2162 = e
-
1,81.10-5
t
( 0,02/3)2
x 0,00135
= e- 5,507.10-4
t
⇒ t = 2781 seg = 46,35 min
b) Si se supone que el medio de enfriamiento es una corriente de agua que circula a 1 m/seg, determinar la tempe-
ratura en el centro y en la periferia al cabo de 1 minuto
Flujo en convección forzada.- Propiedades del agua a 20ºC, para aplicar la ecuación de Whitaker:
ν = 1,006.10-6 m2
seg
; k = 0,579
W
mºK
; Pr = 7
Re =
u F d
ν
=
1 m/seg x 0,04 m
1,006.10 -6
m2
/seg
= 39761
Nu = 2 + (0,4 Re + 0,06 Re2/3
) Pr 0,4 ηF
ηpF
4 , válida en:
3,5 Re 80000
0,7 Pr 380
⎧
⎨
⎩
ηpF se calcula a la temperatura media de la pared: 390ºC = 205ºC ⇒ ηpF =132.10-6 Nseg/m2
ηF se calcula a la temperatura media del fluido 20ºC ⇒ 1004.10-6 Nseg/m2

Nu = 2 + (0,4 39761 + 0,06 x 397612 / 3
) 70 ,4 1004
132
4 = 542,9
h cF =
Nu k
d
=
542,9 x 0,597
0,04
= 8100
W
m2
ºC
Temperatura en el centro y en la periferia de la esfera al cabo de 10 segundos
h cF
k
=
8100
65
= 124,6 ; Bi =
hcF R
k
=
8100 x 0,02
65
= 2,5
Fo =
α t
R 2
=
1,81.10-5 x 10
0,022
= 0,4525 1
⎫
⎬
⎪
⎪
⎭
⎪
⎪
⇒ Condición de contorno isotérmica
Temperatura en la periferia: 20ºC (Por la condición de contorno isotérmica)
TC - TF
T0 - TF
= - 2
n=1
∞
∑ (-1)n
e-λn
2
α t
=
Basta con el primer término del desarrollo
λn =
π n
R
=
π n
0,02
⇒ λ1 = 157,08 ; λ2 = 314,16
=
= - 2 (- e-157,082
x 1,81.10-5
x 10
+ ...) = 2 (e- 4,466
- e- 17,86
) = 0,022986 =
TC - 20
390 - 20
⇒ TC = 28,5ºC
Temperatura en el centro y en la periferia de la esfera al cabo de 60 segundos
Fo = α t
R 2
=
1,81.10-5 x 60
0,022
= 2,715
El proceso pierde la condición isotérmica por lo que con este valor y Bi-1 = 0,4 se obtiene:
TC - TF
T0 - TF
= 2
n=1
∞
∑ e-λ n
2
α t sen (λn R) - (λ n R) cos (λn R)
λn R - sen(λ n R)cos (λn R)}
⇒
TC - 20
390 - 20
= 3,5.10-8
⇒ TC = 20 ºC
*****************************************************************************************
III.39.- Un ladrillo macizo de (5 x 10 x 12) cm se extrae de un horno que se encuentra a 400ºC y se coloca en una
mesa en posición vertical (apoyado sobre la base menor). El coeficiente de película medio para todo el ladrillo
puede suponerse del orden de 25 Kcal/h.m2.ºC. La temperatura del ambiente es de 20ºC.
Datos del material: α = 0,002 m2/hora ; k = 1 Kcal/h.m.ºC
Determinar, al cabo de 1 hora:
a) La temperatura en el centro de la cara apoyada en la mesa y en el centro del ladrillo
b) La temperatura en el centro de la cara del ladrillo opuesta a la apoyada en la mesa
c) La temperatura en el centro de la cara del ladrillo de dimensiones (10 x 12) cm
_________________________________________________________________________________________
RESOLUCIÓN
a) Temperatura en el centro de la cara apoyada en la mesa.- Veamos la condición de contorno:
Disipa calor en las tres direcciones; al duplicar la longitud (por el plano apoyado sobre la base menor, que se convier-
te en un plano adiabático y, por lo tanto, de simetría), parece que podría hacerlo sólo en una, hecho que habrá que
comprobar.
hC/k = 25/1 = 25 ⇒ c.c. convección
Bi =
hC (V/A)
k
= 25 x (0,012/0,082)/1 = 0,366 0,1 ⇒ c.c. convección
⎧
⎨
⎩
Dirección Cx:
Bi =
hcF x
k
=
25(0,05/2)
1
= 0,625 ; Bi-1 = 1,6
Fo =
α t
x2
=
0,002 x 1
(0,05/2)2
= 3,2
⎧
⎨
⎪
⎩
⎪
⇒ Px (0) = 0,23
Dirección Cy:
Bi =
hcF y
k
=
25(0,1/2)
1
= 1,25 ; Bi-1 = 0,8
Fo =
α t
y2
=
0,002 x 1
(0,1/2)2
= 0,8
⎧
⎨
⎪
⎩
⎪
⇒ Py (0) = 0,53

Dirección Cz:
Bi =
hcF z
k
=
25(0,24/2)
1
= 3 ; Bi-1 = 0,333
Fo =
α t
z2
=
0,002 x 1
(0,24/2)2
= 0,1389
⎧
⎨
⎪
⎩
⎪
⇒ Pz (0) = 0,91
TC - TF
T0 - TF
= 0,23 x 0,53 x 0,91 = 0,111 ⇒ TC = 20 + 0,111 (400 - 20) = 62,18ºC
Se observa que el elemento corrector en la dirección Cz es prácticamente la unidad
Temperatura en el centro del ladrillo TA
En la dirección Cx no hay corrección por cuanto el punto A corresponde al plano (y,z).- Parámetro
de posición 0 ; Heysler ⇒ 1
En la dirección Cy no hay corrección por cuanto el punto A corresponde al plano (x,z).- Parámetro
de posición 0 ; Heysler ⇒ 1
En la dirección Cz sí hay corrección: Bi = 3 ; Fo = 0,1389
Parámetro de posición:
z
L
= 0,5 ; Bi-1
= 0,33 ⇒ Heysler 0,83
⎧
⎨
⎩
TA - TF
TC - TF
= 1 x 1 x 0,83 = 0,83 ; TA = 20 + 0,83 (62,18 - 20) = 55ºC
b) La temperatura en el centro de la cara del ladrillo opuesta a la apoyada en la mesa, punto B
En la dirección Cx no hay corrección por cuanto el punto B corresponde al plano (y,z).- Paráme-
tro de posición 0 ; Heysler ⇒ 1
En la dirección Cy no hay corrección por cuanto el punto B corresponde al plano (x,z).- Paráme-
tro de posición 0; Heysler ⇒ 1
z
L
= 1 ; Bi-1
= 0,33 ⇒ Heysler 0,38
⎧
⎨
⎩
TB - TF
TC - TF
= 1 x 1 x 0,38 = 0,38 ; TB = 20 + 0,38 (62,18 - 20) = 36,03ºC
c) Temperatura en el centro D de la cara del ladrillo de dimensiones (10 x 12) cm
En la dirección Cx sí hay corrección: Bi = 0,625 ; Bi-1
= 1,6
x
L
= 1 ; Bi-1
= 1,6 ⇒ Heysler 0,73
⎧
⎨
⎩
En la dirección Cy no hay corrección por cuanto el punto D corresponde al plano (x,z).- Paráme-
tro de posición 0 ; Heysler ⇒ 1
z
L
= 0,5 ; Bi-1
= 0,33 ⇒ Heysler 0,83
⎧
⎨
⎩
TD - TF
TC - TF
= 0,73 x 1 x 0,83 = 0,6059 ⇒ TD = 20 + 0,6059 (62,18 - 20) = 45,56ºC
*****************************************************************************************
III.40.- Un ladrillo macizo de dimensiones (10 x 12 x 5) cm, se extrae de un horno que se encuentra a 400ºC, y se
coloca encima de una mesa, en un ambiente tal que el coeficiente de película vale hc = 40 Kcal/h.m2.ºC, siendo la
temperatura de dicho ambiente de 20ºC. El ladrillo se ha apoyado en la mesa sobre la cara de (10x12) cm y se su-
pone que la cantidad de calor transmitida a la mesa es despreciable.
Determinar al cabo de 1 hora
a) La temperatura superficial en el centro de la cara de (5 x 10) cm.
b) La temperatura de uno de los vértices en contacto con la mesa
Datos del material: α = 0,002 m2/hora ; k = 0,6 Kcal/h.m.ºC
_________________________________________________________________________________________

RESOLUCIÓN
Como a través de la mesa el flujo de calor es cero, se puede considerar
a la superficie del ladrillo en contacto con la mesa como un plano
adiabático, por lo que se puede suponer a la pieza de espesor doble (2
x 5) cm, es decir, se somete una pieza de espesor doble a la acción del
medio exterior, con lo que el ladrillo ficticio será de (10 x 12 x 10) cm,
disipando calor por todas las caras laterales.
a) Temperatura superficial en el centro de la cara de (5 x 10) cm.
Aplicando el método de superposición: Temperatura en el centro O
Dirección Ox,
Bi =
hC Lx
k
=
40(0,12/2)
6
= 4 ⇒ Bi-1 = 0,25
Fo=
α t
Lx
2
=
0,002 x 1
(0,12/2)2
= 0,555
⎧
⎨
⎪
⎩
⎪
⎫
⎬
⎪
⎭
⎪
⇒ Px (0) =
T(0,t) - TF
T0 - TF
= 0,42
Dirección Oy,
Bi =
hC Ly
k
=
40(0,1/2)
6
= 3,33 ⇒ Bi-1 = 0,3
Fo= α t
Ly
2
=
0,002 x 1
(0,1/2)2
= 0,8
⎧
⎨
⎪
⎩
⎪
⎫
⎬
⎪
⎭
⎪
⇒ Py (0) =
T(0,t) - TF
T0 - TF
= 0,35
Dirección Oz,
Bi =
hC Lz
k
=
40(0,1/2)
6
= 3,33 ⇒ Bi-1 = 0,3
Fo=
α t
Lz
2
=
0,002 x 1
(0,1/2)2
= 0,8
⎧
⎨
⎪
⎩
⎪
⎫
⎬
⎪
⎭
⎪
⇒ Pz (0) =
T(0,t) - TF
T0 - TF
= 0,35
luego:
T(0,t) - 20
400 - 20
= 0,42 x 0,35 x 0,35 = 0,05145 ⇒ T(0,t) = 39,55ºC
TEMPERATURA EN EL PUNTO A
Dirección Ox,
Bi = 4 ; Bi-1 = 0,25
Parámetro de posición,
x
Lx
=
Distancia al centro
Semiespesor
=
0,06
0,06
= 1
⎧
⎨
⎩
⎫
⎬
⎭
⇒
T(x,t) - TF
T(0,t) - TF
= 0,3
Dirección Oy,
Bi = 3,33 ⇒ Bi-1 = 0,3
Parámetro de posición,
y
Ly
=
0,025
0,05
= 0,5
⎧
⎨
⎩
⎫
⎬
⎭
⇒
T(y,t) - TF
T(0,t) - TF
= 0,82
Dirección (Oz).- Respecto a esta dirección no hay corrección por cuanto el punto A pertenece al plano (x,y).
(Parámetro de posición 0)
luego:
TA - TF
T(0,t) - TF
=
TA - 20
39,55 - 20
= 0,3 x 0,82 x 1 = 0,246 ⇒ TA = 24,8ºC
De otra forma:
TA - TF
T0 - TF
=
TA - 20
400 - 20
= (0,42 x 0,3) (0,35 x 0,82) (0,35 x 1) = 0,0127 ⇒ TA = 24,8ºC
b) Temperatura de cada uno de los vértices en contacto con la mesa
Dirección Ox:
Bi = 4 ⇒ Bi-1 = 0,25
x
Lx
=
Distancia al centro
Semiespesor
=
0,06
0,06
= 1
⎧
⎨
⎩
⎫
⎬
⎭
⇒
T(x,t) - TF
T(0,t) - TF
= 0,3
Dirección (Oy).- Respecto a esta dirección no hay corrección por cuanto el punto M pertenece al plano (x,z).
(Parámetro de posición ξ/L = 0)
Dirección Oz:
Bi = 3,33 ⇒ Bi-1 = 0,3
z
Lz
=
0,05
0,05
= 1
⎧
⎨
⎩
⎫
⎬
⎭
⇒
T(z,t) - TF
T(0,t) - TF
= 0,35

luego:
TM - TF
T(0,t) - TF
=
TM - 20
39,55 - 20
= 0,3 x 1 x 0,35 = 0,105 ⇒ TM = 22,05ºC
*****************************************************************************************
III.41.- Se tiene una chapa de palastro de 40 mm de espesor, en contacto con un muro de material relativamente
aislante de espesor 30 cm. Propiedades físicas del palastro: k = 43 W/m°C ; α = 1,2.10-5 m2/seg
Propiedades físicas del material del muro: k* = 3 W/m°C ; α = 0,75.10-5 m2/seg
La temperatura inicial del conjunto es de 400°C. Se refrigera la chapa de palastro en forma instantánea de modo
que su temperatura superficial desciende a 40°C
Determinar analítica y gráficamente la temperatura que se tiene en la superficie de contacto y a una distancia de
9,5 cm de dicha zona (en el muro), al cabo de 2 minutos 13 segundos.
Representar el proceso en papel milimetrado
_________________________________________________________________________________________
RESOLUCIÓN
El intervalo que tomamos en el palastro es: Δx1 = 2 cm.- Aplicando Binder-Schmidt se obtiene el tiempo correspon-
diente a cada intervalo:
Fo = 0,5 =
α Δt
Δx1
2
⇒ Δt =
Δx1
2
2 α
=
0,022
2 x 1,2.10 -5
= 16,66 seg
por lo que el intervalo Δx2 del muro de material aislante es: Δx2 = 2 α*Δt = 2 x 0,75.10-5
x 16,66 = 0,0158 m
Para determinar la temperatura en la superficie de contacto se hace una rehubicación de los nodos inmediatos a la
misma, en el palastro o en el material del muro, de forma que el intervalo para estos nodos es:
T0
t + Δt
- T1
t
Δx1
=
T2
t
- T0
t + Δt
k1
k2
Δx2
;
k1
k2
Δx2 = 43
3
x 15,8 = 226,58 mm (En el muro)
Como: 2 min 13 seg = 133 seg ⇒ 133
16,66
= 7,98 = 8 incrementos de tiempo
*****************************************************************************************
III.42.- En un proceso de fabricación de hielo el agua fluye a 5ºC sobre una placa a -10ºC enfriada por un refri-
gerante. Cuando la capa de hielo alcanza los 2 mm de espesor, se raspa para quitarla de la placa, repitiéndose el
proceso. Determinar el tiempo que transcurre entre dos raspaduras.
Datos: Coeficiente global de transmisión de calor del refrigerante a la superficie de contacto entre la placa y el
hielo es de 900 W/m2ºK
Calor latente de fusión del hielo: 335 kJ/kg; khielo a 0ºC = 2,22 W/mºK; ρ hielo = 910 kg/m3
_________________________________________________________________________________________
RESOLUCIÓN.
Se puede considerar solidificación unidimensional; las pérdidas de calor sensible del agua al pasar de -5ºC a 0ºC son
pequeñas comparadas con el calor latente del cambio de estado.

Haciendo un balance de energía en la superficie de separación del
agua y el hielo, el calor cedido por el hielo por conducción más el ca-
lor latente del cambio de estado, tiene que ser igual al calor cedido al
agua (calor sensible), por lo que:
k ∂T
∂x
〉hielo + G ragua-hielo = k ∂T
∂x
〉agua
El calor sensible del agua se desprecia ⇒ k
∂T
∂x
〉agua = 0
El calor sensible del hielo es: k
∂T
∂x
〉hielo =
Tsup. separación agua-hielo- Tplaca
s
khielo
+
1
Uplaca-hielo
Cantidad de hielo producido, por unidad de anchura, en dt: G = - ρ ∂s
∂t
k
∂T
∂x
〉hielo + G ragua-hielo = 0 ;
Tsup. separación agua-hielo - Tplaca
s
k hielo
+
1
U placa-hielo
- ρ
∂s
∂t
ragua-hielo = 0
k ∂T
∂x
〉hielo + G ragua-hielo = 0 ;
ρ ragua-hielo
dt = ( s
k hielo
+ 1
Uplaca-hielo
) ds
Como para: t = 0, s = 0, integrando se obtiene:
ρ ragua-hielo
t =
s2
2 k hielo
+
s
Uplaca-hielo
⇒
0 - (- 10)
910 x 335.10 3
t =
0,0022
2 x 2,22
+
0,002
900
⇒ t = 95 seg
*****************************************************************************************
III.43.- Calcular el tiempo necesario para alcanzar la temperatura de -l°C, punto de congelación, en el núcleo de
un paquete paralelepipédico de pescado, de 4 cm de espesor, siendo su temperatura inicial de +20°C. El enfria-
miento se verificará por inmersión del paquete en salmuera a -20°C.
Se considerará la transmisión de calor unidireccional.
Datos salmuera:
cp= 0,671 Kcal/kg°K; ρ= 1260 kg/m3; η = 0,0119 kg/m.seg; k = 0,457 Kcal/h.m°C ; u = 0,1 m/seg
Datos pescado sin congelar: cp* = 0,85 Kcal/kg°C ; ρ* = 1060 kg/m3 ; k* = 0,445 Kcal/h.m°C
_________________________________________________________________________________________
RESOLUCIÓN
Tiempo para pasar de 20ºC a -1ºC en la superficie
Re =
L u F
ν
=
0,04 m x 0,1 (m/seg)
0,0119 (kg/m.seg)/1260 (kg/m3
)
= 423,5 (reg. laminar)
Pr =
η cF
k =
0,0119 (kg/m seg) x 0,671 (Kcal/kgºC)
0,457(Kcal/hmºC) x 3600
seg
hora = 62,9
Para placa plana:
h cF =
0,664 k ReL Pr1/3
L
=
0,664 x 0,457 435,5 x 62,91/ 3
0,04
= 620 Kcal
hm2 ºC
h cF
k*
=
620
0,457
= 1356,67 (C. c. isotérmica)
T( 0,t) - TF
T0 - TF
=
(-1) - (-20)
20 - (-20)
=
19
40
= 0,475 ⇒ Fo = 0,098
Fo = α* t
L2 = α* = k*
ρ* cpF
* =
0,445
1060 x 0,85 = 0,000494 m2
hora =
0,000494 t
0,042 = 0,31 t = 0,098 ⇒
t = 0,3174 horas
ó 19 min 26 seg
⎧
⎨
⎩
Si en el paquete de pescado se desea continuar el enfriamiento hasta alcanzar la temperatura de -15°C en el cen-
tro del paquete, determinar:
a) El tiempo que tardará en congelarse el pescado .- Una vez se ha logrado -1ºC en el núcleo del paquete, se proce-
de de la siguiente forma: En un instante determinado de la congelación supondremos un plano (2) a la distancia δ que
se encuentra a la temperatura -1ºC y que además está en contacto con el pescado no congelado.
La pared exterior está en contacto con la salmuera a -20ºC. El resto del paquete hasta el plano central (3) se encuentra

a la temperatura -1ºC pero no ha congelado.
Supondremos que en un tiempo dt, el calor dQ disipado a través del plano exterior (1) es capaz de congelar un espe-
sor dd del paquete; así se tendrá:
dQ = U A (Tcong - TF ) dt = U = 1
(1/hcF ) + (δ/k)
=
A (Tcong - TF ) dt
(1/h cF ) + (δ/k )
por cuanto se ha considerado que entre los planos (2) y (3), (L - δ), no hay diferencia de temperaturas y, en conse-
cuencia, no existe conducción de calor entre dichos planos.
0,01 0,02 0,03 0,05 0,07 0,1 0,2 0,3 0,4 0,5 0,7 1
Q/Q0
0,4
0,3
0,2
0,1
0,5
0,2
0,4
0,6
0,8
1,0
Q/Q0
Fo
Por otro lado, la cantidad de calor dQ a extraer, y necesaria para producir la congelación del espesor dδ del paquete,
se puede poner en la forma:
dQ = rlíq-sól dG = rlíq-sól x (A dδ) ρ
en la que:
dG = x (A dδ) ρ, es la masa a congelar en el tiempo dt
x, es el % d agua contenida en el pescado, que es la que hay que congelar, ya que el resto es sólido y no experi-
mentará cambio de estado
rlíq-sól es el calor latente de la congelación del agua, del orden de 80 Kcal/kg
dQ = rlíq-sól (A dδ) ρ =
A ( Tcong - TF ) dt
(1/hcF ) + (δ/k)
⇒ dt =
rlíq-sól (A dδ ) ρ
A ( Tcong - TF )
(1/h cF ) + (δ/k )
=
rlíq-sól (A dδ ) ρ
A ( Tcong - TF )
(
1
h cF
+
δ
k
)
t =
rlíq-sól x ρ
Tcong - TF 0
L/ 2
∫ (
1
h cF
+
δ
k
) dδ =
rlíq-sól x ρ
Tcong - TF
(
δ
h cF
+
δ 2
2 k
)0
L/ 2
=
rlíq-sól x ρ
Tcong - TF
(
L/2
h cF
+
(L/2)2
2 k
) =
=
80 (Kcal/kg) x 0,75 x 1000 (kg/m3
)
(-1) - (-20 )
(
0,02 m
620 ( Kcal/hm2
ºC)
+
0,02 2 m2
2 x 1,5 (Kcal/hmºC)
) = 0,65 horas = 39 min 36 seg
b) El tiempo que se tardará en alcanzar los -15°C en el núcleo, después de congelado el pescado
Datos, Pescado congelado: cp = 0,475 Kcal/kg°C ; ρ = 1000 kg/m3 ; k = 1,15 Kcal/h.m°C
Composición: agua=75% , materia sólida=25% . Se considerará que la transmisión de calor es unidireccional
h cF ( pescado−salmuera) = 620 (Kcal/h m2
ºC) ;
h cF
k
=
620
1,15
= 539,13 (C.c. isotérmica)
T( 0,t) - TF
T0 - TF
=
(-15) - (-20)
(-1) - (-20)
= 0,263 ⇒ Fo = 0,168
Fo = α t
L2
= α = k
ρ cpF
=
1,15
1000 x 0,475
= 0,00242 m2
hora
=
0,00242 t
0,042
= 1,513 t = 0,168 ⇒
t = 0,111 h ó
6 min 40 seg
⎧
⎨
⎩
*****************************************************************************************

Problemas transitorio

Más contenido relacionado

La actualidad más candente (20)

Similar a Problemas transitorio (20)

Último (20)

Problemas transitorio