background image

Z

ALICZENIE 

MECHANIKI PŁYNÓW

 

 

1.

 

Przedmiot mechaniki płynów. Pojęcia podstawowe 

2.

 

Klasyfikacja płynów 

3.

 

Pojęcie lepkości. Miary lepkości. Jednostki 

4.

 

Metody pomiaru ciśnienia. Wakuometry, barometry, manometry. Jednostki ciśnienia 

5.

 

Warunek równowagi. Powierzchnie stałego potencjału. Paradoks hydrostatyczny 

6.

 

Równowaga względna. Powierzchnie stałego potencjału. Przykłady równowagi względnej 

7.

 

Napór na ściany płaskie. Współrzędne środka naporu 

8.

 

Napór na ściany zakrzywione. Współrzędne środka naporu 

9.

 

Zjawisko wyporu. Odkrycie Archimedesa (czas, epoka i miejsce) 

10.

 

Stany stateczności pływania. Metacentrum i odległość metacentryczna 

11.

 

Kinematyka płynów. Cele, zadania, parametry kinematyczne 

12.

 

Metoda Eulera i Lagrange'a w kinematyce płynów 

13.

 

Pochodna wędrowna. Operator Stokesa, przyśpieszenie elementu płynu 

14.

 

Pojęcie cyrkulacji. Interpretacja fizyczna i analogia 

15.

 

Równania toru i linii prądu, rurka i powierzchnia prądu, struga 

16.

 

Przepływy potencjalne. Funkcja prądu i potencjału, jej interpretacja. Potencjał zespolony 

17.

 

Przepływy elementarne. Superpozycja przepływów 

18.

 

Zasada zachowanie masy. Równanie ciągłości 

19.

 

Zasada zachowania pędu. Równanie Eulera i Naviera – Stokesa 

20.

 

Równanie Daniela Bemoulliego dla płynu doskonałego, rzeczywistego i gazów 

21.

 

Formalizm matematyczny, zastosowania 

22.

 

Przygotować wykaz zagadnień opracowanych samodzielnie (poza wykładem i ćwiczeniami) 

 

 

 

 

 

 

background image

1.  Przedmiot mechaniki płynów. Poj

ę

cia podstawowe 

 

Zagadnienia  mechaniki  płynów  wyst

ę

puj

ą

  niemal  we  wszystkich  dziedzinach  techniki, 

mi

ę

dzy innymi: 

-  lotnictwie; 

- maszynach przepływowych (turbiny, spr

ęż

arki, pompy); - przepływowych układach regulacji 

automatycznej; 

- układach hydrauliki siłowej, - itd.. 

  Mianem  płynów  okre

ś

lamy  te  ciała  wyst

ę

puj

ą

ce  w  przyrodzie,  które  w  odró

ż

nieniu  od 

ciał  stałych  zmniejszaj

ą

  swe  kształty  w  sposób  trwały  pod  działaniem  dowolnie  małych  sił 

mechanicznych, je

ś

li działanie tych sił trwa dostatecznie długo. 

Płyny dzielimy na ciecze i gazy. 

Cieczami  nazywamy  te  spo

ś

ród  płynów,  które  zmieniaj

ą

  bardzo  nieznacznie  sw

ą

 

obj

ę

to

ść

 pod działaniem nawet bardzo wielkich sił mechanicznych. Inaczej mówimy, 

ż

e ciecze 

s

ą

 nie

ś

ci

ś

liwe. 

Gazami  nazywamy  te  płyny,  które  pod  działaniem  sił  mechanicznych  zmieniaj

ą

  swe 

obj

ę

to

ś

ci bardzo znacznie. Mówimy równie

ż

ż

e gazy s

ą

 to płyny 

ś

ci

ś

liwe. 

Przedmiotem  mechaniki  płynów  jest  badanie  zjawisk  wyst

ę

puj

ą

cych  podczas  ruchu, 

spoczynku  płynów,  ze  zwróceniem  uwagi  na  oddziaływanie  płynów  na 

ś

cianki  ciał  stałych 

ograniczaj

ą

cych płyn i 

ś

cianki ciał zanurzonych całkowiciew płynie. 

Zadaniem  mechaniki  płynów  jako  dyscypliny  podstawowej  jest  poznanie  praw 

rz

ą

dz

ą

cych  spoczynkiem  i  ruchem  płynów  oraz  stworzenie  podstaw  teoretycznych  

i do

ś

wiadczalnych dla całego szeregu dziedzin specjalistycznych. 

W zale

ż

no

ś

ci od rodzaju u

ż

ytych kryteriów mechanik

ę

 płynów mo

ż

na podzieli

ć

 na:  

1. hydromechanik

ę

, to znaczy mechanik

ę

 cieczy; 

 

- mechanik

ę

 cieczy idealnej; 

 

- mechanik

ę

 cieczy rzeczywistej; 

2.  aeromechanik

ę

, to znaczy mechanik

ę

 gazu. Ponadto mechanik

ę

 płynów dzieli si

ę

 na: 

 

  - statyk

ę

- kinematyk

ę

;  

- dynamik

ę

Statyka  płynów  jest  nauk

ą

  o  zjawiskach  zachodz

ą

cych  podczas  spoczynku  cieczy  

i gazów. 

background image

Kinematyka  płynów  zajmuje  si

ę

  analiz

ą

  ruchu  w  oderwaniu  od  wyst

ę

puj

ą

cych  sił. 

Dynamika płynów, czyli dynamika cieczy (hydrodynamika) i dynamika gazów (aerodynamika), 

jest nauk

ą

 o ruchu cieczy i gazów w powi

ą

zaniu z wyst

ę

puj

ą

cymi siłami. 

Istnieje  równie

ż

  inny  podział  mechaniki  płynów,  a  mianowicie  w  oparciu  o  własno

ś

ci 

dyssypatywne  płynów.  Bior

ą

c  za  podstaw

ę

  wła

ś

ciwo

ś

ci  takie  jak  lepko

ść

  i  przewodno

ść

 

ciepln

ą

  mo

ż

na  podzieli

ć

  mechanik

ę

  płynów  na  tak  zwan

ą

  mechanik

ę

  płynów  doskonałych 

opart

ą

  o  model  płynu  nielekkiego  i  nie  przewodz

ą

cego  ciepła  i  na  mechanik

ę

  płynów 

rzeczywistych. 

Z  podziałem  mechaniki  płynów  wi

ąż

e  si

ę

  poj

ę

cie  hydrauliki.  Jest  ona  najstarszym 

historycznie  działem  mechaniki  cieczy  i  opiera  si

ę

  na  modelu  jednowymiarowego  przepływu 

cieczy  lepkiej  i  nie

ś

ci

ś

liwej.  Model  ten  polega  na  stosowaniu  empirycznych  współczynników 

ujmuj

ą

cych  straty  tarcia  w  przepływie  cieczy.  Przez  mechanik

ę

  płynów  rozumiemy  nauk

ę

zajmuj

ą

c

ą

 si

ę

 badaniem ruchów płynów ( cieczy i gazów b

ę

d

ą

cych pod działaniem sił ). 

W  obliczeniach  stosuje  si

ę

  przybli

ż

one  modele  dotycz

ą

ce  idealnych,  cieczy 

rzeczywistych, gazów pół-doskonałych i rzeczywistych. 

Prekursorami  w  dziedzinie  mechaniki  płynów  byli  Newton,  który  sformułował  pojecie 

lepko

ś

ci,  Euler  jako  pierwszy  napisał  ró

ż

niczkowe  równania  cieczy  i  Archimedes  staro

ż

ytny 

konstruktor,  który  zbudował  maszyn

ę

,  na  jego  cz

ęść

  nazwan

ą

  "  Spiral

ą

  Archimedesa  "  , 

słu

żą

c

ą

 do wymuszania przepływu cieczy. 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

background image

2. Klasyfikacja płynów 

 

Płyn nielepki i nie

ś

ci

ś

liwy zwany idealnym. 

 

Jego definicja wynika wprost z definicji cieczy idealnej 

Mówimy, 

ż

e  płyn  jest  nie

ś

ci

ś

liwy,  je

ś

li  jego  współczynnik 

ś

ci

ś

liwo

ś

ci  k  równa  si

ę

  zero. 

Współczynnik 

ś

ci

ś

liwo

ś

ci oznacza si

ę

 wzorem. 

dp

d

k

σ

=

 

Wyra

ż

a on zmian

ę

 g

ę

sto

ś

ci 

σσσσ

 pod wpływem zmiany ci

ś

nienia p 

Jak  wiemy  ciecze  s

ą

  nieznacznie 

ś

ci

ś

liwe,  to  znaczy,  ich  współczynnik 

ś

ci

ś

liwo

ś

ci  jest 

bardzo  mały,  podobnie  zachowuj

ą

  si

ę

  równie

ż

  gazy  przy  małych  pr

ę

dko

ś

ciach.  Opisuj

ą

przepływ  cieczy,  oraz  przepływ  gazu  przy  małych  pr

ę

dko

ś

ciach,  mo

ż

emy  przyj

ąć

  stał

ą

 

warto

ść

 g

ę

sto

ś

ci. 

Ten typ płynu opisuje równanie Eulera 

gradp

p

F

dt

dv

1

=

 

Definicja lepko

ś

ci znajduje si

ę

 w dalszej cz

ęś

ci rozprawki. 

 

Płyn lepki i nie

ś

ci

ś

liwy 

 

Modelem płynu lepkiego i nie

ś

ci

ś

liwego badamy przepływy w warstwie przy

ś

ciennej.  

dn

dv

 

Warstwa  przy

ś

cienna  charakteryzuje  si

ę

  tym, 

ż

e  gradienty  pr

ę

dko

ś

ci  przepływu,  s

ą

  

w  niej  bardzo  du

ż

e,  oddalaj

ą

c  si

ę

  od 

ś

cianek  w  gł

ą

b  strumienia  pr

ę

dko

ś

ci  te  gwałtownie 

rosn

ą

. Poza warstw

ą

 przy

ś

cienn

ą

 gradienty pr

ę

dko

ś

ci s

ą

 bardzo małe.  

 

υ

υ

2

+

=

gradp

F

dt

dv

 

 

2

2

2

2

2

2

2

dz

v

d

dy

v

d

dx

v

d

+

+

=

υ

 

W  obszarze  warstwy  przy

ś

ciennej  napr

ęż

enia  styczne,  uzyskuj

ą

  znaczne  warto

ś

ci, 

niezale

ż

nie od lepko

ś

ci płynu jest niewielka, ze wzgl

ę

dem na gradienty pr

ę

dko

ś

ci napr

ęż

enia 

styczne s

ą

 pomijalnie małe. Ten typ płynu opisuje równanie: 

 

background image

 

Płyn nielepki i 

ś

ci

ś

liwy 

 

Modelem  tego  płynu  posługujemy  si

ę

  w  dynamice  gazów.  Jest  to  nauka  zajmuj

ą

ca  si

ę

 

przepływami z du

ż

ymi pr

ę

dko

ś

ciami. 

Poniewa

ż

 zachodz

ą

 wyra

ź

ne zmiany g

ę

sto

ś

ci s nale

ż

y uwzgl

ę

dni

ć

 zale

ż

no

ść

 p od ci

ś

nienia  

i temperatury. Tymi zale

ż

no

ś

ciami zajmuje si

ę

 termodynamika. Powi

ą

zania z termodynamik

ą

 

mog

ą

  by

ć

  daleko  bardziej  id

ą

ce  np.:  w  przepływach  w  których  doprowadzamy  lub 

odprowadzamy  ciepło  lub  prac

ę

.  Gazy  s

ą

  płynami  ekspansywnymi,  to  znaczy  zajmuj

ą

  cał

ą

 

woln

ą

 przestrze

ń

, w której s

ą

 zamkni

ę

te. 

 

Płyn lepki i 

ś

ci

ś

liwy 

 

Charakteryzuje si

ę

 lepko

ś

ci

ą

 i niewielk

ą

 

ś

ci

ś

liwo

ś

ci

ą

, jest modelem płynów rzeczywistych 

wyst

ę

puj

ą

cych  w  przyrodzie.  Na  tym  modelu  bazuj

ą

  najbardziej  ogólne  i 

ś

cisłe  rozwi

ą

zania. 

Jednak

ż

e  przy  rozwi

ą

zywaniu  równa

ń

  opisuj

ą

cych  ten  typ  płynu  napotykamy  na  ogromne 

trudno

ś

ci,  głównie  natury  matematycznej,  ze  wzgl

ę

du  na  du

ż

y  stopie

ń

  skomplikowania  tych 

równa

ń

Obserwuj

ą

c  zjawisko  przepływu  płynów  zauwa

ż

amy, 

ż

e  jest  ono  ogromnie 

skomplikowane.  Opisanie  tak  zło

ż

onych  zjawisk  nie  jest  łatwe.  Dlatego  te

ż

  ze  wzgl

ę

du  na 

trudno

ś

ci  w  ich  rozwi

ą

zywaniu  staramy  si

ę

  je  nieco  upro

ś

ci

ć

,  pomijaj

ą

c  te  wielko

ś

ci  które  

w stosunku do pozostałych s

ą

 bardzo małe. 

 

 

 

 

 

 

 

 

 

 

 

 

 

background image

3. Poj

ę

cie lepko

ś

ci. Miary lepko

ś

ci. Jednostki 

 

Hipoteza Newtona 

 

Płyny  rzeczywiste  wykazuj

ą

  zdolno

ść

  przenoszenia  napr

ęż

e

ń

  stycznych,  przy  czym 

napr

ęż

enia 

powstaj

ą

 

mi

ę

dzy 

s

ą

siednimi 

warstwami 

płynu 

poruszaj

ą

cymi 

si

ę

  

z  ró

ż

nymi  pr

ę

dko

ś

ciami.  Napr

ęż

enia  styczne  powstaj

ą

  równie

ż

  pomi

ę

dzy  poruszaj

ą

cym  si

ę

 

płynem  i  ciałem  stałym,  nie  wyst

ę

puj

ą

  natomiast  w  czasie  spoczynku,  lub  w  płynie 

poruszaj

ą

cym  si

ę

  z  wyrównan

ą

  pr

ę

dko

ś

ci

ą

  przepływu  w  całym  przekroju.  Zgodnie  

z  hipotez

ą

  Newtona,  napr

ęż

enia  styczne  i  wyst

ę

puj

ą

ce  mi

ę

dzy  s

ą

siednimi  warstwami,  lub 

elementami  płynu  s

ą

  proporcjonalne  do  przyrostu  pr

ę

dko

ś

ci  w  kierunku  normalnym  do 

kierunku przepływu. Sytuacje t

ą

 przedstawia rysunek numer 1, oraz opisuje poni

ż

szy wzór. 

y

v

=

η

τ

 

gdzie:  

τ

 - napr

ęż

enia styczne; 

η

- współczynnik lepko

ś

ci dynamicznej, zwany lepko

ś

ci

ą

 dynamiczn

ą

 

 

Rys.1. Ilustracja do równania opisuj

ą

cego hipotez

ę

 Newtona 

 

 

 

 

background image

Pojecie lepko

ś

ci dynamicznej 

 

Lepko

ść

  dynamiczna 

η

  jest  funkcj

ą

  temperatury,  ci

ś

nienia  i  rodzaju  płynu.  Zale

ż

no

ść

 

lepko

ś

ci  płynu  od  ci

ś

nienia  jest  nieznaczna,  i  ro

ś

nie  bardzo  wolno  wraz  z  jego  wzrostem,  

a jedynym wyj

ą

tkiem jest woda, która w zakresie poni

ż

ej 32°C maleje ze wzrostem ci

ś

nienia. 

Zale

ż

no

ść

 lepko

ś

ci 

η

 od temperatury jest natomiast bardzo znaczna i zupełnie odmienna dla 

cieczy  i  gazów.  Zale

ż

no

ść

  t

ą

  ilustruje  rysunek  2  .  lepko

ść

  dynamiczna 

η

  cieczy  maleje  ze 

wzrostem  temperatury.  Wynika  z  tego, 

ż

e  w  cieczach  ruch  molekuł  jest  stosunkowo  mało 

intensywny,  wobec  czego  napr

ęż

enia  styczne  powstaj

ą

  głównie  w  skutek  molekularnych  sil 

spójno

ś

ci; przy wzro

ś

cie temperatury rosn

ą

 odległo

ś

ci pomi

ę

dzy molekułami, a zatem malej

ą

 

siły spójno

ś

ci. W gazach, jak mo

ż

emy wnioskowa

ć

 z rysunku jest na odwrót. 

 

 

Rys.2. Zale

ż

no

ść

 lepko

ś

ci od temperatury 

Jednostka lepko

ś

ci dynamicznej 

η

 w układzie SI wynika z równania Newtona: 





s

m

kg

*

1

 

Pojecie lepko

ś

ci kinematycznej 

 

Współczynnik lepko

ś

ci kinematycznej, zwany lepko

ś

ci

ą

 kinematyczn

ą

, okre

ś

lony jest 

wzorem:  

ρ

η

υ

=

 

 

Jednostk

ą

 lepko

ś

ci kinematycznej jest: 

background image

1stokes [1St]=

s

m

2

4

10

 

W u

ż

yciu s

ą

 stosowane jednostki mniejsze, takie jak: 

1centistokes [1cSt] 

Poza układem SI dla okre

ś

lenia lepko

ś

ci u

ż

ywa si

ę

 szeregu innych jednostek, takich 

-  stopie

ń

 Eulera [E]; 

-  sekunda Redwood'a [R sek];  

-  i inne. 

Jednostki te wynikaj

ą

 ze sposobu pomiaru i rodzaju u

ż

ytej aparatury. 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

background image

4. Metody pomiaru ci

ś

nienia. 

Manometry, barometry, wakuometry. Jednostki ci

ś

nienia

 

 

Przyrz

ą

dy do pomiaru ci

ś

nie

ń

 

 

Przyrz

ą

dy  do  pomiaru  ci

ś

nienia  mierz

ą

  nie  jego  warto

ść

  lecz  ró

ż

nice  tej  warto

ś

ci  

i warto

ś

ci ci

ś

nienia odniesienia. Ci

ś

nieniem odniesienia najcz

ęś

ciej jest: 

- pró

ż

nia absolutna (bezwzgl

ę

dna); 

- ci

ś

nienie panuj

ą

ce aktualnie w miejscu i czasie pomiaru. 

Ogólnie  przyrz

ą

dy  do  pomiaru  ci

ś

nie

ń

  nazywaj

ą

  si

ę

  ci

ś

nieniomierzami.  Przyrz

ą

dy  do 

pomiaru  ci

ś

nie

ń

  absolutnych  nazywaj

ą

  si

ę

  barometrami  i  ci

ś

nieniomierzami  ci

ś

nienia 

absolutnego, nad ci

ś

nie

ń

 - manometrami, a podci

ś

nie

ń

 wakuometrami. Poza tym wyró

ż

nia si

ę

 

jeszcze  ci

ś

nieniomierze  (manometry)  ró

ż

nicowe  -  do  pomiaru  ró

ż

nicy  ci

ś

nie

ń

.  Poni

ż

szy 

rysunek przedstawia wy

ż

ej wymienione zale

ż

no

ś

ci w sposób graficzny. 

 

Rys.3. Rodzaje ci

ś

nie

ń

 

Manometry  

Dzielimy je na: 

- hydrostatyczne ( cieczowe ) 

- pr

ęż

ne ( rurkowe lub przeponowe ) 

- przeponowe 

background image

 

Elementem  pr

ęż

nym  tego  manometru  jest  przepona  (0)  płaska  lub  falista.  Pod  wpływem 

ci

ś

nienia p przepona ulega odkształceniu. Odkształcenie to przenosi si

ę

 przez przekładni

ę

(S) 

na wskazówk

ę

 (W) 

Poni

ż

szy rysunek przedstawia zasad

ę

 działania manometru. 

 

 

Metoda  pomiaru manometrem 

Manometr (1) mierzy ci

ś

nienie statyczne ps płynu (wektor pr

ę

dko

ś

ci w płynie jest styczny do 

otworu wlotowego manometru). Je

ż

eli otwór wlotowy manometru jest prostopadły do wektora 

pr

ę

dko

ś

ci  w,  to  manometr  (2)  wska

ż

e  wzrost  ci

ś

nienia.  Wzrost  ten  jest  wywołany  zmian

ą

 

energii  kinetycznej  płynu  o  pr

ę

dko

ś

ci  w  i  g

ę

sto

ś

ci  p  na  energi

ę

  potencjaln

ą

,  jest  ci

ś

nieniem 

kinematycznym i zwyczajowo nazywamy je ci

ś

nieniem dynamicznym 

 

Barometr 

 

Barometr  -  przyrz

ą

d  do  pomiaru  ci

ś

nienia  atmosferycznego.  Pierwszy  tak  zwany  barometr 

rt

ę

ciowy  został  wynaleziony  w  roku  1643  przez  E.  Torricellego  w  zwi

ą

zku  z  jego  badaniami 

nad  ci

ś

nieniem  sferycznym.  Barometr  rt

ę

ciowy  został  ulepszony  w  1665  przez  R.  Hook'ea 

który wprowadził podziałk

ę

 umo

ż

liwiaj

ą

c

ą

 bezpo

ś

rednie odczytywanie wielko

ś

ci. 

 

Poni

ż

szy rysunek przedstawia zasad

ę

 działania barometru. 

background image

 

Barometr 

wykorzystuje 

si

ę

 

do 

pomiaru 

ci

ś

nie

ń

 

barometrycznych.  

Do  zbiorniku  (1)  wstawiona  jest  szklana  rurka  manometryczna  (2),  z  której  uprzednio 

wypompowano  cale  powietrze,  do  stanu  pró

ż

ni  absolutnej.  W  przestrzeni  pomi

ę

dzy  jej 

górnym  zamkni

ę

tym  ko

ń

cem,  a  meniskiem  rt

ę

ci  panuje  pró

ż

nia  Torrcelleogo  (3)  pró

ż

nia 

bezwzgl

ę

dna Poło

ż

enie słupka rt

ę

ci okre

ś

la ci

ś

nienie barometryczne (otoczenia) pb. 

Wynaleziony 

barometr 

odegrał 

podstawowe 

znaczenie 

dla 

rozwoju 

metrologii.  

B.  Pascal  powtarzaj

ą

c  i  kontynuuj

ą

c  badania  Torricellego,  zauwa

ż

ył, 

ż

e  ci

ś

nienie 

atmosferyczne  zale

ż

y  nie  tylko  od  wysoko

ś

ci  miejsca,  w  którym  przyrz

ą

d  si

ę

  znajduje,  ale 

tak

ż

e od stanu pogody. 

 

Wakuometr 

 

Wakuometr jest ci

ś

nieniomierzem słu

żą

cym do pomiaru podci

ś

nie

ń

. Zasada działania i 

budowa  nie  ró

ż

ni

ą

  si

ę

  od  manometrów  pr

ęż

nych  i  hydrostatycznych.  Cz

ę

sto  wykonuje  si

ę

 

manometry  pr

ęż

ne,  które  mog

ą

  słu

ż

y

ć

  do  pomiaru  nadci

ś

nie

ń

  i  podci

ś

nie

ń

.  Nazywamy  je 

mano-wakuometrami. 

 

 

 

 

 

 

 

 

background image

Jednostki ci

ś

nienia. 

 

W  układzie  jednostek  miar  SI  główn

ą

  jednostk

ą

  ci

ś

nienia  jest  Pascal  (  Pa  )  czyli  niuton  na 

metr kwadratowy 

2

2

*

1

1

1

s

m

kg

m

N

Pa

=

=

 

Pascal jest jednostk

ą

 mał

ą

, dlatego w praktyce stosuje si

ę

 megapascal 

1MPa =10

6

 Pa 

oraz bar 

1bar =10

5

 Pa 

Dawniej  manometry  skalowano  w  atmosferach  technicznych,  czyli  kilogramach  siły  na 

centymetr kwadratowy. 

Pa

cm

kG

at

4

2

10

*

80665

,

9

1

1

=

=

 

Cz

ę

sto mierzy si

ę

 ci

ś

nienie za pomoc

ą

 wysoko

ś

ci słupa cieczy 

Przy małych ci

ś

nieniach wzgl

ę

dnych cz

ę

sto mierzy si

ę

 ci

ś

nienie za pomoc

ą

 wysoko

ś

ci słupa 

wody:  

Pa

m

kG

O

mmH

80665

,

9

1

1

2

2

=

=

 

Przy  wi

ę

kszych  ci

ś

nieniach  wzgl

ę

dnych  cz

ę

sto  mierzy  si

ę

  ci

ś

nienie  za  pomoc

ą

  wysoko

ś

ci 

słupa rt

ę

ci ( Hg ) 

Tor  (  Tr  )  jest  jednostk

ą

  ci

ś

nienia  równ

ą

  ci

ś

nieniu  wywieranemu  w  pró

ż

ni  przez  słup  rt

ę

ci  

o wysoko

ś

ci 1 mm i temperaturze 0°C przy normalnym przyspies zeniu ziemskim. 

g

n

 = 9,8066





2

s

m

 

1 Tr = 1 mm Hg = 133,3224 Pa 

Jako normalne ci

ś

nienie fizyczne przyj

ę

to ci

ś

nienie jednej atmosfery fizycznej 

Pn = 1 atm = 760 Tr = 101325 Pa 

 

Wnioski 

Nauka  mechaniki  płynów  skłania  nas  do  rozwa

ż

a

ń

  nad  procesami  zachodz

ą

cymi  obok  nas,  

a  nawet  w  nas.  Układ  krwiono

ś

ny  jest  takim  specyficznym  układem  hydrodynamicznym. 

Układ  krwiono

ś

ny  jest  takim  specyficznym  układem  hydrodynamicznymi,  w  którym  zachodzi 

wiele przemian. Wiemy ze ci

ś

nienie jaki temperatura krwi nie jest stała a co si

ę

 z tym wi

ąż

background image

jej  g

ę

sto

ść

  jest  zmienna.  St

ą

d  te

ż

  wyznaczenie  relacji  krew  -  zastawki  jest  ogromnie 

skomplikowane, ale ciekawe. 

Wchodz

ą

c  w  XXI  wiek  coraz  bardziej  modne  staje  si

ę

  pozyskiwanie  energii  

w  sposób  ekologiczny.  Jednym  ze  sposobów  pozyskiwania  energii  ekologicznej  jest 

budowanie  elektrowni  wodnych,  zarówno  szczytowo  -  pompowych  jak  i  zaporowych. 

Dziedzina  hydroenergetyki  wi

ąż

e  si

ę

 

ś

ci

ś

le  z  mechanik

ą

  płynów.  Projektowanie  tych 

elektrowni  wymaga  bowiem  znajomo

ś

ci  praw  rz

ą

dz

ą

cych  cieczami,  jak  na  przykład:  prawa 

naporu  hydrostatycznego,  i  innych,  bez  których  wykorzystanie  energii  tkwi

ą

cej  

w cieczach byłoby niemo

ż

liwe. 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

background image

5 i 6. Warunek równowagi równowaga wzgl

ę

dna, powierzchnie stałego 

potencjału 

 

W  hydrostatyce  modele  cieczy  lepkiej  i  idealnej  s

ą

  równowa

ż

ne  sobie.  Rozpatruj

ą

równowag

ę

 ciała płynnego poddanego działaniu sił powierzchniowych mo

ż

emy okre

ś

li

ć

 jego 

stan  napi

ę

cia  przy  pomocy  tensora  napr

ęż

enia,  podstawiaj

ą

c  w  nim  zerowe  warto

ś

ci 

napr

ęż

e

ń

 stycznych. 

zx

yz

xy

τ

τ

τ

=

=

 

=

z

y

x

ij

δ

δ

δ

δ

0

0

0

0

0

0

 

i,j=x,y,z    

δ

=

δ

ij

 e

i

 e

j

 

W przypadku cieczy w której nie uwzgl

ę

dniamy sił masowych , tensor napr

ęż

e

ń

 tworzy 

przynale

ż

n

ą

 mu struktur

ę

 , poniewa

ż

 wszystkie składowe tensora s

ą

 sobie równe. 

δ

xx

=

δ

yy

=

δ

zz

=p 

Otrzymali

ś

my  w  ten  sposób  twierdzenie  które  w  hydrostatyce  nosi  nazw

ę

  prawa 

Pascala. Ustalmy element płynu P(x,y,z) za pomoc

ą

 prostopadło

ś

cianu o wymiarach dx,dy,dz 

ś

cianach  zorientowanych  równolegle  do  osi  układu  współrz

ę

dnych,  a  ponadto  funkcja 

ci

ś

nienia p(x,y,z) jest niewiadom

ą

 

Układamy warunki równowagi na kolejne osie układu współrz

ę

dnych 

0

)

(

:

=

+

dydz

dx

dx

dp

p

pdydz

XdV

X

ρ

 

0

)

(

:

=

+

dxdz

dy

dy

dp

p

pdxdz

YdV

Y

ρ

 

0

)

(

:

=

+

dxdy

dz

dz

dp

p

pdxdy

ZdV

X

ρ

 

Po podstawieniu otrzymujemy 

background image

0

=

dx

dp

X

ρ

 

0

=

dy

dp

Y

ρ

 

0

=

dZ

dp

Z

ρ

 

Otrzymali

ś

my w ten sposób równanie równowagi płynu pozostaj

ą

cego w spoczynku lub 

ogólnie warunki równowagi płynów 

Powy

ż

sze równania mo

ż

na zast

ą

pi

ć

 równaniem wektorowym 

+

+

)

(

k

j

i

Z

Y

X

ρ

 

i

dx

dp

 

j

dy

dp

 

k

dz

dp

= 0 

 

(1) 

0

=

gradp

F

ρ

 

 

Skalarne  pole  ci

ś

nie

ń

  mo

ż

na  wyznaczy

ć

  całkuj

ą

c  te  równania,  otrzymamy  wtedy 

szukan

ą

  funkcj

ę

  p(x,y,z).  Je

ś

li  funkcja  ta  jest  stała  to  otrzymujemy  powierzchni

ę

  stałego 

ci

ś

nienia nazywane powierzchniami izobarycznymi. Stwierdzamy ,

ż

e wektor siły obj

ę

to

ś

ciowej 

jest prostopadły w ka

ż

dym punkcie powierzchni izobarycznej prostopadły do niej. 

 

Na podstawie (1) mo

ż

na wyprowadzi

ć

 prawo Pascala : 

dx

dp

 

dy

dp

=

 

0

=

=

dz

dp

 

Je

ś

li (1) pomno

ż

y

ć

 przez przyrost wektora promienia , to otrzymamy  

r

gradpd

r

Fd

=

ρ

 

 

            (2) 

dp

r

Fd

=

ρ

 

Niech  F  b

ę

dzie  sił

ą

  potencjaln

ą

.  Wtedy  musi  istnie

ć

  takie  U(x,y,z)  , 

ż

e  F=-gradU  (3). 

Uwzgl

ę

dniaj

ą

c ( 1 ) i (2) otrzymujemy 

dp

gradUdr

=

ρ

 

0

=

+

dp

dU

ρ

 

Je

ś

li p=const , to mo

ż

emy to scałkowa

ć

 i otrzymamy  

const

p

U

=

+

ρ

 

background image

Jedynie  siły  masowe  s

ą

  w  stanie  wywoła

ć

  równowag

ę

  cieczy  nie

ś

ci

ś

liwej.  

W  szczególnym  przypadku  gdy  p=const  (  dp=0  )  ,  to  pdU=0  ~  U=const.  Czyli  powierzchnie 

stałego ci

ś

nienia s

ą

 równie

ż

 powierzchniami ekwipotencjalnymi. 

 

Wyra

ż

enie w nawiasie jest równe zupełnej dU, st

ą

dp=

ξ

*dU 

Z  równania  tego  wynika, 

ż

e  dla  dp=0  jest  dU=0.  Oznacza  to, 

ż

e  powierzchnie  jednakowego 

ci

ś

nienia  (dp=0,  p=const),  czyli  powierzchnie  izobaryczne,  s

ą

  w  polu  sił  masowych 

jednocze

ś

nie  powierzchniami  stałego  potencjału  (dU=0,  U=const),  czyli  powierzchniami 

ekwipotencjalnymi

Powierzchnie  ekwipotencjalne  z  natury  rzeczy  nie  mog

ą

  si

ę

  wzajemnie  przecina

ć

  

i  ka

ż

da  z  nich  jest  albo  powierzchni

ą

  zamkni

ę

t

ą

,  albo  ko

ń

czy  si

ę

  na 

ś

cianie  zbiornika 

(naczynia).  Swobodne  powierzchnie  cieczy,  czyli  powierzchnie  oddzielaj

ą

ce  cieczy  od  gazu, 

s

ą

 oczywi

ś

cie powierzchniami ekwipotencjalnymi. 

 

Dla  powierzchni  ekwipotencjalnych  (dU=0)  z  równania  wy

ż

ej  wynika, 

ż

e  siły  masowe  przy 

przesuni

ę

ciu wzdłu

ż

 tych powierzchni nie wykonuj

ą

 

ż

adnej pracy (X*dx=Y*dy=Z*dz=0). 

 

Parcie cieczy na powierzchnie 

ś

cian 

 
 

Parcie  jest  to  siła,  jak

ą

  wywiera  ciecz  w  spoczynku  na  dowolnie  zorientowan

ą

 

powierzchni

ę

. Rozwa

ż

ania obejm

ą

 parcie na powierzchnie płask

ą

 i zakrzywion

ą

Powierzchnia  płaska  lub 

ś

ciana  płaska  jest  w  ogólnym  przypadku  nazywana 

płaszczyzn

ą

  pochylon

ą

  pod  k

ą

tem 

α

  wzgl

ę

dem  poziomu  (rys.  4)  za

ś

  w  szczególnych 

przypadkach – płaszczyzn

ą

 pionow

ą

 lub poziom

ą

)

(

dz

z

U

dy

y

U

dx

x

U

dp

+

+

=

ς

)

(

dz

z

U

dy

y

U

dx

x

U

dp

+

+

=

ς

dz

Z

dy

Y

dx

X

dU

*

*

*

+

+

=

background image

 

 

Rys 4. Parcie cieczy na 

ś

cian

ę

 płask

ą

 

 

W  przypadku  zbiornika  (naczynia)  otwartego,  ci

ś

nienie  atmosferyczne  P

a

  nie  jest 

uwzgl

ę

dniane,  gdy

ż

    działa  ono  jednocze

ś

nie  tak

ż

e  na  zewn

ę

trzn

ą

  stron

ę

 

ś

ciany.  

W  przypadku  zbiornika  (naczynia)  zamkni

ę

tego  nale

ż

y  uwzgl

ę

dni

ć

  ewentualn

ą

  ró

ż

nic

ę

 

ci

ś

nie

ń

, jaka istnieje pomi

ę

dzy czynnikami gazowymi wewn

ą

trz i na zewn

ą

trz.  

Parcie  F  na  płaszczyzn

ę

  pochyła  A,  o 

ś

rodku  ci

ęż

ko

ś

ci  S  i  współrz

ę

dnych 

ś

rodka 

ci

ęż

ko

ś

ci Xs, Ys, Zs, wynosi 

 

gdzie 

∫∫

z*dA=z

S

*A – moment statyczny powierzchni A wzgl

ę

dem powierzchni cieczy 

 

p

– ci

ś

nienie hydrostatyczne na gł

ę

boko

ś

ci z

S

 

ś

rodka ci

ęż

ko

ś

ci S. 

 

 

Powierzchnie zakrzywione 

 

Parcie na powierzchnie zakrzywione sprowadza si

ę

 w ogólnym przypadku do skr

ę

tnika, 

tzn.  parcia  wypadkowego  i  pary  sił,  wzgl

ę

dnie  do  dwóch  sił  sko

ś

nych.  Tak  wi

ę

c  utrzymanie 

zakrzywionej powierzchni w równowadze wymaga przyło

ż

enia nie tylko siły przeciwnej parciu, 

ale i pary sił o odpowiednim momencie. 

Znalezienie  parcia  polega  praktycznie  na  znalezieniu  jego  rzutów  na  kierunki  osi 

współrz

ę

dnych.  W  tym  celu  zostanie  przyj

ę

ty  układ  współrz

ę

dnych  prostok

ą

tnych,  

w którym osie X, Y le

żą

 na swobodnej powierzchni cieczy, a o

ś

 z skierowania jest na dół (rys. 

5).  Element  dA  stanowi  cz

ęść

  powierzchni  cylindrycznej  o  tworz

ą

cych  prostopadłych  do 

płaszczyzny yz. 

 

 

 

 

∫∫

∫∫

=

=

=

=

A

A

S

S

A

p

A

z

g

dA

z

g

dA

p

F

*

*

*

*

*

*

*

ξ

ξ

background image

Rys 5. Parcie cieczy na element powierzchni zakrzywionej 

dF=

ξ

*g*z*dA 

 

Paradoks hydrostatyczny 

 

Składowa pozioma F

y

 parcia na powierzchni

ę

 zakrzywion

ą

 A równa si

ę

 zatem parciu na 

powierzchni

ę

 A

y

, która jest rzutem powierzchni A na płaszczyzn

ę

 pionow

ą

. Tak samo oblicza 

si

ę

 poziom

ą

 składow

ą

 parcia na płaszczyzn

ę

 pochyl

ą

. Sposób obliczania poziomej skł

ą

dowej 

parcia nie ró

ż

ni si

ę

 jak wida

ć

, od sposobu obliczania parcia na płaszczyzn

ę

 pionow

ą

gdzie V, m, G – obj

ę

to

ść

, masa i ci

ęż

ar słupa cieczy nad powierzchni

ą

 zakrzywion

ą

 

gdzie 

M

X

=Z

S

*A 

W  przypadku  płaszczyzn  poziomych  (

α

=0),  zale

ż

no

ś

ci  upraszczaj

ą

  si

ę

Ś

rodek  parcia 

le

ż

y w 

ś

rodku ci

ęż

ko

ś

ci płaszczyzny dna naczynia lub zbiornika. 

 

Rys.6. Paradoks hydrostatyczny 

∫∫

∫∫∫

∫∫

=

=

=

=

=

=

Z

Z

A

V

A

Z

Z

Z

G

g

m

V

g

dV

g

dA

z

g

dF

F

*

*

*

*

*

*

ξ

ξ

ξ

background image

a

a

z

z

y

x

z

y

x

p

gz

p

p

c

z

c

gz

p

c

gz

p

gdz

p

d

dz

f

p

d

g

f

f

f

dz

f

dy

f

dx

f

dp

+

=

=

=

+

+

=

=





=





=

=

=

+

+

=

ρ

ρ

ρ

ρ

ρ

ρ

'

0

'

/*

0

1

 

Jak wida

ć

 z powy

ż

szego równania ci

ś

nienie na dnie zbiornika , a wi

ę

c i napór na jego dno 

nie  zale

ż

y  od  kształtu  zbiornika  ,  ale  od  wysoko

ś

ci  słupa  cieczy.  To  zjawisko  okre

ś

la  si

ę

 

mianem paradoksu hydrostatycznego

Moment siły naporu 

Moment siły naporu wzgl

ę

dem dowolnego punktu definiujemy nast

ę

puj

ą

co: 

(

)

*

S

L

r pn dS

= −

                                                                          (4) 

Gdzie  r  jest  promieniem  ł

ą

cz

ą

cym  punkt,  wzgl

ę

dem  którego  liczymy  moment,   

z  elementem powierzchni dS. 

Tak  wi

ę

c  siły    hydrostatyczne  sprowadzaj

ą

  si

ę

  do  naporu    hydrostatycznego   

i momentu  hydrostatycznego. Jak wiadomo ze statyki, taki układ przestrzennych sił daje si

ę

 

sprowadzi

ć

 do wypadkowej siły wtedy i tylko,wtedy, kiedy 

⊥⊥

 L.                                                                     (5) 

Prostopadło

ść

  momentu  i  naporu  jest  oczywista  dla  płaskich 

ś

cian.  Wówczas   

wektorem r

n

 mo

ż

na okre

ś

li

ć

 poło

ż

enie  linii działania siły naporu 

r

n

*N = L.                                                           (6) 

Punkt przebicia 

ś

ciany lini

ą

   działania siły naporu nazywamy 

ś

rodkiem naporu

Je

ż

eli obliczone s

ą

 ju

ż

 N i L to równanie (6) staje si

ę

 równaniem dla trzech składowych 

r

x

, r

y

, r

z

  wektora r

n

 

0

0

0

x

z y

y z

x

z x

y

x z

y

y x

x y

z

z

r

N r

N r

L

N r

r

N r

L

N r

N r

r

L

+

=

+

+

=

+

=

  

(7) 

background image

7 i 8. Napór hydrostatyczny na 

ś

ciany płaskie i zakrzywione Współrz

ę

dne 

ś

rodka naporu 

 

Pole sił masowych  

k

g

f

=

 

∫ ∫

=

=

=

A

pdA

p

A

pd

p

d

dA

n

A

d

dA

,

,

,

0

 

∫ ∫

=

+

=

+

=

pomijamy

p

gzdA

p

dA

p

gz

p

p

gz

p

a

A

a

A

a

,

,

)

(

,

ρ

ρ

ρ

 

∫ ∫

=

dA

n

gz

p

A

ρ

 

k

k

n

j

j

n

i

i

n

n

)

,

cos(

)

,

cos(

)

,

cos(

+

+

=

 

 Po podstawieniu i uporz

ą

dkowaniu otrzymujemy  

k

p

j

p

i

p

p

z

y

x

+

+

=

 

∫ ∫

=

Ax

x

gzdA

p

,

ρ

 

∫ ∫

=

Ay

y

gzdA

p

,

ρ

 

∫ ∫

=

Az

z

gzdA

p

,

ρ

 

Dla p

X

 , p

y

 obliczamy sił

ę

 naporu 

∫ ∫

=

=

Ax

x

gzdA

p

ρ

 

x

sx

Ax

A

gz

zdA

g

ρ

ρ

∫ ∫

=

 

Analogicznie  

y

sy

y

A

gz

p

ρ

=

 

z

sx

 , z

sy

, -współrz

ę

dne 

ś

rodka ci

ęż

ko

ś

ci rzutu powierzchni w kierunku osi x i y.  

A

x

, A

y

 -pola powierzchni rzutów w kierunku osi x i y. 

Warto

ść

 p

Z

 wyliczamy w sposób odmienny. 

∫ ∫

=

=

z

Az

z

gV

zdA

g

p

ρ

ρ

 

V

z

 -obj

ę

to

ść

 słupa cieczy nad poziomym rzutem powierzchni. 

 

 

 

background image

Ś

rodek naporu 

 

]

,

,

[

n

n

n

n

z

y

x

r

=

 

Sumaryczny moment od naporów elementarnych wynosi  

∫ ∫

×

p

d

r

A

 

Moment od wypadkowej 

dA

n

r

gz

p

d

r

A

n

)

(

×

=

×

∫ ∫

ρ

 

 

x

x

x

AX

x

AX

nx

AX

x

nx

x

x

x

x

sx

x

x

AX

x

nx

M

J

M

dA

z

gM

dA

gz

z

dA

gz

gM

z

gzdA

dp

gM

A

gz

p

zdp

p

z

=

=

=

=

=

=

=

∫ ∫

∫ ∫

∫ ∫

∫ ∫

2

2

2

*

,

*

ρ

ρ

ρ

ρ

ρ

ρ

ρ

 

J

x

  - moment bezwładno

ś

ci 

M

x

 - moment statyczny 

 

x

xy

x

Ax

nx

Ax

x

nx

x

Ax

x

nx

M

J

M

yzdA

y

gyzdA

gM

y

ydp

p

y

=

=

=

=

∫ ∫

∫ ∫

∫ ∫

ρ

ρ

;

*

 

 

 

 

 

 

 

 

9. Zjawisko wyporu. Odkrycie Archimedesa (czas, epoka, miejsce) 

background image

 

Je

ż

eli ciało jest zanurzone w płynie (rys.4.4),wypór hydrostatyczny jest okre

ś

lony wzorem: 

=

S

pndS

N

                                                                 (9.1)

 

 

W  takim  przypadku  całk

ę

  zawart

ą

  w  powy

ż

szym  wzorze  mo

ż

na  przekształci

ć

  

na całk

ę

 obj

ę

to

ś

ciow

ą

=

S

gradpd

N

τ

 

(9.2) 

a po uwzgl

ę

dnieniu równania równowagi powy

ż

sze wyra

ż

enie przyjmuje posta

ć

=

=

S

S

d

gk

fd

N

τ

τ

τ

ρ

τ

ρ

                                               (9.3)

 

Ale 

=

)

(

)

(

s

S

d

τ

τ

τ

 

Jest obj

ę

to

ś

ci

ą

 ciała zanurzonego w płynie. St

ą

d siła wyporu 

N= - 

ρρρρ

gk=-G                                                     

(9.4) 

jest  równa  ci

ęż

arowi  cieczy  wypartej  przez  zanurzone  ciało  i  skierowana  jest  przeciwnie  do 

zwrotu siły ci

ęż

ko

ś

ci. Tak wi

ę

c otrzymali

ś

my znane prawo Archimedesa. 

Siła  wyporu  N  przechodzi  przez 

ś

rodek  ci

ęż

ko

ś

ci  obj

ę

to

ś

ci  płynu  zajmowanej  przez 

zanurzone ciało, i wyra

ż

a si

ę

 wzorem: 

N

r

G

r

L

n

n

××××

====

−−−−

××××

====

)

(

                                                     (9.5)

 

 

 

 
 
 
 

 

Archimedes 

background image

 

Jeden  z  najwybitniejszych  greckich  matematyków  i  fizyków  staro

ż

ytno

ś

ci,  odkrył  wiele  praw 

matematycznych  i  fizycznych,  sformułował  wa

ż

ne  zasady  mechaniczne.  Archimedes  urodził 

si

ę

 w 287 roku przed Chrystusem w Syrakuzach; a o jego 

ż

yciu opowiadaj

ą

 Livius, Polibius. 

Nie sprawował 

ż

adnego urz

ę

du, oddaj

ą

c si

ę

 wył

ą

cznie nauce. Przebywał przez pewien czas 

w  sławnej  akademii  Aleksandryjskiej,  jako  ucze

ń

  matematyka  Konona,  z  którym  utrzymywał 

ź

niej  tak

ż

e  korespondencj

ę

.  Zgin

ą

ł  tragicznie  w  roku  212  przy  zdobywaniu  Syrakuz  przez 

Rzymian pod wodz

ą

 Marcellusa ( podczas drugiej wojny Punickiej ) 

 

Prawo Archimedesa 

 

Ka

ż

de  ciało  zanurzone  w  cieczy  traci  pozornie  na  ci

ęż

arze  tyle,  ile  wynosi  ci

ęż

ar  cieczy 

wypartej  przez  to  ciało  tzn.  ciało  zanurzone  w  cieczy  doznaje  ze  strony  tej  cieczy  parcia  do 

góry, równego co do warto

ś

ci ci

ęż

arowi cieczy wypartej przez to ciało. 

 

Z  czasem  prawo  to  uogólniono  na  gazy  (  i  ciała  sypkie  spełniaj

ą

ce  okre

ś

lone  warunki) 

Jak głosi legenda, Hieron II zamówił dla siebie koron

ę

 z czystego złota. Władca nie dowierzał 

jednak    złotnikowi.  Pos

ą

dzał  go  to, 

ż

e  koron

ę

  wykonał  ze  srebra  i  z  zewn

ą

trz  tylko  pozłocił. 

Zwrócił  si

ę

  wtedy  do  przebywaj

ą

cego  na  jego  dworze  Archimedesa,  aby  ten  sprawdził  jego 

przypuszczenie, nie niszcz

ą

c pi

ę

knej korony. 

Archimedes  długo  my

ś

lał  nad  tym  zadaniem,  niestety  bez  skutku.  Zastanawiał  si

ę

  nad 

tym  nawet  w  k

ą

pieli.  Siedz

ą

c  kiedy

ś

  w  wannie  zauwa

ż

ył, 

ż

e  ciała  zanurzone  

w cieczy wydaj

ą

  si

ę

 l

ż

ejsze. W tym momencie przyszło na

ń

 ol

ś

nienie. Z okrzykiem eureka! 

(ła

ć

.  znalazłem)  Archimedes  pono

ć

  wyskoczył  z  wanny  i  w  stroju  mocno  niekompletnym 

pobiegł przez miasto do swego króla, aby mu zakomunikowa

ć

 o rozwi

ą

zaniu problemu. Je

ś

li 

wi

ę

c  wierzy

ć

  legendzie,  to  dzi

ę

ki  zadaniu  króla  Hierona  Archimedes  odkrył  wa

ż

ne  prawo, 

zwane dzi

ś

 prawem Archimedesa, które stanowi podstaw

ę

 teorii pływania ciał. 

Je

ż

eli  rzeczywi

ś

cie  Archimedes  odkrył  to  prawo  w  wannie,  to  trudno  si

ę

  dziwi

ć

ż

e  był 

zaskoczony  prostot

ą

  metody,  jak

ą

  nale

ż

ało  zastosowa

ć

,  aby  rozwi

ą

za

ć

  postawione  mu 

zadanie.  Srebro  ma  bowiem  g

ę

sto

ść

  prawie  dwa  razy  mniejsz

ą

  ni

ż

  złoto.  Fałszywa  korona 

musiałaby  wi

ę

c  mie

ć

  znacznie  wi

ę

ksz

ą

  obj

ę

to

ść

  ni

ż

  korona  z  czystego  złota  o  tej  samej 

masie  i  -  co  za  tym  idzie  -  wypierałaby  wi

ę

cej  wody,  a  zatem  wi

ę

cej  "  traciłaby  pozornie  na 

ci

ęż

arze  ".  Inaczej  fałszywa  korona  byłaby  w  wodzie  znacznie  l

ż

ejsza  ni

ż

  próbka  czystego 

złota o tej  samej masie. 

background image

Prawo  Archimedesa  jest  najwa

ż

niejszym,  lecz  bynajmniej  nie  jedynym  osi

ą

gni

ę

ciem 

Archimedesa  w  zakresie  hydrostatyki.  Jego  badania  i  spostrze

ż

enia  dotycz

ą

ce  warunków 

równowagi  cieczy  i  warunków  pływania  ciał  legły  u  podstaw  rozwoju  tej  dziedziny.  Oprócz 

praw  fizyki  i  matematyki  Archimedes  odkrył  równie

ż

  podstawow

ą

  zasad

ę

  mechaniki  ciał 

płynnych, zwi

ęź

le sformułowana brzmi: ciało zanurzone w cieczy ulega parciu do góry i traci 

pozornie  na  ci

ęż

arze  tyle,  ile  wa

ż

y  wyparta  prze  ze

ń

  ciecz.  Z  tej  zasady  korzysta  si

ę

  przy 

wyznaczaniu  ci

ęż

aru  wła

ś

ciwego,  mianowicie  dzieli  si

ę

  ci

ęż

ar  bezwzgl

ę

dny  przez  obj

ę

to

ść

równ

ą

  liczbowo  ci

ęż

arowi  wypartej  wody.  Wszystkie  przyrz

ą

dy,  słu

żą

ce  do  pomiaru  ci

ęż

aru 

wła

ś

ciwego,  jak  piknometr,  waga  hydrostatyczna,  wolumometr,  a  dla  cieczy  areometr, 

opieraj

ą

 si

ę

 na zasadzie Archimedesa. 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

background image

10. Stan stateczno

ś

ci pływania. Metacentrum i odległo

ść

 metacentryczna 

 

Na ciało zanurzone w cieczy działa ci

ęż

ar ciała G

1

 i wypór W. Je

ż

eli potraktuje si

ę

 wypór 

W jako wypór ciała całkowicie zanurzonego, to mo

ż

liwe s

ą

 trzy przypadki: 

  G

1

<W, siła (w-G

1

) wypiera ciało do góry powoduj

ą

c jego cz

ęś

ciowe wynurzenie; stan 

równowagi  zostaje  osi

ą

gni

ę

ty  wtedy,  gdy  ci

ęż

ar  ciała  b

ę

dzie  równy  wyporowi 

zanurzonej cz

ęś

ci ciała; w tym stanie równowagi ciało pływa 

  G

1

=W, ciało jest całkowicie zanurzone na dowolnej gł

ę

boko

ś

ci 

  G

1

>W, ciało tonie. 

Stateczno

ść

 pływania jest to zdolno

ść

 powrotu ciała pływaj

ą

cego wychylonego ze stanu 

równowagi do pierwotnego poło

ż

enia. 

 

rys.7. Stateczno

ść

 ciała całkowicie zanurzonego: a) równowaga stała, b) chwiejna,  

c) oboj

ę

tna 

 

Pływanie ciała całkowicie zanurzonych. Na ciało całkowicie zanurzone działaj

ą

 dwie siły; 

wypór  W  i  ci

ęż

ar  G

1

  (rys.  wy

ż

ej).  Punkt  S  oznacza 

ś

rodek  ci

ęż

ko

ś

ci  ciała  zanurzonego  

i  w  ogólnym  przypadku  nie  musi  pokrywa

ć

  si

ę

  ze 

ś

rodkiem  wyporu  N,  który  le

ż

y  w 

ś

rodku 

geometrycznym  ciała.  Równowaga  pływania,  jak  wiadomo,  zachodzi  wówczas  gdy  W=  G

1  

i  gdy  W  i  G

1

  le

żą

  wzdłu

ż

  tej  samej  osi  pionowej,  czyli  wzdłu

ż

  osi  pływania.  Mo

ż

liwe  s

ą

  trzy 

przypadki: 

  Punkt S le

ż

y poni

ż

ej punktu N, 

  Punkt S le

ż

y powy

ż

ej punktu N 

  Punkt S i N pokrywaj

ą

 si

ę

 

Mo

ż

na zatem stwierdzi

ć

ż

e: 

1)  równowaga stała – 

ś

rodek ci

ęż

ko

ś

ci S le

ż

y poni

ż

ej 

ś

rodka wyporu N, 

2)  równowaga chwiejna – S le

ż

y powy

ż

ej N, 

background image

3)  równowaga oboj

ę

tna – punkt S i N pokrywaj

ą

 si

ę

Bardziej zło

ż

onym zagadnieniem jest stateczno

ść

 pływania ciał cz

ęś

ciowo zanurzonych. 

Dowolne wychylenie ciała jest, ogólnie bior

ą

c, wypadkow

ą

 trzech przesuni

ęć

 i trzech obrotów 

wzgl

ę

dem osi X, Y, Z., przy czym o

ś

 x jest prostopadła do płaszczyzny. Przy takim poło

ż

eniu 

ś

rodka  ci

ęż

ko

ś

ci  mo

ż

liwe  jest  zachowanie  stateczno

ś

ci  pływania,  co  było  wykluczone  

w przypadku ciała pływaj

ą

cego całkowicie zanurzonego. 

 Ciało  jest  stateczne,  czyli  posiada  równowag

ę

  stał

ą

,  przy  przesuni

ę

ciu  wzdłu

ż

  osi  z. 

Przy  takiej  wymuszonej  zmianie  gł

ę

boko

ś

ci  zanurzenia  zostaje  naruszona  równowaga 

pomi

ę

dzy  ci

ęż

arem  ciała  G

i  wyporem  W,  co  prowadzi  do  zmiany  zanurzenia  i  powrotu  do 

stanu  pocz

ą

tkowego.  Równowaga  oboj

ę

tna  ma  miejsce,  natomiast,  podczas  przesuni

ęć

 

równoległych  do  zwierciadła  cieczy,  czyli  podczas  przesuni

ęć

  wzdłu

ż

  osi  x  i  y  oraz  podczas 

obrotu wokół osi z.. 

Te  rozwa

ż

ania  mo

ż

na  uzupełni

ć

,  wprowadzaj

ą

c  poj

ę

cia  punktu  M.,  zwanego 

metacentrum, czyli punkt przeci

ę

cia linii działania wyporu chwilowego W i pionowej osi ciała 

pływaj

ą

cego.  Odległo

ść

  punktu  M  od 

ś

rodka  ci

ęż

ko

ś

ci  ciała  S  nosi  nazw

ę

  odległo

ś

ci

ą

 

(wysoko

ś

ci

ą

)  metacentrycznej  m.  Dla  informacji  mo

ż

na  poda

ć

ż

e  minimalna  odległo

ść

 

metacentryczna statków wynosi m = 0,5

÷

4,5 m. 

Odległo

ść

  metacentryczna  m.  mo

ż

na  wyrazi

ć

  przez  parametry  geometryczne  ciała 

pływaj

ą

cego. Przy wychyleniu o mały k

ą

ϕ

, wypór chwilowy W jest równy sumie algebraicznej 

wyporu pocz

ą

tkowego W i wyporów W

k

 obj

ę

to

ś

ci klinowych. 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

background image

11. Kinematyka płynów. Cele, zadania, parametry kinematyczne. 

 

DWIE METODY OPISU STANU PŁYNU 

 

Wprowadzaj

ą

c omówion

ą

 we Wst

ę

pie cech

ę

 ci

ą

gło

ś

ci, mo

ż

emy traktowa

ć

 płyn wypełniaj

ą

cy 

rozwa

ż

any  przez  nas  obszar  jako  kontinuum  materialne.  Ka

ż

demu  punktowi  tego  obszaru 

przyporz

ą

dkowujemy pewne jego otoczenie, o wymiarach małych w porównaniu z wymiarami 

ciał  opływanych  lub  wymiarami  obszaru,  a  du

ż

ych  w  porównaniu  z  odległo

ś

ciami  mi

ę

dzy 

molekułami.  Otoczenie  takie  nazywamy  elementem  płynu,  a  przestrze

ń

  traktujemy  jako 

wypełnion

ą

 tymi elementami w sposób ci

ą

gły. Posługiwanie si

ę

 poj

ę

ciem elementu płynu jest 

bardzo  dogodne  —  mo

ż

emy  uto

ż

sami

ć

  go  z  punktem  któremu  przypisane  s

ą

  wszystkie 

parametry  charakteryzuj

ą

ce  płyn  lub  te

ż

  traktowa

ć

  go  jako  sko

ń

czon

ą

  obj

ę

to

ść

  o  cechach 

kontinuum materialnego. 

Stan  płynu  zajmuj

ą

cego  okre

ś

lony  obszar  przestrzeni  mo

ż

emy  opisa

ć

  w  dwojaki  sposób. 

Pierwszy z nich polega na okre

ś

leniu parametrów w ka

ż

dym punkcie przestrzeni zajmowanej 

przez  płyn.  Parametry  te  mog

ą

  si

ę

  zmienia

ć

  zarówno  w  przestrzeni  (mog

ą

  by

ć

  funkcj

ą

 

współrz

ę

dnych przestrzeni), jak i w czasie. Metoda ta nosi nazw

ę

 metody Eulera albo metody 

lokalnej. 

Omawiaj

ą

c  drugi  sposób,  wyobra

ź

my  sobie, 

ż

e  wyodr

ę

bnili

ś

my  element  płynu  

ś

ledzimy  jego  zachowanie  si

ę

  w  czasie.  Znaj

ą

c  wszystkie  interesuj

ą

ce  nas  parametry 

dowolnych  elementów  płynu,  mamy  opisany  stan  płynu  w  całym  zajmowanym  przez  niego 

obszarze.  

Ta metoda nosi nazw

ę

 metody Lagrange^a lub metody w

ę

drownej. 

 

KINEMATYCZNY PODZIAŁ PRZEPŁYWÓW 

 

Je

ż

eli  parametry  opisuj

ą

ce  stan  płynu  s

ą

  niezale

ż

ne  od  czasu  f,  to  d/dt  =  O  i  stan  taki 

nazywamy  stanem  stacjonarnym  (ustalanym).  Je

ż

eli  za

ś

  parametry  te  zale

żą

  jawnie  od 

czasu,  

to przepływ nazywamy przepływem niestacjonarnym (nieustalonym}. 

Z  punktu  widzenia  kinematycznego  b

ę

dziemy  rozró

ż

niali  przepływy  jednowymiarowe,  

w których mamy tylko jedn

ą

 składow

ą

 wektora pr

ę

dko

ś

ci

 

 

przepływy dwuwymiarowe 

background image

 

oraz przepływy trójwymiarowe 

 

W  ka

ż

dym  z  tych  przypadków  składowe  wektora  pr

ę

dko

ś

ci  mog

ą

  by

ć

  funkcj

ą

  dowolnej 

kombinacji  zmiennych  niezale

ż

nych  x,  y,  z,  t.  Tak  na  przykład,  mo

ż

emy  mie

ć

  przepływ 

jednowymiarowy niestacjonarny, w którym 

u

x

 = u

x

(x,f), 

lub przepływ dwuwymiarowy stacjonarny 

u

x

 = u

x

 (x,y)u

y

 = u

y

 (x,y) 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

background image

12. Metoda Eulera i Lagrange’a w kinematyce płynów 

 

 

METODA EULERA 

 

Je

ż

eli  wybrany  zostanie  układ  współrz

ę

dnych  odniesienia,  to  obszar  zajmowany  przez 

płyn  b

ę

dzie  opisany  promieniem  wodz

ą

cym  r.  W  najprostszym  przypadku,  je

ż

eli  układ 

współrz

ę

dnych b

ę

dzie układem kartezja

ń

sktm, to 

 

Wówczas  parametry  opisuj

ą

ce  stan  płynu  b

ę

d

ą

  funkcjami  promienia  wodz

ą

cego  r 

i czasu t. Ka

ż

demu poło

ż

eniu elementu płynu w chwili r, które jest okre

ś

lone współrz

ę

dnymi 

x,  y,  z,  odpowiada

ć

  b

ę

dzie  poło

ż

enie  *„,  y

(i

,  z<>  w  chwili  poprzedzaj

ą

cej  t

t

,.  Zapisujemy  t

ę

 

odpowiednio

ść

 w nast

ę

puj

ą

cej postaci: 

                                          (12.1) 

metodzie Eulera mamy nast

ę

puj

ą

cy zapis dla wektora pr

ę

dko

ś

ci: 

 

co oznacz,a, 

ż

e wektor pr

ę

dko

ś

ci mo

ż

e by

ć

 przedstawiony nast

ę

puj

ą

co: 

            (12.2)                  

lub w skrócie

 

Nie  tylko  składowe  wektora  pr

ę

dko

ś

ci  s

ą

  funkcjami  współrz

ę

dnych  przestrzennych  i  czasu, 

ale i inne parametry, takie jak: ci

ś

nienie, g

ę

sto

ść

, temperatura 

                                              (12.3) 

Je

ś

li posługujemy si

ę

 t

ą

 metod

ą

 opisu parametrów płynu, to w dowolnym miejscu układu 

współrz

ę

dnych  i  w  dowolnym  czasie  powinny  by

ć

  okre

ś

lone  powy

ż

sze  funkcje  czterech 

zmiennych x, y, z, t. 

Istotnym  poj

ę

ciem  dla  metody  Eulera  jest  poj

ę

cie  pochodnej  substancjalnej,  oznaczonej  dla 

dowolnej funkcji f (x, y, z, t) symbolem d//dr. 

Pochodn

ą

 substancjaln

ą

 buduje si

ę

 bior

ą

c za punkt wyj

ś

cia poj

ę

cie ró

ż

niczki zupełnej funkcji 

wielu zmiennych. W tym przypadku

 

background image

k

z

j

y

i

x

+

+

=

                            (12.4) 

W wyra

ż

eniu tym przyrosty dx, dy i dz s

ą

 przyrostami dowolnymi w przestrzeni ryz. Je

ż

eli na 

przyrosty te nało

ż

ymy ograniczenia 

                                                    (12.5) 

co oznacza, 

ż

e s

ą

 one wybierane wzdłu

ż

 kierunku ruchu cz

ą

stki, to wyra

ż

enie (12.4) mo

ż

na 

zapisa

ć

 nast

ę

puj

ą

co: 

                                 (12.6) 

Odnosz

ą

c przyrost d/ do przyrostu czasu dt, z (12.6) otrzymujemy 

                            (12.7) 

Zapis ten mo

ż

e by

ć

 zastosowany do dowolnej funkcji/. Istotny jest jedynie operator typu: 

    (12.8) 

Ze sposobu budowania operatora pochodnej substancjalnej wynika nast

ę

puj

ą

ca interpretacja 

fizyczna poszczególnych wyra

ż

e

ń

d/df oznacza zmian

ę

 danej wielko

ś

ci w czasie z punktu widzenia obserwatora poruszaj

ą

cego 

si

ę

 wraz z elementem płynu, 

d/dt oznacza zmian

ę

 w czasie danej wielko

ś

ci w danym punkcie przestrzeni (przy ustalonym 

x, y, z) — jest to pochodna lokalna,

 

 

oznacza zmian

ę

 danej wielko

ś

ci w przestrzeni w danym 

ustalonym czasie t — jest to pochodna konwekcyjna (cz

ę

sto tak

ż

e zwana adwekcyjn

ą

). 

Tak    wi

ę

c      pochodna    substancjalna  jest      sum

ą

    pochodnej      lokalnej      i    pochodnej 

konwekcyjnej (adwekcyjnej). 

Mo

ż

na  łatwo  zauwa

ż

y

ć

,   

ż

e  pochodna konwekcyjna jest  iloczynem   skalarnym 

wektora pr

ę

dko

ś

ci u i operatora gradient (nabla) grad =

  

                                  (12.9) 

Stosuj

ą

c  operator  ró

ż

niczkowania  substancjalnego  do  składowych  wektora  pr

ę

dko

ś

ci, 

otrzymujemy przyspieszenie substancjalne 

background image

      (12.10; 12.11; 12.12) 

Po  pomno

ż

eniu  powy

ż

szych  równa

ń

  przez  wektory  i,  j,  k  oraz  ich  zsumowaniu  mo

ż

emy 

relacje (1.10)—(1.12) zapisa

ć

 w sposób zwarty 

    (12.13) 

Podobnie, stosuj

ą

c operator (1.8) do innych parametrów, otrzymujemy 

       (12.14; 12.15; 12.16) 

Z  powy

ż

szych  wzorów  wida

ć

,  i

ż

  zmiana  parametru  w  elemencie  płynu  okre

ś

lona  jest 

zmienno

ś

ci

ą

  w  czasie  i  zmienno

ś

ci

ą

  w  przestrzeni,  co  wymaga  znajomo

ś

ci  pola  pr

ę

dko

ś

ci. 

Ponadto  człony  okre

ś

laj

ą

ce  pochodn

ą

  konwekcyjn

ą

  s

ą

  członami  nieliniowymi  typu  iloczynu 

pr

ę

dko

ś

ci i pochodnej danego parametru. 

Nieliniowo

ść

 ta, pojawiaj

ą

ca si

ę

 w podstawowych relacjach mechaniki płynów, jest przyczyn

ą

 

znacznych trudno

ś

ci w rozwi

ą

zywaniu problemów przepływowych. 

 

METODA LAGRANGE'A 

 

W  metodzie  tej  opisujemy  histori

ę

  zmiany  danego  parametru,  zwi

ą

zanego  z  wybranym 

elementem  płynu.  Oznaczamy  wybrany  element  płynu  w  chwili  f

(l

  współrz

ę

dnymi  x

0

,  Y

0

,  Z

0

.  

Z  czasem  b

ę

dzie  si

ę

  zmieniało  jego  poło

ż

enie  i  inne  parametry.  Zapisujemy  ten  fakt 

nast

ę

puj

ą

co: 

                                     (12.17) 

lub w zapisie wektorowym 

 

Dla ci

ś

nienia, g

ę

sto

ś

ci i temperatury mamy 

background image

                                             (12.18) 

Pr

ę

dko

ś

ci  i  przyspieszenia  w  metodzie  Lagrange'a  otrzymujemy  przez  ró

ż

niczkowanie 

współrz

ę

dnych poło

ż

enia elementu płynu wzgl

ę

dem czasu. Tak wi

ę

c dla pr

ę

dko

ś

ci 

                              (12.19) 

lub w zapisie wektorowym 

 

a dla przyspiesze

ń

 

                          (12.20) 

W dalszym ci

ą

gu posługiwa

ć

 si

ę

 b

ę

dziemy przewa

ż

nie metod

ą

 Eulera. 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

13.Pochodna w

ę

drowna. Operator Stokesa, przy

ś

pieszenie elementu 

płynu. 

 

background image

 

Je

ś

li  wezmiemy  pod  uwag

ę

  dowoln

ą

  wielko

ść

  fizyczn

ą

  f,  skalarn

ą

  lub  wektorow

ą

zwi

ą

zan

ą

  z  poruszaj

ą

cym  si

ę

  elementem  płynu,  to  wyra

ż

enie  zmian  tej  wielko

ś

ci  w  czasie 

przy  u

ż

yciu  zmiennych  Eulera  napotka  na  pewne  trudno

ś

ci.  Mamy  wyrazi

ć

  zmiany  f, 

zwi

ą

zane  z  poruszaj

ą

cym  si

ę

po  swym  torze  elementem  płynu,  przy  pomocy  wielko

ś

ci 

zwi

ą

zanych  z  wybranym  punktem  przestrzeni  x,y,z  i  wybran

ą

  chwil

ę

  czasu  t.  Mamy  zatem 

obliczy

ć

 dla elementu płynu 

 

Chc

ą

c  to  wyrazi

ć

  za  pomoc

ą

  operacji  ró

ż

niczkowania  wzgl

ę

dem  x,y,z,t  musimy  x,y,z  w 

powy

ż

szym  wzorzepotraktowa

ć

  jako  równanie  parametryczne  toru  elementu  płynu 

przechodz

ą

cego przez wybranu punkt przestrzeni, a wi

ę

c jako: 

 

Uwzgl

ę

dniaj

ą

c to piszemy 

 

Korzystaj

ą

c ze zwi

ą

zków dla toru elementu płynu otrzymujemy 

 

Przy u

ż

yciu operatora mo

ż

emy ten wzór zapisa

ć

 nast

ę

puj

ą

co 

 

Wyra

ż

enie 

 uwa

ż

amy formalnie za iloczyn skalarny 

 i wektora pr

ę

dko

ś

ci 

 

Wynik tej operacji traktujemy jako operator, którym działamy na dowoln

ą

 wielko

ś

c skalarn

ą

 

np. 

 lub wektorow

ą

 np. 

 

 

Gdy podstawimy 

to równanie oznacza przy

ś

pieszenie elementu płynu 

 

 

Składowe przy

ś

pieszenia w układzie współrz

ę

dnych kartezja

ń

skich s

ą

 okre

ś

lone wzorami 

background image

 

Pochodn

ą

 df/dt dowolnej wielko

ś

ci f nazywamy pochodn

ą

 substancjaln

ą

. Ma to okre

ś

lony 

sens fizyczny, gdy

ż

 d/dt oznacza zmiany dla tego samego poruszaj

ą

cego si

ę

 elementu płynu, 

czyli zmiany zwi

ą

zane z jego „substancj

ą

’’. 

 

We  współrz

ę

dnych  Eulera  zmiany  dotycz

ą

ce  poruszaj

ą

cego  so

ę

  elementu  trzeba 

wyrazi

ć

 przy pomocy pochodnych cz

ą

stkowych wzgl

ę

dem czasu i pkt. Przestrzeni. Pochodna 

substancjalna  jest  sum

ą

  pochodnej  lokalnej 

f  / 

t  wyra

ż

aj

ą

cej  zmiany  wielko

ś

ci  f  w  czasie 

lecz  w  tym  samym  punkcie  przestrzeni  i  pochodnej  konwekcyjnej. 

wyra

ż

aj

ą

cej 

zmiany wielko

ś

ci f przy przej

ś

ciu z pkt. X,y,z do jego najbli

ż

szego otoczenia, co zwi

ą

zane jest 

z polem pr

ę

dko

ś

ci.

 

 

14.Poj

ę

cie cyrkulacji. Interpretacja fizyczna i analogia. 

 
 
 
Cyrkulacja  pr

ę

dko

ś

ci  wzdłu

ż

  odcinka  BC  krzywej  s  w  ogólnej  krzywej  przestrzennej,  jest  to 

całka liniowa z iloczynu skalarnego 

i skierowanego elementu tej lini 

czyli

 

 

Cyrkulacja pr

ę

dko

ś

ci wyra

ż

ona równaniem wcze

ś

niejszym mo

ż

e by

ć

 równie

ż

 przedstawiona 

w nast

ę

puj

ą

cy sposób: 

background image

 

gdzie 

- rzut wektora 

na kierunek 

 

W  przypadku  całki  liniowej  po  krzywej  zamkni

ę

tej  s  ,  czyli  całki 

  cyrkulacja  pr

ę

dko

ś

ci 

 wynosi: 

 

Jako dodatni kierunek obchodu zamkni

ę

tej linii konturowej przy obliczaniu całki przyjmiejmy  

kierunek przeciwny do ruch wskazówek zegara. 

 

Poj

ę

cie cyrkulacji pr

ę

dko

ś

ci wyst

ę

puje w zagadnieniach opływu ciał, profili, łopatek itd. 

Poj

ę

cie to jest podobne do poj

ę

cia pracy z tym, 

ż

e zamiast wektora siły wyst

ę

puje wektor 

pr

ę

dko

ś

ci. 

  

 

 

 

 

Cyrkulacja pr

ę

dko

ś

ci wzdłu

ż

 zewn

ę

trznej lini konturowej równa si

ę

 sumie cyrkulacji wzdłu

ż

 

konturów składowych. Zilustrowano to na rys. gdzie lini

ę

 konturow

ą

 s podzielono odcinkiem 

BD na dwie linie konturowe. Zatem mo

ż

na napisa

ć

 

ż

e: 

 

 

background image

 

gdy

ż

 cyrkulacja wzdłu

ż

 BD znosz

ą

 si

ę

 z uwagi na przeciwne znaki. 

 

W teorii pola dowodzi si

ę

 tw. Stokesa, Które dotyczy zwi

ą

zku mi

ę

dzy całk

ą

 liniow

ą

 i 

powierzchniow

ą

. Dla ró

ż

niczkowalnych pól wektorowych 

ma ono posta

ć

 

 

gdzie 

 oznacza elementarny skierowany odcinek zamkni

ę

tej lini s stanowi

ą

cej lini

ę

 

konturow

ą

 powierzchni A. 

 

Inaczej tw. Stokesa mo

ż

na napisa

ć

       

 

Twierdzenie  formułowane  w  kategoriach  kinematycznych  mówi 

ż

e  cyrkulacja  pr

ę

dko

ś

ci 

wzdłu

ż

 zamkni

ę

tej krzywej s równa si

ę

 strumieniowi wirowa

ń

 (rotacji) przechodz

ą

cemu przez 

powierzchni

ę

 A, której brzegiem jest krzywa s. 

 

15.Równania toru i lini pr

ą

du, rurka i powierzchnia pr

ą

du, struga. 

 

Linia pr

ą

du jest to linia pola wektorowego pr

ę

dko

ś

ci. Je

ś

li mamy pole pr

ę

dko

ś

ci 

 

o składowych w układzie Kartezja

ń

skim 

       

 

oraz  element  linii  pr

ą

du 

o  skł.  Dx,  dy,  dz  to  warunek  równoległo

ś

ci 

  w  ka

ż

dym 

punkcie pola dla dowolnej chwili czasu mo

ż

emy napisa

ć

 w postaci  

 

czyli 

background image

  

 

 

ka

ż

da składowa tego iloczynu musi by

ć

 równa zeru i st

ą

d otrzymujemy równanie lini pr

ą

du 

   

 

W równaniu tym czas wyst

ę

puje jako parametr, od którego zale

żą

 warto

ś

ci Vx, Vy, Vz ale nie 

jest on zmienn

ą

 niezale

ż

n

ą

. W ruchu nieustalonym obraz linii pr

ą

du ma charakter chwilowy, 

zale

ż

ny od warto

ś

ci t. W ruchu ustalonym obraz linii pr

ą

du dla danego przepływu jest 

niezmienny w czasie, tzn. kształt ka

ż

dej linii pr

ą

du dla przepływu ustalonego nie zmienia si

ę

 

w czasie. 

 

Powierzchnia pr

ą

du jest to powierzchnia utworzona z linii pr

ą

du przecinaj

ą

cych 

dowoln

ą

 lini

ę

 l nie b

ę

d

ą

c

ą

 lini

ą

 pr

ą

du. Je

ś

li linia jest zamkni

ę

ta, powierzchnia pr

ą

du 

nazywana jest rurk

ą

 pr

ą

du. 

Zbiór linii pr

ą

du wypełniaj

ą

cych w sposób ci

ą

gły rurk

ę

 pr

ą

du nazywamy włoknem pr

ą

du lub 

strug

ą

 pr

ą

du. 

 

 

 

 

 

Torem elementu płynu nazywamy lini

ę

, po której porusza si

ę

 Element płynu dV traktowany 

jako pkt materialny. Gdy elementarny odcinek toru oznaczymy przez 

o składowych dx, 

dy, dz a elementarny odcinek czasu przez dt, to równanie ró

ż

niczkowe toru jest nast

ę

puj

ą

ce: 

 

Równania skalarne toru elementu wynikaj

ą

ce ze wzoru to: 

background image

 

 

     

 

Czas odgrywa tu rol

ę

 zmiennej niezale

ż

nej. Ostatnie równanie ró

ż

niczkowe toru wyra

ż

ona za 

pomoc

ą

 składowych mal

ą

 posta

ć

 

16.Przepływy potencjalne. 

 

 

Prawie  zawsze  mo

ż

na  potraktowa

ć

  w  przybli

ż

eniu  ka

ż

dy  przepływ  przestrzenny  jako 

przepływ  dwuwymiarowy  (płaski  lub  osiowo  –  symetrycznu).  Takie  uproszczenie  jest 
niezmiernie  korzystne  ze  wzgl

ę

dów  matematycznych,  gdy

ż

  pozwala  stosowa

ć

  bardzo 

dogodn

ą

  i  dobrze  opracowan

ą

  teori

ę

  funkcji  zmiennej  zespolonej.  W  niniejszych 

rozwa

ż

aniach b

ę

d

ą

 omawiane tylko płaski przepływy potencjalne.\ 

 

Mimo, 

ż

e  poj

ę

cie  potencjału  ma  charakter  abstrakcyjny,  to  przy  jego  pomocy 

rozwi

ą

zuje  si

ę

  szereg  wa

ż

nych  zagadnie

ń

  praktycznych.  W  dalszej  tre

ś

ci  b

ę

d

ą

  podane 

przykłady takich rozwi

ą

za

ń

Zale

ż

no

ś

ci podstawowe. 

 

Niektóre podst. poj

ę

cia przepływów potencjalnych, jak potencjał pr

ę

dko

ś

ci 

Φ

 i 

równanie Laplace’a, zostały podane przy opisie pola bezwirowego. 

 

Równania Laplace’a, tj. równanie dla układu płaskiego ma nast

ę

puj

ą

c

ą

 posta

ć

 

background image

Jak wiadomo, wa

ż

n

ą

 cech

ą

 równania Laplace’a jest jego liniowo

ść

, co jest wykorzystywane 

przy superpozycji, czyli nakładaniu przepływów. Zagadnienie superpozycji b

ę

dzie omawiane 

w dalszej cz

ęś

ci tre

ś

ci. 

 

Przyrost potencjału pr

ę

dko

ś

ci mo

ż

e by

ć

 wyra

ż

ony jako ró

ż

niczka zupełna: 

 

 

 

 

17.Przepływy elementarne. Superpozycja przepływów 

 

1. Przepływy elementarne. 

a. Przepływ jednostajny.

 

 

ogólnym 

przypadku, 

jednostajny 

przepływ 

potencjalny 

odbywa 

si

ę

  

pr

ę

dko

ś

ci

ą

 

skierowan

ą

 

pod 

k

ą

tem 

a 

wzgl

ę

dem 

jednej 

osi 

układu,  

np. wzgl

ę

dem osi x. Składowe pr

ę

dko

ś

ci v

x

 i v

y

 wynikaj

ą

 z zale

ż

no

ś

ci trygonometrycznych. 

 

 

Potencjał pr

ę

dko

ś

ci 

Φ

 obliczymy z równania 

Otrzymuj

ą

dy

y

dx

x

d

Φ

+

Φ

=

Φ





Φ

+

Φ

=

Φ

dy

y

dx

x

background image

 

 

Z równa

ń

 

ν

=

x

∂Φ

 

ν

y

y

∂Φ

   jak i rysunku wynika: 

* cos

x

x

ν

α

∂Φ =

 

* sin

y

y

ν

α

∂Φ =

 

st

ą

Funkcj

ę

 pr

ą

du obliczymy z równania 

d

dx

dy

x

y

∂Ψ

∂Ψ

Ψ =

+

 

w

i

edz

ą

c, 

ż

,

y

x

x

y

∂Ψ ∂Φ

∂Ψ ∂Φ

=

=

, otrzymamy  

(

)

* sin

* cos

x

y

ν

α

α

Ψ = −

+

 

Podstawiaj

ą

c do równania 

( )

f x

i

= Φ + Ψ

 otrzymuje si

ę

 warto

ść

 potencjału zespolonego 

( ) (

)

(

)

(

)(

)

* cos

*

sin

*

cos

sin

* *

i

f z

x i y

i x y

x i y

i

z e

α

ν

α

α ν

α

α ν

=

+

+ − +

=

+

=

(17.1)

 

Potencjał zespolony wyra

ż

ony równaniem ( ) ma zatem posta

ć

 

f(z)=a*z 

 

 

 

 

 

(17.2) 

gdzie a jest liczb

ą

 zespolon

ą

Jest to potencjał zespolony przepływu jednostajnego.  

 

b. 

Ż

ródlo płaskie. 

 

Ź

ródło  płaskie  mo

ż

e  by

ć

  dodatnie,  lub  ujemne. 

Ź

ródło  dodatnie  nazywane  jest  po 

prostu 

ź

ródłem,  za

ś

 

ź

ródło  ujemne  nosi  nazw

ę

  upustu.  Potencjał  zespolony  wyra

ż

a  si

ę

 

równaniem 

f(z)=a*lnz                        

 

 

(a) 

gdzie: 

=

2

V

Π

 - liczba rzeczywista

,

 

V-wydajno

ść

 

ź

ródła, lub upustu płaskiego.  

(

) (

)

α

α

ν

α

ν

α

ν

sin

*

cos

*

sin

*

*

cos

*

y

x

dy

dx

+

=

+

=

Φ

background image

W  układzie  współrz

ę

dnych  biegunowych  potencjał  zespolony,  przy  wykorzystaniu 

równania z = r * e

i

ϕ

 przyjmie posta

ć

    

( )

* ln

ln( *

)

(ln *

)

2

2

i

V

V

f z

i

a

z

r e

r i

ϕ

ϕ

= Φ + Ψ =

=

=

Π

Π

             

(17.3)  

Potencjał pr

ę

dko

ś

ci jako cz

ęść

 rzeczywista wynosi 

ln

2

V

r

Φ =

Π

 

Linie 

stałego 

potencjału 

pr

ę

dko

ś

ci 

s

ą

 

okr

ę

gami 

współ

ś

rodkowymi, 

gdy

ż

 

dla  

Φ

 = const jest r = const. 

Funkcja pr

ą

du jako cz

ęść

 urojona wynosi 

ln( *

)

2

i

V

r e

ϕ

Π

 

gdzie 

ϕ

 -k

ą

t(współrz

ę

dna biegunowa). 

Linie  pr

ą

du  (

Ψ

  =  const)  s

ą

  ortogonalne  do  linii 

Φ

  =  const,  a  wi

ę

c  s

ą

  prostymi 

wychodz

ą

cymi z pocz

ą

tku układu współrz

ę

dnych (

ϕ

 = const). 

ź

ródle i upu

ś

cie istnieje tylko pr

ę

dko

ść

 promieniowa 

ν

r

 gdy

ż

 pr

ę

dko

ść

 obwodowa V

u

 

jest równa zero. Zwrot pr

ę

dko

ś

ci 

ν

r

 w 

ź

ródle i upu

ś

cie jest odwrotny, przy czym przyjmuje si

ę

ż

e  w  niesko

ń

czenie  bliskim  otoczeniu  punktu 

ś

rodkowego  płynu  pojawia  si

ę

  lub  przestaje 

istnie

ć

Warto

ść

 pr

ę

dko

ś

ci 

ν

r

 mo

ż

na znale

źć

 w nast

ę

puj

ą

cy sposób. Ró

ż

niczkuj

ą

c równanie (a)  

i wykorzystuj

ą

c równania z=r*e

i

ϕ

 i 

( )

1

*

2

i

df z

V

e

dz

r

ϕ

=

Π

otrzymujemy kolejno 

( )

1

1

*

*

2

2

i

i

df z

V

V

e

e

dz

r

r

ϕ

ϕ

=

=

Π

Π

 

( )

1

*

2

i

df z

V

e

dz

r

ϕ

=

Π

 

Pr

ę

dko

ść

 

ν

r

 jako rzeczywista równania wynosi 

1

*

2

V

Vr

r

=

Π

, 

Pr

ę

dko

ść

 

ν

u

 jest równa zero gdy

ż

 cz

ęść

 urojona pr

ę

dko

ś

ci zespolonej nie istnieje. 

 

 

 

 

background image

c. Wir płaski potencjalny. 

 

W  ruchu  cieczy  nazywamy  kołowym,  płaskim  wirem  potencjalnym(swobodnym)  cyrkulacja 

pr

ę

dko

ś

ci  po  dowolnej  krzywej  nie  obejmuj

ą

cej 

ś

rodka  wiru  jest  równa  zero-st

ą

d  wir 

potencjalny. Cyrkulacja po krzywej jeden raz obejmuj

ą

cej 

ś

rodek wiru jest równa F. 

Dla kołowego wiru płaskiego warto

ść

 potencjału zespolonego wynosi 

f(z)=-i*a*lnz,                                                     (17.4) 

gdzie: a=

2

V

Π

 -liczba rzeczywista, 

Γ

 - cyrkulacja pr

ę

dko

ś

ci. 

W  układzie  współrz

ę

dnych  biegunowych  potencjał  zespolony,  uwzgl

ę

dniaj

ą

c  powy

ż

sze 

warunki przyjmie posta

ć

 

ƒ

( )

+

Φ

=

z

i

Ψ

(

)

(

)

r

i

i

r

i

e

r

i

z

a

i

i

ln

*

2

2

ln

2

*

ln

2

ln

*

*

π

ϕ

π

ϕ

π

π

ϕ

Γ

Γ

=

+

Γ

=

Γ

=

=

   (17.5) 

Potencjał pr

ę

dko

ś

ci jako cz

ęść

 rzeczywista wynosi  

ϕ

π

2

Γ

=

Φ

   

 

Linie  stałego  potencjału  pr

ę

dko

ś

ci  s

ą

  prostymi  wychodz

ą

cymi  z  pocz

ą

tku  układu 

współrz

ę

dnych, gdy

ż

 dla 

Φ

 = const jest 

φ

 = const 

Funkcja pr

ą

du jako cz

ęść

 urojona wynosi : 

r

1

*

2

π

Γ

=

Φ

 

 

Linie  pr

ą

du  s

ą

  okr

ę

gami  współ

ś

rodkowymi,  gdy

ż

  dla  Y  =  const  jest  r  =  const.  

Role 

φ

 i 

Ψ

 s

ą

 dla 

ź

ródła i wiru zamienne. W zwi

ą

zku z tym, pr

ę

dko

ść

 promieniowa jest równa 

zett| a pr

ę

dko

ść

 obwodowa wynosi : 

r

u

1

*

2

π

ν

Γ

=

    

 

 

 

 

(17.6) 

Ogólnie  mo

ż

na  powiedzie

ć

ż

e  zamiana  ról 

Φ

  i 

Ψ

  zachodzi  w  tzw.  Przepływach 

sprz

ęż

onych, którymi s

ą

 w danym przypadku 

ź

ródło i wir płaski.  

 

d. Dipol. 

 

Dipol,  lub 

ź

ródło  podwójne  płaskie  ma  miejsce  wówczas,  gdy  odległo

ść

  

/  pomi

ę

dzy 

ź

ródłem  i  upustem  wynosi  l 

 

  .Obowi

ą

zuje  przy  tym  taka  zale

ż

no

ść

ż

background image

odległo

ść

  l  maleje  w  takim  stopniu,  w  jaki  wzrasta  strumie

ń

  obj

ę

to

ś

ci  V 

ź

ródła  

i  upustu,  czyli  l 

  O,  V 

 

  ,  przy  czym  l*V=M  =  const.  Moment  dipola  M.  zostaje  stały  

i ró

ż

ny od zera. 

Potencjał zespolony dipola wynosi 

z

a

z

f

=

)

(

 

 

 

 

 

 

(17.7) 

gdzie 

π

2

M

a

=

 

Po  uwzgl

ę

dnieniu  równania  z=x+z  *  y  i  pomno

ż

eniu  licznika  i  mianownika  przez  

(x-i * y) otrzymuje si

ę

 : 

)

(

2

)

*

(

)

*

(

2

)

*

(

)

(

2

2

2

2

2

y

x

y

i

x

M

y

i

x

y

i

x

M

z

f

+

=

=

π

π

  

 

 

 

(17.8) 

Potencjał pr

ę

dko

ś

ci jako cz

ęść

 rzeczywista potencjału zespolonego wynosi : 

 

)

(

2

*

2

2

y

x

x

M

+

=

Φ

π

 

 

 

 

 

 

 

(17.9) 

 

Funkcja pr

ą

du jako cz

ęść

 urojona wynosi  

)

(

2

*

2

2

y

x

x

M

+

=

Ψ

π

  

 

 

 

 

 

 

(17.10) 

Z tych zale

ż

no

ś

c

i

 wynika, 

ż

e linie pr

ą

du ( 

Ψ

 = const ) stanowi

ą

 rodzin

ę

 okr

ę

gów stycznych do 

osi  x  w  pocz

ą

tku  układu,  za

ś

  linie 

Φ

 

=

  const  s

ą

  tak

ż

e  okr

ę

gami,  lecz  stycznymi  w  pocz

ą

tku 

układu do osi y. 

Okr

ę

gi   

Φ

 

=  const  i 

Ψ

  =  const  s

ą

  wzajemnie  ortogonalne.  Osi

ą

  dipola  jest  wspólna 

styczna okr

ę

gów  

Ψ

 =  const

.

 Jej zwrot jest zgodny z kierunkiem przepływu. 

 

 

 

2. Superpozycja przepływów. 

 

Wa

ż

n

ą

  cech

ą

  równania  Laplace'a  jest  jego  liniowo

ść

  wykorzystywana  przy 

superpozycji, 

czyli 

nakładaniu 

dwóch, 

lub 

wi

ę

cej 

przepływów 

potencjalnych.  

W przypadku dwóch przepływów mo

ż

na zapisa

ć

 

1

1

1

*

)

(

Ψ

+

Φ

=

i

z

f

 

2

2

2

*

)

(

Ψ

+

Φ

=

i

z

f

  

 

 

                        (17.11) 

background image

sk

ą

)

(

*

)

(

*

)

(

2

2

1

1

z

f

k

z

f

k

z

f

+

=

   

           (17.12) 

gdzie k

1

 , k

2  

- dowolne współczynniki. 

 

 

 

W  wyniku  superpozycji  został  otrzymany  nowy  przepływ  o  potencjale  zespolonym 

f(z), którego cz

ęś

ci

ą

 rzeczywist

ą

 i urojon

ą

 s

ą

 odpowiednio 

Φ

 i 

Ψ

 

2

2

1

1

*

*

Φ

+

Φ

=

Φ

k

k

 

2

2

1

1

*

*

Ψ

+

Ψ

=

Ψ

k

k

    

 

 

 

    (17.13) 

Superpozycja 

mo

ż

by

ć

 

wykonywana 

metod

ą

 

analityczn

ą

,

 

lub 

wykre

ś

ln

ą

.  

W  metodzie  analitycznej  zadanie  sprowadza  si

ę

  do  znalezienia  wypadkowego                                                                         

(sumarycznego)  potencjału  zespolonego,  a  zwłaszcza  cz

ęś

ci  urojonej  tego  potencjału. 

W metodzie wykre

ś

lnej natomiast nale

ż

y nanie

ść

 linie 

Ψ

1

,

Ψ

2

=const, zachowuj

ą

c mi

ę

dzy nimi 

takie  odleg

ł

o

ś

ci,  jakie  wynikaj

ą

  z  równo

ś

ci  strumieni  obj

ę

to

ś

ci  obu  przepływów,  

a wi

ę

c : 

2

1

V

V

=

 

 

 

 

 

 

(17.14) 

Utworzona w ten sposób siatka linii pr

ą

du mo

ż

e by

ć

 oznakowana w sposób umowny. 

 

 

Zagadnienie  superpozycji  sprowadza  si

ę

  wła

ś

ciwie  do  sumowania  wektorów 

pr

ę

dko

ś

ci. 

Szczególne  znaczenie  w  aerodynamice  posiadaj

ą

  przypadki  superpozycji  przepływu 

jednostajnego i 

ź

ródła, lub układu 

ź

ródeł. 

 

 

 

 

 

background image

∫∫

∫∫∫

=

A

V

dV

v

p

div

dA

n

v

p

)

*

(

*

*

*

r

r

r

n

r

∫∫∫

V

dV

dt

t

p

*

18.Zasada zachowania masy. Równanie ci

ą

gło

ś

ci 

 

Równanie ci

ą

gło

ś

ci przepływu w wynika z zasady zachowania masy (niezniszczalno

ś

ci 

materii).      Równanie  ci

ą

gło

ś

ci  odnosi  si

ę

  tylko  do  pól  bez

ź

ródłowych  tzn. 

ż

e  w 

ż

adnym 

punkcie pola masy nie mo

ż

e si

ę

 tworzy

ć

 ani znika

ć

Równanie ci

ą

gło

ś

ci przepływu ma charakter kinematyczny i z tej racji jest identyczny dla  

płynów    nielepkich  i  lepkich.  Warunek  ci

ą

gło

ś

ci  ma  jednak  inny  sens  dla    płynów  

nie

ś

ci

ś

liwych i 

ś

ci

ś

liwych. 

Dla  płynów  nie

ś

ci

ś

liwych  (p=const),  w  obj

ę

to

ś

ci  kontrolnej  musi  by

ć

  zawarta  

w  ka

ż

dej  chwili  tak  sama  masa    płynu,  czyli  przez  powierzchni

ę

  kontroln

ą

  musi  dopływa

ć

  

i  wypływa

ć

  taka  sama  masa  płynu.  Warunek  ten  jest  taki  sam  dla  przepływów  ustalonych  

i nieustalonych. 

W  przepływie  ustalonym  musi  by

ć

  zachowany  powy

ż

szy  warunek,  bo  masa  zawarta 

wewn

ą

trz  powirzchni  kontrolnej  jest  niezmienna  w  czasie.  W  przepływie  nieustalonym 

natomiast,  mog

ą

  mie

ć

  miejsca  lokalne  zmiany  g

ę

sto

ś

ci  płynu,  które  wywołuj

ą

  ró

ż

nic

ę

  

w bilansie dopływaj

ą

cym przez powierzchni

ę

 kontroln

ą

 strumieni masy. 

 

1.   Sformułowanie ogólne 

 

Równanie  ci

ą

gło

ś

ci  zostanie  wyprowadzone  dla  trójwymiarowego  nieustalonego 

przepływu 

ś

ci

ś

liwego,  w  którym  został  wyodr

ę

bniony  obszar  o  obj

ę

to

ś

ci  V  otoczony 

zamkni

ę

t

ą

 powierzchni

ą

 kontroln

ą

 A

Przyrost  masy  na  obszarze  V  w  czasie  dt  mo

ż

e  nast

ą

pi

ć

  na  skutek  dopływu  

z zewn

ą

trz

 

znak minus wynika st

ą

d, 

ż

e przyrost masy ma miejsce wtedy, gdy przepływ odbywa si

ę

 

do 

ś

rodka,  czyli  przeciwnie  do  jednostkowego  normalnego  zewn

ę

trznie  skierowanego 

wektora n. 

Wyra

ż

enie całkowe  jest strumieniem masy. 

 

 

Przyrost masy mo

ż

e, w przypadku nieustalonego ruchu 

ś

ci

ś

liwego, nast

ą

pi

ć

 tylko na 

skutek lokalnego przyrostu g

ę

sto

ś

ci 

 

Z  równo

ś

ci  obydwu  wyra

ż

e

ń

 otrzymuje 

si

ę

 równanie ci

ą

gło

ś

ci w formie całkowej 

∫∫

A

n

dA

v

p

dt

*

*

background image

;

0

)

*

(

)

*

(

)

*

(

=

+

+

+

z

v

p

y

v

p

x

v

p

t

p

z

y

x

∫∫∫

∫∫

=

+

V

A

n

dA

v

p

dV

t

p

0

*

*

∫∫∫

∫∫

=

V

A

n

dA

v

p

dV

v

p

div

)

*

(

)

*

(

r

∫∫∫

=

+

V

dV

v

p

div

t

p

0

)]

*

(

[

r

0

*

=

+

v

div

p

t

p

r

0

)

*

(

=

v

p

div

r

0

=

v

div

r

;

0

=

+

+

z

v

y

v

x

v

z

y

x

                                                            

(a) 

 

Drugi wyraz równania mo

ż

na przekształci

ć

 zgodnie z twierdzeniem Gaussa-Ostrogradskiego. Zgodnie 

z tym twierdzeniem  

 

Po podstawieniu otrzymuje si

ę

 

 

 

(18.1) 

Przy  zało

ż

eniu  ci

ą

gło

ś

ci  v  i  p  (brak  fal  uderzeniowych)  i  z  uwagi  na  dowolno

ść

  obszaru  całkowania 

mo

ż

na zapisa

ć

 to równanie w formie ró

ż

niczkowej, czyli dla elementu dV 

  

 

(b

Równania  (a)  i  (b)  s

ą

  równaniami  ci

ą

gło

ś

ci  przepływu  w  postaci  całkowej  i  ró

ż

niczkowej  w  zapisie 

wektorowym. W zapisie skalarowym równanie (b) przyjmuje posta

ć

 

(c

Która obowi

ą

zuje w kartezja

ń

skim układzie współrz

ę

dnych.  

Równanie (c) mo

ż

na zapisa

ć

 w nast

ę

puj

ą

cy sposób: 

Pierwsze 

cztery 

wyrazy 

stanowi

ą

 

pochodn

ą

 

substancjaln

ą

 

g

ę

sto

ś

ci, 

wyra

ż

enie  

w nawiasie jest diwergencj 

ą

 wektora pr

ę

dko

ś

ci. Zatem równanie (b) przyjmie posta

ć

 

(18.2) 

 

Dla przepływu ustalonego z równania (b) wynika 

                                                                         

(f) 

czyli 

(18.3) 

Dla  cieczy  (p  =  const)  z  równania  (f)  otrzymuje  si

ę

                                     

(18.4)  

czyli w układzie kartezja

ń

skim 

        

(18.5) 

 

;

0

)

*

(

)

*

(

)

*

(

=

+

+

+

z

v

p

y

v

p

x

v

p

t

p

z

y

x

;

0

*

*

*

*

=



+

+

+

+

+

+

z

v

y

v

x

v

p

z

p

v

y

p

v

x

p

v

t

p

z

y

x

z

y

x

0

)

*

(

)

*

(

)

*

(

=

+

+

z

v

p

y

v

p

x

v

p

z

y

x

background image

0

)

*

(

*

=

+

+

v

p

p

v

t

p

r

r

 

Inn

ą

 form

ą

 zapisu wektorowego równania (b) jest 

 

            

(18.6)  

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

19.Zasada zachowania p

ę

du. Równanie Eulera i  Naviera-Stokesa 

 

background image

Równanie zachowania p

ę

du. 

 

Zasad

ę

  zachowania  p

ę

du  zastosujemy  do  obszaru  płynnego  V,  ograniczonego 

powierzchni

ą

 płynn

ą

 A 

Zasad

ę

 t

ę

 sformułujemy nast

ę

puj

ą

co: pochodna p

ę

du płynu zwartego wewn

ą

trz V wzgl

ę

dem 

czasu jest równa sumie sil zewn

ę

trznych działaj

ą

cych na ten obszar. 

Elementarny p

ę

d płynu zwartego w niesko

ń

czenie małej obj

ę

to

ś

ci dV jest równy 

dV

v

*

*

ς

 

P

ę

d  płynu  zwartego  w  całym  obszarze  płynnym  jest  równy  sumie  tych  wszystkich 

p

ę

dów elementarnych, któr

ą

 wyra

ż

amy za pomoc

ą

 całki 

 

(ilustracja do równania zachowania p

ę

du ) 

∫ ∫ ∫

v

dV

v

*

*

ς

 

Siły  zewn

ę

trzne  działaj

ą

ce  na  obszar  V  podzielimy  na  siły  powierzchniowe  i  siły  masowe 

(obj

ę

to

ś

ciowe). 

Siłami  powierzchniowymi  b

ę

dziemy  nazywa

ć

  siły  oddziaływania  płynu  znajduj

ą

cego  si

ę

  na 

zewn

ą

trz obszaru V na płyn znajduj

ą

cy si

ę

 wewn

ą

trz V, przyło

ż

one na powierzchni płynnej A. 

Na  ka

ż

dym  niesko

ń

czenie  małym  elemencie  dA  działa  jednostkowa  siła  powierzchniowa  

wymiarze napr

ęż

enia [

3

m

N

] zwrócona w kierunku obszaru płynnego. Oznaczamy j

ą

 przez pA. 

Linia jej działania tworzy w przypadku poruszaj

ą

cego si

ę

 płynu lepkiego pewien kat ró

ż

ny od 

zera  z  kierunkiem  jednostkowego  wektora  normalnego  n.  Polem  sił  pA,  które  jest  polem 

background image

tensorowym,  zajmiemy  si

ę

  szczegółowo  przy  wyprowadzaniu  równa

ń

  ruchu  płynów  lepkich. 

W  przypadku  płynu  nieruchomego  lub  nielepkiego  jednostkow

ą

  sił

ę

  powierzchniow

ą

 

wyra

ż

amy za pomoc

ą

 ci

ś

nienia i zgodnie z tym ci poprzednio powiedzieli

ś

my  o kierunku siły 

powierzchniowej, oznaczamy - n p. 

Suma  sił  powierzchniowych  działaj

ą

cych  na  cał

ą

  powierzchni

ę

  A  b

ę

dzie  okre

ś

lona 

całk

ą

 

∫ ∫

A

A

dA

*

 

Dla płynu nieruchomego lub płyn

ą

cego, ale nielepkiego, siła powierzchniowa jest równa 

∫ ∫

A

dA

p

n

*

*

 

Siły  masowe  (  obj

ę

to

ś

ciowe  )  s

ą

  to  siły  wywierane  na  ka

ż

dy  element  masy  

q  dV  zawarty  wewn

ą

trz  V  przez  zewn

ę

trzne  pola  sił  ,np.  pole  grawitacyjne, 

elektromagnetyczne 

gdy 

płyn 

jest 

przewodnikiem 

elektryczno

ś

ci, 

od

ś

rodkowe  

i  Coriolisa  -  gdy  płyn  porusza  si

ę

  w  wiruj

ą

cym  układzie  współrz

ę

dnych.  Załó

ż

my, 

ż

e  układ 

współrz

ę

dnych,  w  który  rozpatrujemy  ruch  jest  nieruchomy  lub  porusza  si

ę

  ruchem 

jednostajnym  po  linii  prostej  i  uwzgl

ę

dniamy  tylko  siły  masowe  niezale

ż

nie  od  ruchu  płynu  

w  tym  układzie.  Reprezentantem  pola  sił  masowych  jest  siła  masowa  przypadaj

ą

ca  na 

jednostk

ę

 masy czynnika, czyli mierzymy je w [N/kg] i oznaczamy przez Fm. W przypadku sił 

grawitacyjnych  b

ę

dzie  to  po  prostu  g,  czyli  przyspieszenie  ziemskie.  Suma  sił  masowych 

działaj

ą

cych na obszar płynny wynosi zatem 

∫ ∫ ∫

V

m

dV

F

q

*

*

 

Ostatecznie  zapisujemy  powy

ż

sz

ą

  zasad

ę

  zachowania  p

ę

du  dla  całego  obszaru 

płynnego V, a wi

ę

c w formie całkowej 

∫ ∫

A

A

dA

*

=

∫ ∫

A

dA

p

n

*

*

=

∫ ∫ ∫

V

m

dV

F

q

*

*

 

 

Zasad

ę

 t

ę

 wyrazi

ć

 mo

ż

na równie

ż

 według sformułowania d'Alamberta, które mówi, 

ż

suma sił zewn

ę

trznych i siły bezwładno

ś

ci musi by

ć

 w ka

ż

dej chwili równa zeru. 

Poniewa

ż

 siła bezwładno

ś

ci b

ę

dzie tutaj równa 

∫ ∫ ∫

V

dV

q

dt

dv

*

*

 

To zamiast powy

ż

szego wzoru mo

ż

na zapisa

ć

 

 

Z powy

ż

szych dwóch równa

ń

 wynika, 

ż

e istnieje to

ż

samo

ść

 

background image

∫∫∫

V

dV

p

v

*

*

r

n

r

n

r

dA

p

A

A

*

∫∫

r

dA

p

n

A

*

*

∫∫

r

∫ ∫ ∫

=

V

dV

v

q

dt

d

*

*

∫ ∫ ∫

V

dV

dt

dv

q

*

*

 

Której udowodnienie jest mo

ż

liwe na innej drodze rozumowania. Nale

ż

y podkre

ś

li

ć

 

ż

e wy

ż

ej 

napisana równo

ść

 (to

ż

samo

ść

) obowi

ą

zuje równie

ż

 w przypadku, gdy zamiast v podstawimy 

w niej dowolny inny skalar lub wektor. 

Analogicznie do zasady zachowania p

ę

du mo

ż

emy napisa

ć

 zasad

ę

 zachowania kr

ę

tu. 

∫ ∫ ∫

∫ ∫

∫ ∫ ∫

×

+

×

=

×

V

A

V

m

A

dV

F

r

q

dA

p

r

dV

q

v

r

dt

d

)

(

)

(

*

)

(

 

 

Zasada  zachowania p

ę

du. 

 

Zasad

ę

 zachowania p

ę

du zastosujemy do obszaru płynnego V, ograniczonego powierzchni

ą

 

płynn

ą

 A. Zasad

ę

 t

ę

 sformułujemy nast

ę

puj

ą

co: pochodna p

ę

du płynu zawartego wewn

ą

trz F 

wzgl

ę

dem czasu jest równa sumie sił zewn

ę

trznych działaj

ą

cych na ten obszar. 

Elementarny  p

ę

d  płynu  zawartego  w  niesko

ń

czenie  małej  obj

ę

to

ś

ci  dVjest  równy  

v  *  p  *  dV  .  P

ę

d  płynu  zawartego  w  całym  obszarze  płynnym  jest  równy  sumie  tych 

wszystkich p

ę

dów elementarnych, któr

ą

 wyra

ż

amy za pomoc

ą

 całki 

                            

(19.1) 

 

Siły  zewn

ę

trzne  działaj

ą

ce  na  obszar  V  podzielimy  na  siły  powierzchniowe  i  siły 

masowe(obj

ę

to

ś

ciowe). 

Siłami  powierzchniowymi  b

ę

dziemy  nazywa

ć

  siły  oddziaływania  płynu  znajduj

ą

cego  si

ę

 

wewn

ą

trz  V,  przyło

ż

one  na  powierzchni  płynnej  A.  Na  ka

ż

dym  niesko

ń

czenie  małym 

elemencie  dA  działa  jednostkowa  siła  powierzchniowa  zwrócona  w  kierunku  obszaru 

płynnego. Oznaczymy j

ą

 przez p

A

. Linia jej działania tworzy w przypadku poruszaj

ą

cego si

ę

 

płynu 

lepkiego 

pewien 

k

ą

ż

ny 

od 

zera  

z kierunkiem jednostkowego wektora normalnego  

przypadku  płynu  nieruchomego,  lub  nielepkiego  jednostkow

ą

  sił

ę

  powierzchniow

ą

 

wyra

ż

amy  za  pomoc

ą

  ci

ś

nienia  i  zgodnie  z  tym 

oznaczamy -  .p. 

Suma sił powierzchniowych działaj

ą

cych na cał

ą

 powierzchni

ę

 A b

ę

dzie okre

ś

lana całk

ą

.                                                                                   

 

Dla płynu nieruchomego, lub 

płyn

ą

cego, ale nie lepkiego, siła 

powierzchniowa jest równa  

 

 

background image

∫∫∫

V

dV

m

F

p

*

*

r

dV

p

m

F

dA

p

dV

p

v

dt

d

V

A

V

A

*

*

*

*

*

∫∫∫

∫∫

∫∫∫

+

=

r

r

r

∫∫∫

V

dV

p

dt

v

d

;

*

*

r

dV

p

m

F

dA

p

dV

p

dt

v

d

V

A

V

A

*

*

*

*

*

∫∫∫

∫∫

∫∫∫

+

=

r

r

r

∫∫∫

∫∫∫

=

V

V

dv

dt

v

d

p

dV

v

p

dt

d

*

*

*

*

r

r

(19.3) 

Siły  masowe(obj

ę

to

ś

ciowe)  s

ą

  to  siły  wywierane  na  ka

ż

dy  element  masy  p  *dV  zawarty 

wewn

ą

trz  V  przez  zewn

ę

trzne  pola  sił,  np.  pole  grawitacyjne,  elektromagnetyczne-gdy  płyn 

jest przewodnikiem elektryczno

ś

ci, od

ś

rodkowe i Coriolisa-gdy płyn porusza si

ę

 w wiruj

ą

cym 

układzie współrz

ę

dnych. Załó

ż

my, 

ż

e układ współrz

ę

dnych, w którym rozpatrujemy ruch jest 

nieruchomy,  lub  porusza  si

ę

  ruchem  jednostajnym  po  linii  prostej  i  uwzgl

ę

dniamy  tylko  siły 

masowe 

niezale

ż

ne 

od 

ruchu 

płynu  

w  tym  układzie.  Reprezentantem  pola  sił  masowych  jest  siła  masowa  przypadaj

ą

ca  na 

jednostk

ę

 masy czynnika, czyli mierzymy jaw [N/kg] i oznaczamy przez F

m

. W przypadku sił 

grawitacyjnych  b

ę

dzie  to  po  prostu  g,  czyli  przy

ś

pieszenie  ziemskie.  Suma  sił  masowych 

działaj

ą

cych na obszar płynny wynosi zatem 

                (19.4) 

 

 

 

Ostatecznie  zapisujemy  zasad

ę

  zachowania  p

ę

du  dla  całego  obszaru  płynnego  V  w  formie 

całkowej 

  

(19.5) 

Zasad

ę

 t

ę

 mo

ż

na wyrazi

ć

 równie

ż

 według sformułowania d'Alemberta, które mówi, 

ż

e suma 

sił zewn

ę

trznych i sił bezwładno

ś

ci musi by

ć

 w ka

ż

dej chwili równa zero. 

Poniewa

ż

 siła bezwładno

ś

ci b

ę

dzie tutaj równa

 

         

 

(19.6) 

to zamiast (19.5) mo

ż

na zapisa

ć

 

 

(19.7) 

Z równa

ń

 (19.5) i (19.6) wynika, 

ż

e istnieje to

ż

samo

ść

 

 

(19.8) 

 

To

ż

samo

ść

  (19.8)  obowi

ą

zuje  równie

ż

  w  przypadku,  gdy  zamiast  v  podstawimy  w  niej 

dowolny inny skalar, lub wektor. 

 

2.Równanie Eulera. 

background image

0

*

*

*

*

*

*

=

+

∫∫∫

∫∫∫

∫∫

dV

p

m

F

dA

p

n

dV

p

dt

v

d

V

V

A

r

r

r

∫∫

∫∫∫

=

A

V

dV

gradp

dA

p

n

*

*

*

r

0

*

*

*

=





+

∫∫∫

V

dV

m

F

p

gradp

dt

v

d

p

r

r

0

)

*

(

=

+

v

p

div

dt

v

d

r

r

0

)

*

(

=

+

v

p

div

dt

dp

r

m

F

gradp

dt

v

d

r

r

=

+

ρ

1

 

Zasad

ę

 zachowania p

ę

du w postaci całkowej podaje równanie (19.7). Dla płynu nielepkiego, 

gdy podstawimy p

A

 =-n* p, otrzymamy 

 

 

(19.9) 

Stosuj

ą

c  twierdzenie  Gaussa-Ostrogradzkiego  do  przekształcenia  całki  wyra

ż

aj

ą

cej 

siły powierzchniowe na całk

ę

 obj

ę

to

ś

ciow

ą

, czyli 

 

 

(19.10) 

otrzymujemy 

 

(19.11) 

Ze  wzgl

ę

du  na  dowolno

ść

  w  wyborze  obszaru  całkowania  V,  dla  ka

ż

dego  elementu 

dV funkcja podcałkowa musi si

ę

 zerowa

ć

, czyli dla przepływu o ci

ą

głych zmianach v, p, p 

b

ę

dzie

  

                                                                          

(19.12)  

Równanie (19.12) nazwane równaniem Eulerajest bilansem sił bezwładno

ś

ci, ci

ś

nienia 

i sił masowych dla niesko

ń

czenie małego elementu płynu o obj

ę

to

ś

ci dV podczas ruchu tego 

płynu. 

Równanie  to  ł

ą

cznie  z  równaniem  ci

ą

gło

ś

ci  w  formie  ró

ż

niczkowej  stanowi  zamkni

ę

ty 

układ  równa

ń

  wystarczaj

ą

cy  do  rozwi

ą

zywania  zada

ń

  dotycz

ą

cych  przepływów  cieczy 

nielepkiej. Mamy wtedy cztery niewiadome 

ν

x

ν

y

,

ν

z

, p i cztery równania skalarne. 

Przepływy gazu doskonałego opisujemy za pomoc

ą

 układu

 

  

(19.13)  

 

 

(19.14) 

 

 

 

 

 

 

 

3.Równania Naviera-Stokesa. 

 

background image

∫∫

A

A

dA

*

r

z

z

y

y

x

x

A

n

p

n

p

n

p

p

*

*

*

r

r

r

v

+

+

=

z

p

y

p

x

p

F

p

dt

v

d

p

z

y

x

m

+

+

+

=

r

r

r

r

r

*

*

xz

xy

xx

x

k

j

p

i

p

τ

τ

*

*

*

r

r

r

r

+

+

=

yz

yy

yx

y

k

p

j

i

p

τ

τ

*

*

*

r

r

r

r

+

+

=

zz

zy

zx

z

p

k

j

i

p

*

*

*

r

r

r

r

+

+

=

τ

τ

z

y

x

p

X

p

dt

dv

p

xz

xy

xx

x

+

+

+

=

τ

τ

*

*

z

p

y

x

Z

p

dt

dv

p

zz

zy

zx

z

+

+

+

=

τ

τ

*

*

z

y

p

x

Y

p

dt

dv

p

yz

yy

yx

y

+

+

+

=

τ

τ

*

*

Równanie  zachowania  p

ę

du  płynu  lepkiego  w  postaci  całkowej  dla  sko

ń

czonego  obszaru 

płynnego  V  podałem  jako  równanie  (19.7).  W  przypadku  obszaru  V,  w  którym  wyst

ę

puj

ą

 

ci

ą

głe  zmiany  parametrów  ,  p,pitd.,  równanie  zachowania  p

ę

du  mo

ż

na  wyrazi

ć

  dla 

elementów  obj

ę

to

ś

ci  dV,  czyli  w  formie  ró

ż

niczkowej.  Przej

ś

cie  z  formy  (19.7)  do 

ż

niczkowej  wymaga  przedstawienia  wszystkich  wyst

ę

puj

ą

cych  w  (19.7)  sił  jako  sił 

obj

ę

to

ś

ciowych. W równaniu (19.7) sił

ą

 powierzchniow

ą

 jest  

 

 

Korzystaj

ą

c z zale

ż

no

ś

ci  znanych

    

 

i analizy wektorowej wzorów na przekształcenie całek powierzchniowych na obj

ę

to

ś

ciowe,

 

mianowicie 

mo

ż

emy (19.7) zapisa

ć

 w nast

ę

puj

ą

cej formie

 

 

Z powy

ż

szego równania wida

ć

ż

e wyra

ż

enie typu        jest  

 

charakterystyczne dla siły powierzchniowej, przypadaj

ą

cej na jednostk

ę

 obj

ę

to

ś

ci. 

Ze wzgl

ę

du na dowolno

ść

 obszaru V musi by

ć

 spełniony warunek zerowania si

ę

 funkcji 

podcałkowej, czyli

 

 

 

 

Jest  to  najogólniejsze  równanie  zachowania  p

ę

du  w  formie  ró

ż

niczkowej.  Zwane  ono  jest 

równie

ż

  równaniem  p

ę

du  w  napr

ęż

eniach.  Mo

ż

na  je  przepisa

ć

  jako  trzy  równania  skalarne 

wykorzystuj

ą

c zwi

ą

zki 

 

otrzymujemy 

 

 

 

 

(19.15) 

(19.16) 

(

)

dV

z

p

y

p

x

p

dA

n

p

n

p

n

p

V

z

y

x

A

z

z

y

y

X

X

∫∫∫

∫∫



+

+

=

+

+

r

r

r

r

r

r

*

*

*

0

*

*

*

=



+

+

∫∫∫

dV

z

p

y

p

x

p

m

F

p

dt

v

d

p

V

z

y

x

r

r

r

r

r

background image

dz

z

V

y

z

x

y

y

v

dx

dv

Θ

=

+

x

p

x

r

0

=

=

=

=

y

x

Y

Y

X

Z

ν

ν

ν

ν

Aby  uzale

ż

ni

ć

  wyst

ę

puj

ą

ce  w  tym  równaniu  składowe  napr

ęż

e

ń

  od  parametrów  przepływu  musimy 

skorzysta

ć

  z  hipotez  dotycz

ą

cych  zwi

ą

zku  mi

ę

dzy  napr

ęż

eniami  i  odkształceniami  elementu  płynu, 

które  s

ą

  wyra

ż

one  poprzez  pola  pr

ę

dko

ś

ci.  Hipotezy  te  pochodz

ą

  z  obserwacji  empirycznych. 

Skorzystamy  z  hipotezy  Newtona  mówi

ą

cej, 

ż

e  napr

ęż

enia  w  poruszaj

ą

cym  si

ę

  płynie  s

ą

 

proporcjonalne  do  pr

ę

dko

ś

ci  deformacji  elementu  płynu.  Gdy 

zastosujemy  prosty  newtonowski 

model  przepływu  mi

ę

dzy  dwoma  warstwami  cieczy 

η

  *  i      porównamy  go  z  ogólniejszym 

odkształceniem  postaciowym  elementu  płynu,  to  zauwa

ż

ymy,  model  Newtona  opisuje  szczególny 

przypadek odkształcenia. 

 

 

                        

                       

 

 

 

                                                                                                                Vy  + 

 

            dz                                                                                V

z

=0 

                                            d

ββββ

 

                                                                                          d

αααα

=0 

                                                                 V

y

                                                 

 

 

W modelu Newtona    

i wszystkie wyrazy tensora pr

ę

dko

ś

ci 

deformacji  elementu 

płynu sprowadzaj

ą

 si

ę

 zgodnie ze wzorami 

 

 

                                                                          

                          

do członów                                                                                                                                                            

                                              

X

Y

Z

z

v

y

v

dt

d

d

Θ

=





+

=

+

2

1

2

1

2

β

α

Y

Z

X

x

v

z

v

Θ

=

+

background image

z

v

y

+

2

1

z

v

y

yz

=

*

η

τ

,

*

2

1

*

1

z

v

k

z

yz

=

τ

[ ]

,

2

1

2

1

2

1

2

1

2

1

2

1

z

x

y

x

y

z

y

z

z

d

T

ε

ε

ε

Θ

Θ

Θ

Θ

Θ

Θ

=

[ ]

,

zz

yz

xz

zy

yy

xy

zx

yx

xx

p

p

p

τ

τ

τ

τ

τ

τ

=

Π

[ ]

Π

=

*

n

p

A

r

r

 

 

                            

 

(19.17) 

Hipoteza Newtona mówi w tym przypadku, 

ż

e napr

ęż

enie 

τ

yz 

wynosi 

 
 

 

(19.18) 

Je

ś

li przyjmiemy, 

ż

e napr

ęż

enia zawsze proporcjonalne do pr

ę

dko

ś

ci deformacji, czyli w tym 

przypadku na podstawie (19.17) 

 

 

(19.19) 

gdzie  k

1

-współczynnik  proporcjonalno

ś

ci  mi

ę

dzy  tensorem  deformacji  i  napr

ęż

e

ń

Porównuj

ą

c (18) i (19) widzimy, 

ż

e: 

 

k

1  =

  2

η

                                                              (19.20)  

Je

ś

li  przeniesiemy  ten  rezultat  na  ogólne  przypadki  płynu,  opisanych  członami 

Θ

  tensora 

deformacji 

 

 

 

to porównuj

ą

c tensory 

 

 

 

 

 

 

 

background image

[ ]

,

2

1

2

1

2

1

2

1

2

1

2

1

z

x

y

x

y

x

y

z

x

x

y

x

d

y

v

y

v

x

v

x

v

z

v

y

v

x

v

x

v

T

ε

Θ

Θ

Θ





+

+





+

=

Θ

 

 

 

 

 

background image

 

 

,

*





+

=

y

v

x

v

z

y

yz

η

τ

,

*

+

=

x

v

z

v

z

x

zx

η

τ

,

*





+

=

y

v

x

v

x

y

xy

η

τ

0

=

v

r

p

n

p

A

r

r

*

=

0

v

r

 
otrzymamy: 
 

 
 
 
 
 
 
 

 

 

Wyja

ś

nimy teraz struktur

ę

 napr

ęż

e

ń

 normalnych p

xx

 , p

yy

 , p

zz

 , wyst

ę

puj

ą

cych 

przek

ą

tnej tensora napr

ęż

e

ń

Przyjmiemy, 

ż

e ci

ś

nienie musi by

ć

 zwarte we wzorach okre

ś

laj

ą

cych p

xx

 , p

yy

 , p

zz

 . 

Gdy 

p

ę

dko

ść

 

czynnika               ,      to  napr

ęż

enia  normalne  s

ą

  sobie  równe  i 

równe ci

ś

nieniu, a wi

ę

                                 p

xx

+p

yy

+p

zz

-3p

                                                      

(19.25) 

Znak minus przy ci

ś

nieniu wi

ąż

e si

ę

 z wyja

ś

nionym ju

ż

 zwi

ą

zkiem: dla v=0 jest 

Jednoczenie  suma  napr

ęż

e

ń

  normalnych  na  podstawie  ogólnych 

twierdze

ń

 rachunku tensorowego jest stała i wynosi -3p równie

ż

, gdy   

 

Pr

ę

dko

ś

ci  wydłu

ż

e

ń

  elementu  płynu  E

x

  ,  E

y

  ,  E

z

  musz

ą

  si

ę

  wi

ą

za

ć

  

z  napr

ęż

eniami  normalnymi.  Dla  cieczy,  gdzie  divv  =0,  napr

ęż

enia  normalne 

mo

ż

emy nast

ę

puj

ą

co uzale

ż

ni

ć

 od pr

ę

dko

ś

ci wydłu

ż

e

ń

(19.26) 

Miar

ą

 dodatkowego odkształcenia obj

ę

to

ś

ci gazów jest ró

ż

na od zera warto

ść

 

divv. Przyjmujemy zatem, 

ż

e dla gazów napr

ęż

enia normalne dodatkowo zale

żą

 od 

(19.24) 

x

v

p

k

p

p

X

X

XX

+

=

=

*

2

*

*

1

η

ε

y

v

p

k

p

p

Y

Y

YY

+

=

=

*

2

*

*

1

η

ε

z

v

p

k

p

p

Z

Z

ZZ

+

=

=

*

2

*

*

1

η

ε

background image

 

 

divv,  

a wi

ę

c: 

 

(19.27) 

Sumuj

ą

c stronami równanie (19.27) z uwzgl

ę

dnieniem (19.25) otrzymamy     

3 p = -3 p + k

1

 

divv + 3k

2

 * divv, 

 

 

st

ą

      (19.28)  

Ostatecznie piszemy zwi

ą

zki napr

ęż

e

ń

 normalnych z polem pr

ę

dko

ś

ci dla gazów: 

 (19.29) 

Dla  cieczy  zale

ż

no

ś

ci  te  sprowadzaj

ą

  si

ę

  do  (19.26).  Do  równania  (19.16) 

podstawimy zwi

ą

zki napr

ęż

e

ń

 z polem pr

ę

dko

ś

ci i ci

ś

nie

ń

. Dla pierwszego z równa

ń

 

(19.16) otrzymujemy w wyniku drobnych przekształce

ń

 przy zało

ż

eniu, 

ż

η

 

= const, 

 (19.30) 

st

ą

d: 

v

div

k

k

p

p

X

XX

r

*

*

2

1

+

+

=

ε

v

div

k

k

p

p

Y

YY

r

*

*

2

1

+

+

=

ε

v

div

k

k

p

p

Z

ZZ

r

*

*

2

1

+

+

=

ε

η

3

2

3

1

2

=

=

k

k

v

div

x

v

p

p

X

XX

r

*

3

2

*

2

η

η

+

=

v

div

y

v

p

p

Y

YY

r

*

3

2

*

2

η

η

+

=

v

div

z

v

p

p

Z

ZZ

r

*

3

2

*

2

η

η

+

=

v

div

x

z

x

v

y

x

v

x

v

z

v

y

v

x

v

x

p

X

p

t

v

p

X

X

X

X

X

X

X

r

*

*

3

2

*

*

*

*

2

2

2

2

2

2

2

2

2

2



+

+

+



+

+

+

=

η

η

η

background image

 

 

p

v

η

=

 (19.31)  

gdzie lepko

ść

 kinematyczna   

 

 

Analogicznie mo

ż

na zapisa

ć

 równanie dla kierunku y i z. 

 

 (19.32) 

 

Układ (19.31) i (19.32) mo

ż

emy zapisa

ć

 w postaci wektorowej 

   (19.33) 

 

Jest  to  równanie  zachowania  p

ę

du  płynu  lepkiego, 

ś

ci

ś

liwego  w  formie 

ż

niczkowej,  tzw.  Równanie  Naviera-Stokesa.  Dla  cieczy  odpada  ostatni  człon 

równania  ze  wzgl

ę

du  na  to,ze  divv  =  O.  Dla  płynu  nielepkiego  v  =  O  i  równanie 

(19.33) przechodzi w równanie Eulera (19.12). 

Równanie  (19.33)  jest  najogólniejszym  równaniem  p

ę

du  jakie  znamy  dla 

nieustalonych  trójwymiarowych  przepływów  lepkiego  płynu  newtonowskiego. 

Fizycznie  uzasadnionym  warunkiem  brzegowym,  który  stosuje  si

ę

  przy 

rozwi

ą

zywaniu  równa

ń

  N  -  S  jest  znikanie  pr

ę

dko

ś

ci  na  nieruchomych 

ś

ciankach 

kanałów  lub  ciał  opływanych.  Gdy 

ś

cianka  porusza  si

ę

,  pr

ę

dko

ść

  płynu  na  niej 

równa si

ę

 pr

ę

dko

ś

ci 

ś

cianki. 

Rozwi

ą

zanie  analityczne  równania  (19.33)  jest  mo

ż

liwe  tylko  dla  bardzo 

uproszczonych  jego  postaci.  Metody  numeryczne  i  wykorzystanie  maszyn 

v

div

x

v

v

v

x

p

p

X

dt

dv

X

X

r

*

*

3

1

*

*

1

2

+

+

+

=

v

div

y

v

v

v

y

p

p

Y

dt

dv

Y

Y

r

*

*

3

1

*

*

*

1

2

+

+

=

v

div

z

v

v

v

z

p

p

Z

dt

dv

Z

Z

r

*

*

3

1

*

*

*

1

2

+

+

=

v

graddiv

v

v

v

gradp

p

m

F

dt

v

d

r

r

r

r

*

3

1

*

*

1

2

+

+

=

background image

 

 

cyfrowych  umo

ż

liwiaj

ą

  rozwi

ą

zanie  tego  równania  w  szerszym  zakresie.  Jednak

ż

nawet  metody  numeryczne  nie  zezwalaj

ą

  obecnie  na  teoretyczne  rozwi

ą

zywanie 

dowolnych  przepływów  lepkich  w  oparciu  o  równanie  N  -  S.  Problemy  te  zostan

ą

 

szczegółowiej  omówione  w  dalszym  ci

ą

gu  tego  rozdziału.  Równanie  N-  S  jest 

podstaw

ą

  w  badaniu  zjawisk  podobie

ń

stwa  przepływu  płynów  bez  potrzeby  ich 

rozwi

ą

zywania. 

 

20. Równanie Daniela Bernoulliego

 

 

Przy nast

ę

puj

ą

cych zało

ż

eniach :  

1. Płyn jest nielepki ( 

, = 0 ) 

2. Przepływ jest stacjonarny 

3. Pole sił masowych ma potencjał II ( f = -grad II )  

4. Płyn jest barotropowy (

=

(p) ) 

Równanie za chowania ilo

ś

ci ruchu mo

ż

na zapisa

ć

 nast

ę

puj

ą

co:  

rotu

u

p

P

u

grad

dt

du

×

=

Π

+

+

+

)

(

2

2

 

Trójmian 

Π

+

+

=

)

(

2

2

p

P

u

E

 

nazywamy trójmianem Bernoulliego. 

Jest pi

ęć

 przypadków stało

ś

ci tego trójmianu: 

 

1.  

 

rotu

u

gradE

×

=

 

0

*

=

=

ds

dE

I

gradE

 

I - wektor jednostkowy 

0

=

ds

dE

 

 

2.  

0

*

=

=

dw

dE

I

gradE

w

 

 Iw,- wektor styczny do linii wirowej 

0

=

dw

dE

 

background image

 

 

 

3. 

u

rotu

λ

=

 

 

u

u

rotu

u

λ

×

=

×

 

st

ą

grad E = 0 

 

4. 

rotu=0  

 

grad E = 0 

5. Dla u=0 

const

p

P

E

=

Π

+

=

)

(

 

Stało

ść

  trójmianu  Bernoulliego  wykazana  dla  pi

ę

ciu  przypadków,  stanowi  podstaw

ę

 

do sformułowania równania Bernoulliego. Ma ono posta

ć

E

const

P

u

=

=

Π

+

+

2

2

 

 

background image
background image

ERROR: syntaxerror
OFFENDING COMMAND: --nostringval--

STACK:

(ZALICZENIE Z MECHANIKI P£YNW _Czarno-Bia‡e_Wojnar Brak 

)

/Title 
()
/Subject 
(D:20041124224136)
/ModDate 
()
/Keywords 
(pdfcreator Version 0.8.0)
/Creator 
(D:20041124224136)
/CreationDate 
(root)
/Author 
-mark-