Université Med BOUDIAF Mia Faculté de Techaologies Département de Génie mécanique Cours d’MDF-Approfondie : Master 1 Energétique _Resp-M: Bouaouina-t Chapitre. 3: Dynamique des fluides réels 3-1-Quelques solutions exactes de Navier-Stokes (écoulements paralléles) Les Equations de Navier-Stokes sont des équations aux dérivées partielles non linéaires dont la résolution analytique est généralement impossible. Cependant, pour certaines configurations simples, elles peuvent se réduire a un systéme linéaire dont lintégration permet d’obtenir des solutions exactes. Cas des écoulements paralléles unidirectionnels (une seul composante de vitesse non mulle— trajectoires paralléles des particules fluides) de fluide isovolume, Par ailleurs, les problémes traités seront tels que les conditions aux limites associées aux frontiéres sont simples a appliquer mathématiquement. Ainsi, nous allons essayer de résoudre. {pa P+ v8 (9.1) On P* et p: pression et viscosité dynamique. a Le systéme (9 .1) est compose de quatre équations pour les quatre inconnues P* et ii. Contrairement aux. problémes en dynamique la force -¥P" est a déterminer comme pattie intégrale de la solution, car elle n'est pas donnée, VP + pge; = V(P + pgz) = VP" 3. coulements paralléles 3-2-1-Equations pour les écoulements paralléles en canal > Considézons, un écoulement rectiligne paralléle a Paxe Ox dans le repére cartésien(8,,8,,8,) ‘Uy = Uy(x,y,Z,0) » Les composantes du vecteur vitesse sont: } Uy = 0 u,=0 Pour un fluide isovolume, l’équation dé continuité se réduit a’ oe =0 Les profils de vitesse restent done identiques dune section droite & une autre: 1 = uy(y,z,t)@y En définissant le potentiel de pesanteur © par G=—grad© et la pression étoilée ou dynamique P* par P*=P+pO= P(x)+pgy «il vient pour un fluide iso volume et on négligeant les forces de volume devant celle de viscosité: ty _ . Me er) 7] | ze —_— —> oo z : Il appatait/que *la ae Aloilée ne a évoluer a dans la direction de I’écoulement: P* = P*(x, t). 3-2-2- Equations pour les écoulements paralléles en conduite =n écoulement dans un conduit cylindrique paralléle & l’axe Oz, vewenen Dans un repére de coordonnées cylindriques(é,.,85,2,) = u, = — Les composantes d’un vecteur vitesse sont: {i 0 = u, = u(r, 0,2, t) 5 at Pour un fluide isovolume, I’équation de continuité se réduit dans cette situation a: = 0 a =u,(r,0,0e; Les profils de vitesse sont inchangés d'une section droite a une autre= L’équation de quantité de mouvement en projection sur les trois axes se réduit & w oF 3-2-3- Equations pour les écoulements paralléles en rotation un écoulement plan de rotation autour d’un axe Oz, Dans un repére de coordomnées cylindriques (8,,89,2,) u,=0 Les composantes du vecteur vitesse sont: }Ug=U9(r,0,2,t) u,=0 Pour un fluide isovolume, l’équation de continuité se réduit: oe =0 Les profils de vitesse sont done identiques dans tous les plans azimutaun: % = up(r,z, 0) > L’équation de quantité de mouvement en projection sur les trois axes s’écrit Ainsi P* = P¥(x,0). On peut montrer que ty = uo(1) Gj _ oP La symeétrie eylindrique impose alors 6P*/00= Le mouvement satisfait done le systéme suivant 3. coulement a faible Nombre de Reynolds On s‘intéresse ici aux éconlements oules forees d'inertie (forces daccélération) sont négligeables par rapport aux forces de pression et foree de viscosité, L’équation de Navier-Stokes s’écrit alors f =O equationdecontnulté "On appelle cette équation Péquation de Stokes, VP" = pV" equation dynamique ; >» (Via S+24250 sil’on peut permuter les opérateurs div et V ax ay ae G1) V.(VP") = n02(V.ii) => V2P' =0 On appelle écoulement rampant tout écoulement représente par l'équation (3.1). avec Re << 1 De plus nous supposerons que le nombre de Stokes qui compare le terme instationnaire aux termes fel _ 2 (On va voir des applications assez diverses et plus au moins dizectes de ces équations dans des problémes d’ingénierie Econlement dans un massif poreux, |’écoulement de Stokes autour d’une sphére, sédimentation de particules et lubrification dun palier. visquenx est trés petit: St = «1 3-3-1-La sédimentation : Sédimentation d'une particule sphérique de rayon R animé d’une vitesse up dans un fluide de viscosité pt La réduction des équations de Stokes permet de trouver la vitesse du fluide autour de la particule, done la force exereée par le fluide sur la particule : F = 6mpR. tu,Mg 32, ;Coordonntes sphérique autour d'une sphére Fig 33. > lignes de courant pour un écoulement de Stokes La vitesse de sédimentation ou de chute visquense ou inertielle, Dans le premier cas il y a égalité entre le poids apparent (poids corrigé de la poussée D’Archiméde) et la force de Stokes 4 (Py ~ Pp); AR? = OMAR. Ustokes > Up = Ustokes = 5 Pour que le nombre de Reynolds de chute soit petit et que l'on puisse utiliser la formule de la force de i gery Poser rg) 3-3-2-Lubrification Hydrodynamique Dans les écoulements de Stokes (Re <<1) et dans les écoulements paralléle le terme Inertiel (0.V)0 de I’équation de Navier Stokes est soit négligeable soit exactement nul L’écoulement est laminaire, les lignes de courants sont pratiqnement paralléles aux parois est done presque paralléle a l'axe des x, uU et uy-U8. Si le nombre de Stokes St<<1 on peut négliger les termes instationnaires, tout se passe comme si localement l’écoulement était paralléle. On trouve donc toujours, selon les conditions aux limites de pression et de vitesse , une superposition locale d’un écoulement de Couette et de Poiseuille, En tudiant la projection de Navier Stokes sur x et y , On obtient le systéme 2Rilstoken — Stokes il faut que : Re = 10P_ au, pox dy wa O- et it itr — =0- et on peut montre poy P Remarque ; Cette équation est un compromis entre la condition des écoulements paralléles (6-0 et Re en principe quelonque) et la condition des écoulements de Stokes (Re<<1 et 6 quelconque), igh. Ecoulement ind parle mouvement relat de aise par appt au palin Fig3-S Distribution de @(8)-0(0) pour 0.2mm, bs-0 1mm rvs PIO" kgim La force de pression s’exergant sur le patin pour unité d’envergure, est donnée par mnt GuU (RL yh — ha Ba [par = OP (ne La portance du patin est égale & F'cos@ ~ F car @ <1. Par exemple, pour les données utilisées en figure ELM on trouve F = 159 N.