Problema

Transporte de calor en una pared cilíndrica con transferencia de calor por convección en ambos
lados.
Consideremos una pared cilíndrica de radio interno r0 , radio externo r1 y longitud L . El

coeficiente de transferencia de calor interno es ha y el externo hb ; asimismo las temperaturas

interna y externa del medio son Ta y Tb respectivamente. Encuentre

a) La densidad de flujo de calor en las partes interna y externa

b) Las temperaturas en las interfases solido-fluido
c) El perfil de temperatura

Solución
El esquema es el siguiente

Nota: Para efectos del balance de energía no se requiere visualizar la altura del cilindro

Este tipo de problemas con convección en uno o ambos extremos para una pared de conductividad
constante se aborda del siguiente modo:

Transferencia de calor por convección en la parte interna:

q0 = ha (Ta − T0 ) (1)

Balance de energía en el material solido

d (rqr )
=0 (2)
dr

Integrando la ecuacion (2)

rqr = C = r0 q0 = r1q1 (3)

Es decir, el producto del radio por la densidad de flujo de calor es constante. Por conveniencia
tomamos la siguiente igualdad:
rqr = r0 q0 (4)

De la ecuacion (4) se obtiene

dT 1
−k = r0 q0  (5)
dr r

Integrando (5)

T1 r
r0 q0 1 dr
 dT = −
T0
k r0 r
 (6)

Cuyo resultado es

r0 q0
T0 − T1 =  ln(r1 / r0 ) (7)
k

Transferencia de calor por convección en la parte externa:

q1 = hb (T1 − Tb ) (8)

Ahora sumamos cada una de las ecuaciones del siguiente modo:

De (1)

q0
Ta − T0 = (9)
ha

De (7)

r0 q0
T0 − T1 =  ln(r1 / r0 ) (10)
k

De (8)

q1
T1 − Tb = (11)
hb

Al sumar (9), (10) y (11)

q0 r0 q0 q
Ta − Tb = +  ln(r1 / r0 ) + 1 (12)
ha k hb

En la ecuacion (12) se debe eliminar el valor de q1 mediante la ecuación (3)

 r0
q1 = q0 (13)
r1

Reemplazando (13) en (12)

q0 r0 q0 rq
Ta − Tb = +  ln(r1 / r0 ) + 0 0 (14)
ha k r1hb

De la ecuacion (14) se despeja el valor de la densidad de flujo de calor interno q0

r0−1 (Ta − Tb )
q0 = (15)
1 1 1
+  ln(r1 / r0 ) +
ha r0 k hb r1

La ecuacion (15) permite calcular la densidad de flujo de calor en la superficie interna. Para el
cálculo de la densidad de flujo de calor en la superficie externa se hace uso de la ecuacion (13).

Para el cálculo del perfil de temperatura, previamente debe calcularse las temperaturas T0 y T1

La temperatura T0 se obtiene a partir de (9), es decir:

q0
T0 = + Ta (16)
ha

La Temperatura T1 se obtiene a partir de (11), es decir:

q1
T1 = + Tb (17)
hb

Luego al aplicar el balance de energía y resolviendo de manera similar al problema anterior, se

obtiene:

(T0 − T1 )
T= ln(r / r0 ) + T0 (18)
ln(r0 / r1 )

Problema
Transporte de calor en una pared cilíndrica con transferencia de calor por convección en ambos
lados.
Consideremos una pared cilíndrica de radio interno Ra , radio externo Rb y longitud L . El

coeficiente de transferencia de calor interno es ha y el externo hb ; asimismo las temperaturas

interna y externa del medio son Ta y Tb respectivamente. Considerando que la conductividad

térmica es variable y está dada por k = a + bT , encuentre:

a) El perfil de temperatura
b) El perfil de la densidad de flujo de calor
c) Las temperaturas en las interfases solido-fluido

Datos

Solución

Este tipo de problemas se resuelve usando un procedimiento numérico que consiste en la solución
simultanea del balance de energía y la ley de Fourier.

Balance de energía en el material solido

d (rqr )
=0 (2)
dr

Ley de Fourier

dT
qr = −(a + bT ) (3)
dr

Transferencia de calor por convección en la parte interna

q0 = ha (Ta − T0 ) (4)

Transferencia de calor por convección en la parte externa

q10 = hb (T10 − Tb ) (5)

Derivando la ecuacion (1)

dqr 1
+ qr = 0 (6)
dr r

Discretizando (6)

1
qi − qi −1 + qi = 0 (8)
i

Expandiendo (8)

1
i = 1 q1 − q0 + q1 = 0 con q0 = ha (Ta − T0 )
1

1
i = 2 q2 − q1 + q2 = 0
2

Y así sucesivamente

1
i = 10 q10 − q9 + q10 = 0
10

Discretizando la ecuacion (3)

Ti − Ti −1
qi + (a + bTi ) =0 (9)
r

Expandiendo la ecuacion (9)

T1 − T0
i = 1 q1 + (a + bT1 ) =0
r

T2 − T1
i = 2 q2 + (a + bT2 ) =0
r

Y así sucesivamente

T10 − T9 q
i = 10 q10 + (a + bT10 ) = 0 con T10 = 10 + Tb
r hb

El programa respectivo es el siguiente

f(q1) = q1 - q0 + (1 / 1) * q1
f(q2) = q2 - q1 + (1 / 2) * q2
f(q3) = q3 - q2 + (1 / 3) * q3
f(q4) = q4 - q3 + (1 / 4) * q4
f(q5) = q5 - q4 + (1 / 5) * q5
f(q6) = q6 - q5 + (1 / 6) * q6
f(q7) = q7 - q6 + (1 / 7) * q7
f(q8) = q8 - q7 + (1 / 8) * q8
f(q9) = q9 - q8 + (1 / 9) * q9
f(q10) = q10 - q9 + (1 / 10) * q10

f(T0) = q1 + (a + b * T1) * (T1 - T0) / DR

f(T1) = q2 + (a + b * T2) * (T2 - T1) / DR
f(T2) = q3 + (a + b * T3) * (T3 - T2) / DR
f(T3) = q4 + (a + b * T4) * (T4 - T3) / DR
f(T4) = q5 + (a + b * T5) * (T5 - T4) / DR
f(T5) = q6 + (a + b * T6) * (T6 - T5) / DR
f(T6) = q7 + (a + b * T7) * (T7 - T6) / DR
f(T7) = q8 + (a + b * T8) * (T8 - T7) / DR
f(T8) = q9 + (a + b * T9) * (T9 - T8) / DR
f(T9) = q10 + (a + b * T10) * (T10 - T9) / DR
q0 = ha * (Ta - T0)
Ta = 200
Tb = 80
RA = 0.20
RB = 0.40

DR = (RB - RA) / 10
a = 40
b = -0.1
T10 = q10 / hb + Tb
ha = 50
hb = 120

q1(0) = 100
q2(0) = 100
q3(0) = 100
q4(0) = 100
q5(0) = 100
q6(0) = 100
q7(0) = 100
q8(0) = 100
q9(0) = 100
q10(0) = 100

T0(0) = 80
T1(0) = 80
T2(0) = 80
T3(0) = 80
T4(0) = 80
T5(0) = 80
T6(0) = 80
T7(0) = 80
T8(0) = 80
T9(0) = 80

Resultados

q0 5442.457 T0 91.15086
q1 2721.228 T1 89.39863
q2 1814.152 T2 88.23484
q3 1360.614 T3 87.36442
q4 1088.491 T4 86.66963
q5 907.0761 T5 86.09171
q6 777.4938 T6 85.59712
q7 680.3071 T7 85.16496
q8 604.7174 T8 84.78128
q9 544.2457 T9 84.43634
q10 494.7688 T10 84.12307

Las temperaturas de la interfase solido fluido son 91.15086 oC y 84.12307 oC respectivamente.

a) Balance de energía en coordenadas esféricas

En base a la ecuacion (1), la ecuacion de balance de energía, está dada por

Consideremos una pared esférica tal como se muestra en la figura, en la cual elaboramos el
balance de energía:

Tomando en consideración la ecuacion (2.1) se tiene:

qr r 4 r 2 − qr r +r
4 r 2 = 0 (1)

Llevando al límite se tiene

qr r +r
4 r 2 − qr 4 r 2
lim r
=0 (2)
r →0 4r

A partir de (2) se obtiene

d (r 2 qr )
=0 (3)
dr

Integrando

r 2 qr = C1 (4)

Volviendo a integrar

C1
T =− + C2 (5)
r

Las condiciones de frontera aplicables a este caso son:

r = R1 T = T1 y r = R2 T = T2

Reemplazando las condiciones de frontera se tiene

(T1 − T2 )
C1 = (6)
1 1
−
R2 R1

C1
C2 = T1 + (7)
R1
Reemplazando (6) y (7) en (5)

T = −(T1 − T2 ) 
1/ r − 1/ R1  + T (8)
1/ R2 − 1/ R1  1

La ecuación (8) representa la ecuacion de distribución de temperatura

La densidad de flujo de calor respecto a la superficie interna puede ser obtenida a parir de la ley
de Fourier

dT C1 1 (T1 − T2 ) 1 (T1 − T2 )
qr = −k = −k = −k  = k 2 (9)
r = R1
dr r = R1 r2 r = R1 R1 (1/ R2 − 1/ R1 )
2
R1 (1/ R1 − 1/ R2 )

Del mismo la densidad de flujo de calor respecto a la superficie externa puede ser obtenida a parir
de la ley de Fourier

dT C1 1 (T1 − T2 ) 1 (T1 − T2 )
qr = −k = −k = −k  =k 2 (10)
r = R2
dr r = R2 r2 r = R2 R2 (1/ R2 − 1/ R1 )
2
R2 (1/ R1 − 1/ R2 )

Sin embargo, el flujo de calor es constante

4 R12 (T1 − T2 ) 4 k (T1 − T2 )

Qr = qr 4 R12 = k  = (11)
r = R1 r = R1 R1 (1/ R1 − 1/ R2 ) (1/ R1 − 1/ R2 )
2

4 R22 (T1 − T2 ) 4 k (T1 − T2 )

Qr = qr 4 R22 = k  = (12)
r = R2 r = R2 R2 (1/ R1 − 1/ R2 ) (1/ R1 − 1/ R2 )
2

Se observa que el flujo de calor es constante.

La temperatura media se obtiene mediante (usar las ecuaciones del anexo)

2  R2 R2

   r sen T (r )drd d  (− cos  )
2
r 2T (r )dr
2
0

 T = =
0 0 R1 R1
(13)
2  R2
2 R2

   r sen drd d  (− cos  )
2 2
0
r dr
0 0 R1 R1

R2 R2
3 3 C1
 T =  r T (r )dr = r (− + C2 )dr
2 2
(14)
R2 − R12
2
R1
R2 − R12
2
R1
r

R2
3 3 r2 r3
 T =  r T (r )dr = − +
2
C ( ) (15)
R2 − R12 R22 − R12
2 1
R1
2 3R1
 T  −T1 1 R1 R2  2 R23 
=  R − 3 R 2
+ 2  (16)
T2 − T1 2 ( R23 − R13 )( R2 − R1 ) 
1 2
R1 

La ecuación (16) es la ecuacion de distribución adimensional de temperatura.

Problema
Transporte de calor en una pared esférica con transferencia de calor por convección en ambos
lados.
Consideremos una pared esférica de radio interno r0 , radio externo r1 . El coeficiente de

transferencia de calor interno es ha y el externo hb ; asimismo las temperaturas interna y externa

del medio son Ta y Tb respectivamente. Encuentre

d) La densidad de flujo de calor en las partes interna y externa

e) Las temperaturas en las interfases solido-fluido
f) El perfil de temperatura

Solución

El esquema es el siguiente

Este tipo de problemas con convección en uno o ambos extremos para una pared de conductividad
constante se aborda del siguiente modo:

Transferencia de calor por convección en la parte interna:

q0 = ha (Ta − T0 ) (1)

Balance de energía en el material solido

d (r 2 qr )
=0 (2)
dr
Integrando la ecuacion (2)

r 2 qr = cte = r02 q0 = r12 q1 (3)

Es decir, el producto del radio elevado al cuadrado por la densidad de flujo de calor es constante.
Por conveniencia tomamos la siguiente igualdad:

r 2 qr = r02 q0 (4)

De la ecuacion (4) se obtiene

dT 1
−k = r02 q0  2 (5)
dr r

Integrando (5)

T1 r
r02 q0 1 dr
 dT = −
T0
k
 2
r0
r
(6)

Cuyo resultado es

r02 q0  1 1 
T0 − T1 =  −  (7)
k  r0 r1 

Transferencia de calor por convección en la parte externa:

q1 = hb (T1 − Tb ) (8)

Ahora sumamos cada una de las ecuaciones del siguiente modo:

De (1)

q0
Ta − T0 = (9)
ha

De (7)

r02 q0  1 1 
T0 − T1 =  −  (10)
k  r0 r1 

De (8)

q1
T1 − Tb = (11)
hb
Al sumar (9), (10) y (11)

q0 r02 q0  1 1  q1
Ta − Tb = +  − + (12)
ha k  r0 r1  hb

En la ecuacion (12) se debe eliminar el valor de q1 mediante la ecuación (3)

r02
q1 = 2 q0 (13)
r1

Reemplazando (13) en (12)

q0 r02 q0  1 1  r02 q0
Ta − Tb = +  − + 2 (14)
ha k  r0 r1  r1 hb

De la ecuacion (14) se despeja el valor de la densidad de flujo de calor interno q0

r0−2 (Ta − Tb )
q0 = (15)
1 1 1 1 1
+  − +
ha r0 k  r0 r1  hb r12
2

La ecuacion (15) representa la densidad de flujo de calor en la parte interna de la esfera. Para el
cálculo de la densidad de flujo de calor en la superficie externa se hace uso de la ecuacion (13).

Para el cálculo del perfil de temperatura, previamente debe calcularse las temperaturas T0 y T1

La temperatura T0 se obtiene a partir de (9), es decir:

q0
T0 = + Ta (16)
ha

La Temperatura T1 se obtiene a partir de (11), es decir:

q1
T1 = + Tb (17)
hb

Luego al aplicar el balance de energía y resolviendo de manera similar al problema anterior, se

obtiene:

T = T1 − (T1 − T0 ) 
1/ r − 1/ r0  (18)
1/ r1 − 1/ r0 
T − T1  r0 / r − 1 1/ r * − 1
= o T =
*
(19)
T0 − T1  r0 / r1 − 1  r0 / r1 − 1

Problema
Transferencia de calor por conducción en una pared esférica con conductividad variable
En un recipiente esférico semipermeable de radios R1 y R2 se almacena gas licuado a su

temperatura de saturación. Debido a la diferencia de temperatura entre el gas licuado y el

ambiente, se produce flujo de calor desde la parte externa hacia la interna, de tal forma que parte
del gas licuado se pierde por evapotranspiración. Obtenga la expresión respectiva para el calculo
de la tasa de evaporación, considerando que la conductividad varia linealmente con la
temperatura.

Solución

El balance de energía esta dado por

d (r 2 qr )
=0 (1)
dr

Integrando (1)

r 2 qr = C1 (2)
Reemplazando

dT
kr 2 = C1 (3)
dr

Ahora consideramos la conductividad variable

C1
 (k 0 + k1T )dT = 
r2
dr (4)

Integrando la ecuacion (4)

k1 2 C
k0T + T = − 1 + C2 (5)
2 r

Las condiciones de frontera son:

r = R1 T = T1 y r = R2 T = T2

1  k1 2 
C1 =  k (T − T ) + (T1 − T22 )  (6)
1/ R2 − 1/ R1 
0 1 2
2 

k1 2 C1
C2 = k0T1 + T1 + (7)
2 R1

Reemplazando en (5) se tiene

k1 2 1 1
k0 (T − T1 ) + (T − T12 ) −
2 r R1
= (8)
k1 2 1 1
k0 (T1 − T2 ) + (T1 − T2 )
2
−
2 R2 R1

La ecuacion (8) representa la ecuacion del perfil de temperatura. Ahora calculamos la densidad
de flujo de calor.

dT C1
qr r = R1
= −k =− (9)
dr r = R1 R12

Reemplazando el valor de la constante

dT 1 1  k1 2 2 
qr = −k =   k0 (T1 − T2 ) + 2 (T1 − T2 )  (10)
r = R1
dr r = R1 R1 1/ R1 − 1/ R2
2

Ahora se puede calcular el flujo de calor en la superficie interna que también es la misma que el
de la superficie externa, es decir
 4 R12 1  k1 2 2 
Qr = Qr = 4 R12 qr =   k0 (T1 − T2 ) + 2 (T1 − T2 )  (11)
r = R1 r = R2 r = R1 R12
1/ R1 − 1/ R2

El flujo de calor que ingresa al tanque desde la parte externa, evapora una cantidad proporcional
de gas licuado.

Qr r = R1
= m (12)

 : Calor latente de vaporización del gas licuado

m : Flujo de masa

Finalmente, se tiene

Qr r = R1
m= (13)


Problema
Transferencia de calor por conducción en una pared esférica
En una esfera hueca, tal como se muestra en la figura xxxx se transfiere calor desde la parte
exterior hacia el interior. Calcule la temperatura media en la pared esférica, considerando que la
conductividad varia con la temperatura. Las dimensiones, así como las conductividades en los
puntos 1 y 2 se muestran en la sección de datos

Datos

Solución

A partir de la información podemos obtener una ecuacion de la conductividad en función a la

temperatura mediante el uso de proporciones, dado que la dependencia es lineal.
 T − 200 k − 20
= (1)
250 − 200 15 − 20

De esta ecuación se obtiene

k = 40 − 0.1T (2)

Es decir: a = 40 y b = −0.1

Luego usamos la ecuacion (8) del problema anterior

b 1 1
a(T − T0 ) + (T 2 − T02 ) −
2 r R0
= (3)
b 2 1 1
a(T1 − T0 ) + (T1 − T0 )
2
−
2 R1 R0

La temperatura media se calcula por medio de la siguiente ecuacion, para lo cual debe obtenerse
previamente el perfil de temperatura mediante (3)

2  R1

   r sen T (r )drd d
2
R2
3
 T = =  r T (r )dr
0 0 R0 2
2  R2
(4)
R1 − R02
2

   r sen drd d
2 R1

0 0 R0

El numerador de la ecuacion (4) se obtiene por integración numérica, del siguiente modo

h 
x1 n −1

 f ( x)dx = 
2
y0 + 2i =1
yi + yn 

(5)
x0
Usando la formula (5) y los datos de la tercera columna de la tabla xxx se tiene

R1
0.01
 r T (r )dr = (20 + 2(2.598918 + 3.259821 + 3.984026 + 4.772269) + 5.625000)
2

R0
2

= 0.18427534

3
 T = 0.18427534 = 232.76 0C
0.15 − 0.13
3

Problema
Transporte de calor en una pared esférica con transferencia de calor por convección en ambos
lados.
Consideremos una pared esférica de radio interno Ra , radio externo Rb . El coeficiente de

transferencia de calor interno es ha y el externo hb ; asimismo las temperaturas interna y externa

del medio son Ta y Tb respectivamente. Considerando que la conductividad térmica es variable

y está dada por k = a + bT , encuentre:

d) El perfil de temperatura
e) El perfil de la densidad de flujo de calor
f) Las temperaturas en las interfases solido-fluido

Datos

Solución

Este tipo de problemas se resuelve usando un procedimiento numérico que consiste en la solución
simultanea del balance de energía y la ley de Fourier.

Balance de energía en el material solido

d (r 2 qr )
=0 (2)
dr
Ley de Fourier

dT
qr = −(a + bT ) (3)
dr

Transferencia de calor por convección en la parte interna

q0 = ha (Ta − T0 ) (4)

Transferencia de calor por convección en la parte externa

q10 = hb (T10 − Tb ) (5)

Derivando la ecuacion (1)

dqr 2
+ qr = 0 (6)
dr r

Discretizando (6)

2
qi − qi −1 + qi = 0 (8)
i

Expandiendo (8)

2
i = 1 q1 − q0 + q1 = 0 con q0 = ha (Ta − T0 )
1

2
i = 2 q2 − q1 + q2 = 0
2

Y así sucesivamente

2
i = 10 q10 − q9 + q10 = 0
10

Discretizando la ecuacion (3)

Ti − Ti −1
qi + (a + bTi ) =0 (9)
r

Expandiendo la ecuacion (9)

T1 − T0
i = 1 q1 + (a + bT1 ) =0
r
 T2 − T1
i = 2 q2 + (a + bT2 ) =0
r

Y así sucesivamente

T10 − T9 q
i = 10 q10 + (a + bT10 ) = 0 con T10 = 10 + Tb
r hb

El programa respectivo es el siguiente

f(q1) = q1 - q0 + (2 / 1) * q1
f(q2) = q2 - q1 + (2 / 2) * q2
f(q3) = q3 - q2 + (2 / 3) * q3
f(q4) = q4 - q3 + (2 / 4) * q4
f(q5) = q5 - q4 + (2 / 5) * q5
f(q6) = q6 - q5 + (2 / 6) * q6
f(q7) = q7 - q6 + (2 / 7) * q7
f(q8) = q8 - q7 + (2 / 8) * q8
f(q9) = q9 - q8 + (2 / 9) * q9
f(q10) = q10 - q9 + (2 / 10) * q10

f(T0) = q1 + (a + b * T1) * (T1 - T0) / DR

f(T1) = q2 + (a + b * T2) * (T2 - T1) / DR
f(T2) = q3 + (a + b * T3) * (T3 - T2) / DR
f(T3) = q4 + (a + b * T4) * (T4 - T3) / DR
f(T4) = q5 + (a + b * T5) * (T5 - T4) / DR
f(T5) = q6 + (a + b * T6) * (T6 - T5) / DR
f(T6) = q7 + (a + b * T7) * (T7 - T6) / DR
f(T7) = q8 + (a + b * T8) * (T8 - T7) / DR
f(T8) = q9 + (a + b * T9) * (T9 - T8) / DR
f(T9) = q10 + (a + b * T10) * (T10 - T9) / DR
q0 = ha * (Ta - T0)
Ta = 200
Tb = 80
RA = 0.20
RB = 0.40

DR = (RB - RA) / 10
a = 40
b = -0.1
T10 = q10 / hb + Tb
ha = 50
hb = 120

q1(0) = 100
q2(0) = 100
q3(0) = 100
q4(0) = 100
q5(0) = 100
q6(0) = 100
q7(0) = 100
q8(0) = 100
q9(0) = 100
q10(0) = 100

T0(0) = 80
T1(0) = 80
T2(0) = 80
T3(0) = 80
T4(0) = 80
T5(0) = 80
T6(0) = 80
T7(0) = 80
T8(0) = 80
T9(0) = 80

Resultados

q0 5811.076 T0 83.77847
q1 1937.025 T1 82.55808
q2 968.5127 T2 81.94905
q3 581.1076 T3 81.58405
q4 387.4051 T4 81.3409
q5 276.7179 T5 81.16732
q6 207.5384 T6 81.03718
q7 161.4188 T7 80.936
q8 129.135 T8 80.85508
q9 105.6559 T9 80.78888
q10 88.04661 T10 80.73372

Las temperaturas de las interfases interna y externa son 83.77847 y 80.73372 respectivamente

2.5.4 Balance de energía en un cono divergente-convergente

Consideremos un cono divergente-convergente que se muestra en la figura y además, supongamos
que la transferencia de calor se da solo en la dirección axial y que se conocen tanto las dimensiones
geométricas del cono, así como sus respectivas temperaturas en los extremos.

Para realizar el balance de energía, es necesario transformar este sistema cónico a un sistema de
coordenadas rectangulares, en cuyo campo será resuelto dicho balance.
qx A( x) x − qx A( x) x +x =0 (1)

Llevando al limite

qx A( x) x +x − qx A( x)
lim x
=0 (2)
x →0 x

Por definición, la ecuacion (2) se transforma en:

d
(qx A( x)) = 0 (3)
dx

Dado que varia la densidad de flujo de calor como el area de la sección transversal, se aplica la
derivada del producto.

dqx dA( x)
A( x) + qx =0 (4)
dx dx

El área de transferencia de calor puede ser expresada del siguiente modo

A( x) =  y 2 =  (mx + b) 2 (5)

Puesto que y y x están relacionados por la ecuacion de la recta. Por otro lado, en base a la ley
de Fourier se tiene:

dT
qx = − k (6)
dx

Los términos correspondientes de (4) son obtenidos a partir de (5) y (6)

Para el caso de conductividad constante, se tiene:

dqx d 2T
= −k 2 (7)
dx dx

dA( x)
= 2 m(mx + b) (8)
dx

Reemplazando (7) y (8) en (4)

 d 2T   dT 
 (mx + b)   −k 2  +  −k   2 m(mx + b) = 0
2
(9)
 dx   dx 

d 2T 2m dT
+ =0 (10)
dx 2
(mx + b) dx
 dT
Para resolver (10) hacemos el siguiente cambio de variables: = z , luego (10) se transforma
dx
en:

dz 2m
+ z=0 (11)
dx (mx + b)

La ecuacion (11) es fácilmente integrable pues se trata de una ecuacion diferencial de variables
separables.

ln z = −2 ln(mx + b) + ln C1 (12)

A partir de (12 ) se tiene:

dT
= z = C1 (mx + b) −2 (13)
dx

Integrando (13)

C1
T =− + C2 (14)
m(mx + b)

Condiciones de frontera

x = 0 T = T0 y x = L T = TL

T0 − TL
C1 = (15)
1 1
−
m(mL + b) mb

C1
C2 = T0 + (16)
mb

Reemplazando (16) en (14)

C1 C  1 1 
T =− + T0 + 1 = C1  −  + T0 (17)
m(mx + b) mb  mb m(mx + b) 

Ahora se reemplaza (15) en (17)

T0 − TL  1 1  T0 − TL 1 1  (18)
T= − +T = − +T
1 1  mb m(mx + b)  0 1 1  b (mx + b)  0
− −
m(mL + b) mb (mL + b) b
Reordenando

1 1
−
T − T0 (mx + b) b  mL + b  x (19)
= = ( )
TL − T0 1 1  mx + b  L
−
(mL + b) b

La ecuacion (19), representa la ecuacion de distribución de temperatura.

Respecto a la densidad de flujo de calor, se debe tener en cuenta que su valor varia en cada x por
tanto, haciendo uso de la ecuacion (13) se tiene:

dT T0 − TL
qx x =0
= −k = −k (mx + b)−2
dx 1 1
x =0 −
m(mL + b) mb x =0

Evaluando x = 0

(mL + b)
qx x =0
=k (T0 − TL )
Lb

Problema
Transferencia de calor en un cono divergente-convergente de conductividad variable
Obtenga el perfil de temperatura, así como la densidad de flujo de calor y el flujo de calor en un
cono divergente-convergente, considerando los siguientes datos:
T1=250 T2=80 k0=40 k1=-0.1 L=0.5 m R1=0.20 m R2=0.1 m

Solución

La ecuacion de balance de energía es la siguiente

d
(qx A( x)) = 0 (3)
dx

Derivando se obtiene la siguiente ecuacion

dqx dA( x)
A( x) + qx =0 (4)
dx dx

El área de transferencia de calor puede ser expresada del siguiente modo

A( x) =  y 2 =  (mx + b) 2 (5)

dA( x)
= 2 m(mx + b) (6)
dx
Reemplazando (5) y (6) en (4)

dqx
 (mx + b) 2 + qx (2 m(mx + b)) = 0 (7)
dx

Simplificando se tiene

dqx 2mqx
+ =0 (8)
dx (mx + b)

Luego, la ley de Fick para el caso de conductividad variable

dT
qx = −(k0 + k1T ) (9)
dx

Discretizando la ecuacion (8)

qi − qi −1 2mqi
+ =0 (10)
x (m ix + b)

Ahora expandimos la ecuacion (10)

2mxq1
i = 1 q1 − q0 + =0
(m 1x + b)

2mxq2
i = 2 q2 − q1 + =0
(m  2x + b)

Y así sucesivamente

2mxq9
i = 10 q10 − q9 + =0
(m 10x + b)

Ahora discretizando la ecuacion (9)

Ti − Ti −1
qi + (k0 + k1Ti ) =0 (11)
x

Expandiendo la ecuacion (11)

T1 − T0
i = 1 q1 + (k0 + k1T1 ) =0
x

T2 − T1
i = 2 q2 + (k0 + k1T2 ) =0
x
Y así sucesivamente

T10 − T9
i = 10 q10 + (k0 + k1T10 ) =0
x

Las condiciones de frontera son:

x = 0 T0 = 250 q0 = ? y x = L TL = 80 qL = ?

Para la resolución se requiere conocer q0 o q L ; sin embargo, no se conoce. Por tanto, para este

caso se ha desarrollado el programa en el cual se debe suponer q L hasta lograr que cuando x = L
o
la temperatura sea 80 C (Condición dada en el problema)

El programa desarrollado es el siguiente

f(q0) = q1 - q0 + (2 * m * DX * q1) / (m * 1 * DX + b)
f(q1) = q2 - q1 + (2 * m * DX * q2) / (m * 2 * DX + b)
f(q2) = q3 - q2 + (2 * m * DX * q3) / (m * 3 * DX + b)
f(q3) = q4 - q3 + (2 * m * DX * q4) / (m * 4 * DX + b)
f(q4) = q5 - q4 + (2 * m * DX * q5) / (m * 5 * DX + b)
f(q5) = q6 - q5 + (2 * m * DX * q6) / (m * 6 * DX + b)
f(q6) = q7 - q6 + (2 * m * DX * q7) / (m * 7 * DX + b)
f(q7) = q8 - q7 + (2 * m * DX * q8) / (m * 8 * DX + b)
f(q8) = q9 - q8 + (2 * m * DX * q9) / (m * 9 * DX + b)
f(q9) = q10 - q9 + (2 * m * DX * q10) / (m * 10 * DX + b)

f(T1) = q1 + (k0 + k1 * T1) * (T1 - T0) / DX

f(T2) = q2 + (k0 + k1 * T2) * (T2 - T1) / DX
f(T3) = q3 + (k0 + k1 * T3) * (T3 - T2) / DX
f(T4) = q4 + (k0 + k1 * T4) * (T4 - T3) / DX
f(T5) = q5 + (k0 + k1 * T5) * (T5 - T4) / DX
f(T6) = q6 + (k0 + k1 * T6) * (T6 - T5) / DX
f(T7) = q7 + (k0 + k1 * T7) * (T7 - T6) / DX
f(T8) = q8 + (k0 + k1 * T8) * (T8 - T7) / DX
f(T9) = q9 + (k0 + k1 * T9) * (T9 - T8) / DX
f(T10) = q10 + (k0 + k1 * T10) * (T10 - T9) / DX

m = -(R1 - R2) / L
R1 = 0.2
R2 = 0.1
L = 0.5
DX = L / 10
k0 = 40
k1 = -0.1
b = R1
T0 = 250
q10 = 16595 #Valor asumido

q0(0) = 1000
q1(0) = 1000
q2(0) = 1000
q3(0) = 1000
q4(0) = 1000
q5(0) = 1000
q6(0) = 1000
q7(0) = 1000
q8(0) = 1000
q9(0) = 1000

T1(0) = 100
T2(0) = 100
T3(0) = 100
T4(0) = 100
T5(0) = 100
T6(0) = 100
T7(0) = 100
T8(0) = 100
T9(0) = 100
T10(0) = 100

Resultados

T0
q0
250
3493.684
q1 T1
3904.706 237.952
4392.794 225.3743
4978.5 212.1248
5689.714 198.0386
6565.055 182.9175
7659.231 166.5155
9051.818 148.5185
1.09E+04 128.5135
1.33E+04 105.9399
1.66E+04 80.00944

20000
Título del gráfico

15000

10000

5000

0
0 2 4 6 8 10 12
 Título del gráfico
300

250

200

150

100

50

0
0 2 4 6 8 10 12

2.6 Transporte de calor en paredes compuestas

2.2.1 Pardes compuestas planas

Consideremos una pared plana compuesto de tres materiales tal como ocurre en la industria. El
primer material es un ladrillo refractario, el segundo material es ladrillo aislante y el ultimo
material es generalmente de hierro o acero ademas de aislar ligeramente, se utiliza para conformar
la estructura respectiva, tal como se muestra en la figura.

Consideremos que la conductividad de cada uno de los materiales es constante. En este caso, se
trata de una extensión del problema zzzz

Transferencia de calor por convección en el lado izquierdo de la pared

q0 = ha (Ta − T0 ) (1)

Transferencia de calor por convección en el lado derecho de la pared

q3 = hb (T3 − Tb ) (2)
Ecuacion de transporte de calor en el material solido

dqx
=0 (3)
dx

Integrando (3)

qx = cte = q0 = q1 = q2 = q3 (4)

Aplicando la Ley de Fourier para la pared (I)

dT 
−k 01 = q0 (5)
dx

Integrando la ecuacion (5)

T1 x1
q0
 dT  = −  k
T0 x0
01
dx (6)

Cuyo valor integrado es:

q0
T0 − T1 = ( x1 − x0 ) (7)
k 01

Del mismo modo para la pared (II)

dT 
−k 12 = q0 (8)
dx

Integrando (8)

T2 x2
q0
 dT  = −  k
T1 x1
12
dx (9)

q0
T1 − T2 = ( x2 − x1 ) (10)
k 12

Del mismo modo para la pared (III)

dT 
−k 23 = q0 (11)
dx

Integrando (8)

T3 x3
q0
 dT  = − 
T2 x2
k 23
dx (12)
 q0
T2 − T3 = ( x3 − x2 ) (13)
k 23

Sumando las ecuaciones (1), (2), (7), (10) y (13)

De (1)

q0
Ta − T0 =
ha

De (7)

q0
T0 − T1 = ( x1 − x0 ) (7)
k 01

De (10)

q0
T1 − T2 = ( x2 − x1 ) (10)
k 12

De (13)

q0
T2 − T3 = ( x3 − x2 ) (13)
k 23

De (2)

q3 q0
T3 − Tb = = (14)
hb hb

Cuya suma es

q0 q0 q q q
Ta − Tb = + 01 ( x1 − x0 ) + 120 ( x2 − x1 ) + 230 ( x3 − x2 ) + 0 (15)
ha k k k hb

Luego a partir de (15) se obtiene

Ta − Tb
q0 = (16)
1 ( x1 − x0 ) ( x2 − x1 ) ( x3 − x2 ) 1
+ + + +
ha k 01 k 12 k 23 hb

La ecuacion (16) permite calcular la densidad de flujo ce calor en la pared plana, cuyo valor es
constante.

El flujo de calor se obtiene multiplicando la densidad de flujo de calor por el área de transferencia
de calor, es decir:
 WL(Ta − Tb )
Q = q0WL = (17)
1 ( x1 − x0 ) ( x2 − x1 ) ( x3 − x2 ) 1
+ + + +
ha k 01 k 12 k 23 hb

La ecuacion (17) se puede expresar del siguiente modo:

WL(Ta − Tb ) WL(Ta − Tb )
Q= = (18)
1 ( x1 − x0 ) ( x2 − x1 ) ( x3 − x2 ) 1
+ 01
+ 12
+ 23
+ R
ha k k k hb

 R : Es la suma de resistencias a la transferencia de calor: Dos resistencias a la convección y

tres resistencias a la convección.

WL(Ta − Tb )
Q= = ( R)−1WL(Ta − Tb )  = UWL(Ta − Tb ) (19)
R
U : Coeficiente global de transferencia d calor (incluye a la conducción y convección) y tiene las
mismas unidades que el coeficiente de transferencia de calor por convección, y como puede verse
es la inversa de la suma de las resistencias.

Ahora estos conceptos pueden ser ampliados a estructuras más complejas

El diagrama de resistencias es el siguiente:

En base a lo desarrollado anteriormente, se tiene:

1
R1 = (1)
ha

Para las resistencias 2 y 3 se deben sumar las inversas

1 1 1
= 01 + 01 (2)
R2 R2 R3

Reemplazando se tiene

1 1 1
= + (3)
R2 ( x1 − x0 ) ( x1 − x0 )
k201 k301

( x1 − x0 )
R2 = (4)
k201 + k301

Del mismo modo para las resistencias en paralelo

( x2 − x1 )
R3 = (5)
k + k512 + k612
12
4

x3 − x2
R4 = (6)
k4

1
R5 = (7)
hb

Sumando todas las resistencias se tiene:

(Ta − Tb ) (Ta − Tb )
q0 = = (8)
 R 1 + ( x011 − x001) + 12( x2 −12 x1 ) 12 + x3 − x2 + 1
ha k2 + k3 k 4 + k5 + k 6 k4 hb

La misma que puede ser expresada como

q0 = ( R) −1 (Ta − Tb ) = U (Ta − Tb ) (9)

−1
 1 (x − x ) (x − x ) x −x 1
U =  + 011 001 + 12 2 12 1 12 + 3 2 +  (10)
 ha k2 + k3 k4 + k5 + k6 k4 hb 
En base a la figura, las resistencias en paralelo deben ser adicionadas de forma adecuada para
tener todas las resistencias en serie.

2.2.1 Pardes cilíndricas compuestas

Consideremos una pared cilíndrica compuesta y que la conductividad de cada uno de los
materiales es constante. Para obtener la densidad de flujo de calor, así como el perfil de
temperatura, procederemos de la misma forma que en el caso del problema de paredes planas
Transferencia de calor por convección entre el fluido y la superficie interna del cilindro

q0 = ha (Ta − T0 ) (1)

Balance de energía en el material solido (I)

d (rqr )
=0 (2)
dr

Balance de energía en el material solido (II)

d (rqr)
=0 (2)
dr

Balance de energía en el material solido (III)

d (rqr)
=0 (3)
dr

Integrando las ecuaciones (1), (2) y (3)

rqr = rqr = rqr = r0 q0 = r1q1 = r2 q2 = r3q3 = cte (4)

Para el material (I) se tiene

rqr = r0 q0 (5)

Que es lo mismo que

dT  r0 q0
−k 01 = (6)
dr r

T1 r
r0 q0 1 1
 dT  = −
T0
k 01 r0 r
 dr (7)

Integrando se obtiene

r0 q0
T0 − T1 =  ln(r1 / r0 ) (7)
k 01

Del mismo modo para el material (II)

rqr = r0 q0 (8)

Reemplazando a ley de Fourier e integrando

r0 q0
T1 − T2 =  ln(r2 / r1 ) (9)
k 12

Y para el material (III)

rqr = r0 q0 (10)

Reemplazando a ley de Fourier e integrando

r0 q0
T2 − T3 =  ln(r3 / r2 ) (11)
k 23

Ahora la transferencia de calor por convección en la parte externa

q3 = hb (T3 − Tb ) (12)

Para agrupar adecuadamente, las ecuaciones (1), (7), (9), (11) y (12) deben estar en términos de
q0 ; sin embargo, la ecuacion (12) no lo está, por lo tanto, es necesario utilizar la siguiente
igualdad dada anteriormente.

r0 q0 = r3q3 (13)
 r0
q3 = q0 (14)
r3

Ahora colocamos todos los términos en forma secuencial para sumarlos y eliminar los términos
intermedios, del siguiente modo:

q0
Ta − T0 = (15)
ha

r0 q0
T0 − T1 =  ln(r1 / r0 ) (16)
k 01

r0 q0
T1 − T2 =  ln(r2 / r1 ) (17)
k 12

r0 q0
T2 − T3 =  ln(r3 / r2 ) (18)
k 23

r0 q0
T3 − Tb = (19)
r3hb

Al sumar se tiene

1 ln(r1 / r0 ) ln(r2 / r1 ) ln(r3 / r2 ) 1

Ta − Tb = r0 q0 ( + 01
+ 12
+ 23
+ ) (20)
r0 ha k k k r3hb

1 ln(r1 / r0 ) ln(r2 / r1 ) ln(r3 / r2 ) 1 −1

q0 = r0−1 ( + 01
+ 12
+ 23
+ ) (Ta − Tb ) (21)
r0 ha k k k r3hb

Dado que todos los cálculos se realizaron en función a la densidad de flujo de calor en la superficie
interna, q0 , entonces, el coeficiente global de transferencia U 0 , se denomina coeficiente global

de transferencia de calor basa do en el radio interno

q0 = U 0 (Ta − Tb ) (22)

1 ln(r1 / r0 ) ln(r2 / r1 ) ln(r3 / r2 ) 1 −1

U 0 = r0−1 ( + 01
+ 12
+ 23
+ ) (23)
r0 ha k k k r3hb

Q0 = 2 r0 Lq0 = 2 r0 LU 0 (Ta − Tb ) (24)

Del mismo modo, si tomamos como referencia el radio externo, se obtiene la siguiente ecuacion
de transferencia de calor
 1 ln(r1 / r0 ) ln(r2 / r1 ) ln(r3 / r2 ) 1 −1
q3 = r3−1 ( + 01
+ 12
+ 23
+ ) (Ta − Tb )
r0 ha k k k r3hb

q3 = U 3 (Ta − Tb )

1 ln(r1 / r0 ) ln(r2 / r1 ) ln(r3 / r2 ) 1 −1

U 3 = r3−1 ( + 01
+ 12
+ 23
+ )
r0 ha k k k r3hb

Q3 = 2 r3 Lq3 = 2 r3 LU 3 (Ta − Tb )

Conclusión: Las densidades de flujo de calor interna y externa q1 y q3 no son iguales, debdo a

que el área de transferencia de calor también lo es; sin embargo, el flujo de calor es el mismo (

Q0 = Q3 ), puesto que e fundamento del balance de energía esta soportada en el principio de que
el flujo de calor es constante.

bbbbbbbbbbbbbbbbbbbbbb

2.2.3 Paredes esféricas compuestas

Consideremos una pared esférica compuesta y que la conductividad de cada uno de los materiales
es constante. Para obtener la densidad de flujo de calor, así como el perfil de temperatura,
procederemos de la misma forma que en el caso del problema de paredes planas

Transferencia de calor por convección entre el fluido y la superficie interna del cilindro

q0 = ha (Ta − T0 ) (1)

Balance de energía en el material solido (I)

d (r 2 qr )
=0 (2)
dr
Balance de energía en el material solido (II)

d (r 2 qr)
=0 (2)
dr

Balance de energía en el material solido (III)

d (r 2 qr)
=0 (3)
dr

Integrando las ecuaciones (1), (2) y (3)

r 2 qr = r 2 qr = r 2 qr = r02 q0 = r12 q1 = r22 q2 = r32 q3 = cte (4)

Para el material (I) se tiene

r 2 qr = r02 q0 (5)

Que es lo mismo que

dT  r02 q0
−k 01
= 2 (6)
dr r
T1 r
r02 q0 1 1

T0
dT  = −
k 01 r0 r 2
 dr (7)

Integrando se obtiene

r02 q0  1 1 
T0 − T1 =  −  (7)
k 01  r0 r1 

Del mismo modo para el material (II)

r 2 qr = r02 q0 (8)

Reemplazando a ley de Fourier e integrando

r02 q0  1 1 
T1 − T2 =  −  (9)
k 12  r1 r2 

Y para el material (III)

r 2 qr = r02 q0 (10)

Reemplazando a ley de Fourier e integrando

r02 q0  1 1 
T2 − T3 =  −  (11)
k 23  r2 r3 

Ahora la transferencia de calor por convección en la parte externa

q3 = hb (T3 − Tb ) (12)

Para agrupar adecuadamente, las ecuaciones (1), (7), (9), (11) y (12) deben estar en términos de

q0 ; sin embargo, la ecuacion (12) no lo está, por lo tanto, es necesario utilizar la siguiente

igualdad dada anteriormente.

r02 q0 = r32 q3 (13)

r02
q3 = q0 (14)
r32

Ahora colocamos todos los términos en forma secuencial para sumarlos y eliminar los términos
intermedios, del siguiente modo:

q0
Ta − T0 = (15)
ha

r0 q0  1 1 
T0 − T1 =  −  (16)
k 01  r0 r1 

r02 q0  1 1 
T1 − T2 =  −  (17)
k 12  r1 r2 

r02 q0  1 1 
T2 − T3 =  −  (18)
k 23  r2 r3 

r02 q0
T3 − Tb = (19)
r32 hb

Al sumar se tiene
 1 1  1 1 1 1 1  1  1 1  1
Ta − Tb = r02 q0 ( 2
+ 01   −  + 12   −  + 23   −  + 2 )
r h k
0 a  r0 r1  k  r1 r2  k  r2 r3  r3 hb
(20)

−1
 1 1 1 1 1 1 1 l 1 1 1 
q0 = r0−2  2 + 01  −  + 12  −  + 23  −  + 2  (Ta − Tb ) (21)
 r0 ha k  r0 r1  k  r1 r2  k  r2 r3  r3 hb 

Dado que todos los cálculos se realizaron en función a la densidad de flujo de calor en la superficie

interna, q0 , entonces, el coeficiente global de transferencia U 0 , se denomina coeficiente global

de transferencia de calor basa do en el radio interno

q0 = U 0 (Ta − Tb ) (22)

−1
 1
−2 1 1 1 1 1 1 l 1 1 1 
U 0 = r  2 + 01  −  + 12  −  + 23  −  + 2 
0 (23)
 r0 ha k  r0 r1  k  r1 r2  k  r2 r3  r3 hb 

Q0 = 4 r02 q0 = 4 r02U 0 (Ta − Tb ) (24)

Del mismo modo, si tomamos como referencia el radio externo, se obtiene la siguiente ecuacion
de transferencia de calor

−1
 1
−2 1 1 1 1 1 1 l 1 1 1 
q3 = r  2 + 01  −  + 12  −  + 23  −  + 2  (Ta − Tb )
3 (25)
 r0 ha k  r0 r1  k  r1 r2  k  r2 r3  r3 hb 

q3 = U 3 (Ta − Tb ) (26)

−1
 1
−2 1 1 1 1 1 1 l 1 1 1 
U 3 = r  2 + 01  −  + 12  −  + 23  −  + 2 
3 (27)
 r0 ha k  r0 r1  k  r1 r2  k  r2 r3  r3 hb 

Q3 = 4 r32 q3 = 4 r32U 3 (Ta − Tb ) (28)

Conclusión: Las densidades de flujo de calor interna y externa q1 y q3 no son iguales, debdo a

que el área de transferencia de calor también lo es; sin embargo, el flujo de calor es el mismo (

Q0 = Q3 ), puesto que el fundamento del balance de energía esta soportada en el principio de que
el flujo de calor es constante.
Problema
Espesor óptimo de aislamiento
Se desea evitar pérdidas de calor en una tubería que conduce vapor de agua a la temperatura de
500 oF, para lo cual se aísla dicha tubería con un material aislante cuya conductividad térmica es
25 BTU/h pie oF; la temperatura del ambiente es 77 oF. El radio interno de la tubería es 0.99 pies
y el externo 1 pie; puede estimarse que el coeficiente de transferencia de calor en la parte interna
es de 200 BTU/hr.pie2 oF y el externo 24 BTU/hr.pie2. oF; la conductividad termica de la tubería
es 50 BTU/[Link] oF

Encuentre

a) El espesor óptimo de aislamiento que minimice la perdida de calor

b) El flujo de calor por unidad de longitud de tubería

Solución

Usaremos la ecuacion desarrollada anteriormente

Q0 = 2 r0 LU 0 (Ta − Tb ) (1)

A partir de (1) se tiene

Q0 / L = 2 r0U 0 (Ta − Tb ) (2)

1 ln(r1 / r0 ) ln(r2 / r1 ) 1 −1
U 0 = r0−1 ( + + + ) (3)
r0 ha k 01 k 12 r2 hb

1 ln(r1 / r0 ) ln(r2 / r1 ) 1 −1
Q0 / L = 2 r0 r0−1 ( + + + ) (Ta − Tb ) (4)
r0 ha k 01 k 12 r2 hb

Simplificando la ecuacion (4)

(Ta − Tb )
Q0 / L = 2 (5)
1 ln(r1 / r0 ) ln(r2 / r1 ) 1
( + + + )
r0 ha k 01 k 12 r2 hb

Por el criterio de optimo

d (Q0 / L)
=0 (6)
dr2

Derivando (5) se tiene:

d (Q0 / L) 1 ln(r1 / r0 ) ln(r2 / r1 ) 1 −1  1 r1 1 1 

= −2 (Ta − Tb )( + 01
+ 12
+ )   −   = 0 (7)
dr2 r0 ha k k r2 hb  k 12 r2 h2 r22 

Simplificando la ecuacion (7)

 1 r1 1 1 
 12  −  2 =0 (8)
k r2 h2 r2 

Despejamos el valor del radio critico

k 12
r2 = (9)
h2

Reemplazando datos

25
r2 = = 1.0416 pies (10)
24

Luego, el espesor de aislamiento es:

 = r2 − r1 = 1.0416 − 1 = 0.0416 pies = 0.5 pu lg

Reemplazando datos en (5)

2 (500 − 77) BTU

Q0 / L = = 56690 (11)
1 ln(1/ 0.99) ln(1.0416 / 1) 1 hr pie
( + + + )
0.99  200 50 25 1.0416  24

Problema
Transferencia de calor en un secador rotatorio horizontal
En el secador experimental que se muestra en la figura xxxx pasa aire caliente a 60 oC; el
material del tubo es de vidrio pyrex de 2mm de espesor y su diámetro interno es de 50 mm; el
coeficiente de transferencia de calor interno puede estimase en 18 w/m2 oC y el externo 25
w/m2 oC; la conductividad térmica del vidrio es 0.75W/m oC. Considerando que la temperatura
del medio es 20oC, encuentre la tasa de perdida de calor por unidad de longitud
Solución

Se plantea el siguiente esquema para estimar las pérdidas de calor en el secador de vidrio

A partir de la ecuación (21)

1 ln(r1 / r0 ) 1 −1
q0 = r0−1 ( + ) (Ta − Tb ) (1)
r0 ha k r3hb

Reemplazando datos

1 ln(0.052 / 0.05) 1
q0 = 0.05−1 ( + ) −1 (60 − 20) = 376W / m2 (2)
0.05 18 0.75 0.052  20

El flujo de calor se puede expresar mediante la siguiente ecuacion

Q0 = 2 r0 Lq0 (3)
A partir de la ecuacion (3) se tiene

Q0 / L = 2 r0 q0 = 2  0.05  376.48 = 118.7W / m (4)

2.3. Superficie de enfriamiento

Las superficies de enfriamiento o extensiones o aletas son dispositivos que permiten la

transferencia de calor por conducción dentro de sus límites, así como la transferencia de calor por
convección (o radiación) entre sus límites y sus alrededores. La aplicación más frecuente de las
superficies extendidas es la de incrementar la rapidez de transferencia de calor desde un material
que se calienta por diversas causas (resistencia eléctrica u otro medio), mediante un fluido (aire)
que se encuentra en su cercanía.

En base al diseño, las aletas se pueden clasificar en aletas de área de sección constante y aletas de
área de sección variable.
2.3.1 Balance de energía en una aleta

La realización del balance de energía en una aleta es de suma importancia pues nos permite
relacionar las dimensiones geométricas de la aleta con la tasa de disipación de calor, obteniendo
principalmente el perfil de temperatura en función a las condiciones de frontera que a su vez
depende de las dimensiones de la aleta y el área transversal correspondiente. En forma general, el
balance de energía para cualquier tipo de aleta está dada por:

Velocidad de entrada de  Velocidad de salida de  Velocidad de salida de 

energia por conduccion  − energia por conduccion  − energia por conveccion  = 0 (1)
     

2.3.2 Factor de efectividad

El parámetro más importante a tener en cuenta en una aleta es el factor de efectividad o la

eficiencia de la aleta que permite conocer la tasa de disipación de energía calorífica, cuyo
conocimiento es muy útil en su diseño.

calor disipado por la sup erficie de la aleta

= (2)
calor que se disiparia si la sup erficie se mantuviese a Tw sin var iar h

2.3. Aleta de rectangular

Esta es un tipo de aleta de fácil diseño y construcción motivo por el cual tiene gran aplicación en
el campo industrial. El cuerpo caliente puede tener cualquier forma; sin embargo, la superficie de
enfriamiento es una lámina de longitud L , ancho W y espesor. En este caso el área de la sección
transversal es constante e igual a 2BW (área de transferencia de calor por conducción) , cuya
conductividad se puede considerar constante o variable y 2LW (área de transferencia de calor
por convección: superficie superior y superficie inferior), cuya tasa de transferencia de calor
depende del coeficiente de transferencia de calor por convección.
Asimismo, la temperatura de la superficie caliente (cuerpo del cual se quiere disipar el calor), está

a la temperatura TW constante y la temperatura del medio o temperatura del ambiente igual a Ta

Hagamos el balance de energía para este caso, considerando que hay entrada ÑÑy salida de calor
por conducción por el área normal y una salida adicional por convección por el área lateral del
elemento diferencial de volumen ( 2BW Z ). En base a la ecuacion (1) se tiene:

qz z 2 BW − qz z +z
2 BW − 2hW z (T − Ta ) = 0 (3)

Llevando al limite

qz 2 BW − qz z 2 BW 2hW z (T − Ta )
lim ( z +z
)+ =0 (4)
z →0 2 BW z 2 BW z

dqz h
+ (T − Ta ) = 0 (5)
dz B

La ecuacion (5) es la ecuacion diferencial que describe el comportamiento de la aleta rectangular

en el proceso de disipación de calor.

La resolución para el caso de conductividad térmica constante es de tipo analítico, lo cual

ilustramos a continuación.

Ley de Fourier

dT
qz = − k (6)
dz

Reemplazando (6) en (5)

d 2T h
− (T − Ta ) = 0 (7)
dz 2 kB
Las condiciones de frontera para este caso son:

dT
Aleta de longitud finita: z = 0 T = TW z = L −k = h(T − Ta )
dz

dT
Aleta de longitud semi infinita: z = 0 T = TW z=L =0
dz

Aleta de longitud infinita: z = 0 T = TW z = L T = Ta

Es decir para la condición final hay tres posibilidades, dependiendo de la longitud de la aleta.

METODOS ANALITICOS

Problema
Aleta rectangular 1
Obtenga el perfil de temperatura y el factor de efectividad de la aleta rectangular (longitud finita)
cuyas condiciones de frontera son:
dT
z = 0 T = TW z = L −k = h(T − Ta )
dz

Utilice los siguientes datos:

L = 0.1 m , h = 20 w / m 2 oC , k = 40 w / m 2 oC , TW = 200 oC , Ta = 20 oC , B = 0.02 m

Solución

La ecuacion diferencial de la aleta rectangular es:

d 2T h
− (T − Ta ) = 0 (1)
dz 2 kB

Hacemos el siguiente cambio de variable: u = T − Ta en consecuencia: du = dT . De estos

cambios se obtiene:

d 2u d 2T
= (2)
dz 2 dz 2

Reemplazando (2) en (1)

d 2u h
− u=0 (3)
dz 2 kB

La ecuacion algebraica característica de (3) que se resuelve mediante el método de operadores es:
 h
(m 2 − )u = 0 (4)
kB

A partir de (4) se tiene

1/2
 h 
m =   (5)
 kB 

La solución general de (5) es:

u = C1e mz + C2e − mz (6)

Cuando z = 0 , se tiene

(TW − Ta ) = C1 + C2 (7)

Derivando la ecuacion (6)

du
= C1memz − C2 me− mz (8)
dz

La ecuacion (8) puede ser reescrita del siguiente modo

dT
= C1memz − C2 me− mz (9)
dz

Luego multiplicando (9) por − k se tiene:

dT
−k = −kC1memz + kC2 me− mz (10)
dz

Igualando la ecuacion (10) con la condición de frontera

−kC1me mL + kC2 me − mL = h(TL − Ta ) (11)

A partir de (11) se obtiene la temperatura en el extremo de la aleta TL

C1 mL C
TL − Ta = −k me + k 2 me − mL (12)
h h

Ahora que se conoce TL , se puede igualar con (6) cuando z = L

C1 mL C
−k me + k 2 me− mL = C1emL + C2e − mL (13)
h h
La ecuación (13) se puede reescribir del siguiente modo

−C1 Ae mL + C2 Ae − mL = C1e mL + C2e − mL (14)

km
Siendo A = una constante adimensional, y de (14) se obtiene
h

 Ae− mL − e− mL  ( A − 1) e− mL
C1 = C2  mL 
= C  (15)
 Ae + e  ( A + 1) emL
mL 2

Ahora se resuelve simultáneamente (15) y (7). De (7) se tiene

C1 = (TW − Ta ) − C2 (16)

Igualando (15) y (16)

( A − 1) e− mL
C2  = (TW − Ta ) − C2 (17)
( A + 1) emL

Despejando C 2

(TW − Ta )
C2 = (18)
( A − 1) e − mL
1+ 
( A + 1) e mL

Reemplazando en (17)

 
 1 
C1 = 1 −  (TW − Ta ) (19)
 ( A − 1) e − mL 
 1 +  
 ( A + 1) e mL 

Reemplazando en (6) modificado

T = C1e mz + C2 e − mz + Ta (20)

El factor de efectividad se obtiene mediante:

LW L L

  h(T − T )dydz
a W  (T − Ta )dz  (T − T )dz
a

= 0 0
= 0
= 0
(21)
LW
LW (TW − Ta ) L(TW − Ta )
  h(T
0 0
W − Ta )dydz
Ahora se reemplaza la ecuacion (20) en (21)

 (C e + C2e − mz )dz
mz
1

= 0
(22)
L(TW − Ta )

L
 C1 mz C2 − mz  C1 mL C
 m e − m e  (e − 1) − 2 (e − mL − 1)
= 0
= m m (23)
L(TW − Ta ) L(TW − Ta )

De esta manera queda determinado el factor de efectividad, dado que las constantes C1 y C 2 son

conocidas.

d(Y)/d(Z) = 1
Y(0) = 0
C1 = (Tw - Ta) * (1 - 1 / E)
C2 = (Tw - Ta) / E
E = 1 + ((A - 1) / (A + 1)) * exp(-m * L) / exp(m * L)
T = C1 * exp(m * Z) + C2 * exp(-m * Z) + Ta
FE = ((C1 / m) * (exp(m * L) - 1) - (C2 / m) * (exp(-m * L) - 1)) / (L * (Tw - Ta))
L = 0.1
h = 20
k = 40
Tw = 200
Ta = 20
m = (h/(k * B)) ^ 0.5
B = 0.02
A=k*m/h
Z(0) = 0
Z(f) = 0.1

El factor de efectividad analítico sin considerar la transferencia de calor por convección por la
parte lateral cuya área es 2( L  2 B) , por ser pequeña respecto al área 2( L  W ) , se obtiene

 = 0.9050
Problema
Aleta rectangular solución analítica
Demuestre que la ecuacion de distribución de temperatura, con transferencia de calor por
convección en la sección final, cuyas condiciones están dadas por:
dT
z = 0 T = TW z = L −k = h(T − Ta )
dz

Es igual a:

T − Ta cosh m( L − x) + (h / mk )sinh m( L − x)
=
TW − Ta cosh(mL) + (h / mk )sinh(mL)

Utilice los siguientes datos:

L = 0.1 m , h = 20 w / m 2 oC , k = 40 w / m 2 oC , TW = 200 oC , Ta = 20 oC , B = 0.02 m

Solución

La ecuacion diferencial de la aleta rectangular es:

d 2T h
− (T − Ta ) = 0 (1)
dz 2 kB

Hacemos el siguiente cambio de variable: u = T − Ta en consecuencia: du = dT . De estos

cambios se obtiene:

d 2u d 2T
= (2)
dz 2 dz 2

Reemplazando (2) en (1)

d 2u h
− u=0 (3)
dz 2 kB

La ecuacion algebraica característica de (3) que se resuelve mediante el método de operadores es:

h
(m 2 − )u = 0 (4)
kB

A partir de (4) se tiene

1/2
 h 
m =   (5)
 kB 
La solución general de (5) es:

u = C1e mz + C2e − mz (6)

Cuando z = 0 , se tiene

(TW − Ta ) = C1 + C2 (7)

Cuando z = L

TL − Ta = C1e mL + C2e − mL (8)

Derivando la ecuacion (6)

du
= C1memz − C2 me− mz (9)
dz

La ecuacion (8) puede ser reescrita del siguiente modo

dT
= C1memz − C2 me− mz (10)
dz

Luego multiplicando (9) por − k se tiene:

dT
−k = −kC1memz + kC2 me− mz (11)
dz

Cuando z = L , la transferencia de calor por conducción es igual a la transferencia de calor por

convección.

h(TL − Ta ) = −kC1me mL + kC2 me − mL (12)

A partir de (12) se obtiene la temperatura en el extremo de la aleta (TL − Ta )

km km
(TL − Ta ) = − C1e mL + C2e − mL (13)
h h

Ahora se tienen tres ecuaciones: (7), (8) y (13) con las variables C1 , C 2 y TL . Igualando (8) y

(13)

km km
C1emL + C2e − mL + C1emL − C2e − mL = 0 (14)
h h

km km
(1 + )C1emL + (1 − )C2e− mL = 0 (15)
h h
Por otro lado, a partir de la ecuacion (7)

C1 = (TW − Ta ) − C2 (16)

Reemplazando (16) en (15)

km km
(1 + )((TW − Ta ) − C2 )e mL + (1 − )C2e − mL = 0 (17)
h h

km mL km km
−(TW − Ta )(1 + )e + (1 + )C2emL − (1 − )C2e − mL = 0 (18)
h h h

Despejando el valor de C 2

km mL
(TW − Ta )(1 +
)e
h (TW − Ta )(h + km)e mL
C2 = = − mL
(19)
(1 + )emL − (1 − )e− mL (h + km)e − (h − km)e
mL
km km
h h

− mL
Dividiendo el numerador y denominador por e

(TW − Ta )(h + km)e 2 mL

C2 = (20)
(h + km)e2 mL − (h − km)

A partir de las ecuaciones (16) y (20) se obtiene el valor de C1

(TW − Ta )(h + km)e2 mL

C1 = (TW − Ta ) − (21)
(h + km)e2 mL − (h − km)

(TW − Ta )(h + km)e2 mL − (TW − Ta )(h − km) − (TW − Ta )(h + km)e 2 mL

C1 =
(h + km)e2 mL − (h − km)

Simplificando

(TW − Ta )(h − km)

C1 = − (23)
(h + km)e2 mL − (h − km)

Reemplazando C1 y C 2 en (6)

(TW − Ta )(h − km) (TW − Ta )(h + km)e 2 mL − mz

T − Ta = − e mz
+ e (24)
(h + km)e2 mL − (h − km) (h + km)e 2 mL − (h − km)

(TW − Ta )(h + km)e 2 mL − mz (TW − Ta )(h − km)

T − Ta = e − emz
(h + km)e − (h − km)
2 mL
(h + km)e − (h − km)
2 mL