13º ENCIT 13º ENCIT – Congresso Brasileiro de ... - LEPTEN - UFSC

13º ENCIT – Congresso Brasileiro de Engenharia
e Ciências Térmicas
UFU UFU-Uberlândia
Uberlândia - MG, 05 a 09 de dezembro de 2010
FUNDAMENTOS E APLICAÇÕES DE
TRANSFERÊNCIA DE CALOR COM MUDANÇA DE
FASE: EBULIÇÃO E CONDENSAÇÃO
CONDENSAÇÃO
(4ª aula)
Júlio César Passos
jpassos@emc.ufsc.br
http://energetique-juliocesarpassos.blogspot.com
Universidade Federal de Santa Catarina
Centro Tecnológico - Departamento de Engenharia Mecânica
LEPTEN/Boiling
Laboratórios de Engenharia de Processos de Conversão e Tecnologia de Energia

• Aplicações
Plano da aula
• Aspectos históricos sobre a condensação
• Tipos de condensação
• Condensação em película – modelo de
Nusselt
• Exemplo resolvido
• Condensação em microcanais (resultados
preliminares, obtidos no LEPETEN/Boiling – UFSC)
13º ENCIT Congresso Brasileiro de Engenharia e Ciências Térmicas
UFU-Uberlândia– 05-09/12/2010

Aplicações
• Condensadores com e sem superfícies
microaletadas
indústria de refrigeração,
ciclos de potência (geração
termelétrica)
• Tubos de calor
• Termossifões bifásicos
coletores solares
fornos de cocção
13º ENCIT Congresso Brasileiro de Engenharia e Ciências Térmicas
UFU-Uberlândia– 05-09/12/2010

Aspectos históricos da condensação
1769 - James Watt obtém a patente do
condensador de superfície;
(147 anos)
1916 - Modelo de Nusselt para a condensação em película
(in Nusselt, W., “Die oberflachenkondesation des wasser
dampfes”, Z. ver., 1916)
1930 - Schimidt et al. apresentaram um dos primeiros
trabalhos sobre a condensação em gotas.
13º ENCIT Congresso Brasileiro de Engenharia e Ciências Térmicas
UFU-Uberlândia– 05-09/12/2010

g
q″
TIPOS DE CONDENSAÇÃO
vapor
Líquido
g
q″
vapor
em película em gotas mista
q″
g
vapor
Líquido Líquido
hh gotas > > hpelícula 13º ENCIT Congresso Brasileiro de Engenharia e Ciências Térmicas
UFU-Uberlândia– 05-09/12/2010

TIPOS DE CONDENSAÇÃO
• O tipo de condensação depende do ângulo de contato
líquido líquido-superfície.
superfície.
• Dois tipos:
– Condensação em gotas – maior transferência de
calor, porém sua manutenção é muito difícil.
Ocorre com fluidos não molhantes.
– Condensação em película – modo mais comum de
condensação. Formação de uma película líquida
sobre a superfície. Ocorre com fluidos molhantes.
13º ENCIT Congresso Brasileiro de Engenharia e Ciências Térmicas
UFU-Uberlândia– 05-09/12/2010

CONDENSAÇÃO
EM GOTAS (1)
Condensação sobre uma
superfície vertical
(q=0,12 W/cm W/cm2 ; ΔT=0,05 T=0,05 K)
( h= 24 kW/K.m 2 )
a) início
b) 16 s
c) 32 s
Diâmetro da superfície superfície: : 18mm
in Stephan, K., “Heat Transfer in
Condensation and Boiling”, Springer, 1992, p. 29
13º ENCIT Congresso Brasileiro de Engenharia e Ciências Térmicas
UFU-Uberlândia– 05-09/12/2010

CONDENSAÇÃO EM GOTAS (2)
Condensação sobre uma
superfície vertical
(in J. H. Leinhard, “A Heat Transfer Textbook”,
2nd ed., Prentice-Hall, 1987)
13º ENCIT Congresso Brasileiro de Engenharia e Ciências Térmicas
UFU-Uberlândia– 05-09/12/2010

Codensação em gotas (3)
Condensação sobre uma
superfície coberta
com teflon
(in in Ganzeles, Ganzeles, F.L.A., F.L.A., 2002, 2002,
“Drainage “Drainage and and condensate condensate heat heat
resistence resistence in in dropwise dropwise condensation”,
condensation”,
PhD PhD thesis thesis TuE TuE--Eindhoven Eindhoven--Holanda Holanda)
13º ENCIT Congresso Brasileiro de Engenharia e Ciências Térmicas
UFU-Uberlândia– 05-09/12/2010

Padrões na codensação em gotas (3)
(in Ganzeles, F.L.A., 2002 2002, ,
PhD thesis
TuE TuE-Holanda) Holanda)
13º ENCIT Congresso Brasileiro de Engenharia e Ciências Térmicas
UFU-Uberlândia– 05-09/12/2010

Curva de
Condensação
Comparação:
Condensação
em gotas e
em película
in Massoud Kaviany,
“Principles of Heat Transfer“,
John Wiley & Sons, 2002, p. 610.
13º ENCIT Congresso Brasileiro de Engenharia e Ciências Térmicas
UFU-Uberlândia– 05-09/12/2010

Ângulo de contato e Molhabilidade
(Ângulo de contato estático)
σ LV
SV
Equação de Young
σ = σ + σ cosθ
SV
SL
LV
σ −σ
= σ cosθ
SV
SL
LV
θ
L
σσ σσ
SL
Equilíbrio de forças na linha tripla
13º ENCIT Congresso Brasileiro de Engenharia e Ciências Térmicas
UFU-Uberlândia– 05-09/12/2010
V
S

σ LVII
L
σσ
Força de adesão, F Fad ad
V
SL
S
σ SV
σ LVII
( σ + σ + σ ) − ( σ + )
F =
σ
ad LV LVII
SV LVII
( σ + − )
F = σ σ
ad
LV
SV
SL
V
L
SL
V
σ
LV
S
13º ENCIT Congresso Brasileiro de Engenharia e Ciências Térmicas
UFU-Uberlândia– 05-09/12/2010

σ LVII
c
L
σσ
Força de Coesão, F Fc V
SL
F = 2σ
LV
S
σ LVII
σ V
SL
( σ + σ + σ ) − ( σ + )
F = 2
σ
c LV LVII
SL LVII
L
L
V
SL
σ LV
13º ENCIT Congresso Brasileiro de Engenharia e Ciências Térmicas
UFU-Uberlândia– 05-09/12/2010
S

da Eq. de Young
Força de Espalhamento, F Fe F = F − F
F
F
e
e
=
e
a d
( σ + σ −σ
) − ( σ )
= 2
LV
( σ −σ
−σ
)
SV
SV
LV
c
SL
SL
σ −σ
= σ cosθ
SV
( 1− θ )
SL
LV
Fe = −σ
SL cos = 0 , θ = 0
F < 0, 0 < θ < 90°
Líquido molhante
e
F Líquido totalmente molhante
90° < θ < 180°
Líquido não molhante
= −2σ
LV θ = 180°
e
F Líquido totalmente não molhante
e
LV

Alteração do ângulo de contato
Experimento no LMPT e LEPTEN/UFSC
13º ENCIT Congresso Brasileiro de Engenharia e Ciências Térmicas
UFU-Uberlândia– 05-09/12/2010

O MODELO DE NUSSELT
• Escoamento laminar;
• Propriedades físicas do fluido constantes;
• O sub sub-resfriamento resfriamento do líquido na Camada Limite é
desprezível, no balanço de energia;
• A equação da QM não considera os termos de inércia;
• Vapor em repouso, sem exercer força de arrasto ou
de cisalhamento sobre o líquido;
• A interface líquido líquido-vapor vapor é lisa;
• A transferência de calor, através da película, ocorre
por condução.
13º ENCIT Congresso Brasileiro de Engenharia e Ciências Térmicas
UFU-Uberlândia– 05-09/12/2010

Nusselt (1882 (1882-1957) 1957)
Ernst Kraf Vilhelm Nusselt
Publicou, em 1916,
“Die Oberflachenkondensation des
Wasserdampfes”, Z. Ver. Deut. Ing., Vol.
60, 541, 1916.
13º ENCIT Congresso Brasileiro de Engenharia e Ciências Térmicas
UFU-Uberlândia– 05-09/12/2010

O MODELO DE NUSSELT,
y
Laminar sem
ondulações
Re
δ
≤
30
Laminar com
ondulações
Re
30 ≤ Reδ
≤
Região turbulenta
δ
≥ 1800
1800
2
d u
2
dy
=
Reδ
−
1
µ
l
≤
30
y=0, u=0 ,
( ρ − ρ )g
Condições de contorno
u(y)
=
g
y=δ,
( ρ − ρ )
l
µ
l
v
δ
2
l
du
=
dy
⎡y
⎢
⎢⎣
δ
0
−
1
2
v
⎛
⎜
⎝
y
δ
⎞
⎟
⎠
2
(1)
13º ENCIT Congresso Brasileiro de Engenharia e Ciências Térmicas
UFU-Uberlândia– 05-09/12/2010
⎤
⎥
⎥⎦
(2)

O Modelo de Nusselt (2)
Re
Simplificações na equação de Navier Navier-Stokes Stokes
x
δ
=
CL
ρlum
D
μ
D
hidr
.
l
y
hidr.
dp
dx
=
4.
A 4
= =
P
m
∂u
∂y
.
= ρ g
l
y=δ
v
m D
μ bδ
hidr.
( b.
δ )
= 4δ
b
= 0
l
.
(Dedução da Eq. 1)
⎛
ρl
⎜
⎝
4m
4Γ
= =
μ b μ
∂u
∂x
∂u
⎞
+ v ⎟
∂y
⎠
u = − l
.
l
(3)
b
≈ 0
.
m
b
= Γ
2
∂p
∂ u
+ ρ g + μ 2
∂x
∂y
δ ( x)
( x)
=
∫
δ
0
( x)
ρ u
l
( y)dy
.
3
m gρl
( )
( ρl
− ρv
) δ
= Γ x =
(4)
b
3μl
13º ENCIT Congresso Brasileiro de Engenharia e Ciências Térmicas
UFU-Uberlândia– 05-09/12/2010

x
x + dx
O Modelo de Nusselt (3)
"
q .
sup
( b dx)
dq = hlvd
"
dq = q .
sup
.
m
( b dx)
.
m
.
m+ d m
dΓ
dx
.
=
Após derivar Γ, Eq. (4):
dq = hlvd
q
"
sup
h
lv
dΓ
dx
d
=
=
.
m
k
l
.
m
dq = hlvd
"
dq = q .
( T −T
)
ρ
l
sat
δ.
h
lv
sup
( ρ − ρ )
l
μ
l
v
sup
2
gδ
.
m
( b dx)
dδ
dx
δδ
( x )
⎡ 4 k l lμμ
l s at −
= ⎢
⎣⎢
gρl
ρl
− ρv
( x)
1
4
( T − T ) x x⎤
sup
( ) sup
δ ; x ↑
'
hlv = hlv
correção
h
lv
⎥
⎥⎦
( 1+
0,
68Ja)
Ja: Número de Jakob
13º ENCIT Congresso Brasileiro de Engenharia e Ciências Térmicas
UFU-Uberlândia– 05-09/12/2010

h
L
O Modelo de Nusselt (4)
h
h =
x
kl
δ
( T T )
q = h −
x
"
sup x sat ssup
( ρ − ρ )
3 '
⎡
⎤
4 ( ) sup ⎥ ⎥
⎡ g
⎤
l l v k h l lv
= ⎢
⎢⎣
μl
Tsat
−T
x ⎦
⎡ g ρl
= 0,
943⎢
⎢⎣
μl
1
4
ρ O h local, hx, diminui com x.
Após integrar
( ρ − ρ )
l
( T −T
)
sat
v
sup
k
3
l
L
h
'
lv
⎤
⎥
⎥⎦
1
4
Nu
L
.
T
ρ
m =
v
f
q
h
Propriedades:
'
lv
=
e
⎡ g ρl
= 0,
943⎢
⎢⎣
μlk
T
sat
h
lv
+
2
à
( ρ − ρ )
l
h
l
L
Tsup
T
h
'
lv
sat
( T −T
)
sat
A
( T −T
)
v
sat
h
'
lv
sup
L
3
sup
13º ENCIT Congresso Brasileiro de Engenharia e Ciências Térmicas
UFU-Uberlândia– 05-09/12/2010
=
⎤
⎥
⎥⎦
1
4

Limites para cada regime de escoamento
Re ≤
δ
30
30 ≤ Re ≤ 1800
≤ δ
Re ≥
δ
1800
Re
Re
Re
δ
δ
δ
⎡
⎤
⎢
⎥
⎢κ
( ) ⎥
lL
Tsat
−T
p
= 3,
78⎢
1 ⎥
⎢
2 3 ⎥
' ⎛ν
⎞ l ⎢ μlhlv
⎜
⎟ ⎥
⎢⎣
⎝ g ⎠ ⎥⎦
0,
82
3,
70 ( ) ⎥ ⎥⎤
⎡
⎢
⎢3,
70κ
lL(
Tsat
−T
p ) ⎥
= ⎢ + 4,
8
1 ⎥
⎢
2 3 ⎥
'
⎢
⎥
⎢
⎦⎥
⎢ κl
L Tsat
−T
p
=
+
⎛ν
⎞ l μlhlv
⎜
⎟
⎣ ⎝ g ⎠
⎡
⎢
⎢0,
069κ
lL(
Tsat
−T
p ) 0,
5
= ⎢
Pr −151Pr
1
l
l
⎢
2 3
' ⎛ν
⎞ l
⎢ μlhlv
⎜
⎟
⎢⎣
⎝ g ⎠
3
4
0,
5
⎤
⎥
⎥
+ 253⎥
⎥
⎥
⎥⎦
13º ENCIT Congresso Brasileiro de Engenharia e Ciências Térmicas
UFU-Uberlândia– 05-09/12/2010
4
3

Número de Nusselt modificado:
condensação em película sobre placa vertical
In Incropera et al. ( 2008), LTC, p. 409
13º ENCIT Congresso Brasileiro de Engenharia e Ciências Térmicas
UFU-Uberlândia– 05-09/12/2010

Número de Nusselt :
condensação em película
sobre placa vertical
δ ↑ x
h
x
↑
↓ δδ
↑
In Kern
13º ENCIT Congresso Brasileiro de Engenharia e Ciências Térmicas
UFU-Uberlândia– 05-09/12/2010

λ
Taylor
=
⎡
σ
⎤
2π 3⎢
⎥
⎣(
ρl
− ρv
) g ⎦
1/
2
λTaylor
Condensação em película
sobre um só tubo e sobre
uma coluna com dois tubos
horizontais
horizontais.
13º ENCIT Congresso Brasileiro de Engenharia e Ciências Térmicas
UFU-Uberlândia– 05-09/12/2010

Modelo de Nusselt: condensação em película
para tubo horizontal
h
D
Nu
=
D
0,
729
=
0,
729
⎡ g ρl
⎢
⎣ ⎣⎢
μl
⎡ g
⎢
⎢⎣
ρ
( ρ − ρ )
l
( T −T
)
l
μ k
l
sat
v
sup
( ρ − ρ )
l
l
k
h
3
l
h
D
'
lv
( T −T
)
sat
v
sup
13º ENCIT Congresso Brasileiro de Engenharia e Ciências Térmicas
UFU-Uberlândia– 05-09/12/2010
'
lv
D
3
⎤
⎥
⎦ ⎦⎥
⎤
⎥
⎥⎦
1
4
1
4

hD
Correção para o cálculo de para uma
coluna vertical com tubos horizontais
Com drenagem contínua:
Nu
D,
N
=
0,
729
h
h
⎡ g ρl
⎢
⎢⎣
Nμl
D , N
D , N
=
=
( ρ − ρ )
l
k
( T −T
)
h
N
h
N
sat
Nusselt
1/
4
Com drenagem por gotejamento:
Nusselt
1/
6
v
sup
3
l
h
D
'
lv
⎤
⎥
⎥⎦
Correlação de Kern
13º ENCIT Congresso Brasileiro de Engenharia e Ciências Térmicas
UFU-Uberlândia– 05-09/12/2010
1
4

Exemplo C1: condensação em película
sobre uma placa plana vertical
Vapor d’água saturado a 1 atm condensa sobre uma placa vertical com 2,5 m de altura
e 0,5 m de largura, mantida à temperatura de 54ºC. Calcule a taxa de calor transferido à
placa, a vazão de condensado e o limite x para validade do modelo de Nusselt.
Propriedades da água líquida, à
T
f
ρρ
k
l
Pr
Tsat
+ Tp
=
2
= 973 , 7 kg / m
= 0,
668W
/( mK)
μ = 3,
65x10
c
l
l
pl
σ = 0,
063N
/ m
=
2,
29
−4
3
Pa.
s
= 4195J
/( kgK)
Propriedades à T sat
3
ρ = 0,
596kg
/ m
h
v
lv
= 2257kJ
/ kg
b
Correção de h lv
'
h = h lv
lv
( 1 + 0 , 68 Ja )
( T − T )
c pl sat − p
Ja =
hlv
'
Ja = 0,
0855 → hlv
= 2388,
22kJ
/ kg
, Número de Jakob
L 1ª tentativa, hipótese laminar com ondulações
h
L
⎛ ν
⎜
⎝ g
k
l
2
⎞
⎟
⎠
1
3
=
1,
08
Re
Re
δ
δ
30 Re ≤1800
−
≤ δ
5,
2
Correlação
Recomendada
por
Kutateladze

Re
Exemplo C1 (cont.)
30 ≤ Reδ
≤
Re
δ
1800
⎡
⎢
⎢3,
70κ
lL
sat
= ⎢
⎢
2
' ⎛ν
⎞ l ⎢ μlhlv
⎜
⎟
⎢⎣
⎝ g ⎠
( T −T
)
1
3
p
⎤
⎥
⎥
+ 4,
8⎥
⎥
⎥
⎥⎦
0,
82
2ª tentativa, região turbulenta
δ
= 2425 > 1800
Portanto, deve ser o regime turbulento.
⎡
⎢
⎢0,
069κ
lL(
Tsat
−T
p ) 0,
5
= ⎢
Pr −151Pr
1
l
l
⎢
2 3
' ⎛ν
⎞ l ⎢ μlhlv
⎜
⎟
⎢⎣
⎝ g ⎠
0,
5
⎤
⎥
⎥
+ 253⎥
⎥
⎥
⎥⎦
4
3
=
2941,
2
Re
.
m
δ
=
l
.
4 m
=
μ b
μlb
Re
4
.
δ
= 0,
134kg
/ s
'
q = m h hlv
= 320481
h
L
=
hL =
A
p
q
( T −T
)
sat
5573,
6
W
sup
/( m
W
2 K
)

Exemplo C1 (cont.)
Re ≤
δ
30
, Domínio para o modelo de Nusselt
Re
δ
⎡
⎢
⎢κ
l x sat
= 3,
78⎢
⎢
2
' ⎛ν
⎞ l ⎢ μlhlv
⎜
⎟
⎢⎣
⎢⎣ ⎝ g ⎠
( T −T
)
p
1
3
⎤
⎥
⎥
⎥
⎥
⎥
⎥⎦
⎥⎦
3
4
= 30
3/
4
x = 0,
0338;
x = 0,
011m
Valor local de h cond, em x=0,011m
hx =
hL =
9.
730,
3
W
12.
973,
7
/( m
W
2 K
/( m
)
2 K
)

Condensação em convecção forçada
em mini e microcanais
13º ENCIT Congresso Brasileiro de Engenharia e Ciências Térmicas
UFU-Uberlândia– 05-09/12/2010

Aplicações
Condensador de chiller
http://www.carriercca.com
Multiportas de alumínio extrudado
Utilizadas em microcondensadores
www.hydro.com

Condensação em canais convencionais
• Formulações analíticas permitem cálculo do
coeficiente de transferência de calor por
convecção e queda de pressão para os
regimes de escoamento.
• Forças gravitacionais dominantes sobre
forças de atrito e de tensão superficial.
13º ENCIT Congresso Brasileiro de Engenharia e Ciências Térmicas
UFU-Uberlândia– 05-09/12/2010

Condensação em canais convencionais
• Regimes de escoamento:
Regimes bifásicos normalmente encontrados
em condensação (Stephan, 1992)
13º ENCIT Congresso Brasileiro de Engenharia e Ciências Térmicas
UFU-Uberlândia– 05-09/12/2010

Condensação em microcanais
• Aumento da taxa de transferência de calor e
da queda de pressão.
• A força de tensão superficial torna-se
dominante sobre a aceleração da gravidade e o
atrito.
• Modelos propostos para o cálculo do h e de ΔΔΔΔp
em macrocanais não funcionam muito bem em
microcanais.
13º ENCIT Congresso Brasileiro de Engenharia e Ciências Térmicas
UFU-Uberlândia– 05-09/12/2010

Mapa de padrões de escoamento em
microcanais
Padrões de escoamento para a condensação em microcanais
encontrados por Coleman e Garimella (2003)
R-134a em canais retangulares, circulares e quadrados;
1 mm < D h < 4,8 mm
150 kg/m²s < G < 750 kg/m²s

Condensação em convecção forçada
Mapa de padrões desenvolvido por Coleman e Garimella
para condensação em canais circulares com D=1 mm
13º ENCIT Congresso Brasileiro de Engenharia e Ciências Térmicas
UFU-Uberlândia– 05-09/12/2010

Condensação em convecção forçada do
R-134a, 134a, no interior de microcanais
(baseado nos resultados preliminares do
trabalho de mestrado
de Gil Goss Jr., LEPTEN/Boiling - UFSC)
Goss Jr., G., Marini, S.F., Passos, J.C., Heat Transfer and Pressure
Drop Condensation of R-134a inside Parallel Microchannels, in Proceedings of
the ASME/JSME 2011, 8th Thermal Engineering Joint Conference – AJTEC2011, March 13-17,
2011, Honolulu, Hawaii, USA, 9 pages.
13º ENCIT Congresso Brasileiro de Engenharia e Ciências Térmicas
UFU-Uberlândia– 05-09/12/2010

Condições de teste
• Estudo da condensação convectiva em
microcanais:
– Microcanais paralelos, circulares, D = 0,8 mm.
– R-134a
– 50 kg/m²s < G < 420 kg/m²s
– 6 bar < p < 12 bar
• Dados experimentais x correlações e modelos:
– Coeficiente de transferência de calor por
convecção (h)
– Queda de pressão
13º ENCIT Congresso Brasileiro de Engenharia e Ciências Térmicas
UFU-Uberlândia– 05-09/12/2010

Condições de teste
13º ENCIT Congresso Brasileiro de Engenharia e Ciências Térmicas
UFU-Uberlândia– 05-09/12/2010

R134a
Localização termopares
Seção de teste
R134a
Características:
Composta de 8 microcanais
(tubos de cobre, D = 0,8 mm)
paralelos;
Instrumentada com:
16 termopares;
2 transdutores de pressão Warme;
3 coolers Peltier.

Resultados preliminares
Mass velocity (kg.m -2 s -1 Mass Velocidade velocity Mássica (kg.m (kg/m²s) )
800
700
600
500
400
300
200
Dispersed
Bubble
Flow
Plug/Slug
Flow
Plug/Slug and
Annular Film Flow
Annular Film
and Mist Flow
Annular Film Flow
Mist Flow
100
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Título Vapor de quality vapor
Mapa de padrões
desenvolvido por
Coleman e Garimella
in, Goss Jr., G., Macarini, S. e Passos, J.C. (ASME-2011, Hawaii)
13º ENCIT Congresso Brasileiro de Engenharia e Ciências Térmicas
UFU-Uberlândia– 05-09/12/2010

Resultados
Experimental Heat Transfer Coefficient (W.m -2 K -1 Coeficiente de transferência de calor experimental (W/m²K)
Experimental Heat Transfer Coefficient (W.m )
8000
7500
7000
6500
6000
5500
5000
4500
4000
3500
G=65 kgm -2 s -1 p in =7.0 bar q"=24 kW.m -2
G=95 kgm -2 s -1 p in =9.2 bar q"=29 kW.m -2
G=100 kgm -2 s -1 p in =9.3 bar q"=30 kW.m -2
G=125 kgm -2 s -1 p in =11.2 bar q"=34 kW.m -2
3000
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Título de Vapor vapor quality
13º ENCIT Congresso Brasileiro de Engenharia e Ciências Térmicas
UFU-Uberlândia– 05-09/12/2010

Resultados
Heat Transfer Coefficient (W.m -2 Coeficiente de transferência
Heat Transfer
de
Coefficient
calor experimental
(W.m
(W/m²K)
K)
15000
10000
5000
0
Traviss et al. (1973)
Shah (1976)
Moser et al.(1998)
Soliman (1968)
Dobson and Chato (1998)
Yan and Lin (1999)
Koyama (2003)
Experimental
G= 105 kg.m -2 s -1
q"= 30 - 33 kW.m -2
p in = 9.3 bar
T sat =36 o C
0.3 0.4 0.5 0.6 0.7 0.8 0.9
Título Quality de vapor (-)
13º ENCIT Congresso Brasileiro de Engenharia e Ciências Térmicas
UFU-Uberlândia– 05-09/12/2010

Resultados
Correlação de
Yan e Lin (1999)
Calculated Heat Transfer Coefficient (W.m -2 K -1 Calculated Coeficiente Heat de transferência transferência Transfer Coefficient de calor (W.m calculado calculado)
6000
5000
4000
3000
2000
1000
G=60 kgm -2 s -1
G=65 kgm -2 s -1
G=100 kgm -2 s -1
G=105 kgm -2 s -1
G=130 kgm -2 s -1
+30%
-30%
0
0 1000 2000 3000 4000 5000 6000
Experimental Heat Transfer Coefficient (W.m -2 K -1 Coeficiente de transferência de calor experimental (W/m²K)
)
13º ENCIT Congresso Brasileiro de Engenharia e Ciências Térmicas
UFU-Uberlândia– 05-09/12/2010

Agradecimento aos órgãos financiadores
13º ENCIT Congresso Brasileiro de Engenharia e Ciências Térmicas
UFU-Uberlândia– 05-09/12/2010

Muito obrigado pela atenção.
jpassos@emc.ufsc.br
http://energetique-juliocesarpassos.blogspot.com
www.lepten.ufsc.br
13º ENCIT Congresso Brasileiro de Engenharia e Ciências Térmicas
UFU-Uberlândia– 05-09/12/2010