Post on 13-Oct-2015
description
UNIVERSITATEA DIN BACU FACULTATEA DE INGINERIE
IULIAN FLORESCU
MECANICA FLUIDELOR
NOTE DE CURS PENTRU UZUL STUDENILOR
Editura ALMA MATER Bacu 2007
Tiparul executat sub comanda nr... UNIVERSITATEA din BACU
Str. Spiru Haret nr. 9 Bacu UNIVERSITATEA BACU Aprut n anul 2007
________________________________________________________________________________________________
PREFA Odat cu celelalte discipline tiinifice, mecanica fluidelor s-a dezvoltat rapid n ultimul timp, numeroasele cercetri efectuate lrgind mult cunotinele asupra comportrii fluidelor, ct i a numeroaselor probleme a cror rezolvare depinde de cunoaterea acestora. Paralel a crescut i numrul aplicaiilor n diverse ramuri ale tehnicii moderne, pentru a cror dezvoltare cunoaterea fenomenelor specifice fluidelor a devenit indispensabil. Lucrarea este rezultatul activitii didactice i tiinifice a autorului, profesor doctor inginer n cadrul Catedrei de Energetic, Mecatronic i tiina Calculatoarelor i se bazeaz pe concepia unitar de predare a acestei discipline n toate universitile tehnice din ar. Aceast lucrare ncearc s dea o prezentare a problemelor reprezentative ale disciplinei, precum i modul specific de rezolvare a lor. Lucrarea cuprinde pe ntinderea a 12 capitole probleme ale mecanicii fluidelor i o anex cu aplicaii ale principalelor capitole . Majoritatea capitolelor au un coninut teoretic pronunat cu demonstraii relativ simple i punctate cu exemple tehnice aplicative. Lucrarea se adreseaz n primul rnd studenilor facultilor cu profil mecanic i energetic i are ca scop aprofundarea i consolidarea sub aspect teoretic i aplicativ a cunotinelor legate de echilibrul sau micarea diferitelor tipuri de fluide. Totodat ofer soluii tiinifice pentru alegerea unor subiecte de cercetare aprofundat i este folositoare specialitilor din industriile de profil.
Iulian Florescu
________________________________________________________________________________________________
_______________________________________________________________________
LISTA DE NOTAII
A aria At atmosfera fizic
),,( zyx aaaa vectorul acceleraie a viteza de propagare a loviturii de berbec at atmosfer tehnic b limea deversorului; C coeficientul lui Chzy; viteza absolut n
turbomaini; centru de caren; centrul de presiune
c viteza sunetului; D diametrul E energia total unitar; modulul de
elasticitate (solide) Eu numrul lui Euler F fora
gF fora de greutate
mF fora masic
pF fora de presiune
sF fora de suprafa Fr numrul lui Froude f frecvena
if rezultanta forelor de inerie unitare
mf rezultanta forelor masice exterioare unitare G centrul de greutate g acceleraia gravitaiei H sarcina hidrostatic; energia specific;
energia specific (sarcina) turbomainilor adncimea;
H0 presiunea static Hs nlimea de aspiraie Ht sarcina unui rotor real, fluidul perfect
tH sarcina teoretic a unui rotor ideal cu un numr infinit de pale;
hr energia disipat (pierderea de sarcin) Ix(Iy,Iz) componenta impulsului pe axa Ox (Oy,
Oz); momentul de inerie al suprafeei S fa de Ox (Oy, Oz)
),( kji versorul axei Ox (Oy,Oz) K momentul cinetic k exponentul adiabatic L lungimea; lucrul mecanic l lungimea; lucrul mecanic; lungimea de
amestec (Prandtl); lv lucrul mecanic unitar al forelor de
viscozitate M momentul rezultant Ma numrul lui Mach
0M vectorul moment n raport cu punctul O m masa n vectorul normalei
n exponentul politropic; frecvena; turaia turbomainii
ns rapiditatea n funcie de putere nq rapiditatea n funcie de debit
1n turaia unitar Oxyz triedrul de referin Ox1y1z1 triedrul ataat unui corp n micare P perimetrul udat; p presiunea
np tensiunea unitar de suprafa pa presiunea atmosferic (pat) pcr presiunea critic pd presiunea dinamic pm presiunea manometric pst presiunea static Q debitul volumic; debitul surs punctiform;
debitul QM debitul masic QG debitul gravific
lQ1 debitul unitar q debitul specific; R raza (cilindru, sfer); raza hidraulic;
constanta gazelor perfecte; rezistena la naintare; raza de curbur
R rezultanta forelor exterioare Re numrul lui Reynolds r vectorul de poziie r(,z) coordonata cilindric S suprafaa Sh numrul lui Strouhal s elementul de arc T temperatura absolut; perioada t vectorul tangent la arcul ds t timpul U potenialul forelor masice u(v,w) componenta vitezei pe axa Ox (Oy,Oz)
),( wvu componenta pe Ox (Oy,Oz) a vitezei n micarea medie
ul(vl,wl) pulsaiei componentei vitezei pe Ox(Oy,Oz)
u viteza de transport la turbomaini ),,( wvuV vectorul vitez ),,( wvuV vectorul vitez n micarea medie
),,( llll wvuV vectorul pulsaie al vitezei V viteza medie n seciune v volumul specific (masic) X(Y,Z) coordonata cartezian; componenta forei
masice unitare pe Ox (Oy,Oz); valoarea adevrat a unei msurtori
x(y,z) coordonata cartezian z variabila complex (planul z) W(z.t) potenialul complex coeficientul lui Coriolis;
_______________________________________________________________________
coeficientul de compresibilitate izoterm (modulul de compresibilitate)
circulaia vectorului vitez; intensitatea vrtejului
greutatea specific rugozitatea absolut lungimea caracteristic; grosimea stratului
limit; grosimea peliculei de lubrifiant; grosimea substratului laminar; grosimea (perete);
l grosimea substratului (filmului) laminar modulul de elasticitate (fluide);
coeficientul de viscozitate turbulent (Boussinesq)
coeficientul rezistenei locale viscozitate dinamic; randamentul h randamentul hidraulic v randamentul volumic m randamentul mecanic viteza de deformaie volumic coeficientul lui Darcy; coeficientul
pierderilor de sarcin lineare viscozitate cinematic produsul criterial densitatea tensiunea superficial a lichidului componenta tangenial a tensiunii unitare tensiunea tangenial
),( versorul axei tangeniale n triedrul lui Frenet
0 tensiunea tangenial pe perete funcia de deformaie potenialul vitezelor (planul z); funcia de curent (planul z)
),,( zyx vectorul vrtej viteza unghiular
_____________________________________________________________________________________
CUPRINS
Capitolul I. NOIUNI INTRODUCTIVE 1.1. Generaliti ......................................................................................... ......... 11 1.2. Caracteristicile fizice ale fluidelor ............................................................... 12 1.3. Proprieti fizice fundamentale ale fluidelor................................................ 13 1.4. Proprieti fizice specifice lichidelor .......................................... ................ 16 1.5. Proprieti specifice gazelor ........................................................................ 19 Capitolul II. STATICA FLUIDELOR 2.1. Definiia i obiectul staticii fluidelor .......................................................... 21 2.2. Ecuaiile generale ale hidrostaticii ............................................................. 22 2.3. Legea fundamental a hidrostaticii ...............................................................24 2.4. Ecuaia general a hidrostaticii n cmp gravitaional ..................................25 2.5. Echilibrul relativ al lichidelor .......................................................................26 2.6. Aciunea fluidelor n repaus pe perei solizi ............................................... 27 2.6.1. Aciunea fluidelor n repaus pe suprafee plane ....................................... 27 2.6.2. Aciunea fluidelor n repaus pe suprafee curbe deschise ......................... 30 2.7. Aciunea fluidelor n repaus pe suprafee curbe nchise ............................ . 31 2.7. Plutirea corpurilor ....................................................................................... 34 3. Cinematica fluidelor .................................................................................... 37 3.1. Noiuni specifice ......................................................................................... 37 3.2. Micarea unei particule fluide ..................................................................... 38 3.3. Ecuaia continuitii .................................................................................... 43 4. Dinamica fluidelor ideale ............................................................................. 47 4.1. Ecuaiile lui Euler ........................................................................................ 47 4.2. Relaia lui Bernoulli pentru fluide ideale pe o linie de curent...................... 52 4.3. Relaia lui Bernoulli n micare semipermanent n lungul unei traiectorii ......................................................................... 52 4.4. Relaia lui Bernoulli n micare potenial nepermanent .......................... 54 4.5. Relaia lui Bernoulli pentru cureni cu seciunea finit ............................... 54 4.6. Calculul debitului prin orificii ..................................................................... 56 4.7. Teorema impulsului i teorema momentului cinetic n cazul micrii permanente a fluidelor ....................................................................................... 58 4.8. Teorema impulsului i teorema momentului cinetic aplicate tuburilor de curent n micare permanent ............................................................................. 61 5. Dinamica fluidelor reale ...............................................................................63
_____________________________________________________________________________________ 8
5.1. Micarea laminar a fluidelor reale ............................................................. 63 5.2. Starea de tensiune ntr-un fluid n micare .................................................. 63 5.3. Ecuaiile de micare a fluidelor reale sub forma dat de Cauchy (n componente de tensiuni) ..................................................................................... 65 5.4. Ecuaiile Navier-Stokes pentru micarea laminar a fluidelor reale .................................................................................................................... 66 5.5. Relaia lui Bernoulli pentru o linie de curent, n micarea laminar a fluidelor reale ................................................................................... 70 6. Analiza dimensional i teoria similitudinii ............................................... 73 6.1. Metodele analizei dimensionale .................................................................. 73 6.2. Noiuni despre similitudine ......................................................................... 76 7. Noiuni de teoria hidrodinamic a lubrificaiei ......................................... 78 8. Teoria stratului limit .................................................................................. 85 8.1. Ecuaiile de micare ale stratului limit ...................................................... 86 8.2. Desprindea stratului limit ........................................................................... 89 9. Micarea turbulent a fluidelor reale ......................................................... 91 9.1. Structura micrii turbulente ....................................................................... 91 9.2. Tensiunea tangenial n micarea turbulent ............................................. 92 9.3. Distribuia vitezelor n micarea turbulent ................................................ 94 9.4. Ecuaia Reynolds pentru micarea turbulent a fludelor reale .................... 94 9.5. Calculul pierderilor de sarcin .................................................................... 96 9.6. Conducte netede i conducte rugoase; grosimea stratului laminar ............................................................................................................... 97 9.7. Determinarea coeficientului pierderilor de sarcin liniare .......................... 98 9.8. Calculul pierderilor locale de sarcin .......................................................... 99 10. Curgerea prin orificii i ajutaje ...............................................................101 11.1. Calculul debitului unui orificiu mic, sub sarcin constant .....................101 10.2. Calculul debitului orificiului mare .......................................................... 102 10.3. Calculul debitului orificiului necat ........................................................ 103 10.4. Curgerea sub sarcin variabil, prin orificii situate la baza rezervorului. Timpul de golire al unui rezervor .................................................................... 104 10.5. Curgerea sub sarcin variabil i cu debit afluent constant ..................... 105 10.6. Curgerea sub sarcin variabil, printr-un orificiu necat ......................... 107 10.7. Curgerea prin ajutaje ............................................................................... 108 10.8. Jeturi de fluid .......................................................................................... 111 11. Micri permanente n conducte sub presiune ...................................... 113 11.1. Calculul conductelor compuse n serie ................................................... 114
CUPRINS _____________________________________________________________________________________
9
11.2. Calculul conductelor compuse n paralel ................................................ 115 11.3 Calculul conductelor cu debit continuu i terminal ................................. 116 11.4 Calculul conductelor n sifon ................................................................... 118 12. Micarea nepermanent n conducte sub presiune ............................... 121 12.1. Lovitura de berbec n conducte sub presiune .......................................... 122 12.2. Ecuaiile fenomenului lovitura de berbec ............................................... 123 12.3. Soluiile generale ale ecuaiilor cu derivate pariale ale fenomenului lovitura de berbec ......................................................... 126 Bibliografie ........................................................................................ ............. 129
_____________________________________________________________________________________ 10
Introducere _____________________________________________________________________________________
11
1. INTRODUCERE
1.1. Generaliti
Mecanica fluidelor reprezint o diviziune a Mecanicii teoretice, care studiaz micrile, respectiv repausul fluidelor ideale sau reale, compresibile sau incompresibile, sau interaciunea dintre fluidele n micare sau repaus i corpurile solide cu care acestea vin n contact.
Mecanica fluidelor se mparte n trei pri: statica, cinematica, i dinamica. Statica fluidelor studiaz repausul fluidelor i aciunile exercitate de acestea asupra suprafeelor solide cu care acestea vin n contact. Cinematica fluidelor studiaz micarea fluidelor fr s se in cont de forele care intervin i modific starea de micare. Dinamica fluidelor abordeaz micarea fluidelor considernd forele care intervin i transformrile energetice produse n timpul micrii.
Denumirea de Mecanica fluidelor a aprut relativ recent (n secolul XX) i este atribuit studiului general al micrii i al interaciunii fluidelor cu suprafeele corpurilor solide cu care vin n contact. Iniial cu acest studiu se ocupa Hidraulica, cuvnt care deriv din grecescul hdraulos, provenit din legtura a dou cuvinte: hdor ap i aulos - conduct, i care reflect una din primele probleme practice care a preocupat oamenii.
Aceast tiin a cunoscut o diversificare i dezvoltare n strns legtur cu problemele teoretice privind Aerodinamica (stratul limit, rezistena la naintare, teoria profilurilor aerodinamice), Hidraulica (micarea lichidelor cu suprafa liber, micarea aluviunilor, micarea prin medii poroase) i Dinamica gazelor.
Datorit complexitii fenomenelor aprute n micarea fluidelor reale, a aprut necesitatea experimentrii pe modele n tunele aerodinamice i apoi pe baza teoriei similitudinii s-au extins rezultatele la problemele tehnice care au fost modelate. Rezultatele deosebite acumulate de Mecanica fluidelor au permis ca n ultimul timp s apar noi domenii tehnice precum: transportul pe ap i subacvatic, hidrotransportul, meteorologia , exploatarea modern a zcmintelor etc.
Astzi Mecanica fluidelor este o disciplin mai mult teoretic, care studiaz legile general valabile pentru starea de repaus sau micare a fluidelor. Fenomenele proprii lichidelor, gazelor sau aerului sunt studiate respectiv de hidraulic, termotehnic i aerodinamic sau de alte discipline specifice cum ar fi transferul de cldur, construcii hidrotehnice, construcii aerospaiale .a.
n funcie de condiiile impuse la limit se observ dou feluri de aplicaii ale dinamicii fluidelor i anume:
a) curgerea fluidului n jurul unui corp solid considerat izolat (avioane, automobile, paraute) fenomen la care urmrim n mod preponderent fora necesar naintrii;
b) corpul solid delimitez micarea fluidului (micarea n conducte, canale) la care intereseaz transportul de energie.
Ca o remarc general putem afirma c transportul energiei se realizeaz numai cu ajutorul fluidelor n micare (ap, vapori, petrol, gaz natural, aer comprimat).
Mainile hidraulice i pneumatice, utilizeaz micarea fluidelor n domenii delimitate de frontiere solide, care pot fi i ele n micare. Se deosebesc mainile care
Mecanica fluidelor _____________________________________________________________________________________
12
preiau energia mecanic i o transmit fluidului (pompe, ventilatoare), ct i maini care utilizeaz energia fluidului pentru a crea alte forme de energie (motoare).
S-au diversificat foarte mult domeniile n care se studiaz i utilizeaz micarea fluidelor: hidrotransportul, fenomene de filtrare, poluarea, mecanica suspensiilor, aeroelasticitatea, magnetohidrodinamica i altele.
1.2. Proprieti fizice ale fluidelor
Mecanica teoretic definete dou mari categorii de corpuri materiale: corpuri
solide rigide i solide deformabile i corpuri fluide, care cuprind lichidele i gazele. Fa de cele trei stri de agregare cunoscute: solid, lichid, gazoas, putem aminti i o a patra form: plasma, care se definete ca o stare ale crei proprieti sunt determinate de existena electronilor i ionilor n stare liber.
Fluidele sunt corpuri materiale care se caracterizeaz n primul rnd prin proprietatea de fluiditate. tim c moleculele unui corp au o stare de agitaie nentrerupt, care este funcie de temperatur i c ntre ele exist fore de atracie (de coeziune). La corpurile solide moleculele ocup locuri bine determinate (stabile), n jurul crora execut oscilaii de amplitudine mic funcie de temperatur, aceste corpuri avnd form i volum fix, sub aciunea unor fore exterioare.
Spre deosebire de corpurile solide, lichidele i gazele, sub influena unor fore exterioare relativ mici, pot cpta deformaii orict de mari, astfel nct iau forma recipientului solid n care se gsesc. n consecin lichidele i gazele nu au form proprie i ele se caracterizeaz prin posibilitatea de a deplasa foarte uor particulele din care sunt formate, datorit forelor de coeziune mici.
Fluidele trebuie deosebite de materialele care sunt deformabile (situate ntre fluide i solide) cum ar fi pastele sau metalele topite, de a cror studiu se ocup alt tiin numit reologie.
Lichidele reprezint fluide care sunt practic incompresibile i sub aciunea forelor gravitaionale iau forma vasului n care exist fr a umple acest vas. Spre deosebire de lichide, gazele sunt fluide la care forele de coeziune sunt mult mai mici ca la lichide i care umplu n totalitate recipientul n care se gsesc, oricare ar fi forma i dimensiunea lui.
Mecanica fluidelor studiaz medii continue, omogene i izotrope. Un mediu este continuu i omogen, dac are aceeai densitate n orice punct i este izotrop dac prezint aceleai proprieti n toate direciile. Exist la fluide linii, puncte, sau suprafee de discontinuitate, care prezint condiii specifice la limit.
Ipoteza de continuitate
Omogenitatea i izotropia permit ca proprietile i relaiile stabilite pentru o
particul fluid de dimensiuni mici determinate de condiia neglijrii micrii proprii a moleculelor sau de micarea brownian la gaze s fie valabile pentru tot fluidul. Ipoteza general a continuitii impune pentru mrimile fizice: densitate, vitez, presiune, temperatur, funcii care depind de coordonatele punctului i de timp i care sunt continue cu excepia unor linii, puncte, suprafee de discontinuitate.
n studiul mecanicii fluidelor utilizm diferite modele de fluid, n funcie de ipotezele simplificatoare pentru calcule, cum ar fi: fluid uor (fr greutate), fluid ideal (fr viscozitate), fluid incompresibil, la care volumul unei mase determinante este
Introducere _____________________________________________________________________________________
13
constant, fluid real (compresibil i vscos) , fluide vscoase i incompresibile (lichidele), fluide fr greutate dar compresibile (gazele).
Mecanica fluidelor studiaz fenomenele att cu metode experimentale ct i teoretice, de cele mai multe ori combinndu-le. n studiul teoretic se utilizeaz teoremele generale ale mecanicii (teorema impulsului, teorema momentului cinetic, teorema energiei cinetice, legi de conservare), utiliznd un calcul matematic complex. Metodele experimentale de studiu se aplic pentru verificarea calculelor teoretice, pentru determinarea unor legi generale, determinarea unor corecii utiliznd modele fizice la alte scri, rezultatele extinzndu-se prin similitudine.
1.3. Proprieti fizice fundamentale ale fluidelor a) Densitatea (greutatea specific) Densitatea se definete ca masa unitii de volum:
Vm
Vm
V ddlim
0=
= (1.1) unde: m este masa unui element de volum V.
Admind ipoteza continuitii, densitatea este o funcie continu de coordonatele punctului i de timp: = (x,y,z,t). Densitatea se msoar n kg/m3 i are aceeai valoare n orice punct al fluidului omogen.
Inversul densitii se numete volum specific sau volum masic:
=kgmv
31 (1.2)
Greutatea specific este greutatea unitii de volum:
= 3dd
mN
VG (1.3)
ntre densitate i greutate specific exist relaia: g=
Densitatea variaz funcie de presiune i de temperatur. Pentru lichide variaia n raport cu presiunea poate fi neglijat. Densitatea fluidului scade odat cu creterea temperaturii. Pentru ap densitatea maxim este n jurul valorii de 4oC i are valoarea de 1kg/m3. Variaia densitii apei funcie de temperatur este redat n figura 1.
Fig. 1 Variaia densitii apei n funcie de temperatur
(oC) 0 5 10 15 20 25 30 35 40 45
0,98
0,99
0,97
1
(kg/m3)
Mecanica fluidelor _____________________________________________________________________________________
14
La gaze, variaia densitii i greutii specifice cu temperatura i presiunea este deosebit de mare i este dat de ecuaia de stare.
b) Compresibilitatea izotermic a fluidelor este proprietatea de variaie a densitii
(volumului), sub influena variaiei presiunii. Dac are loc o variaie de presiune p pentru un fluid cu volum V i presiune p, se
produce o variaie relativ de volum VV proporional cu p, dat de relaia:
pVV = (1.4)
unde este coeficient de compresibilitate cubic (m2/N), iar semnul minus arat c unei creteri a presiunii i corespunde o scdere a volumului.
n majoritatea fenomenelor studiate considerm lichidele ca fluide incompresibile. De exemplu apa este de 100 de ori mai compresibil dect oelul. Totui exist fenomene care se studiaz innd cont de compresibilitatea lichidelor, cum ar fi lovitura de berbec sau sonicitatea (teoria sonicitii a fost fondat de Gogu Constantinescu n 1916, un mare savant romn care a trit n Anglia).
Gazele sunt mai compresibile dect lichidele. Se poate neglija compresibilitatea gazelor pentru viteze mai mici de 0,6c (c este viteza sunetului).
Definim modulul de elasticitate cubic ca inversul modulului de compresibilitate cubic:
VPV
dd1 == (1.5)
Relaia poate fi exprimat i funcie de densitatea . Pornind de la condiia ca n procesul de comprimare masa s rmn constant:
V = C sau: dV + V d = 0 Separnd variabilele obinem:
VV dd =
de unde:
d= (1.6)
tiind c viteza de propagatre a sunetului stabilit de Newton este: =c ,
rezult:
p
pc
dd1
dd
== (1.7)
Deoarece dd11 p== i la fluide incompresibile tinde ctre zero, rezult
0dd =
p , deci c , adic o propagare instantanee a sunetului, ceea ce este n
contradicie cu realitatea fizic. Prin urmare, n fenomenele legate de propagarea undelor de presiune n medii fluide este necesar considerarea proprietii de compresibilitate a fluidului.
Introducere _____________________________________________________________________________________
15
Micrile n fluide compresibile pot fi clasificate n funcie de viteza pe care o au fa de viteza sunetului prin numrul lui Mach (Ma), care este adimensional i ne d raportul dintre viteza v i viteza sunetului c (celeritate) n mediul respectiv:
cv=Ma (1.8)
Pentru Ma < 1, micarea este subsonic, iar pentru Ma > 1, micarea este supersonic.
b) Viscozitatea Proprietatea de viscozitate a fost explicat i definit diferit de oamenii de tiin: - Newton a considerat c viscozitatea este o consecin a forelor de coeziune
care reacioneaz la deplasarea relativ a particulelor de fluid. Aceast ipotez nu poate fi valabil pentru gaze, la care distanele intermoleculare sunt mari i forele de coeziune neglijabile.
- Maxwel explica viscozitatea fluidelor, prin capacitatea de a face s apar fore, atunci cnd se produc variaii brute ale formei fluidului.
Putem concluziona c dac fluidul este n micare, n diferite plane de separaie apar fore sau tensiuni tangeniale (fore raportate la aria suprafeei), care se opun variaiei formei volumului considerat, frneaz micarea i modific repartiia vitezelor.
Viscozitatea reprezint mecanismul transmiterii micrii n fluid. Viscozitatea se mai poate defini ca o proprietate comun tuturor fluidelor, prin
care cu fore suficient de mici se pot produce deformaii orict de mari, cu viteze de deformare mici.
Considerm un paralelipiped dreptunghic, cu aria bazei S i nlimea h, figura 2.
Fig. 2. Determinarea forei de viscozitate Suprafaa AIBICIDI alunec cu viteza v fa de baz, fora necesar pentru a
imprima aceast vitez fiind:
hvSF = (1.9)
relaie care a fost determinat experimental i unde este coeficient dependent de natura fluidului i se numete viscozitate dinamic. Stratul aderent la plac are aceeai vitez v cu placa. Atracia dintre acest strat i urmtorul face ca i acesta s fie antrenat cu o vitez mai mic vI , astfel nct diferena creat s produc micarea, .a.m.d.
Al
D
A B
C
Bl
ClDl
All
CllDll
Bll
h
v F
S
Mecanica fluidelor _____________________________________________________________________________________
16
Unitatea de msur pentru coeficientul de viscozitate dinamic este poise
scm
g
n sistemul vechi CGS, sau poisseuille
sm
kg , n sistemul internaional. Tensiunea
tangenial ce apare ntre dou straturi infinit vecine (h dn), (v dv) este:
nv
SF
dd== (1.10)
Aceast tensiune are tendina de a egala vitezele celor dou straturi, deci se opune micrii stratului cu vitez mai mare (are sens opus micrii acestui strat). Fluidele ale cror tensiuni tangeniale de viscozitate n micare laminar sunt date de relaia (1.10), se numesc newtoniene.
Raportul dintre viscozitatea dinamic i densitate se numete viscozitate cinematic:
= (1.11)
Unitatea de msur n sistemul internaional SI, pentru viscozitatea cinematic este m2/s, iar n vechiul sistem CGS este stockes (cm2 /s).
Cteva valori pentru viscozitatea cinematic la temperatura normal pentru diferite fluide sunt prezentate n tabelul 1.
Tabelul 1 Fluidul Ap Benzin Alcool Petrol Ulei Glicerin Aer
m2/s 106 0,6510-6 1,3310-6 2,510-6 1,710-4 8,710-4 7,410-3 Variaia viscozitii cu temperatura este diferit pentru lichide i gaze. La lichide
scade cu creterea temperaturii, iar la gaze crete. 1.4. Proprieti specifice lichidelor Adeziunea este proprietatea ce rezult din atracia dintre moleculele unui fluid i
cele ale suprafeei corpului solid cu care vine n contact. Dac atracia intermolecular a lichidului este mai mic dect cea dintre lichid i perete, atunci lichidul ud peretele (ader la acesta). Ca exemplu putem da apa, care ader la un perete de sticl. n caz contrar spunem c lichidul nu ud pereii (de exemplu mercurul).
La gaze adeziunea este neglijabil. Fig. 3 Suprafa de separaie n echilibru
r2
d1
d2r1
dS1 dS1
dS2
dS2
Introducere _____________________________________________________________________________________
17
Experimental s-a constatat c suprafaa liber a unui lichid tinde s-i micoreze aria.
Aceasta se explic prin faptul c fiecare molecul din suprafaa de separaie este atras de moleculele vecine, iar rezultanta acestor fore este ndreptat spre interior. O mas oarecare de lichid i modific forma astfel ca la un volum dat s fie suprafaa minim (de exemplu: forma sferic a picturilor de ap).Tensiunea superficial , este devinit prin fora ce se exercit tangenial pe unitatea de lungime. Ea modific presiunea n lichid. Pentru a demonstra acest lucru considerm un element din suprafaa de separaie, (un dreptunghi curbiliniu) care are laturile dS1 i dS2 , figura 3.
Se observ c: 1
11
ddrS= i:
2
22 r
dd S=
Proiectnd forele elementare pe normal i tiind c 2
d2
dsin 11 = i
2d
2dsin 22 = , obinem:
=
=
1
12
12
2
21
21
dd2
dd2
dd2
dd2
rSSS
rSSS
(1.12)
Suma acestor dou fore trebuie s echilibreze fora care rezult din diferena de presiune de pe cele dou fee ale suprafeei de separaie p:
+==
1
12
2
2121
ddddddrSS
rSSSSpF
Dup simplificri rezult:
+=
21
11rr
p (1.13) relaie care este cunoscut sub denumirea de formula lui Laplace.
Capilaritatea reprezint consecina tensiunii superficiale pentru tuburi subiri.
Fig. 4 Tub capilar
z
z 2R
r
Mecanica fluidelor _____________________________________________________________________________________
18
Dac introducem n ap un tub de sticl de diametru mic (sub un milimetru) apa ud bine sticla, formnd un menisc concav, iar diferena de presiune este ndreptat n sus, apa urcnd n tub pn ce diferena de nivel z echilibreaz aceast diferen de presiune, figura 4.
( )
+=
2121
11zrr
(1.14) unde: 1 i 2 sunt greutile specifice ale apei, respectiv aerului.
Considernd c suprafaa meniscului are aproximativ forma unei calote sferice, (r1 = r2 = R), relaia 1.14 devine:
( ) ( ) z
Rzr == 2121 cos22 unde R este raza tubului i este unghiul de contact al apei cu tubul. Neglijnd 2 n raport cu 1, putem scrie c:
cos2zR =
iar dac 0, obinem legea lui Jurin ntlnit n fizic:
RgRhz ===
22 (1.15)
Pentru lichide neaderente (mercurul fa de sticl), meniscul este convex iar n tubul capilar se formeaz o denivelare (z < 0).
Studiul fenomenelor capilare prezint importan n studiul fenomenelor de infiltraii, n msurtori efectuate cu aparate ce cuprind tuburi capilare .a.
Absoria i degajarea gazelor Majoritatea gazelor se dizolv (ptrund prin difuziune) n lichide. Absoria este un
proces fizico-chimic i se produce dac concentaia componentelor gazelor este mai mare dect cea corespunztoare lichidului la presiunea i temperatura respectiv. Absorbia crete odat cu creterea presiunii i scade odat cu creterea temperaturii. La presiune i temperatur normal n ap se dizolv 2% gaze n greutate: astfel este posibil viaa florei i faunei acvatice.
Dac componenta fazei absorbante (lichidul) este mai mare dect cea a fazei absorbate (gazul) fenomenul se produce invers i se numete desorbie. Desorbia crete odat cu creterea presiunii.
Dac n anumite poriuni ale unui lichid n micare presiunea scade pn la valoarea presiunii de vaporizare la temperatura dat se produce vaporizarea lichidului nsoit de degajare de gaze dizolvate. Apare fenomenul numit cavitaie, un fenomen care este duntor pentru mainile i instalaiile hidraulice.
Bulele de vapori i de gaz, ajung n zone cu presiuni mai mari i se recondenseaz producnd suprapresiuni i creteri de temperatur mari, precum i zgomote. Aceast ipotez nu a putut explica de ce metalele cu rezisten mecanic mai mic sunt mai rezistente dect metalele cu rezisten mecanic mare (de exemplu: bronzul fa de oel).
O alt ipotez, chimic, explic distrugerea metalelor prin faptul c vaporii i gazele degajate pun n libertate oxigenul atomic, care este foarte activ chimic i corodeaz metalul.
Introducere _____________________________________________________________________________________
19
Ipoteza termochimic pune la baza fenomenului temperaturile foarte mari (mii de grade) care apar.
Exist i o ipotez electric a fenomenului de cavitaie, care se bazeaz pe diferena de potenial dintre bulele de gaz i lichid.
Cavitaia este un fenomen complex, care se aplic prin ansamblul tuturor acestor ipoteze.
Presiunea de vaporizare crete odat cu temperatura unui lichid i este mai mare pentru aceeai temperatur la lichidele volatile, figura 5.
Fig. 5 Variaia presiunii de vaporizare n funcie de temperatur 1.5. Proprieti specifice gazelor
Gazele datorit spaiilor intermoleculare mari sunt fluide mult mai uoare i mai compresibile dect lichidele. Gazele ocup prin expansiune tot volumul disponibil (coeziunea este neglijabil). Datorit compresibilitii accentuate, densitatea variaz mult cu presiunea la o temperatur dat. Ecuaia de stare Clapeyron-Mendeleev se exprim pentru o mas de gaz m dat, de volum V:
p V = mRT sau:
p = RT (1.16) nlocuim pe R cu
obinem:
T=p
unde: este constanta universal a gazelor perfecte
Kmol
J31,8 okk , iar este masa
molar. Gazele care satisfac ecuaia (1.21) se numesc gaze perfecte. Dac se ia n consideraie
expresia densitii mpV
= ecuaia (1.20) se scrie sub forma mpV mRT T= = . (1.17)
n studiul repausului sau al micrii unui gaz perfect (fr frecri sau oc) se deosebesc urmtoarele legi de variaie a densitii n funcie de presiune: izocora =
40
alcool
ap\ t (oC)
Pv
1 at
0,062
0 20 60 80 100
Mecanica fluidelor _____________________________________________________________________________________
20
const., izoterma p = const., adiabata (n care gazele nu schimb cldur cu mediul
exterior) kp = const., unde k este exponentul adiabatic i politropa n
p = const., unde
n este exponentul politropic. Dac ntr-un volum dat se afl un amestec de gaze atunci presiunea parial a
fiecrui gaz este egal cu presiunea pe care ar exercita-o dac ar ocupa singur ntregul volum, iar presiunea total este egal cu suma presiunilor pariale (legea lui Dalton).
Statica fluidelor _____________________________________________________________________________________
21
2. STATICA FLUIDELOR 2.1. Definiia i obiectul staticii fluidelor Statica fluidelor studiaz echilibrul fluidelor i aciunea pe care acestea le exercit
asupra corpurilor solide cu care vin n contact. ntr-un fluid n repaus nu apar fore de viscozitate, ele fiind condiionate de
deplasarea relativ a particulelor. Rezult c n starea de repaus proprietatea de viscozitate nu se manifest i prin urmare relaiile din statica fluidelor sunt stabilite pentru fluidele perfecte (pentru fluide incompresibile i compresibile).
Un fluid n repaus se manifest ca un fluid ideal (nevscos) relaiile fiind valabile i pentru fluide reale (vscoase).
Un fluid n repaus este acionat de dou categorii de fore, care l echilibreaz: forele masice i forele de suprafa.
Forele masice sunt analoage celor ntlnite n mecanica corpurilor rigide i se datoreaz prezenei cmpurilor exterioare. Cele mai obinuite fore masice ntlnite sunt cele de greutate datorate cmpului gravitaional exterior masei de fluid considerate. n cazul unui repaus relativ (fluidul se afl n repaus fa de un sistem de referin mobil care execut o micare accelerat fa de un sistem de referin fix), pe lng forele de greutate apar i forele de inerie.
Forele de suprafa joac rolul forelor de legtur din mecanica rigidului. Pentru un fluid n repaus forele de suprafa sunt fore de presiune fiind compresiuni normale la elementele de suprafa. Pentru a demonstra acest lucru separm cele dou poriuni A i B ale unui fluid aflat n echilibru prin seciunea S figura 6.
Fig. 6 Echilibrul unei particule de fluid Fiecare dintre cele dou poriuni trebuie s fie n echilibru. Considernd echilibrul
poriunii A, putem nlocuim efectul poriunii B printr-o for de suprafa care este o compresiune. O for elementar F este normal la elementul de suprafa S deoarece n cazul n care ar fi nclinat componenta sa tangenial ar strica echilibrul fluidului.
Prin definiie presiunea n punctul M este:
( )SF
SFMp
S ddlim
0== (2.1)
Rezult c presiunea este n funcie de coordonatele punctului, valoarea ei nu depinde de orientarea arbitrar a seciunii S care trece prin punctul M. Presiunea ntr-un punct fiind aceeai pe orice direcie, este o mrime scalar. Pentru a demonstra aceasta detam din fluidul n echilibru tetraedrul OABC aa cum rezult din figura 7.
S M
SFA B
Mecanica fluidelor _____________________________________________________________________________________
22
Fig. 7 Particul de fluid n echilibru Introducem pe feele tetraedrului forele elementare de presiune (compresiuni
normale). Considerm n normala la suprafaa ABC (dirijat spre exteriorul volumului de fluid) i S mrimea acestei suprafee.
Considerm fora masic (elementar) unitar ( )ZYX ,,fm i px, py, pz, pn presiunile n lungul axelor de coordonate i al normalei, ecuaiile de echilibru pe direciile axelor sunt:
( )( )( )
=+
=+
=+
06
cos2
06
cos2
06
cos2
xyzZk,nSpxyp
xyzYj,nSpzxp
xyzXi,nSpyzp
nz
ny
nx
(2.2)
Deoarece:
( ) ( ) ( )2
cos2
cos2
cos xyk,nS,zxj,nS,yzi,nS ===
rezult:333zZpp,yYpp,xXpp nznynx ===
Trecnd la limit, tetraedrul tinznd n toate direciile (dup x,y,z), ctre punctul M obinem:
px = py = pz = pn = p(M) = p(x,y,z) (2.3) Rezultatul este independent de nclinarea feei ABC, deci presiunea ntr-un fluid
n repaus este un scalar. 2.2. Ecuaiile generale ale hidrostaticii
(ecuaiile lui Euler pentru statica fluidelor) Ecuaiile pentru echilibrul fluidelor se obin din anularea rezultantei forelor ce
acioneaz asupra maselor de fluid. Considerm o particul fluid de forma unui paralelipiped desprins dintr-un fluid aflat n echilibru de dimensiuni dx, dy, dz, i densitate , figura 8.
i y
PnS
2zyPx
2xzPy
2yxPz
x
z
A
B
C
jk
Statica fluidelor _____________________________________________________________________________________
23
Rezultanta forelor masice unitare ( ),ZY,X,fm acioneaz n centrul de mas al particulei considerate i are expresia:
zyxffF mmmm ddddd == (2.4) i are componentele:
dFmx = Xdxdydz, dFmy = Ydxdydz, dFmz = Zdxdydz. (2.5) Forele de suprafa sunt fore de presiune datorate aciunii fluidului asupra
particulei considerate. tiind c forele de presiune sunt proporionale cu mrimea suprafeelor elementare considerate i c presiunea este n funcie de coordonatele punctului n spaiu, p = p (x,y,z), atunci forele de presiune pe feele determinate de planele sistemului de referin se pot exprima astfel: pdydz, pdxdz, pdxdy
Fig. 8 Particul de fluid n echilibru La presiunile pe feele opuse se adaug creterile pariale datorate variaiilor
obinute prin deplasarea n cele trei direcii, obinndu-se:
yxzzppz,xy
yppz,yx
xpp ddddddddd
+
+
+
Condiia de echilibru dup axa Ox va fi:
0dddddddd =+
+ zyxXzyx
xppzyp (2.6)
sau:
=+
=+
=+
0dddddd
0dddddd
0dddddd
zyxZzyxzp
zyxYzyxyp
zyxXzyxxp
(2.7)
n mod analog s-a calculat dup Oy i Oz.
Dup simplificare cu dx,dy,dz, ecuaiile difereniale de echilibru dup cele trei axe sunt:
ZzpY
ypX
xp =
==
,, (2.8)
x
zxyypp ddd
+
zyxxpp ddd
+
yxzzpp ddd
+z
O
Ol
C
pdydz
pdxdy pdxdz
yCl Bl
BA
Al
Mecanica fluidelor _____________________________________________________________________________________
24
Cele trei ecuaii exprim condiiile de echilibru ale volumului de fluid considerat, ntre forele de presiune i forele masice. Acestea sunt ecuaiile generale ale hidrostaticii ecuaiile cu derivatele pariale de ordinul I, stabilite de Euler pentru echilibrul fluidului. A gsi condiia de integrabilitate a acestui sistem nseamn a preciza condiiile pe care trebuie s le ndeplineasc fora masic unitar mf (X,Y,Z) pentru ca sub aciunea sa fluidul s rmn n echilibru.
Multiplicm cele trei ecuaii, respectiv prin dx, dy, dz, le adunm i obinem:
( )zZyZxXzzpy
ypx
xp dddddd ++=
+
+ (2.9)
Membrul stng al egalitii reprezint difereniala total a presiunii p, astfel putem scrie:
dp = (Xdx+Ydy+Zdz) (2.10) Ecuaia reprezint expresia variaiei de presiune ntre dou puncte ale unui fluid
situate la distana elementar ds care are proieciile dx, dy, dz. Deoarece membrul stng al egalitii este o diferenial total, expresia are sens
dac paranteza din membrul drept este o diferenial total a unei funcii de x, y ,z n acelai domeniu.
Considerm o funcie scalar U(x,y,z) uniform n domeniul dat. Dac componentele forei masice unitare X, Y, Z pot fi exprimate ca derivate
pariale ale funciei U n raport respectiv cu x, y, z, putem scrie:
zUZ
yUY,
xUX
==
= (2.11) sau: gradUfm = (2.12) unde U reprezint potenialul forelor masice exterioare.
Aceasta reprezint condiia necesar i suficient pentru ca un cmp de fore masice unitare s determine echilibrul unui fluid, adic trebuie s fie potenial. Spre exemplu, cmpul gravitaional este potenial deci sub aciunea sa un fluid poate fi n echilibru.
2.3. Legea fundamental a hidrostaticii
Dac nlocuim n relaia (2.10), expresiile din relaia (2.11), obinem: dp + dU = 0 (2.13)
iar prin integrare: p + U = Ct (2.14) care reprezint legea fundamentaal a hidrostaticii. O remarc general este: dac un fluid se afl n repaus fa de un sistem de referin fix (sau sistem inerial) singurele fore masice care apar n fluid sunt forele de greutate (create de cmpul gravitaional); n cazul unui echilibru relativ apar n plus i forele de inerie.
Consecinele legii fundamentale a hidrostaticii 1) Dac n ecuaia (2.10), considerm dp = 0, (p = ct), se obine Xdx+Ydy+Zdz
= 0 (deoarece 0), care reprezint ecuaia diferenial a suprafeelor de egal presiune, numite i suprafee izobare.
Definim suprafaa echipotenial locul geometric al punctelor n care potenialul forelor masice este constant.
Statica fluidelor _____________________________________________________________________________________
25
Dac fluidul este incompresibil ( = ct), pentru U = Ct se obine p = ct, deci ntr-un fluid n repaus absolut suprafeele echipoteniale sunt i izobare.
2) Din relaia (1.13) rezult c fora masic ce acioneaz asupra unei particule fluide este normal la suprafaa echipotenial (izobar) ce trece prin punctul de aplicaie al forei i este ndreptat n sensul scderii potenialului, deci n sensul creterii presiunii.
3) Suprafeele echipoteniale nu se intersecteaz deoarece presiunea fiind un cmp scalar este unic ntr-un punct dup orice direcie i deci n punctul respectiv nu pot exista dou poteniale egale.
4) O suprafa izobar este i izodens ( = ct); afirmaia este evident pentru un fluid incompresibil la care = ct n toat masa deci i pe suprafaa izobar.
5) Temperatura pe o suprafa echipotenial este constant deoarece p i fiind constante din ecuaia lui Clapeyron-Mendeleev rezult:
ctRpT == . (2.15)
6) Suprafaa de separaie dintre dou fluide imiscibile cu densitile 1 2 este echipotenial; de asemenea i dintre un lichid i un gaz. Considernd c ntre dou puncte infinit vecine ale acestei suprafee diferena de presiune este aceeai, indiferent n ce fluid este calculat: UUp 21 ddd == . Rezult c:
( ) 0d21 = U (2.16) de unde: dU = 0, deci U = ct.
7) Dac forele masice sunt neglijabile n raport cu cele de presiune presiunea n fluid este constant.
n acest caz, dac 0=mf rezult din relaia (2.11) c U = ct, deci conform relaiei (2.14), p = ct (conform primei consecine). Rezult c dac ntr-o zon a fluidului are loc o suprapresiune p aceasta se transmite n toat masa fluidului. Aceast consecin poart numele de principiul lui Pascal.
2.4. Ecuaia general a hidrostaticii n cmp gravitaional Particularizm componentele forei masice unitare pentru cmpul gravitaional X
= 0, Y = 0 i Z = - g. Rezult c: kgfm = (2.17) Integrnd obinem: ( ) +==
)(
dD
CgzzgU (2.18)
Ecuaia fundamental a hidrostaticii devine,CUP =+ CgkP =+ sau
lCzp =+ . Putem scrie relaia sub forma: ctzp =+ (2.19)
n care termenii au dimensiuni de lungime. Considernd un fluid n repaus absolut, figura 9 putem scrie relaia (2.19) pentru
dou puncte ale fluidului A i B:
Mecanica fluidelor _____________________________________________________________________________________
26
BBAA zpzp +=+ tiind c pA= po relaia devine:
pB = p = po + (zA- zB) = po + h (2.20) Se observ c presiunea crete liniar cu adncimea. Unghiul pe care l face
graficul cu verticala este:
==H
Htg (2.21)
deci tg este proporional cu , fiind un coeficient ce rezult n raport de scrile alese.
Fig. 9 Distribuia de presiuni ntr-un lichid n repaus n cazul mai multor lichide imiscibile n repaus ele se aeaz n ordinea
greutilor lor specifice, fig. 10.
Fig. 10 Distribuia de presiuni n cazul a trei lichide imiscibile z = ct rezult p = ct. ntr-un lichid n repaus planele orizontale sunt plane izobare. De aici dou
consecine: - suprafaa liber a lichidelor este orizontal (principiul vaselor comunicante) - suprafaa de separaie dintre dou lichide imiscibile este un plan orizontal. 2.5. Echilibrul relativ al lichidelor Asupra unui lichid aflat n echilibru relativ fa de un sistem de referin ce se
mic accelerat n raport cu un sistem fix, acioneaz fore masice de greutate i fore masice de inerie. Notm forele masice absolute (de greutate) cu ( )ZYXfm ,, i fore de inerie unitare cu ( )iiii ZYXf ,, , rezult:
Ugradff im =+ unde U este potenialul forelor masice, sau:
z = 0
P0 P0 P
P0+ (zA-zB)B
A
H
H
h
zB
zA
h1
h2
h3
1
2
3
1
2
3
P0 P 1 < 2 < 3 P0 + 1 h1
P0 + 1 h1 + 2 h2
P0 + 1 h1 + 2 h2 + 3 h3
Statica fluidelor _____________________________________________________________________________________
27
zUZZ
yUYY
xUXX iii
=+=+
=+ ,, (2.22) Pentru deducerea ecuaiei suprafeei libere unde p = po, din ecuaia de echilibru
CUp =+ rezult: U(x,y,z) = Cl (2.23)
Rezolvarea unei probleme de echilibru relativ necesit parcurgerea urmtoarelor etape:
- se alege un sistem de axe legat de vasul n care se gsete lichidul; - se determin acceleraia (fora de inerie unitar) pentru o particul de fluid, n
funcie de coordonatele sale: ( )zyxafi ,,= i fora de greutate unitar gfm = ;
- se determin potenialul U cu expresia: ( ) ( ) ( ) +++++= zZZyYYxXXU iii ddd ;
- n relaia (2.23), constanta Cl se determin din condiia ca volumul n echilibru absolut s fie acelai cu cel din echilibrul relativ;
- se determin repartiia presiunilor pe pereii vaselor.
2.6. Aciunea fluidelor n repaus pe perei solizi Aciunea fluidelor asupra unui perete solid cu care este n contact se calculeaz
nsumnd forele elementare de presiune: SnpFp dd = (2.24)
unde n este normala la suprafaa elementar dS a peretelui orientat spre fluid. Dac suprafaa este oarecare, forele elementare de presiune sunt oarecare n
spaiu i efectul lor este un torsor format de rezultanta forelor i momentul rezultant n raport cu originea sistemului ales:
=S
p SnpF d (2.25)
=S
o SnxprM d (2.26)
Considernd c se cunoate repartiia de presiuni din starea de repaus a fluidului, vom analiza starea de echilibru absolut.
2.5.1. Aciunea fluidelor n repaus pe suprafee plane Normala n n acest caz este constant i deci forele de presiune sunt paralele i
de acelai sens:
Sp SpnF d (2.27)
Rezultanta este o for normal la suprafa, orientat dinspre fluid spre perete: =
So SprxnM d (2.28)
Mecanica fluidelor _____________________________________________________________________________________
28
Considerm c rezultanta forelor de presiune are punctul de aplicaie n C, numit i centrul de presiune, momentul rezultantei este:
==S
CpClo SprxnFxrM d (2.29)
Conform teoremei lui Varignon pentru sisteme de fore paralele cele dou momente sunt egale: loo MM = (2.30)
rezult:
=S
d
d
Sp
Sprr SC (2.31)
Relaia (2.31) reprezint expresia vectorului de poziie al centrului de presiune.
a. Aciunea unui fluid uor n repaus pe o suprafa plan n acest caz presiunea este constant i rezult:
pSF
pSnSpnF
p
Sp
=
== d (2.32)
G
S
SC rS
Srr ==
d
d (2.33)
Aciunea unui fluid uor n repaus pe o suprafa plan este o for normal la suprafa, de mrime pS, care se aplic n centrul de greutate al suprafeei.
b. Aciunea unui fluid greu n echilibru pe o suprafa plan Considerm un vas cu un perete nclinat i alegem un sistem de axe convenabile
conform figurii 12. n cazul fluidului greu presiunea este p = h
Fig. 12 Aciunea unui lichid pe un perete plan nclinat
===S SS
p SynSynShnF dsindsind (2.34)
C
h
y yl
x
xl
G dS
pF
O
Statica fluidelor _____________________________________________________________________________________
29
unde: S
Sy d este momentul static al ariei S n raport cu axa x i putem scrie:
SySy GS
= d (2.35) deci: ( ) ShnSynF GGp == sin (2.36)
Rezultanta foelor de presiune este normal la suprafa i are valoarea egal cu greutatea unui cilindru de lichid avnd aria seciunii egal cu S i nlimea egal cu adncimea centrului de greutate al suprafeei S.
Poziia centrului de presiune se determin n modul urmtor:
Sy
Syr
Sy
Syr
Sy
Syr
Sh
Shrr
G
S
S
S
S
S
S
SC ====
d
d
d
dsin
dsin
d
d
Rezult coordonatele centrului de presiune C, punctul de aplicaie a forei Fp:
Sy
Sxyx
G
SC
=
d i
Sy
Syy
G
SG
=
d2
(2.37)
unde xyS
ISxy = d este momentul de inerie centrifugal al ariei S i xS
ISy = d2 este momentul de inerie al ariei S fa de axa x.
Rezult: Sy
IySy
Ix
G
xC
G
xyC == (2.38)
Dac alegem convenabil sistemul de axe (axe de simetrie), Ixy = 0 deci i xC = 0. Putem aplica teorema lui Steiner: SyII G
lxx += 2 (2.39)
de unde: Sy
Iyy
G
lx
GC += (2.40) n cazul n care pe suprafaa liber a lichidului se exercit o presiune relativ p1
formulele stabilite sunt valabile dac sistemul de referin se alege cu Ox la nivelul
presiunii atmosferice, deci la nlimea 1
1Ph = .
Dac peretele de suprafa S este orizontal C coincide cu G i presiunea este constant pe toat suprafaa.
n cazul fundului unui rezervor presiunea are mrimea SH, unde S este aria seciunii, H este nlimea lichidului i nu depinde de forma vasului, figura 13.
Fig. 13 Paradox hidrostatic
H S S S S
Mecanica fluidelor _____________________________________________________________________________________
30
Este un paradox hidrostatic deoarece fora de presiune n mod aparent ar trebui s fie egal cu greutatea lichidului din fiecare vas.
2.6.2. Aciunea fluidelor n repaus pe suprafee curbe deschise Torsorul forelor de presiune poate fi echivalat cu un sistem de trei ecuaii:
( )( )( )
=========
xoypp
xozpp
yozpp
SpSknpkFF
SpSjnpjFF
SpSinpiFF
z
y
x
dd,dd
dd,dd
dd,dd
(2.41)
unde: xoyzozyoz SSS d,d,d sunt proieciile elementului de suprafa dS al peretelui respectiv pe planele yoz, zox, xoy.
Rezult:
===xoy
z
zox
y
yoz
xS
xoypS
zoxpS
yozp SpFSpFSpF d;d;d (2.42)
Suportul forei xp
F trece prin centrul de presiune Cl, respectiv yp
F prin Cll i
zpF prin Clll, vectorii de poziie fiind determinai cu relaiile:
===xoy
xoylll
zox
zoxll
yoz
yozi
Sxoy
Sxoy
C
Szox
Szox
C
Syoz
Syoz
C Sp
Spr
rSp
Sprr
Sp
Spr
rd
d
;d
d;
d
d
(2.43)
Forele care nlocuiesc torsorul trec prin centrele de presiune ale proieciilor suprafeei curbe ale peretelui pe planele de coordonate.
a) Aciunea fluidelor uoare n repaus pe perei curbi deschii n acest caz p = ct., deci putem scrie: xoypzoxpyozp pSFpSFpSF zyx === ;; (2.44)
respectiv:
lllllllllllyozl GCGCGS
yoz
C rrrrrS
Sr
r ====
;;
d
(2.45)
Forele echivalente cu torsorul sunt paralele cu axele de coordonate i trec prin centrele de greutate ale proieciilor suprafeei pe planele de coordonate.
b) Aciunea fluidelor grele n repaus pe suprafee curbe deschise Considerm o curb deschis ABCD pe care o proiectm pe planele triedrului
triortogonal Oxyz, figura 14.
Statica fluidelor _____________________________________________________________________________________
31
Fig. 14 Aciunea unui fluid asupra unei suprafee curbe deschise Presiunea fluidului la nlimea z este p = z, n acest caz rezultnd:
===
==
==
VSzSzF
SzSzF
SzSzF
xoylllG
Sxoyp
zoxllG
Szoxp
yozlG
Syozp
xoy
z
zox
y
yoz
x
d
d
d
(2.46)
respectiv:
( )VrV
Vr
S
Szr
r
S
Szrr
S
Szr
r
GS
xoyG
xoy
C
zoxG
zox
CyozG
yoz
C
xoy
lll
xoylll
ll
zoxll
l
l
===
==
d
z
d
z
d;
z
d
S
SSyoz
(2.47)
n concluzie, efectul este un torsor nlocuit nu un sistem de trei fore n general neconcurente i paralele cu axele de coordonate.
Forele orizontale se determin analog cu cele exercitate de un lichid pe perei plani, unde se nlocuiete suprafaa cu proiecia suprafeei curbe.
Fora vertical este egal cu greutatea lichidului cuprins ntre suprafaa peretelui curb i planul xOy i trece prin centrul de grreutate al acestui volum.
2.7. Aciunea fluidelor n repaus pe suprafee curbe nchise tim c proiecia unei suprafee curbe nchise pe un plan este nul. Considerm
planul P cu normala sa N i o suprafa curb nchis S, figura 15. Proiecia suprafeei elementare dS pe un plan este: ( )NnSS p ,dd = Proiecia total este: ===
VSp VNdivSnNS 0dd (2.48)
BA
CD
Al Bl
Cl Dl
All
Dll
Bll
Cll
Alll Blll
ClllDlll
y
z
x
Mecanica fluidelor _____________________________________________________________________________________
32
Fig. 15 Proiecia unei curbe nchise pe un plan
unde s-a aplicat formula lui Gauss de transformare a integralei de suprafa n integral de volum, ( )0=Ndiv .
a) Aciunea fluidelor uoare n repaus pe suprafee curbe nchise innd cont c Syoz = Szox = Sxoy = 0, rezult c: 0===
zyx pppFFF (2.49)
Aciunea fluidului nu este totui nul, el dezvoltnd tensiuni unitare n peretele solid al rezervoarelor. Acest calcul este necesar pentru dimensionarea rezervoarelor.
b) Aciunea fluidelor grele n repaus pe suprafee curbe nchise Considerm o suprafa curb nchis aflat n fluid i un sistem de referin cu
planul x0y situat pe suprafaa liber, figura 16. Fig. 16 Aciunea fluidelor grele asupra suprafeelor curbe nchise innd cont de faptul c proieciile suprafeei curbe sunt nule, rezult: 0==
yx ppFF
Pentru determinarea forei zp
F , considerm suprafaa dat secionat de un plan paralel cu xOy i format astfel din dou suprafee deschise ACB i BDA. Fora de presiune pe poriunea superioar este ndreptat n jos lpF , respectiv pe poriunea
inferioar n sus, llpF .
kVFVVVFFFzz p
cilADB
cilACB
llp
lpp ==== ; (2.50)
Aceast for trece prin centrul de greutate al volumului V.
N
(P)
ds(V)S
n
lpF
llpF
A B
C
D z
y
0 x
Statica fluidelor _____________________________________________________________________________________
33
Un corp scufundat ntr-un lichid este supus unei fore verticale portante egal cu greutatea volumului de lichid dezlocuit de corp i trece prin centrul de greutate al volumului conform principiului lui Arhimede.
Trebuie fcut remarca c aceast for este o for de suprafa (de presiune) i este exprimat prin produsul dintre lichid i volum; nu este o for masic.
Principiul lui Arhimede este valabil i pentru gaze, dar datorit greutii specifice mici a gazelor este sesizat doar la corpuri solide cu greuti specifice mici cum ar fi dirijabilele.
Fora portant poate fi determinat i direct integrnd forele de presiune pe suprafaa nchis S:
==V
p VgradpSnpF ddS
(2.51)
tim c: p = z + po deci: kgradp = rezult: ==
Vp kVVkF d (2.52)
Principiul lui Arhimede rmne valabil i pentru corpuri parial scufundate n fluid, figura 16.
Fig. 16 Corp parial scufundat n lichid
Suprafaa este nchis S = Sl + Sll, So fiind aria delimitat de corp din suprafaa liber a lichidului. Pe cele dou suprafee Sl i Sll acioneaz presiunile p = po + z, respectiv po. ( ) =+=
l ll lS S S Soozop SnzSpnSpnSpnF dddd (2.53)
Integrala S
o Spn d este nul.Pe suprafaa liber a lichidului z = 0, deci i pe So. Putem
aduna la Fp termenul oS
Snz d , obinnd:
kVVgradzSnzF oVSS
p
ool
=== +
dd (2.54) unde Vo este volumul de lichid dezlocuit de corp numit i volum de caren. Principiul lui Arhimede poate fi stabilit i pentru un corp scufundat n dou lichide imiscibile, figura 17, cu greuti specifice 21 < , n repaus.
Presiunea n lichidul cu 1, pentru z < h este: ( ) zpzp o += 11 (2.55) Presiunea n lichidul cu 2, pentru z > h este: ( ) ( )hzhpzp o ++= 212 (2.56) Pentru suprafaa de separaie (z = h), cele dou presiuni sunt egale:
p1(h) = p2(h).
Sll
Sl
S0x
y z
0
Mecanica fluidelor _____________________________________________________________________________________
34
Fig. 17. Aciunea fluidelor grele asupra suprafeelor curbe nchise Rezultanta forelor de presiune ce acioneaz asupra corpului este:
( )[ ] ShzhnSnpSznSnpSpnSpnF
ll
lllllll
S
So
SSo
SSp
d
ddddd
21
121
==
adunm i scdem fora: ohS1 . Rezult: ( )[ ]
+ +=
ol
ollSS SS
p ShzhnSznF dd 211 (2.57) Deci fora portant este suma forelor portante pariale i are modulul egal cu
suma greutilor volumelor de lichid dezlocuite. Dac volumele Vl i Vll au centrele de greutate Gl i Gll, fora portant trece printr-
un punct G situat pe dreapta GlGll, ntre segmente fiind stabilit relaia:
11
22
1VV
FF
GGGG
p
pll
lz
== (2.58)
Rezultatul se poate extinde pentru mai multe lichide imiscibile. 2.8. Plutirea corpurilor Asupra unui corp cufundat ntr-un lichid acioneaz fora de greutate:
Fg = mV unde m este greutatea specific a corpului i fora portant Vc, unde este greutatea specific a lichidului i Vc volumul de caren (volumul de lichid dezlocuit).
Dac m > , corpul se scufund, deoarece Fg > Fp. Dac m = corpul rmne n echilibru indiferent, ceea ce corespunde plutirii scufundate (plutire submarin). Cazul cel mai important este acela cnd m < , deci Fg = Fp, cnd corpul plutete la suprafaa lichidului i creeaz un volum de caren, Vc = V.
Rezult: mV = Vc (2.59) Un astfel de corp liber scufundat parial n lichid se numete plutitor. Dac un vas are anumite ncrcri, acestora le corespunde mai multe carene
(poriunea din corp care este scufundat i creeaz volumul de caren) cuprinse ntre
1 < 2
0
y
z
x
Vl
Vll
Sl
Sll
S0
2
h
Statica fluidelor _____________________________________________________________________________________
35
carena minim pentru vasul gol i carena maxim corespunztoare vasului ncrcat la limit.
Centrul de greutate C al volumului de caren, n care se aplic fora portant se numete centru de caren. Pentru diferite nclinri ale plutitorului, la aceast greutate Fg, carenele au acelai volum, conform condiiei de echilibru, formele lor fiind n general diferite (cu excepia plutitorului omogen sferic) ns poziia centrului de caren se schimb.
Cele dou fore care echilibreaz un plutitor n poziia normal de plutire au acelai suport. n cazul unei nclinri provocate de perturbaii exterioare poziia centrului de caren se modific, din C se mut n Cl, figura 18.a.
Cele dou fore nu au acelai suport i formeaz un cuplu de fore, care tinde s readuc plutitorul n poziia normal, deci plutirea este stabil. n aplicaiile plutirii corpurilor se pune n general problema cutrii condiiilor necesare pentru o plutire stabil n cazul unor perturbaii care nu depesc anumite limite.
Definim planul de plutire planul suprafeei libere a lichidului, iar linia de plutire curba de intersecie a plutitorului cu planul de plutire. Aria figurii plane mrginit de linia de plutire se numete arie de plutire.
Raportm plutitorul la un sistem de axe astfel nct planul xOy s conin centrul de greutate G i centrul de caren C, corespunztor poziiei normale de plutire, axa Oz unind punctele C i G. Axa Oz este ax de plutire i este solidar cu plutitorul.
Axa Oy este de obicei ax longitudinal a plutitorului i se numete ax de nclinaie. Adncimea maxim a plutitorului sub linia de plutire se numete pescaj (h), figura 18. b.
Vericala forei portante care trece prin centrul de caren corespunztor punctului nclinat Cl, intersecteaz axa de plutire n punctul M, numit i metacentru. Cnd plutitorul tinde s revin n poziia vertical, metacentrul tinde ctre un metacentru principal Mo. Poziia metacentrului fa de centrul de greutate al plutitorului depinde de echilibrul plutirii.
Fig. 18 Plutirea corpurilor a. plutitor nclinat; b. plutitor n poziia de echilibru
Echilibrul este stabil atunci cnd metacentrul se afl deasupra centrului de greutate, GM > 0.
M
G C Cl
Fg
Fp GC
Fg
Fp
h
y
x
z
b)a)
Mecanica fluidelor _____________________________________________________________________________________
36
Echilibrul este instabil atunci cnd metacentrul se afl sub centrul de greutate, GM < 0. Echilibrul este indiferent atunci cnd M coincide cu G, GM = 0. Distana dintre punctul metacentric M i centrul de caren, pentru poziia vertical al plutitorului poate fi determinat conform teoremei metacentrului cu relaia:
c
gy
M VI
C = (2.60) unde gyI este momentul de inerie al ariei de plutire calculat fa de axa de nclinaie Oy, care trece prin centrul de greutate al ariei de plutire.
Condiia de stabilitate a plutirii, poate fi scris i altfel:
GM > 0 sau CM CG > 0 sau 0>CGVI
c
gy (2.61)
Cinematica fluidelor _____________________________________________________________________________________
37
3. CINEMATICA FLUIDELOR 3.1. Noiuni specifice Cinematica studiaz micarea fluidelor fr s considere forele care acioneaz
asupra lor determinnd micarea, efectund un studiu geometric al acesteia. De aceea, studiile cinematicii fluidelor sunt valabile att pentru fluide ideale ct i pentru cele reale.
Cinematica fluidelor are la baz ipoteza continuitii, caracterizat de parametrii care sunt funcii de timp i de punct, continue i derivabile. Ca metod de lucru, considerm c masa de fluid este format dintr-un numr foarte mare de particule fluide, analoage punctelor materiale din mecanica solidului.
Studiul cinematic al mecanicii fluidelor const n determinarea traiectoriilor, vitezelor i acceleraiilor particulelor de fluid. Se pot utiliza dou metode.
Prima metod care poart denumirea Lagrange, studiaz micarea fiecrei particule de fluid n raport cu un sistem fix Oxyz, poziia unei particule depinznd de timp i de coordonatele poziiei iniiale xo, yo, zo, corespunztoare timpului iniial to: ( ) ( ) ( )tzyxzztzyxyytzyxxx ,,,;,,,;,,, 000000000 === (3.1)
Componentele vitezei i acceleraiei se determin cu relaiile:
tzw
tyv
txu
==
= ;; (3.2)
22
2
2
2
2
;;tz
twa
ty
tva
tx
tua zyx
==
==
== (3.3)
Aceast metod este mai rar ntlnit i este necesar cnd se studiaz micarea unei particule individualizate.
A doua metod, metoda Euler, studiaz cmpul vitezelor n punctele spaiului ocupat de fluidul n micare i variaia acelor viteze n funcie de timp. Cmpul vitezelor este dat de relaiile:
( ) ( ) ( )tzyxwwtzyxvvtzyxuu ,,,;,,,;,,, === (3.4) sau ( )trVV ,= , unde x, y, z reprezint coordonatele punctelor spaiului i nu ale particulei fluide ca n cazul metodei precedente. Metoda Euler este mai simpl i utilizeaz teoria cmpurilor ca aparat matematic de studiu.
Componentele vitezelor se determin prin derivarea total a funciilor x = x(t), y = y(t), z = z(t):
tzw
tyv
txu
dd;
dd;
dd === (3.5)
Traiectoria particulelor se obine din integrarea ecuaiilor precedente: ( ) ( ) ( )tzyxzztzyxyytzyxxx ooooooooo ,,,;,,,;,,, === (3.6)
unde xo, yo, zo sunt constantele de integrare i reprezint coordonatele particulei la momentul iniial to.
Pentru determinarea cmpului acceleraiilor derivm pe u, v, w, innd cont c sunt funcii de x, y, z i de t, utiliznd regula de difereniere total:
zzuy
yux
xut
tuu ddddd
++
+=
Mecanica fluidelor _____________________________________________________________________________________
38
Pentru componentele acceleraiei obinem urmtoarele expresii:
+
++
==+
++
==+
++
==
zww
ywv
xwu
tw
twa
zvw
yvv
xvu
tv
tva
zuw
yuv
xuu
tu
tua
z
y
x
dddddd
(3.7)
nmulind relaiile cu i , respectiv cu j i k versorii axelor de coordonate i adunnd, obinem:
( ) VVt
VzVw
yVv
xVu
tV
tVa +
=+
++
==dd (3.8)
deci acceleraia este derivata total a vitezei. Se observ c aceast derivat total este
format din acceleraia local t
V i acceleraia convectiv (de antrenare):
( ) VVzVw
yVv
xVu
++
:sau . Acceleraia local reprezint variaia vitezei n puncte fixe n spaiu i este
caracteristic micrilor nepermanente (coordonatele punctelor sunt considerate ca invariabile). n micarea permanent acceleraia local este nul.
Cmpul acceleraiilor se poate exprima n coordonate carteziene i sub alt form echivalent:
+
++
+
==
+
++
+
==
+
++
+
==
xw
zuv
zv
ywvwvu
xtw
twa
zv
yww
yu
xvuwvu
xtv
tva
yu
xvv
xw
zuwwvu
xtu
tua
z
y
x
2dd
2dd
2dd
222
222
222
(3.9)
iar forma vectorial este:
VVrotVgradtVa x
2
2
++= (3.10)
Expresia vectorial pune n eviden partea potenial, 2Vgrad i partea
rotaional, VVrot x a acceleraiei consecutive. Cinematica fluidelor opereaz cu noiuni specifice micrii fluidelor. Curentul de fluid este o mas de fluid n micare. Linia de curent este linia curb care urmrind direcia de curgere este tangent la
vectorii vitez ai particulelor care la un moment dat coincid cu punctele de pe curba respectiv, figura 19. n general linia de curent nu este identic cu traiectoria.
n cazul micrilor nepermanente liniile de curent i modific forma n timp, iar n cazul micrilor permanente, cnd vectorii de vitez au poziii fixe n fiecare punct din spaiu, liniile de curent coincid cu traiectoriile i rmn aceleai n orice moment. n
Cinematica fluidelor _____________________________________________________________________________________
39
VnS n
d
V
general, liniile de curent nu se intersecteaz (o particul nu poate avea dou viteze diferite) dect n cazul unor puncte critice.
Fig. 19 Linie de curent
Ecuaiile difereniale ale liniilor de curent se obin din condiia ca vectorul d r = (dx, dy, dz) s fie paralel cu vectorul vitez, adic:
wz
vy
uxrV ddd:sau0xd === (3.11)
n cazul micrilor permanente vizualizarea traiectoriilor se poate realiza prin introducerea de suspensii fine n fluid.
Tubul de curent este suprafaa tubular format de liniile de curent, care trec la un moment date prin toate punctele unei curbe nchise. n micarea permanent, tubul de curent i pstreaz n timp forma i dimensiunile.
Firul de curent este fluidul din interiorul unui tub de curent elementar (cu seciune transversal foarte mic), care materializeaz o linie de curent.
Seciunea transversal a unui tub de curent (seciunea vie) este suprafaa normal pe toate liniile de curent care o strbat limitat de tub.
Perimetrul udat este lungimea conturului seciunii transversale a unui tub de curent, mrginit de perei rigizi.
Raza hidraulic este raportul dintre aria seciunii transversale i perimetrul udat:
PAR = (3.12)
Debitul unui curent de fluid printr-o suprafa S este fluxul vectorului vitez V, prin aceast suprafa, i reprezint limita raportului dintre volumul V care trece prin suprafaa S ntr-un interval de timp t, cnd aceasta tinde la 0, figura 20.
=== S nSt SVSnVtVQ ddlim
0 (3.13)
Debitul se poate defini i ca volumul de fluid care trece printr-o suprafa n unitatea de timp.
n afara acestui debit numit i debit volumic, se mai definesc debitul masic Qm=Q i debitul gravific Qg=Q.
Fig. 20 Suprafaa S fa de care calculm debitul Circulaia vitezei V de-a lungul unei curbe oarecare este:
==AB
tAB
AB SVsV dd (3.14)
A B C
D1v 3
v
4v
2v
Mecanica fluidelor _____________________________________________________________________________________
40
n cazul n care curba este nchis C, atunci circulaia vitezei poate fi exprimat cu ajutorul unei integrale de suprafa, S fiind o suprafa oarecare care se sprijin pe curba C, (teorema lui Stokes):
==S
SnVrotsV ddC
(3.15)
Vrtejul unei particule de fluid este vectorul definit prin relaia:
wvu
zyx
kji
Vrot
==
21
21 (3.16)
i reprezint viteza unghiular medie de rotaie a particulei n jurul unei axe ce trece prin centrul ei de greutate. Componentele sale sunt:
=
=
=yu
xv
xw
zu
zv
yw
zyx 21;
21;
21 3.17)
Linia i tubul de vrtej sunt noiuni analoage cu cele definite anterior, n care se nlocuiete vectorul .cu V
Tipuri de micri specifice fluidelor Micrile fluidelor se pot clasifica dup diferite criterii. Din acestea menionm: - dup modul de variaie n timp a parametrilor micrii, pot fi permanente
(viteza, presiunea, densitatea sunt constante n timp) i nepermanente, cnd parametrii sunt n funcie de timp;
- dup modul de desfurare a micrii n lungul curentului se deosebesc micri uniforme (liniile de curent sunt paralele i rectilinii) i neuniforme, cnd liniile de curent au o form oarecare, de-a lungul crora vitezele variaz ca mrime i direcie;
- dup cmpul vitezelor se deosebesc micri poteniale, cnd exist o funcie , numit funcie de potenial, astfel nct gradV = (aceste micri se mai numesc i micri irotaionale 0== gradrotVrot ) i micri nepoteniale sau rotaionale;
- dup structura fizic a curgerii unui fluid real sunt micri laminare care se produc la viteze relativ mici i n care straturile se mic paralel unele cu altele i micri turbulente la care particulele diferitelor straturi se amestec ntre ele i se deplaseaz dup traiectorii neregulate.
3.2. Micarea unei particule fluide Micarea unui fluid este mult mai complicat dect micarea unui rigid solid (care
poate fi descompus ntr-o micare de translaie i o rotaie n jurul unei axe instantanee de rotaie) deoarece fiecare particul de fluid sufer o translaie, o rotaie i o deformaie.
Considerm o particul de fluid, figura 21.
Cinematica fluidelor _____________________________________________________________________________________
41
z
x y
rr d+r
V
rdM
rotV
MlV
lV defV
Fig. 21 Deplasarea unei particule fluide n punctul M (x,y,z) viteza are componentele u, v i w, iar ntr-un punct nvecinat M' (x+dx,y+dy,z+dz) viteza V' are componentele:
d d d
d d d
d d d
u u uu u x y zx y zv v vv v x y zx y zw w ww w x y zx y z
= + + + = + + + = + + +
(3.18)
Dac adunm i scdem la prima ecuaie termenii yxv d
21 i z
xw d
21 , putem s-o
scriem sub forma: 1 1d d d2 2
1 1d d2 2
u u v u wu u x y zx y x z x
u w v uz yz x x y
= + + + + + +
(3.19)
n mod analog se obin relaiile pentru v' i w'. Se observ c ultimele dou paranteze sunt componentele vitezei unghiulare de rotaie, deci ele reprezint rotaia particulei. Se noteaz:
=
=
+
=
=
+
=
xua
axv
zua
ayu
xva
xx
xzzx
yxxy
21
21
(3.20)
unde axy = ayx, azx = axz, reprezint viteze specifice de deformare unghiular iar axx i analog ayy, azz reprezint viteze specifice de deformare liniar.
Cu aceste notaii ecuaia (3.19) devine: d d d d dxx xy xz y zu u a x a y a z z y = + + + +
componentele v' i w' se determin n mod analog. Pentru determinarea semnificaiei lui axx se consider un element de fluid liniar
AA', paralel cu Ox, de lungime dx, fig. 22.
Mecanica fluidelor _____________________________________________________________________________________
42
Diferena deplasrilor relative ale captului elementului liniar n intervalul de timp dt reprezint dilatarea sau contractarea acestuia i este:
Fig. 22 Deformarea liniar a unei particule fluide
txxututx
xuu ddddd
=
+ (3.21)
deci, viteza specific de deformaie liniar este:
txxu
txxuaxx dddd
1=
= (3.22) Pentru interpretarea termenilor de forma ayz, azy se examineaz micarea unei
particule de form paralelipipedic, a crei seciune cu un plan paralel cu yOz este dreptunghiul ABCD, figura 23.
Fig. 23 Deformaia unghiular a unei particule fluide ntr-un interval de timp dt, particula se deplaseaz i se deformeaz ocupnd
poziia A'B'C'D'. Dac anulm translaia i rotaia, aducem particula n poziia A''B''C''D''. Deplasarea relativ DD'' se datoreaz diferenei dintre vitezele punctelor A
i D, zzvvv AD d= i are mrimea: d dvDD z t
z = . Analog rezult:
d dwBB y ty
= . n ipoteza unor deplasri mici, deformaia medie a unghiului drept BAD este:
( )1 1 1 d2 2 2
DD BB v w tAD AB z y
+ = + = + (3.23) Deci viteza de deformaie unghiular are expresia:
yzayw
zv
t=
+
=
+21
d21 (3.24)
z
y0
D
A=A''
B''
B
CC''D''
D'
B'
A'
C'
AAl
x
z
y
u
xxuu d+
dx
Cinematica fluidelor _____________________________________________________________________________________
43
( ) tzxyyvv dddd
+
( ) tzyxyuu dddd
+
( ) tyxzyww dddd
+
wdxdydt
udydzdt
vdxdzdt z
x
y
Rezult c axx, ayy, azz reprezint vitezele de deformaie liniar, iar mrimile axy, axz, ayz reprezint viteze specifice de deformaie unghiular. Revenind la expresiile vitezelor u', v', w', putem determina vectorul vitez V nmulind respectiv cu kji ,, i adunnd. Dac considerm funcia scalar:
( )zyazxayxaxaxaxa yzxzxyzzyyxx dd2dd2dd2ddd21 222 +++++= (3.25) vectorul este: xdV V V r grad = + + cu precizarea c gradientul funciei se calculeaz n raport cu dx, dy, dz.
Funcia se mai numete funcie de deformaie, iar cuadrica corespunztoare este un elipsoid de deformaie. Utiliznd relaia vectorial anterioar, putem formula urmtoarea teorem: dac se cunoate micarea particulei fluide M( r ), micarea unei particule vecine Ml( rr d+ ) se compune dintr-o micare de translaie definit de viteza V a punctului M, dintr-o micare de rotaie definit de viteza unghiular Vrot
21=
n jurul unei axe care trece prin M i dintr-o micare de deformaie a cuadricei = ct. cu centrul n M i care trece prin Ml, micare compus dintr-o deformaie liniar definit de mrimile axx, ayy, azz i deformaie unghiular definit de mrimile axy, axz, ayz, figura 21.
Aceast teorem poart numele de teorema lui Cauchy-Helmholtz. 3.3. Ecuaia continuitii Aceast ecuaie este expresia matematic a principiului conservrii masei de fluid
n micare. Se consider cazul general al unui fluid compresibil cu (x,y,z,t) n micare
nepermanent cu V (x,y,z,t). Alegem un volum de fluid de forma unui paralelipiped cu muchiile dx, dy, dz, figura 24. Relaia care exprim continuitatea fluidului se obine egalnd variaia masei de fluid din volumul considerat cu diferena dintre masa care intr n acest volum i masa de fluid care iese din el, n acelai interval de timp.
Fig. 24 Particul fluid paralelipipedic
Mecanica fluidelor _____________________________________________________________________________________
44
La timpul t masa de fluid este dxdydz, iar la timpul t + dt devine zyxt
tdddd
+ . Variaia masei este: tzyx
tm ddddd
= Diferena dintre masa de fluid intrat i cea ieit n intervalul de timp dt,
considernd cele trei direcii este:
( ) ( ) ( ) tzyxzw
yv
xum ddddd
+
+= (3.26)
Din egalarea celor dou expresii rezult: ( ) ( ) ( ) ( ) 00 =+=+++ Vdivtzwyvxut sau (3.27) Dac efectum derivatele produselor u, v, w, obinem alt form a ecuaiei
continuitii:
0sau 0 =+=
+
++ Vdiv
dtdp
zw
yv
xu
dtdp (3.28)
n cazul micrii permanente a unui fluid compresibil ,0=
tp deci se obine:
( ) ( ) ( ) ( ) 0sau 0 ==++ Vdivzwyvxu (3.29) Pentru un fluid incompresibil, = const. n micare permanent sau
nepermanent, ecuaia continuitii are expresia:
0 sau 0 ==+
+ Vdiv
zw
yv
xu (3.30)
adic cmpul vitezei unui fluid incompresibil este solenoidal.
Ecuaia continuitii pentru un tub de curent oarecare Considerm un volum de fluid V(t) n interiorul unui tub de curent limitat de dou
seciuni S1 i S2, respectiv la intrare i la ieire. Masa din interiorul tubului este:
=V
Vm d . n intervalul de timp dt densitatea fluidului crete cu ,d tt
iar masa de
fluid crete cu = V tVtm ddd . Aceast cretere a masei se datoreaz faptului c prin
suprafaa de intrare S1 ptrunde n intervalul de timp dt o mas de fluid mai mare dect cea care iese din acelai interval de timp prin suprafaa de ieire S2, diferena fiind:
tSVSVmS
nS
n dddd2
2
1
1 2211
=
Egalnd cele dou expresii ale variaiei masei rezult:
21
2
2
1
1 2211ddd mm
Sn
Sn
V
QQSVSVVt
== (3.31)
Ecuaia (3.31) se interpreteaz astfel: variaia masei de fluid dintr-un tub de curent n unitatea de timp este egal cu diferena dintre debitul masic intrat n tubul de curent i debitul masic ieit din tub.
Cinematica fluidelor _____________________________________________________________________________________
45
ds
vSdt
( ) tvSs
vS d
+
2
1
n micarea permanent, 0=
t i rezult c:
MS
nS
n QSVSV ==2
2
1
1 2211dd (3.32)
deci n micarea permanent a unui fluid compresibil debitul masic prin orice seciune este constant. Dac fluidul este incompresibil indiferent de felul micrii debitul volumic este constant.
Considerm un tub de fluid elementar limitat de dou seciuni transversale situate la o distan infinit mic ds. ntr-un interval de timp dt prin seciunea l intr masa
VSdt, iar prin seciunea 2 iese masa ( ) tsVSs
VS dd
+ , figura 25.
Fig. 25 Tub de fluid
Variaia masei de fluid n intervalul de timp dt este dat de diferena dintre cele
dou mase, ( ) tsVSs
m ddd = .
Pe de alt parte masa iniial dV = Sds sufer n timpul dt variaia ( ) ( ) tstSt
txSm ddddd
==
Egalnd cele dou expresii diferite ale variaiei masei, rezult:
( ) ( ) 0=+
VSs
St
(3.33) care reprezint ecuaia continuitii pentru un fluid compresibil n micare nepermanent. Dac fluidul este compresibil i micarea este nepermanent ntr-un tub
de curent cu perei solizi, atunci 0=
tS i ecuaia devine:
( ) 0=+
VSst
S (3.34) Dac seciunea este constant n lungul tubului, 0=
sS i rezult:
( ) 0=+
Vst (3.35)
Pentru un fluid compresibil n micare permanent ecuaia ia forma: .const== MQVS (3.36)
Mecanica fluidelor ___________________________________________________________