Rozdz3b, Zbigniew Kosma Podstawy Mechaniki Płynów


Strumień rotacji przez dowolny przekrój poprzeczny s rurki wirowej jest zawsze równy cyrkulacji wzdłuż linii , jest więc stały; stanowi to treść drugiego twierdzenia Helmholtza . Również cyrkulacja wektora 0x01 graphic
wzdłuż dowolnej linii zamkniętej i leżącej na rurce wirowej ma stałą wartość w każdej chwili.

Biorąc pod uwagę otrzymane własności przepływu wirowości przez rurkę wirową wprowadzono pojęcie włókna wirowego , jako granicy, do której dąży rurka wirowa przy zmniejszaniu jej przekroju poprzecznego i zwiększaniu wirowości w jej wnętrzu w taki sposób, że cyrkulacja pozostaje bez zmiany.

Z drugiego twierdzenia Helmholtza wynika ważny wniosek, że rurki wirowe i włókna wirowe nie mogą zaczynać się lub kończyć w obszarze płynu. Mogą one składać się więc z zamkniętych linii wirowych albo też zaczynać się i kończyć na ograniczeniach obszaru (ścianki sztywne, powierzchnia swobodna). Wynik ten możemy potwierdzić również w inny sposób. Obliczając

0x01 graphic
(3.34)

analogicznie do równania (3.20), które mówi o zachowaniu masy w przepływie płynu nieściśliwego, otrzymujemy równanie zachowania wirowości.

ĆWICZENIA

Przykład 3.1. Stacjonarne pole prędkości jest określone składowymi:

0x01 graphic

Wyznaczyć równania linii prądu i torów elementów.

Ponieważ 0x01 graphic
zatem pole prędkości jest płaskie.

Na podstawie wzoru (3.6) piszemy równanie linii prądu

0x01 graphic

z którego po scałkowaniu otrzymujemy

0x01 graphic

Z kolei równania torów elementów (3.7) przyjmują postać:

0x01 graphic

czyli

0x01 graphic

co jest wynikiem poprawnym, gdyż w przepływach stacjonarnych linie prądu i tory elementów pokrywają się.

Przykład 3.2. Wyznaczyć kształt linii prądu i rodzinę torów elementów płynu dla nieustalonego przepływu płaskiego określonego przez następujące składowe prędkości:

.

Całkujemy równanie różniczkowe rodziny linii prądu (3.6)

0x01 graphic

traktując t jako stały parametr. W wyniku otrzymujemy

lub

0x01 graphic

Linie prądu stanowią więc w każdej chwili rodzinę hiperbol i np. w chwili t = 0 przez punkt: przechodzi linia o równaniu 0x01 graphic

Dla określenia torów elementów należy scałkować równania (3.8):

0x01 graphic

a zatem:

0x01 graphic

0x01 graphic

Przykładowo, dla toru, po którym porusza się element płynu znajdujący się w chwili t = 0 w punkcie: jest i następnie mamy

0x01 graphic

Przykład 3.3. Dla nieustalonego przepływu płaskiego określonego polem prędkości:

0x01 graphic

należy wyznaczyć równania: linii prądu, toru elementu płynu i linii wysnutej.

Układ równań różniczkowych linii prądu przyjmuje postać:

0x01 graphic

Po scałkowaniu otrzymujemy równania parametryczne tych linii:

0x01 graphic

gdzie i są stałymi całkowania. W szczególności zakładając, że interesujemy się kształtem linii przechodzącej przez punkt o współrzędnych otrzymujemy , a więc:

0x01 graphic

0x01 graphic

Rys. 3.7

Ostatnie wzory uwidaczniają fakt, że kształt linii przechodzącej przez punkt o współrzędnych zmienia się z czasem. Dla chwili jest ona określona następującymi równaniami parametrycznymi:

Eliminując parametr s otrzymujemy równanie linii prądu, która przechodzi przez punkt o współrzędnych w chwili

0x01 graphic

Linia ta jest przedstawiona na rys. 3.7 łącznie z innymi liniami.

Układ równań toru elementu płynu w rozważanym przypadku przyjmuje postać:

0x01 graphic

Całkując to równanie mamy:

0x01 graphic

Są to równania parametryczne rodziny torów elementów płynu. W szczególności jeżeli interesujemy się torem elementu płynu, który w chwili przechodzi przez punkt o współrzędnych (1, 1), otrzymujemy:

0x01 graphic

Eliminując czas t z tych równań dostajemy równanie toru

0x01 graphic

Tor ten został przedstawiony na rys. 3.7, gdzie widać, że nie pokrywa się on z linią prądu, a jedynie jest do niej styczny w punkcie o współrzędnych (1, 1).

Jak wspomniano w rozdziale 3.2, linia wysnuta jest rozwiązaniem toru elementu płynu. Tak więc rodzina linii wysnutych jest określona równaniami:

0x01 graphic

Przykładowo dla warunków początkowych

, gdy

równania te są następujące:

Są to równania parametryczne linii wysnutej, która przechodzi przez punkt o współrzędnych . Widzimy, że w rozważanym przypadku kształt tej linii zależy od czasu (przepływ jest nieustalony). W szczególności dla mamy:

Eliminując parametr τ otrzymujemy równanie tej linii

,

co przedstawia rys. 3.7.

Należy zwrócić uwagę na fakt, że wszystkie trzy linie nie pokrywają się, mimo iż zostały wyznaczone dla tego samego przepływu. Przyczyną tego jest nieustaloność przepływu.

Przykład 3.4. Powierzchnię ziemi w otoczeniu miasta stanowi płaszczyzna , na której wyidealizowanym źródłem emisji spalin jest komin, znajdujący się w początku układu osi współrzędnych. Pole prędkości wiatru opisują składowe:

0x01 graphic

przy czym V oraz ω mają wartość stałą.

a. Wyznaczyć i naszkicować rodzinę linii prądu wiatru.

b. Określić rodzinę torów poruszania się cząsteczek spalin, które w czasie znajdują się w punkcie o współrzędnych 0x01 graphic
Przedstawić graficznie tory poruszania się kilku cząsteczek spalin, opuszczających komin w różnym czasie 0x01 graphic
Napisać równanie toru poruszania się elementu dymu, który wypływa z komina w chwili 0x01 graphic

a. W równaniu różniczkowym linii prądu wiatru:

0x01 graphic

0x01 graphic

zatem

0x01 graphic

Ogólnym rozwiązaniem otrzymanej zależności jest wyrażenie

0x01 graphic

Stała całkowania jest parametrem rodziny linii prądu, które są prostymi równoległymi, przecinającymi oś x w czasie t pod kątem 0x01 graphic
(rys. 3.8a). Ponieważ dla 0x01 graphic
jest również 0x01 graphic
więc stała

;

stąd linię prądu przechodzącą przez początek układu współrzędnych opisuje równanie

0x01 graphic

a) b)

0x01 graphic

Rys. 3.8

b. Równania różniczkowe torów poruszania się cząsteczek spalin możemy zapisać następująco:

0x01 graphic

0x01 graphic

a po scałkowaniu:

0x01 graphic
,

0x01 graphic
,

Dla czasu oraz współrzędnych 0x01 graphic

0x01 graphic
,

0x01 graphic
,

zatem:

0x01 graphic

0x01 graphic

Po przekształceniu równań i obustronnym podniesieniu do kwadratu otrzymamy:

0x01 graphic

0x01 graphic

skąd po dodaniu stronami jest

0x01 graphic
.

Wprowadzając podstawienia:

0x01 graphic

0x01 graphic
,

otrzymamy

0x01 graphic
.

Z przeprowadzonych obliczeń wynika, że linie torów poruszania się cząsteczek spalin są okręgami o promieniu , których środki znajdują się w punktach o współrzędnych (rys. 3.8b).

Dla toru poruszania się cząsteczki, która w chwili opuszcza komin:

, .

Ponieważ komin znajduje się w początku układu osi współrzędnych, przeto:

, ,

czyli:

wobec tego równanie toru poruszania się elementu dymu będzie miało następującą postać

Przykład 3.5. Wykazać, że dla dowolnej cząstki płynu poruszającej się po powierzchni kuli równanie ciągłości ma postać

0x01 graphic

gdzie to prędkości kątowe cząstki płynu w kierunku określonym kątem θ i kątem ϕ (szerokość i długość kątowa).

Równanie ciągłości we współrzędnych sferycznych jest określone wzorem

0x01 graphic

Ponieważ cząstka płynu porusza się po powierzchni kuli, to składowe prędkości są równe:

0x01 graphic

gdzie R jest promieniem kuli.

Podstawiając te wielkości do równania ciągłości oraz przyjmując mamy

0x01 graphic

Przykład 3.6. Ciecz porusza się w ten sposób, że każda cząstka opisuje okrąg na płaszczyźnie prostopadłej do pewnej osi ze środkiem na niej. Wykazać, że równanie ciągłości upraszcza się do postaci

0x01 graphic

gdzie ω jest prędkością kątową ruchu cząstki, której położenie określa się współrzędnymi: r, ϕ, z.

Zapisujemy równanie ciągłości (3.16) we współrzędnych cylindrycznych

0x01 graphic

Ponieważ cząstki poruszają się stale po tym samym okręgu, prędkości wyrażą się zależnościami:

0x01 graphic

Po uwzględnieniu tych szczególnych wartości składowych prędkości w równaniu ciągłości oraz przyjęciu r = R uzyskujemy

0x01 graphic

Przykład 3.7. Wyznaczyć pole przyspieszeń odpowiadające polu prędkości ok-reślonemu składowymi:

0x01 graphic

Składowe przyspieszenia obliczamy ze wzorów (3.11) i (3.12):

0x01 graphic
0x01 graphic

0x01 graphic

przy czym pierwsze składniki we wzorach na oraz reprezentują składowe lokalnej pochodnej prędkości; następne - składowe konwekcyjne. Łatwo to można sprawdzić obliczając:

0x01 graphic

0x01 graphic

0x01 graphic

Przykład 3.8. Zadane jest pole prędkości płynu:

0x01 graphic
0x01 graphic
0x01 graphic
( a = const) .

Obliczyć:

a) cyrkulację prędkości wzdłuż konturu trójkąta powstałego z przecięcia płaszczyzn układu współrzędnych z płaszczyzną 0x01 graphic

b) składowe wektora wirowości oraz prędkości deformacji objętościowych i po-staciowych.

Biorąc pod uwagę symetrię konturu trójkąta i postać funkcji opisujących pole prędkości (rys. 3.9) otrzymujemy

0x01 graphic

0x01 graphic

0x01 graphic

Rys. 3.9

Składowe wektora wirowości:

0x01 graphic

Prędkości deformacji objętościowych:

0x01 graphic

Prędkości deformacji postaciowych:

0x01 graphic

Przykład 3.9. Prędkości cząstek cieczy są proporcjonalne do ich odległości od osi O x i równoległe do niej tak, że:

0x01 graphic

Określić pole wirów.

Obliczając składowe wektora wiru 0x01 graphic
mamy:

0x01 graphic

0x01 graphic

0x01 graphic
.

Wartość bezwzględna wektora wiru jest równa:

0x01 graphic

a równanie linii wiru ma postać

0x01 graphic

Ponieważ , więc i linie wirowe są liniami płaskimi. Równanie po-wyższe przekształca się wiec następująco

0x01 graphic

stąd jest

0x01 graphic

i następnie po scałkowaniu

0x01 graphic

jest to równanie rodziny okręgów.

70



Wyszukiwarka