Ć w i c z e n i e 5
Wyznaczanie współczynnika Coriolisa
1. Wprowadzenie
Celem ćwiczenia jest doświadczalne określenie strumienia objętościowego,
prędkości średniej oraz współczynnika poprawkowego Coriolisa dla strugi powietrza
przepływającego przez przewód o kołowym przekroju poprzecznym.
Podczas przepływu płynu rzeczywistego przez przewody zamknięte jego lepkość i
związane z nią naprę\
enia styczne powodują niejednorodność rozkładu prędkości w
przekrojach poprzecznych. Prędkość maksymalna występuje w pobli\
u środka
przekroju i w sposób ciągły maleje w kierunku ścianek, osiągając na ich powierzchni
wartość równą zeru.
W obliczeniach technicznych wprowadza się zazwyczaj zało\
enie upraszczające,
polegające na przyjęciu jednorodnego rozkładu prędkości w przekroju poprzecznym,
przy czym charakterystyczna prędkość przyjmowana jest jako równa prędkości
średniej, określonej zale\
nością:
&
V
Uśr = (1)
F
w której:
&
V - strumień objętości przepływu,
F - pole przekroju poprzecznego.
Zało\
enie to pozwala na wykorzystanie w obliczeniach przewodów równania
Bernoulliego w postaci wyprowadzonej dla strugi elementarnej. Odnosząc strumień
energii do strumienia objętości, mo\
emy równanie Bernoulliego (ściślej równanie
zachowania energii) dla przekrojów kontrolnych 1 i 2 strugi elementarnej zapisać
następująco:
2 2
�U1 �U2
+ p1 + � g z1 = + p2 + � g z2 + "pstr1-2 (2)
2 2
lub w tradycyjnym ujęciu, związanym z  cię\
arowym [1] układem jednostek:
2 2
U1 p1 U2 p2
+ + z1 = + + z2 + "hstr1-2 (3)
2g �g 2g �g
Pierwsze wyrazy obydwu stron równań (2) i (3) określają energię kinetyczną w
odpowiednich przekrojach kontrolnych, drugie i trzecie odpowiednio energię ciśnienia
i energię potencjalną wysokości, a wyraz ostatni stratę energii między przekrojami.
Zało\
enie o jednorodnym rozkładzie prędkości w przekroju poprzecznym przewodu
pociąga jednak za sobą konsekwencje w postaci błędnego obliczenia strumienia
energii kinetycznej przenikającego przez ten przekrój.
Rozpatrzmy niejednorodne, ustalone w czasie pole prędkości w przekroju
poprzecznym przewodu o ścianach cylindrycznych.
Dla płynu nieściśliwego (� = const) zakładając, \
e wektory prędkości są normalne do
rozpatrywanego przekroju, strumień energii kinetycznej przenikający przez pole
elementarne dF jest równy:
39
Rys. 1. Szkic do wyznaczania rzeczywistego strumienia energii kinetycznej oraz
strumienia objętościowego przepływu
2
U
&
dE = � U dF , (4)
2
skąd po scałkowaniu otrzymujemy:
�
&
E = (5)
+"+"U 3dF
2
F
Strumień energii wyznaczony w oparciu o prędkość średnią, nazywany dalej
&
pozornym strumieniem energii Ep , wyra\
a się wzorem:
2 3
Uśr Uśr
& &
Ep = �V = �F (6)
2 2
Wykorzystując dla wyznaczenia prędkości średniej definicyjną zale\
ność (1)
przepisaną w postaci:
&
V 1
Uśr = = (7)
+"+"U dF
F F
F
pozorny strumień energii przestawić mo\
na jako:
3
�ł �ł
� 1
& �ł
Ep = (8)
+"+"U dF �ł
2
�ł �ł
2
F
�ł F łł
Nietrudno stwierdzić, \
e wyra\
enia (8) i (5) nie są jednoznaczne i zawsze spełniony
& &
jest warunek E > Ep .
&
Stosunek rzeczywistego strumienia energii kinetycznej E do strumienia
&
obliczonego z prędkości średniej Ep nazywa się współczynnikiem Coriolisa lub
rzadziej współczynnikiem Saint-Venante`a [3]:
� �
+"+"U 3dF +"+"U 3dF
2 2
F F
ą = = (9)
3
�
2
&
�ł �ł
VUśr
� 1
�ł
2
+"+"U dF �ł
2
�ł �ł
2
F
�ł F łł
a wartość jego jest większa od jedności. Energia kinetyczna występująca w równaniu
Bernoulli ego (2) lub (3), a wyra\
ona za pomocą prędkości średniej, winna być
zapisana w postaci:
40
2 2
ą�Uśr ąUśr
lub
2 2g
Rozpatrzmy często spotykany w praktyce przypadek przewodu o kołowym
przekroju poprzecznym, w którym formuje się osiowosymetryczny rozkład prędkości
U = U(r) (rys. 2). Elementarny strumień objętościowy wyniesie:
&
dV = 2Ąr drU (r) (10)
a całkowity strumień objętości:
D / 2
&
V = 2Ą (11)
+"U(r)r dr
0
Prędkość średnia określona jest wówczas wzorem:
D / 2 D / 2
2Ą rU(r)dr 8 rU(r)dr
+" +"
0 0
Uśr = = (12)
ĄD2 D2
4
Rys. 2. Rozkład prędkości w przewodzie o przekroju kołowym
Korzystając z zale\
ności (5) i (6), mo\
emy wyrazić rzeczywisty i pozorny strumień
energii kinetycznej przenikającej przez przekrój kołowy w postaci:
D / 2 D / 2
2
3
&
E = 2Ąr dr � U(r)U (r) = Ą� rU (r)dr (13)
+" +"
2
0 0
3 3
D / 2 D / 2
�ł �ł �ł �ł
� 1 64Ą�
& �ł �ł
Ep = rU(r)dr (14)
+"U(r)2Ąr dr �ł = �ł +"
2
�ł �ł �ł
2
D4 �ł 0
�ł �ł
�ł 0 łł �ł łł
ĄD2
�ł �ł
�ł �ł
4
�ł łł
Współczynnik Coriolisa w przepływie osiowosymetrycznym będzie równy:
D / 2
3
rU (r)dr
+"
&
E D4 0
ą = = (15)
3
&
Ep 64
D / 2
�ł �ł
�ł �ł
rU(r)dr
+"
�ł �ł
�ł 0 łł
41
 Sprawdz
my obecnie, jaką wartość osiąga współczynnik Coriolisa w przepływie
laminarnym i turbulentnym. Rozkład prędkości uzyskany poprzez rozwiązanie
równania Naviera-Stokesa dla przepływu laminarnego w przewodach kołowych,
wykazujący bardzo dobrą zgodność z wynikami doświadczeń, ma postać:
2 2
�ł łł �ł łł
2r 2r
�ł �ł �ł �ł
U(r) = Umax 1- �ł �ł
(16)
�ł śł H" 2Uśr 1- �ł �ł
�ł śł
D D
�ł łł �ł łł
�ł śł �ł śł
�ł �ł �ł �ł
Po wprowadzeniu zale\
ności (16) do wzoru (15) stwierdzić mo\
na, \
e w przepływach
laminarnych, dla których charakterystyczny jest paraboliczny rozkład prędkości,
wartość współczynnika Coriolisa jest stała i równa ą = 2.
Dla przepływu turbulentnego w przewodach o przekroju kołowym promieniowy
rozkład prędkości mo\
emy wyrazić za pomocą doświadczalnej formuły potęgowej
zaproponowanej przez Prandtla i Karmana [1]:
1/ n
2r
U(r) = Umax �ł1- łł (17)
�ł śł
D
�ł �ł
słusznej dla Re < 4 � 106, przy czym wartość n jest zale\
na od liczby Reynoldsa n =
f(Re) i zmienia się w granicach n = 6 � 11. Rozwiązanie całek występujących we
wzorze (15) pozwala uzyskać dla przepływu turbulentnego zale\
ność:
(2n +1)3(n +1)3
ą = (18)
4n4(n + 3)(2n + 3)
z której wynika, \
e w tym typie przepływu współczynnik Coriolisa jest równie\

zale\
ny od liczby Reynoldsa.
Przykładowe wartości n i ą dla kilku liczb Reynoldsa zestawiono poni\
ej:
Re = 4 � 103; n = 6; ą = 1,08
Re = 1,1 � 104; n = 7; ą = 1,06
Re = 3,2 � 106; n = 10; ą = 1,03.
W praktyce podczas obliczania przepływów turbulentnych w przewodach długich
opuszczamy zwykle współczynnik Coriolisa przyjmując, \
e jest on równy jedności.
Niedokładności popełnione przy określaniu strat (liniowych i lokalnych) przewy\
szają
na ogół nieścisłości wynikające z zało\
enia ą = 1.
2. Stanowisko pomiarowe
Schemat stanowiska pomiarowego pokazano na rys. 3. Wentylator promieniowy 1
o regulowanej prędkości obrotowej tłoczy powietrze do komory uspokajającej 2, w
której dzięki systemowi prostownic i siatek strumień powietrza zostaje
ujednorodniony. Z komory przez odpowiednio ukształtowaną dyszę powietrze
przepływa do przewodu 3 o średnicy D. W wybranym przekroju kontrolnym tego
przewodu jest dokonywany pomiar rozkładu prędkości. Rurka Pitota ciśnienia
całkowitego 4 zamocowana w suporcie 5 umo\
liwiającym jej promieniowy przesuw,
mierzy ciśnienie całkowite, natomiast impuls ciśnienia statycznego jest odbierany z
otworka wykonanego w ściance przewodu. Przewody impulsowe ciśnień całkowitego i
statycznego są połączone z ró\
nicowym mikromanometrem pochylnym 6, którego
wskazanie pozwala obliczyć ciśnienie dynamiczne w wybranym punkcie
pomiarowym.
42
Rys. 3. Schemat stanowiska pomiarowego
3. Metodyka pomiarów i obliczeń
Z rozdziału 1 wynika, \
e dla wyznaczenia prędkości średniej i współczynnika
Coriolisa w przepływie osiowosymetrycznym niezbędna jest znajomość rozkładu
prędkości w kierunku promieniowym. Prędkość w danym punkcie określić mo\
na
znając wartość ciśnienia dynamicznego, z zale\
ności:
2 pd
U = , m/s (19)
�
w której �  gęstość przepływającego czynnika, kg/m3.
W układzie pomiarowym opisanym powy\
ej, ciśnienie dynamiczne wyznacza się ze
wzoru:
pd = �m g l i, N/m2 (20)
w którym:
�m - gęstość cieczy manometrycznej, kg/m3,
g - przyspieszenie ziemskie, m/s2,
l - długość słupa cieczy manometrycznej, m,
i - przeło\
enie mikromanometru.
Znajomość doświadczalnie określonego rozkładu prędkości U(r) pozwala na
wyznaczenie wartości całek oznaczonych, występujących w zale\
nościach (11), (12) i
(15):
D / 2
!1 = (21)
+"U(r)r dr
0
D / 2
!2 = (r)r dr (22)
+"U 3
0
przez zastosowanie jednej z metod całkowania przybli\
onego, np. metody prostokątów
zilustrowanej na rys. 4.
Promień przewodu dzielimy na n równych przedziałów o szerokości:
D
"r = , m (23)
2n
43
Rys. 4. Ilustracja metody wyznaczania przybli\
onej wartości całek
a wartość całek ! i ! wyznaczamy ze wzorów:
1 2
n
!1 = "r rk (24)
"U
k
k =1
n
!2 = "r rk (25)
"U 3
k
k=1
w których:
k - numer kolejnego przedziału,
rk - promień środka k-tego przedziału, m,
Uk - prędkość w środku k-tego przedziału, m/s.
Wzór (11) określający strumień objętości przepływu przyjmie postać:
&
V = 2Ą !1, m3/s (26)
Prędkość średnia zgodnie z (12) wyra\
a się jako:
8
Uśr = !1, m/s (27)
D2
a współczynnik Coriolisa określony z zale\
ności (15) przyjmie postać:
D4 !2
ą = (28)
64
(!1)3
Liczbę Reynoldsa charakteryzującą badany przepływ obliczamy z zale\
ności:
Uśr D
Re = (29)
�
w której �  kinematyczny współczynnik lepkości, m2/s.
4. Szczegółowy program ćwiczenia
Po uruchomieniu tunelu aerodynamicznego i ustaleniu warunków jego pracy,
nale\
y dokonać sondowania pola prędkości wzdłu\
promienia przewodu. Dla
uproszczenia póz
niejszych obliczeń współrzędne promieniowe punktów pomiarowych
nale\
y dobrać tak, aby le\
ały one w środkach przedziałów o szerokości "r:
r1 = 0,5"r
44
r2 = (2 - 0,5) "r
& & & & & & & & ......
rk = (k - 0,5) "r
Wyniki wskazań mikromanometru lk wpisujemy w tablicę pomiarowo-obliczeniową i
dla ka\
dego z punktów pomiarowych obliczamy:
- ciśnienie dynamiczne pdk  wzór (20),
- prędkość przepływu Uk  wzór (19),
- iloczyn rkUk,
- iloczyn rkUk3,
a następnie określamy wartości sum:
n n
3
"r Uk ; "r Uk
k k
k =1 k =1
które pomno\
one przez "r pozwalają zgodnie z (24) i (25) obliczyć wartości całek !
1
oraz ! . Korzystając z zale\
ności (26), (27), (28), (29) obliczamy kolejno:
2
&
- strumień objętości przepływu V , m3/s
- prędkość średnią Uśr, m/s
- współczynnik Coriolisa ą,
- liczbę Reynoldsa Re.
Po wykonaniu pomiarów i obliczeń nale\
y sprawdzić, czy uzyskana wartość
współczynnika Coriolisa zawiera się w zakresie charakterystycznym dla liczby
Reynoldsa badanego przepływu.
Literatura:
1. Bukowski J.: Mechanika płynów, PWN, Warszawa 1959
2. Czetwertyński E., Utrysko B.: Hydraulika i hydromechanika, PWN, Warszawa 1969
3. Walden H., Stasiak J.: Mechanika cieczy i gazów w in\
ynierii sanitarnej, Arkady,
Warszawa 1971
45
Tabela pomiarowo-obliczeniowa
D = & & & .m; t = & & & .oC; � = & & .......kg/m3; �m = & & & & kg/m3;
� = & & ..m2/s; n = & & ..; "r = & & & .m;
n n
&
!1 = "r Uk = ..........; !2 = "r U = ..........; V = 2Ą !1 = ..........;
"r "r 3
k k k
k=1 k=1
8 D4 !2 Uśr D
Uśr = !1 = ...........; ą = = .........; Re = = .........
3
64 �
D2 !1
rk lk pdk Uk rkUk rkUk3
L.p.
m m N/m2 m/s m2/s m4/s3
20
"
k=1
46