Fragment notatki:
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 l1 , który
opływany jest strumieniem płynu o charakterystycznej prędkości napływu U1 , jak pokazano
na rys. 9.1a.
a)
b)
U1
U2
l1
l2
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 2 i opływany jest z
prędkością U 2 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
l1
= idem
(9.1)
l2
gdzie l1 i l 2 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:
t1
= idem
(9.2)
t2
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 S1 i S2
w czasach t1 i t 2 , które spełniają warunek (9.2).
Z warunku podobieństwa kinematycznego wynika związek:
S1
U1t1
l
=
= 1 = idem
S2
U2t 2
l2
który po przekształceniu daje warunek:
U1
lt
= 1 2 = idem
U2
l 2 t1
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:
B1 + G1 + P1 + L1 + S1 = 0
i podobnie w przepływie 2 :
B2 + G 2 + P2 + L 2 + S2 = 0
W każdej odpowiadającej sobie parze punktów obydwu przepływów podobnych dynamicznie
spełniona być musi następująca relacja:
B1 G1 P1 L1 S1
=
=
=
=
(9.3)
B2 G 2 P2 L 2 S2
140
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)
U0
c)
U0
Uśr
l
Rys. 9.2.
l
l
U0
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 o , 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 o 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 o
dla wymiarów liniowych l
dla prędkości U o
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
Uy
U
U
Ux = x ; Uy =
; Uz = z
Uo
Uo
Uo
→
→
U
_
p=
(9.4)
U
=
p
po
Uo
_
; ρ=
ρ
ρo
; τ=
t
to
Podstawienie bezwymiarowych zmiennych do równania N-S dla kierunku x daje:
141
2
Uo ∂ U x Uo
∂Ux
∂Ux
∂Ux
+ Uz
+ Uy
+
Ux
=
t o ∂τ
l
∂ζ
∂η
∂ξ
=g
→
X p o 1 ∂ p νU o 2
ν Uo ∂
div U
−
+
∇ Ux +
g lρo ρ ∂ξ
l
3 l 2 ∂ξ
gdzie:
∂2
∂2
∂2
∇
=
+
+
∂ξ 2 ∂η 2 ∂ξ 2
∂
∂
∂
div =
+
+
∂ξ ∂η ∂ξ
Po pomnożeniu obydwu stron równania przez
2
Uo
l
i wykonaniu identycznych przekształceń dla kierunków
bezwymiarowe równanie N-S w postaci wektorowej:
2
l
→
y
oraz
z
otrzymujemy
→
d U gl F
p 1
= 2
− o 2 grad (p ) +
U o t o dτ U o g ρ o U o ρ
(9.5)
→
→
ν
ν
grad div U
+
∇2 U +
lU o
3lU 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)
Uo t o
St
gl
1
=
;
Fr – liczba Froude’a
(9.7)
2
Fr
Uo
po
1 1
=
; M - liczba Macha
(9.8)
2
κ M2
ρo U o
ν
1
=
;
Re – liczba Reynoldsa
(9.9)
lUo
Re
Jeżeli w przepływie 1 występować będą następujące wielkości charakterystyczne:
U 01 ; l1 ; t 01 ; g ; p 01 ; ρ01 ; ν1
a w przepływie 2 będą one równe:
U 02 ; l 2 ; t 02 ; g ; p 02 ; ρ02 ; ν 2
wówczas warunki podobieństwa dynamicznego obydwu tych przepływów będą następujące:
l2
l1
1
=
=
U 01 t 01
U 02 t 02
St
gl1
2
U 01
gl 2
1
=
2
Fr
U 02
p 01
p 02
1 1
=
=
2
2
κ M2
ρ01 U 01
ρ02 U 02
ν1
ν2
1
=
=
l1 U 01
l 2 U 02
Re
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
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:
dL = τ dS
(9.10)
a poszczególne człony aproksymowane być mogą następująco:
dU
∂U
~ µ
τ = µ
∂n
dl
(9.11)
2
S ~ dl
gdzie dl - jest elementarnym wymiarem liniowym.
Podstawienie zal. (9.11) do wz. (9.10) daje następujący wynik analizy wymiarowej:
dL ~ µdUdl
a po sformułowaniu tej zależności dla elementu płynu o wymiarach skończonych:
dL ~ µ U l
(9.12)
Elementarna siła bezwładności zapisana być może jako:
dB = dm ⋅ a
gdzie:
dm ≅ ρ dl3
dU
a ≅
dt
a po przejściu do elementu płynu o wymiarach skończonych:
dB ~ ρ U 2 l 2
(9.13)
Uwzględniając wynik powyższej analizy, liczba Reynoldsa dana wz. (9.9) zapisana być może
jako:
dB
U⋅l
Re =
=
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 ]
b)
a)
U01
U01
1m
5m
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 −6 m 2 / s
[
143
]
wówczas liczba Reynoldsa wynosi:
Re1 = 7.9 ⋅107
Jeżeli model takiego skrzydła o cięciwie (rys. 9.3b) wynoszącej:
l 2 = 1 [m]
będziemy badać w tunelu aerodynamicznym, wówczas dla zachowania podobieństwa sił
oporu (czyli tarcia) konieczne będzie zachowanie identyczności liczby Reynoldsa, tzn.:
U 2 ⋅ l2
Re 2 =
= Re1
ν
skąd otrzymamy:
U 2 = 1180 [m / s]
Oznacza to, że prędkość przepływu w tunelu będzie kilkakrotnie przekraczać prędkość
dźwięku i wyniki takich badań będą całkowicie niewiarygodne. Należy zatem przeprowadzić
te badania w płynie, którego lepkość byłaby dobrana tak, aby zapewnić żądaną wartość liczby
Reynoldsa bez przekraczania prędkości dźwięku.
Dla zilustrowania sensu fizycznego kolejnej liczby podobieństwa zapiszmy
elementarną siłę ciężkości:
dG = dm ⋅ g ~ ρ dl3 ⋅ g
która dla elementu płynu o wymiarach skończonych będzie proporcjonalna do:
dG ~ ρ ⋅ g ⋅ l3
(9.14)
Liczbę Froude’a zapisaną wz. (9.7) przedstawić możemy następująco:
dG
gl
Fr =
=
dB
U2
co oznacza, że stanowić ona będzie kryterium podobieństwa dla zjawisk, w których role
decydującą odgrywają siły ciężkości i bezwładności. Załóżmy, że badania modelowe będą
tym razem dotyczyć oporu kadłuba statku, pokazanego na rys. 9.4.
a)
b)
U2
U1
l1
Rys. 9.4.
l2
Geometria statku a) i jego modelu b).
Dla prawidłowego modelowania oporu tarcia kadłuba konieczne będzie zachowanie równości
liczb Reynoldsa:
U1 l1
U 2 l2
Re1 =
=
= Re 2
ν
ν
co zostało przedyskutowane w przykładzie poprzednim. W trakcie ruchu statku zachodzą
jednak dodatkowe opory, których mechanizm powstawania pokazano na rys. 9.5. Poruszający
się kadłub przyspiesza wodę, która na skutek działania sił bezwładności unosi się ku górze,
pokonując działanie sił ciężkości. Oznacza to, że dla zachowania podobieństwa oporu
falowego konieczne jest zachowanie identyczności liczb Froude’a:
2
U1
U2
2
Fr1 =
=
= Fr2
gl1
gl 2
Prawidłowe zamodelowanie obydwu składowych oporu wymagałoby zatem przeprowadzenia
dwóch eksperymentów, w których należałoby zachować identyczność liczb Reynoldsa (opór
tarcia) i Froude’a (opór falowy).
144
B
G
Rys. 9.5.
Mechanizm powstawania oporu falowego.
Interpretacja liczby Macha:
p
1 1
=
2
κ M2
ρU
wymaga wprowadzenia pojęcia prędkości dźwięku równej:
p
a =
κ
ρ
co pozwala zapisać:
U
(9.15)
M =
a
co oznacza, że liczba Macha jest ilorazem prędkości przepływu i prędkości dźwięku. Będzie
to zatem liczba podobieństwa dla przepływów, w których istotną rolę odgrywają zjawiska
ściśliwości (przepływy okołodźwiękowe i naddźwiękowe).
Liczba Strouhala dana związkiem:
Uo t o
(9.16)
St =
l
jest kryterium podobieństwa dla przepływów nieustalonych, w których przepływ zmienia się
z czasem.
t0
Rys. 9.6.
Przepływ nieustalony w maszynie przepływowej.
Przykład takiego przepływu, pokazany na rys. 9.6 pokazuje, że ważnym parametrem jest tu
okres zmian pola prędkości za okresowo przesuwającymi się łopatkami wirnikowymi. Dla
przepływu takiego badania modelowe wymagać będą zachowania identyczności:
St1 = St 2
i wówczas możemy oczekiwać, że wyniki otrzymane z eksperymentu modelowego będą
wiernie odzwierciedlać charakterystyki obiektu rzeczywistego.
145
... zobacz całą notatkę
Komentarze użytkowników (0)