Projeto Térmico de um

Gerador de Vapor
Dados de projeto:
Vazão de vapor = 50 ton/h
Pressão = 2 MPa (manométrica)
Vapor superaquecido, T = 350 °C

Condições de processo:
Temperatura ambiente: T = 25 °C
Gerador de Vapor Aquatubular
Com superaquecedor, economizador e pré-aquecedor
 Das tabelas de vapor
1-Entalpia do vapor superaquecido 3135 kJ/kg

2-Entalpia do vapor saturado 2800 kJ/kg

3-Entalpia do líquido saturado 920,1 kJ/kg

4-Entalpia do líquido 106,8 kJ/kg

Temperatura de saturação 214,9 °C

T 1

3 2
T sat
4
23/8/2013 19:07 2
s
 Combustível
Carvão do RS
Composição % Composição corrigida
Carbono C 46 41,44%
Hidrogênio H2 4,05 3,65%
Nitrogênio N2 0,29 0,26%
Enxofre S 12,99 11,70%
Oxigênio O2 9,27 8,35%
Cinzas CZ 27,4 24,68%
Água 9,91%
Somatório 100 100,00%

Água H2O 11
 Poder Calorífico
 O
Pci = 33900 C + 141800  H −  + 9200 S − 2440 (W + 9 H )
 8

Pcs = 18819 kJ
kg

Pci = 17776 kJ
kg

23/8/2013 19:07 4
 Cálculo da combustão
Volume de ar teórico
Vart = 8,876 C + 28,238 H + 3,324 S − 3,56 O

3
Nm de ar
Vart = 4,75
Kg de combustível

23/8/2013 19:07 5
 Massa de ar teórico
mart = 11,51 C + 34,298 H + 14,36 S − 4,32 O
Kg de ar
mart = 7,34
Kg de combustível

Excesso de ar
r r
V m
e= =
ar
t
= 1,5 ar
t
V arm ar

23/8/2013 19:07 6
 Volume de ar Real
3
Nm de ar
Varr = 7,13
Kg de combustível

Massa de ar Real
Kg de ar
m = 11,01
r
ar
Kg de combustível

23/8/2013 19:07 7
 Volume de gás teórico
Vgt = 1,853 C + 0,683 S + 0,8 N + 0,79 Vart + (9 H + W )1,24
3
Nm de gás
Vgt = 5,14
Kg de combustível

Massa de gás teórico

mgt = 3,664 C + 8,94 H + 1,998 S + N + 0,77mart + W
Kg de gás
m = 7,83
t
g
Kg de combustível
23/8/2013 19:07 8
 Gás real

3
Nm de gás
Vg = Vg + (e − 1) Var = 7,51
r t t

kg de combustível

kg de gás
m = m + (e − 1) m = 11,50
r
g
t
g
t
ar
kg de combustível

23/8/2013 19:07 9
 Composição dos gases da
combustão (% em volume)
Vgr = 1,864 C + 0,683 S + 0,8 N + 0,79 Vart + K
K (9 H + W )1,24 + (e − 1)Vart

1,864 C
Teor de CO2 = r
= 10,28%
Vg

0,683 S
Teor de SO2 = r
= 1,06%
Vg

23/8/2013 19:07 10
 Vgr = 1,864 C + 0,683 S + 0,8 N + 0,79 Vart + K
K (9 H + W )1,24 + (e − 1)Vart

0,8 N + 0,79 Vart

Teor de N 2 = r
= 50,01%
Vg

Teor de H 2 O=
(9 H + W )1,24
= 7,06%
Vgr

Teor de Ar =
(e − 1) Vart
= 31,63%
r
V g

23/8/2013 19:07 11
 Composição dos gases da
combustão (% em massa)
mgr = 3,664 C + 8,94 H + W + 1,998 S + K
( )
K 0,77 mar + N + e − 1 mar
t t

3,664 C
Teor de CO2 = r
= 13,20%
mg

8,94 H + W
Teor de H 2O = r
= 3,70%
mg

23/8/2013 19:07 12
 mgr = 3,664 C + 8,94 H + W + 1,998 S + K
K 0,77 mart + N + (e − 1)mart

1,998 S
Teor de SO2 = r
= 2,03%
mg

0,77 mart + N
Teor de N 2 = r
= 49,16%
mg

Teor de Ar =
(e − 1)mart
= 31,91%
r
m g

23/8/2013 19:07 13
 Peso específico dos gases da
combustão
mgr 11,50
ρg = r = = 1,53 kg 3
0 °C e 1 atm
Vg 7,51 Nm
 273  kg
a 500 °C ρ g = 1,53.  = 0,54 3
 273 + 500  m

ou
ρ g = V (CO2 ).ρ CO + V (H 2O ).ρ H O + K
2 2

KV ( N 2 ).ρ N 2 + V (O2 ).ρ O2 + V (SO2 ).ρ SO2

14
 Peso específico dos gases da
combustão
1
ρg =
mi (% )
∑
i ρi

1
ρg =
m(CO2 ) m(H 2O ) m( N 2 ) m(O2 ) m(SO2 )
+ + + +
ρ CO 2
ρH O2
ρN 2
ρO
2
ρ SO
2

23/8/2013 19:07 15
 Calor específico dos gases da
combustão

Cp g = ∑ mi (% ).Cpi Cpi em kJ
i
kg.o C

Cp g = ∑ Vi (% ).Cpi Cpi em kJ
i
Nm 3 .o C

23/8/2013 19:07 16
Equações para determinação
das propriedades em
“Geradores de Vapor, Edson
Bazzo”

Viscosidade do SO2 em “Perry

& Chilton, Manual de
Engenharia Química”

Condutividade térmica do SO2:

 10,39  k em W o
k g = µ g  Cp g +  × 1000 m. C
µ em Pa.s
 M g 
 Cp em kJ o
kg. C
23/8/2013 19:07 17
 Condutividade térmica dos gases
da combustão
∑ y i k i (M i )
1
3

kg = i

∑ yi (M i )
1
3

yi = fração molar do componente i

ki = condutividade térmica do componente i
Mi = massa molecular

Equação de Eucken (alternativa)

k em BTU
 2,48  h. ft.o F

k g = µ g Cp g + µ em lb h. ft
 M 
 g  Cp em BTU
lb.o F
23/8/2013 19:07 18
 Viscosidade dos gases da
combustão
Viscosidade dinâmica para uma mistura de gases:

∑ yi µi (M i )
1
2

µg = i

∑ yi (M i )
1
2

yi = fração molar do componente i

µi = viscosidade dinâmica do componente i
Mi = massa molecular

23/8/2013 19:07 19
 Determinação das perdas
P1 – Calor do combustível caído no cinzeiro
Estimada 1,0 %
P2 – Calor sensível nas cinzas
Z .Cp.∆T
P2 = × 100 = 1,96%
Pci
Z = porcentagem de cinzas no combustível = 24,68 %
Cp = calor específico das cinzas ⇒ Cp = 0,3 kcal/kg.°C
∆T = diferença de temperatura entre o combustível na
entrada da fornalha e as cinzas no cinzeiro
Tf - Ta ⇒1050-25= 1025 °C
23/8/2013 19:07 20
 P3 – Perdas de fuligem
Estimada 1,0 %
P4 – Perdas de combustão incompleta
Estimada 0,0 %
P5 – Perdas por irradiação
Estimada 1,0 %
P6 – Perdas de calor sensível
mgr .Cpm .∆T
P6 = ×100 = 16,97%
Pci
mgr = massa real dos gases por unidade de massa do combustível = 11,5 kg gás/kg comb
Cpm = calor específico médio dos gases = 1,028 kJ/kg.°C
∆T = diferença de temperatura entre os gases na saída da chaminé e o ar (T ambiente)
Temperatura na chaminé = 280 °C (arbitrada)

23/8/2013 19:07 21
 Perda total
pt = p1 + p2 + p3 + p4 + p5 + p6 = 21,93%

Rendimento
η = 1 − pt = 78,07%

Calor útil
Q& u = D (hv − ha ) = 42058 kW
hv = entalpia do vapor = 3135 kJ/kg
ha = entalpia da água = 106,8 kJ/kg
23/8/2013 19:07 D = produção de vapor = 13,9 kg/s 22
Pré-aquecimento de ar
Temperatura do ar pré- aquecido = 140 °C (valor
recomendado Hildo Pera)
Tm = 82,5 °C ⇒ Cpar = 0,996 kJ/kg.°C

Consumo de combustível
Q& u
B=
η .[Pci + marr .Cpar .(Tar − Tamb )]

B = 2,83 kg = 10187,57 kg
s h
23/8/2013 19:07 23
 Calor total fornecido
Q& t = B.Pci + B.marr .Cpar .(Tar − Tamb )

Q& t = 53872,80 kW

23/8/2013 19:07 24
 Dimensionamento da fornalha
Área do suporte
B B
Kg = Sg = = 25,47 m 2

Sg Kg

Kg = carga da grelha (kg/m2.h) = 400

B = consumo de combustível (kg/h) = 10187,57
Sg = área do suporte (m2)

23/8/2013 19:07 25
 Volume da fornalha

B ⋅ Pci
Kf = V = 226,36 m 3
V

Kf = carga da fornalha (kJ/m3.h) = 800 000

B = consumo de combustível (kg/h) = 10187,57
Pci = poder calorífico inferior (KJ/kg) =17 776
V = volume da câmara (m3)
23/8/2013 19:07 26
 Dimensões da câmara de combustão

S g = 25,47 m < L × C = 32,2 m

2 2

V = 226,36 m < A × L × C = 241,5 m

3 3

A
Altura - A = 7,5 m
Largura - L = 4,6 m
Comprimento - C = 7,0 m C
L

23/8/2013 19:07 27
 Temperatura dos gases na
câmara de combustão
• Temperatura adiabática (teórica)

Pci + Varr .Cpar (Tar − Tamb )

T f = Tamb +
Vg .Cp g + [Cz ].Cpcz
r

T f = 1523 oC

23/8/2013 19:07 28
• A temperatura dos gases na câmara de combustão
é determinada por processo iterativo
B.Pci + B.Varr .Cpar (Tar − Tamb ) − Q& r − Q& p
T f = Tamb +
B.Vgr .Cp g + B.[Cz ].Cpcz

Q& r  1  de  1 
T p = Tv +  ln  + 
2.π .L.N t  kt  d i  hi .d i 

& (
Qr ≅ σ .ε .S i . T f − T p
4 4
)
23/8/2013 19:07 29
 Passo dos tubos na parede d’água, s=2xde
Número de tubos  2 L + 2C 
Tubos:
N t = INT   = 158 DN = 2½ “
 s 
de = 73,2 mm
hi estimado 5000 W/m .°C 2
di = 62,7 mm
Calor perdido na fornalha
Q p = ( p1 + p2 + p3 + p4 + p5 )× Qt = 2673,27 kW
Superfície irradiada
Comprimento dos tubos da parede d’água = 6 m

Si = f .Lt .de.N t = 62,41 m 2

23/8/2013 19:07 30
Resolvendo o sistema de equações T f = 1161,3 oC

T p = 226,4 o C

Qr = 10328,93 kW

23/8/2013 19:07 31
 Aproveitamento da energia do
T1 T2
T3 combustível
Superaquecedor

Chaminé
Paredes
d’água
Tch
Tf

Fornalha

Pré-aquecedor de ar
Caldeira
23/8/2013 19:07 32
Economizador
Superaquecedor D.(hvsa − hv ) = B.Vgr .Cp g (T f − T1 )
Tm = 1097,9 oC ⇒ Cp g = 1,128 kJ
kg.o C
Qs = 4652,8 kW T1 = 1034,6 oC

Caldeira D.(hv − hla ) − Q& r = B.Vgr .Cp g (T1 − T2 )

A água sai do economizador e entra na caldeira
a 179,9°C, recomendação Hildo Pera

hla = 763,8 kJ Tm = 783,3 oC ⇒ Cp g = 1,097 kJ

kg kg.o C

Qc = 17951,6 kW T2 = 532,0 oC
23/8/2013 19:07 33
Economizador D.(hla − hl ) = B.Vgr .Cp g (T2 − T3 )
hla = 763,8 kJ o ⇒ Tla = 179,9 oC
kg. C

hl = 334,9 kJ ⇒ Tl = 80 oC
kg.o C

Tm = 445,9 oC ⇒ Cp g = 1,063 kJ
kg.o C

Qe = 5956,9 kW T3 = 359,8 oC

23/8/2013 19:07 34
Pré-aquecedor de ar
B.Varr .Cpat .(Tar − Tamb ) = B.Vgr .Cp g (T3 − Tch )

Temperatura de pré-aquecimento do ar ⇒ Tar = 140 °C

Tm = 306,5 oC ⇒ Cp g = 1,047 kJ o
kg. C

Qe = 3630,9 kW Tch = 253,3 oC

Temperatura estimada no rendimento Tch = 280 oC

23/8/2013 19:07 35
 Área de transferência de calor
Superaquecedor

Lado do vapor Lado do gás

Calor trocado = 4652,8 kW Temp de entrada do gás = 1161,3 C
Temp de entrada do vapor = 214,9 C Temp de saída do gás = 1034,6 C
Temp de saída do vapor = 350 C Vazão de gás = 32,6 kg/s
Vazão de vapor = 13,9 kg/s
Temp média do gás = 1097,9 C
Temp média do vapor = 282,45 C
Massa específica gás = 0,305 kg/m3
Massa específica vapor = 8,74 kg/m3
Calor específico vapor = 1,997 kJ/kg.K Calor específico gás = 1,128 kJ/kg.K
Viscosidade dinâmica = 1,88E-05 Pa.s Viscosidade dinâmica gás = 5,082E-05 Pa.s
Condutividade térmica = 3,79E-02 W/m.C Condutividade térmica gás = 0,087 W/m.C

23/8/2013 19:07 36
Largura máxima do superaquecedor é a largura da fornalha
L=4,6 m

gás

Fluxo cruzado
comprimento do superaquecedor 2m
tubos DN= 1,5 pol =
de = 1,9 pol = 48,26 mm
di = 1,61 pol = 40,89 mm
e= 0,145 pol = 3,68 mm
passo paralelo sp= 2 x de = 96,52 mm
passo transvesal st= 2 x de = 96,52 mm
disposição em linha
número de fileiras paralelas ao escoamento = 46 Máx=47
número de fileiras transversais ao escoamento = 4
 Para fluxos cruzados
∆T1 = Tqe − Tf s = 1161,3 − 350 = 811,3
∆T2 = Tqs − Tf e = 1034,6 − 214,9 = 819,7
∆T1 − ∆T2
∆Tml = = 815,5o C
∆T1
ln
∆T2
Q& = U . A.∆Tml 1
U .A =
R1 + R2 + R3

1 de 1
R1 = ln R3 =
π .d i .L.hi di π .d e .L.he
R2 =
23/8/2013 19:07
2.π .L.kt 38
 0 ,14
 µ 
1
Nu = 0,027.R e .P r . 
0 ,8 3
µ 
 p

Coeficiente de película do vapor - hi

Área de escoamento= 0,06 m2
Velocidade do vapor = 26,30 m/s
Reynolds = Re = 4,99E+05
Prandtl = Pr = 9,93E-01
Nusselt = Nu = 9,74E+02 Sieder-Tate
Coeficiente de película - hi = 903,1 W/m2.C

23/8/2013 19:07 39
 Nu = 0,27. f1.R e 0 , 63
.P r 0 , 36

Coeficiente de película do gás - he

Área de escoamento= 4,44 m2
Velocidade do gás = 24,05 m/s
Reynolds = Re = 6,96E+03
Prandtl = Pr = 6,56E-01
Fator de fileiras f1 = 0,96
Nusselt = Nu = 5,87E+01 Zukauskas
Coeficiente de película - hec = 106,3 W/m2.C

23/8/2013 19:07 40
 % H 2O  T 3, 2  T  3, 2  T  0,65 
pw = × pt 0, 4   
Q& c = 9,3. A.( pc .Le ) .
p
100
 −   .  
100  100 T
  p  
%CO2
pc = × pt
100

0, 6  T 
m
 T 
m

Qw = (42 − 76. pw .Le ). A.( pw .Le ) .
& p
 4  s p .st    −  
Le ≅ 0,85.  2  − 1.d e  100  100  
π  d e  

RADIAÇÃO GASOSA
pressão total = 1 atm
pressão parcial CO2 = 0,10 atm
pressão parcial H2O = 0,07 atm
espessura efetiva do gás = 0,17 m
Temperatura da parede = 374,79 C =
Temperatura do gás = 1370,9 K
m= 2,63
Qc r = 6796,71 W/m2
Qw r = 2431,90 W/m2
Qr= 9228,61 W/m2
Coeficiente de película - her = 11,32 W/m2.C
Coeficiente de película - hec = 106,3 W/m2.C 41
Coeficiente de película - he = 117,6 W/m2.C
 &
Q = hi . Ai (Tp − Tm )

23/8/2013 19:07 42
 U = 101,05W
m 2 .o C
Q& s 4652,8
Anec = = = 56,5 m 2
U .∆Tml 101,5 × 815,5

Adisp = π .de.L.N t = π .0,04826.2.(46 × 4 ) = 55,8 m 2

Erro = −1,2%

23/8/2013 19:07 43
7,5 m
Caldeira

4,6

7
ÁREA DE TROCA TÉRMICA - CALDEIRA
Calor trocado = 17951,6 kW
Temp de saturação do vapor = 214,9 C

Lado do gás
Temp de entrada do gás = 1034,6 C
Temp de saída do gás = 532,0 C
Vazão de gás = 32,6 kg/s
Temp média do gás = 783,3 C
Massa específica gás = 0,396 kg/m3
Calor específico gás = 1,097 kJ/kg.K
Viscosidade dinâmica gás = 4,333E-05 Pa.s
Condutividade térmica gás = 0,070 W/m.C

23/8/2013 19:07 45
Fluxo cruzado
comprimento da caldeira = 5m
número passes dos gases = 3
comprimento por passe = 1,67 m
tubos DN= 1,5 pol =
de = 1,9 pol = 48,26 mm
di = 1,61 pol = 40,89 mm
e= 0,145 pol = 3,68 mm
passo paralelo sp= 2 x de = 96,52 mm
passo transvesal st= 2 x de = 96,52 mm
disposição em linha
número de fileiras paralelas ao escoamento = 47 Máx = 47
número de fileiras transversais ao escoamento = 8

Diferença média logarítmica de temperatura

DT1= 819,7 C MLDT = 529,2 C
DT2= 317,1 C

Coeficiente de película - hi = 5000,0 W/m2.C

23/8/2013 19:07 46
 RADIAÇÃO GASOSA
Coeficiente de película do gás - he pressão total = 1 atm
Área de escoamento= 3,78 m2 pressão parcial CO2 = 0,10 atm
Velocidade do gás = 21,76 m/s pressão parcial H2O = 0,07 atm
Reynolds = Re = 9,59E+03 espessura efetiva do gás = 0,17 m
Prandtl = Pr = 6,78E-01 Temperatura da parede = 227,50 C =
Fator de fileiras f1 = 1 Temperatura do gás = 1056,3 K
Nusselt = Nu = 7,57E+01 Zukauskas m= 2,63
Coeficiente de película - hec = 110,0 W/m2.C Qc r = 2947,12 W/m2
Coeficiente de película - he = 117,85 Qw r = 1222,79 W/m2
Qr= 4169,91 W/m2
Coeficiente de película - her = 7,88 W/m2.C

Coefiente global de troca térmica U= 113,51 W/m2.C

Área necessária = 298,87 m2

Área disponível = 285,0 m2

Área necessária = 298,9 m2

Erro = -4,6%
23/8/2013 19:07 47
 Economizador
Calor trocado = 5956,9 kW
Temp de entrada da água = 80,0 C Temp de entrada do gás = 532,0 C
Temp de saída da água = 179,9 C Temp de saída do gás = 359,8 C
Vazão de água = 13,9 kg/s Vazão de gás = 32,6 kg/s
Temp média da água = 129,95 C Temp média do gás = 445,9 C
Massa específica da água = 937,2 kg/m3 Massa específica gás = 0,581 kg/m3
Calor específico da água = 4,66 kJ/kg.K Calor específico gás = 1,063 kJ/kg.K
Viscosidade dinâmica = 1,18E-04 Pa.s Viscosidade dinâmica gás = 3,291E-05 Pa.s
Condutividade térmica = 0,64 W/m.C Condutividade térmica gás = 0,051 W/m.C

Fluxo cruzado
comprimento do economizador 3m largura = 2 m
número de passes nos gases = 1
tubos DN= 1 pol =
de = 1,32 pol = 33,53 mm
di = 1,049 pol = 26,64 mm
e= 0,1355 pol = 3,44 mm
passo paralelo sp= 2 x de = 67,056 mm
passo transvesal st= 2 x de = 67,056 mm
disposição em linha
número de fileiras paralelas ao escoamento = 24
número de fileiras transversais ao escoamento = 20
número de passes nos tubos = 1
Diferença média logarítmica de temperatura
DT1= 352,1 C MLDT = 314,5
DT2= 279,8 C

0 ,14
 µ1
Nu = 0,027.R e .P r 3 . 
0 ,8
µ 
 p

Coeficiente de película do vapor - hi

Área de escoamento= 0,27 m2
Velocidade da água = 0,06 m/s
Reynolds = Re = 1,17E+04
Prandtl = Pr = 8,57E-01
Nusselt = Nu = 4,61E+01 Sieder-Tate
Coeficiente de película - hi = 1111,6 W/m2.C
 Nu = 0,27. f1.R e 0 , 63
.P r 0 , 36

RADIAÇÃO GASOSA
Coeficiente de película do gás - he pressão total = 1 atm
Área de escoamento= 2,41 m2 pressão parcial CO2 = 0,10 atm
Velocidade do gás = 23,19 m/s pressão parcial H2O = 0,07 atm
Reynolds = Re = 1,37E+04 espessura efetiva do gás = 0,12 m
Prandtl = Pr = 6,84E-01 Temperatura da parede = 165,28 C =
Fator de fileiras f1 = 1 Temperatura do gás = 718,9 K
Nusselt = Nu = 9,52E+01 Zukauskas m= 2,60
Coeficiente de película - hec = 145,3 W/m2.C Qc r = 626,22 W/m2
Coeficiente de película - he = 148,22 W/m2.C Qw r = 281,59 W/m2
Qr= 907,81 W/m2
Coeficiente de película - her = 2,89 W/m2.C

U= 125,56 W/m2.C

Área necessária = 150,85 m2 Área disponível = 151,7 m2

Erro = 0,6%
23/8/2013 19:07 50
 Pré-aquecedor de ar
ÁREA DE TROCA TÉRMICA - PRÉ-AQUECEDOR DE AR
Calor trocado = 3630,9 kW
Temp de entrada do ar = 25,0 C Temp de entrada do gás = 359,8 C
Temp de saída do ar = 140 C Temp de saída do gás = 253,3 C
Vazão de ar = 31,2 kg/s Vazão de gás = 32,6 kg/s
Temp média do ar = 82,5 C Temp média do gás = 306,5 C
Massa específica do ar = 1,186 kg/m3 Massa específica gás = 0,721 kg/m3
Calor específico do ar = 0,996 kJ/kg.K Calor específico gás = 1,047 kJ/kg.K
Viscosidade dinâmica = 2,01E-05 Pa.s Viscosidade dinâmica gás = 2,790E-05 Pa.s
Condutividade térmica = 0,03 W/m.C Condutividade térmica gás = 0,043 W/m.C

Fluxo cruzado
comprimento do pré-aquecedor 4,5 m largura = 3m
número de passes nos gases = 3
tubos DN= 2 pol =
de = 2,38 pol = 60,45 mm
di = 2,067 pol = 52,50 mm
e= 0,1565 pol = 3,98 mm
passo paralelo sp= 2 x de = 120,904 mm
passo transvesal st= 2 x de = 120,904 mm
disposição em linha
número de fileiras paralelas ao escoamento = 20 Máx= 24
número de fileiras transversais ao escoamento = 14
número de passes nos tubos = 1
 0 ,14
 µ 1
Nu = 0,027.R e .P r . 
0 ,8 3
µ 
 p
ar - hi
Coeficiente de película do vapor
Área de escoamento= 0,61 m2
Velocidade do ar = 43,35 m/s
Reynolds = Re = 1,34E+05
Prandtl = Pr = 6,72E-01
Nusselt = Nu = 2,99E+02 Sieder-Tate
Coeficiente de película - hi = 170,0 W/m2.C

23/8/2013 19:07 52
 Nu = 0,27. f1.R e 0 , 63
.P r 0 , 36

RADIAÇÃO GASOSA
Coeficiente de película do gás - he pressão total = 1 atm
Área de escoamento= 1,81 m2 pressão parcial CO2 = 0,10 atm
Velocidade do gás = 24,89 m/s pressão parcial H2O = 0,07 atm
Reynolds = Re = 3,89E+04 espessura efetiva do gás = 0,21 m
Prandtl = Pr = 6,84E-01 Temperatura da parede = 171,76 C =
Fator de fileiras f1 = 1 Temperatura do gás = 579,5 K
Nusselt = Nu = 1,83E+02 Zukauskas m= 2,66
Coeficiente de película - hec = 129,7 W/m2.C Qc r = 272,34 W/m2
Coeficiente de película - hec = 131,70 W/m2.C Qw r = 175,66 W/m2
Qr= 448,00 W/m2
Coeficiente de película - her = 2,00 W/m2.C

U= 69,15 W/m2.C

Área necessária = 234,43 m2 Área disponível = 239,3 m2

Erro = 2,1%
23/8/2013 19:07 53
 Tiragem
 V  2
∆p = K  ρ .  0 ,14
 2   µp 
K = 4. f .N f  
 µm 

 
 sp 
 0 , 08  −0,15
de
f = 0,044 + .R emáx
 [0, 43+1,13.(d e s p )] 
 st − d e 
   
  de  
Perda de carga no superaquecedor
rho = 0,305 kg/m3
V= 24,05 m/s f= 0,0541
de = 48,26 mm K= 8,056E-01
sp = 96,52 mm
st = 96,52 mm Dp = 71,02 Pa =
Re máx = 6,96E+03
T= 1097,9 C 7,1 mm C. A
Tp = 374,79 C
Visc m = 5,082E-05 Pa.s
Visc p = 3,041E-05 Pa.s
Nf = 4 número de fileiras

23/8/2013 19:07 55
Perda de carga na caldeira
rho = 0,396 kg/m3
V= 21,76 m/s f= 0,0516
de = 48,26 mm K= 4,580E+00
sp = 96,52 mm
st = 96,52 mm Dp = 429,11 Pa =
Re máx = 9,59E+03
T= 783,3 C 42,9 mm C. A
Tp = 227,50 C
Visc m = 4,333E-05 Pa.s
Visc p = 2,487E-05 Pa.s
Nf = 8 número de fileiras
Np = 3 Número de passagens

23/8/2013 19:07 56
Perda de carga no economizador
rho = 0,581 kg/m3
V= 23,19 m/s f= 0,0489
de = 33,53 mm K= 3,704E+00
sp = 67,056 mm
st = 67,056 mm Dp = 579,12 Pa =
Re máx = 1,37E+04
T= 445,9 C 57,9 mm C. A
Tp = 165,28 C
Visc m = 3,291E-05 Pa.s
Visc p = 2,239E-05 Pa.s
Nf = 20 número de fileiras
Np = 1 Número de passagens

23/8/2013 19:07 57
Perda de carga no pré-aquecedor de ar / lado do gás
rho = 0,721 kg/m3
V= 24,89 m/s f= 0,0418
de = 60,45 mm K= 6,820E+00
sp = 120,904 mm
st = 120,904 mm Dp = 1523,27 Pa =
Re máx = 3,89E+04
T= 306,5 C 152,3 mm C. A
Tp = 171,76 C
Visc m = 2,790E-05 Pa.s
Visc p = 2,265E-05 Pa.s
Nf = 14 número de fileiras
Np = 3 Número de passagens

23/8/2013 19:07 58
 Perda de carga no interior dos tubos
L
 V  2 K = f.
∆p = K  ρ .  dh
 2  f = (1,82 log10 R e − 1,63)
−2

Perda de carga no pré-aquecedor de ar / lado do ar

rho ar = 1,186 kg/m3
V= 43,35 m/s f= 0,0169
di = 52,50 mm K= 1,445E+00
L= 4,5 m
Np = 1 Dp = 1609,79 Pa =
Re máx = 1,34E+05
T= 82,5 C 161,0 mm C. A

23/8/2013 19:07 59
RESUMO PERDA DE CARGA
Superaquecedor = 7,1 mm C. A
Caldeira = 42,9 mm C. A
Economizador = 57,9 mm C. A
Pré-aquecedor / lado gás = 152,3 mm C. A
Pré-aquecedor / lado ar = 161,0 mm C. A

Somatório perda de carga = 421,2 mm C. A =

23/8/2013 19:07 60
 Chaminé
• Tiragem natural
∆p
Hu =
(ρ a − ρ g ) .g
Somatório perda de carga = 421,2 mm C. A = 4212,31 Pa

T ar = 25 C
T gás ch = 253,3 C
∆p
rho ar = 1,415 kg/m3
Hu = = 692 m
rho gás= 0,794 kg/m3
(ρ a − ρ g ) .g

Utilizar tiragem forçada – tipo balanceada

23/8/2013 19:07 61
 Dimensionamento da chaminé
• Diâmetro da chaminé
4.m& g
d= r
π .ρ g .V
Dimensionamento da chaminé
Altura = 15 m
Velocidade do gás na chaminé = 10 m/s

Diâmetro calculado = 2,28 m

23/8/2013 19:07 62
 Perda de carga na chaminé
L
 V  2 K = f.
∆p = K  ρ .  dh
 2  f = (1,82 log10 R e − 1,63)
−2

Perda de carga na chaminé

rho gás = 0,794 kg/m3
V= 10,00 m/s f= 0,0123
di = 2,28 m K= 8,089E-02
L= 15 m
Dp = 3,21 Pa =
Re máx = 7,01E+05
T= 253,3 C 0,3 mm C. A

23/8/2013 19:07 63
 Tiragem da chaminé
∆pc = H u .(ρ a − ρ g ).g

Ft chaminé= 91,32 Pa = 9,1 mm C. A

23/8/2013 19:07 64
 Potência dos ventiladores
m& a .∆pv m& g .∆pv
N= N=
ρ a .η ρ g .η
Acrescentar 20,0% às perdas, para compensar depósitos e perdas localizadas

INSUFLAMENTO
Vence as perdas no circuito de ar
Pré-aquecedor / lado ar Dp = 161,0 mm C. A + Acréscimo= 193,18 mm C. A = 1931,75 Pa
Vazão de ar = 31,2 kg/s
T ar = 82,5 C
rho ar = 1,186 kg/m3
Rendimento do ventilador = 60,0%

Potência= 84594,53 W = 113 HP = 85 kW

23/8/2013 19:07 65
EXAUSTÃO
Vence as perdas no circuito de gás
Superaquecedor = 7,1 mm C. A
Caldeira = 42,9 mm C. A
Economizador = 57,9 mm C. A
Pré-aquecedor / lado gás = 152,3 mm C. A
Chaminé 0,3 mm C. A

Somatório das perdas de carga Dp 260,6 mm C. A + Acréscimo= 260,57 mm C. A = 2605,73 Pa

Vazão de gás = 32,6 kg/s
T gás = 707,3 C
rho gás = 0,426 kg/m3
Rendimento do ventilador = 60,0%

Potência= 331556,4 W = 444 HP = 332 kW

POTÊNCIA DA BOMBA D'ÁGUA

Vazão = 13,9 Kg/s
T= 129,95 C Pot Bomba = 51868,4498 W = 70 HP = 52 kW
rho água= 937,2 kg/m3
Dp = 2,1 Mpa
rendimento 60,0%

23/8/2013 19:07 66
POTÊNCIA DO GERADOR DE VAPOR

P= 42058333 W 56378 HP = 42058 kW

Os ventiladores representam 0,99% da potência gerada

A bomba representa 0,12% da potência gerada

23/8/2013 19:07 67