ROZDZ3A, Zbigniew Kosma Podstawy Mechaniki Płynów


3. Elementy kinematyki pŁynÓw

3.1. Opis ruchu płynu

Przedmiotem kinematyki płynów jest ustalenie ogólnych praw ruchu płynu względem danego układu odniesienia. Zakładając, że ruch płynu jest każdorazowo dany, zajmiemy się analizą kinematycznych wielkości charakteryzujących zachowanie się dowolnie wybranej cząstki płynu.

Dysponujemy dwiema podstawowymi metodami pozwalającymi opisać ruch płynu: metodą Lagrange'a i metodą Eulera.

Metoda Lagrange'a opiera się na badaniu ruchu wybranej cząstki płynu. Jeżeli przez 0x01 graphic
oznaczymy współrzędne dowolnego elementu płynu w chwili początkowej 0x01 graphic
, a przez 0x01 graphic
jego współrzędne kartezjańskie, to położenie elementu w dowolnej chwili t będzie określone równaniami:

0x01 graphic
(3.1)

Składowe wektora prędkości 0x01 graphic
elementu wyznaczają pierwsze pochodne cząstkowe względem czasu:

0x01 graphic
(3.2)

zaś składowe wektora przyspieszenia 0x01 graphic
drugie pochodne cząstkowe względem czasu:

0x01 graphic
(3.3)

0x01 graphic
(3.3cd.)

Stan kinematyczny płynu będzie znany, gdy będą dane równania (3.1) ÷ (3.3) dla wszystkich elementów płynu.

Metoda Eulera polega na badaniu ruchu płynu w określonych punktach przestrzeni i tylko pośrednio zajmuje się ruchem indywidualnego elementu płynu. Istotę tej metody stanowi badanie pól wielkości fizycznych (rozdz. 1.4) i w wyniku jej zastosowania otrzymuje się - najczęściej wykorzystywany - opis ruchu płynu w postaci funkcji (1.16) lub (1.18).

W mechanice płynów rozważamy często zmiany zachodzące w pewnej masie płynu. Wyodrębnioną masę płynu nazywa się obszarem płynnym , jeśli two-rzą ją wciąż te same elementy płynne lub też obszarem kontrolnym , jeśli tworzą ją wciąż te same punkty przestrzenne, natomiast znajdują się w niej, w miarę upływu czasu, coraz to inne elementy płynu.

Zamkniętą powierzchnię ograniczającą obszar płynny nazywamy powierz-chnią płynną; podobnie powierzchnię ograniczającą obszar kontrolny nazywamy powierzchnią kontrolną .

3.2. Graficzna reprezentacja pola prędkości

Niech będzie dane pole wektora prędkości płynu:

0x01 graphic
(3.4)

Linie tego pola wektorowego (styczne w każdym swym punkcie do wektora pola) nazywają się liniami prądu .

Oznaczając element linii prądu przez 0x01 graphic
, równanie linii prądu możemy napisać w postaci iloczynu wektorowego

0x01 graphic
(3.5)

lub też w postaci równoważnej

0x01 graphic
(3.6)

W równaniach (3.5) ÷ (3.6) czas t odgrywa rolę parametru; kształt linii prądu zależy więc od czasu i tylko w ustalonym polu prędkości będzie niezmienny.

Powierzchnia prądu jest to powierzchnia utworzona z linii prądu, przecinających dowolną linię, nie będącą linią prądu. Jeśli ta dowolna linia jest zamknięta, powierzchnia prądu bywa nazywana rurką prądu .

Linia, po której porusza się pojedynczy element płynu nazywa się torem ele-mentu ; jest ona graficzną reprezentacją matematycznego opisu ruchu Lagrange'a.

Równanie torów elementów płynu otrzymujemy wychodząc z definicji wektora prędkości. Oznaczając element toru przez 0x01 graphic
, możemy napisać równanie różniczkowe toru w następującej postaci

0x01 graphic
(3.7)

albo po przekształceniu

0x01 graphic
(3.8)

Między równaniami (3.6) i (3.8) zachodzi istotna różnica; w równaniu (3.8) czas nie jest stałym parametrem, ale zmienną.

Chwilowy obraz linii, łączącej wszystkie cząstki płynu, które wyszły z jednego źródła i przeszły przez ten sam punkt pola nazywa się linią wysnutą . Linię ta-ką wyznaczają elementy płynu odróżniającego się wizualnie od tła dostarczone do pola przepływu w określonym punkcie przestrzeni (np. dym w powietrzu lub atrament w wodzie, wypływające z dyszy do obszaru przepływu płynu); w ustalonym polu prędkości jest ona zarazem linią prądu i torem każdego elementu płynu, przechodzącego przez ten punkt.

3.3. Przyspieszenie elementu płynu

Przyspieszenie elementu płynu, traktowanego jako punkt, jest pochodną prędkości elementu względem czasu, wyraża się zatem wzorem

0x01 graphic
(3.9)

Zgodnie z regułami różniczkowania, różniczka zupełna 0x01 graphic
prędkości 0x01 graphic
0x01 graphic
jest określona następująco

0x01 graphic

i po wykorzystaniu wzorów (3.8) otrzymamy

0x01 graphic
(3.10)

Przy zastosowaniu operatora Hamiltona

0x01 graphic

wzór (3.10) możemy przepisać w postaci

0x01 graphic
(3.11)

gdzie

0x01 graphic
(3.12)

Równość (3.11) orzeka, że przyspieszenie elementu płynu, wyrażone pochod-ną substancjalną prędkości 0x01 graphic
jest sumą pochodnej lokalnej 0x01 graphic
oraz pochodnej konwekcyjnej 0x01 graphic
Nazwy tych pochodnych wynikają z ich sensu fizycznego. Mianowicie pochodna 0x01 graphic
określa zmiany zachodzące w poruszającym się, ale wciąż tym samym elemencie płynu, tzn. zmiany dotyczące jego substancji; pochodna lokalna określa zmiany zachodzące z upływem czasu w stałym punkcie przestrzeni; pochodna konwekcyjna określa zmiany prędkości związane z samym tylko przesunięciem elementu płynu w inne położenie.

Pojęcia pochodnej substancjalnej, lokalnej i konwekcyjnej są pojęciami ogólnymi i mogą być odniesione do dowolnych funkcji (1.18) - opisujących całą klasę wielkości fizycznych. Wobec tego można wprowadzić pojęcie operatora różniczkowania materialnego

0x01 graphic
, (3.13)

nazywanego często operatorem Stokesa .

3.4. Różniczkowe równanie ciągłości

Równanie ciągłości jest podstawowym równaniem mechaniki płynów, wynikającym z zasady zachowania masy i wyrażającym ciągły charakter przepływu.

Załóżmy, że ruch płynu został określony za pomocą pola prędkości (3.4) i pola gęstości płynu

0x01 graphic
(3.14)

0x01 graphic

Rys. 3.1

W przestrzeni wypełnionej poruszającym się płynem wyodrębniamy obszar kon-trolny w kształcie prostopadłościanu o bokach 0x01 graphic
(rys. 3.1).

Masa płynu zawarta w obszarze kontrolnym może się zmieniać z upływem czasu wskutek dwu przyczyn:

- zmiany gęstości płynu,

- dopływu przez ściankę powierzchni kontrolnej.

Jeśli w czasie t gęstość wynosiła 0x01 graphic
to w czasie 0x01 graphic
gęstość będzie równa

0x01 graphic

zatem zmiana masy płynu wynikająca ze zmiany jego gęstości wyniesie

0x01 graphic
(3.15)

Masa płynu przepływającego przez jednostkową powierzchnię w jednostce czasu jest równa iloczynowi jego gęstości i prędkości. Obliczymy masy płynu: wpływającą i wypływającą w czasie 0x01 graphic
do obszaru kontrolnego przez powierzchnie o bokach 0x01 graphic
i 0x01 graphic
(rys. 3.1):

- masa wpływająca

0x01 graphic

- masa wypływająca

0x01 graphic

Różnica tych wielkości jest następująca

0x01 graphic

Jeśli uwzględnimy przepływy przez pozostałe powierzchnie obszaru kontrolnego i porównamy je ze zmianą masy reprezentowaną wzorem (3.15), otrzymamy równanie ciągłości przepływu

0x01 graphic
(3.16)

Przy wykorzystaniu operatora diwergencji równanie (3.16) można zapisać w po-staci

0x01 graphic
(3.17)

a stąd, po wykonaniu różniczkowania i wprowadzeniu pochodnej substancjalnej (3.13), dostajemy

0x01 graphic
(3.18)

Dla ruchu stacjonarnego oraz cieczy równanie ciągłości upraszcza się do następujących postaci:

0x01 graphic
(3.19)

0x01 graphic
(3.20)

0x01 graphic

Rys. 3.2

W niektórych działach stosowanej mechaniki płynów przyjmuje się jednowymiarowość i stacjonarność przepływu jako podstawowe założenia upraszczające i wtedy stosujemy szczególne formy równania ciągłości. W celu uzyskania tych form równania ciągłości rozważymy ruch płynu przez kanał ograniczony ściankami kontrolnymi 0x01 graphic
i 0x01 graphic
(rys. 3.2).

Zdefiniujemy pojęcie średniej prędkości normalnej w przekroju σ

0x01 graphic
(3.21)

Masa płynu zawarta wewnątrz kanału musi być stała, a więc masa płynu dopływającego musi być równa masie płynu odpływającego. Zakładając, że gęstość płynu nie zmienia się w dowolnym przekroju poprzecznym otrzymamy

0x01 graphic
,

tzn, że strumień masy przepływającej przez dowolny przekrój poprzeczny kanału, zwany wydatkiem masowym albo masowym natężeniem prze-pływu , jest stały

0x01 graphic
(3.22)

Dla cieczy zależność (3.22) upraszcza się do postaci

0x01 graphic
(3.23)

w której stała Q nazywana jest wydatkiem objętościowym albo obję-tościowym natężeniem przepływu .

3.5. Struktura pola prędkości płynu

Rozpatrzmy ruch lokalny płynu, rozumiany jako ruch punktów elementu płynu względem bieguna, dowolnie wybranego w rozpatrywanym elemencie.

Z całej masy poruszającego się płynu wycinamy myślowo element objętościowy w kształcie prostopadłościanu o krawędziach 0x01 graphic
(rys. 3.3). Położenie tego elementu odpowiada pewnej chwili t; z upływem czasu element przesuwa się względem układu odniesienia, przy czym może on również zmieniać swój kształt.

Rozważać będziemy ruch punktu M względem bieguna A. Wektor prędkości 0x01 graphic
w punkcie M będzie równy wektorowi prędkości 0x01 graphic
w punkcie A powiększonemu o jego różniczkę zupełną 0x01 graphic

0x01 graphic

Rys. 3.3

0x01 graphic

składowe prędkości w punkcie M będą zatem wyrażać się wzorami:

0x01 graphic

0x01 graphic

0x01 graphic

Do prawych stron powyższych równań dodamy i odejmiemy kolejno następujące wielkości:

0x01 graphic

0x01 graphic

0x01 graphic

i następnie, po przekształceniach, zapiszemy składowe prędkości punktu M za pomocą wzorów:

0x01 graphic
(3.24)

w których dla zwięzłości zapisu wprowadzono oznaczenia:

0x01 graphic
(3.25)

Postaramy się obecnie wykazać, że wielkości (3.25) odpowiadają określonym odkształceniom i obrotowi elementu płynu przedstawionego na rysunku 3.3.

Jeśli różne od zera będą tylko wielkości wtedy boki elementu ulegną wydłużeniu lub skróceniu, a element pozostanie po odkształceniu prostopadłościanem. Oznacza to, że składowe te odpowiadają odkształceniom objętościowym, a zmiana objętości w czasie 0x01 graphic
wyrazi się wzorem

0x01 graphic
(3.26)

Wprowadzając pojęcie prędkości właściwej odkształcenia objętościowego, zdefiniowanej jako stosunek przyrostu objętości (3.26) do objętości początkowej i czasu, w którym ten przyrost nastąpił, otrzymamy

0x01 graphic
(3.27)

Załóżmy obecnie, że znikają wszystkie wielkości (3.25) oprócz . Wynika stąd zależność 0x01 graphic
i ze wzoru (3.24) mamy:

0x01 graphic

0x01 graphic

Element prostopadłościenny odkształcony zgodnie z powyższymi wzorami jest pokazany na rysunku 3.4. Wskutek różnic prędkości prostokąt ABCD odkształci się, przyjmując po upływie czasu 0x01 graphic
kształt równoległoboku AB'C'D'. Bok AB obróci się dookoła osi x o kąt 0x01 graphic
a bok AD o kąt 0x01 graphic
w związku z tym kąt prosty DAB zmniejszy się o wielkość 0x01 graphic
Wykorzystując rys. 3.4 łatwo obliczamy:

0x01 graphic

i następnie otrzymujemy

0x01 graphic
(3.28)

wielkości odpowiadają więc odkształceniom postaciowym elementu płynu.

0x01 graphic

Rys. 3.4

Na koniec załóżmy, że w wyrażeniach określających składowe prędkości (3.24) znikają wszystkie wielkości występujące w nawiasach kwadratowych; element płynu nie doznaje więc ani odkształceń objętościowych, ani też odkształceń postaciowych - jest zatem ciałem sztywnym.

Niech 0x01 graphic
0x01 graphic
będzie wektorem prędkości kątowej obrotu takiego elementu płynu wokół chwilowej osi przechodzącej przez biegun A, a 0x01 graphic
0x01 graphic
wektorem określającym położenie każdego jego punktu względem bieguna. Składowe prędkości każdego punktu elementu, wyznaczone ze wzoru 0x01 graphic
0x01 graphic
wyrażają się więc wzorami:

0x01 graphic

Obliczając następnie wektor 0x01 graphic
otrzymujemy zależność

0x01 graphic
0x01 graphic
(3.29)

z której wynikają oznaczenia różnic pochodnych zastosowane we wzorach (3.25).

Reasumując, stwierdzamy, że wektor prędkości przemieszczenia dowolnego punktu płynu (3.24) składa się w ogólnym przypadku z prędkości postępowej bieguna 0x01 graphic
, prędkości od obrotu elementu jako ciała sztywnego wokół osi przechodzącej przez biegun (wyrazy w nawiasach okrągłych) oraz prędkości deformacji objętościowych i postaciowych elementu (wyrazy w nawiasach kwadratowych); jest to pierwsze twierdzenie Helmholtza .

3.6. Pojęcia i twierdzenia dotyczące p*l wirowych

Ruchem wirowym płynu nazywamy taki ruch, którego pole prędkości spełnia warunek

0x01 graphic
0x01 graphic
(3.30)

Zgodnie ze wzorami (3.24), elementy płynu mogą doznawać obrotów wokół osi własnych, a prędkość kątowa może tworzyć w ogólnym przypadku niestacjonarne pole wektorowe:

0x01 graphic
(3.31)

Wektorowa linia pola prędkości kątowej 0x01 graphic
nazywana jest linią wirową. Jest to linia styczna do wektorów prędkości kątowych elementów płynu, które znajdują się na niej w danej chwili (rys. 3.5).

0x01 graphic

Rys. 3.5

Równanie różniczkowe linii wirowej jest analogiczne do równania linii prądu (3.5) ÷ (3.6):

0x01 graphic
(3.32)

w którym wektor 0x01 graphic
jest elementem linii wirowej.

Analogicznie do pojęcia powierzchni prądu i rurki prądu wprowadza się pojęcie powierzchni wirowej i rurki wirowej, jako powierzchni składającej się z linii wirowych, przecinających dowolną linię zamkniętą, nie będącą linią wirową.

Obliczymy teraz strumień rotacji przez dowolny przekrój poprzeczny rurki wirowej (rys. 3.6). Linia leżąca na rurce wirowej stanowi krawędź powierzchni utworzonej z powierzchni bocznej i przekroju poprzecznego rurki wirowej σ, którego brzegiem jest linia l.

0x01 graphic

Rys. 3.6

Zgodnie z definicją linii wirowych strumień rotacji przez powierzchnię jest równy zeru; zatem z twierdzenia Stokesa otrzymujemy

0x01 graphic
(3.33)

57



Wyszukiwarka