Dynamique Des Fluides Réels
Dynamique Des Fluides Réels
Dynamique Des Fluides Réels
Réalisé par :
Elaziouti Abdelkader
ω Vitesse angulaire
C Couple
𝝏𝒗
Gradient de vitesse
𝝏𝒏
μ Viscosité dynamique.
ν Viscosité cinématique
𝜸̇ Vitesse de cisaillement
μ app Viscosité apparente
K Constante de consistance
n Indice d’écoulement ou de comportement
λ Coefficient de perte de charge linéaire
K' Constante de perte de charge singulière
ε Rugosité moyenne de la conduite.
εr Rugosité relative
ρ Masse volumique du fluide
D Diamètre de la conduite
∆𝐏𝐭
Perte de charge par unité de longueur
𝐋
Re Nombre de Reynolds
Rec Nombre de Reynolds critique
Fr Nombre de Froude
Eu Nombre d’Euler
Ma Nombre de Mach
St Nombre de Stokes
Pr Nombre de Prandtl
Ca Nombre de capillarité ou nombre capillaire
Chapitre I : Viscosité. Lois de comportement
G ère Année Master: Gestion
UMT/FT/ GP/1 G des Déchets Moodule :Dynaamique des fluides
f réelss
I-1Mise en
n évidencee : viscosim
mètre de Couette
C
Le viscoosimètre (F
Figure I-11)est consttitué de deeux cylindrres coaxiaaux l'un est mis en
rotation avec la vittesse anguulaire ω. L'autre est entraîné parr les forcess de viscossité et est
maintennu immobille par l'appplication d'uun couple C. L'espacce entre less deux cylin
ndres est
faible ett rempli duu fluide donnt on veut mesurer
m la viscosité.
On appeelle R1 et R2 les rayyons des cyylindres ett e la difféérence R1 – R2≈e; e est très
petite dee telle façoon que l'on puisse connsidérer la répartitionn des vitessses commee linéaire
entre less deux cylinndres.
On monntre que la force
f de froottement est proportionnelle au gradient dde vitesse soit
s :
(1
1)
C
Cylindre maintenu
immobile
h R1
e R2r
R2 Fluid
de
Cylindree en rotation
Figure I-1
I : Viscoosimètre dee Couette.
F
Figure d fluide placée entree deux plaqques.
I-22 :Couche de
Elaziouti Ab
bdelkaderr
UMT/FT/ GP/1ère Année Master: Gestion des Déchets Module :Dynamique des fluides réels
Remarques
9 Lorsque la vitesse n’est pas trop grande, les couches de fluide ne se mélangent pas. C’est
un écoulement laminaire.
9 Si la vitesse est constante, le mouvement est rectiligne uniforme. Dans ce cas, F est la
force de frottement.
Figure I-3 :La répartition des vitesses de fluide placée entre deux plaques
et
Remarque
9 R1 et R2 ≈ 1 à quelques centimètres et entrefer R2-R1 ≈ quelques millimètres.
cisaillement constant si R1-R2<< R2 (R1-R2)/R2<0,15 entrefer étroit ;
Elaziouti Abdelkaderr
UMT/FT/ GP/1ère Année Master: Gestion des Déchets Module :Dynamique des fluides réels
9 Forme conique de la base du cylindre intérieure pour éviter les effets parasites. Il
existe deux types de rhéomètres de Couette:
9 À cisaillement imposé: on impose la vitesse de rotation et on mesure le couple C.
9 À contrainte imposée: on impose le couple et on mesure la vitesse de rotation
obtenue.
L'expérience de Couettepermet la mesure du coefficient de viscosité par la mesure du
couple C :
(2)
I-2Dimensions et unités
a-Viscosité dynamique
De la relation on en déduit la dimension de la viscosité dynamique.
=M
Elaziouti Abdelkaderr
UMT/FT/ GP/1ère Année Master: Gestion des Déchets Module :Dynamique des fluides réels
b-Viscosité cinématique
Dans de nombreuses formules apparait le rapport de la viscosité dynamique μ et de la
masse volumique ρ. Ce rapport est appelé viscosité cinématique ν.
(3)
2 -1
Dimensions : =L T
2 -1
Unités : Système international (SI), l’unité de la viscosité est m .s .
2 -1 4
En Système CGS, l’unité est le Stockes( symbole : St) :1m .s =10 st
Quelques valeurs typiques
Quelques valeurs caractéristiques de viscosités dynamiques et cinématiques sont récapitulées
dans le tableau I-1.
Elaziouti Abdelkaderr
G ère Année Master: Gestion
UMT/FT/ GP/1 G des Déchets Moodule :Dynaamique des fluides
f réelss
I-3Variatiion de la viscosité
v dyynamique avec la températuree et la pression
Cas dess gaz
Influencce de la prression , (4)
oùlest lee libre parccours moyeen,
ρla massse volumiqque et c la vitesse
v moyenne des moléculess.
En génééral pour lees gaz, la prression infl
flue peu surr la viscosiité.
Influencce de la températuree
Soit μ0 la
l viscositéé à la temppérature T0
(5)
oùK et C sont
s des coonstantes qui
q dépenddent de la nature
n du gaaz
Cas dess liquides
Influencce de la prression
Influencce de la températuree
av
vec
Elaziouti Ab
bdelkaderr
G ère Année Master: Gestion
UMT/FT/ GP/1 G des Déchets Moodule :Dynaamique des fluides
f réelss
a-Fluide newtonien
n
Entre deuxx couches successivees de fluidde en écoullement uniidimensionnnel à la vitesse
v , il
existe des contraintees tangentielles à l’éccoulement qui accélèèrent la coouche la pllus lente et
ralentissennt la coucche la pluus rapide. Par définition d’unn fluide neewtonien, les forcess
visqueusess sont propportionnellees à la difféérence de vitesse
v c’esst-à-dire auu gradient de
d vitesse:
(6)
(7)
Elaziouti Ab
bdelkaderr
UMT/FT/ GP/1ère Année Master: Gestion des Déchets Module :Dynamique des fluides réels
σt
0<n<1
Elaziouti Abdelkaderr
UMT/FT/ GP/1ère Année Master: Gestion des Déchets Module :Dynamique des fluides réels
Il existe une contrainte seuil au-dessous de laquelle il n’y a pas d’écoulementtel que le
dentifrice.
Elaziouti Abdelkaderr
Chapitre II : Dynamique des fluides visqueux
incompressibles :
Equation de Navier Stokes
G ère Année Master: Gestion
UMT/FT/ GP/1 G des Déchets Moodule :Dynaamique des fluides
f réelss
II-11 Écoulem
ments lamin
naires et écoulementtsturbulen
nts
Etudier l’éécoulemennt d’un fluide
fl réel revient à résoudre l'équation fondamen
ntale de laa
dynamiquee (équationn de Navier-Stokes):
∆ (9)
Mais en pratique,
p ceette équatioon ne peutt se résoud
dre analytiqquement qqu’en introd
duisant dess
hypothèsess de travaiil qui perm
mettront dee résoudre cette équaation dans le cadre de
d régimess
d'écoulemeent particuuliers :le réégime lamiinaireet le régime tu
urbulent.
9 Onn dit qu’un écoulemennt est lamiinairelorsq
que le mouvvement dess particules fluides see
faitt de façonrrégulière ett ordonnée.
9 L’éécoulemennt est tu
urbulentlorsque le déplacem
ment est irrégulierr et quee
dessfluctuationns aléatoirres de vitessse se supeerposent auu mouvem
ment moyen
n du fluidee.
Enn régime turbulent,
t on devrra utiliser des relations empiiriques gén
néralemennt
détterminées expériment
e talement.
L’une des premièress analyses de la trannsition d’un bulence est
n régime laminaire vvers la turb
uées par O.
basée surddes observvations d’écoulementts dans une conduitee cylindriqque effectu
Reynolds en
e 1883.
Le montagge expérim
mental schém
matisé sur la figure II-2
I compoorte un réseervoir de liiquide souss
pression débouchant sur un connduitcylinddrique. Un tube mincee permet l’injection de
d colorantt.
Lorsque l’’écoulemennt est laminnaire, le fiilet decolorrant reste mince,
m réggulier et paarallèle à laa
paroi du cylindre.
c E écoulem
En ment turbullent, le collorant estraapidement dispersé. Dans cettee
situation, une
u mesuree de la com
mposante de
d vitesse axiale
a monntre quecellle-ci fluctu
ue de façonn
aléatoire dans
d l’espacce et dans le temps.
Elaziouti Ab
bdelkaderr
G ère Année Master: Gestion
UMT/FT/ GP/1 G des Déchets Moodule :Dynaamique des fluides
f réelss
À partir de
d ces obseervations, les travauxx de Reyn
nolds ont alors
a permis de mon
ntrer que laa
transition du régim
me laminaiire au réggime turbu
ulent (Figgure II-2))n'est pas seulemennt
conditionnnée par le débit maais dépendd d'un enssemble de paramètrees qu'il co
onvient dee
regrouper pour forrmer une quantité sans dim
mension quu'on appeellera « no
ombre dee
Reynolds »Eq. (10).. Ce nombbre prend en
e compte la vitesse moyenne de l’écouleementV, lee
diamètre Dde
D la connduite (ou tout autre paramètree de longueur caractéérisant l'éccoulement),
ainsi que les propriéttés intrinsèèques du fluuide (masse volumiquue ρet visccosité μ) :
(10)
où ρ : masse
m voluumique [ρ ] = ML-3
μ : viscosité
v μ ] = ML-1 T-1
[μ
ν= µ/ρ : viiscosité cinnématique [ ν] = L 2T-1
V : vitesse [v ] = LT-1
D : diamètre [D]
[ =L
10
Elaziouti Ab
bdelkaderr
G ère Année Master: Gestion
UMT/FT/ GP/1 G des Déchets Moodule :Dynaamique des fluides
f réelss
Figu
ure II-3 :L
Les régimess d’écoulem
ment
Remarquee :
En prenannt certainees précautiions pour éviter au maximum
m les pertuurbations (vibrations
( s,
aspérités de
d la condduite...) l'éécoulementt peut restter laminaiire au-delàà deRec, ju
usqu'à unee
certaine lim
mite (imprrécise) qui dépendra essentiellem
e ment du niveau de prrécautions prises
p pourr
repousser l'apparition
l n des turbuulences.
11
Elaziouti Ab
bdelkaderr
UMT/FT/ GP/1ère Année Master: Gestion des Déchets Module :Dynamique des fluides réels
dissipation d'énergie due aux frottements visqueux, ces pertes de charges prendront place
dans la formulation d'une équation deBernoulligénéralisée.
Partons de l'équation de Navier-Stokes(Eq. (9)) établie pour l'écoulement d'un fluide
newtonien incompressible :
∆ Ω
On retrouve ainsi une formulation semblable à celle conduisant pour un fluide parfait à
l'équation de Bernoulli.
De la même manière, le produit vectoriel du vecteur vitesse avec le rotationnel de la vitesse
s'annule lorsque l'équation est considérée sur l'ensemble des points constituant une ligne de
courantEq. (16).
∆ Ω
Ω Ω
Il reste donc que le long d'une ligne de courant, l'écoulement stationnaire d'un fluide de
viscosité non négligeable obéit à l'équation suivanteEq.(17) :
12
Elaziouti Abdelkaderr
UMT/FT/ GP/1ère Année Master: Gestion des Déchets Module :Dynamique des fluides réels
Puis la projection sur chacun des axes d’un repère cartésienEqs. (18-23) :
(18)
• • •
• (19)
(20)
∆ (22)
∆ (23)
∆ (25)
∆ (26)
est la pression totale (également appelée « charge »), notée Pt, , somme
13
Elaziouti Abdelkaderr
UMT/FT/ GP/1ère Année Master: Gestion des Déchets Module :Dynamique des fluides réels
On en déduit que la pression totale n'est fonction que de x, et qu'elle varie en obéissant à
l'équationEq. (28) :
(28)
Par ailleurs, l'écoulement étant conservatif, l'équation de continuité Eq. (29) s'applique dans
sa forme la plus simple:
(29)
, ,
Il s'avère alors que le membre de gauche ne peut dépendre que de xalors que celui de
droite ne peut dépendre que de y et z. La conséquence est donc que les deux
membres sont égaux à une même constante , , . Ceci permet de montrer que la pression
totale Eq. (31) varie linéairement avec x:
Remarque :
9 Pour un fluide parfait dont la viscosité est négligée, la pression totale ne varie pas. En
effet, la constante étant proportionnelle à μ, elle devient nulle, l'équation de Bernoulli
Eq. (32) s'applique le long d'une ligne de courant :
(32)
On dit alors que la charge est constante ou bien encore qu'il n'y a pas de perte de charge.
9 Pour un fluide réel (dont la viscosité n'est pas négligeable), l'écoulement laminaire est
caractérisé par une variation linéaire de la charge avec la distance parcourue. En effet,
puisque les frottements visqueux sont responsables d'une dissipation d'énergie et que la
pression totale n'est autre qu'une énergie volumique, il s'ensuit logiquement que la
pression totale doit diminuer avec la progression de l'écoulement.
14
Elaziouti Abdelkaderr
G ère Année Master: Gestion
UMT/FT/ GP/1 G des Déchets Moodule :Dynaamique des fluides
f réelss
Figure II-4
I : Evoluution de la pression totale (charrge) avec la
l progression de l'éco
oulement
∆
∆
∆ ∆ (33)
régulièresE
Eq. (34):
∆ 0 (34))
∆ (35)
15
Elaziouti Ab
bdelkaderr
UMT/FT/ GP/1ère Année Master: Gestion des Déchets Module :Dynamique des fluides réels
∆ (36)
Expérience de Nikuradze
16
Elaziouti Abdelkaderr
G ère Année Master: Gestion
UMT/FT/ GP/1 G des Déchets Moodule :Dynaamique des fluides
f réelss
rendues arrtificiellem
ment rugueuuses en y collant
c des grains de sable calibbrés. En ch
hangeant laa
taille des grains
g on chhange la ruugosité et on
o trace log
g λ enfoncttion de logg Re.
On appellee ks la taille moyennne des grainns collés et D le diam
mètre de laa conduite. ε /D est laa
rugosité reelative. On remarque sur les abaaques de Nikuradze(F
N Figure II-66 ) que:
9 la rugosité
r n’aa pas d’inffluence en régime lam l nombre de Reynollds critiquee
minaire et le
Rc est indépenndant de laa rugosité.
9 quaand Re>Reecla perte de
d charge dépend
d forttement de ε /D ou du nombre dee Reynoldss
9 quaand R >>>
>Rc la pertte de charge ne dépeend pas duu nombre de Reynollds, elle nee
déppend que dee ε /D.
Figgure II-5 :E
Expériencee de Nikuraadze
17
Elaziouti Ab
bdelkaderr
UMT/FT/ GP/1ère Année Master: Gestion des Déchets Module :Dynamique des fluides réels
(37)
La formule de ColebrookEq. (38) est actuellement considérée comme celle qui traduit le
mieux lesphénomènes d'écoulement en régime turbulent. Elle est présentée sous la forme
suivante :
,
)(38)
√ , √
L'utilisation directe de cette formule demanderait, du fait de sa forme implicite, un calcul par
approximations successives ; on emploie aussi en pratique des représentations
graphiques(abaques).
Pour simplifier la relation précédente, on peut chercher à savoir si l'écoulement
esthydrauliquement lisse ou rugueux pour évaluer la prédominance des deux termes
entreparenthèses dans la relation de Colebrook.
Remarque :
On fait souvent appel à des formules empiriques plus simples valables pour des cas
particuliers et dans un certain domaine du nombre de Reynolds, par exemple :
18
Elaziouti Abdelkaderr
G ère Année Master: Gestion
UMT/FT/ GP/1 G des Déchets Moodule :Dynaamique des fluides
f réelss
, (40)
19
Elaziouti Ab
bdelkaderr
G ère Année Master: Gestion
UMT/FT/ GP/1 G des Déchets Moodule :Dynaamique des fluides
f réelss
∆ (4
41)
où K :L
Le coefficiient de perrte de chargge singulièrre (sans dim
mension). Il dépend de la
natuure de lasinngularité.
v : La
L vitesse moyenne
m i
incidente d fluide arrrivant sur l’obstacle.
du l
Remarquee :
La plupartt des singuularités ne permettennt pas une approche analytique
a simple co
omme cellee
qui vient d'être
d déveeloppée conncernant l'élargissem
ment brusquue. Les coeefficients de
d perte dee
charge sonnt donc lee plus souuvent déduuit de messures ou de
d calculs numériquees, et sonnt
disponiblees sous forrme d'abaqques(Figurre II-8). À titre d'exxemple, la figure III-9 liste dee
façon non exhaustivee quelquess singularittés typiquees avec les coefficiennts de pertee de chargee
associés.
20
Elaziouti Ab
bdelkaderr
G ère Année Master: Gestion
UMT/FT/ GP/1 G des Déchets Moodule :Dynaamique des fluides
f réelss
Figu
ure II-9 :Liiste non exxhaustive de
d quelquees singulariités typiquees
∆ ∑ ∑ (42)
21
Elaziouti Ab
bdelkaderr
Chapitre III : Les écoulements internes
G ère Année Master: Gestion
UMT/FT/ GP/1 G des Déchets Moodule :Dynaamique des fluides
f réelss
III-1 Intrroduction
Par rapporrt aux écouulements exxternes, less écoulemeents internees ont ceci de particulier que less
conditionss aux limittes y sont omnipréseentes. Ils sont
s aborddés dans cee chapitre sous deuxx
aspects complémentaaires : les profils
p de vitesse et les pertes de charge qui condiitionnent lee
dimensionnnement dees circuits fluides.
f
Les écouleements inteernes ont pour
p caracttère particu
ulier l’abseence de connditions au
ux limites à
l’infini, cee qui n’est pas le cas des écouleements exteernes. L’innfluence dees parois see manifestee
donc danss toutes less directionns, sauf évventuellemeent aux seections d’entrée et dee sortie duu
fluide, et en ons doivent
e particulier le dessiin deséchaangeurs ; plus générallement, less dissipatio
être prises en comptee dans les bilansénerg
b gétiques ou
u économiqques.
III-2Écou
ulement paarallèle dan nduitebidimensionnelle
ns une con
On considère une connduite bidiimensionneelle infiniee, d’épaisseeur e=2b faaible, l’axee x coïncidee
avec l’axee de sym
métrie de la conduiite. Le flluide est considéré comme newtonienn,
incompresssible en éccoulement stationnairre (permaneent) selon l’axe x(Figgure III-1)).
∆ (
(44)
∆ (45)
22
Elaziouti Ab
bdelkaderr
UMT/FT/ GP/1ère Année Master: Gestion des Déchets Module :Dynamique des fluides réels
On en déduit que la pression totale n'est fonction que de x, et qu'elle varie en obéissant à
l'équation Eq. (46) :
Le profil des vitesses dans la conduite est donc donné par la résolution de l’équation
différentielleEq. (47) :
′(48)
oùK et K’ sont des constantes d’intégration, elles sont déterminées grâce aux conditions
aux limites sur la vitesse ( une particule de fluide en contact avec une paroi possède la
vitesse de la paroi) soit:u(b) = u(-b) = 0
En insérant ces deux conditions dans l’équation précédente, on obtient un système de deux
équations à deux inconnues facile à résoudre, on obtient Eq. (49):
′
(49)
23
Elaziouti Abdelkaderr
G ère Année Master: Gestion
UMT/FT/ GP/1 G des Déchets Moodule :Dynaamique des fluides
f réelss
q. (48), on trouve
En reportaant dans l’ééquationEq t le profil
p des vitesses
v suivvantEq. (5
50) :
(50))
Le profil des
d vitessess est paraboolique.
∆ ∆
L est la lonngueur de la
l conduitee et P1>P2
Il suffit dee mesurer lees pressionns correspoondantes po
our déterm
miner le graddient de prression
La vitesse moyenne de
d l’écouleement est définie
d Eq. (54) :
parE
(54)
(5
55)
24
Elaziouti Ab
bdelkaderr
UMT/FT/ GP/1ère Année Master: Gestion des Déchets Module :Dynamique des fluides réels
6)
(57)
(58)
III-3Écoulement de Couette
Une autre situation simple est celle de l’écoulement de Couette(Figure III-3). On a les
mêmeshypothèses que précédemment sur la géométrie de la conduite, sur le fluide
etl’écoulement. La différence vient du fait que la paroi supérieure possède unevitesse V0.
Comme dans l’exemple précédent, on néglige les forces de pesanteur.
On obtient bien sûr la même équation différentielle avec des conditions auxlimites qui
changent.
+b
V0
y
0 x
+b
′(59)
25
Elaziouti Abdelkaderr
UMT/FT/ GP/1ère Année Master: Gestion des Déchets Module :Dynamique des fluides réels
(61)
(Poiseuille)(Couette)
(62)
(64)
Le débit total de fluideentre les deux plans est dû à la fois à la vitesse d’entraînement et au
gradient de pression. En désignant par lla largeur de la veine fluide, le débit-volume qv
s’écrit, dans une section SEq. (65):
(65)
26
Elaziouti Abdelkaderr
G ère Année Master: Gestion
UMT/FT/ GP/1 G des Déchets Moodule :Dynaamique des fluides
f réelss
III-4Écou
ulement dee Poiseuillee
(68))
ng de l'axe de la condduite.
traduisant le fait que la vitesse reste consttante le lon
27
Elaziouti Ab
bdelkaderr
UMT/FT/ GP/1ère Année Master: Gestion des Déchets Module :Dynamique des fluides réels
ce qui permet d'en conclure que le vecteur vitesse est purement axial et ne dépend que de
rEq. (70):
(70)
D’autre part, l'équation régissant l'évolution de la pression totale (charge) Pt le long d'un
écoulement laminaire orienté suivant l'axex (voir section précédente) Eq. (71):
∆ (71)
∆ (72)
alorsEq. (73):
(73)
Résolvons cette dernière équation pour établir le profil de vitesse sur une section de la
conduite Eq. (74):
′
′
et donc :
′ "
′
où " sont deux constantes d'intégration dont l'indétermination peut être levée en
utilisant les conditions aux limites.
Ainsi, on sait que la vitesse des particules fluides est nécessairement nulle au contact avec les
parois de la conduite. On a doncEq. (75) :
′ " (75)
Par ailleurs, la vitesse à l'axe de la conduite, enr=o, soit finie, il faut que K’=0alorsEq. (76):
28
Elaziouti Abdelkaderr
G ère Année Master: Gestion
UMT/FT/ GP/1 G des Déchets Moodule :Dynaamique des fluides
f réelss
"
Cette équuation tradduisant lepprofil parabbolique dee u (r)(Fiigure III--5) et que la vitessee
maximale sur l'axe estEq. (78) :
d charges, on défin
Pour caracctériser les pertes de nit la viteesse moyeenne de l'éécoulement
parEqs. (779-80) :
(79)
(80
0)
(81)
et donc Eq
q. (82):
(82)
Elaziouti Ab
bdelkaderr
G ère Année Master: Gestion
UMT/FT/ GP/1 G des Déchets Moodule :Dynaamique des fluides
f réelss
(83)
∆ ∆ ∆ (84)
La différennce de pression totale (ou pertee de chargee) engendrrée par les frottementts visqueuxx
dans une conduite
c dee longueur L=∆xest donc
d ortionnelle à L, mais également au débit et
propo
à la viscossité du fluiddeEq. (85) :
∆ (85))
III-5Écou
ulement daans un espaace annulaaire
30
Elaziouti Ab
bdelkaderr
UMT/FT/ GP/1ère Année Master: Gestion des Déchets Module :Dynamique des fluides réels
∆ (86)
(87)
(88)
(89)
(90)
(91)
En ce qui concerne le débit dû au gradient dPt/dx, les calculs sont un peu pluslaborieux.
Posons provisoirement, pour alléger l’écriture Eq. (93):
Au total, on a doncEq.(95) :
(95)
31
Elaziouti Abdelkaderr
UMT/FT/ GP/1ère Année Master: Gestion des Déchets Module :Dynamique des fluides réels
Au total, après avoir remplacé C par (Eq.(93)), puis groupé et réarrangé un peu les termes, on
a doncEq.(96) :
32
Elaziouti Abdelkaderr
Chapitre IV : Analyse dimensionnelle et
similitude
G ère Année Master: Gestion
UMT/FT/ GP/1 G des Déchets Moodule :Dynaamique des fluides
f réelss
IV-1 Introoduction
D: Diamèètre de la coonduite,
v :Vitessse moyennee del’écoullement(ou Q: Débit volumique)
v ),
ρ :Massee volumiquue du fluidee,
μ : Viscoosité dynam
mique du fluide,
fl
ε :Rugositté moyenne de la connduite.
Par conséqquent, il exxiste une reelation entrre ces différrentes grandeurs Eq. (97):
33
Elaziouti Ab
bdelkaderr
UMT/FT/ GP/1ère Année Master: Gestion des Déchets Module :Dynamique des fluides réels
∆
=f(D,v,ρ,μ,ε)(97)
(nombre de Reynolds)(98)
(rugosité relative)
,
∆ ∆
, , 99
L’analyse dimensionnelle permet aussi de voir que la perte de charge régulière est fonction
uniquement de nombre de Reynolds et de la rugosité relative de la conduite.
Voyons maintenant une description détaillée de la méthode.
IV-2 Théorème de Buckingham( Théorème de )
Si une équation comportant k variables est homogène, elle peut être réduite à une relation
entre (k-r) produit indépendant sans dimension, ou r est le nombre minimal de dimension
requis pour décrire les k variables.
Inventaire des variables
Vitesse v LT-1
Diamètre D L
Longueur L L
Rugosité ε L
34
Elaziouti Abdelkaderr
UMT/FT/ GP/1ère Année Master: Gestion des Déchets Module :Dynamique des fluides réels
Viscosité m ML-1T-1
Masse volumique ρ ML-3
Afin d’illustrer cet énoncée, reprenant l’exemple d’introduction.
∆
Nous avions k=6 variables ( , D,v,ρ,μ,ε),qui nécessitent un minimum de r=3 (trois
(100)
Pour effectuer une analyse dimensionnelle, on doit considérer les étapes suivantes:
1- Dresser la liste de toutes les variables du problème k,
2- Ecrire l’équation aux dimensions de chacun des k variable,
3- Déterminerr, et donc k-r le nombre de produit sans dimension caractérisant le
problème,
4- Parmi les k variables,on choisir un nombre r, qui soient dimensionnellement
indépendant r variables primaires,
5- Former les k-r produits (π) en combinaison les k-rvariable non primaires avec les r
primaires.
6- Formuler la relation entre les k-r produits (π) trouvé.
35
Elaziouti Abdelkaderr
G ère Année Master: Gestion
UMT/FT/ GP/1 G des Déchets Moodule :Dynaamique des fluides
f réelss
4. Chois de
d r=3 varriable prim
maire dimennsionnellem
ment indéppendantes. P
Par exemp
ple h, v, ρ.
Remarrque :On ne
n peut pass choisir à la
l fois h et L
36
Elaziouti Ab
bdelkaderr
UMT/FT/ GP/1ère Année Master: Gestion des Déchets Module :Dynamique des fluides réels
(102)
→ sans dimension
M0L0T0= M1L1T−2Lα1(M1L-3)β1(L1T−1)γ 1= M1+β1L(1+α1−3β1+γ 1) T−2−γ 1
0 1+β1 β1=-1
0=1+α1−3β1+γ1 α1=-2
0=−2−γ 1γ 1=-2
→ sans dimension
M0L0T0= L1Lα2(M1L-3)β2(L1T−1)γ 2= Mβ2L(1+α2 −3β2 +γ 2) T−γ 2
0 β2β2=0
0=1+α2−3β2+γ2 α2=-1
0=−2γ 1 γ2=0
→ sans dimension
∆
= f(D, v, ρ, μ, ε) ,
Avec
37
Elaziouti Abdelkaderr
UMT/FT/ GP/1ère Année Master: Gestion des Déchets Module :Dynamique des fluides réels
, ù (103)
= si
et donc , = ,
En mécanique des fluides, on est souvent amène a manipuler des groupes de variables sans
dimension, appelés nombre adimensionnel ou rapport de similitude . Ces groupes
sont construits en faisant des rapports entre des termes apparaissant dans les équations du
mouvement, ce qui permet de les interpréter physiquement. On distingue ainsi :
Nombre deReynolds
’
(104)
é
avecLune échelle de longueur, vune échelle de vitesse,μla viscosité du fluide, et ρsa masse
volumique. Il exprime le rapport entre les forces d'inertie et les forces de viscositéEq. (104).
38
Elaziouti Abdelkaderr
UMT/FT/ GP/1ère Année Master: Gestion des Déchets Module :Dynamique des fluides réels
Nombre deFroude
’
(105)
é
Remarque : c’est la pesanteur qui est responsable de la forme de lasurface libre : plus Fr est
grand, moins la surface libre a d’effets surl’écoulement, et inversement.
Nombre d’Euler
∆
(106)
′
Nombre de Mach
Le nombre de Machexprime le rapport entre la vitesse caractéristique de l'¶écoulement et la vitesse
de propagation du son caractéristique du fluide. Ce nombre quantité donc les effets de
compressibilitéEq. (107) :
′
(107)
39
Elaziouti Abdelkaderr
UMT/FT/ GP/1ère Année Master: Gestion des Déchets Module :Dynamique des fluides réels
avecuune échelle de vitesse et C la célérité du son (ou célérité des ondes dans l’air). Le
nombre de Mach sert en aérodynamique ç évaluer la compressibilité de l’air. On distingue
ainsi les écoulements supersoniques (M >1) et subsoniques (M <1) ;
Remarque : est la vitesse du son.
Nombre deStrouhal
(108)
Nombre de Stokes
avectple temps de relaxation de la particule (le temps typique de variation de la vitesse quand
on perturbe l’état d’équilibre de la particule) et le temps caractéristique du fluide (l’échelle de
temps sur laquelle le fluide s’ajuste a tout changement de la particule). Ce nombreEq. (109).
sert dans l’étude des écoulements biphasiques (par exemple, une suspension de particules) à
quantifier les effets biphasiques, c’est-a-dire le couplage entre phases.
Lorsque St 1, la phase solide est entièrementgouvernée par la phase fluide tandis que pour
St 1, les deux phases sont découplées.
Nombre de Prandtl
Il exprime le rapport entre la viscosité cinématique et la diffusivité thermique du fluide. Il s'agit
donc d'un nombre caractérisant le fluide. Le nombre de Prandtl Eq. (110). est donc indépendant de
l'¶écoulement.
40
Elaziouti Abdelkaderr
UMT/FT/ GP/1ère Année Master: Gestion des Déchets Module :Dynamique des fluides réels
Avant de procéder à l’analyse complète d’un écoulement, il convienttout d’abord de poser les
hypothèses simplificatrices adéquates.L’évaluation des différents coefficients sans dimension
relatifs àl’écoulement (Reynolds, Froude…) va en effet permettre de simplifierles équations
à résoudre.
w*=w/V z*=z/L
41
Elaziouti Abdelkaderr
UMT/FT/ GP/1ère Année Master: Gestion des Déchets Module :Dynamique des fluides réels
et
avec
42
Elaziouti Abdelkaderr
UMT/FT/ GP/1ère Année Master: Gestion des Déchets Module :Dynamique des fluides réels
43
Elaziouti Abdelkaderr
UMT/FT/ GP/1ère Année Master: Gestion des Déchets Module :Dynamique des fluides réels
α1/α
En pratique, tout n’est pas si simple, dans la mesure où l’on ne peutsouvent pas respecter
simultanément toutes les similitudes.
44
Elaziouti Abdelkaderr
Annexe
G ère Année Master: Gestion
UMT/FT/ GP/1 G des Déchets Moodule :Dynaamique des fluides
f réelss
NB :
• Les unités du système internnational SI, unités dee base, dim
mensions, rèègles d’écrriture.
45
Elaziouti Ab
bdelkaderr
G ère Année Master: Gestion
UMT/FT/ GP/1 G des Déchets Moodule :Dynaamique des fluides
f réelss
B-Alph
habet grecc
C- Préfixees du SI
46
Elaziouti Ab
bdelkaderr
UMT/FT/ GP/1ère Année Master: Gestion des Déchets Module :Dynamique des fluides réels
Bibliographies
47
Elaziouti Abdelkaderr