TOK OBLICZEŃ

Założenia i zakres parametrów

  1. 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

0x08 graphic
- 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

  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.

  1. Obliczenia komory osadczej

    1. Obliczenie prędkości opadania ziarna granicznego

Repg μg

0x08 graphic
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

0x08 graphic
Repgprzybl. =

18 + 0,61 (Arpg)0,5

Ar - liczba Archimedesa obliczana z wzoru:

dpg3 ρgp - ρg) g

0x08 graphic
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!

    1. Bilans masy pyłu dla komory osadczej, wyznaczenie skuteczności odpylania komory

mwi = mw awi gdzie: mw ↔ U - całkowity strumień masy pyłu na

wlocie do komory (unos)

mzi = mwi ηi gdzie: ηi - skuteczność przedziałowa odpylania i-tej

frakcji obliczana z wzoru

ηi = 0x01 graphic
(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.

moi = mwi − mzi

aoi = 0x01 graphic
0x01 graphic

S = 0x01 graphic

η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.

    1. Obliczenie wymiarów komory osadczej

B = H = 0x01 graphic

gdzie: Vg - rzeczywisty strumień objętości gazu, m3/s

wg - założona prędkość przepływu gazu przez komorę, m/s

L = 0x01 graphic

gdzie: upg - prędkość opadania ziarna granicznego, m/s

n = L / H jest to wielkość obliczeniowa, wyznaczana do 2 miejsc po

przecinku, teoretyczna liczba półek w komorze = (n − 1)

h = H / n

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.

    1. Obliczenie wymiarów dyfuzora i konfuzora komory

    1. Obliczenie wymiarów zasobnika pyłu pod komorą i wyznaczenie czasu gromadzenia pyłu

- strumień objętości pyłu zatrzymanego w komorze

Vpz = 0x01 graphic

gdzie: ρp us - gęstość usypowa pyłu = 1/3 gęstości właściwej pyłu

- czas gromadzenia pyłu w zasobniku

τ = 0x01 graphic

gdzie: Vzas - objętość czynna zasobnika (przestrzeni, w której może

być gromadzony pył).

  1. 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)

  1. 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).

  1. Wyznaczenie wymiarów komina

  1. 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

0x08 graphic
wg2 · ρg

gdzie: ∑ Δpm - suma oporów miejscowych ( Δpm = ζ )

2

l wg2 · ρg

0x08 graphic
0x08 graphic
∑ Δpl - suma oporów liniowych ( Δpm = λ )

dz 2

λ - współczynnik tarcia λ = 0,0032 + 0x01 graphic
Re = 0x01 graphic

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.