J. Szantyr WykÅ‚ad nr 25 PrzepÅ‚ywy w przewodach zamkniÄ™tych I Przewód zamkniÄ™ty kanaÅ‚ o dowolnym ksztaÅ‚cie przekroju poprzecznego, ograniczonym liniÄ… zamkniÄ™tÄ…, caÅ‚kowicie wypeÅ‚niony pÅ‚ynem (bez swobodnej powierzchni) Do opisu ruchu pÅ‚ynu w kanaÅ‚ach zamkniÄ™tych stosuje siÄ™ uproszczony model przepÅ‚ywu jednowymiarowego. ZakÅ‚ada siÄ™, że oÅ› kanaÅ‚u jest prawie prosta, a przepÅ‚yw przez przekrój S odbywa siÄ™ z prÄ™dkoÅ›ciÄ… ~ reprezentatywnÄ… , czyli jakÄ…Å› prÄ™dkoÅ›ciÄ… Å›redniÄ… u Najprostszy przypadek: przewód o staÅ‚ym przekroju koÅ‚owym uÅ‚ożony poziomo. PrzepÅ‚yw stacjonarny pÅ‚ynu nieÅ›ciÅ›liwego. Równanie zachowania masy (m-masowe natężenie przepÅ‚ywu): " m ~ ~ ~ (ÁuS)= 0 ÁuS = m = const u = = const "l ÁS D p ~ ~~ "~ ~ Równanie zachowania pÄ™du: (ÁuS)= ÁfS - S - pÄC Dt "l gdzie C obwód przekroju S gdzie C obwód przekroju S pÄ - lepkoÅ›ciowe naprężenia styczne 2 2 "p "p C C 0 = - S - pÄC = - pÄ pÄ dl +"dp = -+" "l "l S S 1 1 Przy staÅ‚ych naprężeniach wzdÅ‚uż kanaÅ‚u o przekroju koÅ‚owym mamy: 4l "p = p1 - p2 = pÄ d Na skutek dziaÅ‚ania siÅ‚ lepkoÅ›ci wzdÅ‚uż kanaÅ‚u wystÄ™puje spadek pÄ ciÅ›nienia wprost proporcjonalny do l i oraz odwrotnie proporcjonalny do d. W przypadku w peÅ‚ni rozwiniÄ™tego przepÅ‚ywu laminarnego, czyli po lw H" 0,03Å" ReÅ" d odcinku poczÄ…tkowym można uzyskać analityczne rozwiÄ…zanie równania Naviera-Stokesa, które prowadzi do wzorów: prÄ™dkość lokalna: "p ( ) ( ) u(r)= (r02 - r2) 0 4µl 4µl prÄ™dkość "p Å" r02 ~ u = Å›rednia: 8µl Wzór na prÄ™dkość Å›redniÄ… można przeksztaÅ‚cić do 2 ~ wzoru Darcy-Weisbacha: l Áu "p = d 2 gdzie współczynnik oporu, lub współczynnik strat liniowych 64 W przepÅ‚ywie laminarnym: = Re W przepÅ‚ywie turbulentnym przez rury hydrodynamicznie 0,3164 = 4 gÅ‚adkie: Re Henri Darcy Julius Weisbach 1803 - 1858 1806 - 1871 W ogólnym przypadku należy r0 öÅ‚ uwzglÄ™dnić wpÅ‚yw = ëÅ‚Re, ìÅ‚ ÷Å‚ k íÅ‚ Å‚Å‚ chropowatoÅ›ci Å›cian: chropowatoÅ›ci Å›cian: Turbulentny profil prÄ™dkoÅ›ci w przewodzie Przypadek trudniejszy: przewód nachylony pod kÄ…tem Ä… Jeżeli zaÅ‚ożymy przepÅ‚yw stacjonarny, to równanie zachowania pÄ™du przyjmie postać: 2 ~ ~ ëÅ‚ öÅ‚ "u 1 "p pÄ C " u p pÄ C ~ ìÅ‚ ÷Å‚ u = fl - - + + gz = - ìÅ‚ ÷Å‚ "l Á "l Á S "l 2 Á Á S íÅ‚ Å‚Å‚ dz gdzie podstawiono skÅ‚adowÄ… siÅ‚y fl = g sinÄ… = -g masowej wzdÅ‚uż l: dl Po scaÅ‚kowaniu pomiÄ™dzy dwoma przekrojami kanaÅ‚u otrzymujemy równanie Bernoulliego dla rzeczywistego przepÅ‚ywu ze stratami: równanie Bernoulliego dla rzeczywistego przepÅ‚ywu ze stratami: 2 ~ ~ ëÅ‚ öÅ‚ ëÅ‚ öÅ‚ u12 p1 u22 p2 pÄ C ìÅ‚ ÷Å‚ ÷Å‚ + + gz1 ÷Å‚ - ìÅ‚ + + gz2 ÷Å‚ = dl +" ìÅ‚ ìÅ‚ 2 Á 2 Á Á S íÅ‚ Å‚Å‚ íÅ‚ Å‚Å‚ 1 albo: 2 ~ ~ u12 p1 u2 p2 + + z1 = + + z2 + hs = H = const 2g Ág 2g Ág hs Daniel Bernoulli gdzie - wysokość strat 1707 - 1783 Wysokość strat możemy podzielić na dwa skÅ‚adniki: -wysokość strat liniowych zwiÄ…zanych z tarciem pÅ‚ynu o Å›cianki przewodu prostoliniowego o staÅ‚ym przekroju, -wysokość strat lokalnych zwiÄ…zanych z obecnoÅ›ciÄ… zaworów, kolan, zwężeÅ„, rozgaÅ‚Ä™zieÅ„ i innych elementów. 2 ~ u l hs = Wysokość strat liniowych wyraża siÄ™ wzorem: 2g d 2 2 ~ ~ ~ ~ u u u12 u22 2 hs =  =  Wysokość strat lokalnych wyraża siÄ™ wzorem: 2g 2g Gdzie Å› jest współczynnikiem strat lokalnych, który może być okreÅ›lony w odniesieniu do prÄ™dkoÅ›ci przed elementem lub prÄ™dkoÅ›ci za elementem. Współczynniki Å› sÄ… okreÅ›lane eksperymentalnie i można je znalezć w odpowiednich tablicach. Poniżej podano kilka przykÅ‚adowych wartoÅ›ci współczynników strat lokalnych. Współczynniki strat lokalnych Rodzaj straty lokalnej Współczynnik straty 2  = 0,5 Wlot ze zbiornika ZaÅ‚amanie przewodu o Ć  = 0,946sin2 Õ 2 + 2,05sin4 Õ 2 2 2 2 ZwiÄ™kszenie przekroju  = (1- A1 A2)  = (A2 A1 -1) Kurek otwarcie 5 stopni Kurek otwarcie 5 stopni  = 0,05  = 0,05  = 31,2 Kurek otwarcie 45 stopni  =10,0 Wlot ssania pompy W przypadku gdy przepÅ‚yw odbywa siÄ™ w przewodach o znacznej Å›rednicy, równanie Bernoulliego powinno być jeszcze uzupeÅ‚nione o współczynnik Coriolisa (lub de Saint-Venanta) Ä… 2 ~ ~ Ä…1 Å"u12 p1 Ä…2 Å"u2 p2 + + z1 = + + z2 + hs = H = const 2g Ág 2g Ág Wynika to z faktu, że energia strumienia niejednorodnego różni siÄ™ od energii obliczonej wedÅ‚ug Å›redniej prÄ™dkoÅ›ci wydatkowej dla tego strumienia. Wobec tego mamy: tego strumienia. Wobec tego mamy: 3 +"u dA A Ä… = 3 ~ u Å" A Współczynnik Ä… jest tym wiÄ™kszy im bardziej niejednorodne jest pole prÄ™dkoÅ›ci przepÅ‚ywu. Gaspard Coriolis Adhemar de Saint Venant 1792 - 1843 1797 - 1886 PrzykÅ‚ad 1: przewodem o zmiennym przekroju przepÅ‚ywa w ciÄ…gu godziny 19600 kg paliwa o gÄ™stoÅ›ci Á=930 kg/m**3 i współczynniku lepkoÅ›ci kinematycznej ½=0,000061 m**2/s. Obliczyć spadek ciÅ›nienia w przewodzie jeżeli wymiary wynoszÄ…: l1 = 5[m],d1 = 50[mm],l2 =10[m],d2 =100[mm],h = 5[m] RozwiÄ…zanie 19600 Q = = 0,00585[m3 s] ObjÄ™toÅ›ciowe natężenie przepÅ‚ywu: 930Å"3600 4Q 4Å"0,00585 c1 = = = 2,98[m s] PrÄ™dkość Å›rednia w części 1: Ä„d12 3,14Å"0,052 4Q 4Å"0,00585 c2 = = = 0,745[m s] PrÄ™dkość Å›rednia w części 2: 2 Ä„d2 3,14Å"0,12 c1d1 2,98Å"0,05 Liczba Reynoldsa w części 1: Re1 = = = 2443 ½ 0,000061 c2d2 0,745Å"0,1 Re2 = = =1221 Liczba Reynoldsa w części 2: ½ 0,000061 Współczynnik strat liniowych w części 1: 0,3164 0,3164 1 = = = 0,045 4 4 Re1 2443 64 64 2 = = = 0,052 Współczynnik strat liniowych w części 2: Re2 1221 Spadek wysokoÅ›ci ciÅ›nienia wywoÅ‚any stratami liniowymi wynosi: 2 2 c1 l1 c2 l2 2,982 5 0,7452 10 h = 1 + 2 = 0,045 + 0,052 = 2,183[m] 2g d1 2g d2 2Å"9,81 0,05 2Å"9,81 0,1 Spadek wysokoÅ›ci ciÅ›nienia wywoÅ‚any stratami w rozszerzeniu przekroju: przekroju: 2 2 2 2 (c - c ) (2,98 - 0,745) (c1 - c2) (2,98 - 0,745) [ ] h = = = 0,255[m] 2g 2Å"9,81 CaÅ‚kowita strata wysokoÅ›ci ciÅ›nienia wynosi: hs =h +h = 2,183+ 0,255 = 2,438[m] Spadek ciÅ›nienia w przewodzie: 2 îÅ‚ - c1 Å‚Å‚ c2 2 "p = ÁgïÅ‚ + hs - hśł = 27245[Pa]H" 0,272[MPa] 2g ðÅ‚ ûÅ‚ PrzykÅ‚ad 2: Syfon o Å›rednicy d i dÅ‚ugoÅ›ci l Å‚Ä…czy dwa zbiorniki, w których powierzchnie cieczy sÄ… oddalone o wysokość h. OkreÅ›lić objÄ™toÅ›ciowe natężenie przepÅ‚ywu wody przez syfon znajÄ…c współczynnik strat liniowych oraz współczynniki strat lokalnych na dopÅ‚ywie i na wypÅ‚ywie. RozwiÄ…zanie Równanie Bernoulliego dla przekrojów 1 i 2: 2 2 c1 c2 + z1 + z2 = + z + "h s 2g 2g Dla przepÅ‚ywu ustalonego mamy: c1 = c2 = c z1 + z2 - z = Co prowadzi do: "h s c2 l ëÅ‚ + Â1 +Â2 öÅ‚ öÅ‚ = +  + ìÅ‚ ÷Å‚ Ponieważ: oraz z1 + z2 - z = h " s "h = c2 ëÅ‚ l 2g d íÅ‚ Å‚Å‚ Otrzymujemy: c2 l ëÅ‚ +Â1 +Â2 öÅ‚ h = ìÅ‚ ÷Å‚ 2g d íÅ‚ Å‚Å‚ 2gh Z tego wyznaczamy Å›redniÄ… c = l prÄ™dkość przepÅ‚ywu: +Â1 + Â2 d NastÄ™pnie obliczamy objÄ™toÅ›ciowe natężenie przepÅ‚ywu przez syfon: 2 Ä„d 2gh Q = l 4 +Â1 + Â2 d PrzykÅ‚ad 3 Z otwartego zbiornika wypÅ‚ywa woda przez przewód o dÅ‚ugoÅ›ci l=200 [m] i Å›rednicy d=100 [mm]. Jaka powinna być wysokość H poziomu wody w zbiorniku aby objÄ™toÅ›ciowe natężenie wypÅ‚ywu z wylotu rurociÄ…gu wynosiÅ‚o Q=40 [l/s]? m m2 Å = 1Å"10 [ ] Å = 1Å"10-6[ ] Dane: h=2 [m] Dane: h=2 [m] s s wlot ze zbiornika  = 0,5 w kolano Âk = 0,2 zawór wylotowy  = 5,0 z RozwiÄ…zanie PrÄ™dkość Å›rednia wypÅ‚ywu z rurociÄ…gu wynosi: 4Q 4Å"0,04 ~ c = = = 5,1[m s] 2 2 Ä„d 3,14Å"(0,1) Liczba Reynoldsa wynosi: ~ cd 5,1Å"0,1 Re = = = 510000 Å 1Å"10-6 PrzepÅ‚yw w rurociÄ…gu jest turbulentny, czyli współczynnik tarcia wynosi: wynosi: 0,3164 0,3164 = = = 0,012 4 4 Re 510000 Wysokość H okreÅ›lamy z równania Bernoulliego: ~ ~ c12 pb c22 pb + + h + H = + + "h s 2g Ág 2g Ág ~ gdzie: c2 = c c1 = 0 AÄ…czna wysokość strat wynosi: 2 ~ ìÅ‚ ÷Å‚ "h = c ëÅ‚Â + 3Å"Âk + + l öÅ‚ s w z 2g d íÅ‚ Å‚Å‚ Czyli niezbÄ™dny poziom wody H w zbiorniku wynosi: 2 ~ c l ëÅ‚1+ Âw + 3Å"Âk +Âw + öÅ‚ H = - h ìÅ‚ ÷Å‚ 2g d 2g d íÅ‚ Å‚Å‚ íÅ‚ Å‚Å‚ Po podstawieniu danych liczbowych: 2 (5,1) ëÅ‚ 200 öÅ‚ H = ìÅ‚ ÷Å‚ - 2 = 39,2[m] 1+ 0,5 + 3Å"0,2 + 5 + 0,012 2Å"9,81íÅ‚ 0,1 Å‚Å‚ Przypadek przewodów o przekroju niekoÅ‚owym lub częściowo wypeÅ‚nionych W przypadku przewodów o przekroju innym niż koÅ‚owy oraz w W przypadku przewodów o przekroju innym niż koÅ‚owy oraz w przypadku przewodów częściowo wypeÅ‚nionych pÅ‚ynem istotnym parametrem jest promieÅ„ hydrauliczny, czyli stosunek pola przekroju strumienia pÅ‚ynu do obwodu zwilżonego: F rh = L z W takich przypadkach liczbÄ™ Reynoldsa obliczamy wedÅ‚ug wzoru: u Å"4rh Re = Å