CHARAKTERYSTYKA I OPIS TECHNICZNY KOTŁA O WYDAJNOŚCI 130 t/h.
Konstrukcja kotła
Kocioł OP-130 jest kotłem parowym wodnorurowym, opalanym pyłem węgla kamiennengo, z naturalną cyrkulacją, wiszącym na konstrukcji stalowej z obmurzem ciężkim opancerzonym. Komora paleniskowa o wymiarach 5905mm x 6015mm wykonana jest ze ścian membranowych w których średnica rury wynosi 57 x 5 mm, rozstaw ich wynosi 75 mm. Ekrany paleniska zawieszone są na belkach górnego rusztu konstrukcji nośnej wspartych na słupach nośnych. W komorze paleniskowej zabudowany jest przegrzewacz naścienny z rur φ 32 pokrywający strop oraz ścianę przednią od poziomu 16.7 m będący pierwszym stopniem przegrzewacza. Drugi stopień przegrzewacza z rur φ 38 zabudowany jest w międzyciągu, a trzeci, grodziowy z rur φ 32 na wylocie z komory paleniskowej. W drugim ciągu kotła zabudowany jest pięciopęczkowy podgrzewacz wody z rur φ 32 . Ostatnią powierzchnią ogrzewalną w kotle są dwa obrotowe podgrzewacze powietrza. Przygotowanie paliwa odbywa się w trzech młynach średniobieżnych.
Przepływ spalin w komorze paleniskowej jest ku górze, skąd przez kanał nawrotny spaliny kierowane są drugim ciągiem ku dołowi. Ciepło wydzielane przy spalaniu pyłu przekazywane jest na drodze promieniowania i częściowo konwekcji powierzchni ogrzewalnej komory paleniskowej. Niewykorzystane w komorze ciepło oddawane jest powierzchniom ogrzewalnym umieszczonym w drugim ciągu kotła.
Pozostała część ciepła zawarta w spalinach oddawana jest do atmosfery powodując tzw. stratę kominową kotła. Spaliny po wyjściu z kotła poprzez kanały spalin przechodzą do filtra workowego, w którym zatrzymywany jest lotny popiół i dalej poprzez wentylatory ciągu do komina.
Przepływ czynnika roboczego
Woda z podgrzewacza wody doprowadzana jest do walczaka skąd trzema
centralnymi rurami opadowymi φ = 355x28 do tzw. pierścienia dolnego i dalej do dolnych komór ekranowych. Mieszanina parowo-wodna z rur ekranowych wchodzi do górnych komór ekranowych i następnie do walczaka gdzie następuje oddzielenie pary od wody. Para nasycona z walczaka dwunastoma rurami φ = 76,1x7,1 doprowadzana jest do górnej komory zabudowanej na poziomie 28,250m
Para z górnej komory ekranowej przepływa rurami stropowymi paleniska 32x4 i dalej jako tzw. przegrzewacz naścienny pionowy do komory dolnej na poziomie
16,380 m zlokalizowanej na przedniej ścianie kotła. Z komory dolnej para jest odprowadzana dwoma rurami φ = 218x18 na bocznych ścianach kotła poprzez schładzacz 1,5 kg/s do komory górnej, a następnie do przegrzewacza konwekcyjnego.
Przegrzewacz konwekcyjny zbudowany jest w układzie szchownicowym z rur φ = 38. Jego podziałka poprzeczna poprzeczna wynosi 150 mm, a wzdłużna 75 mm. Przepływ odbywa pary odbywa się przeciwprądowo.
Między przegrzewaczem konwekcyjnym, a przegrzewaczem grodziowym zastosowano wtrysk wody chłodzącej 0,8 kg/s.
Podgrzewacz grodziowy umieszczony jest na wylocie z komory paleniskowej i składa się z 18 grodzi po 40 rur φ = 31.8 w każdej z czego 10 grodzi pracuje w przeciwprądziu, natomiast 2x4 grodzie we współprądzie. Po przepływie przez przegrzewacz grodziowy para odprowadzana jest do części wysokoprężnej turbiny.
II. OBLICZENIA CIEPLNE
1. Dane obliczeniowe i znamionowe kotła OP-130
Dane obliczeniowe
Skład paliwa
zawartość węgla Cr = 48.8 %
zawartość wodoru Hr = 4.5 %
zawartość tlenu Or = 7.5 %
zawartość siarki Sr = 1 %
zawartość azotu Nr = 0.8 %
zawartość popiołu Ar = 19 5 %
zawartość wilgoci Wr = 18.4 %
1.1.2. Strumień pary pierwotnej D = 36.111 kg/s
1.1.3. Temperatura powietrza za podgrzewaczem tpg = 230 °C
1.1.4. Temperatura powietrza zimnego tpz = 30 °C
1.1.5. Temperatura wody zasilającej t1=200°C
1.1.6. Nadmiaru powietrza w palenisku λ1 = 1.2
1.1.7. Nadmiar powietrza w spalinach wylotowych λ2 = 1.2
Dane znamionowe
1.2.1. Strumień pary pierwotnej D = 36.111 kg/s
1.2.2. Ciśnienie pary przegrzanej p7 = 12.9 MPa
1.2.3. Ciśnienie wody zasilającej p1 = 14.32 MPa
1.2.4. Ciśnienie w walczaku p2 = 13.24 MPa
1.2.5. Temperatura pary przegrzanej t7 = 540°C
1.2.6. Temperatura spalin wylotowych tw =125 °C
1.2.7. Unos popiołu aun = 0.9
2. OBLICZENIA FIZYKOCHEMICZNE PALIWA I SPALIN
2.1.Skład paliwa
rodzaj paliwa - węgiel kamienny
skład masy roboczej paliwa
Cr = 48.8 Hr = 4.5 Or =7.5 Nr = 0.8
Sr = 1 Wr = 18.4 Ar = 18.4
sprawdzenie
Cr+Hr+Or+Sr+Nr+Ar+Wr=100 %
wartość opałowa paliwa - Formuła Mendelejewa
Qwr=4.1868·[81·C r + 246·H r - 26·(O r - Sr) - 6·Wr]
Qwr = 200145.58
2.2. Obliczenia zapotrzebowanie powietrza oraz powstałych składników spalin
teoretyczne zapotrzebowanie powietrza
Vpt = 0.089·(Cr +0.375·Sr)+ 0.265·Hr - 0.0333·Or
Vpt = 5.314
teoretyczna ilość spalin trójatomowych
VRO2 = 0.01866·(Cr+0.375· Sr) VRO2 = 0.918
teoretyczna objętość azotu w spalinach
VN2 = 0.79·Vpt+008·Nr VN2 = 5.044
teoretyczna objętość pary wodnej w spalinach
VH2O = 0.111·Hr+ 0.0124·Wr+ 0.0161·Vpt VH2O = 0.83
teoretyczna objętość tlenu
VO2:=0.21(λ - 1)Vpt VO2 = 0.223
teoretyczna objętość spalin suchych
Vsps:=VN2+VRO2+VO2
Vsps = 6.185
teoretyczna objętość spalin mokrych
Vspm:=VN2+VH2O+VRO2+VO2 Vspm = 7.016
udział objętościowy gazów trójatomowych
rR02 =
rR02 = 0.131
udział objętościowy pary wodnej
rH2O =
rH2O = 0.12
masa spalin
Gg = 1 -
1.306·nVpt Gg = 9.14
3. ENTALPIA SPALIN
dla temperatur t = 100, 200..1600
ciepło właściwe spalin
Csp = VN2· CpN2+ VH2O· CpH20+VRO2·CpRO2 + VO2· CpO2 + 0.01Ar·aun·CpAsh
Isp = Cspt
Ipow = VptCptp
tabela
T |
Isp (λ1,t) kJ/kg |
Isp (λ9,t) kJ/kg |
100 |
9,811· 102 |
1,087·103 |
200 |
1,998· 103 |
2,,202·103 |
300 |
3,022· 103 |
3,344·103 |
400 |
4,083· 103 |
4,516·103 |
500 |
5,17·103 |
5,716·103 |
600 |
6,28· 103 |
6,945·103 |
700 |
7,422· 103 |
8,202·103 |
800 |
8,589· 103 |
9,488·103 |
900 |
9,78· 103 |
1,08·104 |
1000 |
1,099·104 |
1,214·104 |
1100 |
1,223·104 |
1,35·104 |
1200 |
1,349·104 |
1,489·104 |
1300 |
1,476·104 |
1,629·104 |
1400 |
1,605·104 |
1,771·104 |
1500 |
1,735·104 |
1,914·104 |
1600 |
1,865·104 |
2,058·104 |
4. BILANS CIEPLNY KOTŁA
4.1. Straty cieplne
strata niecałkowitego spalania
Sn = 1.5 % [tab.16.7. [2]]
strata niezupełnego spalania
Sg = 0.5 % [rys.16.6. [2]]
strata promieniowania
Sp = 0.3 % [rys.16.4. [1]]
strata wylotowa
Sk =
n) Sk = 5.5 %
4.2. Sprawność kotła brutto
η = 100 - ( Sn - Sg - Sp - Sk)
η = 92.5 %
4.3. Ciepło użyteczne
entalpia pary na wylocie z kotła i7 = 3445
entalpia wody zasilającej i1 = 857.8
ilość pary świeżej Dn = 36.111
ciepło użyteczne
Quz = Dn· (i7-i1) Quz = 9.341·10
W
4.4. Zużycie paliwa w kotle
rzeczywiste zużycie paliwa
B =
·10-3 B = 5.04
teoretyczne zużycie paliwa
Bobl = B·
Bobl = 4.98
5. KOMORA PALENISKOWA
5.1. Obciążenie cieplne objętości komory paleniskowej
objętość komory paleniskowej z danych geometrycznych komory
Vk = (F2+F3+F4)⋅b+Ftr⋅a Vk = 701.6 m
jednostkowe obciążenie cieplne objętości komory paleniskowej
qv =
qv = 1.422⋅10
5.2. Obciążenie cieplne przekroju komory paleniskowej
pole przekroju komory paleniskowej :
Fk = a ⋅ b Fk = 35.5 m2
jednostkowe obciążenie cieplne poziomego przekroju komory paleniskowej
qF =
qF = 2.808 ⋅106
5.3.Określenie teoretycznej temperatury spalin w komorze paleniskowej
założenie wstępne t1 = 1915 oC
ciepło dostarczone do komory paleniskowej
Qpow = λ1⋅Vpt⋅Ctpow⋅tpowg Qpow = 1.965 ⋅ 103
entalpia teoretyczna spalin
It = Qwr⋅
+ Qpow It = 21900
temperatura adiabatyczna t1 =1849 °C
5.4.Wysokość komory paleniskowej
Hk = 20.7 m
5.5. Powierzchnia ścian komory paleniskowej
5.5.1. Ściana przednia
Ftr=
Ftr = 8.52 m
Fp = h1⋅b+ Ftr Fp =124.6 m2
5.5.2. Ściana tylna
Ft = Ftr+ e⋅b + b⋅f Ft =107.6 m2
5.5.3. Ściany boczne
F1 = 0.5⋅g⋅h F1 = 12.26 m2
F2 =
F2 = 17.88 m2
F3 =
F3 = 10.1 m2
F4 = a⋅m F4 = 80.3 m2
Fb = (F1+F2+F3+F4) Fb = 120.55 m2
5.5.4 Okno wylotowe
Fo = j⋅a Fo =33.2 m2
5.5.5. Powierzchnia palników
Fpalnik =17.4 m2
5.5.6. Powierzchnia sufitu Fsufit = 23.6 m2
Całkowita powierzchnia ścian
Fc = 2⋅Fb + Fp + Ft + Fo + Fsufit - Fpalnik Fc = 512.7 m2
5.6. Powierzchnia opromieniowana komory paleniskowej
współczynnik kształtu ekranu dla ścian x1= 1.0
współczynnik kształtu ekranu dla okna x2= 1.0
Fopr = x1⋅(2⋅Fb + Fp + Ft + Fsufit - Fpalnik) + x2⋅Fo
powierzchnia opromieniowana Fopr =512.7 m2
5.7. Temperatura spalin na wylocie z komory paleniskowej
założenie wstępne t” = 1000 ° C
współczynnik koncentracji popiołu w spalinach
współczynnik przechwytu popiołu app = 0.1
udział popiołu unoszonego w spalinach
ala =1- app ala = 0.9
μz =
μz = 0.019
współczynnik osłabienia promieniowania cząsteczkami popiołu
gęstość spalin ρsp = 1.3
średnica cząstki dz = 16 μm
kz =
= 7.495 kz = 7.495
współczynnik osłabienia promieniowania dla gazów trójatomowych
efektywna grubość warstwy promieniującej
s = 3.6⋅
s = 4.926 m
udział objętościowy gazów trójatomowych i pary wodnej
rn = rH2O + rRO2 rn = 0.251
współczynnik osłabienia promieniowania dla gazów trójatomowych
kg =
kg = 4.103
współczynnik osłabienia promieniowania cząsteczkami koksu zawartymi w spalinach
kkoks = 1.0
bezwymiarowe wielkości uwzględniające wpływ koncentracji popiołu
w spalinach
dla węgla kamiennego χ1 = 0.5
przy komorowym spalaniu paliw χ2 = 0.1
całkowity współczynnik osłabienia promieniowania
k = kg ⋅ rn + kz ⋅ μz + kkoks ⋅ χ1 ⋅ χ2 k=1.212
współczynnik efektywności cieplnej ekranów komory paleniskowej ψ
współczynnik zanieczyszczenia powierzchni opromieniowanej
dla węgla kamiennego
ξ = 0.45
współczynnik kształtu ekranu x1 = 1
ψ = x1⋅ξ ψ = 0.45
stopień czerni płomienia
dla kotłów pracujących bez nadciśnienia w komorze p = 0.1 MPa
apl =1-e-k⋅s⋅p apl = 0.4
stopień czerni komory paleniskowej
ap =
ap = 0.645
średnia efektywność cieplna spalin
Vcsr =
Vcsr =
położenie maksymalnej temperatury w komorze paleniskowej
współczynnik poprawkowy Δx = 0.1
względne położenie maksimum temperatury w komorze paleniskowej wysokość położenia osi palników hp = 6.27 m
wysokość komory paleniskowej Hk = 20.7 m
xT =
xT = 0.403
położenie maksymalnej temperatury w komorze paleniskowej
M = 0.59 - 0.5⋅xT M = 0.388
temperatura spalin na wylocie z komory paleniskowej
t'' =
- 273 t'' = 1193 oC
5.8. Sprawdzenie - kryterium Boltzmanna
zużycie paliwa
Bo =
Bo = 0.5
Bo = 0.5 < 10⋅ap=6.4
kryterium Boltzmana zostało spełnione
6. OBLICZENIA FESTONU
6.1 Charakterystyka geometryczna festonu
Powierzchnia festonu Hf = 18 m2
Podziałka poprzeczna S1f = 0.28 m
Podziałka wzdłużna S2f = 0 m
Średnica zewnętrzna rur festonu df = 0.0318 m
Wolny przekrój dla przepływu spalin Fspf = 15.854 m
6.2. Ciepło przejęte przez powierzchnię ogrzewalną
Jednostkowe obciążenie cieplne festonu (założenie wstępne)
qf = 9020
Ciepło przejęte przez powierzchnię ogrzewalną
Qf = qf ·Hf Qf = 1.624·105 W
6.3. Ilość ciepła przejęta przez czynnik
strumień pary
Df =
Df = 9.028
obliczeniowe zużycie paliw Bobl = 4.98
ciśnienie w walczaku pf =13.24 MPa
stopień suchości pary na wyjściu z festonu (założony)
x = 0.2
temperatura pary w festonie (równa temp. na linii nasycenia)
tf = 332.258 °C
entalpia pary na linii wody X=0 if` = 1543
entalpia pary na linii nasycenia X=1 if ”= 2657
entalpia pary wychodzącej z festonu
ixwyj = if` +X·(if” - if`) ixwyj = 1765
entalpia pary wchodzącej do festonu
ixwej = ixwyj- Qf ·
ixwej =1676
ilość ciepła przejęta przez czynnik
Q = 161
6.4. Ilość ciepła oddana przez spaliny przepływająca przez powierzchnię festonu
dossanie powietrza w festonie (założenie) Δnf = 0.01
temperatura zasysanego powietrza tzas = 30°C
entalpia powietrza dosysanego
Ipow :=Vpt·Cp·tzas Ipow = 211.2
entalpia spalin na wlocie na feston
isp1f = Csp·t3 isp1f = 11834
nadmiar powietrza na wylocie z festonu
λ4 = λ3+Δnf λ4 = 1.26
entalpia spalin na wylocie z festonu
isp2f = Isp1f+Δnf ·Ipow-
isp2f=11803.5
Q = φ · ( Isp1 - Isp2 + Δn · Ipow ) Q =31
6.5. Średnia temperatura między spalinami, a czynnikiem grzewczym
temperatura spalin na wlocie na feston t3 = 1.03604·103 °C
temperatura spalin na wylocie z festonu
Csp·t4-Isp2f = 0 t4 = 1.027·103°C
Δ1f = t3-tf
Δ2f = t4-tf
Średnia temperatura
Δtf =
Δtf = 699.487 oC
6.6. Obliczenia współczynnika przenikania ciepła k1
6.6.1 Współczynnik wnikania ciepła dla pęczka w układzie korytarzowym przy
przy omywaniu poprzecznym.
poprawka na ilość rur w kierunku przepływu spalin Cz = 1
średnia temperatura spalin
tfsr =
tfsr = 1032 °C
współczynnik liniowej przewodności cieplnej spalin dla temperatury średniej λ = 0.11182
liczba Prandtla dla spalin Pr = 0.576
współczynnik lepkości kinematycznej ν = 1.736·10-4
względna podziałka poprzeczna
δ1f =
δ1f = 8.805
względna podziałka wzdłużna
δ2f =
δ2f = 0
współczynnik poprawkowy uwzględniający układ geometryczny pakietu
Cs =
Cs =0.063
prędkość przepływu spalin
wf =
wf = 11.018
współczynnik wnikania ciepła dla pęczka w układzie korytarzowym przy poprzecznym omywaniu
αk = 5.155
6.6.2. Współczynnik promieniowania
1) Stopień czerni powierzchni ścian dla palenisk pyłowych aż = 0.8
2) Temperatura zewnętrzna osadu na ściance rury festonu
współczynnik zanieczyszczenia powierzchni ε1 = 0.0043
temperatura pary w festonie (const) tf = 332.3 °C
Obliczona prędkość pary
-średnica wewnętrzna rury dwf = 0.062 m
ilość równolegle włączonych wężownic zrf = 21
przekrój rur włączonych równolegle
Frf = zrf·
Frf = 0.0634 m2
objętość właściwa dla średniej temperatury pary dla przegrzewacza
Vpsrf = 0.01245
Prędkość pary
Wparyf = Df ·
wparyf = 1.773
współczynnik przejmowania ciepła na drodze konwekcji od ścianki do pary
przegrzanej
współczynnik przewodzenia ciepła dla pary przegrzanej
λparyf = 0.417
lepkość kinematyczna dla pary przegrzanej νparyf = 2.716·10-7
liczba Prandtla
Prparyf =
Prparyf = 0.552
współczynnik przejmowania ciepła na drodze konwekcji od ścianki do pary
α2 = 0.0268
α2=4.259·103
Temperatura zewnętrzna osadu na ściance rury
Tsczf = tf+
·qf + 273 Tsczf =646.162 K
3) Stopień czerni strumienia spalin
grubość warstwy promieniującej dla festonu
wysokość A = 0.0761m
podziałka D = 0.28m
długość C = 3.6m
Sgrf =
Sgrf = 0.106 m
ciśnienie w palenisku p = 0.1 MPa
współczynniki koncentracji popiołu
współczynnik przechwytu popiołu app = 0.1
udzia popiołu unoszonego w spalinachł
apl =1-app apl = 0.9
jednostkowa masa spalin Gg = 9.555
współczynnik koncentracji popiołu
μż =
μż =0.0179
współczynnik osłabienia promieniowania cząsteczkami popiołu
kż =
kż =7.373
współczynnik osłabienia promieniowania przez gazy trójatomowe kg1
kg1 =
kg1=30.789
współczynnik osłabienia promieniowania dla strumienia zapylonego
kps = (kg1 · rn + kż · μż) · p · Sgrf kps = 0.079
stopień czerni strumienia spalin
a =1-e-kps a = 0.076
Współczynnik przejmowania ciepła drogą promieniowania
αpr =5.698·10-8·
·a·(tfsr+273)3·
αpr = 16.1713
6.6.3. Współczynnik przenikania ciepła k1
współczynnik uwzględniający stopień omywania powierzchni
dla w > 4
ζ = 0.85
współczynnik wnikania ciepła
α1= ζ · αpr α1 = 13.746
współczynnik przenikania ciepła k1
k1=
k1=12.939
6.7. Sprawdzenie założonej wielkości obciążenia cieplnego festonu
założona wartość qf = 9.02·103
wartość obliczona
qf = k1·Δtf qf = 9.051·103
7. OBLICZENIA PRZEGRZEWACZA GRODZIOWEGO PG
7.1. Charakterystyka geometryczna grodzi
powierzchnia ogrzewalna grodzi H = 318 m2
podziałka poprzeczna S1 = 0.3 m
podziałka wzdłużna S2 = 0.041 m
średnica zewnętrzna rur grodzi dg = 0.0318 m
wolny przekrój dla przepływu spalin Fsp = 24.1 m
7.2. Ciepło przejęte przez powierzchnię ogrzewalną
obciążenie cieplne powierzchni (założenie wstępne) q = 30.8
ciepło przejęte przez powierzchnię ogrzewalną
Q = q⋅H Q = 9.794⋅103 MW
7.3. Ilość ciepła przejęta przez czynnik
entalpia pary wychodzącej z grodzi dla p7 = 12.9 Mpa, t7 = 540 oC
i2 = 3445
entalpia pary wchodzącej do grodzi
i1 = 3173
temperatura pary wchodzącej do grodzi dla t8 = 442.9 oC
Q = 1972
7.4. Ciepło wypromieniowane z komory paleniskowej do grodzi
entalpia spalin dla teoretycznej temperatury spalania
It = 2.19 · 104
rzeczywista entalpia spalin wylotowych z komory paleniskowej
I2 = 1.34 ·104
powierzchnia opromieniowana komory paleniskowej
Fopr = 512.7 m2
obciążenie cieplne powierzchni opromieniowanej
qHopr =
qHopr = 82.5
przekrój okna wylotowego z komory paleniskowej Hokna = 33.6 m2
ciepło wypromieniowane z komory paleniskowej do grodzi
Qpr = 0.7· qHopr·
Qpr = 389.7
7.5. Ilość ciepła oddana przez spaliny przepływające przez powierzchnię
przyrost współczynnika nadmiaru powietrza w grodzi (dossanie przez
nieszczelności) Δn = 0.05
temperatura dosysanego powietrza tzas = 30 oC
entalpia dosysanego powietrza
Ipow = Vpt⋅Cp⋅tzas Ipow = 211.2
entalpia spalin na wlocie na gródź dla
Isp1 = Csp⋅t2 Isp1 = 1.34·104
nadmiar powietrza na wylocie z grodzi
λ3 = λ2 +
λ3 = 1.25
entalpia spalin na wylocie z grodzi
Isp2 =I2 + Qpr + Δn·Ipow -
Isp2 = 1.183·104
temperatura spalin na wylocie z grodzi
Csp·t3 - Isp2 = 0 t3 =1036°C
współczynnik uwzględniający stratę promieniowania
φ = 0.95
ilość ciepła oddana przez spaliny przepływające przez powierzchnię
Q = φ · ( Isp1 - Isp2 + Δn · Ipow ) Q =1575.87
7.6. Średnia temperatura między spalinami, a czynnikiem grzewczym
Δ1 = t2 - t7 Δ1 = 653.24
Δ2 = t3 - t8 Δ2 = 593.14
Δt =
Δt = 622 °C
7.7. Współczynnik przenikania ciepła k1
7.7.1. Współczynnik wnikania ciepła od spalin do grodzi dla pęczka w układzie
korytarzowym przy omywaniu poprzecznym
poprawka na ilość rzędów rur w kierunku przepływu spalin Cz = 1
średnia temperatura spalin
temperatura spalin na wlocie do grodzi t2 = 1193 oC
temperatura spalin na wylocie z grodzi t3 = 1036 oC
tśr =
tśr = 1115 oC
współczynnik liniowej przewodności cieplnej spalin dla temperatury średniej
λ = 0.11912
liczba Prandtla dla spalin dla temperatury średniej Pr = 0.572
współczynnik lepkości kinematycznej ν = 1.915·10-4
względna podziałka wzdłużna
δ1=
δ1 = 9.434
względna podziałka poprzeczna
δ2 =
δ2 = 1.289
poprawka uwzględniająca układ geometryczny pakietu
Cs = [1 + (2 · δ1 -3) · (1-
)3]-2 Cs = 0.777
średni współczynnik nadmiaru powietrza
λsr =
λsr=1.225
prędkość przepływu spalin
w =
w =
współczynnik wnikania ciepła od spalin do grodzi αk
αk = 0.2 · Cz · Cs ·
·
·Pr0.33 αk = 49.673
7.7.2. Współczynnik promieniowania międzygrodziowego
Stopień czerni powierzchni ścian dla palenisk pyłowych aż=0.8
Temperatura zewnętrzna osadu na ścianie rury
współczynnik zanieczyszczenia powierzchni ε1 = 0.0043
średnia temperatura pary
tpsr =
tpsr =
Obliczona średnia prędkość pary
średnica wewnętrzna rury dw = 0.0192 m
ilość równolegle włączonych wężownic zr =180
przekrój rur włączonych równolegle
Frow = zr ·
Frow =
średnie ciśnienie pary w przegrzewaczu
Ppsr =
Ppsr =
prędkość pary
wpary = Dn ·
wpary =
temperatura zewnętrzna osadu na ścianie rury
Tscz = tpsr +(
Tscz =905 K
Stopień czerni strumienia spalin
grubość warstwy promieniującej dla powierzchni grodziowej
wymiary komory utworzonej przez dwie sąsiednie grodzie
wysokość A = 1.958 m
podziałka D = 0.3 m
długość C = 4.4 m
Sgr =
Sgr = 0.442 m
ciśnienie w palenisku (dla komór z nadciśnieniem) p =0.1 MPa
współczynnik koncentracji popiołu w spalinach
współczynnik przechwytu popiołu app = 0.1
udział popiołu unoszonego w spalinach
ala =1- app ala = 0.9
jednostkowa masa spalin Gg = 9.312
współczynnik koncentracji popiołu
μż =
μż = 0.01836
wspłczynnik osłabienia promieniowania cząsteczkami popiołu
kż =
kż =
współczynnik osłabienia promieniowania dla gazów trójatomowych kg1
kg1 =
kg1 = 13.817
współczynnik osłabienia promieniowania dla strumienia zapylonego
kps = (kg1 · rn + kż · uż) · p · Sgr kps = 0.115
stopień czerni strumienia spalin
a =1-e-kps a = 0.144
Współczynnik promieniowania międzygrodziowego
αpr = 5.698·10-8 · a · (tsr+273)3 ·
αpr =40.41
7.7.3. Współczynnik przenikania ciepła dla powierzchni grodziowych
współczynnik uwzględniający stopień omywania powierzchni
dla w > 4
ζ = 0.85
kątowy współczynnik grodzi dla
= 1.3 x = 0.88
współczynnik wnikania ciepła
α1 = ζ·(αpr + αk
α1 = 68.835
współczynnik przenikania ciepła k1
k1 =
k1 = 49.432
7.8. Sprawdzenie założonej wielkości obciążenia cieplnego q
wartość założona q = 3.08 ·104
wartość obliczona
q = k1 · Δt q = 3.07815 ·104
8. PRZEGRZEWACZ KONWEKCYJNY
8.1. Charakterystyka geometryczna przegrzewacza
powierzchnia ogrzewalna przegrzewacza H = 904 m2
podziałka poprzeczna S1 = 0.150 m
podziałka wzdłużna S2 = 0.075 m
średnica zewnętrzna rur Dg = 0.038 m
wolny przekrój dla przepływu spalin Fsp = 15.37 m
8.2. Ciepło przejęte przez powierzchnię ogrzewalną
- jednostkowe obciążenie cieplne przegrzewacza q = 26400
- ciepło przejęte przez powierzchnię ogrzewalną
Q = q⋅H Q = 23.87 MW
8.3. Ilość ciepła przejęta przez czynnik
8.3.1. Strumień masy pary w przegrzewaczu
strumień wtrysku pierwszego
wtr1 = 0.65 wtr 2 = 0.35
Dwtr1= wtr⋅wtr1⋅Dn Dwtr1=1.526
strumień wtrysku drugiego
Dwtr2 = wtr⋅wtr2⋅Dn Dwtr2:=0.822
strumień masy pary w przegrzewaczu
Dprzeg = Dn-Dwtr2 Dprzeg=35.29
Dprzeg1 = Dn - (Dwtr1+Dwtr2) Dprzeg1=33.764
8.3.2. Parametry wody wtryskowej
temperatura wody na wejściu tww1 = t1
tww1= 200°C
entalpia wody na wejściu Iww1 = 857.8
8.3.3. Parametry pary w przegrzewaczu
entalpia pary na wyjściu z przegrzewacza przed wtryskiem wtr2
ip =
ip=3227
entalpia pary na wyjściu z przegrzewacza za wtryskiem wtr2
i1=3173
entalpia pary wchodzącej do przegrzewacza za wtryskiem wtr1
i2 = ip-
i2=2653
entalpia pary wchodzącej do przegrzewacza przed wtryskiem wtr1
i3:=
i3=2734
ciśnienie pary na wejściu do przegrzewacza pwej=13.3 MPa
ciśnienie pary na wyjściu z przegrzewacza pwyj=13 MPa
temperatura pary wchodzącej do przegrzewacza po wtrysku wtr1
twej = 332.612°C
temperatura pary na wyjściu z przegrzewacza przed wtryskiem wtr2
twyj = 461.496°C
Q = 3892
8.4. Ilość ciepła oddana przez spaliny przepływające przez powierzchnię przegrzewacza
dossanie powietrza Δnprzeg = 0.1
temperatura dosysanego powietrza tzas = 30°C
entalpia powietrza dosysanego
Ipow = Vpt⋅Cp⋅tzas Ipow = 211.2
nadmiar powietrza na wylocie z przegrzewacza
λ5=λ4+Δnprzeg λ5 = 1.36
entalpia spalin na wlocie na przegrzewacz Isp1 = 1.18⋅104
entalpia spalin na wylocie z przegrzewacza
Isp= Isp+Δnprzeg ⋅ Ipow-
Isp2 = 7.032⋅103
temperatura spalin na wylocie z przegrzewacza
założenie wstępne t5 = 500
Csp ⋅ t5-Isp2 = 0 t5 = 603.5 °C
Q = φ · ( Isp1 - Isp2 + Δn · Ipow ) Q =4789
8.5. Temperatur średnia między spalinami, a czynnikiem grzewczym
Δ1 = t4 - twyj Δ1= 566
Δ2 = t5 - twej Δ2= 270.847
Δt =
Δt = 400.4 °C
8.6. Współczynnik przenikania ciepła. k 1
8.6.1. Współczynnik wnikania ciepła od spalin do przegrzewacza dla pęczka w układzie szachownicowym przy omywaniu poprzecznym.
poprawka na ilość rur w kierunku przepływu spalin Cz:=1
średnia temperatura spalin
tsr =
tsr = 815.5 °C
współczynnik liniowej przewodności cieplnej spalin dla temperatury średniej
λ = 0.09286
liczba Prandtla dla spalin dla temperatury średniej Pr = 0.597
współczynnik lepkości kinematycznej ν = 1.293 ⋅10-4
względna podziałka poprzeczna
δ1=
δ1 = 3.947
względna podziałka wzdłużna
δ2 =
δ2 = 1.974
δ2prim =
δ2prim= 2.791
ψ =
ψδ= 1.645
poprawka uwzględniająca układ geometryczny pakietu
dla ψ < 1.7 wynosi Cs = 0.34ψ0.1 Cs= 0.357
średni nadmiar powietrza w przegrzewaczu konwekcyjnym
λsr =
λsr= 1.31
prędkość przepływu spalin
w =
w = 9.83
współczynnik wnikania ciepła αk
αk =
αk = 87.836
8.6.2. Współczynnik przejmowania ciepła drogą promieniowania
Stopień czerni powierzchni ścian dla palenisk pyłowych aż = 0.8
Temperatura zewnętrzna osadu na ściance rury
współczynnik zanieczyszczenia powierzchni ε1 = 0.0043
średnia temperatura pary
tpsr =
tpsr = 397°C
Obliczeniowa średnia prędkość pary
średnica wewnętrzna rury dw = 0.030 m
ilość równolegle włączonych wężownic zr = 79
przekrój rur włączonych równolegle Fr = zr ⋅
objętość właściwa pary na wejściu na przegrzewacz konwekcyjny
v1 = 0.01237
objętość właściwa pary na wyjściu z przegrzewacza konwekcyjnego
v2 = 0.02255
objętość właściwa średnia dla przegrzewacza konwekcyjnego
Vpsr =
Vpsr = 0.01746
średnia prędkość pary
Wpary =Dprzeg ⋅
Wpary = 11.035
współczynnik przejmowania ciepła na drodze konwekcji od ścianki do pary
przegrzanej
średnia temperatura pary tpsr = 397.1 °C
średnie ciśnienie pary w przegrzewaczu konwekcyjnym
Ppsr =
ppsr = 13.15 MPa
współczynnik przewodzenia ciepła dla pary przegrzanej
λpary = 0.075
lepkość kinematyczna dla pary przegrzanej νpary = 4.291⋅10-7
liczba Prandtla dla pary
Prpary =
Pr pary= 1.224
współczynnik przejmowania ciepła α2 α2 = 0.0268
α2 = 3.733⋅103
temperatura zewnętrzna osadu na ściance rury
Tzsc = tpsr+
Tzsc = 790.6 °C
Stopień czerni strumienia spalin
grubość warstwy promieniującej dla powierzchni
wysokość A = 3.975 m
podziałka D = 0.15 m
długość C = 3.575
Sgr = 0.9⋅dg⋅
Sgr = 0.305
ciśnienie w palenisku p = 0.1 Mpa
współczynnik koncentracji popiołu
współczynnik przechwytu popiołu aż = 0.1
udział popiołu unoszonego w spalinach
ala = 1-aż ala = 0.9
jednostkowa masa spalin Gg = 9.902
współczynnik koncentracji popiołu μż
μż =
μż = 0.01727
d) współczynnik osłabienia promieniowania cząsteczkami popiołu
kż =
kż = 8.32
e) współczynnik osłąbienia promieniowania przez gazy trójatomowe kg1
kg1 =
kg1 = 20.944
f) współczynnik osłabienia promieniowania dla strumienia zapylonego
kps:=(kgl ⋅ rn+kż ⋅ μż) ⋅ p ⋅ Sgr kps = 0.152
Stopień czerni strumienia spalin
a = 1- e-kps a = 0.141
Współczynnik przejmowania ciepła drogą promieniowania αpr = 5.698⋅10-8 .
αpr = 24.6
8.6.3. Współczynnik przenikania ciepła k1
współczynnik kątowy dla
= 1.973 x = 0.6
współczynnik uwzględniający stopień omywania powierzchni
dla w > 4
ζ = 0.85
współczynnik wnikania ciepła
α1 = ζ⋅(αpr+αk) α1 = 95.555
współczynnik przenikania ciepła
k1 =
k1 = 66.52
8.7. Sprawdzenie wielkości obciążenia powierzchni q
wartość założona q = 2.64⋅104
wartość obliczona
q = k1 ⋅ Δt q = 2.66378⋅104
9. PODGRZEWACZ WODY
9.1. Charakterystyka geometryczna podgrzewacza wody
powierzchnia ogrzewalna podgrzewacza H = 1246 m2
podziałka poprzeczna S1 = 0.150 m
podziałka wzdłużna S2 = 0.065 m
średnica zewnętrzna rur Dg = 0.0318 m
wolny przekrój dla przepływu spalin Fsp = 14.296 m
9.2. Ciepło przejęte przez powierzchnię ogrzewalną
jednostkowe obciążenie cieplne powierzchni przegrzewacza
(założenie wstępne) q = 12720
ciepło przejęte przez powierzchnię ogrzewalną
Q = q⋅H Q = 1.585⋅104 MW
9.3. Ilość ciepła przejęta przez czynnik
strumień masy wody
Dpodg = Dprzeg1 Dpodg = 33.76
9.3.1 Parametry wody wtryskowej
temperatura na wejściu do podgrzewacza wody t1 = 200°C
ciśnienie na wejściu p1= 14.32 MPa
entalpia na wejściu iww = 857.8
ilość ciepła przejęta przez czynnik
Q =
Q= 3174
9.4. Ilość ciepła oddana przez spaliny przepływająca przez powierzchnię
podgrzewacza
przyrost współczynnika nadmiaru powietrza podgrzewaczu
Δn = 0.05
temperatura zasysanego powietrza tzas = 30°C
entalpia powietrza dosysanego Ipow :=211.2
nadmiar powietrza na wylocie z podgrzewacza
λ6:=λ5+Δnprzeg λ6=1.41
entalpia spalin na wlocie podgrzewacza Isp1 = 7032
entalpia spalin na wylocie z podgrzewacza
Isp2 = Isp+Δnprzeg ⋅ Ipow-
Isp2 = 3.86103
temperatura spalin na wylocie z podgrzewacza
założenie wstępne t6 = 500 °C
Csp ⋅t 6-Isp2 = 0 t6:=332.1 °C
entalpia wody na wyjściu z podgrzewacza wody
iwyjw = iwejw+
iwyjw =1326
ciśnienie na wylocie z podgrzewacza wody p3 = 14.1 MPa
temperatura wody na wyjściu z podgrzewacza wody twyj = 297.46 °C
Ilość ciepła oddana przez spaliny przepływająca przez powierzchnię
podgrzewacza
Q = φ · ( Isp1 - Isp2 + Δn · Ipow ) Q =3182
9.5. Temperatura średnia między spalinami, a czynnikiem grzewczym
Δ1 = t5-twyj Δ1 = 306
Δ2 = t6-t1 Δ2 = 132.07
Δt =
Δt = 207 °C
9.4. Współczynnik przenikania ciepła k 1
9.4.1.Współczynnik wnikania ciepła dla pęczka w układzie szachownicowym przy
omywaniu poprzecznym
poprawka na ilość rur w kierunku przepływu spalin Cz = 1
średnia temperatura spalin
tsr =
`tsr = 467.8 °C
współczynnik liniowej przewodności cieplnej spalin dla temperatury średniej
λ = 0.06272
liczba Prandtla dla spalin dla temperatury średniej Pr = 0.632
współczynnik lepkości kinematycznej ν = 6.786 10-5
względna podziałka poprzeczna
δ1:=
δ1 = 4.717
względna podziałka wzdłużna
δ2:=
δ2 = 2.044
δ2prim =
δ2prim = 3.121
ψδ =
ψδ = 1.752
poprawka uwzględniająca układ geometryczny pakietu
dla 1.7 < ψδ <4,5 i dla δ1 >3.0 przyjmuje się Cs = 0.34 ψδ0.2
Cs = 0.3804
średni nadmiar powietrza w podgrzewaczu wody
λsr =
λ = 1.385
prędkość przepływu spalin
w =
w = 7.575
współczynnik wnikania ciepła αk
αk =
αk = 86.999
9.4.2. Współczynnik przejmowania ciepła drogą promieniowania
Stopień czerni powierzchni ścian dla palenisk pyłowych aż = 0.8
Temperatura zewnętrzna osadu na ściance rury
współczynnik zanieczyszczenia powierzchni ε1 = 0.0043
średnia temperatura wody
twsr =
twsr = 248°C
obliczona prędkość pary
średnica wewnętrzna rury dw = 0.0254 m
ilość równolegle włączonych wężownic zr = 80
przekrój rur włączonych równolegle
Frów = zr⋅
Frów = 0.04054 m2
objętość właściwa pary na wejściu do podgrzewacza wody
v1 = 1.14389⋅10-3
objętość właściwa pary na wyjściu z podgrzewacza wody
v2 = 1.37118⋅10-3
objętość właściwa średnia dla podgrzewacza wody
Vwsr =
Vwsr =1.257535
średnia prędkość wody wwody
wwody = Dpodg⋅
wwody = 1.047
współczynnik przejmowania ciepła na drodze konwekcji
średnia temperatura wody twsr = 248.73°C
średnie ciśnienie wody w podgrzewaczu wody
pwsr =
pwody_sr:=14.21 MPa
współczynnik przewodzenia ciepła dla wody λwody:=0.636
lepkość kinematyczna νwody = 1.374⋅10-7
liczba Prandtla dla wody
Prwody =
Prwody = 0.813
współczynnik przejmowania ciepła α2
α2 = 0.0268⋅
α2 = 1.048⋅104
temperatura zewnętrzna osadu na ścianie rury
Tzsc = twsr+
Tzsc = 577.64 K
stopień czerni strumienia spalin
grubość warstwy promieniującej
wysokość A = 5.46 m
podziałka D = 0.15 m
głębokość C = 3 m
Sgr = 0.9⋅dg ⋅
Sgr = 0.323
ciśnienie w palenisku p = 0.1 MPa
współczynnik koncentracji popiołu
współczynnik przechwytu popiołu aż = 0.1
udział popiołu unoszonego w spalinach
ala = 1- aż ala = 0.9
jednostkowa masa spalin dla współczynnika nadmiaru powietrza
Gg = 10.423
współczynnik koncentracji popiołu μż
μż =
μż = 0.01641
współczynnik osłabienia promieniowania cząsteczkami popiołu kż
kż =
kż = 10.753
współczynnik osłabienia promieniowania przez gazy trójatomowe kgl
kgl =
kgl = 25.141
współczynnik osłabienia promieniowania dla strumienia zapylonego
kps = (kgl ⋅ rn+kż ⋅ μż) ⋅ p ⋅ Sgr kps = 0.184
stopień czerni strumienia
a = 1-e -kps a = 0.168
Współczynnik przejmowania ciepła drogą promieniowania αpr
αpr = 5.698⋅10-8 ⋅
αpr = 10.0387
9.4.3. Współczynnik przenikania ciepła k1
współczynnik uwzględniający stopień omywania powierzchni
dla w > 4
ζ = 0.85
współczynnik kątowy dla
= 2.004 x = 0.6
współczynnik wnikania ciepła
α1 = ζ⋅(αpr+αk) α1 = 82.48
współczynnik przenikania ciepła
k1 =
k1 = 60.535
9.5. Sprawdzenie wielkości obciążenia cieplnej powierzchni q
wartość założona q = 12.72⋅103
wartość obliczona
q = k1 ⋅ Δt q = 12.5307 103
10. PRZEGRZEWACZ OPROMIENIOWANY
10.1. Obciążenie cieplne powierzchni opromieniowanej
qHopr = 82.5 ⋅103
10.2. Powierzchnia przegrzewacza opromieniowanego
Hpo = 76.8 m2
10.3. Przyrost entalpii pary w przegrzewaczu opromieniowanym
Δipo =
Δipo = 143.6
11. BILANS KOTŁA - obliczenia sprawdzające
11.1. Przyrost entalpii pary w grodzi
entalpia pary na wejściu do grodzi idg = 3173
entalpia pary na wyjściu z grodzi izg = 3445
Δigr = izg - idg Δigr = 271.2
11.2. Przyrost entalpii pary w przegrzewaczu konwekcyjnym
entalpia pary na wejściu do grodzi idg = 3173
entalpia pary na wyjściu z grodzi izg = 3445
Δipk = izk - idk Δipk = 573.8
11.3. Przyrost entalpii pary w przegrzewaczu opromieniowanym
Δipo = 143.6
11.3. Entalpia pary na wylocie z kotła
iwyl = 3445
11.4. Entalpia pary w walczaku - na linii nasycenia X=1
io = 2657
11.5. Równania sprawdzające
strumień pary w przegrzewaczu konwekcyjnym Dprzeg1 = 33.764
strumień pary pierwotnej Dn = 36.111
iwyl ⋅ Dn - io ⋅ Dprzeg1 = Δigr ⋅ Dn + Δipp ⋅ Dprzeg1 + Δipk ⋅ Dprzeg
iwyl ⋅ Dn - io ⋅ Dprzeg1 =34.6859 MW
Δigr ⋅ Dn + Δipp ⋅ Dprzeg1 + Δipk ⋅ Dprzeg = 34.895 MW
IV LITERATURA
[1] Prof.zw.mgr inż..Teodor Wróblewski, W. Sikorski, K. Rzepa; Warszawa 1973
[2] Prof.magr inż. Piotr Orłowski, W. Dobrzański, E. Szwarc; Warszawa 1972
[3] Prof. E. Kalinowski; Politechnika Wrocławska - Wrocław 1984
„Przekazywanie ciepła”
[4] Dr inż. S. Kruczek „Kotły”
8