ROZDZ8D, Zbigniew Kosma Podstawy Mechaniki Płynów


Ostatnie równanie dowodzi, że:

0x01 graphic
0x01 graphic
(8.75)

Podstawiamy wyrażenie na 0x01 graphic
z (8.75) do (8.73) i otrzymujemy równanie dla funkcji

0x01 graphic

Zatem

0x01 graphic
,

gdzie A i B są dowolnymi stałymi. Jeżeli dla ustalonego y, to i

0x01 graphic

Jeżeli dla ustalonego t, to oraz:

0x01 graphic
,

0x01 graphic

Stąd wynika, że dla ustalonej wartości y prędkość dąży do V, gdy , tzn. w miarę upływu czasu ścianka nadaje danej cząstce cieczy prędkość 0x01 graphic

Przykład 8.7. Należy wykazać występowanie paradoksu Stokesa, uniemożliwiającego konstruowanie przepływów pełzających w płaskich obszarach nieograniczonych, rozważając zagadnienie opływu cylindra o promieniu R jednorodnym strumieniem stacjonarnym cieczy lepkiej (rys. 8.12).

Do rozwiązania zagadnienia użyjemy równania Stokesa ruchu cieczy lepkiej w układzie współrzędnych cylindrycznych. Wprowadzając funkcję prądu, spełniającą równanie ciągłości, określoną związkami:

0x01 graphic
,

0x01 graphic

Rys. 8.12

możemy zapisać równanie Stokesa w postaci równania biharmonicznego

0x01 graphic
(8.76)

Do rozwiązania równania (8.76) zastosujemy metodę rozdzielenia zmiennych przyjmując

0x01 graphic

Po podstawieniu będziemy mieli:

0x01 graphic

0x01 graphic

Rozwiązanie ogólne tego równania różniczkowego

0x01 graphic

ma postać

0x01 graphic
.

Powrót do funkcji f daje równanie różniczkowe

0x01 graphic

i następnie otrzymujemy:

0x01 graphic

0x01 graphic

Wracając do funkcji prądu i wyrażeń dla składowych wektora prędkości uzyskujemy:

0x01 graphic

0x01 graphic

0x01 graphic

gdzie wartości stałych całkowania A, B, C są wyznaczane z warunków brzegowych.

W nieskończenie wielkiej odległości od cylindra funkcja prądu powinna być równa funkcji prądu dla opływu cylindra cieczą doskonałą, co wymaga spełnienia warunku

0x01 graphic

Z tego warunku wynika, że powinny znikać stałe A i B oraz należy przyjąć:
C = U. Jedyna pozostała stała D nie może jednak spełniać równocześnie dwóch warunków znikania na okręgu składowych i wektora prędkości 0x01 graphic

Przykład 8.8. Kula o promieniu a, umieszczona w nieograniczonej przestrzeni wypełnionej cieczą lepką , obraca się z prędkością kątową wokół osi z. Zbadać ruch cieczy wywołany obrotem kuli, jeżeli obrót kuli jest powolny (ω małe).

Jako prędkość charakterystyczną ruchu możemy przyjąć prędkość punktów równika kuli równą 0x01 graphic
wtedy liczba Reynoldsa

0x01 graphic
.

Ponieważ prędkość kątowa jest mała, więc liczba Reynoldsa też będzie mała. Na mocy tego założenia w równaniach ruchu, zapisanych w układzie współrzędnych sferycznych 0x01 graphic
możemy odrzucić ich lewe strony. Tak otrzymane równania ruchu będą spełnione, jeżeli przyjmiemy, że 0x01 graphic
0x01 graphic
prędkość będzie zależeć tylko od r i θ

i powinna spełniać równanie

0x01 graphic
(8.77)

Na powierzchni kuli cząstki cieczy powinny poruszać się z tą samą prędkością liniową 0x01 graphic
jaką mają punkty kuli, stąd mamy warunek graniczny

0x01 graphic
(8.78)

Na mocy warunku (8.78) będziemy poszukiwać rozwiązania równania (8.77) w po-staci

0x01 graphic
(8.79)

Po podstawieniu (8.79) do (8.77) uzyskujemy równanie różniczkowe Eulera

0x01 graphic

którego rozwiązaniem ogólnym jest funkcja

0x01 graphic

Stałe całkowania 0x01 graphic
, wyznaczamy z warunków brzegowych. Otóż, dla powinno być V = 0, stąd 0x01 graphic
Z warunku (8.83) znajdziemy, że 0x01 graphic
więc

0x01 graphic
.

Obliczymy jeszcze wielkość momentu konieczną do podtrzymania ruchu kuli. Naprężenia styczne (sił tarcia) na powierzchni kuli:

0x01 graphic

i moment sił tarcia

0x01 graphic

0x01 graphic

Moment potrzebny do utrzymania ruchu jest więc równy: 0x01 graphic

Przykład 8.9. Dwie okrągłe płaskie równoległe płytki o promieniu R każda, znajdując się w niewielkiej odległości jedna nad drugą, zbliżają się jednostajnie do siebie. Określić ruch warstwy cieczy zawartej między płytkami oraz wyznaczyć siły oporu działające na każdą z płytek.

0x01 graphic

Rys. 8.13

Wprowadzając założenia upraszczające, analogiczne do założeń przyjętych w rozdziale 8.6, układ równań (8.38), zapisany (rys. 8.13a) we współrzędnych walcowych dla składowych wektora prędkości:

0x01 graphic
,

zredukuje się do postaci:

0x01 graphic

0x01 graphic

Ze względności ruchu wynikają warunki brzegowe:

u = w = 0 dla z = 0 ,

u = 0, w = − U dla z = h ,

dla r = R ;

tutaj U jest prędkością górnej płytki (dolną uważamy za nieruchomą), h jest grubością warstwy

Całkując równanie ruchu względem z (tj. w poprzek warstwy) i uwzględniając warunki brzegowe znajdziemy:

0x01 graphic
(8.80)

Po podstawieniu pierwszego z tych równań do drugiego

0x01 graphic
(8.81)

Zakładając z = 0 w górnej granicy całkowania, po wykonaniu działań otrzymamy

0x01 graphic
,

a po rozwiązaniu i uwzględnieniu warunku brzegowego znajdziemy

0x01 graphic
(8.82)

Siła oporu działająca na każdą z płytek równa się

0x01 graphic
.

Uwzględniając zależność (8.82) we wzorach (8.80) i (8.81) wyznaczamy składowe ektora prędkości:

0x01 graphic
.

*

Rozważane zagadnienie można również rozwiązać posługując się układem równań (8.38), zapisanych w układzie współrzędnych prostokątnych (rys. 8.13b):

0x01 graphic

dla składowych wektora prędkości:

0x01 graphic

Warunki brzegowe mają postać:

0x01 graphic
dla

0x01 graphic
dla

dla

Całkując równania ruchu względem z i uwzględniając sformułowane warunki brzegowe znajdziemy:

0x01 graphic

Po podstawieniu dwóch pierwszych równań do trzeciego mamy

0x01 graphic

Założenie z = 0 w górnej granicy całkowania prowadzi do równania Poissona dla ciśnienia

0x01 graphic
.

Poszukując rozwiązania tego równania w postaci funkcji

0x01 graphic

po uwzględnieniu warunku brzegowego: dla 0x01 graphic
otrzymamy

0x01 graphic
,

a następnie wyrażenia dla składowych prędkości:

0x01 graphic

0x01 graphic

Przykład 8.10. Dwie płaszczyzny, tworzące ze sobą bardzo mały kąt 0x01 graphic
przesuwają się względem siebie ze stałą prędkością U (rys. 8.14). Obszar między tymi płaszczyznami wypełnia ciecz lepka. Należy wyznaczyć ciśnienie i siły działające na górną płaszczyznę.

Dla stacjonarnego przepływu jednowymiarowego równanie Reynoldsa (8.61) upraszcza się do postaci

0x01 graphic

0x01 graphic

Rys. 8.14

skąd jest

0x01 graphic

i następnie

0x01 graphic

Po kolejnym całkowaniu dla warunków brzegowych

0x01 graphic

i dla grubości filmu olejowego określonego związkiem

0x01 graphic

otrzymujemy rozkład ciśnienia

0x01 graphic

oraz nieznany parametr 0x01 graphic

0x01 graphic

Znając ciśnienie możemy obliczyć siłę działającą na górną płytkę

0x01 graphic

gdzie współczynnik siły nośnej jest określony wzorem:

Różniczkując wyrażenie dla składowej prędkości (8.59a) obliczymy jeszcze na-prężenie styczne

0x01 graphic

i siłę tarcia działającą na płytkę

0x01 graphic

250



Wyszukiwarka

Podobne podstrony:
ROZDZ8C, Zbigniew Kosma Podstawy Mechaniki Płynów
ROZDZ8B, Zbigniew Kosma Podstawy Mechaniki Płynów
ROZDZ0, Zbigniew Kosma Podstawy Mechaniki Płynów
ROZDZ9C, Zbigniew Kosma Podstawy Mechaniki Płynów
ROZDZ11A, Zbigniew Kosma Podstawy Mechaniki Płynów
ROZDZ5A, Zbigniew Kosma Podstawy Mechaniki Płynów
ROZDZ7B, Zbigniew Kosma Podstawy Mechaniki Płynów
ROZDZ5B, Zbigniew Kosma Podstawy Mechaniki Płynów
ROZDZ5C, Zbigniew Kosma Podstawy Mechaniki Płynów
ROZDZ7A, Zbigniew Kosma Podstawy Mechaniki Płynów
ROZDZ10B, Zbigniew Kosma Podstawy Mechaniki Płynów
ROZDZ7D, Zbigniew Kosma Podstawy Mechaniki Płynów
ROZDZ4B, Zbigniew Kosma Podstawy Mechaniki Płynów
ROZDZ3A, Zbigniew Kosma Podstawy Mechaniki Płynów
ROZDZ12B, Zbigniew Kosma Podstawy Mechaniki Płynów
ROZDZ2B, Zbigniew Kosma Podstawy Mechaniki Płynów
ROZDZ11B, Zbigniew Kosma Podstawy Mechaniki Płynów
ROZDZ2A, Zbigniew Kosma Podstawy Mechaniki Płynów
ROZDZ7E, Zbigniew Kosma Podstawy Mechaniki Płynów

więcej podobnych podstron