Trabalho de Fenômenos - Pesquisa a acadêmica sobre calor, condução, radiação e convecção

Prévia do material em texto

UNIVERSIDADE DE PERNAMBUCO 
Curso de Graduação em Engenharia Elétrica 
 
 
 
 
 
 
 
 
Thiago Brito Cassimiro da Silva 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
Pesquisa a acadêmica sobre calor, condução, radiação e convecção 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
Recife 
2018 
 
Sumário 
1. DEFINIÇÕES ............................................................................................................. 3 
1.1. Calor; .................................................................................................................. 3 
1.2. Condução; .......................................................................................................... 3 
1.3. Radiação ............................................................................................................. 5 
1.4. Convecção .......................................................................................................... 6 
2. DEDUÇÃO DAS EQUAÇÕES GERAIS ............................................................... 6 
2.1. Transmissão de Calor por Condução ................................................................. 6 
2.1.1. Coordenadas Cartesianas: .................................................................................. 6 
2.1.2. Coordenadas Cilíndricas: ................................................................................... 9 
2.2. De Transmissão de Calor por Radiação ........................................................... 12 
2.3. Transmissão de Calor por Convecção .............................................................. 14 
2.3.1. Equação da Continuidade: ................................................................................ 14 
3. BIBLIOGRAFIA ...................................................................................................... 32 
 
 
 
 
 
 
 
1. DEFINIÇÕES 
 
1.1. Calor; 
 
Segundo definição do dicionário Aurélio, calor é “... a forma de energia que se transfere de um 
sistema para o outro em virtude duma diferença de temperatura existente entre os dois, e que se 
distingue das outras formas de energia porque, como o trabalho, só se manifesta num processo 
de transformação”.1 Desta definição de calor pode-se observar que o calor nada mais é do que 
uma forma de energia. A maneira com que esta energia altera as propriedades (dependentes ou 
independentes) de um sistema no estado de equilíbrio é escopo do estudo da Termodinâmica 
Clássica. Já os efeitos que ocorrem durante o processo da transmissão da energia em forma de 
calor é escopo da Transferência de Calor. Dentre estes efeitos, exalta-se a variação da taxa 
temporal de transmissão de calor (variável não considerada na termodinâmica). Os principais 
princípios da Termodinâmica devem estar bem assimilados a fim de se absorver o conteúdo 
deste trabalho da forma mais efetiva. 
 
 
1.2. Condução; 
 
“A condução é um processo pelo qual o calor flui de uma região de temperatura mais alta para 
outra de temperatura mais baixa, dentro de um meio (sólido, líquido ou gasoso) ou entre meios 
diferentes em contato físico direto”.4 Essa explicação abrange tanto a apresentação da Segunda 
Lei da Termodinâmica, quando se diz que a transmissão de calor parte de uma região de 
temperatura mais alta para outra de temperatura mais baixa, quanto a definição específica do 
processo de transmissão de calor por Condução: “...em contato físico direto.”. Portanto, a 
transmissão de calor por Condução ocorre quando corpos em diferentes temperaturas estão 
literalmente “encostados” um no outro. A energia (calor) do corpo de temperatura mais alta 
agita as moléculas do corpo de temperatura mais baixa, fazendo com que a energia cinética 
 
média das moléculas deste último se eleve, aumentando, assim, sua energia interna. 
Conseqüentemente, a temperatura do corpo que está “recebendo” a energia em forma de calor 
se eleva até o estado de equilíbrio. Para ilustrar este fenômeno, imagina-se um bule com água 
fervendo ao fogão. O fogo aquece o bule, o qual, por condução, aquece a parcela de líquido que 
está em contato direto com o mesmo. Esta transmissão de calor por condução é a única maneira 
de que o calor pode ser transmitido entre corpos sólidos opacos. Já em meios líquidos, a 
condução também apresenta grande importância, embora esteja, quase sempre, relacionada com 
outros meios de transmissão de calor. 
Segundo definição do cientista francês J.B.J. Fourier, em 1882, a quantidade de calor 
transmitida por condução segue a seguinte lei: 
 
 
q  kA 
dT
 
k 
dx 
 
 
 
 
 
Na formulação acima, k representa a condutividade térmica do material, A representa a área da 
seção através da qual o calor flui por condução (medida perpendicularmente à direção do 
fluxo), e dT/dx representa o gradiente de temperatura na seção. Nesta formulação, toma-se 
como convenção a direção de aumento na coordenada x como fluxo de calor. Sabendo-se que, 
pela segunda lei da termodinâmica, o calor flui da região de maior temperatura para a região de 
menor temperatura, deve-se adotar o sinal negativo para o produto acima, conforme mostra a 
equação. 
Como se pode observar pelo positivo balanço de unidades da fórmula acima, qk é medido em 
quantidade de calor por unidade de tempo. Usualmente, esta grandeza é expressa em 
kilocalorias por hora, ou kcal/h. 
 
1.3. Radiação 
“A radiação é um processo pelo qual o calor é transmitido de um corpo a alta temperatura para 
um de mais baixa quando tais corpos estão separados no espaço, ainda que exista vácuo entre 
eles”. Por esta definição, vê-se que não há necessidade de um contato físico entre os corpos 
para que a energia (na forma de calor) seja transmitida entre eles. Ao calor transmitido desta 
forma dá-se o nome de calor radiante. Esta forma de energia se assemelha 
fenomenologicamente à radiação da luz, diferindo-se apenas nos comprimentos de onda. A 
transmissão do calor radiante ocorre na forma de quanta (porções discretamente definidas) de 
energia. Para ilustrar este fenômeno, tenta-se imaginar como a energia solar é transmitida para 
os demais astros (a Terra, por exemplo), sendo que há praticamente vácuo entre eles. Nas nossas 
aplicações práticas, o calor radiante tem vital importância no projeto de, por exemplo, uma 
caldeira. Além da energia que é transmitida do combustível queimado para as paredes da 
caldeira, existe também uma parcela de calor radiante que é transmitido para seus demais 
componentes. 
 
 
 
1.4. Convecção 
“A convecção é o processo de transporte de energia pela ação combinada da condução de calor, 
armazenamento de energia e movimento de mistura. A convecção é importante principalmente 
como mecanismo de transferência de energia entre uma superfície sólida e um líquido ou gás”. 
Em um fluido, onde a mobilidade das partículas é grande, as partículas aquecidas pelo contato 
direto com a superfície sólida tendem a migrar para locais onde as temperaturas são mais baixas. 
Esta movimentação de partículas acarreta uma transferência de energia de uma posição para a 
outra, caracterizando a transmissão de calor por convecção. 
 
 
2. DEDUÇÃO DAS EQUAÇÕES GERAIS 
 
 
2.1. Transmissão de Calor por Condução 
 
Aqui serão deduzidas analiticamente as fórmulas da Transmissão de Calor por Condução nas 
coordenadas Cartesianas e Cilíndricas: 
 
2.1.1. Coordenadas Cartesianas: 
 
Considere-se um pequeno elemento material de um corpo sólido, na forma de paralelepípedo, 
de volume dxdydz, conforme mostra a Figura 1. 
 
 
 
 
 
 
 
 
Figura 1: Elemento materialde 
corpo sólido, na forma de 
paralelepípedo, a partir do qual 
serão deduzidas as equações da 
transmissão de calor por 
condução para coordenadas 
cartesianas. 
 
 
 
 
 
 
 
 
 
O balanço de energia deste elemento pode ser calculado da seguinte forma: 
 
Calor que entra no 
elemento por 
unidade de tempo 
 
 
+ 
Calor gerado no 
elemento por 
unidade de tempo 
 
 
= 
Calor que Sai do 
elemento por 
unidade de tempo 
 
 
+ 
Variação da 
energia 
interna com 
o tempo 
 
 
 
 
Algebricamente, esta equação pode ser escrita da seguinte forma: 
 
(qx  qy  qz )  
dq 
 (dx  dy  dz)  (q 
dt 
 
 
xdx  qy dy  qz dz 
)  c (dx  dy  dz)  
T
 


 
 
Na equação acima,
 dq
 
 
é o calor gerado por unidade de tempo e por unidade de 
dt 
volume, T é a temperatura do corpo, θ é o tempo, c é o calor específico do material 
e ρ é sua densidade. 
 Aplicando a equação C.2 a este problema, o calor transmitido por condução 
para dentro do corpo na direção x, mostrado na Figura 1, pode ser escrito como 
segue: 
 
q  (k 
T 
)  dy  dz 
x 
x 
 
Aplicando esta mesma equação C.2 ao calor que deixa o corpo, por condução, na direção x, 
pode-se escrever: 
qxdx 
 [(k 
T 
) 
x 
 
(k 
T 
)dx]d 
x x 
y d z Eq. D.a1.4 
 
Fazendo a subtração destes dois termos, temos: 
 
 
(qx 
 
)  (q 
 
 
xdx 
)  [ 


x 
(k 
T 
)]d d 
x 
x
 
y d z 
Eq. D.a1.5 
 
Aplicando, de forma análoga, a mesma equação C.b.1 para as direções y e z, e aplicando os 
resultados assim obtidos na equação D.a1.2, podemos facilmente chegar à seguinte equação: 
 
(k 
T 
) 
 
 
(k 
T 
) 
 
 
(k 
T 
)  
dq 
 c 
T 
 
Eq. D.a1.6 
x x y y z z dt 
 

A equação D.a1.6 é a equação geral, obtida de forma analítica, para o problema de transmissão 
de calor por condução em três dimensões. Observe que esta equação pode ser bastante 
simplificada se assumido que a condutividade térmica do material, k, é uniforme, e se 
assumido que o calor específico, c, e a densidade, ρ, forem independentes da temperatura. A 
equação resultante, nestas condições, é chamada de equação de Fourier. Se, além disto, for 
admitido regime permanente, o termo da direita se anula. Adicionando- se esta condição à 
equação de Fourier, chega-se na equação de Poisson. Se, ainda mais, não for considerada a 
geração interna de calor, a equação D.a1.6 assume sua forma mais simplificada possível, à 
qual se da o nome de equação de Laplace: 
 2T 
 
 2T 
 
 2T 



x 2 y 2 z 2 = 
0 
 
Para o problema unidimensional apresentado na secção C.a, sendo a segunda derivada da 
temperatura em relação à posição igual a zero, pode-se dizer que a primeira derivada da 
temperatura em relação à posição é constante.
Portanto: 
dT 
 const., 
dx 
 
conforme sugerido na equação C.b.1. 
 
 
 
 
2.1.2. Coordenadas Cilíndricas: 
Considere-se o sistema de referência em coordenadas Cilíndricas apresentado na Figura 2 
 
 
 
 
 
 
 
 
 
 
 
Figura 2: Representação de 
um ponto no espaço a partir 
das coordenadas cilíndricas. 
 
 
 
Para este sistema de coordenadas, valem as seguintes relações, comparando-se com o sistema 
de coordenadas cartesiano: 
 
 
 
x  r cos

y  r sin 


Eq. D.a2.1 
 
Eq. D.a2.2
z  z Eq. D.a2.3 
 
Portanto, sendo a temperatura T = T(x,y,z) para as coordenadas cartesianas, valerá a 
expressão T = T(rcosΦ, rsinΦ, z) para as coordenadas cilíndricas mostradas na Figura 2. 
Sendo assim, pode-se dizer que os termos independentes das coordenadas x e y apresentados 
na Eq. D.a1.6 são mantidos da mesma maneira para coordenadas cilíndricas e 
cartesianas. Já os dois primeiros termos desta mesma expressão, 
 2T 
k 
x 2 
 2T 
e k 
y 2 
, devem 
 
ser trabalhados. Como a constante k é independente do sistema de coordenadas, ela será 
desconsiderada aqui. 
Inicialmente, será derivada a temperatura em relação a r, no sistema de 
coordenadas cilíndricas. Aplicando-se a regra da cadeia a este problema, tem-se: 
 
T 
 ( 
T 
)( 
x 
)  ( 
T 
)( 
y 
)  ( 
T z 
 
 
Eq. D.a2.4 
r x r y r z 
)( 
r 
)
 
 
Sendo 
x 
 cos , 
r 
y 
 sin 
r 
e 
z 
 0 , tem-se o que segue: 
r 
 
T 
 ( 
T 
) cos  ( 
T 
) sin 
  

Eq. D.a2.5 
r x y 
 
Aplicando a segunda derivada de T em função de r, tem-se: 
 
 2T 
r 2 
 
 
[( 
T 
r x 
) cos  ( 
T
 
y 
) sin  ] 
 
Eq. D.a2.6 
 
Portanto, 
 
 2T 
r 2 
 
 2T 
x 2 
 
cos
2
   2 
 2T 
 
 
xy 
(cos sin ) 
 2T 
y 2 
sin 
2 

Eq. D.a2.7 
Esta equação será guardada, por enquanto, da maneira em que está apresentada. Agora, 
realizando o mesmo procedimento para a derivada de T em relação a Φ, 
tem-se: 
 
T 
 ( 
T 
)( 
x 
)  ( 
T 
)( 
y 
)  ( 
T 
)( 
z 
) 
 
 
 
Eq. D.a2.8 
 x  y  z 


Sendo 
x 
 r sin  , 

y 
 r cos , e 

z 
 0 , tem-se o que segue: 

 
 
T 
 ( 
T 
)(r sin )  ( 
T 
)(r cos) 
 
 
Eq. D.a2.9 
 x y 
 
Agora, fazendo-se a segunda derivada de T em relação a  , tem-se: 
 
 
 2T 
 2 
 
 [( 
 2T 
x 2 
 
)(r 2 sen
2)  ( 
 2T 
xy 
 
)(r 2 sen cos)  ( 
T
 
x 
)(r cos)] 


Eq. D.a2.10 
 2T 
[( 
y 2
 
)(r 
2
 cos
2
 )  ( 
 2T 
xy 
)(r 2 sen cos)  ( 
T
 
y 
)(rsen)] 
 
Somando-se os dois termos da direita da Eq. D.a2.10, tem-se (ver Eq. D.a2.5): 
 
( 
T 
)(r cos)  ( 
T 
)(rsen)  r 
T
 
 
 
Eq. D.a2.11 
x y r 
 
Com a Eq. D.a2.11 e fazendo as devidas contas, pode-se reescrever a Eq. D.a2.10 como 
segue: 
 
 2T 
 2 
 2T 
( 
x 2 
 
)(r 2 sen
2)  2  ( 
 2T 
xy 
 
)(r 2 sen cos)  ( 
 2T 
y 2 
 
)(r 
2
 
 
cos
2
 
)  r 
T
 
r 
 
Eq. D.a2.12 
 
Dividindo-se a Eq. D.a2.12 por r2, e somando-a à Eq. D.a2.7, tem-se o que segue: 
 
1  2T 
 
 2T 
 
 2T 
( 
 
)  ( 
 2T 
) 
1 T 
 
 
 
Eq. D.a2.13 
r 
2  2 r 2 x 2 y 2 r r 
 
Então, pode-se reescrever a Eq. D.a2.13 como segue: 
 
 2T 
 
 2T 
 
1  2T 
 
 2T 
 
1 T 
 
x 2 y 2 r 2  2 r 2 r r Eq. D.a2.14 
 
Aplicando-se a Eq. D.a2.14 na Eq. D.a1.6, pode-se chegar na equação geral de transmissão 
de calor por condução em coordenadas cilíndricas, como segue: 
 2T 
k 
r 2 
 k(
1
) 
T
 
r r 
 k( 
1
 
r 2 
 2T 
) 
 2 
 2T 
k 
z 2 
 dq 
dt 
 c 
T
 

Eq. D.a2.15 
 
Esta equação coincide com aquela encontrada na literatura. 
 



 
 
2.2. De Transmissão de Calor por Radiação 
 
Como citado anteriormente, a natureza da maneira com a qual o calor radiante se transmite de 
um corpo para outro ainda não é totalmente conhecida. Existem hipóteses que dizem que o 
calor radiante é transmitido da mesma forma das ondas eletromagnéticas, porém existem 
outras hipóteses que dizem que o calor radiante é transmitido por fótons. 
Cada corpo, dependendo de suacomposição química e de outros parâmetros, possui um certo 
poder de radiação. Como já citado, os corpos negros, ou irradiadores perfeitos, são aqueles 
corpos capazes de emitir a máxima quantidade de radiação para uma dada temperatura. Este 
conceito é teórico, e estabelece um limite máximo para a radiação de calor, conforme Segunda 
Lei da Termodinâmica. A quantidade de energia de radiação por unidade de tempo e por 
unidade de área transmitida por um corpo é chamada de poder emissivo espectral (ou 
monocromático). Este termo sugere que a quantidade de energia transmitida por um corpo na 
forma de radiação térmica varia com o comprimento de onda. De fato, a relação que fornece 
o poder emissivo espectral de um corpo negro é a seguinte: 
 
A 
B 
 
Eb (T )  
5 
(e T  1) Eq. D.b.1 
 
 
 
Na equação acima, Ebλ representa o poder emissivo monocromático de um corpo negro, como 
função da temperatura, em kcal/hm2µm, A é uma constante que vale 3,22 × 108 kcalµm4/hm2, 
B é uma constante que vale 1,4388 µmK, λ é o comprimento de onda emitido pelo corpo, em 
µm, e T é a temperatura absoluta do corpo, em Kelvin. 
Como apresentado na secção C.c.1, o calor total por unidade de tempo emitido por radiação 
por um corpo pode ser calculado da seguinte forma: 
E (T )  
qr 
b 
A
  T 
4
 Eq. D.b.2 
 
O produto, σT4, desta equação pode ser calculado a partir da integração da equação D.b.1 em dλ de 0 a 
∞. 
Para as aplicações na engenharia, é comum precisar calcular a quantidade de calor radiante 
emitida por um corpo em uma determinada faixa de comprimentos de onda. Para tanto, basta 
integrar a equação D.b.1 em dλ, definida entre o intervalo de comprimentos de onda com o qual 
se deseja trabalhar. 
Existe um valor de comprimento de onda para a qual o valor de Ebλ(T) é máximo, em uma dada 
temperatura. A relação entre este comprimento de onda, batizado de λmax, e a temperatura do 
corpo pode ser dada pela equação empírica que segue: 
max  
2898 
T 
Eq. D.b.3 
 
Nesta equação, λ aparece em µm e T em K. 
 
ara facilitar a resolução dos problemas que envolvem transmissão de 
calor por radiação, geralmente baseia-se nos seguintes gráficos: 
1) Ebλmax × λ,T; 
 
2) Ebλmax × T; 
 
3) Ebλ/Ebλmax × λT; 
 
 
 
2.3. Transmissão de Calor por Convecção 
A Transmissão de Calor por Convecção não é puramente um modo de transmissão de calor, 
mas sim uma combinação de transmissão de calor com movimentos mecânicos de fluidos. 
Assim, a fim de se entender mais a fundo os conceitos abordados nesta seção, é importante 
conhecer bem alguns aspectos abordados em mecânica dos fluidos. Em particular, os conceitos 
de escoamento em regime laminar e turbulento, assim como os conceitos de camada limite e 
separação da mesma, têm vital importância no entendimento de Transmissão de Calor por 
Convecção. Os aspectos aqui abordados foram obtidos do livro de Fox e McDonald, p. 12 a 28, 
e 326 a 347. 
A equação C.d.1 apresentada anteriormente é frequentemente usada pelos engenheiros para o 
cálculo da Transmissão de Calor por Convecção, porém ela nada mais é do que a definição do 
coeficiente de transmissão de calor por convecção, hc. A fim de se entender melhor a 
transmissão de calor por convecção, serão aqui abordados os métodos de avaliação do 
coeficiente de transmissão de calor por convecção e será explicado o papel das aletas neste 
método de transmissão de calor. 
 
2.3.1. Equação da Continuidade: 
Segue abaixo a Equação da Continuidade: 
 
 
u 
  
v 
  
w 
  

   
 


Eq. D.c.1.ii.1 
x y z t 
 
Esta equação também é conhecida como equação da conservação da massa. Nela, sugere-se que 
o gradiente de velocidade multiplicado pela massa especifica do fluido se iguala à variação da 
massa especifica com relação ao tempo. Para fluidos em escoamento permanente (i.e. ρ = 
ρ(x,y,z)), o termo do lado direito desta formulação se anula. Para fluidos incompressíveis (i.e. 
ρ(x,y,z) = constante), pode-se simplificar esta equação ainda mais, restando apenas que o 
gradiente da velocidade do fluido é nulo. Se o escoamento, assumindo todas as hipóteses acima, 
é bidimensional, esta formulação pode ser rescrita 
como segue: 
 
u 
 
v 
 0
 
 
 
 
Eq. D.c.1.ii.2 
x y 
 
 
Esta situação é bem simplificada, porém sua análise pode ser extrapolada para inúmeras 
situações reais. Sendo assim, para se obter a solução matemática da camada limite, 
consideraremos fluidos escoando nas condições citadas acima. 
a) Equação da Força de Cisalhamento: 
 
A definição de um fluido Newtoniano segue abaixo: 
 
   
du
 
dy 
 
 
Eq. D.c.1.ii.3 
 
Esta relação sugere que a força de cisalhamento causada pelo escoamento de um fluido é diretamente 
proporcional à derivada da sua velocidade em relação à posição perpendicular ao escoamento. A 
constante de proporcionalidade, µ, é uma propriedade do fluido, chamada viscosidade dinâmica, de 
unidades comumente usadas 
 
kgf  s 
.
 
m2 
 
 
 
b) Equação da Quantidade de Movimento aplicada a uma Placa Plana: 
 
Aplicando-se as relações acima sobre uma placa plana e usando-se um pouco de álgebra, 
pode-se facilmente chegar na seguinte relação: 
 (u 
u
 
x 
 
v 
u 
y 
)  
 2u 
 
 
y 2 
 
p
x Eq. D.c.1.ii.4 
 
Onde p é a pressão exercida pelo fluido. 
 
Resolvendo-se esta equação simultaneamente à equação da continuidade, pode-se 
obter a espessura da camada limite e a força de atrito na parede da placa plana. Para tal, 
defini-se a função de corrente. 
 
c) Função de Corrente: 
 
A função de corrente, ψ (x,y) é definida de forma que sejam válidas as duas relações 
abaixo: 
u  
 
Eq. D.c.1.ii.5 
y 
 
 
v   
 
Eq. D.c.1.ii.6 
x 
 
Observe que esta equação satisfaz automaticamente a equação da continuidade. 
 
Para dar sequência ao estudo, introduz-se uma nova variável, η, como segue: 
 
  y 
 
Eq. D.c.1.ii.7 
 
 
Então, pode-se dizer o seguinte: 
 
  f () 
 
Eq. D.c.1.ii.8 
 
 
Onde f(η) designa uma função de corrente adimensional. 
 
Obtendo-se as componentes de velocidade em função de f(η), e exprimindo ∂u/∂x, 
 
∂u/∂y, ∂2u/∂y2, podemos facilmente chegar à equação diferencial ordinária, de terceira 
ordem e não-linear, como segue: 
f () 
d 
2 
[ f ()] 
 
 
d 2 
d 
3 
[ f ()] 
2 
d 3 
 0 Eq. D.c.1.ii.9 
 
Sujeita a três condições de contorno em η = 0 e η = ∞. 
Após análise dos fatores podemos concluir que o coeficiente de Transmissão de Calor por 
Convecção é: 
h  
q 
(T  T )  0,332 
k 
Re 


Eq. D.c.1.ii.21 
cx 
A 
 s 
x 
L 
 
E assim extrai-se o valor do coeficiente de Convecção. Na prática, as propriedades físicas 
usadas na formulação acima variam com a Temperatura. Entretanto, nesta análise, as 
propriedades foram tomadas como contantes. Os resultados da Eq. D.c.1.ii.21 são 
satisfatórios se considerado o valor da temperatura média entre a temperatura da parede 
sólida e do fluido ao longe. 
3 Pr 
(
u ) 
vx 
u 
vx 

 
3. BIBLIOGRAFIA 
 
Chapman, Alan J., Fundamentals of Heat Transfer, 1984, Macmillan Publishing Company, 
New York. 
Ferreira, Aurélio Buarque de Holanda, Novo Dicionário Aurélio, Editora Nova Fronteira, Rio 
de Janeiro, Brasil. 
Fox, Robert W. e Alan T. McDonald, Introdução à Mecânica dos Fluidos, 1998, LTC, Rio de 
Janeiro, Brasil. 
Harrington, R. L, Marine Engineering, 1992, The Society of Naval Arquitecture and Marine 
engineering, Jersey City. 
Kreith, F. eBohn, MS. Princípios de Transferência de Calor, 2003, Editora Edgard Blücher, 
São Paulo. 
Kreith, F. e Bohn, MS. Princípios de Transferência de Calor, 1977, Editora Edgard Blücher, 
São Paulo. 
Li, K. W., e A. P. Priddy, Power Plant System Design, 1985. John Wiley & Sons Inc., New 
York City, USA. 
Moran, Michael J. e Howard N. Shapiro, Fundamentals of
 Engineering Thermodynamics, 2000, John Wiley & Sons Inc., New 
York City, USA. 
Pohlhausen, E., Der Wärmeaustausch swischen festen Körpern und Flüssigkeiten mit kleiner 
Reibung und Kleiner Wärmeleitung, 1921, ZAMM, Vol. 1. 
Stewart, James, Cálculo, Vol. II, 1999, Ed. Pioneira, 4ª. Edição, São Paulo, Brasil.