Polycopie MDF App Bonne-Version Preliminaire

1
Ministère de l’Enseignement Supérieur et de la Recherche Scientifique
Université de Batna 2 - Mustapha Benboulaı̈d
Faculté de Technologie
Département de Génie Mécanique
Support de cours
Destiné au parcours Master ERM / EN
Mécanique des fluides approfondie

partie I
Dr. Nabil BESSANANE
Année Universitaire : 2019 - 2020

2
Table des matières
1 Rappel d’analyse vectoriel 8
1.1 Les systèmes de coordonnées courants . . . . . . . . . . . . . . . . . . . . . 8
1.2 Champs scalaire et vectoriel . . . . . . . . . . . . . . . . . . . . . . . . . . 9
1.2.1 Champ scalaire : . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
1.2.2 Champ vectoriel : . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
1.3 Comment utiliser Nabla ? . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
1.4 La signification du Gradient . . . . . . . . . . . . . . . . . . . . . . . . . . 11
1.5 La signification de la divergence . . . . . . . . . . . . . . . . . . . . . . . . 11
1.6 La signification du rôtationnel . . . . . . . . . . . . . . . . . . . . . . . . . 12
1.7 La signification du Laplacien . . . . . . . . . . . . . . . . . . . . . . . . . . 13
1.8 Quelques relations importantes . . . . . . . . . . . . . . . . . . . . . . . . 14
1.9 Applications : . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
2 La physique des fluides 16
2.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
2.1.1 Échelle mésoscopique . . . . . . . . . . . . . . . . . . . . . . . . . . 17
2.1.2 Champ des vitesses dans un fluide . . . . . . . . . . . . . . . . . . 17
2.2 Équation fondamentale de l’hydrostatique . . . . . . . . . . . . . . . . . . 18
2.2.1 Forces de pression d’un fluide en repos . . . . . . . . . . . . . . . . 19
2.2.1.1 Cas de contrainte normale . . . . . . . . . . . . . . . . . . 19
2.2.1.2 Cas de contrainte tangentielle . . . . . . . . . . . . . . . 21
3
TABLE DES MATIÈRES 4
2.2.1.3 Variation de force de pression . . . . . . . . . . . . . . . . 22
2.2.2 Forces de pression d’un fluide en mouvement (uniforme) . . . . . . 23
2.2.3 Forces de pression d’un fluide en mouvement (quelconque) . . . . . 24
2.2.4 Application de de l’équation de BERNOULLI (approche statique) . 27
2.2.4.1 Effet de gravitation sur un liquide . . . . . . . . . . . . . 27
2.3 Théorèmes de transport . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
2.3.1 Volumes et surfaces de contrôle . . . . . . . . . . . . . . . . . . . . 31
2.3.2 Théorème de transport . . . . . . . . . . . . . . . . . . . . . . . . 31
3 Équations d’Euler (fluides parfaits) 33
3.1 Rappel de cinématique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
3.1.1 Approche de Lagrange . . . . . . . . . . . . . . . . . . . . . . . . . 33
3.1.2 Approche de Euler . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
3.2 Dynamique des fluides parfaits . . . . . . . . . . . . . . . . . . . . . . . . 35
3.2.1 Notion de dérivée particulaire : . . . . . . . . . . . . . . . . . . . . 35
3.2.2 Équation d’Euler : . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
3.2.3 Équation de BERNOULLI : ( autres approches) . . . . . . . . . . 37
3.2.4 Écoulement irrôtationnel : . . . . . . . . . . . . . . . . . . . . . . . 39
4 Équation de continuité (conservation de matière) 41
4.1 Courant de matière . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
4.2 Bilan de matière (conservation de masse) . . . . . . . . . . . . . . . . . . . 43
4.3 Flux de matière d’un volume fini "V" . . . . . . . . . . . . . . . . . . . . 44
4.4 Équation de continuité . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
5 Équations de Navier-Stokes (fluides réels) 47
5.1 Généralités et définitions . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
5.2 Régime d’écoulement . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
5.3 Tenseur de contraintes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48

5.4 Équation de Navier-Stokes . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
5.4.1 Notion du tenseur de déformation . . . . . . . . . . . . . . . . . . . 49
5.4.2 Décomposition du tenseur de déformation . . . . . . . . . . . . . . 50
5.4.3 Équation de Navier-Stokes . . . . . . . . . . . . . . . . . . . . . . . 51
5.5 Pertes des charges (application aux fluides réels) . . . . . . . . . . . . . . . 54
5.6 Pertes de charge linéaires . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
5.6.1 Application (mesure de la viscosité) : . . . . . . . . . . . . . . . . . 55
5.7 Pertes de charge singulières . . . . . . . . . . . . . . . . . . . . . . . . . . 56
5.7.1 Coefficient de contraction : . . . . . . . . . . . . . . . . . . . . . . 57
5.7.2 Diagramme de Moody . . . . . . . . . . . . . . . . . . . . . . . . . 58
6 Analyse dimensionnelle et similitude 63
6.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63
6.2 Procédure de l’analyse dimensionnelle . . . . . . . . . . . . . . . . . . . . 63
6.3 Théorème de Vaschy-Buckingham . . . . . . . . . . . . . . . . . . . . . . . 64
6.3.1 Exemple d’application : . . . . . . . . . . . . . . . . . . . . . . . . 65
6.4 La signification du rôtationnel . . . . . . . . . . . . . . . . . . . . . . . . . 66
6.5 Similitude . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
6.5.1 Similitude géométrique . . . . . . . . . . . . . . . . . . . . . . . . 68
6.5.2 Similitude cinématique . . . . . . . . . . . . . . . . . . . . . . . . 68
6.5.3 Similitude dynamique . . . . . . . . . . . . . . . . . . . . . . . . . 69
6.6 Régimes d’écoulement et paramètres de similitude . . . . . . . . . . . . . 70
6.6.1 Nombre de Reynolds . . . . . . . . . . . . . . . . . . . . . . . . . . 70

Avant Propos
Le contenu de la première partie de cours MDF Approfondie est destiné aux étudiants
de spécialités énergétique et énergies renouvelables en mécanique. Il s’agit de la matière
de la mécanique des fluides appliquée, il regroupe quelques rappels et notions de bases,
suivi de détailles plus approfondies en ce qui concerne le développement des équations de
bilans ; conservation de masse et de quantités de mouvement. Où l’étudiant peut connaitre
l’établissement des équations essentielles sur les quelles est fondée la mécaniques des
fluides ; tout en partant du 2nd principe de Newton avec les hypothèses adéquates pour
déterminer l’équation de la continuité, les équations d’Euler, les équations de Navier-
Stokes, l’équation de Bernoulli généralisée et voir aussi ses différentes formes appliquées
selon le cas à étudier.
Ce polycopie cours est réparti à six chapitres et est présenté d’une manière méthodologique
et simplifié du point de vue pédagogique. Le plan de ce cours est présenté comme suit :
1er chapitre : rappel sur les opérateurs mathématiques et des notions de l’analyse vectoriel
pour la mécanique des fluides.
2ème chapitre : introduction à la physique des fluides.
3ème chapitre : (équation de Euler) rappel cinématique et approche Euler

Bernoulli pour
la dynamique des fluides parfaits.
4ème chapitre : (équation de continuité), détermination de l’équation de bilan de matière

(conservation de masse)
5ème chapitre : (équations de Navier-Stokes), détermination des équations de moments,

Notion de tenseur de déformation et de contraintes (conservation de quantité de mouve-
ment) et applications des pertes de charges dans le calculs de conduites (fluides réels et
Bernoulli généralisé).
Remerciements : pour les relecteurs : Dr L.Messaoudi, Pr M.Si–Ameur, Pr M. Aksas,

Pr A.H.Benmachiche, Dr G.Mebarki, Dr A.Benderradji.
Ce polycopie de cours n’est plus soumis à aucun droit et l’auteur a renoncé entièrement
pour que son contenu soit largement diffusées et que le bénéfice puisse se répandre et
l’étudiant de spécialité en profite autant que possible !
La gestion typographique du texte français a été réalisée avec Lyx et LATEX à l’aide du
package french.sty de Bernard Gaulle.
Dans l’attente de l’achèvement d’une deuxième partie de ce cours (MDF approfondie) qui
suivra dans le futur et qui va regrouper les volets suivants :
- Similitudes et analyse adimensionnel (adimensionnement).
- Écoulements cisaillés et théorie de la couche limite.
- Écoulements turbulents et introduction aux modèles de turbulence.

Chapitre 1
Rappel d’analyse vectoriel
Ce chapitre présente l’essentiel des notions mathématiques portant sur les opérateurs
vectoriels et tensoriels et la notation indicielle, qui est largement utilisée dans le cours de
Mécanique des Fluides.
1.1 Les systèmes de coordonnées courants
Pour repérer la position d’un point M dans l’espace à trois dimensions, il est nécessaire
d’introduire trois axes non coplanaires. Le repérage peut alors être réalisé par l’introduc-
tion de (voir figure ci-dessous) :
- trois distances : coordonnées cartésiennes,
- deux distances et un angle : coordonnées cylindriques,
- une distance et deux angles : coordonnées sphériques.
8
CHAPITRE 1. RAPPEL D’ANALYSE VECTORIEL 9
1.2 Champs scalaire et vectoriel
Un champ est un outil physique qui donne, pour un point de l’espace, une valeur d’une
grandeur physique. Autrement dit, le champ établit une correspondance entre une position
de l’espace et une valeur prise par la grandeur physique étudiée.
Pour commencer, nous allons définir un champ scalaire ou vectoriel à partir de notions que
tout étudiant connait. Soit un trièdre orthonormé (→−e x, →
−
e y, →
−
e z ) et M un point de l’espace,
−−→
de coordonnées (x, y, z) ; d’où le vecteur position de M s’écrit : OM = X → −e x +Y →
−
e y +Z →
−ez
1.2.1 Champ scalaire :
Un champ scalaire donne la valeur d’une grandeur physique scalaire en fonction de la po-
sition, pour un temps donné. exemple de : température, densité, pression...etc. La fonction
f(M) est dite fonction scalaire de point M ou champ scalaire si :f (M ) = f (x, y, z).
On parle de champ scalaire lorsque la grandeur physique est un nombre (réel). La température
et la pression d’une zone sont décrits par des champs scalaires. Il existe différentes
manières de représenter un champ scalaire selon son application : coloriage de contours,
équipotentielles, graphe 3D...
1.2.2 Champ vectoriel :
Pour un moment donné, un champ vectoriel, ou champ de vecteurs, donne la valeur d’une
grandeur vectorielle en fonction de la position M. Autrement dit, chaque point de l’espace
est mis en correspondance avec un vecteur. Comme pour les champs scalaires, un champ
vectoriel peut être étudié dans un espace 2D, ou 3D et dans le temps. Par exemples, le
champ de la pesanteur, le champ magnétique, en mécanique des fluides le mouvement de
gaz ou de liquides est étudié en utilisant des champs de vitesse.
→
−
Le vecteur vitesse V (M ) est dit fonction vectorielle de point M ou champ vectoriel si :
→
−
V (M ) = →
−u (x, y, z)→
−
ex+→−v (x, y, z)→
−
ey+→ −
w (x, y, z)→
−
e z.
Un champ vectoriel établit un lien entre une position de l’espace est une grandeur physique
vectorielle. Les champs de vitesse en sont un exemple. Le champ de pesanteur est un
champ vectoriel uniforme localement. Les champs électriques et magnétiques sont d’autres
exemples de champ vectoriels.
1.3 Comment utiliser Nabla ?
On l’appelle Nabla ou aussi l’opérateur d’Hamilton chez les mathématiciens, les opérateurs
mathématiques (Gradient, Divergence et rôtationnel) peuvent être appliqués mathématiquement
et utilisés sur des champs scalaires et vectoriels à travers l’opérateur différentiel Nabla,
pour retrouver : le Gradient, la divergence, le rôtationnel, Laplacien scalaire et Laplacien
vectoriel.
Pour présenter le concept très utile de l’opérateur Nabla, il faut insister de connaitre
le principe du champ scalaire et champ vectoriel (définis précédemment), pour ne plus
confondre entre les deux notions intuitives.
— Champ scalaire : soit une  fonction ϕ, en coordonnées cartésiennes son champ
ϕ : R3 −→ R



scalaire est définie par :

 M (x, y, z) −→ ϕ(x, y, z) = ϕ(M )

exemples : température, pression ; énergie potentiel...etc.
→
−
— Champ vectoriel : soit un vecteur A (M ), mathématiquement on défini le champ
→
−
A : R3 −→ R3





 



 
 A (x , y , z )
 
 x 1 1 1
vectoriel de ce vecteur par :


→
− →
− → 
r (x, y, z) −→ A (−r)=




  Ay (x2 , y2 , z2 )

 


 


  Az (x3 , y3 , z3 )

c’est dans le cas d’un état stationnaire ; c’est à dire, qui ne dépend pas de t (temps).
exemples : champ des vitesses, champ des forces, champ magnétique, champ électrique
...etc.
Donc, pour trouver les notions de gradient, divergence et rôtationnel, on va utiliser

l’opérateur différentiel Nabla. C’est le symbole triangle tète en bas avec flèche ∇ ,
on le note formellement dans le système de coordonnées cartésiennes par :


 ∂

 ∂x
→
− 
∇= ∂
∂y



 ∂

∂z
Remarque : Il est alors évident que pour avoir les composantes de cet opérateur dans
les différents systèmes de coordonnées, il suffit de reprendre celles du gradient (ci-après)
et d’enlever la fonction ϕ .
1.4 La signification du Gradient

−−→
— Le gradient : grad est une quantité vectoriel et est appliqué seulement pour
un scalaire et non pas sur un vecteur. on dit le gradient d’un champ scalaire.
 
∂ϕ
 ∂x 
−−→ →
− 
∂ϕ

pour trouver le gradient d’une fonction ϕ, on écrit : gradϕ = ∇ϕ=
 ∂y
 ; comme un

 
∂ϕ
∂z
produit scalaire entre Nabla et la fonction scalaire ϕ.
Il faut noté ! Le gradient donne donc des informations sur la direction, sur le sens de la
variation de la fonction f mais aussi sur l’importance de cette évolution.
1.5 La signification de la divergence

→
−
— La divergence : div A est une quantité  scalaire
 et ellemême appliquée à un
∂
 ∂x   Ax 
→
− →
− →
− → −  ∂   •  A  = ∂Ax + ∂Ay + ∂Az ;
 
vecteur A , on écrit : div A = ∇ • A = 
 ∂   y  ∂x ∂x ∂x
   
∂
∂z
Az
→
−
comme un produit scalaire entre Nabla et le vecteur A .
L’expression de la divergence dans les différents systèmes de coordonnées est présentée

dans la suite :
→
− ∂A
Coordonnées cartésiennes : div A = ∂A y ∂Az
∂x + ∂x + ∂x
x
→
− ∂(rAr ) 1 ∂Aθ ∂Az
Coordonnées cylindriques :div A = 1r ∂r + r ∂θ + ∂x
→
− ∂ (r2 Ar ) 1 ∂(sinθAθ ) 1 ∂Aϕ
Coordonnées sphériques :div A = r12 ∂r + rsinθ ∂θ + rsinθ ∂ϕ
→
−
Remarque : en coordonnées cartésiennes, il faut utiliser ∇ . Avec les coordonnées
→
−
sphériques, compte tenu des multiples implications de ∇ sur la base, il sera préférable
d’utiliser le théorème d’Ostrogradski-Green (voir section 3.3). Dans le cas des coordonnées
cylindriques, les deux possibilités sont à peu près équivalentes.
1.6 La signification du rôtationnel
— Le rôtationnel : est une quantité vectoriel et est appliqué sur un vecteur

→
− −→→−
A (très utile
 en 
mécanique
 des fluides
  et électromagnétisme),
 on écrit : rot A =
∂ ∂Az ∂Ay
   Ax   ∂y − ∂z 
∂x
→
− →
−      
∇∧A =

∂ ∧
∂y   A y
 =  − ∂Az + ∂Ax 
  ∂x ∂z 
     
∂ ∂Ay ∂Ax
∂z
Az − ∂y
 ∂x 
∂Az ∂Ay
 ∂y − ∂z 
−→ →
− 
∂Ax ∂Az 

— Coordonnées cartésiennes :rot A = 
 ∂z − ∂x 
 
∂Ay ∂Ax
− ∂y
 ∂x 
1 ∂Az ∂(rAθ )
 r ∂θ − ∂z 
−→→− 
∂A ∂A

— Coordonnées cylindriques :rot A = 
∂z
r
− ∂r
z 

 
1 ∂(rAθ )
r ∂r
− ∂A ∂θ
r
 
1 ∂(sinθAϕ ) ∂Aθ
 rsinθ ∂θ
− ∂ϕ 
−→ →
− 
1

1 ∂Ar ∂(rAϕ )

— Coordonnées sphériques :rot A = 
 r sinθ ∂ϕ − ∂r


 
1 ∂(rAθ ) ∂Ar
r ∂r
− ∂θ
Remarque : une remarque fondamentale consiste à dire que le caractère rotationnel
n’est pas lié à la courbure des lignes de champ. Il ne suffit pas que ces lignes de champ
soient rectilignes pour que l’écoulement soit irrotationnel bien au contraire ! Cet opérateur
nous donne donc des informations sur le caractère localement tourbillonnaire du champ
vectoriel.
1.7 La signification du Laplacien
— Le Laplacien : il faut faire attention ! Il y a le Laplacien scalaire et le

Laplacien vectoriel. il est noté par le symbole triangle tête en haut 4 .
1. Le Laplacien 
scalaire
 :Le Laplacien
  d’une fonction scalaire ϕ est : 4ϕ =
∂ϕ ∂ ∂ϕ
 ∂x   ∂x   ∂x 
−−→ 
∂ϕ
 
∂
 
• ∂ϕ
 ∂ ∂ϕ
∂

∂ϕ

∂ ∂ϕ

div(gradϕ) = div 
 ∂y
=
  ∂   ∂y
=
 ∂x ∂x
+ ∂y ∂y ∂z ∂z
=
     
∂ϕ ∂ ∂ϕ
∂z ∂z ∂z
∂2ϕ ∂2ϕ ∂2ϕ ∂2ϕ ∂2ϕ ∂2ϕ →
− →− →−
∂x2
+ ∂y 2
+ ∂z 2
=⇒ 4ϕ = ∂x2
+ ∂y 2
+ ∂z 2
= ∇ • ∇ϕ = ∇ 2 ϕ = 4ϕ
→
−
2. Le Laplacien  vectoriel
 : Le Laplacien d’un champ vectoriel (vecteur A )
 4Ax 
→
−→−  
est : 4 A =   4Ay , à noter que Ax , Ay etAz ce sont des champs scalaires. il

 
4Az
est bien appliqué en mécanique des fluides dans le cas des calculs des équations de
Navier-Stokes.
−−→

∂ 2 Ax ∂ 2 Ax ∂ 2 Ax
4Ax = div(gradAx ) = + +



 ∂x2 ∂y 2 ∂z 2
−−→

∂ 2 Ay ∂ 2 Ay ∂ 2 Ay
donc, on écrit : 4Ay = div(gradAy ) = ∂x2
+ ∂y 2
+ ∂z 2

 4Az = div(−−→

∂ 2 Az ∂ 2 Az ∂ 2 Az

gradAz ) = + +

∂x2 ∂y 2 ∂z 2
→
−
on obtient
 à la fin le Laplaciend’un vecteur A , qui s’écrit sous cette forme de vecteur :
∂ 2 Ax ∂ 2 Ax ∂ 2 Ax
 ∂x2 + ∂y2 + ∂z2 
→
−→−  2
∂ Ay ∂ 2 Ay ∂ 2 Ay 

4A =   ∂x2 + ∂y2 + ∂z2 
 
∂ 2 Az ∂ 2 Az ∂ 2 Az
∂x2
+ ∂y2 + ∂z2
Remarque : la majorité des étudiants se trompent dans l’écriture du Laplacien de
 
∂ 2 Ax
 ∂x2 
→
− →
−→− 
∂ 2 Ay

vecteur A , ou il mette : 4 A = 
 ∂y 2
 ; qui est évidemment faut !

 
∂ 2 Az
∂z 2
Remarque très importante : il faut bien noter qu’on ne peut pas utiliser un pseudo
opérateur Nabla ”∇” en coordonnées cylindriques, cylindropolaires ou sphériques. Atten-
tion ! Ça ne marche pas et ça ne fonctionne que en coordonnées cartésiennes, car c’est le
seule systèmes de coordonnées avec lequel on peut trouver ces opérateurs sus-cités.
1.8 Quelques relations importantes

−−→
— div rotA = 0 : ceci est évident en écrivant cette relation avec l’opérateur nabla :
−−→ → − → − → −
div rotA = ∇ • ∇ ∧ A =0. Cette égalité implique que si la divergence d’un
→
− →
− →
− −−→
champ vectoriel B est égale à 0, nous pourrons mettre B sous la forme B = rotA
.
−−−−−−
−−
→−→ → −
— rot(gradf ) = 0 : en utilisant l’opérateur nabla, cette égalité est elle aussi évidente :
−−−−−−
−−
→−→ →
− → −
rot(gradf ) = ∇Λ ∇f . En mécanique des fluides, si le rotationnel du champ des
vitesses est nul, on parle d’écoulement irrotationnel et l’on introduit la fonction
→
− −−−−−−→
scalaire potentielle des vitesses ϕ définie par V = grad(−ϕ) .
1.9 Applications :
Exercice 1 : Soit deux points M et P de coordonnées respectives M(x,y,z) et P (xP , yP , zP )

.
−−→ −−→
1. Calculer →
−
r = P M et r = P M .

→
−
2. Calculer ∇( 1r ) au voisinage de M.
→
− −
→
3. Calculer ∇ × ( rr3 ) au voisinage de M.
→
− −
→
4. Calculer ∇ • ( rr3 ) au voisinage de M.
Réponses :
−−→ −
1. Le vecteur position P M :→
r = (x − xP )→
−
e x + (y − yP )→
−
e y + (z − zP )→
−
ez
−−→ −
et le module du vecteur position P M :→
r = (x − xP )→
−
e x + (y − yP )→
−
e y + (z − zP )→
−
ez
→
− (x−x )→
−
e +(y−y )→
−
e +(z−z )→
−
e →
−
2. ∇ 1r = − = − rr3
P x P y P z
2 2 2 3/2
[(x−xP ) +(y−yP ) +(z−zP ) ]
3. En tenant compte du résultat précédent et sachant que le rotationnel d’un gradient est
nul, on obtient le produit vectoriel :
→
− −→
→−
∇ × rr3 = 0
4. Le produit scalaire est :
→
− −→

∇ • rr3 = 0
Exercice 2 :
Soit :f (x, y) = x2 + y 3
1. Calculer la divergence de f : (∇f )
2. Calculer Laplacien de f : (4f )

Chapitre 2
La physique des fluides
On appelle Mécanique l’étude des déplacements et des déformations des corps au cours du
temps, y compris l’étude des conditions qui entrainent ces mouvements. Nous considèrerons
ici la Mécanique au sens restreint où n’interviennent ni changements d’état physique, ni
transformations chimiques (vaporisation, cavitation, combustion. . .). Ce chapitre contient
des rappels de base en mécanique des fluides, liés à la physique des fluides et ses concepts
fondamentaux.
2.1 Introduction
La dynamique est la partie de la Mécanique des fluides qui étudie (sans expliciter la
variable température T) les mouvements ou le repos dans leurs rapports avec les forces
qui les engendrent. La cinématique fournit le cadre spatiotemporel dans lequel sont décrits
les mouvements dans l’espace euclidien à 3 dimensions. La cinétique se construit à partir
de la cinématique en introduisant la notion de masse.
La Mécanique des fluides (MDF) étudie la physique des fluides, afin de décrire un fluide
en repos ou en mouvement.
Dans la physique du solide : on s’intéresse à la distance : d (M1 M2 )= Cste

→
− −→ → − −→
=⇒ V A .AB= V B .AB
on peut tirer un champ de vitesse pour solide, c’est la projectivité des vitesses (translation,
rotation ou les contrainte).
16
CHAPITRE 2. LA PHYSIQUE DES FLUIDES 17
→
− →
− →
− −→
où : V B = V A + Ω ∧AB.
Pour un fluide : les molécules peuvent se déplacer librement , les unes par rapport aux
autres, et donc la distance d(M1 M2 ) peut varier à n’importe quel instant.
On ne peut pas suivre les particules une par une (impossible !), en thermodynamique
(1 mole → N =6.1023 ) d’un système en équilibre - dit volume de contrôle (VC) - est
décrit par des grandeurs moyennes intensives “ρ, x, P, T ” et extensives “volume, masse”.
Cela n’est pas suffisant dans le cas de l’hydrodynamique des fluides, car elles sont non
homogènes sur l’ensemble du VC, et les différentes propriétés du fluide ne sont plus les
mêmes dans les différents points du volume de contrôle.
2.1.1 Échelle mésoscopique
Entre l’échelle microscopique et macroscopique, il y a l’échelle mésoscopique (échelle ponc-

tuelle). C’est une échelle située entre les deux échelles micro et macro.
Définition : l’échelle mésoscopique ponctuel est une échelle inférieure aux aux échelles
sur les quelles les grandeurs moyennes varient suffisamment grand pour décrire l’échelle
moyen et sa variation est très petite d’un point matériel à l’autre.
2.1.2 Champ des vitesses dans un fluide
Pour un petit élément de fluide à l’échelle mésoscopique, le champ des vitesses en un point
M est présenté par : →
−
v (M ) = h→
−
v i i.
En chaque point M, il y a une vitesse spécifique →

−
v i lié à ce dernier.
Soit deux points matériels M et M’ de masse m et de volume élémentaire V suffisamment

petit (V=0).
En peut définir des champs des différent paramètres :
— Champ des masses volumiques, ρ(M ) = dm

dV
— Champ des températures T(M), d’après les lois de la statistique à l’équilibre d’un
système, la description des particules ne change pas .
— Champ des vitesses → −v (M ) : la distribution de la vitesse obéie à la probabi-
1
2
mvx2
−
lité :P (vx ) = K exp kB T
x2
− 2σ
c’est une Gaussienne de forme : e 2
, la densité de distribution de vitesse est unique
à l’équilibre et la seule différence est la largeur de la gaussienne 12 mv 2 i=h 23 kB T . La

température est le seule et unique paramètre qui détermine la vitesse.
Champ de pression : P= 31 nmv 2 ou P = nkT , c’est la pression exercée par un fluide sur
une paroi immergée dans un fluide.
→
− →
−
La 2nd loi de Newton : m ddt
vi
= f i , à un instant t il y a des milliard de molécule
qui viennent frapper la paroi.
−→
On a : dF = P.dS.→
−
n ⊥ , avec : P = nkB T et n = N
V
= nombre de particules
volume
où : P V = N kB T =⇒ P V = nNA kB T =⇒ P V = nRT .
Une fois qu’on a compris qu’il y a une distribution de vitesse pour un système à l’équilibre,
caractérisée par une grandeur de T et avec le calcul de la pression, on trouve P V = nRT .
Cela est à l’échelle moléculaire (càd microscopique).
Lorsqu’on a un large système, on n’est pas obligé de suivre chaque molécule, car dès que
le système est suffisamment vaste (largement grand) et à l’équilibre, alors la distribution
des vitesses est déterminée de la distribution dite de Maxwell Boltzmann.
P→− →
−
Pression en un point du fluide à l’équilibre ⇒ F ext = 0 (c’est l’état statique du fluide).
2.2 Équation fondamentale de l’hydrostatique

→
− →
− →
−
En général, lorsqu’il s’agit d’un système solide en équilibre, on a : (→
−
v = 0 et Ω = 0
→
−
Tandis que pour le cas d’un fluide en équilibre, on a : (→
−
v (M ) = 0 ; ∀M =⇒→
−
vA =
→
− →
− −−→
v G + Ω ∧ GM
D →
−
E D→
− E →
− P→
−
càd. m. ddtv i = f i =⇒M. d dt
vG
= f ext
Soit un point matériel M de fluide (élément de volume dV = dxdydz). En plus des forces
standards (poids) qui s’exerce sur l’élément il y a aussi les forces de pression. Dans un
récipient, il y a une infinité de particule fluide.
P→− →
−
Au repos, le bilan de forces : f ext = 0 , le fluide est à l’équilibre statique et il est
soumis à des forces standards :
— poids : dm→
−
g =ρdV →
−
g = ρdxdydz →
−
g
Figure 2.1 – Bilan de forces sur un élément de fluide dV
— autres forces éventuelles (magnétique, inertie)

— forces de pression : comme les forces exercées sur chaque parois solide.
2.2.1 Forces de pression d’un fluide en repos
2.2.1.1 Cas de contrainte normale
Pour un élément de volume dV = dxdydz en équilibre, les projections des forces exercées
sur cet élément par rapport aux axes donnent :
→
− →
− →
− →
− →
−
OX : d F x = F x − F x+dx = Px .dydz. i − Px+dx .dydz. i avec : Px+dx = Px + dPx
→
− →
− ∂f
d’où : d F x = −dPx .dydz. i =⇒ Fx = − ∂P
∂x
dxdydz ; car :dfn = ∂n
dn
de même pour les autres axes :

→
− →
− →
− →
− →
−
OY : d F y = F y − F y+dy =⇒d F y = −dPy .dxdz. j =⇒ Fy = − ∂P
∂y
dydxdz
→
− →
− →
− →
− →
−
OZ : d F z = F z − F z+dz =⇒d F z = −dPz .dxdy. k =⇒ Fz = − ∂P ∂z
dzdxdy
→
− →
− →
− ∂P →
−
h i
alors, la variation de force de pression est écrite : d F pression = − ∂P
∂x
i + ∂P
∂y
j + ∂z
k dxdydz
−−→ − →
− →
− →
−
on a : dP = gradP.d→
r avec :d→
−
r = dx i + dy j + dz k
→
− −−→
donc : d F pression = −gradP.dV ; c’est l’effet de variation des forces de pression.
→
−
et pour les forces volumique (gravitationnelle), on a : d F gravité = ρ.→
−
g .dV
P→
− →
−
Donc, à l’équilibre : f ext = 0
→
− →
− →
−
=⇒ d F pression + d F gravité = 0
−−→ →
−
=⇒−gradP.dV + ρ.→
−
g .dV = 0
−−→
=⇒ gradP = ρ.→
−
g
Figure 2.2 – forces de vérin hydrostatique
Vérin hydrostatique :
Soit un dispositif de vérin hydraulique, le système est en équilibre et les pistons sont
maintenu dans le même niveau. Si on applique une force F1 de 100N sur le piston P1,
quelle doit être la charge possible à soulever par le vérin. On donne le rapport de diamètres
(S2 /S1 =10).
F1 : La force de pression au point 1 en [N].
F2 : La force de pression au point 2 en [N].
S1 : La surface pressée 1 en [m2 ].
S2 : La surface pressée 1 en [m2 ].
F1 F2
S1
= S2
Dans cet exemple, les forces de pression de variation importante sont celles agissent dans
la direction verticale. Faisant un bilan de forces sur (OZ) pour chacun des deux cotés du
système.
On a à l’équilibre du vérin : F2 + Patm S2 − P2 S2 = 0 ⇒ F2 + P S2 = 0
et à l’équilibre de la pompe : F1 + Patm S1 − P1 S1 = 0 ⇒ F1 + P S1 = 0

F2 F1
à l’équilibre de tout le système, on a : z1 = z2 et P1 = P2 d0 où : S2
= S1
, donc : F2 =
F1 (S2 /S1 ) = 100 × 102
=⇒ F2 = 10KN . (c’est la charge qui peut être soulevée).

Poussée d’Archimède :
On prend l’exemple d’un corps complètement émergé dans un fluide, ce dernier reçoit de
la part de ce fluide une force verticale appelée (poussée d’Archimède), dirigée vers le haut
dont l’intensité est égale au poids du volume du fluide déplacé (ce volume est donc égal
au volume immergé du corps), voir figure ci-dessous (1.2).
Figure 2.3 – Poussée d’Archimède
La différence d pression est : FB − FA = ρf .Vimmergé .g
La force de poussée est donnée comme suit : FBA = ρf .Vimmergé .g

P→
−
à l’équilibre, on a : F ext = 0 =⇒ FB − FA − P = 0 ⇒ Π.A − P = 0
i
avec : P = m.g. = ρ.V.g
Alors : Π.A − ρ.V.g = 0 , d’où : Π.A = ρ.V.g
2.2.1.2 Cas de contrainte tangentielle
On prend l’expérience de Couette, le cas d’un fluide entre deux plaques parallèles de même
surface A, distantes de h (voir Figure 2.3), l’une fixe et l’autre est mobile. La plaque mobile
se déplace avec une vitesse horizontale U et la force qui cause son déplacement est F.
Il est montré expérimentalement que suite à l’effet de la force F, il résulte une contrainte
de cisaillement dans le fluide entre les deux plaques. Cette contrainte (tension) est opposée
au mouvement de la plaque et du fluide, elle est exprimée par : τ = FA mN2 .

N
Remarque : L’unité de cette contrainte (de cisaillement) est du Pascal : m2
= P a, il
s’agit d’une force F par rapport à une surface S, seulement cette force est est de direction
tangentielle et elle n’est pas normale comme pour la pression ordinaire notée par P.
Figure 2.4 – Expérience de viscosité
Il est prouvé aussi que la contrainte de cisaillement est proportionnelle au gradient de

∂u
vitesse u suivant la distance h entre les deux plaques : τ ∼ ∂y
Autrement dit, la contrainte de cisaillement est proportionnelle au taux de déformation

et la constante de proportionnalité est le coefficient de viscosité dynamique (µ) :
∂u
τ =µ (2.1)
∂y
Fs N.s kg
Contrainte de cisaillement
Un réarrangement de l’éq. (1) donne : µ = taux de cisaillement
= As
dU = τ dU
dy
; m2
= m.s ou
dy
[P a.s]
La viscosité cinématique (ν) est définie comme étant : la viscosité dynamique sur la
h 2i
masse volumique, son unité est [m²/s]. Elle est exprimée comme suit : ν = µρ ms
2.2.1.3 Variation de force de pression
Pour les gaz (fluide compressible) et dans le cas de l’air (atmosphère), la pression varie
avec l’altitude, et obéie à la loi des Gaz Parfait. (figure 2.4)
Mg
P (z) = P0 exp−( RT )z

Mg
− z0
tel que : P0 = exp RT0
, pour l’atmosphère à : T0 = 300K
Tendis que pour les liquides (fluide incompressible), la pression varie avec profondeur et
obéie à la loi de l’hydrostatique. (figure 2.4)
Figure 2.5 – variation de pression atmosphérique et sous marine [6]
Figure 2.6 – Liquides en mouvements uniformes
P (z) = P0 + ρgh
avec : h = z − z0 est la différence de profondeur et P0 = P (z0 ).
2.2.2 Forces de pression d’un fluide en mouvement (uniforme)
L’étude de déformation de la surface libre d’un liquide en mouvement uniforme (à vitesse
constante). Par exemple, les deux cas (connus), celui du réservoir (contenant un liquide)
→
−
en translation à vitesse linéaire constante ( V = Cste) et le second cas est le vortex forcé
qui tourne à vitesse angulaire constante (ω = C ste ). Dans cette catégorie de problème, il
est nécessaire d’appliquer la loi fondamentale de l’hydrostatique. (figure 2.5)
Figure 2.7 – Conduite avec tube en U
2.2.3 Forces de pression d’un fluide en mouvement (quelconque)
L’étude de mouvement quelconque d’un fluide (ex. conduites), nécessite parfois la mesure
de pression locale ou totale à des endroits bien spécifiés. Pour cela dans l’ingénierie, les
instruments de mesure souvent employés utilisent des techniques basées sur l’application
de la loi fondamentale de l’hydrostatique, par exemples : les baromètres en U, à cuvette,
à siphon, les tubes de Pitot et les l’effet de Venturi...etc.
Bien sûr, ça nous ramène à introduire la notion de conservation d’énergie dans un volume
de contrôle (V.C.). C’est à dire, on va ajouter un 3ème terme à l’équation fondamentale
de l’hydrostatique et qui deviendra l’équation (très connue) de Bernoulli.
Dans ce genre de problème, on fait appel à la fois à l’équation de Bernoulli et aussi

à l’équation de l’hydrostatique (au niveau des sondes de mesure de pression). quelques
exemples d’application sont présentés ci-après.
Application 1 : conduite avec tube en U
De l’huile de densité 0,75 circule à travers une conduite horizontale (figure 2.7), et fait
monter le mercure dans le manomètre en "U". Calculer la valeur de la hauteur h, si la
pression mesuré au point A est égale à 1, 4kg/cm2 . On donne : z= 0,825 m.
Réponse :
On a : PB = PC (points B et C même niveau et même fluide).
En appliquant la loi fondamentale de l’hydrostatique entre les points A-B et C-D, on

Figure 2.8 – Principe de la sonde de Pitot
trouve :
=⇒ PA + ρ.g.(z + h) = PD + ρm .g.h
=⇒ h.g.(ρm + ρ) = PA − PD + ρ.g.z
PA −PD +ρ.g.z
=⇒ h = g.(ρm +ρ)
1.4×105 −101325+750×10×0.825
A.N. :=⇒ h = 10.(113600+750)
=⇒ h = 1.14m
Application 2 : sonde de Pitot tube en U
La sonde de Pitot permet de mesurer la vitesse d’un écoulement d’air autour d’une struc-
ture, par exemple un avion.
A l’entrée du tube, en A, la vitesse de l’air est nulle. Au point B, la vitesse de l’air est
sensiblement égale à sa valeur non perturbée par l’instrument. L’équation de Bernoulli
nous donne :
=⇒ PA = PB + ρ2 v 2
pρ
donc, la vitesse de l’air est : v = 2
(PA − PB ) m/s
L’équation de Bernoulli permet de comprendre le principe de divers instruments et ap-

pareils tels que trompe à eau, bec Bunsen, pistolet à peinture, carburateur, etc. . . En
revanche, elle ne permet pas d’expliquer les pertes de charge dans les tuyaux, car le fluide
est supposé non visqueux.
Application 3 : deux conduites avec tube en U
Soit deux conduites de sections SA et SB , qui contiennent de l’eau aux pressions respectives
de 2,80 et 1,40 bar. on demande de calculer la dénivellation h de mercure du manomètre
différentielle. On donne : x + y = 2m. La densité du mercure est : d = 13, 57.
Figure 2.9 – Différence de pression entre deux conduites
En appliquant l’équation de l’hydrostatique entre les points A et C, on trouve :
PC − PA = ρeau .g.(ZA − ZC )=⇒ PC = PA + ρeau .g.(h + x)
En appliquant l’équation de l’hydrostatique entre les points C et D, on trouve :
PD − PC = ρmercure .g.(ZC − ZD ) =⇒ PD = PC + d.ρeau .g.(−h)
=⇒ PD = PA + ρeau .g.(h + x) + d.ρeau .g.(−h)
=⇒PD = PA + ρeau .g.x + ρeau .g.h(1 − d)
En appliquant l’équation de l’hydrostatique entre les points D et B, on trouve :
PB − PD = ρeau .g.(ZD − ZB ) =⇒ PB = PD − ρeau .g.y
=⇒ PB = PA + ρeau .g.(x + y) + ρeau .g.h(1 − d)
PB −(PA +ρeau .g.(x+y)) 2,8×105 −(1,4×105 +1000×9,81×2)

=⇒ h = ρeau .g.(1−d)
= 1000×9,81×(1−13,57)
AN : h = 0, 95m.
2.2.4 Application de de l’équation de BERNOULLI (approche
statique)
Le passage de la loi fondamentale de l’hydrostatique à l’équation de Bernoulli, s’effectue

selon le principe de conservation d’énergie. Tout en considérant la variation de l’énergie
cinétique, qui est (toujours) égale à la somme des travaux des forces extérieures.
δW = −dEP + δwf orces pression =⇒ 12 v 2 + ρg z+ Pρ =C te

P
dEc =
et l’énergie totale d’un fluide statique est conservé tout point de ce fluide.
càd. Etotale = Ecinétique + Epotentiele + Elocale = C te
2.2.4.1 Effet de gravitation sur un liquide
Application de vidange d’un réservoir (formule de Torricelli) Considérons un

réservoir de grande section, rempli d’un liquide qui s’écoule à travers un orifice de section
faible (devant celle du réservoir) et situé à une hauteur sous la surface libre (voir figure).
La pression atmosphérique P0 s’exerce à la fois sur la surface libre et sur le jet à la sortie
de l’orifice. En supposant que le liquide est incompressible et non visqueux, il est possible
d’utiliser l’équation de Bernoulli pour déterminer la vitesse de vidange du réservoir.
Figure 2.10 – Vidange de réservoir
Considérons un point A de la surface libre et un point B situé dans le jet. Il existe une
ligne de courant passant par ces deux points et le long de laquelle on peut appliquer
l’équation de Bernoulli :
PA + ρgZA + 12 ρUA2 = PB + ρgZB + 21 ρUB2
Compte tenu du rapport de section entre le réservoir et l’orifice, on peut y négliger

UA devant UB (UA << UB ), car la surface du récipient est grande vis-à-vis de celle de
la section l’orifice : UA ≈ 0 . Par ailleurs, le point A étant situé sur la surface libre,
on peut admettre que : PA =PB pour de petites différences d’altitude. Dans le domaine
du jet où se trouve le point B, les lignes de courant sont rectilignes et parallèle et par
conséquent la pression motrice P est uniforme sur une même section. En négligeant les
variations d’altitude au sein du jet, il s’ensuit que la pression statique est également à
peu près uniforme sur une même section. Or, par continuité de la pression à l’interface
jet-atmosphère, la pression atmosphérique règne donc en tout point du jet (ce résultat
peut être généralisé à toute situation dans laquelle un liquide s’écoule sous la forme d’un
jet libre).
Appliquant Bernoulli entre A et B : PA + ρgZA + 12 ρUA2 = PB + ρgZB + 12 ρUB2
Comme : PA =PB et UA ≈ 0, l’équation de Bernoulli prend la forme suivante :
⇒ P0 + ρgZA + 0 = P0 + ρgZB + 12 ρUB2
=⇒ 12 ρUB2 = ρg (ZA − ZB ) = ρgH
qui, après simplification, permet d’écrire la formule non-uniformede Torricelli : UB =

√
2gH
Remarque :
La vitesse d’écoulement du liquide et le temps de vidange sont indépendants de sa masse

volumique du liquide ; le seul paramètre dont ils dépendent est donc la hauteur séparant
l’orifice de la surface libre du liquide.
Exercice d’application (Cas d’un Siphon) : On considère un siphon de diamètre

d=10 mm alimenté par un réservoir d’essence de grandes dimensions par rapport d et
ouvert l’atmosphère. On suppose que :
- le fluide est parfait.
- le niveau du fluide dans le réservoir varie lentement.
- L’accélération de la pesanteur g = 9.81m.s2 .
- le poids volumique de l’essence : ω = 6896N/m3 .
- H = ZA –ZS = 2, 5m.
Figure 2.11 – Siphon
1) En appliquant le Théorème de Bernoulli entre les points A et S, calculer la vitesse

d’écoulement VS dans le siphon.
2) En déduire le débit volumique qV .
3) Donner l’expression de la pression PB au point B en fonction de h, H, et Patm . Calculer

PB pour h = 0.4m.
4) h peut-elle prendre n’importe quelle valeur ? Justifier votre réponse.
Solution :
PS U2 PA U2
1) Appliquant l’équation de Bernoulli entre les points S et A : ω
+ZS + 2ωS = ω
+ZA + 2ωA
on a : PS = PA = Patm , UA = 0 et ZA –ZS = H
√ √
US = 2gH, A.N. : US = 2 × 9.81 × 2.5 = 7m/s
2 2
2) Le débit volumique : qv/S = US ×S = US × πd4 A.N. : qv/S = 7× 3.14×(0.01)
4
=5.5×104 m3 /s =
0.55l/s
2
UB US2
PB PS
3) Théorème de Bernoulli entre B et S : ω
+ ZB + 2ω
= ω
+ ZS + 2ω
Or : US = UB , ZB − ZS = (H + h) et P s = P atm
PB = Patm – ω(H + h) A.N. : PB = 105 – 6896 × (2, 5 + 0, 4) = 80001, 6P a = 0, 8bar
4) Non ! Il faut que PB > 0

Patm 105
Équivaut à : h < ω
− H, A.N. : h < 9.81×700
− 2.5 = 12m.
Exercice d’application (tube de Pitot) : Le tube de Pitot, représenté ci-dessus,

permet de mesurer la vitesse (ou le débit) d’un fluide dans une canalisation de section S.
Le fluide circulant dans la canalisation est animé d’une vitesse U et possède une masse
volumique ρ1 alors que le tube en U est rempli d’un liquide de masse volumique ρ2 .
Enoncé
Figure 2.12 – Tube de Pitot
Le fluide est incompressible (ρ = Cste)et parfait (µ = 0) . - L’écoulement est permanent

(∂U/∂x = 0).
Nous pouvons donc appliquer la relation de Bernoulli sous la forme :
P + ρgZ + 21 ρU 2 = Cste
— L’ouverture B est dirigée face au jet et constitue un point d’arrêt (la vitesse du
fluide en B est nulle, UB = 0).
— ZA = ZB
— UB = UB 0 = UA0 = 0
— La vitesse dans le tube est nulle (pas de mouvement), ce qui nous laisse écrire en
appliquant le principe fondamental de la statique des fluide ans cette partie du
système :

PB 0 − PA0 = ρ2 gh





PA0 − PA = ρ1 g (ZA − ZA0 )



 PB 0 − PB = ρ2 g (ZB − ZB 0 )

PA + ρ1 gZA + 12 ρ1 UA2 = PB + ρ1 gZB + 12 ρ1 UB2
=⇒ PB = PA + 12 ρ1 UA2
PA0 − PA = ρ1 g (ZA − ZA0 )
PB 0 − ρ2 g (ZB − ZB 0 ) = PA0 + 21 ρ1 UA2 − ρ1 g (ZA − ZA0 )

r
ρ2
=⇒ UA = 2hg ρ1
−1
2.3 Théorèmes de transport
2.3.1 Volumes et surfaces de contrôle
Un volume de contrôle est un volume imaginaire sur lequel au procède au bilan intégral
d’une grandeur physique comme la masse, la quantité de mouvement ou encore l’énergie.
On appelle surface de contrôle l’enveloppe délimitant un volume de contrôle.
Le fluide peut entrer et sortir d’un volume de contrôle c’est-à-dire traverser la surface de
contrôle qui peut elle même être fixe ou mobile.
Le volume de contrôle peut être sous différente formes. Il peut être fixe (constant),
matériel à volume variable (masse constante), ou encore arbitraire dont le déplacement et
la déformation sont différents de ceux du fluide qui le traverse.
2.3.2 Théorème de transport
En Mécanique des Fluides, l’évolution des grandeurs matérielles est analysée à l’aide
d’équations intégrales de bilan sur des domaines fluides macroscopiques. Le transport
de ces grandeurs dans l’écoulement est explicité en suivant le mouvement ; il est par
conséquent nécessaire d’établir l’expression de la dérivée particulaire d’une intégrale vo-
lumique.
Généralement, un volume de contrôle arbitraire est noté par V a(t) limité par l’enveloppe
fermée Sa(t). Soient f (→
−
x , t) une fonction scalaire continue et dérivable et I(t) son intégrale
sur le volume V a.
I(t) =
f (→
−
x , t)dV
Va
Le taux de variation dans le temps de l’intégrale I(t) est donné par :
df →
−
d
dt
I(t) = d
dt
f dV = dt
dV + f(V a • →
−
n )dS —————-(*)
Va Va (t) Sa (t)
→
−
où V a désigne la vitesse locale de la surface de contrôle et →
−
n la normale extérieure. La
relation précédente (*), est connue sous le nom de la règle de Leibnitz, s’interprète comme
suit :
Taux de variation de Intégrale de la variation Flux de f (→

−
x , t)
l’intégrale de f (→
−
x , t) = temporelle de f (→
−
x , t) + à travers
sur le volume V a(t) sur le volume V a(t) l’enveloppe Sa(t)
Chapitre 3
Équations d’Euler (fluides parfaits)
3.1 Rappel de cinématique
La cinématique des fluides est l’étude des mouvements de fluides sans tenir compte des
forces exercées sur le fluide. Dans la cinématique, on s’intéresse uniquement à l’étude du
champ des vitesses. on trouve deux approches différentes, l’une est de Lagrange (approche
Lagrangienne) et l’autre est de Euler (approche Eulérienne).
3.1.1 Approche de Lagrange
Le point de vue Lagrangien, est de suivre l’élément du fluide dans son mouvement, càd.
suivre sa position à tout moment ; i.e. connaitre les coordonnés du point M(x,y,z) à chaque
instant t .
→
− d→
−r →
−
r (t + dt) − →
−
r (t)
v = = lim (3.1)
dt dt→0 dt
une
→
− d→
−v d2 →
−
r
a = = (3.2)
dt dt2
−−→ −
→ −
→ −
→ −
→ − →
OM (t) =→
−
r (t,→
−
r 0 )=⇒→
−
v (t,→
−
r 0 )= dtd (→
−
r (t, →
−
r 0 )) = lim r (t+dt, r 0 )−dtr (t)− r (t, r 0 )
dt→0
33
CHAPITRE 3. ÉQUATIONS D’EULER (FLUIDES PARFAITS) 34
Figure 3.1 – Approche de Lagrange

Les lignes de trajectoires représentent les chemins suivis par les particules fluides dans
l’espace avec l’évolution du temps, comme l’explique la figure 3.2.Théorèmes de
transport
Figure 3.2 – Lignes de trajectoires
3.1.2 Approche de Euler
Du point de vue Eulérien, on s’intéresse au champ des vitesses de la particule fluide

(élément de volume) en mouvement →
−
v (→
−
r , t). L’équation de Euler détermine l’évolution
des vitesses dans un volume de fluide. Donc, la question est comment déterminer cette
équation ?
Figure 3.3 – Approche Eulérienne
Les lignes de courant représentent des profiles de l’espace (courbes), décrivent le fluide
Figure 3.4 – lignes de courant
en mouvement et qui, à tout instant, possède en tout point une tangente parallèle à la
vitesse des particules de ce fluide.
3.2 Dynamique des fluides parfaits
3.2.1 Notion de dérivée particulaire :

−
→
v (−
→
r +d−
→
r ,t+dt)−−
→
v (−
→
on a l’accélération d’un élément de fluide : →
− d−
→v r ,t)
a = dt
= lim dt
dt→0
Si on a une fonction f(x,y,z,t), on sait que : f (x+dx, y+dy, z+dz, t+dt) = f (x, y, z, t)+ ∂f
∂x
dx+
∂f ∂f ∂f
∂y
dy + ∂z
dz + ∂t
dt
→
− −
→ −
→ − → −→
v (−
→
1 ∂−
h →
∂−→ ∂−→ ∂−
→
h −
∂→ ∂−
→ ∂−
→
i
a = lim v ( r +d r ,t+dt)−
dt
r ,t)
=lim dt
v
∂x
dx + v
∂y
dy + v
∂z
dz + v
∂t
dt = lim v dx
+ v dy
+ v dz
+ ∂
dt→0 dt→0 dt→0 ∂x dt ∂y dt ∂z dt
=⇒→
−a = lim ∂∂xv vx + ∂∂yv vy + ∂∂zv vz + ∂∂tv =⇒→
h −
→ −
→ −
→ −
→ − →
− −
→
i
∂ ∂ ∂
a = vx ∂x + vy ∂y + vz ∂z v + ∂∂tv
dt→0
→ − −
=⇒→
−
a = →
−
v .∇ → ∂−
→v
v + ∂t
→ − →
d−
→ ∂−
→
— Dérivée particulaire : v
dt
= v
∂t
+ v .∇ −
→
− v
C’est la variation temporaire d’une particule fluide, du point de vue Lagrangien, en la

suivant dans son mouvement,
avec : d et →
−
r est le vecteur position.
d→
−v ∂→
−v → − −
d0 où : = + →
−
v .∇ →v (3.3)
dt ∂t
Lagrangien Eulerien
Donc, l’approche de dynamique des fluides parfaits est basée sur l’étude des écoulements
de fluides, tout en considérant les forces extérieurs sans prise en compte des forces de frot-
tement ( ou forces visqueuses). Les écoulements parfaits sont gouvernés par les équations
connues sous le nom des Équations de Euler .
3.2.2 Équation d’Euler :
d−
→v →
−
Pour une particule fluide de volume dV et de masse dm=ρdV =⇒ ρdV dt
=Σ f ;
→
−
f : forces extérieures qui influent sur la particule fluide.
→
∂−

v →
→
− − →
− →
−
=⇒ ρdV + v .∇ v = Σ f (3.4)
∂t
Il y a les forces volumiques (poids) et les forces de pression.

→
− →
− −−→
poids : d f poids =ρdV →
−
g et pression :d f pression = −dV gradP
on peut traiter le cas particulier des forces de poids et de pression, donc on obtient :
→ − − →
− → − →
−
ρdV ∂∂tv + →
−
→ −v .∇ → v =Σ f =d f poids +d f pression
→ − − −−→ → − − → −−→
=⇒ρdV ∂∂tv + ρdV →
−
→ −
v .∇ → v =ρdV →−g −dV gradP =⇒ρ ∂∂tv + ρ →
−
→ −
v .∇ → v =ρ−g −gradP
∂→
−v → − − 1 −−→
=⇒ + →
−
v .∇ →v = →
−
g − gradP (3.5)
∂t terme de gravité ρ terme pression
terme temporaire terme convectif
(6)99KC’est la dite équation de Euler (1755).

→ −
Soit une fonction G ; sa dérivée particulaire est : dG
dt
= ∂G
∂t
+ →
−
v .∇ G
Lagrangien Eulerien
→ −
=⇒ dG
= ∂G
+ →
−
v .∇ G
dt ∂t
terme temporaire terme local transport ou terme convectif
Pour le transport d’une grandeur physique dans l’espace, il y a deux façons présentées
par :
→ le terme d’advection (de convection ou convectif).

→ le terme de diffusion (de conduction).
Remarque : Dans la dynamique des fluides, le terme de diffusion est lié aux forces de
frottement entre les particules et qui dépend de la viscosité du fluide, appelé aussi le terme
visqueux. Ce terme n’apparaitra pas ici avec Euler), il s’agit du cas d’écoulements réels,
qui seront présentés avec les équations de Navier-Stokes dans le dans la chapitre suivant.
3.2.3 Équation de BERNOULLI : ( autres approches)
Approche Bernoulli
Euler :
En partant du principe de thermodynamique, par l’application de la loi de conservation

d’énergie sur un volume de contrôle :
→
− →
− →
−
dEc = ΣδW f ext =δW f poids + δW f pression =ρgdz − dP
h →
d−v →
−i → −
et d’autre part, on a : m dt = Σ f .d r
avec : d→
−
r le vecteur de déplacement
→
− →
− →
− →
−
=⇒m ddtv .d r = Σ f .d r
→
− −
=⇒m ddtr .d→
v = ΣδW
=⇒m→
−
v .d→
−
v = ΣδW
1 →
m−

sur l’axe (OZ) : =⇒d 2
v2 = ΣδW
→
−
et on a : ΣδW = δWpoids + δWf.pression = m→ −g d→
−
r − V ∇P.d→
−
r
d’où, on obtient sur l’axe (OZ) : d 21 ρV v 2 = −ρV gdz − V dP

dz
.dz
1

=⇒ d 2
ρv 2 =−ρgdz − dP
après intégration, on aura : 12 ρv 2 =−ρgz − P + C
⇒ 12 ρv 2 + ρgz + P = C 99Kl’équation de Bernoulli sur l’axe OZ.

on dérive l’équation de Bernoulli, par rapport à (OZ) : =⇒ d( 21 v 2 ) + ρ1 dP − gdz = 0
→
−
h i
1 ∂ 2 1 ∂P
⇒ 2 ∂z
(v ) + ρ ∂z
− g .dz = 0
−→
−
g .d→
−
h i
1 ∂ 1 ∂P

=⇒ 2 ∂−
→
r
(v 2 ) + ρ ∂−
→
r
r =0
h −−→
1 −−→ →
− .d→
−
i
1 2
=⇒ 2
grad (v ) + ρ
gradP − g r =0
−−→ −−→
=⇒ 12 grad (v 2 ) + ρ1 gradP − →
−
g =0
−−→ → − − −→→
d’après le théorème de Green Ostrogradski : 12 gradv 2 = →
−
v .∇ → v − rot−
v ∧→
−
v
→ − − −→→ −−→
alors ; on peut écrire : →
−v .∇ →v − rot−v ∧→ −v + ρ1 gradP − →
−

g =0
→ − → −→− → −−→
→
−
⇒ v .∇ −v − rot→
v ∧−
v =→
−
g − ρ1 gradP
−→− →
−
pour un fluide irrôtationnel rot→
−
→
v = 0 et instationnaire ( ∂∂tv 6= 0) :
∂−
→
→ − → −−→
=⇒ v
∂t
+ v .∇ −
→
− v =→
−
g − ρ1 gradP 99Kc’est l’équation d’Euler
Approche Euler
Bernoulli :
−
→
Maintenant, on prend la démarche inverse et partant de l’équation d’Euler , on a : ∂∂tv +
→ − − −−→
→
−v .∇ → v =→ −g − gradP
ρ
→ − − −−→ −→−
et en utilisant le théorème de GREEN OSTROGRADSKI : → −v .∇ →v = 12 gradv 2 + rot→v∧
→
−
v
on aura :
∂−
→ −−→ −→− → −−→
⇒ v
∂t
+ 12 gradv 2 + rot→
v ∧−
v =→
−
g − gradP
ρ
−→− →
−
pour un écoulement incompressible stationnaire ( ∂∂tv = 0) et irrôtationnel (rot→
−
→
v = 0 ),
l’équation précédente devient :
−−→ −−→
⇒ 12 gradv 2 = →
−
g − gradP
ρ
et le long d’une ligne de courant, on a :d→

−
r=→
−
v dt ; tel que (→
−
r ) est le vecteur déplacement
d’une particule quelconque d’un fluide incompressible. d’où, on peut écrire :
h −−→ −−→ P →
− .d→
−
i
1 2
=⇒ 2
gradv + grad( ρ
) − g r =0
−→
−
g .d→
−
h i
1 ∂ ∂ P
=⇒ 2 ∂−
→
r
v2 + ∂−
→ ( )
r ρ
r =0
→
−
h i
1 ∂ 2 1 ∂P
après la projection sur l’axe (OZ) : 2 ∂z
v + ρ ∂z
− g .dz = 0
=⇒ d( 12 v 2 ) + ρ1 dP − gdz = 0 ; et en intégrant :
=⇒ 12 ρv 2 + P − gρz = C ste 99K équation de Bernoulli sur l’axe OZ

0 dx

avec : →
− r = dy =⇒projection sur OZ : →
g = 0 ; d→
− −
g .d→
−

r V g.dz

g dz

3.2.4 Écoulement irrôtationnel :

−→− →
−
Si : rot→
v = 0 partout , alors :
1
2
ρv 2 + P − gρz = C ste (écoulement uni-directionnel)
si : un écoulement est irrôtationnel initialement (µ = 0) viscosité négligeable, il restera

irrôtationnel dans tout le fluide.
Exercice 1 :
On considère un siphon de diamètre d = 10 mm alimenté par un réservoir d’essence de

grandes dimensions par rapport à d et ouvert à l’atmosphère. On suppose que :
- le fluide est parfait.
- le niveau du fluide dans le réservoir varie lentement.
- l’accélération de la pesanteur g = 9.81m/s2 .
- le poids volumique de l’essence : ω = 6896N/m3 .
- H = ZA − ZS = 2, 5m.
1) En appliquant le Théorème de Bernoulli entre les points A et S, calculer la vitesse

d’écoulement vs dans le siphon. 2) En déduire le débit volumique qv .
3) Donner l’expression de la pression PB au point B en fonction de h, H, ω et Patm . Faire

une application numérique pour h = 0.4m.
4) h peut elle prendre n’importe quelle valeur ? Justifier votre réponse.
Réponse :
Chapitre 4
Équation de continuité (conservation

de matière)
La matière a une structure discontinue et la notion de milieu continu est une pure no-
tion schématique . Elle consiste à admettre que la masse et toutes ses propriétés sont
réparties continûment dans le matériau (ce qui n’exclut pas les discontinuités aux in-
terfaces). Bien entendu ce schéma ne prétend représenter que les phénomènes macro-
scopiques ou mésoscopique dont les échelles caractéristiques sont très grandes devant la
distance intermoléculaire moyenne. Comme il n’est pas question d’ignorer complètement
les phénomènes dont le siège est à l’échelle moléculaire (comme celui de la diffusion), ceux-
ci devront être représentés à travers une description macroscopique (ou mésoscopique) de
leurs conséquences à grande échelle. Le concept du continuum présente l’immense avan-
tage d’autoriser le calcul différentiel et intégral dont les outils sont présentés au premier
chapitre d’analyse vectoriel.
4.1 Courant de matière
Soit un flux de matière qui travers une surface élémentaire dS. La quantité de matière
−→ −
→
traversant dS pendant un temps dt est : dm = ρqV = ρV.dt.dS = ρ→ −
v. dS.dt= ρ→
−v. →
−
n ⊥ .dS.dt
−
→
⇒ dm =ρ→
−
v. dS.dt
→
−
Par définition : j = ρ→
−
v est la densité de courant de matière (ou de masse).
41
CHAPITRE 4. ÉQUATION DE CONTINUITÉ (CONSERVATION DE MATIÈRE) 42

dm − −
→ →
donc, dt /dS
= j .dS
Figure 4.1 – Flux de matière
La quantité de matière qui travers dS (ou dite aussi l’intensité de masse), par analogie au
=→ − −→ − −
→ →
courant électrique : dq

dt /dS
j .dS99Kϕ/dS = j .dS (flux).
− −
→ →
−
→
I/S = j .dS = →
−
v. dS (4.1)
S S
Figure 4.2 – Courant de matière

= →

− −→
exemple d’une cExercice 1 : analisation, le flux de matière est conservé : I/S1 =qm1 = dm
dt /S1
j .dS
S
qm1 = qm2 [kg/s] : le débit massique est conservé le long d’un tube de lignes de courant
(figure 4.3), entre l’entrée (section S1) et la sortie (section S2), c’est similaire à une
conduite réelle.
Par analogie au circuit électrique, le courant est conservé le long du circuit (figure 4.4).
Figure 4.4 – Analogie au circuit électrique

Enoncé
Figure 4.3 – Conservation de débit d’une conduite
Un écoulement stationnaire et incompressible, le flux de matière est conservé, comme dans

le cas du courant le long d’un dans un fil électrique, et donné par : I/S1 =I/S2 =I/Si .
4.2 Bilan de matière (conservation de masse)
Soit un élément de volume dxdydz : dm = ρdV = ρdxdydz
Figure 4.5 – Bilan de matière suivant ox
→
−
Pendant dt la quantité de matière entrante dans l’élément suivant (ox) est : j (x) .dS.→
−
n x .dt =ρ.→
−
u .dydz.→
−
n
→
−
La quantité qui sorte suivant (ox) est : − j (x+dx) .dS.→
−
n x .dt =−ρ.→
−
u .dydz.→
−
n x .dt
Le bilan de matière sur suivant trois axes :

→
− →
−
(ox) : j (x) .dydz.→
−
n x .dt = − j (x+dx) .dydz.→
−
n x .dt = 0
→
− →
−
(oy) : j (y) .dxdz.→
−
n y .dt = − j (y+dy) .dxdz.→
−
n y .dt = 0
→
− →
−
(oz) : j (z) .dydx.→
−
n z .dt = − j (z+dz) .dydx.→
−
n z .dt = 0
→
− →
− − → − − →
− −
on a : j = j x .→
n x + j y .→
n y + j z .→
nz
−
→ −
→ −
→ −
→ −
→ −
→

j (x+dx) − j (x) →
− j (y+dy) − j (y) →
− j (z+dz) − j (z) →
−
alors : dm = −dxdydz.dt dx
nx+ dy
ny + dz
nz
−
→ −
→ −
→
→
j (x) →
− j (y) →
− j (z) →
− dm − →
− d(ρdV )
→
− →
−
⇒dm = −dV.dt dx
nx + dy
ny + dz
nz ⇒ dt
= −dV ∇. j ⇒ dt
= −dV ∇. j
→
− →
−
⇒ d(ρ)
dt
= −div j = −div ρ U
Donc :
dρ →
−
+ div ρ U = 0 (4.2)
dt
C’est la dite équation de continuité ou bien dite : équation de conservation de masse .
4.3 Flux de matière d’un volume fini "V"
Soit un volume fini de fluide "V" et de surface d’enveloppe Σ(V).
Donc : dM
=
∂ρ dV ; avec : M = ρdV
dt ∂t
V V
Par ailleurs, le bilan de matière donne :
! →
−
Flux de matière entrant à travers l’enveloppe Σ : dM
dt
= j(→
−
r , t).d S in
!
Σ(V )
→
− → →
−
Flux de matière sortant à travers l’enveloppe Σ : dM = j (−
! →−j .d→−S = div→−j .dV
dt
r , t).d S out

Σ(V )
D’après le théorème de Green Ostrogradski :

Σ(V ) V
→
−
dρ
Si le fluide est incompressible (ρ=Cste ) ⇒ dt
+ div j = 0
→
−
⇒ div ρ U = 0
→
−
⇒ div U = 0
Courant à travers la surface de l’enveloppe fermée (V) est :
! − →
→ − →
−
I/Σ = j .d S = div j .dV = 0
Σ(V ) V
Application :
La loi des nœuds ( Σ de tout ce qui entre = Σ de tout ce qui sort) à travers la surface de
l’enveloppe :
→
− ! − →
→ − ! − →
→ − − →
→ − →
− →
−
div j = 0 ⇒ j .d S out = 0=⇒ j .d S out = − j .d S f low + j .d S f low +
→
− →
−
Σ(V ) Σ(V ) S1 S2
j .d S f low
S3
⇒ 0 = I1 + I2 + I3 ⇒ I1 = I2 + I3 (c’est la loi de nœuds en électricité).
Par analogie à la loi des nœuds d’électricité, la masse est aussi conservé dans une bifur-
cation de conduite, qui est considéré comme un nœud.
Selon la loi des nœuds la masse est conservée et le débit aussi, donc la somme des débit
entrant est égale à la somme des débits sortant dans le cas de la figure 3.6.
Figure 4.6 – Application de la loi des nœuds
→
−
Exemples de divergence "div( U ) " d’écoulement de fluide : la divergence représente la
variation de débit entre les lignes de courant pour un écoulement d’un fluide incompressible
(figure ci-après).
écoul. uniforme écoul. laminaire source puits
→
− →
− →
− →
−
div( U ) = 0 div( U ) = 0 div( U ) > 0 div( U ) < 0
Figure 4.7 – Exemples de divergence
4.4 Équation de continuité
Rappel : L’équation de continuité traduit le principe de conservation de masse.
Soit un élément de volume dV = dxdydz, sa masse est dm = ρdV = ρdxdydz.
Si la matière est conservée, sa variation dans le temps est nulle, pendant un temps dt on
dm
a: dt
=0
→
−→−
alors : ∇ U = dm
dt
= ∂m
∂t
dt + ∂m
∂x
dx + ∂m
∂y
dy + ∂m
∂z
dz
on a : dm
dt
= 0 ⇒0 = ∂m
∂t
dt + ∂m
∂x
dx + ∂m
∂y
dy + ∂m
∂z
dz
∂m
⇒0 = ∂t
+ ∂m∂x
u + ∂m
∂y
v + ∂m
∂z
w

⇒ ∂m
∂t
= − ∂m
∂x
u + ∂m
∂y
v + ∂m
∂z
w ; avec : m = ρdV

⇒ ∂ρdV
∂t
= − ∂ρdV
∂x
u + ∂ρdV
∂y
v + ∂ρdV
∂z
w

⇒ ∂ρ
∂t
dV = − ∂ρu ∂x
+ ∂ρv
∂y
+ ∂ρw
∂z
dV

∂ρ ∂ρu ∂ρv ∂ρw
=− + + (4.3)
∂t ∂x ∂y ∂z
∂ρ →
− →− P
⇒ ∂t
= − ∇ ρ U + ρqv
avec : qv > 0 cas de source et qv < 0 cas puits.

∂− →
→ − P
— écoulement permanent (stationnaire) : = 0⇒ ∇ ρ U = ρqv
∂t
→
− →− →
−→−
— écoulement incompressible : ρ = C ste −→ ∇ ρ U = ρ qv =⇒ ∇ U = qv .
P P
P →
−→−
— écoulement conservatif et incompressible : qv = 0 =⇒ ∇ U = 0 (ni puits, ni
source)
Exercice 1 : Soit un écoulement dont le potentiel des vitesses est donné par :
−−→ →
−
ϕ(x, y) = x2 − 2y − y 2 , avec : →
−
v = gradϕ = ∇ϕ
— Démontrer que l’écoulement est : bidimensionnel, permanent et incompressible
(vérifiant l’équation de continuité) ?
Exercice 2 : L’équation de la continuité pour l’écoulement permanent isovolume est-elle
vérifier ? si les composantes de la vitesse sont les suivantes pour les deux cas (a) et (b) :

u = 2x2 − xy + z 2





a) v = x2 − 4xy + y 2



 w = −2xy − yz + y 2


u = (2x − 3y)t





b) v = (x − 2y)t



0


Chapitre 5
Équations de Navier-Stokes (fluides

réels)
5.1 Généralités et définitions
L’étude des écoulements réels (visqueux) est gouvernée par l’équation dite équation de
Navier-Stokes :
→
−
dU →
− →
−
ρ = − ∇P + µ4 U + ρ→
−
g 99K équationdeN avier − Stocks (5.1)
dt
En pratique, la solution analytique de cette équation est impossible. Généralement, on

pose des hypothèses simplificatrices et on procède aux solutions numériques basées sur des
méthodes numériques (différences finies, volumes finis, éléments finis...etc). Dans toutes les
situations, les forces de frottements ont un grand rôle, ou la viscosité n’est plus négligeable
(fluide réel), ce qui fait la différence avec un fluide parfait.
5.2 Régime d’écoulement
Dans la pratique, parmi les effets de viscosité, on retrouve le changement de la nature et du

comportement d’écoulement avec l’accélération (vitesse) sous l’effet des forces visqueuses.
47
CHAPITRE 5. ÉQUATIONS DE NAVIER-STOKES (FLUIDES RÉELS) 48
Ce qui caractérise ce qu’on l’appel régime d’écoulement. O. Reynolds a réalisé une étude
systématique du régime d’écoulement en fonction des différents paramètres intervenant
(débit, viscosité, dimension de géométrie ... etc). Il a montré que le régime dépend d’une
U ρ`
grandeur sans dimension Re = µ
, appelé Nombre de Reynolds , qui regroupe les
grandeurs dimensionnelles suivantes : la vitesse d’écoulement (U), la masse volumique
(ρ), la viscosité dynamique (µ) et la longueur caractéristique de la géométrie (`). D’un
point de vue sens physique, ce nombre adimensionnel, caractérise le rapport des forces
d’inertie sur les forces visqueuses.
5.3 Tenseur de contraintes
On a un point M qui appartient à une surface dS, la force de surfacique (de pression),
→
−
s’exprime : d F = −→
−
n P dS, qui est normale à la surface dS.
→
− →
−
Pour une contrainte quelconque, on a : d F = − T n dS.
→
− →
− →
−
Soit la norme de cette surface : →
−
n ⊥ = nx i + ny j + nz k
→
−
et dans ce cas : T n = Tx →
−
n x + Ty →
−
n y + Tz →
−
nz
Les projections de cette contrainte surfacique sur les surfaces ⊥aux axes, soient :
→
− →
− →
− →
−
— la projection de dS sur le plan ⊥à X est : T x = σxx i + τyx j + τzx k
→
− →
−
— a projection de dS sur le plan ⊥à Y est : Équation des moments T x = τxy i +
→
− →
−
σyy j + τzy k
→
− →
− →
− →
−
— a projection de dS sur le plan ⊥à Z est : T x = τxz i + τyz j + σzz k
→
− →
− →
− →
− − →
− →
−
Alors : T n = Tx → −
n x + Ty →
−
n y + Tz →
−
n z = (σxx i + τyx j + τzx k )→n x + (τxy i + σyy j +
→
− − →
− →
− →
− −
τzy k )→
n y + (τxz i + τyz j + σzz k )→nz
  
 σxx τyx τzx  nx 
→
−  n =⇒→ −
T n = T .→
−
  
on peur écrire sous forme matricielle : T n =  τ
 xy σ yy τzy  y  n
  
τxz τyz σzz nz
99K c’est le tenseur des contraintes.
d’où : on peut le décomposer par la sommation de deux matrices, un tenseur sphérique

et un autre de trace nulle.
   
0
 1 0 0   σxx τyx τzx 
0   0  
 0 1 0  et T =  τxy σyy τzy →c’est le tenseur des
=⇒T = α.I + T teq :I=    0 
   
0
0 0 1 τxz τyz σzz
contraintes visqueuses.

0
σxx = α + σxx





0 0 0 0
avec : σyy = α + σyy et σxx + σyy + σzz = 0 = −3α + σxx + σyy + σzz


 σzz = α + σ 0


zz

1
σxx + σyy T + σzz ⇒ α = 31 T r T

⇒α= 3
→
− 0 0
On écrit : T n = T .→
−
n = α.I.→
−
n + T .→
−
n = α.→
−
n + T .→
−
n
→
− →
− 0 0
donc : d F = T n .dS = α.→
−
n .dS + T .→
−
n .dS = −P.→
−
n .dS + T .→
−
n .dS (où : le 1er terme
présente les forces de pression hydrostatique).
0
d’où : P = 31 T r T et T r(T ) = 0 avec : α = −P c’est la pression hydrostatique.
5.4 Équation de Navier-Stokes
5.4.1 Notion du tenseur de déformation
Il faut voir aussi le tenseur de déformation G ?
soit un élément de fluide dV en déplacement de M vers M 0 , le champ de vitesses

dérive de ce qui suit :
→
−0 →− →
− →
− − →
− − →
− −
V = V + d V = V (→ r + d→−r ) = V (→r ) + V (d→r ) et on écrit :

u0 = u + ∂u dx + ∂u dy + ∂u dz



 ∂x ∂y ∂z

v 0 = v + ∂x
∂v ∂v
dx + ∂y dy + ∂v
∂z
dz —— Équation des moments



 w0 = w + ∂w dx + ∂w dy + ∂w dz

∂x ∂y ∂z
      
0 ∂u ∂u ∂u
 u   u   ∂x ∂y ∂z   dx 
      
⇒
 v 0  = 
  v +
 
∂v
∂x
∂v
∂y
∂v
∂z
  dy 
 
      
w0 w ∂w
∂x
∂w
∂y
∂w
∂z
dz
on note :  
∂u ∂u ∂u
 ∂x ∂y ∂z 
 
G= ∂v ∂v ∂v  (5.2)
 ∂x ∂y ∂z 
 
∂w ∂w ∂w
∂x ∂y ∂z
G : est la matrice Jacobienne ou le Jacobien
→
− − →
− −
=⇒ V (→
r + d→
− r ) + G.d→
r ) = V (→ −
r
On appelle G le tenseur de déformation, qui représente quatre types de déformation

(translation pure, élongation/contraction, déformation angulaire et rotation pure), comme
présenté sur la figure ci-contre.
Figure 5.1 – Cas de déformations
5.4.2 Décomposition du tenseur de déformation

=
on écrit : G = e + w
 
∂u 1 ∂u ∂v 1 ∂u ∂w
 ∂x 2 ∂y
+ ∂x 2 ∂z
+ ∂x
=  

1 ∂u ∂v ∂v 1 ∂v
avec : e = 
 2 ∂y
+ ∂x ∂y 2 ∂z
+ ∂w
∂y
 c’est un tenseur symétrique

 
1 ∂u ∂w 1 ∂v ∂w ∂w

2 ∂z
+ ∂x 2 ∂z
+ ∂y ∂z
 
1 ∂u ∂v 1 ∂u ∂w
 0 − 2
−∂y ∂x 2 ∂z ∂x
 

et w =  1 ∂u ∂v 1 ∂v ∂w  est un tenseur anti-symétrique.
 2
− ∂y
+ ∂x
0 2 ∂z
− ∂y 
 
1 ∂u ∂w 1 ∂v ∂w

2 ∂z
− ∂x 2 ∂z
− ∂y
0
 
 0 G (d→
−
r ) − Ωz Ωy 
 
On note : w = 
 Ωz 0 −Ωx 

 
−Ωy Ωx 0
  
 0 −Ωz Ωy   dx 
teq : w.d→
−   
r =
 Ωz 0 −Ωx 
  dy 
 
  
−Ωy Ωx 0 dz
 
 0 −Ωz dy Ωy dz 
 →
⇒ w.d→
− =− →
−

r =
 Ωz dx 0 −Ωx dz  Ω ∧dr
 
−Ωy dx Ωx dy 0
Remarque : On définie le mouvement et la déformation d’une particule fluide en termes

de simple translation, élongation/contraction, déformation angulaire et rotation pure, tout
→
−
en développent l’expression de variation de vitesse (d V ) :
→
− − →
− →
= V (→ →
− − = →− →
−
r ) + G.d
| {z r} = V ( r ) + e.d r + w.d r
−
→
dV
→
− − →
− − →
−
⇒ V (→
r + d→
−
r ) = V (→ →
− →
−
=
r) + |e.d
{zr} + |Ω ∧
{zd r}
| {z }
translation déf ormation rotation pure
5.4.3 Équation de Navier-Stokes
Par définition d’un fluide est dit Newtonien , on a le tenseur de contraintes visqueuses
0 =
T qui dépend linéairement du tenseur de déformation e .
=

∂ui ∂uj
e est symétrique −→ eij = 21 ∂x j
+ ∂xi
⇒ eij = eji
0
On a la relation entre les éléments du tenseur T et e
pour un fluide isotrope :
σij0 = 2µeij + µ0 (exx + eyy + ezz ) δij 99K (∗)
avec, δij : symbole Kronecker et δi=j = 1 et δi6=j = 0
µ : viscosité dynamique.
µ0 : viscosité de dilatation.
Les éléments diagonaux du tenseur eij ne sont que la divergence de la vitesse,

3 3
∂ui →
−
c’est à dire : e11 + e22 + e33 = eii = ∂u ∂v ∂w
P P
∂x
+ ∂y
+ ∂z
= ∂xi
= ∇V
i=1 i=1
→
−
pour un fluide incompressible : ∇ V = 0,
donc l’équation (*) devient : σij0 = 2µeij
car : exx + eyy + ezz = 0

0 =
Alors, les éléments de T et e} sont directement proportionnels et le coefficient de pro-
portionnalité est la viscosité dynamique du fluide.
0 =
d’où : T = 2µ e
D’après la seconde loi de la dynamique (Newton) et avec l’hypothèse de fluide Newtonien,

incompressible et isotrope :
→
− →
− →
− 0 −
→
on a : Σ F ext = m.→
−
a ⇒ −dV ∇P + dV ∇T + ρdV →
−
g = ρdV dV
dt
→
− →
−= −
→
⇒ − ∇P + 2µ ∇ e + ρ→
−
g = ρ ddtV
→
−=
et comme ∇ eest la notation compacte du tenseur de déformation, qui est égale à :
! " #
→
− = P P ∂eij →
.− .→
−

P P 1 ∂ ∂Vi ∂V
∇e = ∂xj
ei = 2 ∂xj ∂xj
+ ∂xji ei
i j i j
! !
→
−=
.→
− .→
− .→
−
PP ∂ 2 Vi ∂ 2 Vj
P P ∂ 2 Vi P ∂Vi P ∂Vj
Enoncé⇒ ∇ e = 12 ∂x2
+ ∂xi ∂xj
ei = 1
2 ∂x2j
ei + 1
2 ∂xi ∂xj
ei
j
i j i j i j
!
.→
−
1
P ∂Vi P ∂Vj
le second terme est nul : 2 ∂xi ∂xj
ei = 0 ;
i j
car la divergence de vitesse est nulle pour un écoulement incompressible :
P ∂V
( ∂xjj = ∇Vj = 0).
j
et le premier terme est égale à Laplacien de vitesse :
P ∂ 2 Vi
∂x2j
= ∆Vj
j
!
→
−=
.→
−
ei = 21 (∆Vj ) .→
−
1
P P ∂ 2 Vi P
donc : ∇ e = 2 ∂x2j
ei
i j i
et pour un fluide incompressible : ∇Vj = 0, alors :
→
− = 1P →
−
∇ e = 2 (∆Vj ) .→
−
ei = 21 ∆ V
i
Donc, l’équation de la 2nd loi fondamentale de la dynamique s’écrira :

→
− →
− −
→
− ∇P + µ∆ V + ρ→
−
g = ρ ddtV
→
− →
− −
→ →
−→− →
−
− ∇P + µ∆ V + ρ→
−
g = ρ ∂∂tV + ρ V ∇ V 99KC’est l’équation de Navier-Stokes
−
→ −
→ →
−→− →
−
dV ∂V
avec : dt
= ∂t
+ V ∇ V
Donc, la projection de l’équation de Navier-Stokes sur les trois axes, on obtient le système
d’équations en coordonnées cartésiennes :
 
∂2u ∂2u ∂2u



 ρ ∂u
∂t
+ρ u ∂u
∂x
+ v ∂u
∂y
+ w ∂u
∂z
= − ∂P
∂x
+µ ∂x2
+ ∂y 2
+ ∂z 2




 
∂2v ∂2v ∂2v
ρ ∂v
∂t
+ ρ u ∂v
∂x
+ v ∂v
∂y
+ w ∂v
∂z
= − ∂P
∂y
+ µ ∂x2
+ ∂y 2
+ ∂z 2
(5.3)

 2 

2 2
 
ρ ∂w + ρ u ∂w
+ v ∂w
+ w ∂w
= − ∂P
+ µ ∂ w
+ ∂ w
+ ∂ w
− ρg

 

∂t ∂x ∂y ∂z ∂z ∂x2 ∂y 2 ∂z 2
Le système d’équations de Navier-Stokes en coordonnées cylindriques,
avec : →
−
u = (ur , vθ , w) et 1 ∂(rur )
r ∂r
+ 1 ∂vθ
r ∂θ
+ ∂uz
∂z
=0
 h 2 i 
vθ 2 2 2 ∂vθ
ρ ∂t + ρ ur ∂r + vθ r ∂θ − r + w ∂z = − ∂r + µ ∂∂ru2r + 1r ∂u
∂ur
∂ur 1 ∂ur ∂ur ∂P
− ur2r + r12 ∂∂θu2r + ∂∂zu2r −
 r 
r2 ∂θ

 ∂r



 2 

2 2
ρ ∂v
∂v 1 ∂vθ ∂vθ ∂ vθ
ur 1 ∂P
+ 1r ∂v − vrθ2 + r12 ∂∂θv2θ + ∂∂zv2θ + 2 ∂ur

∂t
θ
+ ρ ur ∂r
θ
+ v θ r ∂θ
− r
+ w ∂z
= − r ∂θ
+ µ ∂r2 ∂r
θ
r2 ∂θ

 2 

1 ∂2w ∂2w
 vθ ∂w

ρ ∂w ∂w ∂w ∂P ∂ w 1 ∂w

+ ρ u + + w = − + µ + + + − ρg

 

∂t ∂r r ∂θ ∂z ∂z ∂r 2 r ∂r 2
r ∂θ 2 ∂z 2
(5.4)
Le système d’équations de Navier-Stokes en coordonnées sphérique,
avec :→
−
u = (ur , uθ , uϕ ) et ∂ur
∂r
+ 2ur
r
+ 1 ∂uθ
r ∂θ
+ vθcotθ
r
+ 1 ∂uϕ
rsinθ ∂ϕ
=0
Sur l0 axe r :
u2ϕ
h i
uθ ∂ur uϕ ∂ur u2θ
ρ ∂u
∂t
r
+ ρ ur ∂u
∂r
r
+ r ∂θ
+ rsinθ ∂ϕ
− r
− r
= − ∂P
∂r
+ fr
h i
∂ 2 ur 2 ∂ur 2ur 1 ∂ 2 ur cotθ ∂ur 1 ∂ 2 ur 2 ∂vθ 2uθ cotθ 2 ∂uϕ
+µ ∂r2
+ r ∂r
− r2
+ r2 ∂θ2
+ r2 ∂θ
+ r2 sin2 θ ∂ϕ2
− r2 ∂θ
− r2
− r2 sinθ ∂ϕ
Sur l0 axeθ :
u2ϕ cotθ
h i
uϕ ∂uθ
ρ ∂u
∂t
θ
+ρ ur ∂u
∂r
θ
+ ur uθ
r
+ uθ ∂uθ
r ∂θ
+ rsinθ ∂ϕ
− r
= − 1r ∂P
∂θ
+ fθ
h i
∂ 2 uθ 2 ∂uθ uθ 1 ∂ 2 uθ cotθ ∂uθ 1 ∂ 2 uθ 2 ∂ur 2cosθ ∂uϕ
+µ ∂r2
+ r ∂r
− r2 sin2 θ
+ r2 ∂θ2
+ r2 ∂θ
+ r2 sin2 θ ∂ϕ2
+ r2 ∂θ
− r2 sin2 θ
∂ϕ
Sur l0 axeϕ :
h i
ρ ∂u
∂t
ϕ
+ ρ ur ∂u
∂r
ϕ
+ ur uϕ
r
+ uθ ∂uϕ
r ∂θ
+ uθ uϕ cotθ
r
+ uϕ ∂uϕ
rsinθ ∂ϕ
1 ∂P
= − rsinθ ∂ϕ
+ fϕ
h i
∂ 2 uϕ 2 ∂uϕ uϕ 1 ∂ 2 uϕ cotθ ∂uϕ 1 ∂ 2 uϕ 2 ∂ur 2cosθ ∂uθ
+µ ∂r2
+ r ∂r
− r2 sin2 θ
+ r2 ∂ϕ2
+ r2 ∂θ
+ r2 sin2 θ ∂ϕ2
+ r2 sinθ ∂ϕ
+ r2 sin2 θ ∂ϕ
5.5 Pertes des charges (application aux fluides réels)
Il existe deux types de pertes de charges, les pertes de charges linéaires liées à la longueur
et les pertes singulières liées aux changement de formes du circuit de l’écoulement du
fluide ( variation des section de passage du fluide).
5.6 Pertes de charge linéaires
Les pertes de charge linéaires s’expriment par l’équation (1) donnée Équation des moments
en fonction du coefficient de pertes de charge linéaire (λ) :
v2

D
JL12 = −λ (5.5)
2g L
avec :
D : Diamètre hydraulique,
L : Longueur de canalisation,
v : vitesse moyenne d’écoulement,
g : accélération de pesanteur,
λ :coefficient de perte de charge linéaire, qui est lié aussi au coefficient de frottement (f :
friction factor), donné par l’équation ci-après :
λ
f= 4
2τ ∂u
et f = ρv 2
avec : τ = µ ∂n n :normale à la direction d’écoulement.
λn’est fonction que du nombre de Reynolds et de la rugosité relative de la paroi (ε/D),

qui peut être obtenue du diagramme de Moody.
On distingue différent cas :
- Si Re <2300 : Cas du régime laminaire, appelé aussi droite de Poiseuille, où λn’est
fonction que du nombre de Re :
Re
λ= 64
- si 2300<Re>10000 : cas d’écoulement turbulent et paroi lisse, où λÉquationdesmomentsne

dépend encore que de Re et on utilise la droite de Blasius :
Équation des momentsλ = 0.316 × Re−0.25
√
ou encore l’équation de Van Karman : √1 = 2log Re λ − 0.8
λ
qui se résout par méthode itérative (Newton) avec : λ0 = 0.02 comme solution initiale.
- Si Re>10000 : Cas de régime d’écoulement turbulent et rigoureux, alors λne dépend plus
du nombre de Reynolds. D’après le diagramme de Moody, en très peu de cas la courbe
ε/D est presque une droite horizontale. Pour les calculs industriels o, utilise la relation
de Blench :
λ = 0.790 Dε 99K équation de Blench

p
Remarque : pour les conduites expérimentales de rugosité uniforme, on utilise la relation

de Karman Nikuradse :
√1 D

λ
= 2log 2ε
+ 1.74
Loi générale du coefficient de pertes de charge en régime turbulent et paroi lisse ou rigou-
reuse, ainsi que pour le transitoire. Colebrook et White ont regroupé les lois de Prandtl
et Karman, en proposant l’équation suivante :

√1 = −2log 2.51
√ + ε
λ Re λ 3.71D
5.6.1 Application (mesure de la viscosité) :
Dans un large réservoir, on place un liquide de masse volumique ρ. On fixe au fond du

récipient un long tube fin horizontal (longueur L rayon R) ; les forces de viscosités sont
prépondérantes ici, l’écoulement est très lent (goutte à goutte). La mesure du débit permet
de calculer la viscositéµ.
En appliquant Bernoulli, dans le réservoir (les forces de viscosité n’interviennent pas là)
on obtient :
Figure 5.2 – Mesure de la viscosité
4P = ρgh − 12 ρv 2 ≈ ρgh.
πR4 ρgh πR4 ρgh

La loi de poiseuille donne : Q = 8µ L
⇒µ= 8Q L
par exemple, pour l’huile de masse volumique ρ = 860kg/m3 , h = 0.3m, R = 2mm, L =

0.15m et on mesure un débit : Q = 53ml/min. D’où :
4
π (2.10−3 ) 860×9.81×0.3
µ= 8×(53.10−6 /60) 0.15
= 0.12 [kg/m2 .s] = 0.12 [P a.s]
5.7 Pertes de charge singulières
Ce type de pertes de charge est lié à des accidents ponctuels (changement de forme de
tuyauteries), par exemple l’élargissement et le rétrécissement brusques ou progressifs,
coudes, formes U, T et Y , vannes, robinets, clapets ...etc.
Elles s’expriment à partir d’un coefficient de pertes de charge adimensionnel noté par
e et donnée par la formule suivante :
2
JS12 = −e v2g
avec :
e : coefficient des pertes de charge singulière,
v : vitesse moyenne d’écoulement à la section de passage considérée,
g : accélération de pesanteur,
Le coefficient e est déterminé explicitement pour quelques cas (ou bien par formule
empiriques), ou encoure par des abaques. on peut retenir quelques exemples des cas ci-
après :
2
S1
- Élargissement brusque : e = 1 − S2
2
S2
- Rétrécissement brusque : e = 0.45 × 1 − S1
- Tube plongeant dans un bac :e = 1
- Coude 90° : e = 0.8
- Vanne ouverte : e = 1.2
5.7.1 Coefficient de contraction :
L’évaluation du débit de vidange nécessite de prendre en compte la contraction du jet.

Comme le montre l’encart de la figure (..) , la vitesse évaluée précédemment correspond à
une section plus faible que celle de l’orifice. En conséquence, le débit volumique de vidange
s’obtient en calculant :qv = UB × SB
où SB est la section du jet où les lignes de courant peuvent être considérées rectilignes et
parallèles. On peut ainsi définir un coefficient de contraction Cc = SB /S, lequel dépend
essentiellement du type de paroi ainsi que du profil de l’orifice dans la paroi. On peut
alors reformuler le débit de vidange ainsi : :qv = UB × S×C c
Exemple : Les figures ci-après donnent de manière non exhaustive quelques valeurs ty-
piques du ce coefficient de contraction.
Figure 5.3 – Coefficient de contraction [5]
Figure 5.4 – Diagramme de Moody[5]
5.7.2 Diagramme de Moody
Présentation et lecture du diagramme :
Utilisation en régime laminaire :
Le coefficient se lit directement à partir de la droite 64/Re, ( voir figure ci-dessous).

Le coefficient se lit directement à partir de la droite 64/Re
Figure 5.5 – Régime laminaire
Utilisation en régime turbulent :
Le coefficient se lit directement à partir de la zone délimitée par les deux profils en rouge,
( voir figure ci-dessous).
On calcule la rugosité relative ε et on sélectionne la courbe correspondante 0.02 ou 5.10−4

(ici, dans la figure ci-après). On détermine le nombre de Reynolds et on lit à l’intersection
de la courbe et de la verticale
— On voit qu’au-delà de la courbe Complète turbulence, le coefficient ne dépend

plus que de la rugosité et est indépendant du nombre de Reynolds.
— La ligne Smooth Pipe correspond à la limite du diagramme en régime turbu-
lent : les conduites ne sont plus rugueuses sur cette ligne.
On peut se demander quelle zone du diagramme est intéressante pour les écoulements
habituels dans les conduites horizontales : Si on fixe Re à 2000 (valeur critique laminaire
– turbulent), on peut calculer les vitesses critiques v *, au-delà desquelles le régime est
turbulent :
Figure 5.6 – Régime turbulent - a
Figure 5.7 – régime turbulent - b

V > v ∗ = 2000.ν/D
où : υ est la viscosité cinématique du fluide.
On voit bien que les vitesses critiques sont très inférieures aux vitesses usuellement ren-
contrées donc les régimes seront toujours turbulents.
Exercice d’application :
Un pipe-line de diamètre d = 25cm est de longueur L est destiné à acheminer du pétrole

brut dune station A vers une station B avec un débit massique qm = 18kg/s.
Les caractéristiques physiques du pétrole sont les suivantes :
- masse volumique ρ = 900kg/m3 ,
- viscosité dynamique µ = 0, 261P a.s.
On suppose que le pipe-line est horizontal.
1) Calculer le débit volumique qv du pétrole.

→
−
2) Déterminer sa vitesse d’écoulement V .
3) Calculer le nombre de Reynolds Re.
4) Quelle est la nature de l’écoulement ?

5) Calculer la valeur du coefficient des pertes de charge linéaire.
6) Exprimer la relation de Bernoulli entre A et B.
Préciser les conditions d’application et simplier.
7) Déterminer la longueur L maximale entre deux stations A et B à partir de laquelle la

chute de pression (4P = PA − PB ) dépasse 3 bar.
Réponse :
Chapitre 6
Analyse dimensionnelle et similitude
6.1 Introduction
Les méthodes analytiques (ou modèles) utilisées pour résoudre les problèmes de mécanique
des fluides ne sont pas toujours satisfaisantes dans la pratique (souvent des simplifications
sont nécessaires ou des analyses détaillés couteuses sont indispensables). La méthode alter-
native est d’utiliser l’expérience et d’en déduire des corrélations applicables à l’ensemble
des cas de problèmes du même type. Comme les conditions des expériences effectuées au
laboratoire peuvent être différentes d’un cas réel à l’autre (par exemple, les dimensions
de la maquette d’un avion et du prototype, les fluides différents...).
La solution adéquate est d’utiliser une méthode d’analyse dimensionnelle qui permet de
d’obtenir des corrélations sans dimensions, applicables dans des conditions pratiques de
façon quasi-universelle. C’est à dire ces corrélations sans dimensions peuvent être ap-
pliquées dans des conditions dynamiques similaires à celle dans les quelles elles ont été
établies en utilisant, par exemple un fluide différent.
6.2 Procédure de l’analyse dimensionnelle
1. établir une liste de variable pertinentes du problème considéré dépendantes ou

indépendantes.
2. Identifier le nombre et la forme des paramètres adimensionnels.
63
CHAPITRE 6. ANALYSE DIMENSIONNELLE ET SIMILITUDE 64
Une de ces technique consiste à appliquer le théorème de Vaschy-Buckingham dit aussi

théorème Pi.
6.3 Théorème de Vaschy-Buckingham
Énonce : Une équation physique complète de la forme générale :
f (q1 , q2 , ..., qn ) = 0
où qi représente les n variables physiques choisies pour décrire le problème, exprimées en
terme de j unités physiques interdépendantes, peut être réécrite sous la forme :
F (Π1 , Π2 , ..., Πk ) = 0
où Πi sont les nombres sans dimension construit à partir des qi par k = n − j équations
de la forme :
= q1m1 q2m2 ...qnmn

Q
i
où les mi sont des constantes.
En pratique soit :
n : le nombre de paramètres indépendants du problème.
j 0 : le nombre de dimensions de base (M, L, T, θ...etc) exprimant les unités de n paramètres.
j : le nombre de dimensions de base correspondant aux paramètres répétés qui ne forme

pas un produit.
Q
k : le nombre de produit indépendants permettant de décrire le problème, soit : k =
n − j.
Les étapes de démarche à suivre sont donc :
1. Lister les n paramètres du problème.
2. Exprimer les dimensions de chaque paramètre en utilisant les dimensions de base

(M, L, T, θ), compter le nombre des dimensions de base utilisés j 0 dans l’ensemble
de ces paramètres.
3. Trouver le nombre j en posant initialement j=j’. En suite chercher les paramètres

Q
répétés qui ne forme pas un produit . Si c’est impossible réduire j un par un et
répéter la procédure.
Q
4. Choisir j paramètres répétés qui ne forment pas le produit .
Q
5. En considérant les paramètres non répétés, un par un, former les produits en
mettant ensemble les paramètres répétés et non répétés. En déduire algébriquement
Q
les puissances de chaque paramètre répété pour construire les produits sans
dimension.
Q
6. Écrire la combinaison des produits trouvés et les mettre sous la forme :Πk =
F (Π1 , Π2 , ...Πi ).
6.3.1 Exemple d’application :
considérant le problème d’un fluide visqueux dans un conduit cylindrique. On cherche la

chute de pression sans dimension en fonction des autres paramètres sans dimension.
1- Les paramètres du problème sont :
4P : chute de pression.
V : vitesse moyenne de l’écoulement.
D : diamètre de conduit.
l : longueur du conduit.
ρ : masse volumique du fluide.
µ : viscosité dynamique du fluide.
ε : rugosité de la parois du conduit.
Donc, on a 4P est un paramètre dépendant et le reste (V, D, l, ρ, µ, ε) sont des paramètres

indépendants. alors, le nombre des paramètres du problème est n=7.
2- La liste des paramètres et leurs dimensions est présentée dans le tableau ci-dessous :
Paramètres 4P V D l ρ µ ε
dimensions M L−1 T −2 LT −1 L L M L−3 M L−1 T −1 L
⇒le nombre de dimension de base est : j 0 = 3.

3- Choisir : j = j 0 = 3 et les paramètres répétés comme (V, D, ρ). Cependant ces pa-
Q
ramètres ne peuvent pas former un produit sans dimension (aucune combinaison ne
peut éliminer la masse de la densité et le temps de la vitesse).
Q
4- V, D, ρ sont les paramètres répétés qui ne forme pas le produit (étape identique à
3).
Donc, k = n − j = 7 − 3 = 4. Alors, il existe 4 produit indépendants.

Q
5- construction des produits :
−1 a c 1
a b c 1
) (L)b (M L−3 ) (M L−1 T −1 )
Q
1 = V D ρ µ = (LT
NB : Choisissons µ comme paramètre non répété.

Q
Pour que le produit 1 soit sans dimension, il faut que la puissance de chaque dimension
de base soit égale à zéro.
Ainsi, on a :
pour M : c + 1 = 0 ⇒ c = −1
pour L : a + b − 3c − 1 = 0 ⇒b = −1
pour T : −a − 1 = 0 ⇒ a = −1
on a donc : 1 = ρVµD ,
Q
Q ρV D
on peut l’écrire sans perdre la généralité comme : 1 = µ
= Re ; qui est le nombre de
Reynolds.
* Choisir l comme paramètre non répété (l : longueur de conduite et non unité de base).
−1 a c
a b c 1
) (L)b (M L−3 ) (L)1
Q
2 = V D ρ l = (LT
Alors,
pour M : c = 0
pour L : a + b − 3c + 1 = 0 ⇒b = −1
pour T : −a = 0
l
Q
on a donc : 2 = D
: le rapport longueur / diamètre
6.4 La signification du rôtationnel
, c’est un facteur de forme.

* Choisir ε comme paramètre non répété :

−1 a c
a b c 1
) (L)b (M L−3 ) (L)1
Q
3 = V D ρ ε = (LT
d’où : a = c = 0 et b = −1
alors : 3 = Dε : c’est la rugosité relative.

Q
* Choisissons finalement 4P comme paramètre non répété :

−1 a c 1
a b c 1
) (L)b (M L−3 ) (M L−1 T −2 )
Q
4 = V D ρ 4P = (LT
pour M : c + 1 = 0 ⇒ c = −1
pour L : a + b − 3c − 1 = 0 ⇒b = 0
pour T : −a − 2 = 0 ⇒ a = −2
Q 4P
alors : 4 = ρV 2 : c’est le coefficient de pression.

4P ρV D l ε
ou Cp = f Re, ¯l, ε̄

6- combinaison des paramètres : ρV 2 = f µ
, ,
D D
le coefficient de chute de pression d’un écoulement visqueux dans une conduite ayant une
rugosité est une fonction de :
- nombre de Reynolds.
- longueur sans dimension du conduit.
- rugosité sans dimension du conduit.
6.5 Similitude
Considérons deux problèmes différents par exemples deux écoulements visqueux distinctes
(ρ1 , µ1 , l1 , V1 ) et (ρ2 , µ2 , l2 , V2 ) concernant le même phénomène (ici il s’agit de la viscosité).
Les problèmes sont dits semblables si leurs équations sans dimension (équations réduites)
sont identiques et si les conditions aux limites et initiales des problèmes exprimées sont
identiques.
Conséquence : deux problèmes différents sont dits semblables si les paramètres sans di-
mension ou paramètres de similitude sont les mêmes, deux à deux, et si les conditions
aux limites réduites sont les mêmes.
On distingue trois types de similitude en mécanique des fluides : similitude géométrique,

similitude cinématique et similitude dynamique.
Figure 6.1 – profil d’aile
6.5.1 Similitude géométrique
L’échelle linaire est la même, les dimensions linéaires du modèle réduit correspond à
celles du prototype définie pzr un facteur d’échelle constant. Exemple : considérant le
profil d’aile.
rm Lm
si :SFk est le facteur d’échelle :SFk = rp
= Lp
=...
Dans ce type de similitude on conserve les mêmes angles, les mêmes directions d’écoulement,
la même orientation par rapport à l’environnement (angle d’incidence identique).
6.5.2 Similitude cinématique
L’échelle de vitesse est la même, les vitesses aux points correspondants sur le modèle et
le prototype sont dans la même direction. Elles sont différentes seulement par un facteur
d’échelle constant β.
Conséquence : les écoulements ont le même motif de lignes de courant et le même

régime d’écoulement.
exemple : écoulement autour d’un objet solide.
Les vitesses en chaque point sont proportionnelles par un facteur d’échelle constant β.
Figure 6.2 – Écoulement autour d’objet solide

Vp∞ = βVm∞ ; Vp1 = βVm1 ; Vp2 = βVm2
6.5.3 Similitude dynamique
La similitude dynamique est réalisée si le rapport des forces aux points homogènes sur le
modèle et le prototype est identique à un facteur d’échelle constant.
Exemple : écoulement à travers une porte d’écluse permettant à un bateau de passer d’un
fleuve ç un bassin ou l’inverse.
La similitude dynamique se produit dans les conditions suivantes :
1. Écoulement incompressible sans surfaces libres : Rem = Rep
2. Écoulement incompressible avec surfaces libres : Rem = Rep et F rm = F rp
3. Écoulement incompressible : Rem = Rep ; M am = M ap et γm = γp .
4. Écoulement avec tension surfacique :Rem = Rep ; W em = W ep .

6.6 Régimes d’écoulement et paramètres de simili-
tude
6.6.1 Nombre de Reynolds
Considérons un écoulement de fluide visqueux le long d’une plaque plane de longueur L.

On suppose les forces de volumes négligeables devant les forces visqueuses.
L’expérience montre que si la parois est imperméable, le champ de vitesse est tangent à
la paroi et que la vitesse diminue près de la paroi et s’annule sur la paroi si celle-ci est
immobile. On suppose que le champ de vitesse est parallèle à l’axe des X et la viscosité
dynamique est fonction de la vitesse : µ = U (y).
L’équation de bilan de quantité de mouvement s’écrit dans ce cas :
∂ 2u

−
→ −−→ →
− →
−
ρ ddtV = −gradP + ρ f + µ x
|{z} ∂y 2
:≈0 | {z }
f orce de viscosité
par unité volume
Adimensionnons cette équation de la manière suivante :
u v P
Lu = U0
,v = U0
,P = P0
x = Lx , y = Ly , t = t
t0
où U0 , P0 , t0 et L sont choisis comme grandeurs de référence.

Bibliographie
Livres :
[1] Sakir Amiroudine & Jean-Luc Battaglia, Mécanique des fluides , 2nd édition, DU-
NOD, Paris 2014
[2] Jean-Francois Sini, cours de MDF , HAL Id : cel-00356205, 2009
Supports de cours :
[3] Laid MESSAOUDI, MDF Polycopie de cours, Université de Batna2, 2016
[4] Christophe Ancey, Mécanique des fluides, introduction , Notes de cours : version
12.2 du 2 juin 2016
[5] M’hamed BERIACHE, MÉCANIQUE DES FLUIDES I, Université de Chlef, 2016
[6] Philippe Marty, MECANIQUE DES FLUIDES : Maste 1 GDP ,Université Joseph
Fourier, Grenoble, 2012
[7] Daniel Huilier, Écoulements en conduites , polycopie de cours 2010
[8] Riadh BEN HAMOUDA, NOTIONS DE MECANIQUE DES FLUIDES, Centre de

Publication Universitaire, 2008
[9] METICHE MEHDI, MÉCANIQUE DES FLUIDES APPROFONDIE , Centre Uni-

versitaire de Béchar, 2004
Sites internet :
schémas tirés de sites :

[10] https ://pbs.twimg.com/media/EUcF1HiWkAASM6K.png
[11] https ://encrypted-tbn0.gstatic.com/images ?q=tbn :ANd9GcQP52OH9F2qPX7BZ AlkS8k0i5h7yqMlfWk2Q&usqp=CAU
71

Polycopie MDF App Bonne-Version Preliminaire

Transféré par

Droits d'auteur :

Formats disponibles

Polycopie MDF App Bonne-Version Preliminaire

Transféré par

Informations du document

Description originale:

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Polycopie MDF App Bonne-Version Preliminaire

Transféré par

Droits d'auteur :

Formats disponibles

1

Ministère de l’Enseignement Supérieur et de la Recherche Scientifique

Université de Batna 2 - Mustapha Benboulaı̈d

Département de Génie Mécanique

Destiné au parcours Master ERM / EN

Mécanique des fluides approfondie

Dr. Nabil BESSANANE

Année Universitaire : 2019 - 2020

1 Rappel d’analyse vectoriel 8

1.1 Les systèmes de coordonnées courants . . . . . . . . . . . . . . . . . . . . . 8

1.2 Champs scalaire et vectoriel . . . . . . . . . . . . . . . . . . . . . . . . . . 9

1.2.1 Champ scalaire : . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9

1.2.2 Champ vectoriel : . . . . . . . . . . . . . . . . . . . . . . . . . . . 9

1.3 Comment utiliser Nabla ? . . . . . . . . . . . . . . . . . . . . . . . . . . . 10

1.4 La signification du Gradient . . . . . . . . . . . . . . . . . . . . . . . . . . 11

1.5 La signification de la divergence . . . . . . . . . . . . . . . . . . . . . . . . 11

1.6 La signification du rôtationnel . . . . . . . . . . . . . . . . . . . . . . . . . 12

1.7 La signification du Laplacien . . . . . . . . . . . . . . . . . . . . . . . . . . 13

1.8 Quelques relations importantes . . . . . . . . . . . . . . . . . . . . . . . . 14

2 La physique des fluides 16

2.1.1 Échelle mésoscopique . . . . . . . . . . . . . . . . . . . . . . . . . . 17

2.1.2 Champ des vitesses dans un fluide . . . . . . . . . . . . . . . . . . 17

2.2 Équation fondamentale de l’hydrostatique . . . . . . . . . . . . . . . . . . 18

2.2.1 Forces de pression d’un fluide en repos . . . . . . . . . . . . . . . . 19

2.2.1.1 Cas de contrainte normale . . . . . . . . . . . . . . . . . . 19

2.2.1.2 Cas de contrainte tangentielle . . . . . . . . . . . . . . . 21

2.2.1.3 Variation de force de pression . . . . . . . . . . . . . . . . 22

2.2.2 Forces de pression d’un fluide en mouvement (uniforme) . . . . . . 23

2.2.3 Forces de pression d’un fluide en mouvement (quelconque) . . . . . 24

2.2.4 Application de de l’équation de BERNOULLI (approche statique) . 27

2.2.4.1 Effet de gravitation sur un liquide . . . . . . . . . . . . . 27

2.3 Théorèmes de transport . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31

2.3.1 Volumes et surfaces de contrôle . . . . . . . . . . . . . . . . . . . . 31

2.3.2 Théorème de transport . . . . . . . . . . . . . . . . . . . . . . . . 31

3 Équations d’Euler (fluides parfaits) 33

3.1 Rappel de cinématique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33

3.1.1 Approche de Lagrange . . . . . . . . . . . . . . . . . . . . . . . . . 33

3.1.2 Approche de Euler . . . . . . . . . . . . . . . . . . . . . . . . . . . 34

3.2 Dynamique des fluides parfaits . . . . . . . . . . . . . . . . . . . . . . . . 35

3.2.1 Notion de dérivée particulaire : . . . . . . . . . . . . . . . . . . . . 35

3.2.2 Équation d’Euler : . . . . . . . . . . . . . . . . . . . . . . . . . . . 36

3.2.3 Équation de BERNOULLI : ( autres approches) . . . . . . . . . . 37

3.2.4 Écoulement irrôtationnel : . . . . . . . . . . . . . . . . . . . . . . . 39

4 Équation de continuité (conservation de matière) 41

4.1 Courant de matière . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41

4.2 Bilan de matière (conservation de masse) . . . . . . . . . . . . . . . . . . . 43

4.3 Flux de matière d’un volume fini "V" . . . . . . . . . . . . . . . . . . . . 44

4.4 Équation de continuité . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45

5 Équations de Navier-Stokes (fluides réels) 47

5.1 Généralités et définitions . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47

5.2 Régime d’écoulement . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47

5.3 Tenseur de contraintes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48

5.4 Équation de Navier-Stokes . . . . . . . . . . . . . . . . . . . . . . . . . . . 49

5.4.1 Notion du tenseur de déformation . . . . . . . . . . . . . . . . . . . 49

5.4.2 Décomposition du tenseur de déformation . . . . . . . . . . . . . . 50

5.4.3 Équation de Navier-Stokes . . . . . . . . . . . . . . . . . . . . . . . 51

5.5 Pertes des charges (application aux fluides réels) . . . . . . . . . . . . . . . 54

5.6 Pertes de charge linéaires . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54