STATIKA FLUIDA
Ovaj deo mehanike fluida se bavi dejstvom sila na fluid u stanju mirovanja. Budući da se fluid ne kreće, sile (težina fluida, spoljni pritisak) ne vrše rad, već utiču samo na iznos njegove energije. Stoga kažemo da statika fluida proučava energetsko stanje fluida u miru . Ukupna mehanička energija fluida u stanju mirovanja sastoji se od potencijalne energije i energije pritiska, pritiska, koje su povezane izrazom, poznatim kao osnovna jednačina hidrostatike. hidrostatike. Smisao ove fundamentalne relacije može se objasniti razmatranjem sistema prikazanog na slici 3.
Slika 3. Otvoreni sud u kome tečnost miruje nalazi se na određenoj visini iznad podloge 0-0. U tečnosti su označena dva preseka (i (i ii j ), ), paralelna s podlogom. Sve tačke preseka i nalaze i nalaze se na istoj visini (hi) i zato imaju istu potencijalnu energiju – presek i je i je jedna ekvipotencijalna ekvipotencijalna ravan. ravan. Iznad preseka i se nalazi određena količina tečnosti, čija težina, zajedno sa atmosferskim pritiskom (Pa), rezultuje u pritisku (Pi) na ravan i : Pi
= Pa +
Gi S
= Pa +
V i ⋅ γ S
= Pa +
S ⋅ z i ⋅ γ S
= Pa + zi ⋅ γ
gde su Gi i Vi ukupna težina, odnosno zapremina tečnosti iznad ravni i , tečnosti, a S površina poprečnog preseka suda.
(1)
γ je specifična težina
Čestica fluida, mese m, težine G i zapremine V , koja se nalazi u preseku i , ima ukupnu mehani čku energiju: E i
= E potencija ln a + E pritiska = m ⋅ g ⋅ hi + Pi ⋅ V
(2)
STATIKA, KINEMATIKA, DINAMIKA FLUIDA, ISTICANJE FLUIDA
40
Kada ukupnu energiju čestice svedemo na jedinicu njene težine, dobijamo: E i
= hi +
G
Pi
(3)
γ
Jednakim razmatranjem za presek j dobija se: E j G
= h j +
P j
(4)
γ
Kada u izraze (3) i (4) unesemo izraz (1): E i G
= hi +
Pa
+ zi ⋅ γ γ
= hi + zi +
Pa
;
E j
γ
G
= h j +
Pa
+ z j ⋅ γ γ
= h j + z j +
Pa
(5)
γ
oučavamo da su oni me đusobno jednaki: hi
+ zi = h j + z j
(6)
odakle sledi da je ukupna jedinična mehanička energija fluida jednaka za preseke i i j . S obzirom na to da su preseci i i j proizvoljno izabrani, isti zaključak važi i za bilo koji presek u sudu. Prema opštem obliku osnovne jednačine hidrostatike: E G
=h+
P
= const .
(7)
γ
ukupna energija fluida u miru je po jedinici njegove težine konstantna , iako odnos njegove potencijalne energije i njegove energije pritiska zavisi od položaja u odnosu na podlogu. Komentar: Podloga na slici 3 predstavlja tzv. referentnu ravan, tj. ekvipotencijalnu ravan u odnosu na koju merimo potencijalnu energiju ta čaka posmatranog sistema. (U mehanici fluida ekvipotencijalna ravan je uvek paralelna zemljinoj površini.) Apsolutni iznos potencijalne energije tela u gravitacionom polju nas ne zanima jer uvek razmatramo samo razliku u potencijalnoj energiji izme đu dve tačke (dva preseka sistema). Referentna ravan se može postaviti proizvoljno ( čak i iznad suda na slici), ali je najpraktičnije da se postavi u najnižu ta čku sistema. U razmatranom primeru tome bi odgovaralo dno suda. Komentar: Svaki član osnovne jedna čine hidrostatike predstavlja energiju po jedinici težine fluida (J/N ). Pošto džul predstavlja N·m, ona se može izraziti i kao N·m/N . Ako formalno skratimo njutn u brojiocu i imeniocu, osta će samo jedinica dužine – m. Zbog toga u inženjerskoj praksi pomenuti članovi osnovne jedna čine hidrostatike nose naziv visine. Dakle, član potencijalne energije se naziva visina potencijalne energije ili geometrijska visina, ili geodetska visina, a član energije pritiska visina stati čk og pritiska ili jednostavnije visina pritiska. Uprkos tome, nikako ne treba zaboraviti da je pomenuto skraćivanje jedinica samo formalno i da se suštinski uvek radi o energiji po jedinici težine fluida.
Jednu od primena osnovne jedna čine hidrostatike nalazimo kod pijezometra (πιεσις – pritisak; μετρον – mera): prave, vertikalne, na vrhu otvorene staklene cevi pomo ću koje se može meriti statički pritisak u fluidu. U tačkama A i B horizontalnog cevovoda, ispunjenog fluidom (slika 4), postavljeni su pijezometri. Usled pritiska Pa, odnosno Pb u ta čkama A i B, tečnost se u pijezometrima diže do
STATIKA, KINEMATIKA, DINAMIKA FLUIDA, ISTICANJE FLUIDA
41
odgovarajućih nivoa a i b: zbir atmosferskog pritiska i hidrostatičkog pritiska tečnosti u svakom od pijezometara održava ravnotežu odgovarajućem pritisku u osi cevovoda.
Slika 4. Osnovna jednačina hidrostatike važi za svaki od dva pijezometra: P A γ
+ h A =
Pa γ
+ (h A + z A ) ;
P B γ
+ h B =
Pa γ
+ (h B + z B )
(8)
odakle se dobija: z A
=
P A
− Pa γ
=
P A γ
−
Pa
;
γ
z B
=
P B
− Pa γ
=
P B γ
−
Pa
(9)
γ
Ukoliko nas zanima razlika pritisaka u tačkama A i B, ona se može dobiti direktno iz razlike visina nivoa u pijezometrima:
Δ z = z A − z B =
P A γ
−
P B
(10)
γ
Razlika nivoa tečnosti u pijezometrima brojčano je jednaka razlici visina pritiska u ta čkama A i B, što i predstavlja izvesno opravdanje za upotrebu termina "visina". Kao najjeftiniji ure đaj za merenje pritiska, pijezometar se najčešće koristi za praćenje njegove promene duž nekog hidrauličkog sistema. Komentar: Pijezometrom se ne može meriti pritisak gasovitih fluida, što je očigledno iz same njegove konstrukcije. Komentar: Ista razlika nivoa u pijezometrima ne zna či automatski i istu razliku pritisaka – ova druga zavisi i od specifične težine fluida. Važi i obrnuto: pri proticanju razli čitih fluida cevovodom ista razlika pritisaka će dati drugačiju razliku nivoa u pijezometarskim cevima. Komentar: Kada se radi o protočnim, osnosimetričnim sistemima, kao što je cevovod sa slike 4, sve visine (pijezometarska visina, visina u odnosu na referentnu ravan) se ra čunaju od težišta preseka sistema; u prikazanom primeru to je osa cevovoda.
42
STATIKA, KINEMATIKA, DINAMIKA FLUIDA, ISTICANJE FLUIDA
KINEMATIKA FLUIDA
Ova oblast mehanike fluida proučava brzinsko polje u fluidu ne zalaze ći u uzroke njegovog kretanja, odnosno u sile koje pri tome dejstvuju. Najznačajniji njeni rezultati sadržani su u jednačini kontinuiteta, koja predstavlja zakon o održanju mase primenjen na mehaniku fluida. Prema slici 5, na kojoj je prikazan segment horizontalnog cevovoda, jednačina kontinuiteta tvrdi da su maseni protoci fluida u presecima 1 i 2, međusobno jednaki: mτ 1
=m 2
(11)
τ
Drugim rečima, pri proticanju fluida od preseka 1 do preseka 2 masa fluida ne može nestati neznano kud, niti se ni iz čega može stvoriti, što i jeste tvrdnja zakona o održanju mase. Kada se radi o nestišljivim fluidima, izraz (11) važi uvek. On važi i za stacionarno strujanje stišljivih fluida. Međutim, ako stišljivi fluid struji nestacionarno, što zna či da mu se protok s vremenom menja, ta promena će se uvek pre osetiti u preseku 1 nego u preseku 2. U ovom slučaju se jednačina kontinuiteta potvr đuje tek ako se izjednače prosečni maseni protoci, mereni u dužem intervalu vremena.
Slika 5 Izraz (11) se može dodatno razviti ako fluid struji pod izotermskim uslovima, tj. ako su temperature u presecima 1 i 2 jednake. U tom slučaju je i gustina fluida duž cevovoda konstantna ( ρ 1 = ρ 2 ), pa iz: ρ 1
⋅V 1 = ρ 2 ⋅ V 2 τ
τ
dobijamo
V τ 1
= V 2 τ
(12)
Dakle, za izotermsko strujanje fluida, pored konstantnosti masenog protoka, važi i konstantnost zapreminskog protoka duž celog cevovoda.
STATIKA, KINEMATIKA, DINAMIKA FLUIDA, ISTICANJE FLUIDA
Uzevši u obzir da je V τ
43
= w ⋅ S , gde je w - brzina fluida, a S - površina poprečnog preseka
cevovoda, prema slici 5 sledi: w1 ⋅ S1
= w2 ⋅ S 2
w1
i
w2
=
S2 S1
(13)
Brzina strujanja u pojedinim presecima sistema obrnuto je proporcionalna površini tih preseka . Ako je presek cevovoda kružnog oblika, kao na slici 5, izraz (13) se može i dalje razviti: w1 w2
⎛ d ⎞ = ⎜⎜ 2 ⎟⎟ ⎝ d 1 ⎠
2
(14)
gde su d1 i d2 - prečnici odgovarajućih preseka. Izrazi (13) i (14) imaju veoma široku primenu u mehanici fluida uprkos ograničenjima u pogledu stacionarnosti i izotermnosti strujanja, odnosno stišljivosti fluida.
44
STATIKA, KINEMATIKA, DINAMIKA FLUIDA, ISTICANJE FLUIDA
DINAMIKA FLUIDA
Dok je statika fluida proučavala energetsko stanje fluida u mirovanju, dinamika fluida se bavi energetskim stanjem fluida u kretanju. Svako telo koje se kreće pored drugih vidova energije sadrži i odre đeni iznos kinetičke energije. Isto tako, i čestica fluida mase - m, težine - G i zapremine - V , koja se nalazi na visini - h iznad referentne ravni, ukoliko je izložena pritisku - p i kreće se brzinom - w , ima ukupnu mehani čku energiju: E = E potencija ln a
+ E pritiska + E kineticka
E = m ⋅ g ⋅ h + p ⋅ V +
m⋅w
(15)
2
2
Svođenjem ukupne energije čestice na jedinicu njene težine, dolazimo do Bernulijeve jednačine: E G
= h+
p
+
γ
w2 2g
= const .
(16)
koja tvrdi da je ukupna mehanička energija fluida u kretanju konstantna i da se sastoji od potencijalne energije, energije pritiska i kinetičke energije. Bernulijeva jednačina, kako se vidi, predstavlja proširenje osnovne jednačine hidrostatike, a obe su izraz zakona o održanju energije, primenjenog na mehaniku fluida. Slično slučaju kada je fluid mirovao, sve njegove čestice i pri kretanju održavaju istu ukupnu energiju, dok se menja samo odnos izme đu pojedinih njenih vidova. Fluid može sticati brzinu (kinetičku energiju) na račun padanja s neke visine (smanjenja potencijalne energije) ili usled dejstva pritiska (pada energije pritiska). Na primer, u sistemu sa slike 6 na narednoj strani, ukupna energija fluida u presecima 1 i 2 je jednaka: h1 +
p1 γ
+
w1
2
2g
= h2 +
p2 γ
+
w2
2
2g
(17)
ali se već pomoću jednačine kontinuiteta može zaključiti da su njegove brzine, znači i kinetičke energije fluida u dva preseka, različite.
STATIKA, KINEMATIKA, DINAMIKA FLUIDA, ISTICANJE FLUIDA
45
Slika 6 Pošto je cevovod horizontalan (potencijalna energija fluida je u oba preseka ista), gornji izraz se svodi na: p1 γ
+
w1
2
2g
=
p2 γ
+
w2
2
2g
i dalje
p1 γ
−
p2 γ
=
w2
2
2g
−
w1
2
2g
odakle se vidi da je protičući od preseka 1 do 2 fluid zadobio određeni iznos kinetičke energije, ali na račun pada energije pritiska, tako da je ukupna energija ostala konstantna. Ovaj primer daje povoda da se uoči važno pravilo: na mestima u strujnom sistemu gde fluid poveć ava brzinu, dolazi do pada pritiska, i obratno! Iako to nije posebno naglašeno, gornje razmatranje se odnosilo na idealni fluid, dakle takav koji struji bez trenja. Kod realnih fluida, međutim, jedan deo energije se pri strujanju troši na savlađivanje sila trenja i prelazi u toplotu, tako da je ukupan iznos mehani čke energije u preseku 2 manji nego u preseku 1. Isto tako, između dva preseka može postojati crpka koja će fluidu stalno predavati određeni iznos energije u jedinici vremena. Kada se sve uzme u obzir, Bernulijeva jednačina dobija oblik:
⎛ p1 w1 2 ⎞ ⎛ p2 w2 2 ⎞ ⎜ h1 + + ⎟ + E − f 1−2 = ⎜ h2 + + ⎟ ⎜ ⎟ ⎜ ⎟ γ γ 2 2 g g ⎝ ⎠ ⎝ ⎠
(17)
koji se naziva Energijska jednačina. Izraz (17) je najopštiji oblik energetskog bilansa realnog fluida pri strujanju i stoga ima najširu primenu. Komentar: Kao kod osnovne jedna čine hidrostatike, tako i u Bernulijevoj, odnosno u Energijskoj jedna čini svaki član predstavlja energiju po jedinici težine fluida. Zbog formalnog skra ćivanja jedinica ovi članovi u inženjerskoj praksi nose naziv "visina". Uporedo s visinom pritiska i visinom potencijalne 2 energije koje su pomenute u STATICI FLUIDA, sada za član kinetičke energije (w /2g) imamo naziv visina brzine ili visina dinami čk og pritiska, za energiju crpke (E) – visina dizanja crpke, a za član energije gubitaka (f 1-2) – visina gubitaka. Recimo uzgred da visina dizanja crpke i fizi čki predstavlja visinu do koje crpka može da potisne te čnost, a da se visina gubitaka u odre đenim slučajevima može neposredno predstaviti razlikom nivoa u pijezometarskim cevima postavljenim u dva uo čena preseka.
46
STATIKA, KINEMATIKA, DINAMIKA FLUIDA, ISTICANJE FLUIDA
ISTICANJE FLUIDA
Fluide se unutar pogona može transportovati na razli čite načine, pa ipak, ako to dopu[ta tehnološki proces, transport pod dejstvom gravitacije je sigurno najekonomičniji. Na primer, vrlo je pogodno ako je u mlekari sabirni rezervoar za mleko tako postavljen da se iz njega mleko može slobodnim isticanjem distriburati u sve faze obrade. Zbog toga se problematika isticanja tečnosti iz rezervoara obrađuje kao posebna operacija. Osnovna pitanja koja se ovde postavljaju odnose se na činioce koji određuju brzinu, odnosno protok tečnosti pri isticanju, kao i na vreme pražnjenja rezervoara. Proces isticanja može biti stacionaran i nestacionaran. Prvi slučaj je jednostavniji i stoga pogodan za razumevanje principa isticanja, dok drugi više odgovara realnim uslovima.
STACIONARNO ISTICANJE Uslovi stacionarnosti isticanja biće zadovoljeni ukoliko se parametri koji odre đuju brzinu isticanja ne menjaju tokom procesa, što zna či da će i brzina isticanja biti konstantna. Takvi uslovi se ostvaruju samo ako se u rezervoar (slika 7) stalno doliva te čnost u takvoj količini da njen gornji nivo ostane neizmenjen. Opisana situacija je u praksi veoma retka, ali se mnoge realne situacije mogu približno svesti na nju. Na primer, ako se iz rezervoara u kratkom vremenskom intervalu istače samo srazmerno mala koli čina tečnosti, može se smatrati da to praktično neće uticati na opadanje nivoa u njemu. Ta pretpostavka je još opravdanija ako je površina popre čnog preseka rezervoara velika.
Slika 7 Kod isticanja pod dejstvom zemljine teže dolazi do stalnog pretvaranja potencijalne energije u energiju pritiska i energije pritiska u kinetičku energiju. Stoga se za opisivanje isticanja može iskoristiti energijska jednačina. Prema slici 7, presek 1 je odabran u ravni gornjeg nivoa tečnosti, a presek 2 u ravni bočnog otvora pri dnu rezervoara kroz koji tečnost ističe, jer je pritisak (Pa) u ta dva preseka poznat. Referentna ravan je postavljena tako da potencijalna energija preseka 2 bude jednaka nuli, što uprošćava razmatranje. Shodno ovome, za realni fluid važi: h1 +
P1 γ
+
w1
2
2g
− f 1−2 = h2 +
P2 γ
+
w2
2
2g
STATIKA, KINEMATIKA, DINAMIKA FLUIDA, ISTICANJE FLUIDA
Prema slici 7 je h1 = H ; h2
=0 i
P1
= P2 = Pa , pa sledi: Pa
H +
+
w1
γ
što sređivanjem daje:
47
w2
2
2
2g
− w1 2
2g
− f 1−2 = 0 +
Pa γ
+
w2
2
2g
= H − f 1−2
(18)
Brojilac na levoj strani predstavlja razliku kvadrata brzina fluida u otvoru i u ravni nivoa. Pošto su prema jednačini kontinuiteta brzine fluida obrnuto proporcionalne površinama poprečnih preseka kroz koje fluid protiče, a poprečni presek rezervoara (1) je mnogo veći od poprečnog preseka otvora (2), onda važi: S1 >> S1 i w1 << w2 , pa je w2
2
− w12
2g
≈
w2
2
(19)
2g
Na primer, ako je prečnik rezervoara 1 m (100 cm), a pre čnik otvora 5 cm, što odgovara nekoj realnoj situaciji, onda je: 2
w1
2
5 ⎞ 2 2 = ⎛ ⎜ ⎟ ⋅ w2 = 0,0025 ⋅ w2 ⎝ 100 ⎠
To znači da zanemarivanjem brzine fluida u ravni nivoa činimo grešku od 0,25%, što je u tehničkim proračunima dopušteno. Uvođenjem aproksimacije (19) u izraz (18) i njegovim naknadnim preuređivanjem dobijamo brzinu isticanja realnog (r ) fluida kroz otvor:
(w2 )r =
2 g ( H − f 1−2 )
(20)
Uz isti postupak, za idealni fluid (id ) bi se dobio izraz:
(w2 )id =
(21)
2 gH
Član energije gubitaka sada ne postoji jer je idealni fluid neviskozan – struji bez trenja.
Poredeći izraze (20) i (21) uočavamo da pri istom nivou u rezervoaru realni fluid postiže manju brzinu u otvoru jer se deo potencijalne energije tečnosti utrošio na savlađivanje trenja. Odnos brzine isticanja postignute u realnim uslovima u odnosu na idealne naziva se brzinskim koeficijentom:
ϕ
=
(w2 )r = (w2 )id
1−
f 1−2 H
<1
(22)
Njega nije jednostavno izračunati u svakom konkretnom slu čaju, već se koriste neke usvojene prosečne vrednosti, dobijene na osnovu brojnih laboratorijskih eksperimenata. Tako za ve ćinu
48
STATIKA, KINEMATIKA, DINAMIKA FLUIDA, ISTICANJE FLUIDA
tečnosti po osobinama bliskih vodi, dakle i za vodene rastvore, srednja vrednost brzinskog koeficijenta iznosi 0,97, što znači da se pri isticanju prosečno oko 3% energije utroši na trenje. Posle ovoga, brzina isticanja realnog fluida prema izrazu (22) postaje:
(w2 )r = ϕ ⋅ (w2 )id = ϕ ⋅
2 gH ≈ 0,97 2 gH
(23)
Slika 8 Ako ovako izračunatu brzinu pomnožimo površinom popre čnog preseka otvora, dobi ćemo protok znatno veći od protoka koji bi bio izmeren u praksi. uzrok ovome je pojava kontrakcije mlaza. Na slici 8, koja predstavlja okolinu otvora rezervoara sa slike 7, vidi se da tečnost u otvor ulazi konvergentno. Zbog sopstvene inercije, čestice tečnosti slede konvergentne putanje i iza otvora, što ima za posledicu suženje popre čnog preseka mlaza na veličinu Sk. Brzina w2 (potpuna konverzija potencijalne energije u kinetičku) zapravo se javlja u ovom, najužem preseku, tako da se tačna vrednost protoka dobija iz: V τ
= w2 ⋅ S k
(24)
Položaj i veličinu preseka Sk je, nažalost, teško odrediti u pojedinačnoj konkretnoj situaciji, pa se kao i kod brzinskog koeficijenta koristi prose čna vrednost tzv. koeficijenta kontrakcije mlaza: ψ
=
S k S2
= 0,65
(25)
Kombinovanjem izraza (23), (24) i (25), dobija se: V τ
= (w2 )r ⋅ S k = ϕ ⋅
2 gH ⋅ψ ⋅ S 2
= μ ⋅ S 2 ⋅
2 gH
(26)
Proizvod brzinskog koeficijenta (ϕ) i koeficijenta kontrakcije mlaza (ψ) se naziva koeficijentom isticanja (μ); njegova vrednost je za tečnosti bliske vodi 0,62. Izraz (26) više ne sadrži nijednu veli činu koja opisuje osobine tečnosti, već samo geometrijske parametre: visinu nivoa tečnosti i površinu poprečnog preseka otvora. To navodi na zaključak da će, bez obzira na osobine, sve te čnosti pod istim uslovima isticati istim protokom. Ispravnije im protokom. Kod tečnosti koje se znatno bi bilo reći: sve tečnosti bliske vodi isticaće sli čn
STATIKA, KINEMATIKA, DINAMIKA FLUIDA, ISTICANJE FLUIDA
49
razlikuju od vode (naročito po viskozitetu), i koeficijent isticanja će se znatno razlikovati od navedene prosečne vrednosti, pa će i ostvareni protoci biti drugačiji. Komentar: Energijsku, odnosno Bernulijevu jedna činu u praksi treba posmatrati kao izraz koji povezuje energiju fluida (potencijalnu, kineti čki energiju i energiju pritiska) u dva preseka sistema. Preseke treba tako birati da se na osnovnu poznatih parametara može izra čunati onaj nepoznati. Komentar: Izrazi (23) i (26) imaju tipičan oblik, koji se veoma često sreće u tehnici. Kao što je ve ć ranije naglašeno, teorijsko razmatranje nekog fenomena se obi čno prvo izvodi na jednostavnijem modelu (u ovom slu čaju, na modelu idealnog fluida), pa se potom dobijeni rezultat prilagođ ava realnoj situaciji. Kod isticanja se brzina u otvoru najpre izra čunava za idealni fluid, a zatim se uvo đenjem korekcionih koeficijenata (brzinskog, kontrakcije mlaza i isticanja), dobijenih na osnovu obimnog eksperimentalnog materijala, prilagođava realnom stanju.
NESTACIONARNO ISTICANJE U primeru, prikazanom na slici 8, stacionarnost uslova je održavana stalnim nadokna đivanjem istekle tečnosti. U realnim situacijama (slika 9) to nije slučaj, pa će se njen nivo smanjivati, što će posledično izazvati i stalno smanjenje protoka kroz otvor.
Slika 9 Kada je rezervoar iz koga tečnost ističe cilindričnog oblika (S1=const.), metodama diferencijalnog i integralnog rač una se može dobiti vreme potrebno da se nivo tečnosti spusti sa zadate početne vrednosti hp na zadatu krajnju vrednost hk: τ
2 ⋅ S1
= μ
⋅ S2 ⋅
2g
(
h p
−
hk
)
(31)
I za druge oblike rezervoara kod kojih se presek S 1 menja tokom vremena, na primer, za sferni ili rezervoar u obliku položenog cilindra, dobijaju se slični, ali složeniji izrazi. Ne zalazeći u detalje metoda diferencijalnog i integralnog računa pomoću kojih se takva rešenja dobijaju, istaknimo da je na osnovu izraza (31) mogu će predvideti vreme istakanja određene zapremine tečnosti iz rezervoara sa slike 9, što je s gledišta prakse od jedinog značaja.