Universidad Nacional Experimental Politécnica

“Antonio José de Sucre”

Vicerrectorado Puerto Ordaz
Departamento de Ingeniería Mecánica
Sección de Termofluidos

Problema Tipo Para el Primer Parcial de Transferencia de Calor

1. Un reactor nuclear consiste en un material radiactivo R1 = 100 mm

R 2 = 160 mm
encapsulado en dos cilindros de acero inoxidable. El T 2
hi = 3,5 kW/m2°C
espesor de los cilindros de acero se pueden T1 Ti = 600 °C
despreciar. En la parte interior del cilindro circula k he = 20 W/m2K
Calentador
vapor de agua a 600 °C que genera un coeficiente de Te = 27°C
2
convección de 3,5 kW/m °C, mientras que en la parte R1
exterior hay aire atmosférico que se encuentra a 27 °C hi,
y produce un coeficiente de convección de 20 Ti
2 R2
W/m K. El reactor genera calor por unidad de
volumen a razón de 3 MW/m3 y tiene una he, Te
conductividad térmica de 57 W/m°C. Se pide: (a)
Dibuje el circuito térmico equivalente [2 ptos]; (b)
Determine la temperatura máxima dentro del reactor y las temperaturas en la pared interior y exterior del
mismo [3 ptos]; (c) Calcule el flujo de calor por unidad de longitud del reactor que absorbe el vapor de
agua y el que se pierde en el ambiente [3 ptos].
Datos:
R1 = 100 mm R2 = 160 mm
hi = 3,5 kW/m2°C Ti = 600 °C Q& i T1
he = 20 W/m2K Te = 27°C Ti
q” = 3 MW/m3 k = 57 W/m°C Ri 1
2πR1Lhi
Respuesta:
T máx Q&
El circuito térmico del reactor se muestra en la figura siguiente.
1
Para calcular la temperatura máxima se debe determinar la distribución Re 2π R2 Lhe
de temperatura dentro del elemento radiactivo, por lo cual se debe
Te
resolver la ecuación diferencial de conducción en coordenadas T2
Q& e
cilíndricas, bajo las siguientes consideraciones:
- Estado estacionario.
- Temperatura uniforme en la superficie interior (T1) y exterior (T2).
- Generación de calor uniforme.
- Propiedades del elemento radiactivo son constante.
- La longitud del cilindro es mucho más grande que el diámetro externo del elemento radiactivo, por lo
que se puede suponer que el flujo de calor es radial.
La ecuación diferencia de conducción bajo las consideraciones anteriores se reduce a:
1 d ⎛ dT ⎞ q"
⎜r ⎟+ = 0 (i)
r dr ⎝ dr ⎠ k
Utilizando separación de variable para integrar por primera vez se tiene:
 d ⎛ dT ⎞ q" dT q" c
∫ dr ⎜⎝ r dr ⎟⎠dr = − ∫ k rdr ⇒
dr
=− r+ 1
2k r
(ii)

Integrando nuevamente
q" 2
T =− r + c1 ln r + c2 (iii)
4k
Aplicando las condiciones de borde
q" 2
r = R1 ⇒ T = T1 ⇒ T1 = − R1 + c1 ln R1 + c2 (iv)
4k
q" 2
r = R2 ⇒ T = T2 ⇒ T2 = − R2 + c1 ln R2 + c2 (v)
4k
Restando (iv) – (v),

(T1 − T2 ) − q" (R22 − R12 )

c1 = 4k (vi)
⎛ R1 ⎞
ln⎜⎜ ⎟⎟
⎝ R2 ⎠
Sustituyendo (vi) en (iv) se tiene que:
⎡ ⎤
⎢ (T − T ) − q" R 2 − R12
q" 2 ⎢ 1 2 4k 2
( )⎥
c2 = T1 + R1 − ⎥ ln R (vii)
4k ⎢ ⎛R ⎞ ⎥ 1

⎢ ln⎜⎜ 1 ⎟⎟ ⎥
⎣⎢ ⎝ R2 ⎠ ⎦⎥
La intensidad de flujo de calor (calor por unidad de área) es:
⎛ ⎞
dT
⎜ k (T − T ) − q" R 2 − R 2
q" ⎜ 1 2 4
2 1 ( ) ⎟⎟ 1
q = −k ⇒ q = r −⎜ ⎟r (viii)
dr 2 ⎜ ⎛ R1 ⎞ ⎟
ln ⎜⎜ ⎟⎟
⎜ ⎟
⎝ ⎝ R2 ⎠ ⎠
El flujo de calor hacia el exteriori es
⎡ ⎛ ⎞ ⎤
⎢ ⎜ k (T − T ) − q" R 2 − R 2
⎜ 1 2 ( ) ⎟ ⎥
⎟ 1⎥
∫ q(R2 )dA = q(R2 )2πR2 L Q& e = ⎢ R2 − ⎜
q" 2 1
Q& e = ⇒ 4
⎟ R ⎥ 2πR2 L ⇒
⎢2 ⎛R ⎞
2πR2 L
⎢ ⎜ ln⎜⎜ 1 ⎟⎟ ⎟ 2⎥
⎜ ⎟ ⎥
⎣⎢ ⎝ ⎝ R2 ⎠ ⎠ ⎦

⎛ ⎞
Q& e ⎜ 1 2 4
2 1(
⎜ k (T − T ) − q" R 2 − R 2 ) ⎟
⎟
= πq" R22 − 2π ⎜ ⎟ (ix)
L ⎜ ⎛R ⎞ ⎟
⎜ ln⎜⎜ 1 ⎟⎟ ⎟
⎝ ⎝ R2 ⎠ ⎠

i
Note que el vector intensidad de flujo de calor y el vector normal de la superficie exterior del tubo son paralelos.
El flujo de calor hacia el interiorii es
⎡ ⎛ ⎞ ⎤
⎢ ⎜ k (T − T ) − q" R 2 − R 2
⎜ 1 2 ( ) ⎟ ⎥
⎟ 1⎥
∫ q(R1 )dA = −q(R1 )2πR1 L Q& i = − ⎢ R1 − ⎜
2 1
q"
Q& i = − ⇒ 4
⎟ R ⎥ 2πR1 L ⇒
⎢2 ⎛ R1 ⎞
2πR1L
⎢ ⎜ ln ⎜⎜ ⎟⎟ ⎟ 1⎥
⎢⎣ ⎜ ⎟ ⎥
⎝ ⎝ R2 ⎠ ⎠ ⎦

⎛ ⎞
&
Qi 2
4
(
⎜ k (T − T ) − q" R 2 − R 2
⎜ 1 2 1 ) ⎟⎟
= 2π ⎜ ⎟ − πq" R1
2
(x)
L ⎜ ⎛ R1 ⎞ ⎟
ln ⎜⎜ ⎟⎟
⎜ ⎟
⎝ ⎝ R2 ⎠ ⎠
El calor total generado por reactor cilíndrico es:
Q& Q& e Q& i Q&
Q& = Q& e + Q& i ⇒ =
L L
+
L
⇒
L
(
= πq" R22 − R12 ) (xi)

De circuito térmico se tiene que:

T − Te Q& e
Q& e = 2 ⇒ = 2πR2 he (T2 − Te ) (xii)
1 L
2πR2 Lhe

T −T Q& i
Q& i = 1 e ⇒ = 2πR1hi (T1 − Ti ) (xiii)
1 L
2πR1 Lhe
Igualando (ix) y (xii)
⎛ ⎞
(
⎜ k (T − T ) − q" R 2 − R 2
⎜ 1 ) ⎟⎟
⎟ = 2πR2 he (T2 − Te )
2 2 1
πq" R2 − 2π ⎜
2 4 ⇒
⎜ ⎛R ⎞ ⎟
⎜ ln⎜⎜ 1 ⎟⎟ ⎟
⎝ ⎝ R2 ⎠ ⎠
Sustituyendo valore numéricos y simplificando
761,9974T1 − 782,1036T2 + 85407,1821 = 0 ⇒ T1 − 1,0264T2 = −112,0833 (xiv)
Igualando (x) y (xiii)
⎛ ⎞
⎜ 1 (
⎜ k (T − T ) − q" R 2 − R 2 ) ⎟
⎟
⎟ − πq" R1 = 2πR1hi (T1 − Ti )
2 2 1
2π ⎜ 4 2

⎜ ⎛R ⎞ ⎟
⎜ ln⎜⎜ 1 ⎟⎟ ⎟
⎝ ⎝ R2 ⎠ ⎠
− 2961,112298T1 + 761,997405T2 + 1381631,136 = 0 ⇒ − T1 + 0,2573T2 = −466,5919 (xv)
Sumando (xiv) y (xv) se tiene
− 0,7691T2 = −578,6752 ⇒ T2 = 752,41°C

ii
Note que el vector intensidad de flujo de calor y el vector normal de la superficie interior del tubo son anti-paralelos.
Sustituyendo T2 en (xiv
T1 − 1,0264(752,41) = −112,0833 ⇒ T1 = 660,19°C
Con los valores de T1 y T2 se calculan las constantes c1 y c2.
c1 = 632,938°C y c2 = 2249,162°C
La temperatura máxima se obtiene haciendo
dT q" c 2kc1 2(57W / m°C )632,938°C
=0 ⇒ − rmáx + 1 = 0 ⇒ rmáx = = ⇒
dr r = rmáx 2k rmáx q" 3x10 6 W / m 3

rmáx = 0,155m(155mm )
Luego la temperatura máxima es:
q" 2
Tmáx = − rmáx + c1 ln rmáx + c2 ⇒ Tmáx = 753,04°C
4k
El flujo de calor al exterior e interior se calcula mediante las ecuaciones (xii) y (xiii) respectivamente
Q& e Q& e
= 2πR2 he (T2 − Te ) ⇒ = 14585,23W / m
L L
Q& i Q& i
= 2πR1hi (T1 − Ti ) ⇒ = 132364,72W / m
L L

2. La Figura muestra un chip (microprocesador) de computadora que esta empotrado en un elemento

cerámico (aislado térmicamente). Para mejorar la disipación de calor, se le coloca al chip por la cara no
empotrada, una placa de aluminio con un banco de 225 aletas recta cilíndricas de 75 mm de longitud y 5
mm de diámetro. La temperatura del aire alrededor de la aleta es de T∞ = 24 °C. En la unión entre el
Chip y la placa de aluminio hay una resistencia de
contacto que se puede estimar en 8x10-5 m2K/W. El
chip disipa 100 W a plena potencia. Se pide: (a)
Dibuje el circuito térmico equivalente [1 ptos]. (b) L
Determine la eficiencia de cada aleta y la eficiencia Empotramiento B B = 100 mm
Cerámico H = 100 mm
global del banco de aleta [2 ptos]; (c) Calcule el Lc = 5 mm
Lc
coeficiente de convección que debe general el La = 20 mm
Chip
sistema de enfriamiento con aire para garantizar L = 75 mm
D = 5 mm
que la temperatura máxima dentro del chip no kc = 1,2 W/mK
supere los 80 °C [3 ptos]; (d) Cuanto sería la ka = 165 W/mK
Rc'' = 8x10-5 m2/WK
temperatura máxima en el chip, si se elimina la H kc ka
placa con el banco de aleta y el coeficiente de
convección es igual al calculado en el punto Tmáx T1 T2 Tb
anterior [2 ptos]. Considera que el flujo de calor por T∞ = 27 °C
D
'' Tmáx
conducción es unidimensional para todos los Rc La
h=?
elementos del sistema y que puede usar el modelo
de aleta adiabática corrigiendo la longitud de la aleta.
Datos:
N = 225 aletas L = 75 mm (0,075 m) D = 5 mm (0,005 m)
T∞ = 24 °C Rc = 8x10-5 m2K/W Q& = 100 W
B = 100 (0,1 m) H = 100 mm (0,1 m) Lc = 5 mm (0,005 m)
La = 20 mm (0,02 mm) kc = 1,2 W/m°C ka = 165 W/m°C
Tmáx = 80 °C
Respuesta:
La Figura muestra el circuito térmico equivalente del sistema.
1
(BH − N π
4 D 2 )h
Tmáx T1 T2 Tb
Q&
Rchip Rc'' La T∞
BH ka BH

1
Nη f A f h

Al considerar al chip como una pared plana en estado estacionario, con propiedades del material constantes y
que el flujo es unidimensional, se tiene que la ecuación de conducción es:
d 2T q′ dT q′ q′ 2
2
=− ⇒ = − x + c1 ⇒ T =− x + c1 x + c2 (i)
dx kc dx kc 2k c
Aplicando las condiciones de borde
dT q′
x =0⇒ =0 ⇒ 0=− (0) + c1 ⇒ c1 = 0
dx x =0 kc

q′ 2 q′L2c
x = Lc ⇒ T = T1 = 0 ⇒ T1 = − Lc + c2 ⇒ c2 = T1 +
2k c 2k c
La distribución de temperatura es:
q′ 2
T = T1 +
2k c
(
Lc − x 2 ) (ii)

La temperatura máxima ocurre en x = 0, entonces:

q′L2c
Tmáx = T1 + (iii)
2k c
El calor generado, considerando una generación uniforme, es:
Q&
Q& = ∫∫∫ q ′dV = q ′V ⇒ Q& = q′Lc BH ⇒ q′ = (iv)
V
Lc BH
Sustituyendo (iv) en (iii)
 Q& Lc
Tmáx = T1 + (v)
2k c BH
O también
T −T
Q& = máx 1 (vi)
Lc
2k c BH
De (vi) se obtiene que:
Lc
Rchip = (vii)
2k c BH
La temperatura Tb se obtiene del circuito térmico
Tmáx − Tb ⎛ Lc R′′ La ⎞
Q& = ⇒ Tb = Tmáx − Q& ⎜⎜ + c + ⎟⎟ (viii)
Lc R′′ La ⎝ 2k c BH BH k a BH ⎠
+ c +
2k c BH BH k a BH
Sustituyendo valores numéricos
Lc 0,005m
= = 0,2083°C / W
2k c BH 2(1,2W / m°C )(0,1m )(0,1m )

Rc′′ 8 x10 −5 m 2 K / W
= = 0,008 K / W = 0,008°C / W
BH (0,1m )(0,1m )
La 0,02m
= = 0,0121°C / W
k a BH (165W / m°C )(0,1m )(0,1m )
Luego
Tb = 57,16 °C
La eficiencia de una aleta de sección transversal constante, en éste caso una aleta cilíndrica, es:
tanh (mLe )
ηf =
mLe
Donde
4h D
m= y Le = L +
ka D 4
Luego
⎛⎛ D ⎞ 4h ⎞⎟
tanh⎜⎜ ⎜ L + ⎟
⎝⎝ 4 ⎠ k a D ⎟⎠
ηf =
⎛ D ⎞ 4h
⎜L+ ⎟
⎝ 4 ⎠ ka D
Sustituyendo valores numéricos
 ⎛⎛ 0,005 ⎞ 4h ⎞
tanh⎜⎜ ⎜ 0,075 + ⎟ ⎟
ηf = ⎝ ⎝ 4 ⎠ (165 )(0, 005 ) ⎟
⎠ ⇒ ηf =
(
tanh 0,1679 h )
⎛
⎜ 0,075 +
0,005 ⎞
⎟
4h ( 0,1679 h ) (ix)

⎝ 4 ⎠ (165)(0,005)
Por otra parte, del paralelo del circuito térmico se tiene:

Q& =
Tmáx − Tb
[( )
= BH − π4 ND 2 + Nη f A f h(Tb − T∞ ) ]
[(BH − ]
(x)
)
4 ND h + Nη f A f h
π 2 −1

Donde
⎛ D⎞
A f = πDLe = πD⎜ L + ⎟ = 1,1977 x10 −3 m 2
⎝ 5⎠
Sustituyendo los valores numéricos y la ecuación (ix) en (x) se tiene:
⎡
(
100 = ⎢ 5,5821x10 −3 + 1,60502
tanh 0,1679 h ⎤
) (
⎥ h(57,16 − 24)
) ⇒
⎣ h ⎦
3015,68 = 5,5821h + 1605,02 h tanh 0,1679 h ( ) (xi)
La ecuación se resuelve numéricamente y se obtiene que:
h = 12,16 W/m2°C
Sustituyendo h en (ix) se tiene que la eficiencia de cada aleta es:
η f = 0,8995(89,95% )
La eficiencia global del banco de aleta es:
Q& Q&
ηo = = ⇒
[
AT h(Tb − T∞ ) (BH − π4 ND 2 )h + NπDL h(Tb − T∞ ) ]
100W
ηo = ⇒
[((0,1m)(0,1m) − ) ](
4 225(0,005m ) + 225π (0,005m )(0,075m ) 12,16W / m °C (57,16 − 24 )°C
π 2 2
)
η o = 0,9163(91,63% )
El circuito térmico para el caso de que no se instale el banco de aletas es el mostrado en la siguiente Figura
Tmáx T1
Q& T∞
Rchip 1
BHh
De éste circuito térmico se tiene que:
Tmáx − T∞ ⎛ Lc 1 ⎞
Q& = ⇒ Tmáx = T∞ + Q& ⎜⎜ + ⎟⎟ ⇒
Lc 1 ⎝ 2k c BH HBh ⎠
+
2k c BH HBh
 ⎛ 0,005m 1 ⎞
Tmáx = 24°C + 100W ⎜⎜ + ⎟⎟
(
⎝ 2(1,2W / m°C )(0,1m )(0,1m ) (0,1m )(0,1m ) 12,16W / m °C
2
) ⎠
Tmáx = 867,20°C
 La figura muestra la sección
transversal de un horno. El k1 = 0,13 w/mºK
k2 = 13 – 0,01T [w/mºK]
refractario tiene una conductividad con T en ºC
térmica K1 = 0,13 w/mºK, y la k1 k2=f(T) Ti = 1000 ºC
hi = 600 w/m2ºK
pared externa, la conductividad Te = 27 ºC
térmica varía con la temperatura he = 15 w/m2ºK

siguiendo la relación K2 = 13-

T1 T2 T3
0,01T [w/mºK], donde T se
hi , Ti
expresa en ºC. Considerando que
las propiedades del fluido externo he , Te
e interno son uniforme, determine:
(a) Dibuje el circuito térmico
equivalente. (b) La temperatura en L1 = 200 mm L2 = 100 mm

cada una de las interfases de los

materiales (T1, T2, T3). (c) El flujo
de calor que circula a través de la
pared del horno.

DATOS:
K1 = 0,13 w/mºK
K2 = ko-BT w/mºK
ko = 13 w/mºK
B= 0,01 w/m^2ºK
L1 = 0,2 m Ti T1 T2 T3 Te
L2 = 0,1 m L1 1
Ti = 1000 ºC
1 R2
hi = 600 w/m^2ºK
k1A he A
hi A
Te = 27 ºC
he = 15 w/m^2ºK

Aplicando la Ley de Fourier para la sección de la pared que tien k variable (k2)

L T3 T3
⎡B 2 ⎤
= − (k o − BT ) ∫T (BT − k o )dT ⇒ ⎢⎣ 2 T − k oT ⎥⎦ T = qL 2
2
dT dT
q = −k2 = ctte ⇒ q ∫ dx =
dx dx 0 2 2

Luego

q=
1 ⎡B 2
⎢
L2 ⎣ 2
( ) ⎤
T3 − T 22 − k o (T3 − T 2 )⎥ =
1 ⎡
⎢ ko −
B ⎤
T 2 + T3 ⎥ (T 2 − T3 ) ( )
⎦ L2 ⎣ 2 ⎦

( )
Haciendo B
k2 = ko − T2 + T3 (i )
2
Se tiene: k
q = 2 ⎛T + T ⎞
⎜ 2 ⎟
L 2 ⎝ 3 ⎠
El calor total es: Ak & k2 ⎛
& = qA = 2 ⎛T + T ⎞ ⇒ Q = ⎞
Q ⎜ 2 ⎟ ⎜T2 + T3 ⎟ (ii)
L2 ⎝ 3 ⎠ A L2 ⎝ ⎠
Por otra parte, del circuito se tiene:

Q& (T3 − Te ) Q& (Ti − T2 )

= (iii) = (iv)
A 1 / he A 1 / h1 + L1 / k1
Igualando (iii) y (iv) y depejando T2

T 2 = T1 − (1 / h1 + L1 / k 1 )h e (T3 − T3 ) (v )

Para resolver el problema se sigue del siguiente algoritmo:

1) Se supone una temperatuta T3
2) Con la ecuación (v) se determina T2
3) Se calcula k2 con la ecuación (i)
4) se calcula Q/A con la ecuación (ii)
5) se calcula Q/A con la ecuación (iii)
6) Se calcula el error = Q/A(ec-ii)-QA(ec.iii)
7) se repite desde (1) hasta (7), interpolando en las dos ultimas temperatura calculadas, hasta
que el error sea nulo
8) Se calcula T1 con la siguiente ecuación

1 ⎛ Q& ⎞
T1 = T i − ⎜⎜ ⎟⎟ ( vi )
hi ⎝ A⎠
Los resultados de la iteración se muestra en la tabla siguiente:

Iteración T3 [ºC] T2[ºC] k2 [w/mºK] (Q/A) I (Q/A) IV Error

1 65,000 122,127 12,0644 6892,0007 570,0000 6322,0007
2 67,000 75,923 12,2854 1096,2343 600,0000 496,2343
3 67,170 71,987 12,3042 592,7074 602,5554 -9,8480
4 67,167 72,064 12,3038 602,5204 602,5057 0,0147
5 67,167 72,064 12,3038 602,5057 602,5057 0,0000

Q/A = 602,506 w/m

T1 = 998,996 ºC
T2 = 72,064 ºC
T3 = 67,167 ºC
La figura muestra un reactor nuclear, cuyo
elemento combustible (torio) es de forma de un
cilindro hueco de radios interno R1 = 200 mm y T3
Refractario
radio externo R2 =400 mm, y tiene una T2 Torio
conductividad térmica de k2 = 25 w/mºK y una
T1
generación de calor uniforme de 10E8 w/m^3 . El
torio está empotrado entre dos placas de acero R0
especial cuyo espesor se pueden considerar R1 hi,
despreciable. La parte exterior del torio se cubre Ti
con un aislamiento refractario, para evitar la R2 k1
perdida de calor al exterior. El recubrimiento
refractario tiene un espesor de 150 mm y una R3 k2 he, Te
conductividad térmica de k3 = 0,15 w/mºK. En la k3
parte interior de torio, soldado al tubo de acero
interior, se montan un sistema de aletas de acero
(k1 = 15,2 w/mºK), rectas de sección transversal
constante, espesor 4 mm y longitud 100 mm. La
soldadura de la aleta con el tubo interior tiene una
resistencia de contacto de 4E-6 ºKm^2/w . El aire
externo tiene una temperatura de Te = 27 ºC y un
coeficiente de convección de he = 10 w/m2ºK,
mientras que el vapor de agua se encuentra a una
temperatura de Ti = 600 ºC y un coeficiente de
convección de hi =600 w/m2ºK. Se pide: (a) Dibuje
el circuito térmico equivalente. (b) Determine las
temperaturas T1, T2, T3. (c) Determine la máxima
temperatura que ocurre en el torio y el punto donde
ocurre. (d) Calcule el flujo de calor total que genera
el torio, el calor que se pierde en el aire y el calor
que absorbe el agua.

DATOS:
R0 = 0,1 m Te = 27 ºC
R1 = 0,2 m he = 10 w/m^2ºK
R2 = 0,4 m q' = 4,00E+05 w/m^3
R3 = 0,55 m Rc'' = 4,00E-06 m^2ºK/w
k1 = 15,2 w/mºK N= 16
k2 = 25 w/mºK t= 0,004 m
k3 = 0,15 w/mºK
Ti = 600 ºC
hi = 600 w/m^2ºK

Primeramente se obtendrá una ecuación

para la conducción dentro del torio.

Consideraciones:
* Estado estacionario
* k2 es constante
* Temperatura en la superficie interior y exterior Q& aire
se consideran uniforme T2 T3 Te
* Generación de calor es uniforme Rt ,e ln (R3 / R 2 ) 1
Q& 2π k 3 L 2πR3 Lh e
Tmáx

T1

Q& agua Rc''

Tb 1
R t ,i Ti
NtL N η f A f hi

Luego, integrando la ecuación de conducción en forma cilídrica 1

se tiene: (2πR1 − Nt )Lhi

1 d ⎛ dT ⎞ q' dT q' 2 dT q' C

⎜r ⎟+ ⇒r =− r + C1 ⇒ =− r+ 1 (i )
r dr ⎝ dr ⎠ k 2 dr 2k 2 dr 2k 2 r
q' 2
T =− r + C1 ln r + C 2 (ii )
4k 2
q ' R12
r = R1 ⇒ T = T1 ⇒ T1 = − + C1 ln R1 + C 2 (iii )
4k 2
q ' R22
r = R2 ⇒ T = T2 ⇒ T2 = − + C1 ln R2 + C 2 (iv)
4k 2
Restando (iii) y (iv) y despejando C1 y luego se obtiene C2:

(T1 − T2 ) + q'
4k 2
(
R12 − R 22 ) q ' R 22
C1 = ⇒ C 2 = T2 + − C 1 ln R 2
ln ⎛⎜ 1 ⎞
R 4k 2
⎟
⎝ R2 ⎠
Haciendo
q'
4k 2
(
R12 − R 22 )
(T1 − T 2 ) + D
D= ⇒ C1 =
ln ⎛⎜ R1 ⎞
⎟ ln ⎛⎜ R1 ⎞
⎟
⎝ R2 ⎠ ⎝ R2 ⎠
q ' R 22
B= ⇒ C 2 = T2 + B − C1 ln R 2
4k 2

La temperatura máxima se obtiene de la siguiente forma:

dT q' C 2C1 k 2
=− rmáx + 1 = 0 ⇒ rmáx =
dr r = rmáx 2k 2 rmáx q'
q' 2 q' 2 ⎛r ⎞
Tmáx = − rmáx + C1 ln rmáx + C 2 = − rmáx + C1 ln ⎜⎜ máx ⎟⎟ + T2 + B ⇒
4k 2 4k 2 ⎝ R2 ⎠
q' 2 ⎛r ⎞
Tmáx − T2 = − rmáx + ln ⎜⎜ máx ⎟⎟ + B
4k 2 ⎝ R2 ⎠
El Flujo de calor por unidad de área es:
Q&
= r − 2 1 ⇒ Q& aire = q( R2 )(2πR2 L) ⇒ aire = πq' R22 − 2πk 2C1
dT q' kC
q = −k 2 (v)
dr 2 r L
Por otra parte, como la generación de calor es uniforme. Se tiene:
&
2 (1
L
)
& = q'V = πq' R 2 − R 2 L ⇒ Q = πq' R 2 − R 2
Q 2 1 ( ) (vi)

Q/L = 150796,4474 w/m

Luego Q& Q& agua Q& aire Q& agua Q& Q& aire
= + ⇒ = − (vii)
L L L L L L
Del circuito electrico se tiene:
T2 − Te Q T2 − Te
Q a ire = ⇒ a ire = (viii)
ln ⎛⎜ 3 ⎞ ln ⎛⎜ 3 ⎞
R L R
⎝ R 2 ⎟⎠ 1 ⎝ R 2 ⎟⎠ 1
+ +
2πLk 3 2πR3 Lhe 2πk 3 2πR3 he

T1 − Ti
Qagua = −1
⇒
⎡⎛ R '' ⎞
−1
⎤
⎢⎜ c + 1 ⎟ + (2πR1 − Nt )Lhi ⎥
⎢⎜⎝ NtL Nη f A f hi ⎟
⎠ ⎥
⎣ ⎦
Qagua T1 − Ti
= −1
(ix)
L ⎡⎛ ⎞
−1
⎤
⎢⎜ '' ⎟ ⎥
⎢⎜ Rc + 1 ⎟ + (2πR1 − Nt )hi ⎥
⎢⎜ Nt ⎛A ⎞ ⎟ ⎥
⎢⎜ Nη f ⎜ f ⎟ hi ⎟ ⎥
⎝ L ⎠
⎢⎣⎝ ⎠ ⎥⎦
Para evaluar la eficiencia de la aleta se utilizará el concepto de aleta adiabática con la correccón
en la punta. Entonces:

L c = (R 1 − R 0 ) + t / 2 Lc = 0,102 m
Af
A f = 2 LLc ⇒ = 2 Lc Af/L = 0,204 m
L
2 hi La eficiencia es muy pequena, por lo
m = m = 140,488 1/m que no justificaría instalarla en la
k 1t
parte interior del tubo, ya que se
tgh (mL c ) aumenta el costo y no se mejora
ηf = nf = 0,070 6,98% significativamente la transferencia de
mL c
calor.

Evaluando numéricamente las constantes y las resiestencia equivalentes de las ecuaciones (viii) y (ix)
B = 640,000
D = 692,494

ln ⎛⎜ 3 ⎞
R
⎟
⎝ R2 ⎠ + 1
= 0,3668
2π k 3 2π R 3 h e
 −1
⎡⎛ ⎞
−1
⎤
⎢⎜ '' ⎟ ⎥
⎢⎜ Rc + 1 ⎟ + (2πR − Nt )h ⎥ = 0,0012
⎢⎜ Nt ⎛ A ⎞ ⎟ 1 i
⎥
⎢⎜ N η f ⎜ f ⎟ hi ⎟ ⎥
⎢⎣⎝ ⎝ L⎠ ⎠ ⎥⎦

Como C1 y C2 son funciones de la temperatura se debe resolver el problema de forma iterativa.

El algoritmo requerido es el siguiente:
1) Se supone una temperatura T2
2) Con la ecuación (viii) se calcula Qaire/L
3) Con la ecuación (vii) se calcula Qagua/L
4) Con la ecuación (ix) se clacula T1
5) Con T1 y T2 se evalúan las constantes C1
6) Con la ecuación (v) se calcula Qaire/L
7) Se calcula el error = (Qaire/L)(ec-viii)-(Qaire/L)(ec-v)
8) se repite desde (1) hasta (8), interpolando en las dos ultimas temperatura calculadas, hasta
que el error sea nulo

T2 [ºC] Qaire/L (ec-viii) Qagua/L T1 C1 Qaire/L (ec-v) Error

1000,000 2652,475 148143,972 774,064 1018,451 41084,015 -3,843154E+04
1500,000 4015,515 146780,932 772,462 1742,109 -72587,922 7,660344E+04
1167,043 3107,847 147688,600 773,528 1260,215 3107,847 2,683009E-11
1167,043 3107,847 147688,600 773,528 1260,215 3107,847 -3,228706E-11

Luego se tiene que:

T1 = 773,528 ºC
T2 = 1167,043 ºC
Qaire/L = 3107,847 w/m 2,06%
Qagua/L = 147688,600 w/m 97,94%
Q/L = 150796,447 w/m
Q aire T3 − Te 1 ⎛ Q aire ⎞
= ⇒ T3 = Te + ⎜ ⎟ ( x)
L 1 2πR3 he ⎝ L ⎠
2πR3 he
T3 = 116,933 ºC
Rmáx = 0,397 m (397) mm
Tmáx = 1176,928 ºC
El calor que llega al agua por la sección sin aleta es:
Q agua ,b T1 − Ti
=
L 1 Qagua,b/L = 124173,9055 w/m 84,08%
(2π R1 − Nt )hi
El calor que va al agua y pasa por el sistema de aleta es:
Q agua , f Q agua Q agua ,b
= − Qagua,f/L = 23514,695 w/m 15,92%
L L L