Mechanika PÅ‚ynów Dr Tomasz Wajman Zespół Maszyn Wodnych i Mechaniki PÅ‚ynów Instytut Maszyn PrzepÅ‚ywowych PA Pokój 110, tel. 42 631 23 60, (24 54) E-mail: tomasz.wajman@p.lodz.pl Literatura " Zbyszko Kazimierski, Podstawy mechaniki pÅ‚ynów i metod komputerowej symulacji przepÅ‚ywów , Aódz, 2004. " Zbyszko Kazimierski, ZdzisÅ‚aw Orzechowski, Laboratorium z mechaniki pÅ‚ynów , Aódz, 2001. " ZdzisÅ‚aw Orzechowski, PaweÅ‚ Wewiórski, Ćwiczenia " ZdzisÅ‚aw Orzechowski, PaweÅ‚ Wewiórski, Ćwiczenia audytoryjne z mechaniki pÅ‚ynów , Aódz, 1999. " Zbyszko Kazimierski, ZdzisÅ‚aw Orzechowski, Mechanika pÅ‚ynów , Aódz, 1986. PÅ‚yny PÅ‚yny ciaÅ‚a posiadajÄ…ce zdolność do zmiany swych ksztaÅ‚tów pod dziaÅ‚aniem dowolnie maÅ‚ych siÅ‚ zewnÄ™trznych, jeÅ›li dziaÅ‚ajÄ… one dostatecznie dÅ‚ugo. PÅ‚yny - różniÄ… siÄ™ od ciaÅ‚ staÅ‚ych tym, że nie posiadajÄ… tzw. PÅ‚yny - różniÄ… siÄ™ od ciaÅ‚ staÅ‚ych tym, że nie posiadajÄ… tzw. sztywnoÅ›ci postaciowej. Wynika to z ich struktury czÄ…steczkowej i maÅ‚ych wewnÄ™trznych siÅ‚ spójnoÅ›ci w porównaniu z ciaÅ‚ami staÅ‚ymi. PÅ‚yny dzielimy na ciecze i gazy. PÅ‚yny ciecze i gazy Ciecze - pÅ‚yny, które zmieniajÄ… nieznacznie swojÄ… objÄ™tość pod dziaÅ‚aniem siÅ‚ zewnÄ™trznych (nawet bardzo dużych) - ciecze sÄ… nieÅ›ciÅ›liwe. WypeÅ‚niajÄ… część zbiorników, równÄ… objÄ™toÅ›ci wprowadzonej cieczy, tworzÄ… powierzchnie swobodne, oraz krople w wyniku dziaÅ‚ania napiÄ™cia powierzchniowego. dziaÅ‚ania napiÄ™cia powierzchniowego. Gazy - pÅ‚yny, które zmieniajÄ… Å‚atwo swÄ… objÄ™tość pod dziaÅ‚aniem zewnÄ™trznych siÅ‚ duża Å›ciÅ›liwość. WypeÅ‚niajÄ… caÅ‚e objÄ™toÅ›ci naczyÅ„, w których siÄ™ znajdujÄ…. Procesy sprężania i rozprężania gazów podlegajÄ… znanym przemianom termodynamicznym. Mechanika PÅ‚ynów " Zjawiska wystÄ™pujÄ…ce podczas ruchu i spoczynku pÅ‚ynów wraz z oddziaÅ‚ywaniem pÅ‚ynów na Å›cianki ciaÅ‚ staÅ‚ych zanurzonych w pÅ‚ynie. " Poznanie praw rzÄ…dzÄ…cych ruchem pÅ‚ynów (szczególnie przepÅ‚ywów turbulentnych). przepÅ‚ywów turbulentnych). Mechanika PÅ‚ynów Statyka Kinematyka Dynamika PÅ‚yn nieruchomy Ruch Ruch "! siÅ‚y wzglÄ™dem naczynia OÅ›rodek ciÄ…gÅ‚y. Element pÅ‚ynu. OÅ›rodek ciÄ…gÅ‚y - pÅ‚yn wypeÅ‚nia przestrzeÅ„ w sposób ciÄ…gÅ‚y, bez pustych obszarów typowych dla struktury molekularnej i zwiÄ…zanych z tym zjawisk mikroskopowych. SformuÅ‚owane zostanÄ… prawa makro-mechaniki w odniesieniu do tzw. elementu pÅ‚ynu. Element pÅ‚ynu dV jest to bardzo maÅ‚a ilość pÅ‚ynu, której wymiary liniowe sÄ… dużo mniejsze maÅ‚a ilość pÅ‚ynu, której wymiary liniowe sÄ… dużo mniejsze od wymiarów opÅ‚ywanych ciaÅ‚ lub kanałów, a jednoczeÅ›nie o wymiarach dużo wiÄ™kszych od Å›redniej drogi swobodnej molekuÅ‚ lub ich temperaturowej amplitudy drgaÅ„. Elementy pÅ‚ynu zawierajÄ… ogromnÄ… liczbÄ™ molekuÅ‚. (1 µm3 powietrza zawiera w warunkach normalnych okoÅ‚o 27mln molekuÅ‚) Zasady zachowania Model oÅ›rodka ciÄ…gÅ‚ego pozwala na okreÅ›lenie makroskopowych wÅ‚asnoÅ›ci pÅ‚ynu jako funkcji przestrzeni oraz czasu i potraktowanie ich jako pól, co stwarza możliwoÅ›ci zastosowania ogólnych twierdzeÅ„ teorii pola w mechanice pÅ‚ynów. Prawa rzÄ…dzÄ…ce zachowaniem siÄ™ pÅ‚ynów zostanÄ… sformuÅ‚owane w oparciu o podstawowe zasady fizyki tj.: " zasadÄ™ zachowania masy, " zasadÄ™ zachowania pÄ™du i momentu pÄ™du (krÄ™tu), " zasadÄ™ zachowania energii. Parametry termodynamiczne i kinematyczne Podstawowymi parametrami termodynamicznymi sÄ…: Á, " gÄ™stość " ciÅ›nienie p, " temperatura bezwzglÄ™dna T, " energia wewnÄ™trzna u, " energia wewnÄ™trzna u, " entalpia h. dla cieczy Á = const. dla gazów idealnych: p = ÁRT, u = cvT, h = cpT r v Podstawowy parametr kinematyczny: prÄ™dkość WÅ‚asnoÅ›ci dyssypatywne pÅ‚ynów Podczas ruchu pÅ‚ynów w ukÅ‚adach adiabatycznych, energia mechaniczna pÅ‚ynu zmniejsza siÄ™ przechodzÄ…c w energiÄ™ cieplnÄ… (w sposób nieodwracalny). Proces taki nazywamy dyssypacjÄ… (rozpraszaniem) energii mechanicznej. mechanicznej. Procesy dyssypacji sÄ… zwiÄ…zane z takimi makroskopowymi wÅ‚asnoÅ›ciami pÅ‚ynu jak: " lepkość, która powoduje zjawiska tarcia (naprężenia styczne), " przewodnoÅ›ciÄ… cieplnÄ…, która powoduje przepÅ‚yw energii z obszarów o wyższej temperaturze do obszarów o temperaturze niższej. Statyka cieczy Statyka cieczy Statyka cieczy zajmuje siÄ™ zagadnieniami równowagi cieczy, a także obliczaniem siÅ‚ wywieranych przez nieruchomÄ… ciecz na Å›ciany zbiorników oraz na ciaÅ‚a nieruchome w niej zanurzone. Równowaga ta może mieć charakter bezwzglÄ™dny lub wzglÄ™dny. W obu przypadkach elementy cieczy nie zmieniajÄ… swego poÅ‚ożenia wzglÄ™dem siebie i Å›cian naczynia, lecz w przypadku równowagi bezwzglÄ™dnej ciecz jest nieruchoma wzglÄ™dem ziemi, a w przypadku równowagi wzglÄ™dnej ciecz wraz z naczyniem znajduje siÄ™ w ruchu. Równanie równowagi Element pÅ‚ynu bÄ™dÄ…cy w spoczynku musi znajdować siÄ™ pod dziaÅ‚aniem caÅ‚kowicie równoważących siÄ™ siÅ‚. Równowaga musi zaistnieć miÄ™dzy siÅ‚ami powierzchniowymi i siÅ‚ami masowymi. dV = dx dy dz r r r r j j j j Y Y Y Y dx Xr Xr r r Rozważamy skÅ‚adowÄ… X k i k i Z Z dy r SiÅ‚y powierzchniowe Fx p = f (x) p "p p + dx "p "x p(x + dx) = p + dx dz "x Równanie równowagi SiÅ‚a masowa siÅ‚a, której wartość jest proporcjonalna do masy ciaÅ‚a na które dziaÅ‚a: Fx = X dm = XÁ dV = XÁ dx dy dz r r Fx = i Fx r r r r j j j j Y Y Y Y dV = dx dy dz dx Xr Xr r r k i k i Z Z dm = Á dV dy r Fx p "p p + dx "x dz Równanie równowagi "p ëÅ‚ p dy dz - p + dxöÅ‚dy dz + Á X dx dy dz = 0 ìÅ‚ ÷Å‚ "x íÅ‚ Å‚Å‚ "p ëÅ‚ öÅ‚ - + Á X dx dy dz = 0 ìÅ‚ ÷Å‚ "x íÅ‚ Å‚Å‚ r r "p "p Á X = Ò! i Á X = i Á X = Ò! i Á X = i Fx = XÁ dx dy dz Fx = XÁ dx dy dz r r r r j j j j "x "x "x "x Y Y Y Y dx Xr Xr r r k i k i r df r r r Z Z "p "p "p dy i + j + k = "p r "x "y "z Fx r r r r p "p Fm = i X + j Y + k Z p + dx "x r dz Á Fm = "p Równanie równowagi r r "p = Á Fm r r r r r r r r "p "p "p "p = i + j + k Fm = i X + j Y + k Z "x "y "z "p "p "p "p "p "p = Á X , = Á Y, = Á Z = Á X , = Á Y, = Á Z "x "y "z "p "p "p dx + dy + dz = Á(Xdx + Ydy + Zdz) "x "y "z dp = Á (X dx + Y dy + Z dz) Różniczka ciÅ›nienia Równanie manometryczne "p "p "p dx + dy + dz = Á(Xdx + Ydy + Zdz) "x "y "z Najczęściej spotykanym polem siÅ‚ masowych w zagadnieniach równowagi cieczy jest pole grawitacji ziemskiej. X X X = 0, Y = 0, Z = g Y Y r pa a g g Z Z Z Z p p "p "p dp "p "p dp p pa = = = Á = = 0, = Á g "x "y dz p = Á g z + C dla z = 0, p = pa Ò! C = pa H pa + Á g H p = pa + Á g z z z Równowaga wzglÄ™dna cieczy Gdy w poruszajÄ…cym siÄ™ naczyniu z cieczÄ… elementy cieczy sÄ… nieruchome wzglÄ™dem siebie i Å›cianek naczynia możemy stosować równania równowagi cieczy uzupeÅ‚niajÄ…c siÅ‚y masowe o nowy rodzaj tych siÅ‚ zwiÄ…zany z ruchem naczynia. Wykorzystamy to w dwóch przypadkach: " postÄ™powego, jednostajnie przyspieszonego ruchu naczynia, " ruchu obrotowego. PostÄ™powy ruch naczynia Wyznaczenie wartoÅ›ci dp = Á (X dx + Y dy + Z dz) ciÅ›nienia X = -a, Y = 0, Z = -g Z Z Y Y dp = -Á (a dx + g dz) z r X X g r r p pa a a p = -Á (a x + g z) + C p = -Á (a x + g z) + C r - a dla x = 0, z = z1, p = pa z1 r C = pa + Á g z1 g x p = pa + Á [g (z1 - z)- ax] PostÄ™powy ruch naczynia Wyznaczenie równania dp = Á (X dx + Y dy + Z dz) powierzchni swobodnej powierzchnia swobodna Ò! dp = 0 X dx + Y dy + Z dz = 0 Z Z Y Y X = -a, Y = 0, Z = -g z r X X g r r p pa a a a dx + g dz = 0 a dx + g dz = 0 Ä… r ax + gz = C - a z1 r g dla x = 0, z = z1 Ò! C = gz1 x Ä… a a z = z1 - x tgÄ… = g g Ruch obrotowy naczynia wirówki Wyznaczenie ciÅ›nienia p dp = Á (X dr + Z dz) pa 2 z X = rÉ , Z = -g r 2 dp = Á (É r dr - g dz) z1 2 2 za pa ëÅ‚ öÅ‚ É r ìÅ‚ ÷Å‚ p = Á - g z + C ìÅ‚ ÷Å‚ a 2 íÅ‚ Å‚Å‚ dla r = 0, z = za, p = pa Ò! 0 Ò! C = pa + Ágza r p F m -g 2 2 îÅ‚ r É p = pa + Á + g (za - z)Å‚Å‚ ïÅ‚ śł 2 ðÅ‚ ûÅ‚ Ruch obrotowy naczynia wirówki Wyznaczenie równania powierzchni swobodnej pa z 2 dp = Á (É r dr - g dz)= 0 r z1 2 r É dr - g dz = 0 za 2 2 2 2 r É r É a - g z = C 2 dla r = 0, z = za Ò! C = -Ágza 0 r 2 2 F m -g r É z = za + 2g Ruch obrotowy naczynia wirówki Wyznaczenie równania powierzchni swobodnej 2 2 r É z z = za + pa 2g zR ObjÄ™tość pod paraboloidÄ… obrotowÄ… jest o poÅ‚owÄ™ mniejsza od z1 objÄ™toÅ›ci walca cylindrycznego, w który ta paraboloida jest wpisana za za 1 1 ( ) ( ) Ä„ R2 (zR - za ) = Ä„ R2 (z1 - za ) 2 1 (zR + za ) = z1 2 2 0 za = z1 - R2 É 4g r R 2 ëÅ‚ öÅ‚ É R2 2 ìÅ‚ - ÷Å‚ z = z1 + ìÅ‚r ÷Å‚ 2g 2 íÅ‚ Å‚Å‚