UNIVERSITATEA POLITEHNICA BUCURESTIFACULTATEA ENERGETICA

PROIECTCurgerea printr-un canal cu ngustarea seciunii

STUDENT: Neagu CosminGRUPA: MS7

Curgerea printr-un canal cu ngustarea seciunii

Dinamica fluidelor studiaz micarea fluidelor i interaciunea acestora cu corpurile solide, innd seama de forele care determin starea de micare i de transformrile energetice produse n timpul micrii. La curgerea fluide, fenomen datorat viscozitii fluidelor i interaciunii fluidelor cu contururile solide.Pentru caracterizarea micrii unui fluid este necesar cunoaterea: distribuiei vitezelor, distribuiei presiunii, distribuiei temperaturii n masa de fluid.Aceti parametri depind de o serie de factori ca: forma i dimensiunile spaiului de curgere, debitul fluidului, cmpul de fore care acioneaz asupra fluidului, etc.Clasificarea micrii fluidelorDup condiiile de variaie n timp a parametrilor locali: curgere staionar curgere nestaionarDup condiiile de variaie n spaiu a parametrilor locali: curgere unidimensional(unidirecional) curgere bidimensional (plan sau axial-simetric) curgere tridimensional (spaial)Dup condiiile de contact cu suprafeele solide care delimiteaz spaiul de curgere: curgere sub presiune curgere cu suprafaliber Dup natura cmpului vectorial al vitezelor: curgere irotaional(potenial) curgere rotaionalDup mecanismul curgerii: curgere laminar curgere turbulent.Modelarea curgeriiTITLE 'Curgerea irotationala printr-un canal cu ingustarea sectiunii' { the problem identification }COORDINATES cartesian2 { coordinate system, 1D,2D,3D, etc }VARIABLES { system variables }fi {potentialul vitezei}! SELECT { method controls }DEFINITIONS { parameter definitions }lx=5 ly=5vx0=4 {viteza de intrare}p0=10^5 {Presiunea atmosferica }ro=10^3 { Densitate}vx=dx(fi) { Componenta vitezei pe directia x}vy=dy(fi) { Componenta vitezei pe directia y}v=vector(vx,vy) vm=sqrt( vx^2+ vy^2) { Viteza }p=p0+ 1/2*ro*(vx0^2-vm^2) { Presiunea }div_v=dx( vx)+ dy( vy) ! INITIAL VALUESEQUATIONS { PDE's, one for each variable }

dxx( fi)+ dyy( fi)=0! CONSTRAINTS { Integral constraints }BOUNDARIES { The domain definition } REGION 1 { For each material region } START (0,ly) value(fi)=-vx0line to (0,-ly) natural(fi)=0line to(lx,-ly) natural(fi)=0line to(lx,-ly/2)natural(fi)=0line to(2*lx,-ly/2) value(fi)=0line to(2*lx,ly/2)natural(fi)=0line to(lx,ly/2) natural(fi)=0line to(lx,ly) natural(fi)=0line to close! TIME 0 TO 1 { if time dependent }MONITORS { show progress }PLOTS { save result displays }contour( fi)vector( v) contour( vm) painted contour( p) painted END

Vectorii prezeni n figur arat repartiia veitezei.

Se poate observa c viteza la intrare este mult mai mic fa de cea la ieire. Viteza minim este prezent n colurile unde se ngusteaz seciunea iar maximul se nregistreaz n colurile n care seciunea se liniarizeaz.

Se poate observa distribuia presiunii, aceasta prezint o variaie mic ntre intrare i ieire. Valoarea cea mai mic este semnalat n aceleai puncte n care este nregistrat viteza maxim.