TOK OBLICZEŃ
Założenia i zakres parametrów
Ziarna mają kształt kulisty,
Ziarna nie zderzają się ze sobą ani ze ściankami komory,
Przepływ gazu jest jednorodny i jednakowy w całym przekroju komory,
Prędkość gazu w przewodach 15-18 (20) m/s
Prędkość gazu w komorze osadczej 0,5 - 1,5 m/s
Średnica ziarna granicznego 60 lub 40 µm
Przeliczenie wszystkich parametrów na warunki rzeczywiste
- temperatura
Trz = t + 273, K
gdzie: t - temperatura, oC;
- strumień objętości gazów
o o
Vrz = Vo Trz / To
gdzie: Vo - strumień objętości gazów w warunkach normalnych, m3/h lub m3/s
To = 273 K
- stężenie pyłu
S = So To/Trz
gdzie: So - stężenie pyłu w warunkach normalnych, g/m3
(1 + x) p
- gęstość gazu w warunkach normalnych ρgo =
(0,622 + x) Rw To
gdzie: x - zawartość wilgoci w gazach (podana w temacie), kg/kgpow. suchego
Rw - stała gazowa pary wodnej, Rw = 461,5 kJ/(kg K)
- gęstość gazów w warunkach rzeczywistych
ρrz = ρgo To/Trz
- dynamiczny współczynnik lepkości
μrz = μo (To + C) /(Trz + C) (Trz/To)1,5
gdzie: μo = 17,08 x 10-6 kg/(m s)
C = 112
- określenie składu granulometrycznego pyłu na wlocie do odpylacza (udziału
masowego poszczególnych frakcji ziarnowych pyłu ai, na podstawie
podanego w temacie udziału masowego skumulowanego Aj).
Lp |
dg μm |
Zakres średnic μm |
dś μm |
Aj % |
ai |
1 2 3 … n-1 n |
< 2 < 5 < 10 ... < 150 > 150
|
0 - 2 2 - 5 5 - 10 … 100 - 150 > 150 |
1 3,5 7,5 … 125 150 |
5 11 22 … … 100 |
0,05 0,06 0,11 … … … Σ = 1 |
Obliczenie niezbędnej skuteczności odpylania całej instalacji
ηc = (U - E) / U
gdzie: U, E - odpowiednio unos i emisja, g/s, obliczane z wzorów:
U = Spw Vg E = Spo Vg
gdzie: Spw i Spo - stężenie pyłu odpowiednio na wlocie i wylocie instalacji, g/m3
Vg - strumień objętości gazów, m3/s
Uwaga: obie wielkości podawane są w tych samych warunkach - albo normalnych albo
rzeczywistych.
Obliczenia komory osadczej
Obliczenie prędkości opadania ziarna granicznego
Repg μg
upg =
dpg ρg
gdzie: dpg - założona średnica ziarna granicznego, zwykle 60 μm (40 μm),
μg,, ρg - dynamiczny współczynnik lepkości i gęstość gazu, w
warunkach rzeczywistych
Repgprzybl. - przybliżona liczba Reynoldsa (dla ziarna granicznego), obliczana wg wzoru:
Arpg
Repgprzybl. =
18 + 0,61 (Arpg)0,5
Ar - liczba Archimedesa obliczana z wzoru:
dpg3 ρg (ρp - ρg) g
Ar =
μg2
ρp - gęstość pyłu, kg/m3,
g - przyspieszenie ziemskie.
Ostateczną liczbę Reynoldsa określa się na podstawie Reprzybliżonej:
- dla ruchu laminarnego opadania ziarna Rep ≤ 2 Re = Ar / 18
- dla ruchu przejściowego 2 > Rep ≤ 500 Re = 0,152 Ar0,715
- dla ruchu burzliwego Rep > 500 Re = 1,74 Ar0,5
i tę wartość wstawia się do wzoru na prędkość opadania ziarna!
Bilans masy pyłu dla komory osadczej, wyznaczenie skuteczności odpylania komory
Strumień masy pyłu i-tej frakcji na wlocie do komory
mwi = mw awi gdzie: mw ↔ U - całkowity strumień masy pyłu na
wlocie do komory (unos)
Strumień masy pyłu i-tej frakcji zatrzymany w komorze
mzi = mwi ηi gdzie: ηi - skuteczność przedziałowa odpylania i-tej
frakcji obliczana z wzoru
ηi =
(stosunek prędkości opadania ziarna i-tej frakcji do prędkości
opadania ziarna granicznego)
Prędkości opadania ziaren wszystkich frakcji obliczane są dla średniej średnicy ziarna danej frakcji.
Strumień masy pyłu i-tej frakcji na wylocie komory (opuszczający komorę)
moi = mwi − mzi
Udział masowy i-tej frakcji ziarnowej pyłu opuszczającego komorę
aoi =
Stężenie pyłu w gazie opuszczającym komorę
S =
Całkowita skuteczność komory osadczej
ηk = ∑ (awi ηi)
W pracy należy zamieścić przykład obliczeń dla wybranej frakcji pyłu (korzystnie dla frakcji o średnicy ziarna poniżej granicznego), pozostałe wyniki zestawić w tabeli.
Nr frakcji |
Rozmiar ziaren µm |
dśri
µm |
Ari |
Repi |
Rei |
upi
m/s |
ηi |
awi |
ηi awi |
mwi
g/s |
mzi
g/s |
moi
g/s |
aoi |
1 2 3 4 5 … n |
0 - 2 2 - 5
|
|
|
|
|
|
0 0 |
|
0 0 |
|
0 0 |
|
|
|
|
|
|
|
|
|
|
∑=1 |
∑=ηko |
∑=U |
∑ |
∑ |
∑=1 |
Prędkości opadania ziaren < 5 µm są na tyle małe, że można założyć skuteczność przedziałową odpylania tych frakcji ziarnowych równą 0.
W zaznaczonych kolumnach należy zsumować obliczone wartości.
Obliczenie wymiarów komory osadczej
Obliczenie szerokości i wysokości komory
B = H =
gdzie: Vg - rzeczywisty strumień objętości gazu, m3/s
wg - założona prędkość przepływu gazu przez komorę, m/s
Obliczenie długości komory bezpółkowej
L =
gdzie: upg - prędkość opadania ziarna granicznego, m/s
Obliczenie liczby przestrzeni międzypółkowych
n = L / H jest to wielkość obliczeniowa, wyznaczana do 2 miejsc po
przecinku, teoretyczna liczba półek w komorze = (n − 1)
Obliczenie wysokości przestrzeni międzypółkowej
h = H / n
Wyznaczenie skorygowanych rzeczywistych wymiarów komory:
Lrz = L / n (1,1 - 1,15) dolicza się 10-15% rezerwy
Hrz = H + (n − 1) 0,002 zwiększa się o grubość blachy półek (2 mm
każda)
Brz = B + ∑przestrzeni zsypowych szerokość pojedynczej przestrzeni
zsypowej 40 lub 50 mm.
Obliczenie wymiarów dyfuzora i konfuzora komory
Obliczenie wymiarów zasobnika pyłu pod komorą i wyznaczenie czasu gromadzenia pyłu
- strumień objętości pyłu zatrzymanego w komorze
Vpz =
gdzie: ρp us - gęstość usypowa pyłu = 1/3 gęstości właściwej pyłu
- czas gromadzenia pyłu w zasobniku
τ =
gdzie: Vzas - objętość czynna zasobnika (przestrzeni, w której może
być gromadzony pył).
Wyznaczenie niezbędnej skuteczności odpylacza w drugim stopniu
- skuteczność odpylaczy połączonych szeregowo
ηc = 1 − (1 − ηI) (1 − ηII)
gdzie: ηc - skuteczność całkowita odpylania wyliczana z wartości unosu i emisji
ηI - skuteczność odpylania odpylacza w I stopniu (komory osadczej)
Dobór odpylacza II stopnia
W zależności od wymaganej skuteczności będzie to:
- bateria cyklonów (ηII < 89%)
- multicyklony (ηII = 90-95%) lub
- filtr włókninowy pulsacyjny (dla ηII > 95%).
Multicyklony i filtry dobiera się w zależności od strumienia odpylanych gazów oraz całkowitej skuteczności odpylania, wyznaczanej ze skuteczności przedziałowej, określanej na podstawie nomogramów, i udziału masowego danej frakcji ziarnowej wg wzoru:
ηcyk = ∑ (awi ηi)..
Dobór cyklonów na podstawie materiałów dostarczonych przez prowadzących wymaga przeanalizowania 5-6 wariantów, wg załączonej tabeli:
Nr wariantu |
Liczba cyklonów n |
Średnica nominalna D m |
Prędkość gazu wg m/s |
Skuteczność
ηcyk |
Opory przepływu Δp Pa |
|
|
|
|
|
|
Skuteczność baterii cyklonów (ηcyk) oblicza się dla każdego wariantu analogicznie jak dla komory osadczej: ηcyk = ∑ (awi ηi).
Wyznaczenie wymiarów komina
Średnice komina dobiera się tak aby prędkość przepływu gazów mieściła się w zakresie 10-12 m/s,
Wysokość komina - przyjmuje się o 5 m większą od najwyższego punktu instalacji (komora osadcza lub urządzenie w II st. odpylania).
Dobór wentylatora
Przed doborem wentylatora należy wykonać schematy (rzut z góry i widoki boczne, tak aby pokazać każdy element instalacji), w celu wyznaczenia oporów przepływu. Całkowite opory przepływu (Δpc) oblicza się wg wzoru:
Δpc = ∑ Δpm + ∑ Δpl + ∑ Z − Δpkom
wg2 · ρg
gdzie: ∑ Δpm - suma oporów miejscowych ( Δpm = ζ )
2
l wg2 · ρg
∑ Δpl - suma oporów liniowych ( Δpm = λ )
dz 2
λ - współczynnik tarcia λ = 0,0032 +
Re =
dz - średnica zastępcza przewodu,
l - długość odcinka przewodu (można sumować długości przewodów o takich
samych wymiarach)
∑ Z - suma oporów na urządzeniach; dla komory osadczej z dyfuzorem i konfuzorem
przyjąć 100 Pa
Δpkom - ciąg kominowy wspomagający pracę wentylatora obliczany z wzoru:
Δpkom = h · g (ρpow - ρg)
gdzie: h - geometryczna wysokość komina,
g - przyspieszenie ziemskie,
(ρpow - ρg) - różnica gęstości powietrza na wylocie komina (można przyjąć
1,25 kg/m3) i gazów odlotowych.
Wentylator dobiera się wg katalogu na obliczone opory przepływu (zwiększone o 10%) oraz rzeczywisty strumień spalin.