mpw_9.pdf

(195 KB) Pobierz
Microsoft Word - roz-dziewiec.doc
9.
Kryteria podobieństwa przepływów.
Analiza podobieństwa jest jednym z najczęściej stosowanych narzędzi w mechanice
płynów i w rozdziałach poprzednich wielokrotnie używaliśmy tej metody, chociaż nie
wprowadziliśmy jeszcze tego pojęcia. Przykładowo, wykresy współczynników strat tarcia i
strat lokalnych sporządziliśmy w funkcji liczby Reynoldsa, przyjmując intuicyjnie, że wartość
Re a nie U uzasadnia przyjęcie do obliczeń danej wartości λ lub ξ . Najważniejszym
powodem stosowania analizy podobieństwa jest złożoności opisu ruchu płynu i duży nakład
obliczeń numerycznych potrzebnych do rozwiązania równań N-S czy też Eulera. Podobnie
pracochłonny i kosztowny jest w mechanice płynów eksperyment i stąd też celowym jest
stosowanie metod analizy, pozwalających przenosić uzyskane numerycznie lub
eksperymentalnie rozwiązania na inne przypadki, jeżeli tylko spełnione będą kryteria
podobieństwa tych przepływów. Przykładowo, załóżmy że znamy rozwiązanie opisujące pole
prędkości wokół profilu aerodynamicznego o charakterystycznej skali liniowej l , który
opływany jest strumieniem płynu o charakterystycznej prędkości napływu U , jak pokazano
na rys. 9.1a.
a)
b)
U 1
U 2
l 1
Rys. 9.1. Przepływ dla którego znane jest rozwiązanie a) i przepływ podobny b) dla
którego poszukujemy sposobu przeniesienia rozwiązania.
Wiedząc, że przepływ z rys. 9.1b charakteryzuje się skalą liniową l i opływany jest z
prędkością U interesuje nas, czy możliwe jest przeniesienie znanego rozwiązania na ten
przypadek i w jaki sposób należy dokonać przeskalowania istniejącego rozwiązania.
Zagadnienie to jest właśnie przedmiotem analizy podobieństwa.
9.1. Klasyfikacja kryteriów podobieństwa.
Analiza podobieństwa przepływów oparta jest o kryteria, których spełnienie
powoduje, że rozwiązanie z jednego przypadku może być w odpowiedni sposób przeniesione
na inny przypadek. Rozwiązaniem może być kształt linii prądu (pole przepływu) a także
współczynnik oporu opływanego ciała, czy też współczynnik straty lokalnej zaworu.
Traktując dwa przepływy jako dwa pola fizyczne będziemy mogli zastosować kryteria
podobieństwa mechanicznego zakładające, że jeżeli w każdej odpowiadającej sobie parze
punktów dwóch pól fizycznych, stosunki charakteryzujących je jednorodnych wielkości
fizycznych będą jednakowe, wówczas pola te będą podobne .
Podstawowymi jednostkami fizycznymi wykorzystywanymi w analizie podobieństwa są
długość, czas i siła a stosunki tych wielkości są ilorazami podobieństwa tychże wielkości.
Najprostszym rodzajem podobieństwa jest podobieństwo geometryczne zakładające, że
między dwoma polami fizycznymi spełniony jest warunek stałości odpowiadających sobie
wymiarów liniowych:
139
146615935.011.png
l
1
=
idem
(9.1)
l
2
gdzie l i l są charakterystycznymi wymiarami obydwu pól fizycznych.
Bardziej zaawansowanym rodzajem podobieństwa mechanicznego jest podobieństwo
kinematyczne które wymaga, aby oprócz warunku (9.1) spełniony był też wymóg stałości
ilorazów charakterystycznych czasów:
t
1
=
idem
(9.2)
t
2
Podobieństwo kinematyczne zakłada zatem, że w obydwu przepływach elementy płynu w
odpowiadających sobie punktach przebywają odcinki podobnych geometrycznie dróg
S i
S
w czasach t i t , które spełniają warunek (9.2).
Z warunku podobieństwa kinematycznego wynika związek:
S
1
=
U
1
t
1
=
l
1
=
idem
S
U
t
l
2
2
2
2
który po przekształceniu daje warunek:
U
1
=
l
1
t
2
=
idem
U
l
t
2
2
1
wykazujący, że podobieństwo kinematyczne jest podobieństwem pól prędkości.
Najważniejszym rodzajem podobieństwa mechanicznego jest podobieństwo
dynamiczne wymagające aby w odpowiadających sobie punktach dwóch przepływów
podobnych kinematycznie, odbywających się w konfiguracjach podobnych geometrycznie
spełnione były warunki stałości ilorazów podobieństwa wartości sił działających na elementy
płynu przy zachowaniu identyczności kierunków działania odpowiadających sobie sił. Jeżeli
zatem dwa przepływy mają być podobne dynamicznie, wówczas wieloboki sił działających na
odpowiadające sobie elementy płynu muszą być podobne. Oznacza to z kolei, że ilorazy
odpowiednich sił składowych muszą być identyczne i jednakowe muszą być kierunki
działania tych sił.
9.2. Bezwymiarowe równanie ruchu.
Jak wynika z analizy przeprowadzonej w rozdz. 3, w przepływach płynu
rzeczywistego występuje pięć rodzajów sił, które w dalszej części rozważań oznaczać
będziemy następująco:
- siły bezwładności B
- siły ciężkości G
- siły ciśnieniowe P
- siły lepkości L
- siły ściśliwości S .
Równania ruchu zgodnie z prawem Newtona wyrażają następujący warunek równowagi:
(wszystkich sił) = 0
który musi być spełniony w każdym z punktów przepływu. Jeżeli rozważać będziemy dwa
przepływy oznaczone odpowiednio 1 i 2 , wówczas w każdym z punktów przepływu 1
spełniony będzie warunek:
B
1
+
G
1
+
P
1
+
L
1
+
S
1
=
0
i podobnie w przepływie 2 :
++
W każdej odpowiadającej sobie parze punktów obydwu przepływów podobnych dynamicznie
spełniona być musi następująca relacja:
B
2
G
2
P
2
+
L
2
+
S
2
=
0
B
1
=
G
1
=
P
1
=
L
1
=
S
1
(9.3)
B
G
P
L
S
2
2
2
2
2
140
146615935.012.png 146615935.013.png 146615935.014.png
stanowiąca warunek całkowitego podobieństwa dynamicznego obydwu przepływów. Jeżeli
między tymi przepływami ma być spełniony związek (9.3), wówczas nasuwa się wniosek, że
rozwiązanie równania ruchu winno być wyrażone w taki sposób, aby przez jego odpowiednie
przeskalowanie możliwe było uzyskanie rozwiązań dla wszystkich dynamicznie podobnych
przepływów. Najprostszym sposobem uzyskania takiego rozwiązania jest przekształcenie
równania ruchu do postaci bezwymiarowej, w której wszystkie wielkości fizyczne będą
odniesione do charakterystycznych wymiarów.
a)
b)
c)
U
0
U
0
U
śr
l
l
l
U 0
Rys. 9.2.
Przykłady charakterystycznych skal dla wybranych przepływów.
Na rys. 9.2 pokazano przykłady przepływów z zaznaczonymi charakterystycznymi
wymiarami liniowymi l i charakterystycznymi prędkościami U , które nazywamy skalami
przepływu. Dla opływu profilu aerodynamicznego (rys. 9.2a) charakterystyczną skalą liniową
jest cięciwa profilu a dla opływu kolumny (rys. 9.2b) jej średnica. W obydwu przypadkach
charakterystyczną skalą prędkości U jest prędkość przepływu niezakłóconego, określona w
miejscu, gdzie linie prądu nie są zdeformowane obecnością opływanego ciała. W przepływie
płynu rzeczywistego w kanale (rys. 9.2c) charakterystyczną skalą prędkości jest prędkość
średnia, natomiast skalą liniową jest średnica kanału.
Jeżeli przyjmiemy następujące skale:
- dla czasu t
- dla wymiarów liniowych l
- dla prędkości U
wówczas zmienne występujące w równaniach ruchu będziemy mogli przekształcić do
następującej postaci bezwymiarowej:
ξ
=
x
;
η
=
y
;
ζ
=
z
l
l
l
U
=
U
x
;
U
=
U
y
;
U
=
U
z
x
y
z
U
U
U
o
o
o
U
(9.4)
U
=
U
o
p
=
p
;
ρ
_
=
ρ
;
τ
=
t
p
ρ
t
o
o
o
Podstawienie bezwymiarowych zmiennych do równania N-S dla kierunku x daje:
141
_
146615935.001.png 146615935.002.png 146615935.003.png 146615935.004.png 146615935.005.png
U
o
U
x
+
U
2
o
U
U
x
+
U
U
x
+
U
U
x
=
x
y
z
t
l
o
=
g
X
p
o
1
p
+
ν
U
o
2
U
+
ν
U
o
div
U
x
g
l
ρ
ρ
l
3
l
2
o
gdzie:
2
2
2
2
=
+
+
2
2
2
div
=
+
+
Po pomnożeniu obydwu stron równania przez
U o
i wykonaniu identycznych przekształceń dla kierunków y oraz z otrzymujemy
bezwymiarowe równanie N-S w postaci wektorowej:
l
l
d
U
gl
F
p
1
=
o
grad
()
p
+
U
t
d
τ
U
2
o
g
ρ
U
2
o
ρ
o
o
o
(9.5)
+
ν
2
U
+
ν
grad
div
U
lU
3
lU
o
o
Dla podobieństwa dynamicznego przepływów wystarczy identyczność współczynników przy
odpowiednich członach równania, co daje następujące liczby podobieństwa:
l
=
1
;
St – liczba Strouhala
(9.6)
U
t
St
o
o
gl
2
o
=
1
;
Fr – liczba Froude’a
(9.7)
U
Fr
p
o
=
1
1
;
M - liczba Macha
(9.8)
ρ
U
2
o
κ
M
2
o
ν
=
1
;
Re – liczba Reynoldsa
(9.9)
lU o
Re
Jeżeli w przepływie 1 występować będą następujące wielkości charakterystyczne:
U
01
;
l
1
;
t
01
;
g
;
p
01
;
ρ
01
;
1
a w przepływie 2 będą one równe:
U ρ
wówczas warunki podobieństwa dynamicznego obydwu tych przepływów będą następujące:
02
;
l
2
;
t
02
;
g
;
p
02
;
02
;
2
l
1
=
l
2
=
1
U
t
U
t
St
01
01
02
02
gl
1
=
gl
2
=
1
U
2
01
U
2
02
Fr
p
01
=
p
02
=
1
1
ρ
U
2
01
ρ
U
2
02
κ
M
2
01
02
ν
1
=
ν
2
=
1
l
U
l
U
Re
1
01
2
02
W praktyce niemożliwe jest jednoczesne spełnienie wszystkich tych warunków i zmuszeni
jesteśmy ograniczyć się do korzystania z podobieństwa częściowego . Należy wówczas wybrać
z kryteriów danych wz. (9.6) ÷ (9.9) ten warunek, który dotyczy najbardziej istotnej siły
142
146615935.006.png
występującej w przepływie. Konieczna jest tu jednak wiedza dotycząca interpretacji liczb
podobieństwa i zagadnienie to zostanie omówione w rozdziale następnym.
9.3. Sens fizyczny liczb podobieństwa.
Zastosujmy analizę wymiarową dla zjawisk tarcia zachodzących przy przepływie
płynu wzdłuż ściany. Elementarna siła lepkości wynosi:
dS
dL τ=
(9.10)
a poszczególne człony aproksymowane być mogą następująco:
τ
=
µ
U
~
µ
dU
n
dl
(9.11)
S
~
dl
2
dL µ
a po sformułowaniu tej zależności dla elementu płynu o wymiarach skończonych:
l
~
dL µ
~
U
(9.12)
Elementarna siła bezwładności zapisana być może jako:
a
dB
=
dm
gdzie:
dm
ρ
dl
3
a
dU
dt
a po przejściu do elementu płynu o wymiarach skończonych:
2
dB ρ (9.13)
Uwzględniając wynik powyższej analizy, liczba Reynoldsa dana wz. (9.9) zapisana być może
jako:
~
U
2
l
Re
=
dB
=
U
l
dL
ν
co oznacza, że liczba Reynoldsa jest ilorazem siły bezwładności i lepkości stanowiąc
kryterium podobieństwa w przepływach, gdzie te dwie siły są dominujące. Stwierdzenie to
zilustrujemy na przykładzie badań modelowych pokazanych na rys. 9.3, gdzie obiektem
rzeczywistym jest profil skrzydła samolotu o cięciwie wynoszącej (rys. 9.3a):
[ ]
l =
5
m
a)
b)
U
01
U
01
5m
1m
Rys. 9.3.
Geometria profilu rzeczywistego a) i jego modelu b).
Jeżeli samolot porusza się z prędkością:
[ ] [ ]
U
=
850
km
/
h
=
236
m
/
s
w powietrzu, którego lepkość wynosi:
ν
=
15
10
[ ]
m
2
/
s
143
gdzie dl - jest elementarnym wymiarem liniowym.
Podstawienie zal. (9.11) do wz. (9.10) daje następujący wynik analizy wymiarowej:
dUdl
146615935.007.png 146615935.008.png 146615935.009.png 146615935.010.png
Zgłoś jeśli naruszono regulamin