ROZDZ8B, Zbigniew Kosma Podstawy Mechaniki Płynów


Strumień ciepła doprowadzonego do elementu płynu może być skutkiem reakcji chemicznych, przewodnictwa cieplnego, promieniowania, jonizacji i innych procesów. Uwzględnimy tylko strumień ciepła dostarczony drogą przewodnictwa cieplnego; jego zmiana w kierunku osi x (rys. 8.5) wynosi

0x01 graphic

Składowa gęstości strumienia ciepła jest określona zgodnie z prawem Fouriera wzorem

0x01 graphic
(8.23)

gdzie jest współczynnikiem przewodzenia ciepła, zależnym - w ogólnym przypadku - od temperatury i ciśnienia.

Analogicznie otrzymuje się zmiany strumienia ciepła w pozostałych kierunkach osi współrzędnych, co pozwala na zapisanie strumienia ciepła 0x01 graphic
w postaci równania

0x01 graphic
(8.24)

Wynikiem połączenia wszystkich uzyskanych rezultatów jest równanie zachowania energii

0x01 graphic

0x01 graphic
(8.25)

które - po podstawieniu zależności dla naprężeń (8.4), (8.5), (8.11) i wykorzystaniu oznaczeń (3.25) - można przedstawić następująco

0x01 graphic

0x01 graphic
(8.26)

Równanie (8.26) możemy przekształcić w sposób czysto formalny i wprowadzić entalpię (7.7). Po wykorzystaniu tożsamości

0x01 graphic
(8.27)

oraz równania ciągłości (3.18)

0x01 graphic
(8.28)

mamy

0x01 graphic

0x01 graphic
(8.29)

Mnożąc kolejne równania ruchu (8.13) przez składowe prędkości 0x01 graphic
oraz dodając je stronami uzyskamy równanie określające zmianę energii wewnętrznej płynu 0x01 graphic
Jeśli równanie to odejmiemy następnie stronami od równania energii (8.25), to otrzymamy równanie

0x01 graphic
, (8.30)

w którym jest składnikiem reprezentującym dyssypację energii mechanicznej

0x01 graphic

0x01 graphic

0x01 graphic
(8.31)

tzn. tę część pracy sił powierzchniowych, która nie ulega przekształceniu na energię kinetyczną płynu, ale na jego energię wewnętrzną, powodując podwyższenie temperatury płynu.

Równanie (8.30), przy wykorzystaniu zależności (8.28) i (7.7), można zapisać również w postaci

0x01 graphic
(8.32)

Z kolei zajmiemy się zmianami entropii zachodzącymi w elemencie gazu. Opierając się na związkach (7.8) i (7.11) mamy

0x01 graphic
(8.33)

Z porównania zależności (8.32) i (8.33) otrzymamy wzór

0x01 graphic
(8.34)

z którego wynika, że zmiany energii płynu zależą wyłącznie od efektów dyssypatywnych, tzn. od lepkości i przewodności cieplnej. Oznacza to, że entropia poruszającego się elementu płynu jest stała, jeśli płyn jest nielepki i nie przewodzi ciepła.

8.4. Podstawowe zagadnienie mechaniki p*yn*w

Podstawowe zagadnienie mechaniki płynów polega na wyznaczeniu ruchu płynu w sąsiedztwie zadanego ciała i sprowadza się do rozwiązania układu równań różniczkowych wynikających z praw zachowania, z odpowiednimi warunkami brzegowymi i początkowymi.

Układ równań dla płynu lepkiego i ściśliwego, doskonałego w sensie termodynamicznym, składa się z równania ciągłości (3.18), trzech równań ruchu (8.13), równania zachowania energii (8.26) i układu równań (1.13) ÷ (1.15). Możliwe jest zatem wyznaczenie sześciu niewiadomych funkcji: 0x01 graphic
- przy założeniu, że znane są zależności wyznaczające lepkość płynu μ i przewodność cieplną λ - przykładem takiej zależności jest wzór (1.8).

Warunki początkowe , formułowane tylko dla niestacjonarnych ruchów płynu, charakteryzują stan ruchu płynu oraz jego stan fizyczny w pewnej chwili uznanej umownie za chwilę początkową badanego zjawiska. Dla musimy więc określić wartości każdej spośród funkcji niewiadomych, np.:

0x01 graphic
(8.35)

Warunki brzegowe odnoszą się natomiast do brzegu 0x01 graphic
obszaru przepływu i określają wartość każdej funkcji niewiadomej lub pochodnej tej funkcji w każdym punkcie brzegu obszaru, w dowolnej chwili

W odniesieniu do prędkości płynu warunki brzegowe postulują zazwyczaj „przyklejanie się” płynu do powierzchni ciał stałych (rozdz. 1.3); jeśli powierzchnia brzegowa jest nieruchoma wtedy

0x01 graphic
(8.36)

Warunki brzegowe dla temperatury T mogą dotyczyć albo samej temperatury, albo jej gradientu, zależnie od tego czy powierzchnia ciała stałego ma zadaną temperaturę, czy też znany jest strumień ciepła transportowanego przez rozważaną powierzchnię. Oprócz tego zwykle zakładamy jednorodność przepływu w nieskończoności; prędkość, ciśnienie, gęstość i temperatura w nieskończoności powinny więc spełniać następujące warunki:

0x01 graphic
(8.37)

8.5. Formy opisu ruchu cieczy lepkiej

Cieczą lepką nazywać będziemy płyn newtonowski o stałej gęstości, stałej lepkości μ oraz stałej przewodności cieplnej λ. Przy tych założeniach układ równań, składający się z równania ciągłości (3.20) i równania Naviera-Stokesa (8.17):

0x01 graphic
(8.38)

można rozwiązać niezależnie od równania energii (8.26). Oznacza to, że pole temperatury T jest wyznaczane dopiero po określeniu pola prędkości 0x01 graphic
i pola ciśnienia p. Temperatura odgrywa więc rolę podrzędną w ruchu cieczy lepkiej, nie wpływa bowiem na prędkość i ciśnienie, co jest konsekwencją przyjęcia stałości μ oraz λ.

Rozwiązywanie układu równań (8.38) dla zmiennych fizycznych: 0x01 graphic
napotyka na szereg trudności, gdyż równanie ciągłości ma istotnie odmienną budowę od równania Naviera-Stokesa. Z tego też względu układ równań (8.38) jest zastępowany często innymi równaniami równoważnymi albo też stosowane są metody obliczeniowe oparte na wykorzystaniu zmiennych Lagrange'a.

Działając operatorem rotacji na wirowość (3.29)

0x01 graphic

i wykorzystując równanie ciągłości otrzymamy równanie

0x01 graphic
(8.39)

które może być użyte zamiast równania ciągłości. Jeśli pole sił masowych jednostkowych 0x01 graphic
jest potencjalne, to taka sama operacja zastosowana do równania Naviera-Stokesa zezwala na uzyskanie równania Helmholtza

0x01 graphic
(8.40)

którego lewa strona wynika z przekształcenia zależności (4.4) za pomocą operatora rotacji, po podstawieniu: 0x01 graphic
(przykł. 4.8)

0x01 graphic

Równania (8.39) i (8.40) tworzą układ równań, który musi być uzupełniony nieznanymi warunkami brzegowymi dla wektora wirowości 0x01 graphic
- wyznaczanymi w trakcie obliczeń. Ponadto, w otrzymanym układzie równań, przybliżone wartości składowych prędkości mogą nie spełniać równania ciągłości. Tę niedogodność układu (8.39) ÷ (8.40) można usunąć wprowadzając potencja* wektorowy 0x01 graphic
0x01 graphic
zdefiniowany wzorem

0x01 graphic
(8.41)

równanie ciągłości jest w tym przypadku spełnione tożsamościowo, ponieważ

0x01 graphic

Istniejąca pewna dowolność w określaniu 0x01 graphic
- wzór (12.31) - pozwala na przyjęcie założenia, że 0x01 graphic
tworzy pole solenoidalne

0x01 graphic
(8.42)

i wtedy z (8.41) mamy

0x01 graphic
(8.43)

Równania (8.40) i (8.43) tworzą zamknięty układ równań dla wektorów 0x01 graphic
i 0x01 graphic
po uprzednim wykorzystaniu zależności (8.41).

Inne możliwości zmodyfikowania układu równań (8.38) mogą polegać na wykorzystaniu zamiast równania ciągłości równania Poissona dla ciśnienia, otrzymanego jako wynik działania operatora diwergencji na obie strony równania Naviera-Stokesa, lub też na przyjęciu koncepcji sztucznej ściśliwości, która może być zastosowana przy wyznaczaniu stacjonarnych przepływów cieczy lepkiej. Zakładając zmienność gęstości według zależności

0x01 graphic
(8.44)

gdzie 0x01 graphic
jest gęstością cieczy nieściśliwej, - małym parametrem, możemy równanie ciągłości (3.18) zapisać następująco

0x01 graphic
(8.45)

przy pominięciu członów 0x01 graphic
Biorąc pod uwagę dodatkowo związek (7.18), mały parametr ε możemy zastąpić również prędkością c rozchodzenia się nieskończenie słabych zaburzeń; otrzymamy

0x01 graphic
(8.46)

*

Wielkie trudności matematyczne związane z wyznaczaniem rozwiązań pełnego układu (8.38) lub też układów równoważnych, zachęcały do badania szczególnych klas przepływów, dla których pewne składniki równania Naviera-Stokesa stają się pomijalnie małe, tak że można je odrzucić i uzyskać tą drogą równania zlinearyzowane, łatwiejsze do rozwiązania.

Jedna z takich uproszczonych form została zaproponowana przez Stokesa dla przypadku powolnego ruchu ciała w cieczy o dużej lepkości, gdy siły spowodowane lepkością są znacznie większe od sił bezwładności - wobec czego te ostatnie można pominąć. Zakładając 0x01 graphic
i ruch stacjonarny otrzymamy równanie Stokesa

0x01 graphic
(8.47)

W zagadnieniu opisywanym układem równań (8.38a) i (8.47) tkwi jednak swego rodzaju wewnętrzna sprzeczność. Okazuje się bowiem, że odrzucone z założenia siły masowe w równaniu rządzącym przepływem, stają się - w dostatecznie dużej odległości od ciała opływanego - większe od sił lepkości; ich odrzucenie jest więc usprawiedliwione tylko w obszarach bliskich ścianek ciała stałego. Ponadto okazuje się (przykł. 8.6), że zagadnienie to dla płaskiego opływu ciała w ogóle nie ma rozwiązania - paradoks Stokesa (przy założeniu, że w obszarze rozwiązania nie ma żadnych osobliwości).

Paradoks Stokesa nie wystąpi, jeśli w równaniach wyjściowych zachowane będą pewne składniki całkowitej siły masowej, które w dużej odległości od ciała są porównywalne z siłami lepkości. Tego typu linearyzację równania Naviera-Stokesa zaproponował Oseen postulując alternatywne równanie, nazywane równaniem Oseena , rządzące przepływem cieczy o dużej lepkości

0x01 graphic
(8.48)

gdzie 0x01 graphic
jest prędkością strumienia jednorodnego, równoległego do osi x.

Przybliżenie nieliniowe równań Naviera-Stokesa, oparte na koncepcji tzw. warstwy przyściennej, jest omawiane w rozdziale dziewiątym.

*

Postać przedstawionych równań opisujących przestrzenny ruch cieczy lepkiej ulega znacznemu uproszczeniu przy ograniczeniu się do rozwiązania płaskiego ruchu cieczy lepkiej

Stosunkowo najmniejszemu uproszczeniu ulega układ (8.38), który w formie rozwiniętej - po pominięciu pola 0x01 graphic
sił masowych jednostkowych - przedstawia się następująco:

0x01 graphic
(8.49)

Równanie Poissona dla ciśnienia otrzymamy po zróżniczkowaniu równania (8.49b) względem x oraz równania (8.49c) względem y, dodaniu ich stronami i odpowiednim wykorzystaniu równania ciągłości (8.49a)

0x01 graphic

0x01 graphic
(8.50)

Dzięki znikaniu składowych potencjału wektorowego 0x01 graphic
i składowych rotacji prędkości 0x01 graphic
w kierunku osi x i w kierunku osi y, bardzo dużemu uproszczeniu ulegają równania (8.40) i (8.43) - otrzymujemy:

0x01 graphic
(8.51)

gdzie a jest funkcją prądu (6.7).

Płaski, stacjonarny ruch cieczy lepkiej może być więc opisany - w szczególności - układem równań:

0x01 graphic
(8.52)

lub też równoważnym quasi-liniowym równaniem czwartego rzędu dla funkcji prądu

0x01 graphic
(8.53)

8.6. Hydrodynamiczna teoria smarowania

Hydrodynamiczna teoria smarowania objaśnia zjawiska występujące w cienkiej warstwie oleju - zwanej filmem olejowym - wypełniającej szczelinę pomiędzy dwiema gładkimi, poruszającymi się względem siebie powierzchniami. Zjawiska takie obserwowane są w łożyskach ślizgowych maszyn wirnikowych i tłokowych, w przekładniach zębatych i wielu innych węzłach tarcia różnych maszyn - pracujących w mgle olejowej.

0x01 graphic

Rys. 8.6

W celu wyprowadzenia podstawowych równań hydrodynamicznej teorii smarowania obieramy układ współrzędnych 0x01 graphic
- w taki sposób, żeby powierzchnie ograniczające przepływ (0x01 graphic
na rys. 8.6) mogły być opisane równaniami:

0x01 graphic
(8.54)

Opis taki umożliwia określenie wielkości szczeliny między powierzchniami za pomocą różnicy tych funkcji

0x01 graphic
(8.55)

oraz wprowadzenie pomocniczego układu współrzędnych

0x01 graphic
(8.56)

Niech powierzchnie ograniczające przepływ: 0x01 graphic
poruszają się z prędkościami, odpowiednio,  i .

Ze względu na małą grubość warstwy smarnej (rzędu mikrometrów) występuje w niej laminarny ruch oleju, ale panujące w obszarze styku ciał bardzo wysokie ciśnienia (rzędu megapascali) powodują konieczność uwzględnienia zarówno ściśliwości oleju, jak i silnej zależności jego lepkości od ciśnienia. Przepływ oleju w szczelinie między współpracującymi elementami maszyn musi więc być opisany najogólniejszymi równaniami mechaniki płynów: równaniem ciągłości (3.16), równaniami Naviera-Stokesa (8.13) i równaniem energii (8.26). Ponieważ jednak grubość filmu olejowego jest bardzo mała w porównaniu z powierzchniami węzłów tarcia, możemy wprowadzić założenia upraszczające - wynikające z faktu, że zmiany parametrów hydrodynamicznych w kierunku normalnym do powierzchni  i powierzchni są dużo większe od zmian w kierunkach stycznych. Przy tych założeniach równania Naviera-Stokesa przepisujemy w postaci:

0x01 graphic
(8.57)

i następnie pomijamy wszystkie siły masowe jako małe w porównaniu z siłami ciśnieniowymi i siłami tarcia:

0x01 graphic
(8.58)

Po uśrednieniu temperatury i lepkości oleju w każdym przekroju poprzecznym warstwy smarnej, możemy obliczyć składowe prędkości w kierunku osi x i y:

0x01 graphic
(8.59)

gdzie 0x01 graphic
są składowymi wektorów prędkości, z jakimi poruszają się powierzchnie (rys. 8.6).

Scałkujemy teraz uproszczone równanie ciągłości (3.16) w granicach od 0x01 graphic
0x01 graphic
do

0x01 graphic
(8.60)

227



Wyszukiwarka