FIZYKA WODY sciaga, studia, Politechnika Warszawska, I rok mgr Inżynieria Środowiska, Fizyka wody


Założenia ośrodka ciągłości: jego właściwości nie zmieniają się przy przechodzeniu do objętości nieskończenie małych.

Założenie ciągłości odkształceń: dwie cząsteczki, które się stykały w miarę deformacji, obciążeń nadal się stykają, jedynie mogą się zdeformować.

Metoda Lagrange'a jest analizą wędrowną, polega na badaniu właściwości fizycznych( np.: prędkość, gęstość, ciśnienie, temperatura) i opisie poruszającego się płynu. W tej metodzie należy wybrać cząstkę np. poprzez opisanie jej położenia w czasie t0 otrzymujemy wtedy: x0, y0, z0 a następnie śledzi się dana cząstkę w kolejnych chwilach t. Otrzymujemy w ten sposób opis ruchu płynu w postaci funkcji współrzędnych położenia początkowego i czasu:0x01 graphic
= f (x0, y0, z0, t) gdzie 0x01 graphic
jest współrzędną położenia lub własnością płynu w chwili t.. Współrzędne położenia elementu płynu w kolejnych chwilach czasowych, wyrażające parametryczne równanie toru cząsteczki określają równania:

x =x(x0,y0, z0, t); y= y(x0, y0, z0, t); z = z(x0,y0,z0, t)

Składowe prędkości poszczególnych elementów płynu wynoszą:

Vx = ∂x/∂t; Vy = ∂y/∂t; Vz = ∂z/∂t

a składowe przyspieszenia:

ax = ∂2x/∂t2; ay = ∂2y/∂t2; az = ∂2z/∂t2

Inne parametry stanu ruchu płynu np. gęstość i ciśnienia, wyraża się zgodnie ze wzorem:

0x01 graphic
= f (x0, y0, z0, t) a więc:

ρ = ρ (x0, y0, z0, t); p = p (x0, y0, z0, t)

0x08 graphic
Metoda Eulera - analiza lokalna

Opisuje zmiany wielkości fizycznej płynu z poszczególnych punktów przestrzeni, w której odbywa się ruch. Vx = f (R, t)

Vx = fx(x, y, z, , t) Vy = fy(x, y, z, , t)

Vz = fz(x, y, z, , t)

0x01 graphic
0x01 graphic
-przyśpieszenie cząsteczki wyrażone za pomocą pochodnej substancjalnej

Zależność między L i E:

0x01 graphic

Linia prądu - linia pola wektorowego definiowana dla pewniej chwili. W każdym miejscu jest styczna do wektora prędkości w tym miejscu.

Tor cząstki - linia zakreślona przez cząsteczkę płynu w ruchu.

Linia prądu = tor cząstki dla ruchu ustalonego dla zmiennego w czasie jest inaczej.

Smuga - miejsce geom. Zajmowane przez cząstki, które przeszły przez określony punkt w przestrzeni.

Wzory ściśliwości i rozszerzalności.

Ściśliwość - to zdolność do zmiany gęstości płynu w danej temperaturze przy zmianie ciśnienia. Uwzględnia zmianę gęstości płynu przypadającą na jednostkę przyrostu ciśnienia dla stałej temperatury ośrodka: β = (1/ρ)(dρ/dp)

Rozszerzalność - zdolność do zmiany gęstości płynu przy stałym ciśnieniu, temperaturze pod wpływem zmiany temperatury. Wyraża względną zmianę gęstości przypadającą na jednostkę przyrostu temperatury przy danym ciśnieniu:

α = - (1/ρ)(dρ/dT)

Równanie stanu:

dρ /ρ = (1/ρ)(dρ/dp) dp + (1/ρ)(dρ/dT) dt =

βdp - αdT

zmiana funkcji (gęstości) pod wpływem ciśnienia i temperatury

Wzory Helmholtza:

Vx= Vxo +∂Vx/∂x ξ +∂Vx/∂y η+ ∂Vx/∂z ζ

Vy= Vyo +∂Vy/∂x ξ +∂Vy/∂y η+ ∂Vy/∂z ζ

Vz= Vzo +∂Vz/∂x ξ +∂Vz/∂y η+ ∂Vz/∂z ζ

przemieszczenie obrót zmiana czasu

∂Vx/∂z ζ+ ∂Vy/∂z ζ- deformacja

∂Vx/∂z ζ- ∂Vy/∂z ζ- obrót

Co to znaczy Vx/x, Vy/y , Vz/z

założenia: Vz= 0, ∂ /∂z =0, ∂Vy/∂x= 0, ∂Vx/∂y =0

∂Vx/∂x - zmiana prędkości w kierunku ruchu, droga jaką przebyła cząsteczka względem P(0,0)

0x01 graphic

AA'=∂Vx/∂x* ξ*dt +∂, BB'=∂Vy/∂y*η*dt. Cząsteczka nie przemieściła się, lecz zmieniła swoje wymiary, nie zmieniła kształtu

Prędkość wydłużenia względnego

εx = ∂Vx/∂x, εy = ∂Vy/∂y, εz = ∂Vz/∂z

Szybkość wydłużenie względnego = wydłużenie względne/czas

Wydłużenie względne = wydłużenie/ wartość początkowa

Względna prędkość objętości:

∂Vx/∂x ξ+∂Vy/∂y η+∂Vz/∂z ζ

Względna szybkość zmiany przyrostu objętości:

∂Vx/∂x +∂Vy/∂y +∂Vz/∂z =div V(wektor) = F

Jeżeli div V jest ujemne to cząstka zmniejsza swoją objętość

Co to znaczy Vx/y, Vy/x

założenia: Vz=0, ∂ /∂z= 0 , ∂Vx/∂x=0, ∂Vy/∂y=0

0x01 graphic

AA'=∂Vy/∂x *ξ*dt , BB'=∂Vx/∂y*η*dt. Przekątna cząstki zmieniła krawędź, cząstka uległa obróceniu.

dβ=(dα1- dα2)/2=dt(∂Vy/∂x+ ∂Vx/∂y)/2

Prędkość kątowa obrotu cząstki dookoła osi z

ωz=(∂Vy/∂x - ∂Vx/∂y)/2= Ωz*1/2= 1/2*(rot V wektor)z

(rot V wektor)z - wektor rotacji

Ωz - szybkość z jaka przesuwają się boki po przekątnych

Szybkość, z jaka przesuwają się boki do osi cząsteczki

θ z= (dα1+ dα2)/ (2dt) - szybkość kątowa dwóch ścianek względem siebie

1/2θ z= (∂Vy/∂x+ ∂Vx/∂y)/2

θ z+Ωz=∂Vy/∂x-∂Vx/∂y+∂Vy/∂x+∂Vx/∂y

Vy/x=1/2(θ z+Ωz) deformacja kątowa

Vx/y= 1/2(θ z-Ωz) szybkość obrotu

Równanie Eulera:

0x01 graphic

Założenia: płyn jest nie lepki ν = 0

dv/dt - siły bezwładności

Fmj - siły masowe jednostkowe

1/ρ grad p - siły powierzchniowe normalne

Ruch turbulentny - duża wartość Re, energia przekazywana z ruchu podstawowego do pulsacyjnego jest większa niż ciepło wytworzone na skutek tarcia. W efekcie przepływ traci stateczność, zaburzenia wzmagają się i rozprzestrzeniają się.

Wartość przeciętna prędkości:

0x01 graphic

Wartość przeciętna:

0x01 graphic

Dla ustalonych warunków brzegowych pole parametrów jest stałe.

Pulsacja:

0x01 graphic

z def: 0x01 graphic

Przeciętna wartość pulsacji jest równa zero.

Intensywność turbulencji:

0x01 graphic

0x01 graphic

Drgania niewielkie, ale bardzo szybkie.

Turbulencja izotropowa:

0x01 graphic
Niezależna od kierunku.

Równanie ruchu Reynoldsa

Jest to równanie Reynoldsa dla składowej x, dotyczącego przepływów turbulentnych.

Dla ruchu turbulentnego oprócz tarcia wynikającego z lepkości pojawia się tarcie turbulentne:

0x01 graphic

gdzie(-ρ) odpowiada za turbulentne naprężenia styczne.

Ε - współczynnik lepkości burzliwej turbulentnej

Naprężenie to działa na powierzchni xy, więc jest prostopadłe do osi z.

Założenia teorii Prandtla.

Przepływ płynu charakteryzuje się dużymi wartościami liczby Re, siły lepkości nie odgrywają istotnych zmian i można je pominąć. Ze wzrostem Re przepływ przybiera charakter płynu nielepkiego.

Założenia:

- ruch płaski, nie ma prędkości w trzecim kierunku oraz zmian prędkości

- prędkości w kierunku x są dużo większe od prędkości w kierunku y.

- zmiany czegokolwiek w kierunku y są o wiele większe niż w kierunku x.

Teoria ta określa związek między naprężeniami τ turbulentnymi a prędkością uśrednioną.

W strefie laminarnej - τL > τT

W strefie przejściowej - τL = τT

W strefie turbulentnej - τL < τT

0x01 graphic

gdzie Vx i Vy to uśrednione prędkości w kierunku odpowiednich osi

0x01 graphic

Teoria Prandtla mówi o związku pomiędzy naprężeniami turbulentnymi a prędkością uśrednioną

Wzór na przeciętną wartość prędkości:

0x01 graphic

dla teorii Prandtla

0x01 graphic

Naprężenia styczne turbulentne

0x01 graphic

w teorii Prandtla

0x01 graphic

Prawo zachowania masy:

0x01 graphic

1 - zmiana masy w obszarze uwarunkowana zmianą gęstości płynu

2 - ilość masy doprowadzana lub odprowadzana z ośrodka

W objętości kontrolnej, zmiana masy w ośrodku równa jest masie doprowadzanej lub odprowadzanej z obszaru.

Prawo zachowania pędu:

0x01 graphic

1-zmiana pędu

2-impuls sił masowych

3-impuls sił powierzchniowych

Suma wszystkich sił plus siła bezwładności jest równa zeru.

Prawo zachowania energii bilansowe równanie

Energia kinetyczna zmienia się w wyniku pracy wykonanej przez siły masowe i powierzchniowe

0x01 graphic
0x01 graphic

1-wszy człon - zmiana energii kinetycznej w czasie

2-gi człon -ilość ciepła w objętości, zmiana energii wewnętrznej na jednostkę masy

0x01 graphic

strumień ciepła, ciepło

0x01 graphic

-praca sił skalarnych = wykonana praca

Równanie zachowania energii mechanicznej

0x01 graphic

Energia kinetyczna zmienia się w pracy wykonanej przez siły masowe i powierzchniowe

Równanie zachowania(bilansowe) ciepła

0x01 graphic

RÓWNANIE NAVIERA-STOKESA

Równanie analizy ruchu płynu lepkiego w postaci wektorowej:

dV /dt = Fmj - 1/ρ gradp + ν ∇2V + ν/3*grad* divV

1 - siła bezwładności;

2 - działanie sił masowych;

3 - różnica sił normalnych, działanie sił powierzchniowych normalnych

4 - siła tarcia stycznego, szybkość poślizgu;

2,3,4 - efekt sił tarcia;

5- dodatkowe siły tarcia związane ze zmianą objętości cząstki, siły tarcia styczne związane ze zniekształcaniem się płynu ściśliwego

Równanie opisuje przepływ płynów przy Założeniach: płyn ma własność lepkości (płyn Newtowski), pola prędkość są ciągle i różniczkowalne

Przypadki szczególne (modyfikacje N-S)

(krzyżyk przed oznacza wektor):

1)płyn nieściśliwy (*=const => div#υ = 0)

d#υ/dt = #Fmj-1/ρ gradp+mi∇2#V

2)płyn nie lepki (mi=0)

d#υ/dt = #Fmj-1/ρ gradp - równanie Eulera

3)płyn nieruchomy albo w ruchu jednostajnym

(brak ruchu względnego) dv/dt = 0

0 = #Fmj-1/ρ gradp - równanie statyki

Modyfikacja

d#υ/dt = #Fmj-1/ρ gradp N-S mówi o tym ze gdy płyn jest nielekki to stosowane jest równanie zmodyfikowane

Założenia: 1) mi=0 , 2) ruch ustalony , 3) płyn barotropowy ro=f(p) , Fmj = grad - siły masowe powodują potencjał

0x01 graphic



Wyszukiwarka