Część 1.Wyznaczenie krytycznej liczby Reynoldsa, przy której następuje przejście z ruchu laminarnego w turbulentny.
Wstęp teoretyczny
W przewodach i w bocznicach kopalnianej sieci wentylacyjnej, mogą zachodzić dwa rodzaje przepływów powietrza: laminarny i turbulentny.
Przepływ laminarny odbywa się warstwowo, przy czym oddzielne warstwy powietrza przesuwają się równolegle do osi przewodu nie mieszając się między sobą. Ze względu na lepkość największa prędkość przepływu jest odnotowana w środku wyrobiska. W kopalniach przepływ laminarny może zachodzić w szczelinach górotworu, otamowanych zrobach i polach pożarowych, podsadzce suchej itp.
W przepływie turbulentnym elementy płynu poruszają się w sposób nieuporządkowany i po bardzo zawiłych torach, wskutek czego powstają ciągłe chaotyczne zaburzenia przepływu.
Dla stwierdzenia, z jakim przepływem mamy do czynienia wyznacza się liczbę Reynoldsa:
(1)
gdzie:
- prędkość średnia powietrza, m/s,
- średnica ekwiwalentna (zastępcza, równoważna) wyrobiska górniczego, przy czym
(2)
- lepkość kinematyczna powietrza kopalnianego; = 15⋅10-6 m2/s,
A - pole przekroju poprzecznego wyrobiska, m2,
B - obwód wyrobiska, m.
W praktyce rozróżnia się pierwszą (dolną) krytyczną liczbę Reynoldsa Rekr1 (charakteryzuje utratę stateczności laminarnego przepływu powietrza, dla wyrobisk górniczych Rekr1 = 1000÷1500 ) i drugą (górną) liczbę Reynoldsa Rekr2 (po przekroczeniu której występuje statecznie burzliwy przepływ powietrza, Rekr2=50000÷80000). Przy przepływie powietrza przez długi gładki przewód kołowy krytyczna liczba Reynoldsa wynosi ok 2300.
Liczba oporu, a co za tym idzie, dyssypacja energii, jest zależna od liczby Reynoldsa i chropowatości względnej przewodu. W przedziale liczb Re odpowiadającym przepływowi laminarnemu nie zależy od chropowatości ścian wyrobiska, dyssypacja energii jest wtedy liniowo zależna od prędkości przepływu powietrza.. W przypadku przewodów chropowatych dla liczb Re z przedziału od około 4000 (przy dużych ) do około 630000 (przy małych) zależy od i Re. Dla wyższych Re nie zależy od Re, dyssypacja energii jest proporcjonalna do przepływu w drugiej potędze.
Przepływy dzielą się ponadto na stacjonarne (niezależące od czasu) i niestacjonarne (zależne od czasu)
Przebieg ćwiczenia
W przekroju II (D=292mm) umieszczono rurkę dymną i powoli zwiększano prędkość przepływu powietrza, za pomocą termoanemometru zmierzono prędkość w przekroju I (D=100mm), przy której w strużce dymu w przekroju II pojawiają się turbulencje, co oznaczało przekroczenie krytycznej liczby Reynoldsa..
Dane techniczne stanowiska pomiarowego:
Rys. 1. Schemat stanowiska
Dane techniczne stosowanych przyrządów
Termoanemometr
Wyniki pomiarów
Pomiar prędkości wykonano 3 razy:
wmax1_I =1,6 m/s
wmax2_I =1,6m/s
wmax3_I =1,6m/s
stąd: wmsr_I =1,6m/s
Tok obliczeń
Obliczono przybliżoną wartość liczby Reynoldsa, następnie lgRe, z nomogramu odczytano stosunek prędkości średniej do maksymalnej i obliczono prędkość średnią w przekroju I:
lgRe=4,028
wn_I/wmax_I =0,78
wn_I =0,78wmax_I =1,248m/s
Obliczono strumień objętości powietrza (jednakowy dla rur o przekroju I i II, co wynika z prawa ciągłości przepływu) i wyznaczono prędkość średnią w rurze o przekroju II:
Obliczono krytyczną liczbę Reynoldsa:
Dyskusja błędów i wnioski
Chwila przejścia z ruchu laminarnego laminarnego w turbulentny była wyznaczana na podstawie subiektywnej oceny dymu, co obarcza pomiar błędem - otrzymany wynik był wyższy od teoretycznej wartości krytycznej liczby Reynoldsa (Rekr=2300). Prawdopodobnie turbulencje pojawiły się w strużce dymu już wcześniej, jednak były one na tyle małe, że nie zostały zauważone.
Część II. Wyznaczenie rozkładu prędkości w przekroju poprzecznym rury o średnicy 292
Wstęp teoretyczny
Rozkład prędkości w ruchu laminarnym stanowi paraboloidę, co spowodowane jest charakterem tego przepływu - tarcie między ściankami przewodu a warstwą powietrza płynącą najbliżej ścianki powoduje jej hamowanie, a ta z kolei wskutek lepkości hamuje dalej położone warstwy. W ruchu turbulentnym górę nad siłami lepkości biorą siły bezwładności, profil prędkości ulega spłaszczeniu. Początkowo w pobliżu ścianek przewodu występuje warstwa laminarna, w której występują duże różnice prędkości. Grubość tej warstwy w miarę wzrostu liczby Reynoldsa zmniejsza się, by całkowicie zaniknąć przy w pełni rozwiniętym przepływie turbulentnym.
Przebieg ćwiczenia
Zmierzono prędkość przepływu powietrza w 15 punktach, poczynając od środka przekroju, kończąc na brzegu rury.
Dane techniczne stanowiska pomiarowego
Rys. 2. Schemat stanowiska do badania profilów prędkości powietrza
Dane techniczne przyrządów
Prędkość mierzono termoanemometrem
Wyniki pomiarów
Tabela nr 1. Wyniki pomiarów cz. II
l.p. |
w [m/s] |
1 |
1 |
2 |
1 |
3 |
1 |
4 |
1 |
5 |
1 |
6 |
0,9 |
7 |
0,9 |
8 |
0,8 |
9 |
0,8 |
10 |
0,8 |
11 |
0,8 |
12 |
0,8 |
13 |
0,8 |
14 |
0,8 |
15 |
0 |
Na podstawie wyników pomiarów sporządzono dołączony do sprawozdania wykres.
Wnioski i dyskusja błędów
Najwyższą prędkość odnotowano w środku przekroju, jednak różnice prędkości w różnych punktach nie są zbyt duże, gdyż przepływ ma charakter turbulentny. Pomierzone prędkości często powtarzają się, następnie następuje skok prędkości - jest to spowodowane niską dokładnością urządzenia, w rzeczywistości przejścia między różnymi prędkościami (dla ruchu turbulentnego - średnimi) są ciągłe.
Część III. Wyznaczenie depresji naturalnej w oczku sieci wentylacyjnej
Wstęp teoretyczny
Depresja naturalna jest to praca wszystkich czynników naturalnych (takich jak zmiany cieplne w prądzie powietrza, zmiana jego składu, działanie wiatru w atmosferze zewnętrznej, porywanie cząstek powietrza przez spadającą wodę) odniesiona do 1kg lub 1m3 powietrza. Największe znaczenie ma depresja cieplna - powstająca wskutek zmian cieplnych w powietrzu kopalnianych.
Dane techniczne stanowiska pomiarowego.
Rys. 3. Schemat stanowiska pomiarowego
Dane techniczne stosowanych przyrządów
Termoanemometr AM-4204
A.Wyznaczenie depresji naturalnej metodą pomiaru wydatku przy zatrzymanym wentylatorze
Najpierw obliczono się opór sieci wentylacyjnej - włączono wentylator, zmierzono strumień powietrza przepływający przez sieć wentylacyjną i spiętrzenie wentylatora, na podstawie których wyznaczono opór sieci. Następnie wyłączono wentylator i włączono grzałkę; po ustabilizowaniu się temperatury wyznaczono strumień objętości powietrza wywołany działaniem depresji naturalnej.
Wyniki pomiarów:
wmax=3,7 m/Max
l=0,898mm (przynachyleniu 1/15)
Po włączeniu grzałki:
t=39,06°C
w=0,5 m/s
Tok obliczeń
Obliczono podciśnienie w punkcie 6 (
), równe spiętrzeniu całkowitemu wentylatora (
):
Obliczono gęstość powietrza:
Obliczono przybliżoną wartość liczby Reynoldsa, następnie lgRe, z nomogramu odczytano stosunek prędkości średniej do maksymalnej i obliczono prędkość średnią:
lgRe=4,857
wn /wmax =0,83
wn =0,83wmax =3,071m/s
Obliczono spiętrzenie całkowite wentylatora (
):
Obliczono strumień objętości powietrza:
Obliczono wartość oporu sieci:
Obliczono prędkość średnią:
Re=9730
lgRe=3,99
wsr/wmax=0,78
wsr=0,39 [m/s]
Obliczono strumień powietrza:
Obliczono depresję:
B.Wyznaczenie depresji naturalnej Metodą W. Budryka
Po ustabilizowaniu się temperatury powietrza w kilku punktach zmierzono temperaturę powietrza i jego podciśnienie w stosunku do otoczenia, zmierzono również parametry powietrza wlotowego: temperaturę na termometrze suchym, wilgotnym i ciśnienie statyczne bezwzględne. Wyznaczono gęstość powietrza dla wszystkich punktów, wyniki naniesiono na wykres p=f( 1/ρ). Wielkość pola otrzymanej figury jest miarą depresji naturalnej.
Wyniki pomiarów:
Tabela nr 2. Wyniki pomiarów cz. 3B
Nr punktu |
1 |
2 |
3 |
4 |
5 |
6 |
l [m] (nach. 1/25) |
|
0,042 |
0,069 |
0,082 |
0,082 |
0,082 |
wmax [m/s] |
|
3,6 |
3,6 |
3,6 |
3,5 |
3,6 |
Ts [K] |
300,15 |
303,05 |
301,65 |
305,15 |
309,05 |
310,45 |
Z [m] |
125,2 |
122,87 |
122,76 |
123,25 |
124,09 |
125,1 |
Ts w laboratorium: 27°C
Tw w laboratorium: 17,1 °C
Ciśnienie w laboratorium: 748,9 Tr = 98778,57 Pa
Tok obliczeń:
Obliczono podciśnienie w punktach 2-6
Obliczono ciśnienie barometryczne:
Obliczono ciśnienie bezwzględne statyczne powietrza:
Obliczono gęstość powietrza:
Tabela nr 3. Wyniki obliczeń cz. 3B
Nr punktu |
1 |
2 |
3 |
4 |
5 |
6 |
pc |
|
13,51 |
22,20 |
26,39 |
26,38 |
26,38 |
pb |
98778,57 |
98804,53 |
98805,88 |
98800,14 |
98790,69 |
98779,65 |
p |
98778,57 |
98791,01 |
98783,68 |
98773,76 |
98764,31 |
98753,27 |
ρ |
1,147 |
1,136 |
1,141 |
1,128 |
1,113 |
1,108 |
1/ρ |
0,872 |
0,881 |
0,877 |
0,887 |
0,898 |
0,902 |
Wykres 1. Zależność p=f( 1/ρ).
Pole figury zaznaczonej na wykresie obliczono za pomocą programu AutoCAD, wynosi ono 0,4557, stąd h=0,456 Pa
C. Wyznaczenie depresji naturalnej metodą bezpośrednią
Pomierzono różnicę ciśnień na szczelnej tamie (7-8) zatrzymującej w sieci przepływ powietrza wywołany działaniem depresji naturalnej.
Wyniki pomiarów
l=0,014m przy nachyleniu 1/25
Tok obliczeń
Obliczono różnicę ciśnień:
D. Obliczenie depresji naturalnej metodą zredukowanego przepływu
Dokonano pomiaru strumieni objętości powietrza (za pomocą termoanemometru umieszczonego w środku rurociągu) oraz spiętrzenia wentylatora przy zasuwie całkowicie zamkniętej i częściowo zamkniętej.
Wyniki pomiarów
Tabela nr 4. Wyniki pomiarów w cz 3D
Położenie |
1 |
2 |
L (nach 1/25) |
0,085 |
0,025 |
nach |
0,04 |
0,04 |
wmax |
3,6 |
1,8 |
Tok obliczeń
Metodą punktową obliczono średnią prędkość powietrza
Obliczono strumień powietrza dla obu ustawień zasuwy
Obliczono opór i depresję naturalną rozwiązując układ równań
Tabela nr 5. Obliczenia cz 3D
Położenie zasuwy |
1 |
2 |
Re |
70 080 |
35 040 |
lgRe |
4,8456 |
4,5446 |
wn/wmax |
0,83 |
0,81 |
wm |
2,988 |
1,458 |
V |
0,200 |
0,098 |
pc |
27,35028 |
8,0442 |
h=1,238 [Pa]