Transmissão de Calor – Convecção Forçada

P.J. Oliveira
Departamento Engenharia Electromecânica, UBI,
Dezembro 2014

1. Introdução
Na convecção forçada o fluxo de calor qɺc entre uma parede aquecida à temperatura Tw
e um fluido mais frio (temperatura média, longe da parede, T∞ ) é determinado a partir
da lei de Newton do arrefecimento:
qɺc = h(Tw − T∞ ) [W/m2]
O número de Nusselt é o parâmetro adimensional principal, fornecendo o coeficiente
convectivo que permite calcular o calor transferido:
hL
u =
k
com significado físico U% T%
fluxo por convecção h∆T
u = = . .
fluxo por condução pura k ∆T / L qc
Tw
em que
h - coeficiente convectivo [W/(m2 K)]
L - dimensão característica [m]
k - condutibilidade térmica do fluido [W/(m K)]
∆T - diferença de temperaturas, entre a parede e o fluido longe da parede [K]

O movimento do fluido em convecção forçada é gerado por uma fonte de energia

externa (bomba, ventilador, etc) que permite definir uma escala de velocidade típica:
U ∞ - velocidade longe da parede [m/s]
Esta velocidade aparece num novo número adimensional, o número de Reynolds,
importante em convecção forçada (e de facto em toda a mecânica de fluidos), definido
como a razão entre as forças de inércia (convectivas) promotoras do movimento, e as
forças de viscosidade que se opõem ao movimento:
forças de inércia ρUU ρUL ρU ∞ L
Re = = = ⇒ Re =
forças de viscosidade µU / L µ µ
O escoamento é laminar quando o número de Reynolds é inferior a um valor
crítico, Recr , que depende da geometria, e transita para o regime turbulento quando Re é
superior a esse valor crítico.

A viscosidade é definida através da lei de Newton para fluido viscoso, que relaciona a
força por unidade de área (tensão, τ = F / A ) que uma camada de fluido exerce sobre a
que lhe está adjacente, a uma distância ∆y , quando a diferença de velocidades entre as
duas é ∆u :
∆u ∂u
τ =µ ⇒τ = µ
∆y ∂y
O coeficiente de viscosidade µ [Pa.s] = [kg/(m.s)] é constante (a uma determinada
temperatura) para um fluido newtoniano (como o ar, ou a água). Mede a facilidade do
deslizamento entre camadas de fluido quando são submetidas a uma determinada força
de corte (tangencial). Se a temperatura aumentar, a viscosidade dos gases aumenta

-1-
(maior tendência para choques moleculares) e a viscosidade dos líquidos diminui
(menor ligação entre moléculas).

Camada limite. Qualquer fluido usual adere a uma parede sólida que o contenha. Num
escoamento sobre uma superfície sólida estacionária (imóvel), existe uma fina camada
de fluido na qual as velocidades variam
abruptamente, desde a velocidade imposta a infinito y
U ∞ pela força motriz que gera o movimento, até a U%
U%
velocidade nula na parede resultante da condição de
não escorregamento ( u = 0 para y = 0 ). Toma-se x camada limite
laminar
turbulenta
como a distância segundo o movimento, a que δ(x)
corresponde uma componente de velocidade u , e y
0 x
a distância normal à parede. Esta zona de elevados xcr
gradientes de velocidade ∂u / ∂y designa-se como
camada limite (conceito introduzido por Prandtl). O escoamento fora da camada limite é
potencial (o fluido comporta-se como ideal, e a viscosidade é irrelevante); é somente
dentro da camada limite que a viscosidade se torna importante. A espessura da camada
limite δ é definida como a distância da parede para a qual a velocidade fica 1 %
diferente da velocidade a infinito: δ = y para u = 0.99U ∞ . Por causa da condição de
não escorregamento, o fluxo de calor convectivo é igual ao fluxo por condução
calculado exactamente sobre a parede: qɺc = qɺcond , w = − k (∂T / ∂y ) w .

Camada limite térmica. Quando a superfície sólida (parede) é aquecida (ou arrefecida
relativamente à temperatura do fluido a infinito, longa da parede), existe uma camada
limite de temperaturas semelhante à camada limite de velocidades. Nessa camada de
fluido a temperatura varia muito rapidamente desde a temperatura imposta na parede até
à temperatura do fluido não perturbado, longe da parede. A espessura dessa camada
limite térmica é indicada por δ t ( y para T − Tw = 0.99(T∞ − Tw ) ), e verifica-se que a
razão entre as duas espessuras é função do número de Prantdl,
µc p µ/ρ ν
Pr = = =
k k / (ρcp ) α
em que
ρ - massa volúmica [kg/m3]
c p - capacidade térmica a pressão constante [J/kg K]
ν - viscosidade cinemática [m2/s]
α - difusividade térmica [m2/s]

Fisicamente o número de Prandtl representa

ν difusão da quantidade de movimento (velocidades)
Pr = =
α difusão do calor (temperatura )
Quando:
Pr ≫ 1 : o campo de velocidade desenvolve-se rapidamente (óleos viscosos);
Pr ≪ 1 : o campo de temperatura desenvolve-se rapidamente (metais líquidos);
Pr ≈ 1 : os campos de velocidade e temperatura desenvolvem-se
simultaneamente (gases).

-2-
Valores típicos do número de Prandtl:

Fluido Pr
Metais líquidos 0.004 – 0.03
Gases 0.7 – 1.0
Água 1.7 – 13.7
Fluidos orgânicos 5 – 50
Óleos 50 – 100 000
glicerina 2000 – 100 000

No caso de camada limite sobre placa plana, demonstra-se que:

3
δ 
Pr ≅  
 δt 

2. Escoamentos exteriores (camada limite)

2.1 Placa plana de comprimento L (geral): u = CRem Pr n Recr = 5 × 105

y T%
U%
a. Convecção forçada sobre placa isotérmica ( Tw )
u(y) T(y) δ (y)
t
δh
#usselt médios ( u = hL / k ; Re = U ∞ L / ν ): Tw x

u = 0.664 Re1/ 2 Pr1/3 laminar, Re < 5 × 105 ; Pr ≥ 0.6

u = 0.037 Re 4/5 Pr1/3 turbulento, 5 × 105 < Re < 107 ; 0.6 ≤ Pr ≤ 60

u = ( 0.037 Re 4/5 − 871) Pr1/3 parte laminar, resto turbulento, 5 × 105 < Re

#usselt locais ( u x = hx x / k ; Rex = U ∞ x / ν ):

u x = 0.332 Rex1/ 2 Pr1/3 laminar, Rex < 5 × 105 ; Pr ≥ 0.6 (Nota: h = 2hx )

u x = 0.0296 Rex 4/5 Pr1/3 turbulento, 5 × 105 < Rex < 107 ; 0.6 ≤ Pr ≤ 60
L
1
(Nota: h = ∫ hx dx = 54 hx )
Lo
Coeficientes de atrito médios ( C f = τ w / ( 12 ρU ∞2 ) ):
1.328
Cf = laminar, Re < 5 × 105 ; Pr ≥ 0.6
Re1/ 2
0.074
Cf = turbulento, 5 × 105 < Re < 107 ; 0.6 ≤ Pr ≤ 60
Re1/5
C f = 0.074 / Re1/5 − 1742 / Re parte laminar e restante parte turbulenta

b. Convecção forçada sobre placa com fluxo de calor ( qɺ w ) imposto:

u x = 0.453Re1/x 2 Pr1/3 laminar, Re < 5 × 105

-3-
u x = 0.0308 Rex0.8 Pr1/3 turbulento, Re > 5 × 105

No caso laminar, a diferença de temperaturas média entre a parede e o fluido é dada por:
∆T = Tw − T∞ = (qɺw L / k ) / (0.6795Re1/2 Pr 1/3 ) ⇔ u = 0.6795 Re1/2 Pr 1/3

2.2 Convecção forçada sobre cilindro ( Re = U ∞ D / ν , D -diâmetro; Recr = 2 × 105 )

Nusselt médio ( u = hD / k ):
4/5
0.62 Re1/2 Pr 1/3   Re 5/8 
u = 0.3 + 1 +    RePr > 0.2
(1 + (0.4 / Pr ) ) 2/3 1/ 4
  28 200  
(Churchill e Bernstein; propriedades a T f = 0.5(T∞ + Tw ) )
Coeficiente de resistência (“drag”) médio (cilindro: diâmetro - D; comprimento - L):
6.8 1.96 0.0004 Re
CD = 1.18 + 0.89 + 0.5 − 10 −4 < Re < 2 × 105
Re Re 1 + 3.64 ×10 −7 Re 2
A força de resistência é: FD = CD Af 12 ρU 2 (área frontal do cilindro: A f = LD ).

2.3 Convecção forçada sobre esfera de diâmetro D ( Re = U ∞ D / ν )

Nusselt médio ( u = hD / k ):
1/ 4
 µ∞ 
u = 2 + ( 0.4 Re 1/ 2
+ 0.06 Re 2/3
) Pr 0.4
  3.5 < Re < 80 000 ; 0.7 < Pr < 380
 µw 
(Whitaker; propriedades a T∞ )
Coeficiente de resistência (“drag”) médio (F.A. Morrison)
24 2.6( Re / 50) 0.411( Re / 263 000)−7.94 Re0.80
CD = + + −8.00
+ 10 −4 < Re < 106
Re 1 + ( Re / 50) 1.52
1 + ( Re / 263 000) 461 000
A força de resistência é: FD = CD Af 2 ρU (área frontal da esfera: Af = π D 2 / 4 ).
1 2

2.4 Convecção forçada sobre corpos de várias formas ( Re = U ∞ L / ν )

Fórmulas aproximadas, com u = hL / k = CRe m Pr 1/3 ; L- dimensão transversal
Geometria Fluido L C m Validade Re
Cilindro Gás ou D 0.989 0.330 0.4 – 4
circular líquido 0.911 0.385 4 – 40
0.683 0.466 40 – 4000
0.193 0.618 4000 – 40 000
0.027 0.805 40 000 – 400 000
Quadrado Gás a 0.102 0.675 5000 – 100 000
(lado a)
Quadrado, 45º Gás 0.246 0.588 5000 – 100 000
2a
Hexágono Gás 0.153 0.638 5000 – 100 000
2a/ 3
(alinhado,a-lado)
Hexágono, 45º Gás 2a 0.160 0.638 5000 – 19 500
0.0385 0.782 19 500 – 100 000
Placa vertical Gás H 0.228 0.731 4000 – 15 000
Elipse (alinhada Gás b 0.248 0.612 2500 – 15 000
b- eixo menor)
(Propriedades a T f = 0.5(Tw + T∞ ) )

-4-
 3. Escoamentos interiores

Em escoamentos interiores a velocidade típica usada na definição do número de

Reynolds é a velocidade média dentro da conduta, obtida por divisão do caudal
volumétrico ( Vɺ [m3/s]) pela área da secção transversal ( A [m2]):
U = Vɺ / A [m/s]
Para fluido com massa volúmica ρ [kg/m3], o caudal
U Tw u(r)
mássico numa conduta é dado por:
δ D
mɺ = ρ AU [kg/s] .
m=ρAU x
O escoamento no interior dum tubo circular de diâmetro camada Lh
D = 2 R ocorre no regime laminar (as partículas de limite
fluido seguem linhas paralelas bem definidas) quando L
Re = ρUD / µ < 2300 (laminar)
e no regime turbulento (as linhas instantâneas de escoamento são irregulares e a
velocidade local varia aleatoriamente ao longo do tempo) se
Re > 4000 (turbulento).
Entre estes dois valores de Re o escoamento está num regime de transição, entre laminar
e turbulento:
2300 < Re < 4000 (transição).

Para um tubo redondo, de raio R e área transversal A = π R 2 , o perfil de velocidades no

regime laminar, em escoamento completamente desenvolvido (a velocidade já não
depende da coordenada axial, ao longo do tubo) tem forma parabólica,
u = 2U 1 − (r / R ) 2 
r
com a velocidade máxima a ocorrer no eixo ( r = 0 ), Tw Tw
igual a: u(r)
U D U
U 0 = 2U x
Ti T(r) Te
e velocidade média:
∆p R 2
U = (− dp / dx) R 2 / 8µ ⇒U =
L 8µ
em que − dp / dx é o gradiente de pressão (constante) que faz mover o fluido.

O fluido que circula no interior de uma conduta pode, idealmente, ser aquecido por duas
formas:
- a parede da conduta está quente, a uma temperatura constante, Tw ;
- um fluxo de calor constante é imposto na parede qɺ w .

A potência calorífica recebida pelo fluido, entre a entrada e a saída da conduta, é dada
por um balanço de energia (1ª Lei da Termodinâmica):
ɺ p (Te − Ti ) [W]
Qɺ = mc
em que a temperatura média do fluido à entrada é Ti , à saída é Te , e c p é a capacidade
térmica específica do fluido a pressão constante. Esta potência calorífica transferida
pode também ser obtida através de aplicação dos princípios de transmissão de calor,
usando um coeficiente convectivo h (assumido como constante):
Qɺ = Ap h∆T [W]

-5-
em que a área de permuta, Ap , é em geral dada por Ap = ∫ Pdx ( P =perímetro molhado
da conduta) e, no caso particular de um tubo redondo, Ap = 2π RL , com L -
comprimento do tubo. Como a temperatura média do fluido vai aumentando ao longo
do tubo, à medida que é aquecido, a diferença de temperaturas na expressão anterior
( ∆T ) tem de assumir um valor médio.
#o caso de Tw = Cte , uma hipótese grosseira é utilizar a diferença entre a temperatura
na parede e a temperatura global média do fluido ( Tm = 0.5(Ti + Te ) ), ∆T = ∆Tm :
∆Tm = Tw − Tm
e um hipótese mais correcta será usar uma média
T
logarítmica ∆T = ∆Tln :
Tw
∆Ti − ∆Te ∆Te
∆Tln =
 ∆T  ∆Ti Te
ln  i 
 ∆Te  Ti
com:
0 x
∆Ti = Tw − Ti e ∆Te = Tw − Te . L
Demonstra-se que esta segunda hipótese, da média logarítmica, resulta teoricamente da
situação em que as propriedades do fluido são constantes, o coeficiente h não varia ao
longo do tubo, e o perímetro molhado da conduta é constante. Nestas condições a
temperatura à saída do tubo é dada por:
Te = Tw − (Tw − Ti ) e p p
− A h / mc
ɺ

Se o parâmetro TU ≡ Ap h / mc
ɺ p (TU = número de unidades de transferência) for
grande (tipicamente maior do que 5) o valor de Te vem praticamente igual à temperatura
na parede Tw .
A temperatura média do fluido numa secção qualquer do tubo ( Ap ( x) = Px ) é dada por
uma expressão semelhante:
T ( x) = Tw − (Tw − Ti ) e
− ( Ph / mc
ɺ p )x
.
Estas fórmulas são válidas para o caso de aquecimento do fluido ( Tw > T , como na
figura acima), ou arrefecimento do fluido ( Tw < T ).

#o caso de qɺ w = Cte , a diferença de temperaturas entre a parede e o fluido é constante

ao longo da conduta:
qɺ w = h (Tw − T ) = Cte ⇒ (Tw − T ) = Cte T
o que implica que Tw e T aumentam linearmente: Twe
Tw
T ( x) = Ti + qɺw Px / mc
ɺ p Twi ∆T Te
e
T
Te = Ti + qɺw Ap / mc
ɺ p i

em que a área de permuta é dada pelo perímetro da secção da 0 L x

conduta multiplicada pelo seu comprimento, Ap = PL . As temperaturas na parede, à
entrada e saída duma zona de desenvolvimento completo, são dadas por:
Twi = Ti + qɺ w / h e Twe = Te + qɺw / h .

-6-
3.1 Convecção forçada dentro de tubos lisos

a. Laminar, em desenvolvimento térmico:

1/3 0.14
D   µb  −1 L 1
u = 1.86  RePr    Pr > 0.5 ; Re < 2300 ; Gz = < 0.05
L   µw  D Re Pr
(Sieder e Tate; propriedades a Tm = 0.5(Ti + Te ) )

b. Laminar desenvolvido, em condutas com secções diversas( Re = UDH / ν ):

Geometria Tw = Cte. qw = Cte. Factor atrito f
Cilindro u=3.66 u=4.36 64.00/Re
Hexágono u=3.35 u=4.00 60.20/Re
Quadrado u=2.98 u=3.61 56.92/Re
Rectângulo 2:1 u=3.39 u=4.12 62.20/Re
3:1 u=3.96 u=4.79 68.36/Re
4:1 u=4.44 u=5.33 72.92/Re
6:1 u=5.14 u=6.05 78.80/Re
8:1 u=5.60 u=6.49 82.32/Re
infinito u=7.54 u=8.24 96.00/Re
Elipse 2:1 u=3.74 u=4.56 67.28/Re
4:1 u=3.79 u=4.88 72.96/Re
8:1 u=3.72 u=5.09 76.60/Re
16:1 u=3.65 u=5.18 78.16/Re
Triângulo 10º u=1.61 u=2.45 50.80/Re
isósceles 30º u=2.26 u=2.91 52.28/Re
60º u=2.47 u=3.11 53.32/Re
90º u=2.34 u=2.98 52.60/Re
120º u=2.00 u=2.68 50.96/Re

c. Turbulento, desenvolvido: u = 0.023Re0.8 Pr n Re > 10 000 ; 0.7 ≤ Pr ≤ 160

com: n = 0.4 aquecimento; n = 0.3 arrefecimento

(Dittus e Boulter; propriedades a Tm = 0.5(Ti + Te ) )

Tubos não cilíndricos, diâmetro hidráulico: DH = 4 A / P e Re = ρ DH U / µ

com: A − área da secção transversal; P − perímetro molhado.

3.2 Tubos rugosos – Em escoamento laminar a rugosidade não tem influência na

transferência de calor ou no atrito. Em escoamento turbulento, a rugosidade da parede,
superfícies artificialmente corrugadas ou aplicação de alhetas internas, faz aumentar o
coeficiente convectivo h . Nesta situação, de regime turbulento, o número de Nusselt
deve obter-se do coeficiente de atrito, usando a analogia de Chilton-Colburn:
u
1/3
= 12 C f ( Re ≥ 4000 )
RePr
com C f dado abaixo (Moody ou Haaland).

3.3 Comprimentos de desenvolvimento, hidrodinâmico ( Lh ) e térmico ( Lt ):

Laminar: Lh / D = 0.05Re ; Lt / D = 0.05Re × Pr
Turbulento: Lh = Lt = 10 D

-7-
 L1
3.4 Coeficientes de atrito ( ∆p = f 2
ρU 2 ; f = 4C f ; C f = τ w / 12 ρU 2 )
D
O coeficiente de atrito C f não é mais do que a tensão de corte na parede ( τ w )
adimensionalizada com a energia cinética média do escoamento. Permite obter a perda
de pressão (perda de carga) que ocorre numa conduta ∆p = pi − pe , e assim calcular a
potência de bombagem necessária para mover o fluido: Wɺ = Vɺ ∆p = mɺ ∆p / ρ [W].
b

Tubo, laminar: f = 64 / Re ou C f = 16 / Re (Hagen-Poiseuille, teórico)

Tubo liso, turbulento: f = 0.184 / Re0.2 ; C f = 0.046 / Re0.2 (Moody; Re ≥ 4000 )
Tubo liso, turbulento: C f = 0.079 / Re1/ 4 (Blasius; 3000 ≤ Re ≤ 105 )
1  6.9  k / D 1.11 
Tubo rugoso, turbulento: = −3.6 log10  +   (Haaland)
Cf  Re  3.71  

Rugosidades absolutas (Moody; Massey)

Material k (mm)
Aço com rebites 1 – 10
Betão 0.3 – 3
Réguas madeira 0.2 – 1
Ferro fundido 0.25
Aço galvanizado 0.15
Ferro revestido a asfalto 0.12
Aço comercial ou ferro 0.045
forjado
Tubo extrudido 0.0015

-8-
Anexo (material avançado)

Dedução da diferença média logarítmica

Considerando que o fluido está a aquecer ao longo de uma conduta cuja parede está à
temperatura uniforme Tw , superior à temperatura média do fluido ( T = T ( x) ) numa
secção qualquer x , temos que a transferência de calor entre a
T
parede e o fluido numa porção de conduta entre x e x + dx
é: Tw temp. parede
. ∆Te
dQɺ = ( Pdx)h (Tw − T ) dQ
∆Ti Te
com P a representar o perímetro molhado da conduta e h o T(x)
coeficiente convectivo local. Este calor faz aumentar a Ti
entalpia (e a temperatura) do fluido:
dQɺ = mc
ɺ p dT 0 x x+dx L x
em que mɺ é o caudal mássico de fluido na conduta (constante, por conservação de
massa) e c p a capacidade térmica específica do fluido a pressão constante. Igualando
estas duas expressões para a taxa de calor diferencial, temos:
d (T − Tw ) Ph
ɺ p dT = Ph (Tw − T ) dx
mc ⇒ =− dx
(T − Tw ) mc
ɺ p
Designando T ′ = T − Tw e integrado entre a entrada da conduta ( x = 0 , T = Ti ) e a saída
da conduta ( x = L , T = Te ), temos:
e
dT ′
e
Ph  T −T  e
Ph Ph Ah
∫i T ′ = −∫i mc dx ⇒ ln  e w  = − ∫ dx = − L=− p
ɺ p  Ti − Tw  i
mc
ɺ p mc
ɺ p mc
ɺ p
Na integração assumiu-se que o perímetro é constante (para tubo redondo de diâmetro
D, P = π D ), e o coeficiente convectivo ou é também constante, ou será baseado num
L
1
valor médio, h ≡ hm = ∫ hdx . A área de permuta é Ap = PL . O resultado acima
L0
permite desde logo obter a temperatura do fluido à saída da conduta:
 T −T  Ah T −T  Ap h 
ln  e w  = − p ⇒ e w = exp  − 
 mc
 Ti − Tw  mc
ɺ p Ti − Tw  ɺ p
ou
 Ap h 
Te = Tw + (Ti − Tw ) exp  −  .
 mc
 ɺ p
Uma integração de 0 até x fornece ainda a temperatura numa secção qualquer da
conduta:
 Ap h x 
T ( x) = Tw + (Ti − Tw ) exp  −  .
 mcɺ L
 p 
A equação que dá a transmissão de calor global em toda a conduta, e define a diferença
média de temperaturas, é:
Qɺ = Ap h∆T
e o balanço global de energia (calor transferido para o fluido aumenta a sua entalpia) é:
ɺ p (Te − Ti )
Qɺ = mc

-9-
Dividindo uma pela outra, e usando o resultado acima para o logaritmo das diferenças
de temperatura entre a parede e o fluido na entrada e saída, temos:
mc
ɺ p (T − T )
∆T = (Te − Ti ) ⇒ ∆T = e i
Ap h  T −T 
ln  i w 
 Te − Tw 
finalmente, com ∆Ti = Tw − Ti , ∆Te = Tw − Te , ∆Ti − ∆Te = (Tw − Ti ) − (Tw − Te ) = Te − Ti
∆Ti − ∆Te
∆T = ∆Tln = .
 ∆Ti 
ln  
 ∆Te 

Esta expressão é válida para o caso em que o fluido está a ser arrefecido (neste caso
∆Ti = Ti − Tw e ∆Te = Te − Tw ) e continua a ser válida quando a temperatura da parede
também varia ( ∆Ti = Twi − Ti e ∆Te = Twe − Te ). Pode escrever-se trocando a ordem das
parcelas (ou seja, é indiferente saber se o fluido está a entrar ou sair por i ou e):
∆T − ∆Ti
∆Tln = e .
 ∆Te 
ln  
 ∆Ti 

-10-
Exemplo 1- Escoamento exterior sobre placa plana. Ar atmosférico a 20 ºC e 1 m/s
escoa-se sobre parede quadrada de 1 m x 1 m. A temperatura da parede é 40 ºC.
Calcular a taxa de calor transferido da parede para o ar.

A temperatura de filme (média entre temperatura da parede e do fluido), usada como

referência para o cálculo das propriedades do ar, é:

Tf = 1
2 (Tw + T∞ ) = 12 ( 40 + 20 ) = 30 ºC

e das tabelas para ar, obtém-se:

ρ = 1.164 kg/m3; k = 0.02588 W/m.K; µ = 1.872 ×10−5 kg/m.s;

ν = 1.608 ×10−5 m2/s; Pr = 0.7282 .

O número de Reynolds baseado no comprimento da placa (nesta caso, quadrada) é

U∞L 1× 1
Re = = = 62 179
ν 1.608 ×10 −5
o que corresponde a escoamento laminar na camada limite ( Re < 5 × 105 ). O coeficiente
de atrito médio, entre o bordo de ataque da placa e o seu comprimento, é

1.328 1.328
Cf = 1/ 2
= = 0.00532
Re 62 1791/2

pelo que a tensão de corte na parede vem

τw
Cf = ⇒ τ w = C f 12 ρU ∞2 = 0.00532 × (0.5 × 1.164 ×12 ) = 0.00310 N/m2
1
2 ρU ∞2

e a força de arrastamento, que o ar a 36 km/h aplica tangencialmente sobre a placa, será

dada pelo produto da tensão e a área superficial:

Fw = τ w A = 0.00310 × (1× 1) = 3.10 × 10 −3 N.

O número de Nusselt médio obtém-se da correlação:

u = 0.664 Re1/ 2 Pr1/3

que fornece

u = 0.664 × 62 1791/ 2 × 0.72821/3 = 149.0

O coeficiente convectivo calcula-se de:

u = hL / k ⇒ h = u × k / L = 3.855 W/m2.ºC

e a taxa de transferência de calor por convecção, da parede para o ar atmosférico,

-11-
 Qɺ = hA (Tw − T∞ ) = 3.855 × 1× (40 − 20) = 77 W.

Exemplo 2- Escoamento interior em conduta quadrada. Uma conduta de secção

quadrada, 15 cmx15 cm, e comprimento 10 m, transporta 0.10 m3/s de ar quente e serve
para aquecer o espaço envolvente num sótão. O ar é admitido na conduta a 85 ºC e a
temperatura das paredes da conduta é aproximadamente constante, a 70 ºC. Calcular o
calor transferido por unidade de tempo, a temperatura do ar à saída da conduta e a
potência do ventilador. (Cengel)

Assumindo uma temperatura média entre o ar e a parede de 80 ºC (para evitar fazer

interpolações nas tabelas), têm-se as seguintes propriedades físicas do ar:

ρ = 0.9994 kg/m3; k = 0.02953 W/m.K; c p = 1008 J/kg.ºC;

ν = 2.097 × 10−5 m2/s; Pr = 0.7154 .

O caudal volumétrico à entrada permite calcular a velocidade média:

Vɺ = UA m3/s ⇒ U = Vɺ / A = 0.10 / (0.15 × 0.15) = 4.444 m/s

e o número de Reynolds vem:

UDH 4.444 × 0.15

Re = = = 31 791
ν 2.097 × 10−5

correspondendo a escoamento turbulento na conduta ( Re > 4000 ). O diâmetro

hidráulico é dado pelo lado da secção quadrada:
A 0.15 × 0.15
DH = 4 = 4 = 0.15 m.
P 4 × 0.15
O caudal mássico, necessário adiante, é obtido multiplicando o caudal volumétrico pela
massa volúmica baseada na temperatura de entrada (onde o caudal foi medido):
mɺ = ρiVɺ = 0.9856 × 0.10 = 0.09856 kg/s
com ρi = 0.9856 kg/m3 a 85ºC.
Tratando-se de um processo de arrefecimento (o ar está mais quente do que as paredes
da conduta), a correlação de Dittus e Boulter é:

u = 0.023Re0.8 Pr 0.3 ⇒ u = 0.023 × 31 7910.8 × 0.71540.3 = 83.16

e o coeficiente convectivo vem:

h = u × k / L = 83.16 × 0.02953 / 0.15 = 16.37 W/m2.ºC

Verifica-se que se trata de um valor superior aos habitualmente encontrados com

convecção natural, como seria esperado. A taxa de transferência de calor por convecção
é dada por:

Qɺ = hAp ∆Tln

-12-
necessitando do cálculo da temperatura média logarítmica:

∆Ti − ∆Te
∆Tln =
 ∆T 
ln  i 
 ∆Te 
com ∆Ti = Ti − Tw = 85 − 70 = 15 ºC e ∆Te = Te − Tw .

A temperatura média do ar à saída da conduta obtém-se de:

Te = Tw − (Tw − Ti )e− TU = 70 − (70 − 85)e−0.9888 = 75.58 ºC

em que o número de unidades de transferência, para uma área de permuta de
Ap = 4aL = 6 m2, é definido como
TU = Ap h / mc
ɺ p = 6 × 16.37 /(0.09856 × 1008) = 0.9888 .

Assim, ∆Te = 75.58 − 70 = 5.58 ºC, e

15 − 5.58
∆Tln = = 9.526 ºC
 15 
ln  
 5.58 
pelo que o taxa de transferência de calor é:
Qɺ = hAp ∆Tln = 16.37 × 6 × 9.526 = 935.6 W
ou
ɺ p (Ti − Te ) = 0.09856 × 1008 × (85 − 75.58) = 935.8 W.
Qɺ = mc

Podia-se voltar a fazer os cálculos com a propriedades obtidas à temperatura média

global Tm = 0.5(Ti + Te ) = 0.5(85 + 75.58) = 80.3 ºC mas os resultados iriam ser
praticamente idênticos.

Cálculo da perda de carga: o coeficiente de atrito, da expressão de Blasius, é dado por:

C f = 0.079 / Re1/4 = 0.079 / 317910.25 = 0.005916
pelo que a tensão média na parede interior da conduta vem:
τ w = 12 ρU 2C f = 0.5 × 0.9994 × 4.4442 × 0.005916 = 0.05840 N/m2.
A perda de pressão correspondente é dada por
L 1 10
∆p = 4C f 2
ρU 2 = 4 × 0.005916 × × 0.5 × 0.9994 × 4.444 2 = 15.6 N/m2
DH 0.15
que também poderia ser obtida do balanço de forças global: A∆p = Apτ w . O ventilador
usado para mover o ar quente terá de produzir um aumento de pressão superior a esta
perda de carga. A sua potência é calculada de:

WɺVentilador = Vɺ ∆p = 0.1× 15.6 = 1.6 W.

-13-