Capítulo 3

CONDUÇÃO DE CALOR 1-D,

REGIME PERMANENTE CALOR

∂ ! ∂T $ ∂ ! ∂T $ ∂ ! ∂T $ ∂T
Parede plana # k & + # k & + # k & + q! = ρ c p
∂ x " ∂ x % ∂ y " ∂ y % ∂z " ∂z % ∂t

T¥1 d " dT % dT
$ k ' = 0 ⇒ qx = −kA = cte
dx # dx & dx
Ts1
qx ⇒ T( x) = ax + b
Ts2 *T(0) = Ts1 Ts2 − Ts1
CC : + ⇒ b = Ts1 a=
,T(L) = Ts2 L
T¥2
Ts2 − Ts1
T( x) = x + Ts1
x=0 x=L L
k
qx = − (Ts2 − Ts1 ) A
Balanço de energia L
"
numa superfície: qcond "
- qconv =0
q x = h1 A(T¥1 - Ts1 ) q x = h2 A(Ts 2 - T¥ 2 ) 1

€
Parede plana, sem fonte
qx
L q x = h1 A(T¥1 - Ts1 ) = (T¥1 - Ts1 )
T¥1 h1 A
k A
Ts1 qx = (Ts1 - Ts 2 ) qx
= (Ts1 - Ts 2 )
qx
Ts2 L ( k A) / L
qx
k T¥2 q x = h2 A(Ts 2 - T¥ 2 ) = (Ts 2 - T¥ 2 )
h2 A
x=0 x=L
q x = UA(T¥1 - T¥ 2 ) qx
= (T¥1 - T¥ 2 )
UA
U = Coeficiente global de
transferência de calor

æ 1 L 1 ö÷
somando q x çç + + ÷ = (T¥1 - T¥ 2 )
è h1 A kA h2 A ø
2
1 æ 1 L 1 ö÷
ç
=ç + +
UA è h1 A kA h2 A ÷ø
 Resistência Térmica - analogia com resistência elétrica

qx T¥1 Ts1 Ts2 T¥2

1/h1A L/kA 1/h2A

Resistência térmica à condução Resistência térmica à convecção

Ts1 − Ts 2 L q = hA(Ts − T∞ )
Rt ,cond ≡ =
qx kA T −T 1
Rt ,conv ≡ s ∞ =
q hA
Circuito equivalente
€ Resistência Total
qx T¥1 T∞2
€ T∞1 − T∞2 1 L 1
Rtot = = + +
1/(UA) qx h1 A kA h2 A
qx=qc1=qk=qc2 3
 Resistência térmica à radiação
Ts − Tamb 1
Rt ,rad = =
q rad hr A
hr ≡ εσ (Ts + Tamb )(Ts2 + Tamb
2
)
(
q rad = esA Ts4 - Tamb
4 )
(Ts4 - Tamb
4 ) = (T - T
s amb )(Ts + Tamb )(Ts + Tamb )
2 2
€
Paredes compostas em série
T∞1 − T∞2 T∞1 − T∞2
qx = =
T¥1 ∑ Rt 1/ h1A + LA / kA A + LB / kB A + LC / kC A + LD / kD A + 1/ h2 A
qx
T¥2 Coeficiente Global de Troca de Calor
€
qx = UAΔT
A B C D 1 1
U≡ =
Rtot A 1/ h1 + LA / k A + LB / kB + LC / kC + LD / k D + 1/
4
h2
 Paredes em paralelo

B
qx A D
Exercício: Obtenha a resistência térmica total
C para a parede plana composta pelos
materiais A, B, C, D

L2/(kBAB)
qx qx
Rtot , serie = å Rt
qx Ts1 Ts2 qx1 Ts3 Ts4
qx2 1 1
L1/(kAAA) L3/(kDAD) =å
Rtot , parelelo Rt
L2/(kCAC)
5
Resistência térmica de contato
-Provoca queda de temperatura na interface entre 2 paredes
-Ocorre principalmente devido à rugosidade das superfícies
-Valores são obtidos experimentalmente

A B
A B T A - TB
qx TA TB Rt , contato º
q ¢x¢

6
 Exemplo
Um fabricante deseja projetar um forno auto-limpante, que utiliza
uma janela composta, separando a cavidade do forno do ar ambiente.
A janela consiste de 2 plásticos (LA=2LB, kA=0,15W/mK e
kB=0,08W/mK). Durante o funcionamento, a parede do forno e a
temperatura do ar interno valem 400 0C, enquanto o ar ambiente esta
a 25 0C. Os coeficientes convectivos interno e externo, e o de
radiação (hi, he e hr) valem 25 W/m2K. Qual a espessura mínima da
parede (L= LA+LB) necessária para garantir que a temperatura na
superfície externa da janela não ultrapasse 50 0C?

A B

Forno
7
 qci,n
T¥1
qkA qkB qc,e
T¥1=400 Tsin TsA Ts,e T¥2
T¥2=25 1/(hiA)
A B qr,i,n
T¥1 LA/(kAA) LB/(kBA) 1/(heA)
T¥1
Forno

1/(hrA)
hi= he = hr = 25 W/m2K L= LA+LB=?
LA=2LB, kA=0,15W/(mK) e kB=0,08W/(mK) Para Ts,e < 50 oC

T - T¥ 2 50 - 25 W
q x = qc , in + q r , in = q k , A = q k , B = qc , e qc , e = se q ¢x¢ = = 625
1 /( he A) 1 / 25 m2

(T¥1 - Tsin )
q x = qc , in + q r , in = (hi A + hr A )(T¥1 - Tsin ) =
1 /[ hi A + hr A]
T¥1 - T¥ 2
qx =
1 /[ hi A + hr A] + L A /( k A A) + LB /( k B A) + 1 /( he A)

400 - 25 375
q ¢x¢ = = = 625
1 /[ 25 + 25] + 2 LB / 0,15 + LB / 0,08 + 1 / 25 0,06 + 25,833 LB

375 375
- 0,06 - 0,06
625 L = L A + LB = 3 LB = 625 = 0,06271 m
LB = = 0,020903m 25,833 8
25,833
 Sistemas radiais: % ∂T
q"= −k ' eˆ r +
1 ∂T
eˆθ +
∂T (
eˆ z *
& ∂r r ∂θ ∂z )
Paredes cilíndricas 1 ∂ % ∂ T ( 1 ∂ % ∂T ( ∂ % ∂ T ( ∂T
'kr * + 2 ' k * + 'k * + q˙ = ρc p
Cilindro: regime permanente, r ∂r & ∂r ) r ∂θ & ∂θ ) ∂z & ∂z ) ∂t
sem geração

1 d " dT % €
$kr ' = 0
r dr # dr &
dT dT
q r = −kA = −k (2πrL) = q "r (2πrL)
dr dr   
função de r
independe de r

Para k constante :
T(r) = a lnr + b
*T(r1 ) = Ts1
CC : +
,T(r2 ) = Ts2
Ts1 − Ts2 " r % 2πLk(Ts1 − Ts 2 )
⇒ T(r) = ln$ ' + Ts2 qr =
ln(r1 /r2 ) # r2 & ln(r2 /r1 )
ln(r2 /r1 )
Rt ,cond = 9
2πkL
Paredes cilíndricas em série

10
Exercício: Raio crítico de isolamento em sistemas radiais
Um tubo de cobre de parede fina transporta refrigerante a baixa temperatura Ti. O
ambiente externo está a T∞ . Calcule a resistência térmica total para um tubo de 10
mm de diâmetro e espessuras de isolamento: 0, 5, 10, 20 e 40 mm. O coeficiente de
troca de calor é 5 W/m2K, e o comprimento do tubo é 10 m.

T∞,h r
ri

Ti k

Aumento de r (ou da espessura de isolante): resistência a condução

cresce, e resistência a conveção cai
Hipóteses: Regime permanente; condução 1-D radial; resistência térmica da
parede do tubo desprezível
11
 qx Ti T∞ T∞,h
ri r
RCD RCV

Δ$ $( − $* quanto maior R menor q e Ti k

!= =
∑ &' &+, + &+. menor perda de energia
2
/0(2 ) <=> 9 9
&+, = 3 cresce com r =0: − =0
5678 <: 5678: 5678:?
&+. =
9
cai com r C
56:8; → AB: =
ℎ
r-ri r RT
(mm) (mm) (K/W)
0 5 0,637
5 10 0,518
10 15 0,530
20 25 0,593
Existe uma espessura crítica de isolamento,
40 45 0,706
abaixo da qual q cresce com o aumento de r, e acima dela
q cai com o aumento de r: rcr=k/h 12
 & ∂T 1 ∂T 1 ∂T )
Paredes esféricas q"= −k ( eˆ r +
' ∂ r r ∂θ
eˆθ +
r sin θ ∂φ
eˆφ +
*
1 ∂ & 2 ∂T ) 1 ∂ & ∂T ) 1 ∂ & ∂T ) ∂T
( kr + + ( k + + ( k sin θ + + ˙
q = ρc p
r 2 ∂r ' ∂r * r 2 sin 2 θ ∂φ ' ∂φ * r 2 sin θ ∂θ ' ∂θ * ∂t

qr+dr dT q r2 dr Ts 2
€ q r = −k ( 4 πr ) 2
⇒ r ∫ 2
= ∫ k(T)dT
qr dr 4π r1 r Ts1

Para k constante :
4 πk (Ts1 − Ts 2 )
qr =
1/r1 −1/r2
(Ts1 − Ts2 ) 1 &1 1)
Rt ,cond = = ( − +
qr 4 πk ' r1 r2 *

13
Resumo

14
Exemplo: Reservatório esférico de nitrogênio com vazamento.
Calcular a taxa de transferência de calor para o nitrogênio e
o fluxo de massa que evapora. Considere que a temperatura
da parede interna é igual a do Nitrogênio.
ri=0,25 m
Parede ro=0,275 m
isolada
Pó de silica
K=0,0017 W/mK

Ar
T=300K q Nitrogênio líquido
T=77 K
h=20W/m2K r=804 kg/m3
hlf=2x105 J/kg (calor latente vaporização)

15
16
17
Condução 1D, regime permanente, com geração de energia:

Eq. condução: ∂ ! ∂T $ ∂ ! ∂T $ ∂ ! ∂T $ ∂T
# k & + # k & + # k & + q! = ρ c p
∂ x " ∂ x % ∂ y " ∂ y % ∂z " ∂z % ∂t

Para condutividade constante, e geração uniforme:

d 2T q˙ q˙ 2
2
+ = 0 ⇒ T(x) = − x + ax + b
dx k 2k
$T(x = −L) = Ts1
CC : %
&T(x = L) = Ts2
q˙L2 ' x 2 * Ts 2 − Ts1 x Ts 2 + Ts1 18
⇒ T(x) = )1− 2 , + +
2k ( L + 2 L 2
Obs.: - O fluxo de calor qx (=-kAdT/dx) é função de x
- Se as temperaturas das paredes são iguais, o problema é
simétrico:
! 2 " x2 %
qL
T (x) = $1− 2 ' + Ts
2k # L &
! 2
qL
Tmax em x=0: Tmax = T0 = + Ts
2k

- No plano de simetria, dT/dx=0 ®não há fluxo de calor

Þsuperfície adiabática
2
" %
Temperatura adimensional: θ (x) = T (x) − T0 = $ x ' = ( x * ) x
2 *
x =
Ts − T0 #L& L

Balanço na superfície, para dT

−k = h(Ts − T∞ )
CC de convecção: dx x= L
19
 1 ∂ ! ∂T $ 1 ∂ ! ∂T $ ∂ ! ∂T $ ∂T
Sistemas radiais: # kr & + 2 # k & + # k & + q! = ρ c p
r ∂ r " ∂ r % r ∂θ " ∂θ % ∂ z " ∂ z % ∂t
cilindro longo
1 d " dT % q˙
$r ' + = 0
r dr # dr & k
T¥,h qr q˙
T(r) = − r 2 + a lnr + b
4k
)T(r0 ) = Ts
q˙ + q˙r02 " r 2 %
CC : * dT ⇒ T(r) = − $1− 2 ' + Ts
+ dr =0 4k # r0 &
, r= 0
€ Temperatura adimensional :
Ts
2
T(r) − Ts " r %
θ (r) = = 1− $ '
T0 − Ts # r0 &
Balanço global :
q˙πr02 L = h2πr0 L(Ts − T∞ )
q˙r0
Ts = T∞ +
2h
20
Exemplo: Tubo sólido longo, isolado termicamente, com
condutividade da parede igual a k
e geração de calor uniforme na parede igual a q˙ (W/m3).
O tubo tem raio interno r1 e externo r2.
1. Obter a distribuição de temperaturas na parede do tubo,
considerando que a temperatura da superfície externa
é igual a Ts2. €
2. Calcule o fluxo de calor por unidade de comprimento,
na superfície interna do tubo.

21
22
23
24
Condução de Calor em Superfícies Extendidas (Aletas):

Balanço de calor:
q x = q x +dx + dq conv
dT dq
q x = −kAc q x +dx = q x + x dx
dx dx
dT d $ dT '
⇒ q x +dx = −kAc − k & Ac )dx
dx dx % dx (
dq conv = hdAs (T − T∞ )
d $ dT ' h dAs
⇒ & Ac )− (T − T∞ ) = 0
dx % dx ( k dx
ou
d 2T 1 dAc dT $ 1 h dAs '
+ −& )(T − T∞ ) = 0
dx 2 Ac dx dx % Ac k dx (

Equação geral energia

€ para superfícies extendidas
25
Aletas de seção reta constante Þ Ac=cte As=Px
ÞdAc/dx=0 e dAs/dx=P
Então a equação de energia fica:

d 2T # hP &
2
−% ((T − T∞ ) = 0
dx $ kAc '

Def.: q(x) = T(x) - T¥ Þ dq/dx = dT/dx

Então€ a eq. energia fica:
d 2θ 2 hP
2
− m θ =0 m2 =
dx kAc
sol. :
θ (x) = C1e mx + C2e −mx

26
CC: x=0 T=Tb q= Tb - T¥ º qb
x=L: são possíveis 4 CC diferentes
Caso A: convecção
dT
hAc (T(L) − T∞ ) = −kAc
dx x= L
ou
dθ
hAcθ (L) = −kAc
dx x= L

Assim,
θ cosh(m(L − x)) + (h /mk)sinh(m(L − x))
=
θb cosh(mL) + (h /mk)sinh(mL)
dT
q f = q b = −kAc = −kAc
dθ Calor perdido
dx x= 0 dx x= 0 Pela aleta
sinh(mL) + (h /mk)cosh(mL)
q f = hPkAc θ b
cosh(mL) + (h /mk)sinh(mL)
27
Caso B: convecção desprezível na extremidade da aleta
Þdq/dx½x=L =0

θ cosh(m(L − x))
=
θb cosh(mL)
dT dθ
q f = q b = −kAc = −kAc
dx x= 0 dx x= 0

q f = hPkAc θ b tanh(mL)

Caso C: temperatura prescrita na extremidade da aleta

€ Þ q(L)= qL
θ θ L /θ b sinh(mx) + sinh(m(L − x))
=
θb sinh(mL)
cosh(mL) − θ L /θ b
q f = hPkAc θ b
sinh(mL)
28
 Caso D: Aleta muito longa
L®¥, qL®0
θ
= e −mx
θb
q f = hPkAc θ b

Exemplo: Barra longa exposta ao ar ambiente

€

T¥=25 0C
1. Determine a distribuição de temperatura na
Tb=100 0C h= 100 W/m2K barra e as perdas de calor pela barra,
considerando 2 materiais: aço (k=14W/mK)
D=5mm
e cobre (k=398 W/mK).
2. Estime o comprimento que a barra deve ter
para que a hipótese de barra infinita seja
satisfatória. 29
Solução:
Hipóteses: Reg. permanente, condução 1D, propriedades ctes,
radiação desprezível, L
CC: caso D
Então:

! = !# + !% − !# ' ()*
m=
+, 0/2
=
3+ 0/2 Cobre Aço
-./ -4
m 14,2 m-1 75,6 m-1
Calor perdido pela barra:
qf 8,3 W 1,6 W
56 = ℎ89:; !% − !#

30
Variação da temperatura na aleta para os dois materiais

Cobre Aço

m 14,2 m-1 75,6 m-1

qf 8,3 W 1,6 W

x(mm)

Para verificar a validade da hip L vamos comparar esta solução

(CC caso D) com a solução que teríamos se usássemos a CC do caso
B (convecção desprezível na extremidade)
31
Calor perdido pela barra:
CC caso D CC caso B
!" = ℎ%&'( )* − ), !" = ℎ%&'( )* − ), tanh(mL)

Assim, para que o resultado obtido seja próximo ao do caso B,

tanh(mL)≥0,99 , ou (mL) ≥ 2,65
Logo, Lcr=2,65/m
Laço Lcobre
Para que a hipótese seja satisfatória:
Para o cobre, L ≥ 190 mm
Para o aço, L ≥ 40 mm

x(mm)

32
 Desempenho das aletas: as aletas aumentam a troca de calor
crescendo a área superficial de troca. Porém, a aleta é também
uma resistência adicional a troca de calor original.
Efetividade da aleta
qf troca de calor na aleta
εf = =
hAc,bθ b troca calor que existiria sem aleta

• ef deve ser o maior possível. ef < 2 o uso da aleta não se justifica

€ • Para o caso D: ef =(kP/hAc)1/2
ÞMaior k ® maior ef
ÞMaior P/Ac (aletas mais finas) ® maior ef
ÞMenor h (gases) ® maior a utilização de aletas
Þ ef máximo quando L®¥
33
 Resistência térmica da aleta

θb
Rt , f ≡
qf
Rt ,b 1
εf = Rt ,b ≡ (= resistência à convecção na base exposta)
Rt , f hAc,b

Eficiência da aleta
€
qf qf
ηf ≡ =
q max hA f θ b
fluxo por convecção caso toda
a aleta estivesse a Tb

€
34
Eficiência global de superfície: caracteriza o desempenho de um
conjunto de aletas dispostas sobre a superfície

qt NA f
η0 =
q max
= 1−
At
(1− η f )

= NA f + Ab

q t = η f NhA f θ b + hAbθ b

35