2.

CONDUCCIÓN
UNIDIMENSIONAL EN ESTADO
ESTABLE.
 a) LA PARED PLANA

En flujo estable con fuente no distribuida

de energía.

Fluido Fluido
Ts 2  Ts1
Caliente Ts1 frío C2  Ts1 y C1  entonces
Ts2
L
T∞1,h1 T∞2,h2 Ts 2  Ts1
T ( x)  x  Ts1
L
 dT kA
d  dT  dT Q x  kA  (Ts1  Ts 2 )
 k   0, si k  Cte ;  C1 dx L
dx  dx  dx
 k
 T ( x)  C1 x  C2 Q"x  (Ts1  Ts 2 )
L
y con T (0)  Ts1 y T ( L)  Ts 2
 RESISTENCIA TÉRMICA
Haciendo una analogía con el sistema eléctrico:
Re → Resistencia eléctrica; V → Voltaje; I → Intensidad de corriente eléctrica
Rt → Resistencia térmica; T → Temperatura; Q → Calor.

V1  V2 L Ts1  Ts 2 L
Re   ; Rtcond  
I A 
Qx kA
Ts  T 1
En convección : Rtconv  

Q hA
En el circuito anterior
 T1  Ts1 Ts1  Ts 2 Ts 2  T 2
Qx   
1 L 1
h1 A kA h2 A
 CONTINUA RESISTENCIA TÉRMICA

En tér min os de (T1  T 2 )

 (T1  T 2 )
Qx 
Rtot

Rtot  1 L  1
h1 A kA h2 A

Que representa la resistencia de un circuito de resistencias en serie.

T∞1 Ts1 Ts2 T∞2 T∞1 T∞2



Qx 1/h1A L/kA 1/h2A Qx Rtot
 PARED COMPUESTA

Para una pared compuesta.

Si; U = Coeficiente global de transferencia de
calor, se define:

A B C 
Ts1 T2 Q x  UAT

T3 Ts2
1
U
T∞1, h1 T∞4;h4 Rt A
fluido fluido 1
caliente x frío 
1 h1  LA k A  LB k B  LC kC  1 h4
 T1  T 4
Qx  ; donde T 1
 Rt Rtot   Rt  

1 L L L 1 UA
 Rt   A  B  C  Qx
h1 A k A A k B A kC A h4 A
 T T T T
Q x  1 s1  s1 2    
1 h1 A LA k A A
 CIRCUITO TÉRMICO EN PARALELO

Una pared compuesta como se muestra Se puede representar como:

A B D RA Req RD
T1 T2
Donde.
C
RB RC
Req 
El circuito térmico es RB  RC
RB
RA RD Y también:
T1 T2
Rtot  RA  Req  RD
RC
 RESISTENCIA DE CONTACTO
Cuando se tienen dos superficies en contacto, debido a sus irregularidades, se presentan
secciones en donde se tienen caídas de temperatura entre estas dos superficies y por lo
tanto, una resistencia térmica llamada resistencia de contacto ( R”tc). El valor de esta
resistencia depende de la presión a que están unidas esta dos superficies, su material y del
tipo de fluido entre estas irregularidades.

A B
TA  TB
R"tc  
Q"x TA

Q" x TB
R”tc depende de:

* Acabado superficial
* Presión de contacto. x
* Fluido entre irregularidades

RA R”tc RB
Ejemplo 2.1. El vidrio de ventana trasera de un carro de vidrio de 4 mm espesor, es
desempañada por una resistencia eléctrica en su superficie interna. Determine la potencia
eléctrica por unidad de área de la ventana para mantener una temperatura de 15 0C. La
temperatura interior es de 25 0C y hi = 10 w/m2k, al exterior -10 0C y he = 65 w/m2k

SE CONOCE: Temperatura deseada vidrio y ANÁLISIS. El circuito Térmico:

condiciones interior y exterior de un carro. T∞i Tsi T∞e
SE BUSCA: Potencia por unidad de área para

mantener esa temperatura deseada. 1/hiA L/kA 1/heA
SE ASUME:,Flujo unidimensional, estado estable
Q " 

Propiedades constantes, radiación y resistencia Q"e

de película despreciables.
ESQUEMA. Ti  Tsi  Tsi  Te
 Q"e 
1 L  1
Aire interior Aire del ambiente hi k he
 T T T T
Td Q"e  si e  i si  1270w / m 2
vidrio L  1 1
k he hi
T∞i  T∞ 
hi Qe" he Nota : con Q"e  0; Tsi  4.60 C

Propiedades: Vidrio a 300 0K; k = 1.4 w/m 0K

 b) SISTEMAS RADIALES
Un problema común es tener un cilindro hueco Con las condiciones de que:

T r1   Ts1 y T r2   Ts 2

cuyas superficies interior y exterior están a
fluidos de diferentes temperaturas.

L d2 Ts1  C1 ln r1  C2 y Ts 2  C1 ln r2  C2
d1
Ts1  Ts 2  r 
T r   ln    Ts 2
ln  1   r2 
fluido fluido r
caliente frío  r2 
T∞1 h1 Ts1 Ts2 T∞2 h2
 2Lk (Ts1  Ts 2 )
Qr  por lo que
ln  r2 
En estado estable y sin generación.

1 d  dT   r1 
 kr 0 ln  r2 
r dr  dr   r1  1
Rtcond  y Rtconv 
Si k = Cte, integrando dos veces: 2Lk h 2rL

T r   C1 ln r  C2
 CILINDRO HUECO COMPUESTO (TUBO)
Un tubo con dos capas de otros materiales Otra forma:

T3 Ts4
T1  T 4
Qr   UA(T1  T 4 )
T2 B C T∞4, h4 Rtot
r1 A r2 r3
r4 U 1 r1  r2  r1  r3
1
 r1  r4  1 r1
 ln   ln    ln   
T∞1, h1 h1 k A  r1  k B  r2  kC  r3  h1 r4
Ts1
Considerando el concepto de resistencia
Se cumple que : U1 A1  U 2 A2  U 3 A3  U 4 A4
térmica en sistemas radiales, se puede

  Rt 
deducir la ecuación del calor radial como: 1

 T1  T 4
Qr 
r  r  r 
ln  2  ln  3  ln  4 
  1  2  3
1 r r r 1
2r1 Lh1 2k A L 2k B L 2kC L 2r4 Lh4
 EL RADIO CRÍTICO

Cuando se usa un aislante en un cilindro o un tubo, se reducen las pérdidas de Calor, se incrementa la
resistencia de conducción. También se tiene el efecto de incrementar el área transferencia de calor por
convección reduciendo la resistencia exterior de la película. Estos efectos se deben cuando al variar el
radio exterior del aislante. Considerando un tubo con una capa de aislante.
 T1  T 3
Qr 
r  r 
ln  2  ln  a 
  1  2
1 r r 1
2r1 Lh1 2kt L 2k a L 2ra Lh3
T∞1h1 T∞3 , h3
kt y k a Conductividades tubo y aislante ra
r2 L
Suponiendo que T1 = T2 = T∞t y que h1 y kt son muy grandes
T1 T2T3
 T1  T 3
Qr  r1
r 
ln  a 
 r2   1 r
2k a L 2ra Lh3
derivando respecto a ra e igualando a cero se obtiene :
k
ra  rcrítico  a
h 3
 LA ESFERA HUECA

Aplicando este método a una esfera Si la Rt se define como la diferencia de

Hueca, para un volumen de control Temperaturas dividida por la razón de calor.
Diferencial, la conservación de la energía 
r2 T
s2
requiere que.  
Qr dr
Q Q
r r  dr
4 r r 2  T k (T )dT ; k  Cte
1 s1

  4k (Ts1  Ts 2 )
Q r dr Qr 
1 1
r     

Ts1 Ts2  r1   r2 
Qr 1 1 1
Rtcond    
dr 4k  r1 r2 
1
Rtconv 
En estado estable, unidireccional sin 4r22 h
generación de energía.
 T 2  T1
 dT dT Q
Q r  kA  k (4r 2 ) 1 1 1  1 
dr dr      2 

 
con Q r  Q(r ) y Cte
4k  r1 r2   4r2 h 
 Ejemplo 2.2. Se tiene un tubo de vapor de diámetro exterior de 120 mm y aislado con silicato de
calcio con 20 mm. Las temperaturas Ts1 = 800 0K y Ts2 = 490 0K. Encuentre el calor radial / m.

SE ASUME: Condición de estado estable, unidimensional y k = Cte

PROPIEDADES: k = 0,089 w/m K.
DIAGRAMA:
Ts2 ANÁLISIS
Ts1

 Qr 2k (Ts1  Ts 2 )
Q´r  
L ln  DD2 
Vapor  1
 2 (0.089)(800  490)
Q´r   603w / m
 
0.16
ln 0.12

D1 = 0.12 m
D2 = 0.16 m COMENTARIO: El calor transferido fuera de
la superficie es disipado a los alrededores
por convección y radiación.
 c) CONDUCCIÓN CON GENERACIÓN DE ENERGÍA TÉRMICA

La pared plana.

q → Energía uniforme Gen / Vol T∞1 ; h1 T∞2 ; h2
Si k = Cte

2
d T q
2
 0
dx k Ts1


q
q 2
T  x  C1 x  C2 ; T ( L)  Ts1 ; T ( L)  Ts 2 Ts2
2k

Ts 2  Ts1 q 2 Ts1  Ts 2
C1  ; C2  L 
2L 2k 2
x

qL2  x 2  Ts 2  Ts1 x Ts1  Ts 2
T ( x)  1  L2   
2k   2 L 2
 CASO ESPECIAL
dT
Cuando: Ts1 = Ts2 = Ts Note que en x = 0 dx
 0 no hay transferencia de
x 0
 Calor a través del plano, puede representarse por una
2
qL  x 2  superficie adiabática. En x = L
T ( x)  1  L2   Ts
2k    dT
k  h(Ts  T ) consideran do T ( x)
qL2 dx
T (0)  T0   Ts xL
 
2k 2 2
dT qL 2x qL
T ( x)  T0  x   ( 2 ) ; Ts  T 
  dx 2k L xL h
Ts  T0 L

-L x L
T∞ h q T∞ h Qcond Qconv
Qcond Qconv
T0 T0 q
Ts Ts Ts T∞ h

x L
 CASO DE SISTEMAS RADIALES CON GENERACIÓN TÉRMICA

El cilindro. El modelo matemático es: Evaluando en r = 0; fluido frío

T∞ ,h Ts

1 d  dT  q
r   0 T(r = 0) = T0 Qr
r dr  dr  k r0

T (r )  Ts r
dT q
  r 2  C1  1   
r T0  Ts  r0  L
dr 2k
q
T (r )   r 2  C1 ln r  C2 Relacionando Ts a la temperatura del fluido frío T∞
4k
dT 
CI :  0 ; T (r0 )  Ts
dr r 0 q(2r02 L)  h(2r0 L)(Ts  T )

q 2 
C2  Ts  r0 ; C1  0 q r0

4k Ts  T 
q r02  r 2  2h 
T (r )  1  2   Ts q r0 q r02  r 2 
4k  r0  T (r )  T   1  2 
2h 4k  r0 
 Ejemplo 2.3. Se tiene un conductor de cobre calibre 12 (2.33 mm de diámetro). La
resistividad del cobre es de 1.73 x 10-8 Ώm, la conductividad térmica es de 380 w/mK y el
coeficiente de transferencia de calor de 10 10 w/m2 k. Determine la ecuación en función
de la corriente eléctrica de la diferencia de temperaturas máxima y del ambiente.

Para un cilindro la temperatura T( r ) tiene Se puede calcular el radio crítico si se forra

su valor máximo en el centro, cuando r = 0 el conductor con un material que tenga por
  ejemplo una ka = 0.11 w/mK.
q r0 q r02
T (r  0)  Tmax  T   k a 0.11
2h 4k Rcrítico    11mm
 2
i R i e
2
h 10
pero : q  2 e 
 
r0 L r02 2
i 2  e  2k  Es interesante evaluar ΔT para este caso
Tmax  T  T  2 
1  
4r0 k  hr0  del conductor aislado.

T 

i 2 1.73 x10 8  
 
2380  


4 2 1.165 x10 3 
2
1

380  10 1.165 x103  

T  0.1079i 2 en C
0
 d) ANÁLISIS DE ALETAS

Se usan aletas para incrementar el área de Haciendo el balance de energía:

contacto del fluido enfriador y así no
incrementar “h” por aumento de potencia.  dT d  dT 
Q x  dx  kAc  k  Ac dx
dx dx  dx 
dQconv

Qx dx
d Q conv  h.dAs (T  T ); dAs  área dif .
dAx Ac(x)
d  dT  h dAs
 Ac  (T  T )  0
dx  dx  k dx
x Qx+dx
d 2T  1 dAc  dT  1 h dAs 
     (T  T )  0
  
Q x  Q x  dx  d Q conv 2
dx  Ac dx  dx  Ac k dx 
 dT
Q x  kAc
dx
Es la ecuación generalizada de una aleta
Ac  Sección

  d Qx
Q x  dx  Q x  dx
dx
 ALETAS DE SECCIÓN UNIFORME

Cuando se tienen aletas como en el Def.

Diagrama Qconv
fluido d dT
T(0) = Tb ; T∞ → fluido T∞ , h  ( x)  T ( x)  T ; 
dx dx
Tb t
Ac = Cte Ac d 2 hP
As = Px Qf w  m 2  0 ; m 2 
dx kAc
x
L P = 2w+2t
 ( x)  C1 mx  C2   mx
P → Perímetro Qconv Ac= wt
d Con  (0)  Tb  T   b
As → área de base a “x” L P = πd
Qf Ac =πd2/4 y con x  L ; Se tiene :
dAc
 0; el mod elo queda
dx
d 2T hP
 (T  T )  0
dx 2 kAc
 CASO (A) Convección en el filo de la aleta

El calor fluye por conducción en la aleta y pasa Se nota que el gradiente de temperatura
a convección en su filo como muestra la figura decrece con “x” por la pérdida continua de
calor por convección en caras de la aleta.
dT
hAc [T ( L)  T ]   kAc
dx xL
Qconv
dT   dT d
h ( L )   k Q b  Q f  kAc  kAc
dx xL Tb dx x 0 dx x 0
 b  C1  C2 Qb = Qf  kAc
dT
 hAc [T ( L)  T ]
dx xL
 h 

Senh.mL   Cosh.mL
Q f  hPkAc  b  mk 
 h 
Cosh.mL    Senh.mL
h(C1 mL  C2   mL )  km(C1 mL  C2   mL )  mk 

Resolviendo para C1 y C2 Q f   hT ( x)  T dAs   h ( x)dAs
Af Af

 h 
Cosh.m( L  x)    Senh.m( L  x)
  mk  Af → Área total de aleta incluyendo el filo
 de la aleta.
b  h 
Cosh.mL    Senh.mL
 mk 
 OTROS CASOS DEL ANÁLISIS DE LA ALETA

CASO (B). Si la convección en el filo del aleta CASO ( C). Θ(L) = θL

es despreciable, se trata como adiabático.

d   L  b Senh.mx  Senh( L  x)
0 
dx xL
b Senh.mL

 Cosh.m( L  x)
 Cosh.mL   L  b 
 Q f  hPkAc  b
Senh.mL
b Cosh.mx

Q f  hPkAc  bTanh.mL CASO ( D ). L → ∞ ; θL → o


   mx
NOTA. Para usar los resultados del análisis
del CASO (A), se tiene que en la práctica es
b
válido si (mL) < 2.65. Si (mL) ≥ 2.65 se puede

usar la aproximación infinita.
Qf  hPkAc  b
 EJEMPLO 2.4. Una barra de bronce de 0.1 m largo y 0.005 m diámetro, se extiende
horizontalmente de una fundición a Tb = 200 0C.La barra está en el ambiente a T∞ = 200
C y h = 30 w/m2 K. ¿ cual es la temperatura de la barra a 0.025, 0.050 y 0.1 m ?. Bronce
a 110 0C; k = 133 w/mK

Diagrama. L = 0.1 m
Evaluando.
Aire a T∞ y h x1 = 0.025 m,
x2 = 0.050 m
Cosh.mL  2.04; Senh.mL  1.78
 h   30 
Tb d     0.0168
 mk   133x0.005 
x1 x2 L y con :  b  200  20  180
x :
12
  Cosh.m( L  x)  0.0168Seh.m( L  x)
12     (180)
 hP  
12

m  h d    4 h 2 .07
   kd 
 kAc   d 2   
 k   
  4  X( Cosh. Senh. θ T(0 )
12 m)
 4 x30  1
m(L-x) m(L-x)
   13.43m con mL  1.34 X 1.55 1.19 136.5 156.5
133 x 0 .005  1

 h  X2 1.24 0.725 108.9 128.9

Cosh.m( L  x)    Senh.m( L  x)
  mk  b
 h  L 1.00 0.00 87.0 107.0
Cosh.mL    Seh.mL
 mk 
 RENDIMIENTO DE ALETAS

Rendimiento de una aleta Eficiencia de una aleta “ηf”.

εf → Efectividad. Relación de la transferencia Af → Área de la superficie de la aleta.
de calor de la aleta a la razón de calor  

transferido si no existiera la aleta. εf > 2 para Qf Qf

f  
justificar las aletas. 
Q max hA f  f

En caso ( D ) 1 Aleta recta, área transversal uniforma y filo

 kP  2
 f    adiabático.
 hAc  MTanh.mL Tanh.mL
f   ; 0   f  1; L  
El rendimiento se puede evaluar en términos hPL b mL
de resistencia térmica.
Filo adiabático, sección recta o cilíndrica
b 1
Rtcond  
y Rtconv  t
Qf hAcb Lc  L   Secc recta
2
Rtcond d
f  Lc  L   Secc cilíndrica
Rtconv 4

Acb → Área de sección transversal de aleta Q f  MTanh.mLc ; M  hPkAc . b
en su base. Tanh.mLc
f 
mLc
 RENDIMIENTO DE ALETAS II

Errores con la aproximación despreciables si: Aleta sección transversal no-uniforme:

Caso de secc. anular; Ac = 2πrt, varía con “r ”

ht hd Reemplazando “r” por “x” en Ec. de calor.
ó  0.0625
k 2k
d 2T 1 dT 2h
Si w  t  P  2 w y   (T  T )  0
dr 2 r dr kt
 
1 1

 hP  As  2 r 2  r12  Área de la Sup

 Lc    Lc
2
2h 2

mLc  
 kAc   kt  2
con m 2  y   T  T
 L2 
1 kt
multiplica ndo por  c1  y con Ap  Lc t d 2 1 d
 2   2
 0
 Lc  dr 2
r dr
m

 2h  3 solución :  (r )  C1 I 0 (mr )  C2 K 0 (mr )

mLc    Lc2
 kA 
 p I 0 y K 0 funciones de Bessel orden cero
primera y segunda clase
 RENDIMIENTO DE ALETAS III

Eficiencia secc transversal no uniforme

 I (mr ) K1 (mr2 )  K 0 (mr ) I1 (mr2 )

 0
 b I 0 (mr1 ) K1 (mr2 )  K 0 (mr1 ) I1 (mr2 ) t
 dT d
Q f  kAcb  k (2r1t ) r2c = r2 + t/2
dr r  r1 dr r  r1
r1 L Lc = L + t/2
K1 (mr1 ) I1 (mr2 )  I1 (mr1 ) K1 (mr2 ) Ap = Lct
 (2kr1t m b )
K 0 (mr1 ) I1 (mr2 )  I 0 (mr1 ) K1 (mr2 ) r2

Qf 2r1 K1 (mr1 ) I1 (mr2 )  I1 (mr1 ) K1 (mr2 )
f  
2h(r22  r12 ) b m(r22  r12 ) K 0 (mr1 ) I1 (mr2 )  I 0 (mr1 ) K1 (mr2 )

1
Rtaleta   Re sistencia térmica de aleta
hA f  f
 EFICIENCIA DE SUPERFICIE GLOBAL
Se tienen “N” aletas en un equipo térmico, la eficiencia de superficie global es:

 
Qf Qf
o  

Q max hA f  b

A f  Área de aleta  Área porción exp uesta

AT  NAf  Ab

N  Número de aletas
AT  ÁreaTotal

Q T  N f hA f  b  hAb b

NAf
o  1  (1   f )
AT
 Problema: Vapor de agua fluye por tubo Dext = D1 = 3 Cm a T = 120 0C. El tubo tiene aletas
circulares de Al (k = 180 w/m 0C) de D2 = 6 Cm y espesor, t = 2 mm. El espacio entre aletas
es de 3 mm por lo que son 200 aletas/m. El aire exterior está a T∞ = 25 0C, h = 60 W/m2 K.
Determine el incremento de la TC del tubo/m por la adición de las aletas.

Análisis: Con estos datos en la gráfica de eficiencia para esta aleta:

η = 0.96
Si no se tienen aletas:
Asa = πD1L = π (0.03)(1) = 0.0942 m2 Aaleta  2 (r22c  r12 )  2 (0.0312  0.0152 )  0.004624m 2
 
 Q aleta   aleta Q max   aletahAaleta (Tb  T )
Q sa  hAsa (Tb  T )  60(0.0942)(120  25)
 0.96(60)(0.004624)(120  25)  25.3W
 537W
La TC en parte libre de aletas es:
Para aletas circulares sujetas a un tubo en
una gráfica se tiene: Alibre  D1S   (0.03)(0.003)  0.000283m 2

L  12 ( D2  D1 )  12 (0.06  0.03)  .015m Q libre  hAlibre (Tb  T )  60(0.000283)(120  25)
r2 c  r2  2t  0.03  0.002
2
 0.031m  1.6W
Lc  L  2  0.015  2  0.016m
t 0.002

Ap  Lc t  0.016(0.002)  3.2 x10 5 m 2 Se tienen 200 aletas / m 

  

r2 c 0.031 Q T  200(Q aleta  Q libre )  200(25.3  1.6)  5380W

  2.07
r2 0.015   
Q incremento  Q T  Q sa  5380  537  4843W
3
3 h 60 
 L 2
 (0.016) 2  2.07 QT 5380
c
kAp 80(3.2 x10 5 )  aleta  
  10
Q sa 537