Problemas Conducción Estacionario
Problemas Conducción Estacionario
Problemas Conducción Estacionario
270°K, con un coeficiente de película de 40 W/m2 °K. El lado opuesto de la pared está en contacto con el aire
en calma a 330°K, y coeficiente de película de 10 W/m2 °K. Calcular el calor transmitido por unidad de área
y unidad de tiempo.
_________________________________________________________________________________________
RESOLUCION
q = Tint - Text = Tint - Text
3 R1 + R2 + R3
∑ Ri
hci k hce i=1
R1 = 1 = 1 = 0,025 (°K/W)
A hci A 40 x 1
e
Ti Te
R2 = e = 0,1 = 0,143 (°K/W)
kA 0,7 x 1
Ti To
Tpi Tpe R 3=1/h ceA
1 = 1 = 0,1 (°K/W)
R3 =
R 1=1/hci A R 2=e/kA
hce A 10 x 1
Calor transmitido por unidad de superficie y unidad de tiempo,
q 330 - 270
= = 224 W
A 0,025 + 0,143 + 0,10 m2
*****************************************************************************************
Una pared plana grande, tiene un espesor de 0,35 m; una de sus superficies se mantiene a una tempe-
ratura de 35°C, mientras que la otra superficie está a 115°C. Unicamente se dispone de dos valores de la
conductividad térmica del material de que está hecha la pared; así se sabe que a 0°C, k = 26 W/m°K y a
100°C, k = 32 W/m°K. Determinar el flujo térmico que atraviesa la pared, suponiendo que la conductividad
térmica varía linealmente con la temperatura
_________________________________________________________________________________________
RESOLUCION
Tpe 100
k dT (1 + 0,004 T) dT 2 2
(100 - 80) + 0,004 100 - 80
k=
Tpi
= 80
= 2 = 1,36 Kcal
Tpe - Tpi 100 - 80 100 - 80 h.m.°C
Tpe - Tpi
q=2πkL
r = 2 π x 1,36 ( Kcal ) x 1 m 100 - 80 = 246,56 Kcal
ln e h.m.°C ln 4 hora
ri 2
Temperatura T del tubo en un diámetro correspondiente a, d = 6 cm
2 π x 1,36 ( Kcal ) x 1 m 100 - T = 246,56 Kcal ; T = 91,7 °C
h.m.°C ln 4 hora
3
b) En el supuesto de trabajar con el valor de “k” correspondiente a cada punto del tubo, se puede suponer un
valor de “k” de la forma,
k + kT (1 + 0,004 x 100) + (1 + 0,004 T)
k = Tpe = = 1,2 + 0,002 T
2 2
2 π (1,2 + 0,002 T)(100 - T)
246,56 = ; 0,002 T2 + T - 108,72 = 0 ; T = 91,84°C
ln 4
3
Con esta temperatura de 91,84ºC habría que volver a reiniar los cálculos.
*****************************************************************************************
Un tubo de diámetro de = 0,5 metros, cuya emitancia superficial vale ε=0,9, que transporta vapor de
agua, posee una temperatura superficial de 500°K. El tubo está localizado en una habitación a 27°C, y el
coeficiente de transmisión de calor por convección entre la superficie del tubo y el aire de la habitación se
puede considerar igual a hC = 20 W/m2 °K.
Calcular
a) La conductancia superficial unitaria combinando radiación y convección
b) El calor disipado por unidad de tiempo y por metro de longitud del tubo
_________________________________________________________________________________________
RESOLUCION
a) El tubo se puede considerar como un cuerpo emisor, rodeado por un cuerpo negro que es la habitación; ade-
más hay que tener presente la convección. Por lo tanto, la conductancia global será,
h = hC + hR
hC = 20 W
m 2 °C
QR = σ ε A (T4pared - T4est ) = hR A (Tpared - Text )
5,67 x 10-8 W x 0,9 x 1 m2 (5004 - 3004 )°K4
σ ε A (T4pared - T4ext ) m 2 °K4 W
hR = = = 13,88
A (Tpared - Text ) (500 - 300)°K m 2 °C
h = hC + hR = 20 + 13,88 = 33,88 W
m2 °K
b) Pérdida de calor por unidad de tiempo y por metro de longitud de tubo,
Q = π de L h (Tpared - Text ) = π x 0,5 x 1 x 33,88 (500 - 300) = 10.650 W
*****************************************************************************************
Por una tubería de aluminio circula vapor de agua a 110°C. La tubería tiene un coeficiente de conduc-
tividad térmica k = 185 W/m°K, un diámetro interior di = 10 cm, y un diámetro exterior de = 12 cm. La tube-
ría está situada en una habitación en la que la temperatura ambiental del aire es de 30°C, siendo el coefi-
ciente de transferencia térmica convectiva entre la tubería y el aire hc=15 W/m2 °K.
Determinar la transferencia de calor para los siguientes casos,
a) La tubería no se encuentra aislada
b) La tubería se encuentra aislada y, para ello, se recubre con una capa de aislante de 5 cm de espesor,
k1 = 0,20 W/m°K. Se admitirá que es despreciable la resistencia convectiva del vapor.
_________________________________________________________________________________________
RESOLUCION
a) Como el h Cvapor es muy elevado, su resistencia convectiva será muy pequeña (despreciable), y podemos con-
siderar que la temperatura interior de la pared coincide con la temperatura del vapor.
Se observa que la presencia del aislamiento reduce la pérdida de calor en un 70%. En ambos casos se podía
haber despreciado la resistencia térmica de la tubería de Al sin perder mucha exactitud en el cálculo de la trans-
ferencia de calor por unidad de tiempo.
*****************************************************************************************
En la parte exterior de una caldera existe una temperatura de 80ºC. Para evitar accidentes se cons-
truye un muro que aisla la caldera del medio exterior, que se encuentra a 25ºC.
Las dimensiones de este muro son, Longitud 18 m; Altura = 8,50 m; Espesor = 0,40 m
Los coeficientes de transmisión de calor son, hCaire-pared interior = 8 Kcal/m2 .h.ºC
hCaire-pared exterior = 20 Kcal/m2 .h.ºC
La conductividad térmica del muro, k = 0,70 Kcal/m.h.ºC
Régimen estacionario.II.-12
Determinar
a) El número de calorías perdidas al exterior a través del muro
b) Se recubre exteriormente la pared con un material aislante, tanto por la parte que mira a la caldera,
como la que da al medio exterior; su espesor es de 1 cm., y su conductividad térmica k*= 0,06
Kcal/m.h.ºC.¿Cuál será ahora la cantidad de calor cedida al exterior?
_________________________________________________________________________________________
RESOLUCION
a) Número de calorías perdidas al exterior a través del muro
Q= S Text. caldera - Tmedio exterior = (18 x 8,5) m2 80 - 25 = 11.273,7 Kcal
1 +∑ e + 1 1 + 0,4
+ 1
hora
hC(aire-pared interior) k hC(aire-pared exterior) 8 0,70 20
b) Calor cedido al exterior si se recubre exteriormente la pared con un material aislante, tanto por la parte
que mira a la caldera, como la que da al medio exterior, de espesor 1 cm., y conductividad térmica k*= 0,06
Kcal/m.h.ºC.
b) El coeficiente global de transmisión de calor U, basado en la superficie exterior de la pared Ae, es de la for-
ma,
Ue Ae = 1 ; Ae = 2 π re L ; L = 1 metro
ln re
1 + ri
2 π ri L hcF 2 π L k
Ue = 1 = r 1 = 1 = 62,69 W
2 π re L x ln r e e r e r e 0,02 0,02 m2 °K
2 π re L + ri + ln r
ri hcF k i
+ ln 2
0,015 x 300 0,5 1,5
2 π ri L hcF 2πLk
Régimen estacionario.II.-13
Comprobación: q = Ue Ae (TF - TpF ) = 62,69 x 2 π x 0,02 (100 - 36,53) = 500 W
m
*****************************************************************************************
Un conductor eléctrico de 1 mm de diámetro, se recubre con una capa de un aislante plástico de espesor
e = 2 mm, k1 = 0,5 W/m°C. El hilo está rodeado de aire con una temperatura ambiental TF = 25°C y un
coeficiente de convección h c = 10 W/m2 °K siendo la temperatura del conductor de 100°C.
Determinar
a) El calor disipado por unidad de tiempo, con aislamiento y sin él.
b) La temperatura exterior del aislamiento
Se supondrá que la temperatura del hilo no se ve afectada por la presencia del aislamiento.
_________________________________________________________________________________________
RESOLUCION
a.1) Calor transferido por unidad de longitud, con aislamiento,
r0 = rnúcleo + espesor aislamiento = (0,5 x 10-3)+ (2 x 10-3) = 2,5 x 10-3 metros
10 x {2,5 x 10-3 }
N° de Biot: Bi = hc0 r0 = = 0,05 < 0,1 (R.T.D.)
k1 0,5
Como Bi < 1, la presencia del aislamiento hará aumentar la transferencia de calor del hilo.
0,5 W
rC = k1 =
m°K
= 0,05 m
hc0 10 W
m2°K
2 2
T = TF + E L {L - x + k } = Tmáx - E x = T1 + E (L2 - x2 )
k 2 2 L hcF 2k 2k
1 m x 16.000 m
3.600 seg
Re = Lνu = = 225.430
-6 m2
17,46 10 seg
x
q Te - Ti 35 - 22
= = = 13 = 21,70 Kcal
A R1 + R2 + R3 + R4 + R5 1 + 0,1 0,15 0,0125 0,5989 h m2
+ + + 1
12,55 1,1 0,6 0,4 6,8
Temperatura de la pared revocada con yeso,
q
= 6,8 (T pi - 22) = 21,70 Kcal ; Tpi = 25,2°C
A h m2
*****************************************************************************************
Régimen estacionario.II.-16
II.14.- Sea una pared compuesta de dos capas de 20 cm de espesor de forma que uno de sus lados está com-
pletamente aislado, a la temperatura de 50°C.
Sabiendo que las conductividades térmicas de las placas son de 15 y 20 W/m°C y que la primera genera
calor a razón de 1000 W/m3 , determinar las temperaturas en la unión de las paredes y en la pared exterior
_________________________________________________________________________________________
RESOLUCION
a) Temperaturas en la unión de las paredes y en la pared exterior
Régimen permanente
Pared aislada térmicamente; todo el calor sale por la otra pared
Primera pared,
2
Distribución de temperaturas: T1 (x) = - E x + C1 x + C2
2k
Flujo de calor,
∂T
q1 = - k = k E x - k C1 = E x - k C1
∂x k
Condiciones,
Tx=0 = 50°C ; 50°C = - 0 + 0 +C2 ; C2 = 50°C
Para {x = 0} ⇒
∂T
( )x=0 = 0 ; 0 - k C1 = 0 ; C1 = 0
∂x
2
T1 (x) = - E x + 50
2k
Segunda pared,
*****************************************************************************************
Régimen estacionario.II.-17
II.15.- Una tubería de hierro de d e = 102 mm y d i = 92 mm, k = 50 Kcal/h.m°C de una instalación de calefac-
ción, conduce agua a 90°C, atravesando un local cuya temperatura es de 15°C. La tubería está aislada con
una coquilla de material aislante de e = 25 mm y k* = 0,04 Kcal/h.m°C, envuelta en carbón asfáltico de 5
mm de espesor y k´ = 0,12 Kcal/h.m°C.
hcF = 1000 Kcal/h.m2 °C; hce = 8 Kcal/h.m2 °C
Determinar,
a) La pérdida horaria de calor por metro lineal de tubería
b) El coeficiente global de transmisión de calor U.
c) Las temperaturas superficiales de la tubería aislada y del interior de la tubería
d) Comparar los resultados obtenidos con los correspondientes a la pared desnuda
_________________________________________________________________________________________
RESOLUCION
a) Pérdida de calor horaria por metro lineal de tubería
2 π (90 - 15)
= = 39,04 Kcal
ln 51 ln 76 ln 81 h.m
1 + 46 + 51 + 76 + 1
0,046 x 1000 50 0,04 0,12 0,081 x 8
b) Coeficiente global de transmisión de calor U, referido a la sección exterior 2 π r4 L
U(r=r4) = 1 =
r4 + r4 ln r2 + r4 ln r3 + r4 ln r4 + 1
r1 hcF k12 r1 k23 r2 k34 r3 hce
= 1 = 1,023 Kcal
0,081 x ln 51 0,081 x ln 76 0,081 x ln 81 h.m2 °C
0,081 46 51 76 1
+ + + +
0,046 x 1000 50 0,04 0,12 8
c) Temperatura de la pared interior-fluido caliente,
q
= 2 π r1 hcF (TF - TpF ) = 39,04 Kcal
L
39,04 39,04
TpF = TF - = 90°C - = 89,86°C
2 π r1 hcF 2 π x 0,046 x 1000
Temperatura de la pared exterior-fluido exterior:
q
= 2 π r4 hce (Tpe - Te) = 39,04 Kcal
L
39,04 39,04
Tpe = Te + = 15°C - = 24,6°C
2 π r4 hce 2 π 0,081 x 8
x
1 2 π (90 - 15)
U= r = = 7,92 Kcal
r2 ln r2 0,051 0,051 x ln 51 h.m2 °C
r2 + 1 + 1 + 46 + 1
r1 hcF k12 hce 0,046 x 1000 50 8
Régimen estacionario.II.-18
q
L 190,4
T1 = TF - = 90 - = 89,34°C
2 π r1 hcF 2 π x 0,046 x 1000
q
L 190,4
T2 = Te + = 15 + = 89,27°C
2 π r2 hce 2 π x 0,051 x 8
*****************************************************************************************
II.16.- Una esfera hueca de radio interior ri y radio exterior re está calentada eléctricamente por la pared
interior a razón de q0 W/m2 . Por la pared exterior se disipa el calor a un fluido que se encuentra a TF,
siendo hcF el coeficiente de transmisión de calor por convección-radiación y “k” la conductividad térmica
del sólido. Determinar las temperaturas Tpi y TpF sabiendo que ri = 3 cm, r e = 5 cm, h cF = 400 W/m 2 ºC, TF =
100ºC, k = 15 W/mºC, q0 = 105 W/m2 .
_________________________________________________________________________________________
RESOLUCION
Régimen estacionario.II.-19
a) Distribución de temperaturas
∞ π n (b-y)
Sh 1 - (-1)n
T = 2 T0 ∑ a sen π an x =
n=1 Sh π n b π n
a
∞ π n (0,5 - y) ∞
Sh
= 2x 100 ∑ a 2 sen π n x (con: n = 1, 3, 5, 7,...) = 127,3 ∑ ShSh{π(πn (0,5 - y)} sen(π n x)
n=1 Sh π n b πn a
n=1
n x 0,5) n
a
1 (m/seg) x 0,05 m
Re = u νd1 = = 62500 (turbulento)
0,001 Kg/seg.m)
1,25 x 103 (Kg/m3 )
Kg
0,001 seg.m x 0,3 Kcal
η cp Kg.ºC seg
Pr = = x 3600 = 5,4
k 0,2 Kcal hora
h.m.ºC
Como este fluido se calienta, la ecuación de Dittus Boelter toma la
forma,
Nu = 0,023 x Re0,8 x Pr0,4 = 0,023 x 625000,8 x 5,4 0,4 = 310
Régimen estacionario.II.-20
310 x 0,2 Kcal
h.m.ºC
hcF = Nu k = = 1240 Kcal
d1 0,05 m h.m2 .ºC
b) Coeficiente global de transmisión referido a la superficie externa, si el coeficiente de película exterior es
de 10 Kcal/h.m2 .°C,
2 π k´ (T3 - T2 ) 2 π k (T2 - T1 )
q = 2 π r3 hce (Te - T3 ) = r = = 2 π r1 hcF (T1 - TF) =
ln r3 ln rr2
2 1
2 π (Te - TF)
= = 2 π r U (Te - TF)
1 + 1 ln r3 + 1 ln r2 + 1
r3 hce k´ r2 k r1 r1 hcF
Si se toma como referencia la superficie exterior (3), se tiene,
q = 2 π (r U)3 (Te - TF) ⇒ U3 = 1 =
1 + r3 ln r3 + r3 ln r2 + r3
hce k´ r2 k r1 r1 hcF
= 1 = 0,2815 W
1 0,1 0,1 0,1 0,028 0,1 m2 ºC
+ ln + ln +
11,627 0,04 0,025 46,5 0,025 1441,8 x 0,025
U*3 = 1 = 1 = 0,3059 W
1 + r3 ln r3 1 0,1 0,1 m2 ºC
hce k´ r2 + ln
11,627 0,04 0,028
Flujo térmico,
q = 2 π r3 U*3 (Te - T2 ) = 2 π x 0,1 x 0,3059 x {20 - (-40)} = 11,53 W
m
d) Temperatura de la superficie exterior
q 11,53
q = 2 π r3 hce (Te - T3 ) ⇒ T3 = Te - = 20ºC - = 18,42ºC
2 π r3 hce 2 π x 0,1 x 11,627
e) Error cometido en la lectura al situar el termopar en (2) y no en (1)
q = T2 -r T1 = T1 - TF =
r2
T2 - TF
ln r2 1 ln r
1 hcF 2 π r1 1 + 1
2πk 2πk hcF 2 π r1
r
ln r2 ln 28
T2 - TF = q { 1 + 1 } = 11,53 x { 25 + 1 } = 0,137ºC
2πk hcF 2 π r1 2 π x 46,5 1441,8 x 2 π x 0,025
*****************************************************************************************
Régimen estacionario.II.-21
II.19.- Mediante los balances de energía necesarios, hallar la
matriz de temperaturas correspondientes a los puntos A, B, C,
(1) y (2), de la sección transversal de una chimenea, por cuyo
interior circulan gases a 500ºC, estando la pared exterior a
30ºC, sabiendo que el coeficiente de película de los gases es de
50 W/m 2 .ºC, y la conductividad térmica del material de la chi-
menea, k = 0,5 W/m.ºC
____________________________________________________
RESOLUCION
T + 2 TB + T3
TA = 1
4
TB = T 2 + T A + TC + T3
4
TC = 2 T B + 2 T 3 = TB + T3
4 2
T1 + 2 TB + T3 - 4 T A = 0
T2 + T A + TC + T3 - 4 TB = 0
TB + T3 = 2 TC
Bi = hcF ∆x =
50 x 0,5
= 50
k 0,5
Cálculo de T1 ,
2 T2
+ TA - (2 + Bi) T1 + Bi TF = 0 ; T2 + TA - (2 + 50) T1 + (50 x 500) = 0
2
T2 - 52 T1 + TA + 25000 = 0
Cálculo de T2 ,
T1 - 53 T2 + 2 TB + 25000 = 0
T1 + 0 + 0 + 2 T B - 4 TA = - 30 T1 = 494,84ºC
0 + T2 + TC - 4 TB + TA = -30 ⇒ T2 = 489,41ºC
0 + 0 - 2 TC + TB + 0 = - 30 TA = 242,15ºC
T1 - 53 T2 + 0 + 2 TB + 0 = - 25000 TB = 221,80ºC
- 52 T1 + T2 + 0 + 0 + T A = - 25000 TC = 125,94ºC
*****************************************************************************************
Régimen estacionario.II.-22
II.20.- Determinar la distribución de temperaturas en régimen estacionario y las transferencias de calor por
unidad de tiempo, procedentes de las cuatro superficies del cuerpo bidimensional de sección cuadrada indi-
cado en la figura, k = 1 W/mºC, con E = 0. Dos de los límites son isotermos a TB = 200°C y TD = 100°C res-
pectivamente, un tercero está aislado y el cuarto transfiere energía por convección a un fluido con hC = 50
W/m2 °C y TF = 50°C. Dimensiones(20 x 20) cm.
a) Por el método numérico ; b) Por el método de relajación ; c) Por técnicas matriciales ; d) Por el método
de iteración
RESOLUCION
a) Método numérico
Dividimos el cuerpo de forma que los nudos se nume-
ren del 1 al 9 tal como se indica en la figura
Las celdillas son cuadradas con ∆x = ∆y = 10 cm
Los nudos 1, 2, 3, tienen una temperatura igual a
TB = 200ºC
Los nudos 7, 8, 9, tienen una temperatura igual a
TD = 100ºC
Los nudos 4, 5, 6, tienen temperaturas desconocidas
hcF ∆x 50 x 0,1
Bi = = = 5
k 1
T5 - 7 T4 = - 400
El nudo 6 está situado sobre el contorno aislado por lo que, q3→6 + q5→6 + q9→6 = 0
k ∆x T3 - T6 + k ∆y T5 - T6 + k ∆x T9 - T6 = 0 ; T3 + T9 + T5 - 2 T6 = 0
2 ∆y ∆x 2 ∆y 2
200 + 100 + T5 - 2 T6 = 0 ; T5 - 2 T6 = - 150
2
T4 + T6 - 4 T5 = - 300 T4 = 75,5ºC
T5 - 7 T4 = - 400 ⇒ T5 = 128,7ºC
T5 - 2 T6 = - 150 T6 = 139,4ºC
Transferencia de calor por unidad de tiempo y por unidad de espesor del sólido desde la superficie A al fluido,
qA = qF→1 + qF→4 + qF→7 = hCF ∆y {TF - T1 + (TF - T4 ) + TF - T7 } =
2 2
= ( 50 x 0,1) { 50 - 200 + (50 - 75,5) + 50 - 100 } = - 627 W
2 2 m
El signo (-) indica que el calor se desprende del sólido
Régimen estacionario.II.-23
Transferencia de calor por unidad de tiempo y por unidad de espesor del sólido en la superficie B,
T -T T -T T - T6 ∆y
qB = q1→4 + q2→5 + q3→6 + q 1→F = k ∆x { 1 4 + 2 5 + 3 } + h cF (T1 - TF) = 538,8 W
2 ∆y ∆y 2 ∆y 2 m
Transferencia de calor por unidad de tiempo y por unidad de espesor del sólido en la superficie C, (aislada),
qC = 0
Transferencia de calor por unidad de tiempo y por unidad de espesor del sólido en la superficie D,
T7 - T4 T8 - T5 T9 - T6 ∆y
qD = q7→4 + q 8→5 + q 9→6 + q 7→F = k ∆x { + + } + h cF (T7 - TF) = 88,8 W
2 ∆y ∆y 2 ∆y 2 m
Comprobación global,
qneta = qA + qB + qC + qD = - 627,5 + 538,8 + 0 + 88,8 = 0,1 W
m
que es una muy buena aproximación.
b) Técnicas de relajación
T5 - 7 T4 + 400 = R 4
T4 + T6 - 4 T5 + 300 = R 5
T5 - 2 T6 + 150 = R 6
Temperaturas que se pueden presuponer inicialmente
T4 = 80ºC ; T5 = 100ºC ; T6 = 150ºC
de acuerdo con los límites extremos del problema, T F = 50ºC y TB = 200ºC
Cálculo de los residuos,
R4 = T5 - 7 T4 + 400 R4 = 100 - (7 x 80) + 400 = - 60ºC
R5 = T4 + T6 - 4 T5 + 300 R5 = 80 + 150 - (4 x 100) + 300 = 130ºC
R6 = T5 - 2 T6 + 150 R6 = 100 - (2 x 150) + 150 = - 50ºC
Corrección de temperaturas.- Tomamos el mayor residuo en valor absoluto R5 y aumentamos T5 , por ejemplo
en 35ºC, T5 = 135ºC, de forma que el residuo R5 cambie de signo, por lo que los nuevos residuos serán,
Tomamos el mayor residuo en valor absoluto R4 y modificamos T4 en una cantidad que haga R4 positivo, por
ejemplo disminuyendo T4 en 4ºC, T4 = 76ºC ,
y así sucesivamente,
T4 = 76ºC R4 = 3
DisminuimosT6 en 10ºC: T5 = 135ºC R5 = - 24
T6 = 140ºC R6 = 5
T4 = 76ºC R4 = - 4
DisminuimosT5 en 7ºC: T5 = 128ºC R5 = 4
T6 = 140ºC R6 = - 2
Régimen estacionario.II.-24
T4 = 75ºC R4 = 3
DisminuimosT4 en 1ºC: T5 = 128ºC R5 = 3
T6 = 140ºC R6 = - 2
T4 = 75ºC R4 = 3,8
Aumentamos T5 en 0,8ºC: T5 = 128,8ºC R5 = - 0,2
T6 = 140ºC R6 = - 1, 2
T4 = 75,55ºC R4 = - 0,05
Aumentamos T4 en 0,55ºC: T5 = 128,8ºC R5 = 0,35
T6 = 140ºC R6 = - 1, 2
etc.
c) Técnicas matriciales
Los nudos situados en las superficies isotermas a TB y T D tienen temperaturas conocidas, quedando por deter-
minar las temperaturas de los 15 nudos restantes
∆x
k
1
∆x/2
2 0 hc 3
Fluido
T2 + T3 T2 + T3
+ T1 - 2 T0 = 0 + T1 + (Bi) TF - (2 + Bi) T0 = 0
2 2
Nudo (1)
TB
TB + T6 + T + 125 - 4,5 T = 0
2 1
2
1
200 + T6 + T + 125 - 4,5 T = 0
2
2 1
2
TF
T6 + T - 4,5 T = - 225
2 1
6 2
Régimen estacionario.II.-25
Nudo (2)
TB
T1 + TB + T3 + T7 - 4 T2 = 0
T1 + 200 + T3 + T7 - 4 T2 = 0
1
2 3 T1 - 4 T2 + T3 + T7 = - 200
TF
6 7
Nudo 3 : T2 - 4 T3 + T4 + TB = - 200
Nudo 4 : T3 - 4 T4 + T5 + T5 = - 200
Nudo 5,
TB
TB + T10
Aislada + T4 - 2 T5 = 0
2
200 + T10 + T - 2 T = 0
4 5
2
T10
T4 - 2 T5 + = - 100
2
Nudo 7 : T2 + T6 - 4 T7 + T8 + T12 = 0
Nudo 8 : T3 + T7 - 4 T8 + T9 + T13 = 0
T5 + T15
Nudo 10 : + T9 - 2 T10 = 0
2
T10
Nudo 15 : + T14 - 2 T15 = - 50
2
Muchos de los elementos de la matriz A son cero, concentrándose los no nulos en las proximidades de la diago-
nal principal.
Los elementos de la matriz B son las constantes de los segundos miembros.
Temperaturas
(1) = 88,57 (2) = 138,9 (3) = 158,56 (4) = 166,3 (5) = 168,38 (6) = 69,3 (7) = 108,5 (8) = 129,00
(9) = 138,28 (10) = 140,9 (11) = 68,1 (12) = 96,83 (13) = 110,69 (14) = 116,9 (15) = 118,68
c) Técnicas de iteracción Para los datos del primer caso (método numérico)
Régimen estacionario.II.-26
Nudo 4 : T4 = 400 + T5
7 7
Nudo 5 : T5 = 300 + T4 + T6
4 4
T5
Nudo 6 : T6 = 150 +
2 2
T4 = 80ºC
supondremos que: T5 = 100ºC
T6 = 150ºC
130,66
T4 = 400 + = 75,77ºC
7 7
77,75 + 150
2ª etapa de iteracción: T5 = 75 + 40,18 = 128,98ºC
4
128,98
T6 = 75 + = 139,49ºC
2
128,98
T4 = 400 + = 75,57ºC
7 7
3ª etapa de iteracción: 75,57 + 139,49
T5 = 75 + = 128,77ºC
4
128,77
T6 = 75 + = 139,38ºC
2
128,77
T4 = 400 + = 75,538ºC
7 7
4ª etapa de iteracción: 75,57 + 139,38
T5 = 75 + = 128,737ºC
4
128,737
T6 = 75 + = 139,368ºC
2
128,737
T4 = 400 + = 75,534ºC
7 7
5ª etapa de iteracción: 75,538 + 139,368
T5 = 75 + = 128,726ºC
4
128,737
T6 = 75 + = 139,368ºC
2
128,726
T4 = 400 + = 75,532ºC
7 7
75,534 + 139,368
6ª etapa de iteracción: T5 = 75 + = 128,725ºC
4
128,726
T6 = 75 + = 139,363ºC
2
*****************************************************************************************
Régimen estacionario.II.-27
II.21.- Una placa sólida de sección rectangular de dimensiones (70 x 100) cm2 y espesor infinito, tiene en
régimen permanente las temperaturas en sus caras tal como se indica,
(x = 0; T = 50°C), (x = 70 cm; T = 150°C), (y = 0; T = 250°C), (y = 100 cm ; T = 200°C)
Determinar,
a) La distribución de temperaturas
b) Valor de la temperatura en el centro geométrico de la sección
_________________________________________________________________________________________
RESOLUCION
a) La configuración inicial se puede descomponer en otras más simples para poder aplicar el método de super-
posición.
A B C
b) Temperatura en el centro; hay que tener en cuenta que, x= 35; y= 50; x´= 50; y´= 35
Resolviendo se obtiene, TA = 20,03 - 0,076 + 0,0005 - ... = 19,95
TB = 26,7 - 0,101 + 0,0007 - ... = 26,6
TC = 38,17 - 1,57 + 0,104 - 0,0083 + 0,0007 -... = 36,699
Tcentro = 50 + 19,95 + 26,6 + 36,699 = 133,25
*****************************************************************************************
II.22.- Sea un cilindro macizo de 50 cm de radio y 50 cm de altura; la cara superior se encuentra a una tem-
peratura doble que la superficie lateral, y la cara inferior se encuentra aislada térmicamente; en estas condi-
ciones se sabe que la temperatura en el centro es de 120°C.
Determinar
a) La temperatura en la superficie lateral
b) La temperatura en la base aislada térmicamente
_________________________________________________________________________________________
RESOLUCION
Al tener el cilindro una base aislada térmicamente, esto equivale a duplicar su longitud, y estudiar la distribu-
ción de temperaturas en un cilindro de longitud (0,50 x 2 = 1 m), y temperatura en las bases 2 TA.
Régimen estacionario.II.-28
En estas condiciones, el nuevo cilindro se puede descomponer en otros dos, uno a temperatura uniforme TA, y
otro con temperatura en las bases TA, y temperatura 0 en la superficie lateral. Haciendo uso de la gráfica corres-
pondiente, teniendo en cuenta que, L/R= 1/0,5 = 2 y z/L = 0,25/1 = 0,25, resulta que t/T = 0,46; por lo tanto,
1,0
t/T
0,8 L/R=1,0
Calentado en
dos bases
0,6
1,5 L/R=1,0
Calentado en
0,46 una base
1,5
0,4
2
2
3
0,2 4
3
4
0
0 0,2 0,4 0,6 0,8 z/L 1,0
Es como hallar la temperatura TC´ en el centro de la siguiente configuración, que a su vez se descompone en
otras dos, una de temperatura uniforme 82,19°C, y otra con temperatura en las bases de 82,19°C y temperatura
lateral 0°C.
Por lo tanto, como, z/L = 0,5/1 = 0,5, L/R = 2, resulta, t/L = 0,28
Tc´= 82,19 + 0,28 x 82,19 = 105,2°C
*****************************************************************************************
Régimen estacionario.II.-29
II.23.- Un cilindro macizo de 1 m de longitud y 0,5 m de radio se calienta de tal modo que su base superior
posee una temperatura doble que la de su base inferior; la superficie lateral se encuentra a 100°C. En estas
condiciones se ha medido la temperatura en el centro del cilindro obteniéndose un valor de 120°C. ¿Cuál es
la temperatura de las bases?
_________________________________________________________________________________________
RESOLUCION
Por el principio de superposición descomponemos el problema en otros más sencillos,
Mediante la utilización de gráficas para cilindros,
L = 1 = 2 ; centro: z = 0,5
R 0,5 L
Cilindro calentado en ambos extremos: t = 0,28
T
Cilindro calentado en un extremo: t = 0,14
T
Régimen estacionario.II.-30
Para el cilindro de base superior a 100°C, L/R = 1 ; z/L = 0,5, la gráfica, para el caso de cilindro calentado en
una base, proporciona,
t/T = 0,385 ; t = 0,385 x T = 0,385 x 100 = 38,5°C
Para el cilindro de base inferior a -50°C, L/R = 1 ; z/L = 0,5, la gráfica, para el caso de cilindro calentado en
una base, proporciona,
t/T = 0,385 ; t = 0,385 x T = 0,385 x (-50) = -19,25°C
Al añadir los 150°C restados en un principio y superponer las dos soluciones, se obtiene,
Tcentro = 38,5 - 19,25 + 150 = 169,25°C
*****************************************************************************************
II.25.- Determinar la transferencia de calor por unidad de tiempo a través de la viga de la figura, si T1 =
500°C , T2 = 200°C y k = 70 W/m°K.
_________________________________________________________________________________________
RESOLUCION
En la viga se tiene que, M = 13 ; N = 4 ; ∆Tglobal = T1 - T2 = 500 - 300 = 300ºC
Como para determinar N sólo se ha considerado la mitad de la viga, la transferencia de calor en la unidad de
tiempo a través de la unidad de espesor de la misma es,
Q = 2 k N (T1 - T2 ) = 2 x 70 W x 4 x 300ºC = 12.923 W/m
M m.ºK 13
*****************************************************************************************
II.26.- Un tanque esférico de diámetro d=0,5 metros contiene materiales radiactivos y va enterrado en la tie-
rra a una profundidad de 1,25 metros, medidos desde el centro de la esfera hasta la superficie del suelo. La
superficie del tanque se mantiene a una temperatura uniforme T1 =100°C, consecuencia del proceso de
radiación, mientras que la superficie de la tierra está a una temperatura de 15°C.
Sabiendo que la conductividad térmica de la tierra es de 0,8 W/m°C, determinar el calor generado en el tan-
que como consecuencia de los procesos radiactivos que acontecen en su interior.
_________________________________________________________________________________________
RESOLUCION
T2
4 π r1 4 π x 0,25
k F= r = z >> ri = = 3,49 metros
z
1- 1 0,25
Ti
2z 1 -
ri
2 x 1,25
Q = F k (Ti - T2 ) = 3,49 m x 0,8 W (100 - 15)ºC = 237,3 W
m.ºC
*****************************************************************************************
Régimen estacionario.II.-31
II.27.- La sección recta de las paredes de un horno tiene la forma
que se indica en la figura, en la que sus dimensiones exteriores son
(2,8 x 2,8 m2 , el espesor de las paredes L=40 cm, temperatura de las
superficies interiores 1100°C, temperatura de las paredes exteriores
100°C, conductividad térmica k = 0,1 Kcal/h.m°C.
Determinar,
a) La distribución de temperaturas
b) El calor disipado al exterior
_________________________________________________________________________________________
RESOLUCION
Cálculo de la resistencia térmica,
Considerando un cuadrante y aplicando la ecuación de paredes en ángulo
recto,
R = 1 = Lx = Ly =
L2 + L2 Ly Ly
x y 2 x y 2 Lx 2
k[ + + ln ] + arc tg ( ) + arc tg
Ly Lx π 4 Lx Ly π Lx Ly π Lx
1 1 h.m.ºC 1 h.m.ºC
= = = 1,7989 =
k (2 x - 2 ln 2 + 1) 0,1 (2 1 - 2 ln 2 + 1) Kcal 0,5558 Kcal
Lx π 0,4 π
y L2 + L2y L
Eje y = T1 - T0 = + 1 [ln x + y arc tg (Lx )]
Ly π 4 L2x Lx Ly
b) Calor disipado al exterior
Régimen estacionario.II.-32
Los puntos en contacto con la superficie exterior están a 100ºC
T0-0 = T0-1 = T0-2 = T0-3 = .... = T0-14 = 100ºC
Los puntos en contacto con la superficie interior están a 1100ºC
Td-0 = Td-1 = Td-2 = Td-3 = .... = Td-10 = 1100ºC
Asu vez,
T1-a = T-1-a ; T1-b = T-1-b
Nudos interiores: TA + TB + TC + TD = 4 Ti
(11´)
q0-1 = k (T1´ - 100) = q1´-1´ = q2´-2´ = q3´-3´ ...
(10)
Si los cuadriláteros se han construido de forma que la suma de los
dos lados opuestos sean iguales, se cumple,
≅ ∆x
(11´) (1´2´) (2´3´)
≅ ≅
(10) (1´1´) (2´2´) ∆y
por lo que,
q0-1 = k (1100 - 100) (1 + 1 + 1 + 1) = 4000 k
Régimen estacionario.II.-33
METODO DEL FACTOR DE FORMA
0,4
a = 1 m ; a > e/5 F = A1 + 2,16 (a + b + c) + 1,2 e
e
e = 0,4 m ; e = = 0,08 m ; b = 2 m ; b > e/5 ⇒
5 5 (a,b,c) > e/5
c = 2 m ; c > e/5
T0 a
T0
T1
c
e
T1
C = 0,0266
Nu = C Ren Pr1/3 = = 0,0266 x (150000)0,805 x 60,0266 = 710 (general)
n = 0,805
C = 0,0266
Nu = 1,11 C Ren Pr0,42 = = 1,11 x 0,0266 x (150000)0,805 x 60,42 = 920 (líquidos)
n = 0,805
Régimen estacionario.II.-34
Tomaremos, Nu = 920
hcF d 920 x 0,01
920 = ; hcF = = 467,7 W
k 1,2 m2 .ºK
b) Temperatura de la pared del contenedor exterior e interior
Se considerará conducción monodimensional radial, viniendo los resultados por unidad de longitud del cilindro
1 m. La energía que generan los productos radiactivos se disipa al exterior por convección, atravesando las
paredes del contenedor.
El calor generado por los residuos radiactivos por unidad de longitud del cilindro es,
Q = E V = E x π r2 = 105 W x π x 0,5 2 m2 = 78.500 W m
m3
Temperatura exterior del contenedor TpF
TpF - Tagua TpF - 5ºC
Q= = = 78500 W
m
1 1
2 π re hcF 2π x 0,6 m x 467,7 W
m2 ºK
TPF = 49,5ºC
Temperatura interior del contenedor Tpi
Régimen estacionario.II.-35