COURS MDF

Cours MDF Aroua Aouadi Enit 2022-2023

 I. INTRODUCTION

Cours MDF Aroua Aouadi Enit 2022-2023

 Mécanique des fluides : Domaine d’applications

Presque tous les domaines sont concernés !

Génie civil Génie hydraulique

▪ Bâtiments ▪ Station de pompage/Pompes hydrauliques
▪ Ponts ▪ Turbine
▪ Digue de protection ▪ Écoulements à surface libre
▪ Conception et dimensionnement des ouvrages GC ▪ Écoulements en charge (dans les canaux)
▪ Architecture et urbanisme… ▪ Hydrologie
▪ Hydraulique maritime
▪ Océanographie

▪ Barrages (conception…)
▪ Assainissement

3
Cours MDF Aroua Aouadi Enit 2022-2023
 Mécanique des fluides : Domaine d’applications

Presque tous les domaines sont concernés !

Environnement Aéronautique Génie énergétique Génie nucléaire

▪ Dispersion atmosphérique ▪ Avions (mesure de vitesse…) ▪ Éolienne ▪ Réacteur nucléaire
▪ Pollution ▪ Soufflerie ▪ Domaine pétrolier
▪ Turbomachines
▪ Transfert thermique
(refroidissement, climatisation)
Génie des procédés Météorologie Physique fondamentale
▪ Prédiction ▪ Magma
▪ Procédés de traitement des eaux
▪ Prévision de temps et de climat ▪ Plasmas
▪ Station d’épuration ▪ Formation des étoiles
▪ Réacteur chimique ▪ Réchauffement climatique
▪ Écoulement granulaires
▪ Turbulence
Génie naval et côtière

4
Cours MDF Aroua Aouadi Enit 2022-2023
 Mécanique des fluides : Domaine d’applications

Presque tous les domaines sont concernés !

Électromécanique

Domaine médical
▪ Hémodynamique ▪ Écoulements sanguins ▪ Mesure de tension artérielle

Cours MDF Aroua Aouadi Enit 2022-2023

 I. RAPPELS MMC

Cours MDF Aroua Aouadi Enit 2022-2023

 Eléments d’analyse vectorielle
Produit scalaire
Le produit scalaire des deux vecteurs A et B s’écrit comme :
En coordonnées cartésiennes
A.B = A x B x + A y B y + A z B z
En coordonnées cylindriques
A.B = A r B x + A  B y + A z B z
Exemple d’application
Soit les deux vecteurs : u ( 6; −5 ) et v ( −2;9 )
u.v = 6  (−2) + (−5)  9 = − 57
Remarques :
Un vecteur A d’ordre 1 scalaire un vecteur B d’ordre 1 résulte un scalaire d’ordre 0.
A i Bi = A.B 7
Cours MDF Aroua Aouadi Enit 2022-2023
 Eléments d’analyse vectorielle
Produit scalaire
Norme / module d’un vecteur
A.A = A 2 = A x A x + A y A y + A z A z A= A 2x + A 2y + A 2z

Soit O et M deux points du vecteur A et O et N deux points du vecteur B tels que A = OM et B = ON. La
norme du vecteur A correspond à la distance OM et la norme du vecteur B correspond à la distance ON.

Le produit scalaire des deux vecteurs A et B s’écrit comme le produit

du module A fois le module B : A.B = ABcos
α est l’angle entre les deux vecteurs A et B.

2 2
A.A = A  A  cos (A, A) = A  cos 0 = A
2
A.A = A = A
2

8
Cours MDF Aroua Aouadi Enit 2022-2023
 Eléments d’analyse vectorielle
Produit vectoriel
Le produit entre deux vecteurs A et B donne un autre vecteur :

( ) ( )
Ax B = A y B z − A z B y i x + ( A z B x − A x B z ) i y + A x B y − A y B x i z

Ax B =

9
Cours MDF Aroua Aouadi Enit 2022-2023
 Eléments d’analyse vectorielle
Produit vectoriel
Le nouveau vecteur résultant du produit vectoriel AxB est normal au plan contenant A et B.
Le produit vectoriel AxB représente le produit du module du vecteur A fois le module du vecteur B, le
tout multiplié par le sinus de l’angle formé entre A et B. le vecteur résultant de AxB est à la fois
perpendiculaire au vecteur A et au vecteur B

A x B = ABsin 

Le produit vectoriel du vecteur A avec lui-même résulte un vecteur nul

Ax A = 0
10
Cours MDF Aroua Aouadi Enit 2022-2023
 Eléments d’analyse vectorielle
Gradient d’un scalaire
Le vecteur gradient noté par grad ou par le symbole nabla , s’écrit comme :
  
grad =  = e x + e y + e z
x y z
Le gradient d’un champ scalaire définit un champ vectoriel :

En coordonnées cartésiennes   
 x 
 
 f f f   
grad f = f = ex + ey + ez Soit f =
x y z y  f
 
  
 z 
 
(ex, ey, ez) sont les 3 vecteurs unitaires pour les coordonnées cartésiennes (x, y, z).
11
Cours MDF Aroua Aouadi Enit 2022-2023
 Éléments d’analyse vectorielle
Gradient d’un scalaire
Le gradient d’un champ scalaire définit un champ vectoriel :  
En coordonnées cylindrique  r 
 
f 1 f f 1  
grad f =  f = er + e + ez Soit  f =   f
r r  z r 
 
  
 z 
 
(er, e, ez) sont les 3 vecteurs unitaires pour les coordonnées cartésiennes (r, , z).
 
En coordonnées sphérique  r 
 
f 1 f 1 f  = 1  
grad f =  f = er + e + e Soit f   f
r r  r sin  r 
 
 1  
 r sin  
(er, e, e) sont les 3 vecteurs unitaires pour les coordonnées cartésiennes (r, , ).  
12
Cours MDF Aroua Aouadi Enit 2022-2023
 Éléments d’analyse vectorielle
Divergence d’un vecteur
La divergence d’un vecteur peut être exprimée comme le produit scalaire de  et d’un vecteur A.
La divergence d’un champ vectoriel définit un champ scalaire.
En coordonnées cartésiennes     
  x   Ax 
   
 A  A y  Az   
div A =  . A = x + + =   .  Ay  Avec A = A x ex + A y ey + A z ez
x y z y  
   
   
  z   A z 
En coordonnées cylindriques  
1  ( r A r ) 1  A  A z
div A =  . A = + + Avec A = A r er + A eθ + A z e z
r r r  z
En coordonnées sphériques
1  ( r 2 Ar ) 1  ( sin A  ) 1  A
div A =  . A = 2 + + Avec A = A r e r + A  e θ + A  e 
r r r sin  r sin  13
Cours MDF Aroua Aouadi Enit 2022-2023
 Éléments d’analyse vectorielle
L’opérateur Laplacien
Le Laplacien peut être calculé comme le produit scalaire de  avec lui-même, il sera noté dans ce cas  2.
Comme il peut être défini en combinant les opérateurs de divergence et de gradient sur un champ scalaire,
noté dans ce cas . On distingue :
Le Laplacien scalaire : s’applique à un champ scalaire et produit un scalaire :
2 a 2 a 2 a
Laplacien ( a ) = div ( grad ( a ) ) = a =  a = 2 + 2 + 2
2

x y z
Le Laplacien vectoriel : s’applique à un vecteur et produit un vecteur :
 2
A  2
A  2
A
Laplacien ( A ) = div ( grad ( A ) ) =  =  = + +
2
A A
x 2
y 2
z 2
 2Ax 2Ax 2Ax 
  x 2 +  y 2 + z 2 
 
A =  A =
2   2
A  2
A  2
Ay 

y
+ y
+ 2 
 x y z 
2 2

 2A  2
A  2
Az 
 z
+ z
+ 2 

14
Cours MDF Aroua Aouadi   x 2
 y 2
 z  Enit 2022-2023
 Éléments d’analyse vectorielle
L’opérateur Laplacien
En coordonnées cylindrique
Le Laplacien scalaire :

1   a  1  2 a  2 a
Laplacien ( a ) = div ( grad ( a ) ) =  a =  a =
2
 r r  + 2 + 2
r r   r 
2
z

Le Laplacien vectoriel :

 2 A r 2 A    2 A  2 A  
=
A =  A  2A r − − +
2  r 
e  A  − − 2  
e + (  2
Az ) ez
 r 2
r  
2
 r 2
r  
2

15
Cours MDF Aroua Aouadi Enit 2022-2023
 Éléments d’analyse vectorielle
Rotationnel d’un vecteur
Le rotationnel s’applique à un champ vectoriel résultant à son tour un champ vectoriel. Il est exprimé comme
le produit vectoriel  avec le vecteur A, noté x A .
Coordonnées cartésiennes
 A  Ay  A  A  Ay  Ax
rot A =  x A = ( z
− ) ex + ( x
− z
) ey + ( − ) ez
y z z x x y
Deux propriétés importante et faciles à vérifier sont le rotationnel d’un gradient ainsi que la divergence
d’un rotationnel sont nuls.
 x ( a) = 0
 .( x A) = 0
Coordonnées cylindriques
1 A
 A A A 1  (r A )  A r
rot A =  x A = ( z
− ) er + ( r
− z
) e + ( − ) ez
r  z z r r r 
Coordonnées sphériques
1  (sin  A )  A  1 1  A  (r A ) 1  (r A )  A r
rot A =  x A = ( − ) er + ( r
− ) e + ( − ) e
r sin    r sin   r r r  16
Cours MDF Aroua Aouadi Enit 2022-2023
 Éléments d’analyse vectorielle
Récapitulatifs
grad a =  a = ( a /  x,  a /  y,  a /  z)

div A = . A =  A x /  x +  A y /  y +  A z /  z)

rot A =  x A = ( A z /  y −  A y /  z ,  A x /  z −  A z /  x ,  A y /  x −  A x /  y)

2 a 2 a 2 a
Laplacien ( a ) = div ( grad ( a ) ) =  a = a = 2 + 2 + 2
2
x y z

 2
A  2
A  2
A
Laplacien ( A ) = div ( grad ( A ) ) =  A =  A = + +
2

x 2 y2 z 2

17
Cours MDF Aroua Aouadi Enit 2022-2023
 Éléments d’analyse vectorielle
Relations de base
Soit p et q deux scalaires et u et v deux vecteurs.
 .(p v) = ( v.  ) p + p (. v)
 .(p q) = q ( p) + p (q)
 .(u v) = (u.  ) v + v (. u)
div (rot u) =  .( x u) = 0
rot (grad p) =  x ( p) = 0
 x (p v) = (p) x v + p ( x v)
div (u x v) = rot u . v + rot v .u =  .(u x v ) = ( x u). v + ( x v ).u

 x (u x v) = u ( . v) − v ( . u) + ( v . ) u − (u . ) v
 x ( x v) =  ( . v) −  2 v) = rot (rot v) = grad (div v) − v

18
Cours MDF Aroua Aouadi Enit 2022-2023
 Éléments d’analyse vectorielle
Formules de Green
Formule du gradient

 grad f dV =  f n dS
V S

Formule d’Ostogradski

 div v dV =  v.n dS

V S

Formules de Stockes-Ampère

 rot v.n dS = 
S C
v . dl
Avec
S : la surface.
V : le volume.
C : le contour.
dl : élément du contour et tangent au contour en tout point.
n dS = dS (où n est le vecteur normale à la surface orienté vers l’extérieur).
19
Cours MDF Aroua Aouadi Enit 2022-2023
 Éléments de calcul tensoriel
Notion de tenseur

Dans un espace de dimension d (avec d = 1 ; 2 ou 3), on considère une quantité physique, décrite par
un tenseur de rang n, noté : a i jk...
Les indices n (i, j, k,…) varient entre 1 et d.
Le nombre de composantes de ce tenseur est égal à : dn

- Un tenseur d’ordre ou de rang 0 : zéro indices, c’est un scalaire de d 0 = 1 composante. Exemples : la masse
(M); la température (T)…

- Un tenseur d’ordre ou de rang 1 : 1 indice, il s’agit d’un vecteur. Soit le vecteur de la vitesse u, de
composantes u i avec i = 1, 2 et 3. On a donc d1 = 31 = 3 composantes.

20
Cours MDF Aroua Aouadi Enit 2022-2023
 Éléments de calcul tensoriel
Notion de tenseur
- Un tenseur d’ordre ou de rang 2 : 2 indices, il s’agit d’une matrice de d 2 = 3 2 = 9 composantes. On peut citer
u i
le gradient d’un vecteur dans 3 directions : G i j = . Il y aura donc 9 combinaisons de dérivées spatiales pour
x j
le gradient de la quantité G i j avec i = 1, 2, 3 et j = 1, 2, 3. Ex : fluctuations de vitesses.

u x /  x u x /  y u x / z 
G 11 G 12 G 13 
   
   i j  u y /  x
G =  G  = u y /  y u y / z  G = G 21 G 22 G 23 
u /  x  
 z u z /  y u z / z  G 31 G 32 G 33 

Un autre exemple de tenseur de rang 2, c’est le tenseur des contraintes  i j , défini comme la ième composante
de la force appliquée par unité de surface normale à n j.

- Un tenseur d’ordre ou de rang 3 : 3 indices, on aura d 3 = 33 = 27 composantes. Soit a i jk , i =1, 2, 3; j = 1,

2, 3 et k = 1, 2, 3.
- Un tenseur d’ordre ou de rang n : n indices, on aura d n composantes.
21
Cours MDF Aroua Aouadi Enit 2022-2023
 Éléments de calcul tensoriel
Quelques propriétés de tenseur
Tenseur symétrique
Un tenseur est symétrique pour 2 indices i et j s’il est égal à sa transposée
a i jk = a jik
Tenseur antisymétrique
Un tenseur est antisymétrique pour 2 indices i et j s’il est égal à l’opposé de sa transposée
a i jk = − a jik
Ceci n’est possible que si les termes diagonaux de a sont nuls : a 11 = a 22 = a 33 = 0
Tout tenseur de rang 2 peut s’écrire comme la somme d’un tenseur symétrique et d’un autre antisymétrique

a = a sym + a asym
a sym 1
(
= a + aT
2
)
Avec
a asym 1
2
(
= a − aT ) 22
Cours MDF Aroua Aouadi Enit 2022-2023
 Éléments de calcul tensoriel
Notations et conventions
Convention d'Einstein ou convention de sommation sur les indices répétés
n
a 1 b1 + a 2 b 2 + ........ + a n b n = a b
i =1
i i = a i bi

Indice muet
Tout indice répété deux fois dans le même monôme (terme) est dit indice muet ou sommé et donc on lui fera
prendre les valeurs 1, 2, 3 et on le fera la somme sur les trois dimensions de l’espace. La lettre utilisée pour
l’indice pour écrire la somme n’a pas d’importance dans ce cas et on peut la substituer par n’importe quelle
lettre : a i b i = a j b j = a k b k = a n b n......
Indice libre
Tout indice qui n’apparaissent qu’une fois dans chaque terme d’une équation est appelé un indice libre ou indice
franc, non sommé.
Exemple : a iki b j = S k jl T lnn
Les indices i, l, n sont muets (ou sommés), tandis que les indices k et j sont libres.
23
Cours MDF Aroua Aouadi Enit 2022-2023
 Éléments de calcul tensoriel
Notations et conventions
Convention d'Einstein ou convention de sommation sur les indices répétés
Autre Exemple
 a 11 a 12 a 13 
 
A =  a 21 a 22 a 23 
 
 a 31 a 32 a 33 
3
S1 = a 11 + a 12 + a 13 = a
j =1
1j = a 1j

3 3 3
S 2 = a 21 + a 22 + a 23 = a
j =1
2j = a 2j Si = a i1 + a i2 + a i3 = a = aij ik = a i j = a ik
j =1 j =1
3
S 3 = a 31 + a 32 + a 33 = 
j =1
a 3j = a 3j Indice libre : i
Indices muets : j et k

24
Cours MDF Aroua Aouadi Enit 2022-2023
 Éléments de calcul tensoriel
Notations et conventions
Symbole de Kronecker
Le symbole de Kronecker est un tenseur unité (ou tenseur identité), noté  i j , tel que :
 i j = 1 si i = j

 i j = 0 si i  j
Il s’agit d’un tenseur symétrique qui correspond à la matrice unité (ou matrice identité) :
1 0 0 
  = I = 0 1 0 
 
0 0 1
L’égalité x = I x correspond à x i =  i j x j
 i j permet de modifier la lettre utilisée pour l’indice d’un tenseur.

En utilisant la convention d’Einstein, on peut écrire :

 ii = 11 +  22 +  33 = 3 = d
25
Cours MDF Aroua Aouadi Enit 2022-2023
 Éléments de calcul tensoriel
Notations et conventions
Symbole Levi-Civita
On définit le tenseur antisymétrique de rang 3, noté  i jk (dit tenseur de Levi-Civita ou de permutation
circulaire) tel que :
0 si deux indices égaux

 i jk = +1 si ( i, j, k ) est une permutation paire de (1, 2, 3 )
−1 si i, j, k est une permutation impaire de 1, 2, 3
 ( ) ( )

- Les permutations paires : (1, 2 3), (2, 3, 1) et (3, 1, 2).

- Les permutations impaires : (2, 1, 3), (1, 3, 2) et (3, 2, 1).

Ce tenseur peut se représenter comme un cube de 33 = 27 éléments, dont seuls 6 sont non nuls. Ce tenseur
est utile pour former le produit vectoriel c = a  b se note :
c i =  i jk a j b k
En effet, on peut vérifier que l'on a bien :
c1 = 123 a 2 b 3 + 132 a 3b 2 = a 2 b 3 − a 3b 2
26
Cours MDF Aroua Aouadi Enit 2022-2023
 Éléments de calcul tensoriel
Produit tensoriel
Le produit tensoriel d’un tenseur d’ordre m et un tenseur d’ordre n est un tenseur d’ordre m+n. Par
exemple :
• Le produit tensoriel de deux vecteurs est un tenseur d’ordre 2 :
BC = A Bi C j = A i j

• Le produit tensoriel de deux tenseurs d’ordre 2 est un tenseur d’ordre 4 :

BC = A Bi jC kl = A i jkl
Quelques propriétés :

( B  C ).X = B ( C.X )
( A  B).( C  D ) = ( B.C )( A  D )
27
Cours MDF Aroua Aouadi Enit 2022-2023
 Éléments de calcul tensoriel
Contraction d’un tenseur
La contraction d’un tenseur d’ordre m et d’un tenseur d’ordre n, donne un nouveau tenseur d’ordre
m + n − 2.

B.C = A Bi j C kl = Bi j C jl = A il

Ainsi, le produit scalaire : a i b i (de rang 0) est la contraction du produit tensoriel : a i b j (de rang 2).

La contraction d’un tenseur de rang 3, a i jk , résulte un vecteur a i ji de rang 1 (n – 2 = 3 – 2 = 1), on écrit :

a i ji =  ik a i jk

28
Cours MDF Aroua Aouadi Enit 2022-2023
 Éléments de calcul tensoriel
Double contraction
Un produit doublement contracté est de rang n-4.
La double contraction d’un tenseur d’ordre m et d’un tenseur d’ordre n, donne un nouveau tenseur
d’ordre m + n − 4.
B C = A Bi j C kl = Bi j C ji = A ii = A

( ) ( ) ( )(
B  C = Bi j e i  e j  C kl e k  e l = Bi j C kl e i  e j  e k  e l )
= Bi j C kl ( e j .e k ) ( e i . e l )
 jk  il

= Bi j C ji
29
Cours MDF Aroua Aouadi Enit 2022-2023
 Éléments de calcul tensoriel
Double contraction

Autrement ( ) ( )
B  C = Bi j e i  e j  C kl e k  e l = Bi j C kl e i  e j  e k  e l( )( )
= Bi j C kl ( e i .e k )( e j. e l )
 ik  jl

= Bi j C i j

3 3 3
B  C = Bi j C i j =  B
i =1 j =1
ij Ci j =  (B
i =1
i1 C i1 + Bi2 C i2 + Bi3 C i3 )

= ( B11 C11 + B12 C12 + B13 C13 )

+ ( B 21 C 21 + B 22 C 22 + B 23 C 23 )
+ ( B31 C 31 + B32 C 32 + B33 C 33 ) 30
Cours MDF Aroua Aouadi Enit 2022-2023
 Éléments de calcul tensoriel
trace d’un tenseur
La trace d’un tenseur d’ordre 2 est la somme de ses termes diagonaux

Tr A = A i j  i j = A ii
Avec  i j = d = 3

31
Cours MDF Aroua Aouadi Enit 2022-2023
 II. CINÉMATIQUE DES
FLUIDES

32

Cours MDF Aroua Aouadi Enit 2022-2023

 Introduction

La cinématique des fluides repose sur la notion d’espace et de temps. En effet, elle décrit l’évolution d’un
milieu fluide dans l’espace-temps sans tenir compte ni de la nature des forces ni des lois et des causes
résultant le mouvement. L’espace et le temps sont indépendants. L’écoulement est décrit principalement
par des vitesses en différents points de l’espaces et à des instants de temps différents.

L’évolution cinématique des fluides est décrite soit par la méthode Lagrangienne soit la méthode
Eulérienne.

La description Lagrangienne s’avère peu adaptée pour décrire le mouvement d’un milieu fluide. C’est
plutôt la méthode eulérienne qui est la mieux appliquée pour suivre la cinétique des fluides.

33
Cours MDF Aroua Aouadi Enit 2022-2023
 Qu’est ce qu’un fluide ?

Liquide/gaz
• Pas de forme propre
• S’écoule si on lui applique une force
• Prend la forme du récipient
• Mobile : Aucune résistance aux déformations : ensemble
de molécules très mobiles les unes par rapport aux autres et
se déplacent sous le moindre effort

Solide
• Possède une forme propre
• Rigide : Résistance aux déformations aux cisaillements : la
déformation croît linéairement avec la contrainte appliquée
jusqu’à rupture 34
Cours MDF Aroua Aouadi Enit 2022-2023
 Qu’est ce qu’un fluide ?

• Molécules rapprochées et ordonnées • Espace entre les molécules

• Forces de cohésion intermoléculaires très • Molécules sont plus espacées encore plus grand
grandes qui permettent au solide de maintenir • Forces intermoléculaires plus faibles • Forces intermoléculaires de
sa forme initial et de résister aux déformations • Solide désordonné ou un gaz dense cohésion extrêmement petites

35
Cours MDF Aroua Aouadi Enit 2022-2023
 Le fluide est un milieu continu

36
Cours MDF Aroua Aouadi Enit 2022-2023
 II.1. Description Eulérienne - Description Lagrangienne
II.1.1. Description Lagrangienne
Suivre la position et la vitesse de chaque particule matérielle en fonction de temps (suivre la trajectoire
des particules à chaque instant à partir de sa position d’origine).

Cette approche est utilisée en mécanique des solides où il y a intérêt de suivre le mouvement de point
matériels (mécanique de point). Cependant, de point de vue expérimental, il est difficile de l’appliquer en
mécanique des fluides

37
Cours MDF Aroua Aouadi Enit 2022-2023
 II.1. Description Eulérienne - Description Lagrangienne
II.1.1. Description Lagrangienne
II. 1.1.1. Trajectoire de particule
La trajectoire est le chemin parcouru au cours du temps par une particule identifiée à un instant initial en
une position initiale. La trajectoire est relative à une même particule et à des instants différents.

 X (t) 
Une particule est représentée par ses coordonnées cartésiennes dans un trièdre : OP  Y (t) 
 
 Z (t) 
 
Avec OP = R (t) est le rayon de la particule à l’instant t.

38
Cours MDF Aroua Aouadi Enit 2022-2023
 II.1. Description Eulérienne - Description Lagrangienne
II.1.1. Description Lagrangienne
II. 1. 1. 1. Trajectoire de particule
Le champ lagrangien des vitesses est l’ensemble des vitesses des particules situées à la position initiale R 0
à l’instant t0 (t= 0).

La vitesse d’une particule est tangente à sa trajectoire. Par définition, c’est la dérivé du vecteur position
de cette particule :

V (P) =
dOP dR
dt
=
dt
= V R0 , t ( )

L’accélération d’une particule s’écrit donc :

dV d 2 OP
A = ==
dt dt 2

39
Cours MDF Aroua Aouadi Enit 2022-2023
 II.1. Description Eulérienne - Description Lagrangienne

II. 1. 2. Description Eulérienne

Contrairement à l’approche Lagrangienne, au lieu de suivre la trajectoire de la particule, l’observateur

fixe un point et il s’intéresse à connaître à chaque instant la vitesse des différentes particules qui au cours
du temps vont passer par ce point. On s’intéresse donc à un point fixe qui voit défiler le fluide.

Les coordonnées d’espace et le temps sont décorrélés.

Cette approche est utilisée en mécanique des fluides et en mécanique des milieux continues à grandes
déformations.

40
Cours MDF Aroua Aouadi Enit 2022-2023
 II.1. Description Eulérienne - Description Lagrangienne
II. 1. 2. Description Eulérienne
A chaque instant t, l’écoulement est décrit par un champ de vitesse.
La vitesse V ( M , t ) représente la vitesse de la particule P passant par le point M (x, y, z) à l’instant t.
Cette vitesse change dans le temps.

La particule dans le point M1 est différente de la particule dans le point M2.

41
Cours MDF Aroua Aouadi Enit 2022-2023
 II.2. Relation entre les méthodes Eulérienne et Lagrangienne
II. 2. 1. Dérivée matérielle
De point de vue mathématique, les deux méthodes sont équivalentes. À partir d’une méthode, on peut obtenir l’autre.
La méthode Lagrangienne : La méthode Eulérienne :
définit le vecteur position : dx dy dz définit le champ de vitesse :
Par dérivation u = , v= , w =
x = F1 ( x 0 , y 0 , z 0 , t ) dt dt dt u = f 1(x , y , z , t)
y = F2 ( x 0 , y0 , z 0 , t ) v = f 2 (x , y , z , t)
z = F3 ( x 0 , y 0 , z 0 , t ) x = u dt w = f 3 (x , y , z , t)
y = v dt Par intégration
z = w dt
Dérivées

Temps =?
t
Espace d
=?
dt
Les deux D
=?
Dt
42
Cours MDF Aroua Aouadi Enit 2022-2023
 II.2. Relation entre les méthodes Eulérienne et Lagrangienne
II. 2. 2. Dérivée particulaire - Dérivée matérielle
Considérons la fonction F ( x , y , z , t )

La différentielle totale de F s’écrit

F F F F
dF = dt + dx + dy + dz
t x y z
F F F F
dF = dt + ( i x + iy + i z ) . (dx i x + dy i y + dz i z )
t x y z

F dR = V dt
F
En indicielle v i dt
xi
F
dF = ( + V .  F) d t Pour différencier, on note la dérivée matérielle : D
t Dt D 
= ( + V . )
Dt t
dF F
= ( + V .  F)
dt t 43
Cours MDF Aroua Aouadi Enit 2022-2023
 II.2. Relation entre les méthodes Eulérienne et Lagrangienne
II. 2. 3. Récapitulations

Dérivée particulaire = Dérivée totale = Dérivée matérielle

D d 
= = ( + V . )
Dt dt t

D d  
Forme indicielle = = ( + vi )
Dt dt t xi

44
Cours MDF Aroua Aouadi Enit 2022-2023
 II.2. Relation entre les méthodes Eulérienne et Lagrangienne
II. 2. 4. Dérivée particulaire et accélération
Ce champ d’accélération est donné en description eulérienne par :
dV V
a= = + (V .  ) V
dt t

Forme indicielle d vi vi vi

ai = =( + vj )
dt t x j
Accélération locale Accélération spatiale convective

u u u u
ax = + u + v + w
t x y z
v v v v
ay = + u + v + w
t x y z
w w w w
az = + u + v + w
t x y z
45
Cours MDF Aroua Aouadi Enit 2022-2023
 II.2. Relation entre les méthodes Eulérienne et Lagrangienne
II. 2. 5. Comparaison entre les descriptions Eulérienne et Lagrangienne

46
Cours MDF Aroua Aouadi Enit 2022-2023
 II.3. Étude graphique des écoulements
II. 3. 1. Ligne de courant
Une ligne de courant à un instant t donné correspond à la courbe ou la trajectoire à laquelle le
vecteur vitesse est tangent en chaque point.

Comme le vecteur vitesse, la ligne de courant varie dans le temps et dans l’espace (ligne à t = t0  à t = t1).

Puisque en tout point il n’y a qu’un seul vecteur vitesse avec une direction donnée, il ne peut y avoir
d’intersections de lignes de courant.
47
Cours MDF Aroua Aouadi Enit 2022-2023
 II.3. Étude graphique des écoulements
II. 3. 1. Ligne de courant
En considérant un déplacement élémentaire dOM , les lignes de courant sont déterminées à partir
de la condition suivante :
V( M, t )  dOM = 0

De point de vue mathématique, pour déterminer les lignes de courant, il suffit de résoudre à un instant
t fixe le système d’équations différentielles suivant :

dx dy dz
= =
u ( x , y, z, t ) v ( x , y, z, t ) w ( x , y, z, t )

On obtient les équations de la ligne de courant en intégrant par rapport à x, y, z (variables), t étant
considéré comme une constante.

48
Cours MDF Aroua Aouadi Enit 2022-2023
 II.3. Étude graphique des écoulements
II. 3. 2. Trajectoire
Une trajectoire est un cheminement décrit au cours du temps par une particule quelconque. Une
trajectoire considère une même particule et à des instants différents.

 dx
 dt = U ( x , y, z, t )

Les équations de la trajectoire s’écrivent :  dy = V ( x , y, z, t )
 dt
 dz
 = W ( x , y, z, t )
 dt 49
Cours MDF Aroua Aouadi Enit 2022-2023
 II.3. Étude graphique des écoulements
II. 3. 3. Comparaison : Ligne de courant – Trajectoire
La trajectoire considère une même particule à des instants différents alors que les lignes de
courant sont relatives à des particules différentes au même instant.

Trajectoire de la particule P
Ligne de courant à t0

50
Cours MDF Aroua Aouadi Enit 2022-2023
 II.3. Étude graphique des écoulements
II. 3. 3. Comparaison : Ligne de courant – Trajectoire
En écoulement permanent, le champ de vitesse ne varie pas dans le temps. Les lignes de courant et
trajectoires passant par le même point se confondent parce que le champ de vitesse ne dépend
que de l’espace.

51
Cours MDF Aroua Aouadi Enit 2022-2023
 II.3. Étude graphique des écoulements
II. 3. 4. Tube de courant
Un tube de courant est l’ensemble des lignes de courant s’appuyant au même instant sur un contour
fermé quelconque. Un tube de courant est défini de telle sorte qu’il n’y ait pas d’intersection entre les
lignes de courant.
Le vecteur vitesse est tangent en tout point au tube de courant. Il en résulte que le fluide ne
traverse pas la paroi d’un tube de courant, et donc le débit qui traverse une section droite d’un tube
de courant reste constant.

52
Cours MDF Aroua Aouadi Enit 2022-2023
 II.3. Étude graphique des écoulements
II. 3. 5. Ligne d’émission
Une ligne d’émission relative à un point E représente l’ensemble des points ou de positions occupées par des
différentes particules de fluides ayant passé à un instant antérieur t par le même point E (point d’émission).
Ces lignes d’émission sont obtenues expérimentalement par injection continue d’un traceur (colorant)

Dans un écoulement stationnaire, trajectoires, lignes de courant et lignes d’émission sont confondues.
53
Cours MDF Aroua Aouadi Enit 2022-2023
 II.4. Déformation des milieux continus
Dans son mouvement, la particule fluide peut changer de position,
modifier sa forme et son orientation. Le mouvement de la particule
fluide peut alors être décrit comme une combinaison de translation,
rotation cisaillement et de compressions d’éléments rigides.

Ci-dessous, un exemple de déformation totale

54
Cours MDF Aroua Aouadi Enit 2022-2023
 II.4. Déformation des milieux continus
II. 4. 1. Tenseur de taux de déformation

Soit M1 et M2 deux points voisins d’une même particule fluide :

M1(x, y, z) et M2(x + dx, y + dy, z + dz)

La distance entre M1 et M2 est : dr = M1M2. Cette distance n’est pas conservée lors de son évolution
entre les deux points à cause de la déformation.

On note les vecteurs vitesse respectifs aux deux points tels que : V M1 (u, v, w) et V M2 (u ', v ', w ')

V M2 = V( r + dr ) = V M1 + dV = V( r ) + dV

55
Cours MDF Aroua Aouadi Enit 2022-2023
 II.4. Déformation des milieux continus
II. 4. 1. Tenseur de taux de déformation
Par définition le développement au premier ordre des composantes
de vecteur vitesse (puisque élément infinitésimal) :

 u u u V M2 = V( r + dr ) = V M1 + dV = V( r ) + dV
 u ' = u + dx + dy + dz
x y z

 v v v
 v ' = v + dx + dy + dz
 x y z
G est le tenseur de taux de déformation
 w w w
 w ' = w + dx + dy + dz  u u u 
 x y z  x y z 
 
 v v v 
V( r + dr ) = V( r ) + G .dr V =G = 
Avec x y z 
 
V( r + dr ) = V( r ) +  V .dr  w w w 
 x y  z 

56
Cours MDF Aroua Aouadi Enit 2022-2023
 II.4. Déformation des milieux continus
II. 4. 2. Termes diagonaux du tenseur de taux de déformation
Si on considère que les termes non-diagonaux sont nuls, le tenseur de taux de déformation G s’écrit :

 u   u
 x 0 0 
 u ' = u + dx
  x On s’intéresse à un écoulement bidimensionnel

 v   v d’une particule fluide dans un plan (x, y) afin
G = 0 0  v ' = v + dy d’examiner la signification physique des
x y
   termes diagonaux
 w   w
0 0
 z  w ' = w + dz
  z

57
Cours MDF Aroua Aouadi Enit 2022-2023
 II.4. Déformation des milieux continus
II. 4. 2. Termes diagonaux du tenseur de taux de déformation

58
Cours MDF Aroua Aouadi Enit 2022-2023
 II.4. Déformation des milieux continus
II. 4. 2. Termes diagonaux du tenseur de taux de déformation

u /  x et  v /  y ont été vérifié, respectivement, comme le taux d’élongation par unité de longueur
dans la direction de x et la direction de y.

Ainsi l’accroissement de la longueur du rectangle ABCD est calculé comme :

d ( AB ) A 'B' − AB dx + u / x dx dt u
= = = dt
AB AB dx dx 59
Cours MDF Aroua Aouadi Enit 2022-2023
 II.4. Déformation des milieux continus
II. 4. 2. Termes diagonaux du tenseur de taux de déformation

60
Cours MDF Aroua Aouadi Enit 2022-2023
 II.4. Déformation des milieux continus
II. 4. 2. Termes diagonaux du tenseur de taux de déformation

61
Cours MDF Aroua Aouadi Enit 2022-2023
 II.4. Déformation des milieux continus
II. 4. 2. Termes non-diagonaux du tenseur de taux de déformation
On considère maintenant que les termes diagonaux sont nuls, le tenseur de taux de déformation G s’écrit :

 u u   u u
0 y z  u ' = u +  y dy +  z dz On s’intéresse à un écoulement bidimensionnel
  
 v v   v v d’une particule fluide dans un plan (x, y) afin
G =  0 v ' = v + dx +
z 
dz d’examiner la signification physique des
x x z
   termes non-diagonaux tel que :
 w w   w w u v
 x y 0   w ' = w + dx + dy et sont non nuls
    x  y y x

62
Cours MDF Aroua Aouadi Enit 2022-2023
 II.4. Déformation des milieux continus
II. 4. 2. Termes non-diagonaux du tenseur de taux de déformation

63
Cours MDF Aroua Aouadi Enit 2022-2023
 II.4. Déformation des milieux continus
II. 4. 2. Termes non-diagonaux du tenseur de taux de déformation

64
Cours MDF Aroua Aouadi Enit 2022-2023
 II.4. Déformation des milieux continus
II. 4. 2. Termes non-diagonaux du tenseur de taux de déformation

65
Cours MDF Aroua Aouadi Enit 2022-2023
 II.4. Déformation des milieux continus
II. 4. 2. Termes non-diagonaux du tenseur de taux de déformation

66
Cours MDF Aroua Aouadi Enit 2022-2023
 II.4. Déformation des milieux continus
II. 4. 3. Récapitulatif sur le tenseur du taux de déformation G :
• Termes diagonaux : Taux d’élongation (causant une élongation ou une compression) : u ,  v ,  w
dx dy dz
• Termes non-diagonaux :
Déformation angulaire pure (symétrique) :  =  u ,  =  v
dt dy dt dx
   v u 1 v u 1 v u 1
Rotation pure (anti-symétrique) : =− = = − = ( − ) = ( − ) =  V  z
dt dt dx dy 2 dx dy 2 dx dy 2

Comme tout tenseur au second ordre, le tenseur de taux de déformation G se décompose en parties
symétrique (e) et antisymétrique () : G = e +  tel que :

V( r + dr ) = V( r ) + dV = V( r ) + G .dr

dV = G .dr = e.dr + ω .dr

67
Cours MDF Aroua Aouadi Enit 2022-2023
 II.4. Déformation des milieux continus
II. 4. 3. Récapitulatif sur le tenseur du taux de déformation G :
e : Partie symétrique représentant le taux de
déformation pure (élongation et déformation angulaire)
G=e+ avec
 : Partie antisymétrique représentant le tenseur de
taux de rotation pure

 u u u   u 1  u v 
+
1  u w   
+
1  u v 
−
1  u w  
−
  0 
 x y z  x 
2   y  x   
2  z x    
2   y  x  
2   z  x  
  
 v v v   1  u v  v 1  v w   1  u v  1  v w  
G =  =   + + +  − + − 
2   y  x   
2   z  y    2   y  x   
 0
x y z  y 2  z  y  
 
 w w w   1  u w  1  v w  w
  
1 u w  1  v w 

 x  + +   − + − + 
 z   2   z  x       
  0 
 y 2  z  y  z   2  z x  2  z  y  

symétrique e antisymétrique 
68
Cours MDF Aroua Aouadi Enit 2022-2023
 II.4. Déformation des milieux continus
II. 4. 3. Récapitulatif sur le tenseur du taux de déformation G :
 1  u v  1  u w  
0 − −

 2   z  x  
2   y  x  
 0 − z y 
 1  u v  1  v w    
ω =  − +  0  −   =  z 0 − x 
 2  y x  2  z  y  
− y x 0 
 1  u w   
1  v w 
  − +   − +  0 
 
2  z  x  2  z  y  

69
Cours MDF Aroua Aouadi Enit 2022-2023
 II.4. Déformation des milieux continus
II. 4. 3. Récapitulatif sur le tenseur du taux de déformation G :
Finalement :
G=e+
Symétrique Antisymétrique

u i
Gi j = = e i j + i j
x j

1 u i u j 1 u i u j
ei j = ( + ) i j = ( − )
2 x j x i 2 x j x i

1 1
e = ( u + ( u) T ) ω = ( u − ( u) T )
2 2

70
Cours MDF Aroua Aouadi Enit 2022-2023
 II.5. Fonction courant
II. 5. 1. Fonction courant d’un écoulement plan incompressible
La vorticité  comme la vitesse V est un champ vectoriel. La ligne de courant de la vorticité est en tout point
tangente à la direction du vecteur vorticité. Elle est appelée ligne de vortex.

On a montré que la divergence du rotationnel d’un vecteur est nulle :  . ( x A) = 0

Pour un écoulement incompressible :  .V = 0

D’un point de vue mathématique, on peut définir un champ vectoriel A pour un écoulement permanent
incompressible comme : V =  x A

 u = A z /  y − A y / z

V =  x A =  v = A x / z − A z /  x
 w = A /  x − A /  y
 y x

71
Cours MDF Aroua Aouadi Enit 2022-2023
 II.5. Fonction courant
II. 5. 1. Fonction courant d’un écoulement plan incompressible
Considérons un écoulement 2-D dans le plan (x, y) où u et v sont fonction de x, y et de t (mais pas z) et w=0.
L’écoulement est invariant suivant z.
 / z = 0

 u = A z /  y − A y / z  u = A z /  y
 
V =  x A =  v = A x / z − A z /  x V =  v = − A z /  x
 w = A /  x − A /  y w = 0
 y x 
Cette équation montre que dans le plan (x,y) la vitesse est définie en tout point à l’aide d’une fonction scalaire Az (x, y)
dépendant de la position, Az (x, y).
On définit ainsi La Fonction Courant de l’écoulement incompressible par :  ( x, y ) = A z ( x, y )

 
 u =
La fonction courant  permet de déduire le champ de vitesse bidimensionnel tel que: y

 v = − 
 x
72
Cours MDF Aroua Aouadi Enit 2022-2023
 II.5. Fonction courant
II. 5. 1. Fonction courant d’un écoulement plan incompressible
II. 5. 1. 1. Ligne de courant
La fonction courant  (x, y) définit les lignes de courant de l’écoulement. En écoulement plan, le long d’une ligne
de courant, on a la relation différentielle :

 = constante
 
d = dx + dy = 0
x y

Ligne de courant  
= −v =u
x y

− v dx + u dy = 0

dy v
=
dx u

73
Cours MDF Aroua Aouadi Enit 2022-2023
 II.5. Fonction courant
II. 5. 1. Fonction courant d’un écoulement plan incompressible
II. 5. 1. 1. Ligne de courant
 = constante   dy v
d = dx + dy = 0 − v dx + u dy = 0 =
x y dx u
Dans ce cas, la pente dy/dx de la courbe  = constante est la même que celle de la vitesse v/u. Ceci permet de
déduire que  = constante est une ligne de courant.
Si on connait la fonction courant, on peut tracer les lignes de courant de l’écoulement.

74
Cours MDF Aroua Aouadi Enit 2022-2023
 II.5. Fonction courant
II. 5. 1. Fonction courant d’un écoulement plan incompressible
II. 5. 1. 2. Débit volumique
La différence entre les valeurs de  sur deux lignes de courant différentes est égale au débit de fluide entre les deux
lignes de courant caractérisées par des fonctions courant constantes. Le débit volumique entre deux point M et M’ est :
dq v = u dy − v dx

 
dq v = dx + dy = d
x y

Entre deux points 1 et 2 reliant les deux points M et M’ sur deux lignes de
courant différentes de constantes  1 et  2 . Le débit volumique s’écrit :

2 2
qv = qv =  d =  2 − 1
1 1

75
Cours MDF Aroua Aouadi Enit 2022-2023
 II.5. Fonction courant
II. 5. 1. Fonction courant d’un écoulement plan incompressible
II. 5. 1. 2. Débit volumique
La différence entre les valeurs de  sur deux lignes de courant différentes est égale au débit de fluide entre les deux
lignes de courant caractérisées par des fonctions courant constantes. Le débit volumique entre deux point M et M’ est :

dq v = u dy − v dx

 
dq v = dx + dy = d
x y

Entre deux points 1 et 2 reliant les deux points M et M’ sur deux lignes de courant différentes dans le plan (x,y). Le
débit volumique s’écrit :

2 2
qv = qv =  d =  2 − 1
1 1

76
Cours MDF Aroua Aouadi Enit 2022-2023
 II.6. Théorème de transport de Reynolds

Le théorème de transport de Reynolds a pour but d’établir une relation ou un lien entre la méthode Eulérienne
basée sur un volume de contrôle versus la méthode Lagrangienne qui s’appuie sur un volume matériel.
Autrement il faut trouver une loi ou une formulation mathématique qui nous permet d’exprimer de point de vue
Eulérien l’évolution d’un volume au cours du temps encore plus si ce volume bouge avec le mouvement de fluide, il
est une fonction du temps.
Pour établir l’expression du taux de variation dans le temps (entre t et t + dt) de la propriété f sur le volume du fluide, qui
bouge avec le mouvement du fluide on considère un volume fixe (volume de contrôle), délimité par une surface (surface de
contrôle) à travers laquelle on pourra comptabiliser le flux de f :

dF f
dt
=  t dV +  f v . n dS
Vc Sc

Taux de variation dans le Flux à travers la surface du

temps sur le volume fixe Vc volume de contrôle Sc
77
Cours MDF Aroua Aouadi Enit 2022-2023
III. LOIS DE CONSERVATION

78

Cours MDF Aroua Aouadi Enit 2022-2023

 III.1. Conservation de la masse

En absence de réactions chimiques, la masse de la matière se conserve. Ceci signifie que dans un volume V, le
fluide peut s’accumuler ou il peut y avoir un échange avec l’extérieur mais pas en créer ni en détruire.

m =
 V
 dV

dS Le déplacement ente t et t + dt : v.n dt

Le volume d’un fluide s’écoulant à travers une surface dS sur une durée dt :
S V
v.n dS dt

Le débit volumique : volume de fluide s’écoulant par unité de temps à travers l’élément de surface : q v =  S
v.n dS

Le débit massique : masse de fluide s’écoulant par unité de temps à travers l’élément de surface : q m = 
S
 v.n dS

79
Cours MDF Aroua Aouadi Enit 2022-2023
 III.1. Conservation de la masse
III. 1. 1. Principe de conservation
III. 1. 1. 1. Bilan d’une grandeur en général
Bilan d’une grandeur G dans un volume V

Flux entrant
Flux sortant
de G
de G
V

Destruction
Production de G
de G

dG
= Flux entrant - Flux sortant + Production - Destruction
dt

Bilan temporel Bilan spatial

80
Cours MDF Aroua Aouadi Enit 2022-2023
 III.1. Conservation de la masse
III. 1. 2. Conservation de la masse
III. 1. 2. 1. Forme intégrale
En absence de réactions chimiques et de flux thermique, la masse se conserve. dm Flux entrant
= - Flux sortant
Ce que signifie qu’il y a ni production ni la destruction : dt

Bilan temporel Bilan spatial


Entre deux instant t et t +dt infiniment voisins (bilan temporel) :  ( x, y, z; t + dt ) −  ( x, y, z; t ) = (développement au
1er ordre de Taylor).
t

 
d   dm   dm 
 =  dV = dV =   −  dt  1
dt V(t) V(t)  t  dt  entrant  sortant

D’autre part :  dm 

 dm 
−   =  v .n dS 2
 dt sortant  dt entrant S

 

1 et 2 dV +  v .n dS = 0
V(t)  t S 81
Cours MDF Aroua Aouadi Enit 2022-2023
 III.1. Conservation de la masse
III. 1. 2. Conservation de la masse
III. 1. 2. 1. Forme intégrale
dF f
En utilisant le théorème de transport de Reynolds :
dt
=  t dV +  f v . n dS
V
c Sc

On identifie : f ( x, y, z, t ) =  ( x, y, z, t )

Et m =
 V
 ( x, y, z, t ) dV

 
dm 
= dV +  v .n dS = 0
dt Vc  t Sc

82
Cours MDF Aroua Aouadi Enit 2022-2023
 III.1. Conservation de la masse
III. 1. 2. Conservation de la masse
III. 1. 2. 2. Forme locale

On transforme l’intégrale de surface en intégrale de volume en utilisant le théorème de divergence d’Ostrogradski :

 Sc
 v .n dS =
Vc
 . ( v) dV


dm    dV → 0 
=   t +  . ( v)  dV = 0 +  . ( v) = 0
dt Vc   t
Équation de continuité

 d
 + v . ( ) +  ( . v ) = 0 +  ( . v ) = 0
t dt

d
dt
83
Cours MDF Aroua Aouadi Enit 2022-2023
 III.1. Conservation de la masse
III. 1. 2. Conservation de la masse
III. 1. 2. 3. Formulation de conservation de la masse sous des conditions d’écoulements (forme locale)

Forme différentielle de la conservation de la masse :
t
+ . ( v ) = i  q Vi

Écoulement permanent : = 0  . ( v ) =   q Vi
t i


Écoulement incompressible :  = cste = 0  . ( v ) = v . ( ) +  ( . v )   . ( v ) =  ( . v )
t
.v = i q Vi
Écoulement conservatif : i  q Vi =0

Écoulement conservatif permanent :  . ( v ) = 0

Écoulement conservatif incompressible :  . v = 0

84
Cours MDF Aroua Aouadi Enit 2022-2023
 III.1. Conservation de la masse
III. 1. 2. Conservation de la masse
III. 1. 2. 4. Taux de dilatation volumique

  D
+  . ( v) = + v .   +  . v = +  .v
t t Dt

1 D
= − .v
 Dt

1 DV
= .v
V Dt

85
Cours MDF Aroua Aouadi Enit 2022-2023
 III.1. Conservation de la masse
III. 1. 2. Conservation de la masse
III. 1. 2. 5. Forme conservative
L'équation de conservation de la masse permet une formulation dite "conservative" de la dérivée matérielle :

D   
 =  + ( v .  )   =    + . (  v) −     +  . ( v ) 
Dt  t  t t 
En utilisant l’équation de conservation de la masse :

D  
 = +  . (  v)
Dt t

Forme indicielle des composantes cartésiennes d’un champs au second ordre :

D i j   i j     i j 
 =  + vk i j  = + (  i j v k )
Dt  t x k  t x k

86
Cours MDF Aroua Aouadi Enit 2022-2023
 III.2. Conservation de la quantité de mouvement
III. 2. 1. Forces exercées sur un fluide
Soit un volume V fermé par une surface S. On considère des petits éléments de surface dS et des petits éléments
volume dV

Dans le volume de fluide V, on distingue deux types de forces dans les milieux continus :
• Les forces volumiques qui s’exercent sur chaque particule fluide interne au volume V ou sur chaque élément de
volume dV (poids, force d’inertie, magnétique…)

• Les forces surfaciques ou de contact qui s’exercent sur la frontière S du volume V ou sur chaque élément de
surface dS (pression, frottement visqueux).

87
Cours MDF Aroua Aouadi Enit 2022-2023
 III.2. Conservation de la quantité de mouvement
III. 2. 1. Forces exercées sur un fluide
III. 2. 1. 1. Forces volumiques
Les forces volumiques sont exprimées à partir de champs de densités massiques : d F =  f dV

Poids :
Le poids est une force volumique. Dans ce cas, la densité de force f représente l’accélération de la pesanteur (f = g).

dm g =  dV g
dV
V

 dV g
p=
òòòV g dV

On peut également citer les forces externes et les forces mutuelles résultant des processus d’interactions identifiées de la
physique, gravitation, champs électrique, électromagnétique. Ces forces sont toujours formulées en une densité massique.

88
Cours MDF Aroua Aouadi Enit 2022-2023
 III.2. Conservation de la quantité de mouvement
III. 2. 1. Forces exercées sur un fluide
III. 2. 1. 2. Forces surfaciques ou de contact
Les forces surfaciques sont représentées via une densité superficielle Tn : d F = t n dS

Les forces de contact ne dépendent que du point et de l'orientation de l'élément de surface considéré en ce point.
Exemples de forces surfaciques : pression, contrainte de cisaillement. Les dimensions de la force surfacique sont celles
d’une pression (Pa = N. m-²).

Pression :
S

n
dS
M dFp= -pn dS
dS
Fp = òò dFp Fp = òò
S
-pn dS
S
M n
V

89
Cours MDF Aroua Aouadi Enit 2022-2023
 III.2. Conservation de la quantité de mouvement
III. 2. 2. Equations intégrales de quantité de mouvement
La quantité de mouvement contenue dans un volume de contrôle D est exprimée comme :
M =   v dV
D

En appliquant le principe fondamental de la dynamique, la dérivé matérielle de la quantité de mouvement est égale
à la somme des forces extérieures. Ce principe implique que toute variation de la quantité de mouvement M est due
à l’action d’une force appliquée par le milieu extérieur et réciproquement.

=  Fext
DM
Dt
Les forces extérieures représentent la somme des contributions des forces volumiques représentées par une
intégrale de volume et des forces surfaciques (ou de contact) représentées par une intégrale de surface
correspondant au flux du tenseur des contraintes à travers un élément de surface.

DM D
Dt
=
Dt   v dV =  t n dS +   f dV
D S D

90
Cours MDF Aroua Aouadi Enit 2022-2023
 III.2. Conservation de la quantité de mouvement
III. 2. 2. Equations intégrales de quantité de mouvement
DM D
Dt
=
Dt   v dV =  t n dS +   f dV
D S D

L’application du théorème de transport de Reynolds à la dérivée totale de la quantité de mouvement du fluide, permet
d’exprimer la conservation de quantité de mouvement comme :

D  v
Dt   v dV =   t dV +   v ( v .n) dS =  t n dS +   f dV
D D S S D

 v
   t dV = −   v ( v .n) dS +  t n dS +   f dV
D S S D

Taux d’accumulation de Flux de QDM rentrant Forces extérieures

QDM dans le volume D à travers la surface S s’exerçant sur D

91
Cours MDF Aroua Aouadi Enit 2022-2023
 III.2. Conservation de la quantité de mouvement
III. 2. 2. Equations intégrales de quantité de mouvement
III. 2. 2. 1. Forces surfaciques : Tenseur des contraintes
 v  v
  t dV = −   v ( v .n) dS +  t n dS +   f dV    t dV +   .  vv dV −   f dV =  t n dS
D S S D D D D S

v 
 ( t
+v
t
+ v  .  v +  ( v.  ) v −  f ) dV =  t n dS
D S

v 
 ( ( t
+ ( v.  ) v) + v (
t
+  .  v) −  f ) dV =  t n dS
D S

v
En considérant l’équation de conservation de la masse :  ( ( t
+ ( v.  ) v) −  f ) dV =  t n dS
D S

  ( a −  f ) dV =  t n dS =  T.n dS
D S S

92
Cours MDF Aroua Aouadi Enit 2022-2023
 III.2. Conservation de la quantité de mouvement
III. 2. 2. Equations intégrales de quantité de mouvement
III. 2. 2. 2. Equations aux dérivées partielles exprimant la conservation de quantité de mouvement
 v
 (
t
+  .  vv −  f −  . T)dV = 0
D

Dv v  v
a =  = +  ( v . ) v = +  . ( vv ) =  . T +  f
Dt t t

Dvi  vi  vi  vi  ( vi v j )  i j p  ij
 ai =  = +  vj = + = +  fi = − + + fi
Dt t x j t x j x j  xi x j

T = −p + 

i j = − p i j + i j

93
Cours MDF Aroua Aouadi Enit 2022-2023
 III.3. Conservation du moment cinétique
III. 3. 1. Equation intégrale du moment cinétique
Le moment cinétique stocké dans un volume de contrôle D par rapport
à un point O s’exprime comme :

M (O, D) =  OM  ( v) dV =  x  ( v) dV

D D

Le principe de conservation du moment cinétique stipule que toute variation du moment cinétique par rapport au point
O est due à l’action du moment d’une force extérieure exercée sur D :

D M (O, D)
Dt
=  Moment O (Fext )

94
Cours MDF Aroua Aouadi Enit 2022-2023
 III.3. Conservation du moment cinétique
III. 3. 1. Equation intégrale du moment cinétique
En utilisant le pseudo-tenseur d’orientation ɛ , l’équation de conservation du moment cinétique devient :
Rappel : (x  v) i = i jk x j v k

D M (O, D) D
Dt
=  D Dt
(x   v ) dV =  x  tn dS +  x   f dV
S D

D
 D Dt ( i jk x j vk ) dV =  i jk x j t n k dS +   i jk x j f k dV
S D

D M (O, D) D
Dt
= 
D Dt
(x   v ) dV =  x  T.n dS +  x   f dV
S D

D

D Dt ( i jk x j vk ) dV =  i jk x j km n m dS +   i jk x j f k dV
S D

95
Cours MDF Aroua Aouadi Enit 2022-2023
 III.3. Conservation du moment cinétique
III. 3. 1. Equation intégrale du moment cinétique
En appliquant le théorème de la divergence :
D
D Dt ( i jk x j vk ) dV =  i jk x j km n m dS +   i jk x j f k dV
S D

D 
 D Dt
( i jk x j v k ) dV =  i jk
x m
(x j km ) dV +   i jk x j f k dV
D D

On développe le premier membre à droite de l’équation :

 km x j k m

 i jk
x m
(x j km ) dV = 
i jk x j
 xm
dV + i jk  k m
 xm
dV =   i jk x j xm
dV +  i m k  k m dV
D D D D D

 jm

D Dv k Dx j Dv k
D’autre part, on note que : (i jk x j v k ) = i jk x j + i jk v k = i jk x j
Dt Dt Dt Dt
96
Cours MDF Aroua Aouadi Enit 2022-2023
 III.3. Conservation du moment cinétique
III. 3. 1. Equation intégrale du moment cinétique
Dv k k m
  ( i jk xj
Dt
−  i jk x j
xm
− i m k  k m −  i jk x j f k ) dV = 0
D

 Dv k k m 
 i jk x j ( Dt −  x −  f k ) dV − i jk  k j dV  = 0
 m 
D

En considérant la loi fondamentale de la dynamique ainsi l’équation de conservation de quantité de mouvement ,

l’équation locale de la conservation du moment cinétique se réduit à :

 i jk  k j dV = 0  i jk  k j = 0
D

Ce qui implique que les éléments du tenseur des contraintes vérifient la propriété :

k j =  jk pour j k

 le tenseur des contraintes est symétrique.

97
Cours MDF Aroua Aouadi Enit 2022-2023
 III.4. Conservation de l’énergie
III. 4. 1. Energie – Travail - Chaleur
III. 4. 1. 1. Energie

L’énergie totale d’un système matériel, notée e0 est la somme de son énergie cinétique (ec) et de son énergie interne
(e). La forme intégrale de l’énergie totale s’écrit comme la somme de densité massique de l’énergie interne (e) et de
la densité massique de l’énergie cinétique (v²/2) :

v²
e0D =  e0 dV =   (e +
2
) dV
D D

Le premier principe de la thermodynamique énonce que l'énergie d'un système matériel est égale à la somme du
travail des forces extérieures et de la chaleur échangée avec le milieu extérieur.

98
Cours MDF Aroua Aouadi Enit 2022-2023
 III.4. Conservation de l’énergie
III. 4. 1. Energie – Travail - Chaleur
III. 4. 1. 2. Travail des forces
Le système reçoit un travail produit par les forces surfaciques et les forces volumiques au cours d’un intervalle de temps
Dt. Le travail peut être exprimé en termes de puissances par l’équation suivante :

DWeD
Dt
=  ( v .T).n dS +   f .v dV
S D

Avec la puissance des forces extérieures : PFext = Fext . v

Ainsi la puissance des forces extérieures dans un volume D : PeD =  ( v .T).n dS +   f .v dV
S D

99
Cours MDF Aroua Aouadi Enit 2022-2023
 III.4. Conservation de l’énergie
III. 4. 1. Energie – Travail - Chaleur
III. 4. 1. 3. Echange de chaleur
Le transfert de chaleur se fait essentiellement par 3 modes :
- La conduction : c’est un mode d’échange par contact où le transfert d’énergie se fait par diffusion moléculaire
ou diffusion (agitation) thermique de proche en proche.
- Le rayonnement : c’est mode d’échange à distance, résultant d’interaction d’ondes électromagnétiques avec la
matière.
- La convection : c’est un mode mixte, l’échange de la chaleur se fait par le mouvement de fluide.

Les échanges radiatifs (par rayonnement) sont évalués par une densité volumique r représentant le taux d’apport
(ou de retrait) d’énergie interne par unité de masse et par unité de temps (quantité de chaleur échangée par
rayonnement avec le milieu extérieur par unité de temps et par unité de masse) . Ainsi, une partie D du milieu reçoit
ou émet une puissance calorifique échangée par rayonnement :
DQ rD
Dt
=   r dV
D

100
Cours MDF Aroua Aouadi Enit 2022-2023
 III.4. Conservation de l’énergie
III. 4. 1. Energie – Travail - Chaleur
III. 4. 1. 3. Echange de chaleur
Les échanges par conduction sont représentés par le flux d’un vecteur courant de chaleur q (ou densité de flux de
chaleur) : la conduction entraîne un échange avec l'extérieur d’une puissance calorifique :

DQcD
Dt
=−  q.n dS
S

Les échanges par convection sont représentés par le terme convectif de la dérivée particulaire de
l’énergie totale e0. La dérivée particulaire de l’énergie totale s’écrit :

De0D   v²  v²
Dt
=  t   (e + )  dV +
2    (e +
2
) v.n dS
D S

101
Cours MDF Aroua Aouadi Enit 2022-2023
 III.4. Conservation de l’énergie
III. 4. 2. Forme intégrale de la conservation d’énergie
La dérivée matérielle de l’énergie totale d’une partie D du milieu (volume D) est égale à la somme de la puissance
des forces extérieures et de la puissance calorifique échangée avec l’extérieur :

D  v² 
Dt    (e + )  dV =
2   ( v .(T.n) dS +   f .v dV +   r dV −  q.n dS
D S D D S

  v²  v²

t   (e + 2 )  dV = −   (e + 2 ) v.n dS +  ( v .(T.n) dS +   f .v dV +   r dV −  q.n dS
D S S D D S

102
Cours MDF Aroua Aouadi Enit 2022-2023
 III.4. Conservation de l’énergie
III. 4. 3. Forme locale de la conservation d’énergie

  v²  v²
t    (e + )  dV = −
2    (e + 2
) v.n dS +  ( v .(T.n) dS +   f .v dV +   r dV −  q.n dS
D S S D D S

En utilisant le théorème de divergence :

  v²   v² 
 t   (e + )  dV +
2   .  (e + ) v  dV =
 2    .(T.v) dV +   f .v dV −   .q dV +   r dV
D D D D D D

Le bilan local s’écrit donc en fonction des dérivées partielles comme suit :

D v²   v²   v² 
 (e + ) =   (e + )  +  .   (e + ) v =  .(T.v) +  f .v −  .q +  r
Dt 2 t 2   2 

103
Cours MDF Aroua Aouadi Enit 2022-2023
 III.4. Conservation de l’énergie
III. 4. 5. Equation de bilan d’énergie cinétique et d’énergie interne
La projection de l’équation de la loi fondamentale de la dynamique sur la vitesse :
Dv
v = ( .T).v +  f .v
Dt
On peut écrire le premier membre à gauche de l’équation dans le cas d’écoulement incompressible et en utilisant la
dérivée eulérienne de la vitesse comme :
Dv v  v² v² 
v =  . v +  ( v .  v) . v = ( ) + . ( v ) = (ec ) + . (e c v)
Dt t t 2 2 t
dv d(v² / 2)
NB : v =
dt dt
Le premier membre à droite de l’équation, ( .T).v , peut être développé comme :

( .T).v =  . (T.v ) − T: ( v )

D v²  v² v²
  ( ) = ( ) +  . ( v) =  . (T.v ) − T: ( v ) +  f . v
Dt 2 t 2 2
104
Cours MDF Aroua Aouadi Enit 2022-2023
 III.4. Conservation de l’énergie
III. 4. 5. Equation de bilan d’énergie cinétique et d’énergie interne
1 1
( . v ) = e + ω =
( v + ( v) T ) + ( v − ( v) T )
2 2
1  1 
T: ( v ) = T:  ( u + ( u) T )  + T:  ( u − ( u) T )  = T:e + T:ω = T:e
2  2 
symétrique Anti symétrique

0
l’équation de bilan de l’énergie cinétique s’écrit donc comme

D v²
 ( ) =  . ( v .T) − T: (e) +  f . v
Dt 2
En intégrant sur un volume de contrôle D :
D v²
Dt  (
2
) dV =  v . (T .n) dS −  T: (e) dV +   f . v dV
D S D D

105
Cours MDF Aroua Aouadi Enit 2022-2023
 III.4. Conservation de l’énergie
III. 4. 5. Equation de bilan d’énergie cinétique et d’énergie interne
Bilan de l’énergie totale
D v²
Dt   (e +
2
) dV =  v . (T .n) dS +   f . v dV +   r dV −  q.n dS
D S D D S

Bilan de l’énergie cinétique

D v²
Dt   ( ) dV =
2  v . (T .n) dS −  T: (e) dV +   f . v dV
D S D D

Bilan de l’énergie interne

Le bilan d’énergie interne est obtenu en retranchant le bilan d’énergie cinétique du bilan d’énergie totale :

D
Dt   e dV =   r dV −  q.n dS +  T: (e) dV
D D S D

106
Cours MDF Aroua Aouadi Enit 2022-2023
 III.4. Conservation de l’énergie
III. 4. 5. Equation de bilan d’énergie cinétique et d’énergie interne
Équation locale de l’énergie totale
D v²
 (e + ) =  . (T . v) +  f . v +  r −  .q
Dt 2

Équation locale de l’énergie cinétique

D v²
 ( ) = . (T . v) − T: (e) +  f . v
Dt 2

Équation locale de l’énergie interne

De
 =  r −  .q + T: (e)
Dt

107
Cours MDF Aroua Aouadi Enit 2022-2023
 III.5. Entropie et second principe

Le second principe de la thermodynamique introduit la température et l’entropie : En tout point d’un milieu continu
on peut associer un champ scalaire, positif ou nul, noté T(x,t) : T est la température absolue de la particule qui
passe en x à t. A chaque partie du milieu on associe l’entropie SD, grandeur additive d’ensembles, (à t fixé), et sa
densité massique s (M,t) telle que :
SD = 
D
 s dV

Selon le second principe, La dérivée matérielle de l’entropie d’une partie D du milieu est égale à la puissance
calorifique par degré échangée avec l’extérieur :
DS D 1 r
Dt
=− S T
q .n dS + 
D

T
dV

En introduisant le taux de production d’entropie par unité de masse s (associé au phénomène d’irréversibilité),
l’inégalité fondamentale peut se mettre sous la forme :
DS D 1 r
Dt
=− ST
q .n dS + 
D
 dV +
T 
s dV
D
s  0

108
Cours MDF Aroua Aouadi Enit 2022-2023
 III.5. Entropie et second principe
III. 5. 1. Bilan d’entropie
DS D 1 r
Dt
=− S T
q .n dS + 
D

T
dV + 
D
s dV s  0

En variable d’Euler :

 1 r
t 
D
 s dV = − 
S
 s v .n dS − 
S
T
q .n dS + 
D

T
dV +  
D
s dV s  0

L’équation de bilan local d’entropie s’écrit :

Ds s q r
 =  ( + v . s) = −  . +  + s s  0
Dt t T T

109
Cours MDF Aroua Aouadi Enit 2022-2023
 III.6. lois constitutives et positions de problème
III. 6. 1. Position de problème
Considérons les équations de l’entropie et de l’énergie interne :
Ds s q r
 =  ( + v . s) = −  . +  + s s  0
Dt t T T

De
 =  r −  .q + T: (e)
Dt

En faisant la soustraction entre les deux équations ci-dessus et en éliminant le rayonnement, on peut exprimer la
production d’entropie sous la forme :

 Ds De p 1 1
s = (T − −  .v) − 2 q .  T + τ : e  0
T Dt Dt  T T

110
Cours MDF Aroua Aouadi Enit 2022-2023
 III.6. lois constitutives et positions de problème
III. 6. 2. Hypothèse de l’état local
L’énergie interne est exprimée en fonction de l’entropie comme : de = Q − p dV  de = T ds − p dV

En utilisant le théorème de conservation de la masse et les notions de la thermodynamique, la divergence de la vitesse

1 DV 1 D
v peut exprimée comme un taux de dilatation en fonction du volume V :  . v = = −
V Dt  Dt

L'hypothèse de l'état local (voir cours) conduit à formuler la relation de Gibbs en terme des dérivées matérielles :
De Ds
 = T − p . v
Dt Dt

Compte tenu des équations ci-dessus, le taux de production d’entropie s’écrit :

 Ds De p 1 1 1 1
s = (T − −  .v ) − q .  T + τ :e  0 s = − q .  T + τ :e  0
T Dt Dt  T 2
T T 2
T

111
Cours MDF Aroua Aouadi Enit 2022-2023
 III.6. lois constitutives et positions de problème
III. 6. 2. Hypothèse de l’état local
Le taux de production d’entropie (équation ci-dessous) indique que l'irréversibilité dans le fluide résulte des processus
de conduction thermique et du travail de déformation du milieu par les contraintes visqueuses.
1 1
s = − 2
q .  T + τ :e  0
T T
Le taux de production d'entropie peut s'exprimer en fonction du taux de dissipation intrinsèque et de dissipation
thermique : T  s = 1 +  2

Avec 1 = τ : e  0 La dissipation intrinsèque

1
2 = − q . T  0 La dissipation thermique
T

112
Cours MDF Aroua Aouadi Enit 2022-2023
 III.6. lois constitutives et positions de problème
III. 6. 3. Fluide classique : lois de Newton et de Fourrier
III. 6. 3. 1. Loi de Newton
Le fluide est un milieu dont le tenseur des contraintes visqueuses, appelé aussi le tenseur de cisaillement, s'exprime en
fonction du tenseur des vitesses de déformations.
Le fluide classique ou fluide Newtonien est un milieu de constitution homogène et isotrope pour lequel le tenseur des
contraintes visqueuses  est une fonction linéaire de taux de déformation e : Loi de Newton
2
τ = 2  e +  d . v δ Avec d = ( − )  0
3
 est la viscosité dynamique et d est la viscosité volumique de dilatation ou deuxième viscosité rencontrée dans le cas
des fluides compressibles.
En considérant l’équation de conservation de la masse, la loi de Newton se réduit dans le cas de fluide incompressible à:
τ = 2e
III. 6. 3. 2. Loi de Fourrier
La loi de Fourier pour un milieu homogène et isotrope est une loi linéaire exprimant le vecteur courant de chaleur en
fonction du gradient de température :
q = −T Où  est le coefficient est le coefficient de conduction de la chaleur
113
Cours MDF Aroua Aouadi Enit 2022-2023
 III.6. lois constitutives et positions de problème
III. 6. 3. Fluide classique : lois de Newton et de Fourrier
En tenant compte de la loi de Newton : τ = 2  e +  d . v δ , la dissipation intrinsèque : 1 = τ : e  0 devient :
On sait que : e : δ =  .v 1 = τ : e = 2  e : e +  d (  . v ) 2

1
Compte tenu de la loi de Fourrier : q = −   T, la dissipation thermique  2 = − q .  T  0 devient :
T
 2 =  ( T) 2

La viscosité dynamique  est responsable de la dissipation par cisaillement.

La viscosité volumique d est associée à une dissipation liée à la déformation volumique du fluide.
 
On introduit la viscosité cinématique  :  = et la diffusivité thermique k : k =
  Cp

Le nombre de Prandtl compare la propriété de diffusion de quantité de mouvement du milieu à la propriété de

diffusion thermique : k
Pr =

114
Cours MDF Aroua Aouadi Enit 2022-2023
 III.6. lois constitutives et positions de problème
III. 6. 3. Fluide classique : lois de Newton et de Fourrier : lois de comportements
6 équations pour ij (loi de Newton) et 3 équations pour le vecteur courant de chaleur q (loi de Fourrier)

III. 6. 4. Equations d’état :

3 équations d’état : voir le cours chapitre 3

115
Cours MDF Aroua Aouadi Enit 2022-2023
 III.7. Retour sur les équations du fluide classique
III. 7. 1. Equations du fluide classique à propriétés physiques constantes
En utilisant les lois constitutives, à propriétés physiques constantes, les équations de bilans de masse, de quantité de
mouvement et d’énergie interne deviennent :

Conservation de la masse +  . ( v) = 0
t
Conservation de la quantité de mouvement
Dv v  v
a =  = +  ( v . ) v = +  . ( vv ) = −  p +  .τ +  f
Dt t t
Dv v  v
a =  = +  ( v . ) v = +  . ( vv ) = −  p +  .  2  e +  d  . v  +  f
Dt t t

Or ( )
 .  2  e +  d . v  =  .  ( v + ( v) T ) +  (d  . v )

D’autre part, on sait que  . ( v + ( v) T ) =  2 v +  ( . v ) =  v +  ( . v )

116
Cours MDF Aroua Aouadi Enit 2022-2023
 III.7. Retour sur les équations du fluide classique
III. 7. 1. Equations du fluide classique à propriétés physiques constantes
Conservation de la quantité de mouvement
v  v
+  .( vv ) = −  p +   v + (  +  d )  (. v ) +  f
Dv
a =  = +  ( v . ) v =
Dt t t
Dv v  v
Dans le cas de fluide incompressible :  a =  = +  ( v . ) v = +  . ( vv ) = −  p +   v +  f
Dt t t

Conservation de l’énergie interne

De
 =   2 T +  r − p  . v + 2  e: e + d ( . v ) 2
Dt

117
Cours MDF Aroua Aouadi Enit 2022-2023
IV. STATIQUE DES FLUIDES

118

Cours MDF Aroua Aouadi Enit 2022-2023

 IV. 1. Statique des fluides : définitions et applications
IV. 1. 1. Définitions
Un fluide est en état statique si il n’y a pas de mouvement relatif entre ses particules :
❑ Fluides au repos (immobile).
❑ Fluides uniformément accélérés.

De point de vue mécanique :

❑ Il n’y a pas de contraintes dues aux frottements entre particules.

❑ Les forces mises en jeu sont d’une part les forces de surface dues à la pression. D’autre part, le fluide est soumis
à son poids (force volumique).

119
Cours MDF Aroua Aouadi Enit 2022-2023
 IV. 1. Statique des fluides : définitions et applications
IV. 1. 2. Applications

Barrage
Réservoir de stockage

Sous-marin
Stabilité d’un bateau
120
Cours MDF Aroua Aouadi Enit 2022-2023
 IV. 2. Equation fondamentale de la statique
IV. 2. 1. Formulation de l’équation fondamentale de la statique
A partir de l’équation de conservation de quantité de mouvement :
v  v
+  .( vv ) = −  p +   v + (  +  d )  (. v ) +  g
Dv
a =  = +  ( v . ) v =
Dt t t

Etat statique : Il n’y a pas de mouvement relatif entre les particules fluides (fluide immobile, au repos) :

 a = −  p + g

a x
 a x , a y , a z sont des constantes.
Avec a = a y

a z

121
Cours MDF Aroua Aouadi Enit 2022-2023
 IV. 3. Statique des fluides incompressibles
IV. 3. 1. Cas d’un fluide au repos dans un champs de pesanteur
a = −  p + g
Le fluide est immobile (au repos) :

a =0  0 = − p +  g Loi de l’hydrostatique pour un fluide immobile

0

 p = g   p = 0
−  g

p p p
=0 =0 = −g  p ( x, y, z ) = p(z)
x y z

dp
= −g Equation différentielle à résoudre pour un fluide immobile
dz (pour connaître la pression en tout point)

p +  g z = C te 122
Cours MDF Aroua Aouadi Enit 2022-2023
 IV. 3. Statique des fluides incompressibles
IV. 3. 1. Cas d’un fluide au repos dans un champs de pesanteur

p +  g z = C te

C te = p 0 +  g z s0 p 0 désigne la pression atmosphérique à la surface libre à la côte z s0 à x = 0.

p = p 0 −  g (z − z s0) = p 0 +  g h

123
Cours MDF Aroua Aouadi Enit 2022-2023
 IV. 3. Statique des fluides incompressibles
IV. 3. 1. Cas d’un fluide au repos dans un champs de pesanteur

p +  g z = C te
Conséquences
Les surfaces isobares, c'est-à-dire d'égale pression, sont des plans horizontaux.
En statique des fluides (pas de mouvement), la pression P ne dépend que de l’altitude z.

124
Cours MDF Aroua Aouadi Enit 2022-2023
 IV. 3. Statique des fluides incompressibles
IV. 3. 1. Cas d’un fluide au repos dans un champs de pesanteur

p +  g z = C te

Conséquences

pA +  g zA = pB +  g zB

p B − p A =  g (z A − z B ) =  g h

h pB  pA et zA  zB
Avec

p B = p A +  g (z A − z B ) = p A +  g h

p Bas = p haut +  g h

125
Cours MDF Aroua Aouadi Enit 2022-2023
 IV. 3. Statique des fluides incompressibles
IV. 3. 1. Cas d’un fluide au repos dans un champs de pesanteur
p +  g z = C te

Conséquences

p
On pose : =h
g

pB p
− A = h B − h A = zA − zB
g g

126
Cours MDF Aroua Aouadi Enit 2022-2023
 IV. 3. Statique des fluides incompressibles
IV. 3. 1. Cas d’un fluide au repos dans un champs de pesanteur
Conséquences
Les surfaces de séparation de liquides non miscibles sont des plans horizontaux, ce sont des surfaces d'égales
pressions appelées aussi surfaces isobares.
Si la masse liquide est constituée par plusieurs liquides de densités différentes, ces liquides se superposent par ordre
de densité croissante. Les surfaces de niveau toujours perpendiculaire à la direction de la gravité donc horizontales.
Par conséquent, les surfaces de séparation des différents liquides sont toujours horizontales.

127
Cours MDF Aroua Aouadi Enit 2022-2023
 IV. 3. Statique des fluides incompressibles
IV. 3. 2. Cas d’un fluide en accélération uniforme dans un champs de pesanteur
a = −  p + g

−  a x

 p = −  a y

−  a z −  g

p p p
 = −  ax = −  ay = −  az −  g
x y z

128
Cours MDF Aroua Aouadi Enit 2022-2023
 IV. 3. Statique des fluides incompressibles
IV. 3. 2. Cas d’un fluide en accélération uniforme dans un champs de pesanteur
IV. 3. 2. 1. Cas d’un fluide en accélération uniforme unidirectionnelle dans un champs de pesanteur
Le fluide s’accélère suivant la direction de x, exemple : un chariot avec une accélération longitudinale suivant l’axe (x)
a x −  a x
a = −  p + g  
a = 0  p = 0
0 −  g
 
p p p
= −  ax =0 = − g  p ( x, y, z ) = p(x, z)
x y z

p p C’est l'équation différentielle à résoudre pour connaître la pression en tout

 dp = dx + dz = −  a x dx −  g dz point du fluide en accélération uniforme. Les isobares (p(x, z) = Cte) sont
x z
des droites : les surfaces isobares sont des plans.

p 0 désigne la pression atmosphérique à la

 p = −  ax x −  g z + C
te
p = p 0 −  a x x −  g (z − z s0) surface libre à la côte z à x = 0.
 s0
te
= p 0 +  g z s0 a
C z = z s = z s0 − x x L’équation de la surface libre est obtenue en
g
posant p = p 0 (z = z s , x)
129
Cours MDF Aroua Aouadi Enit 2022-2023
 IV. 4. Statique des fluides compressibles
IV. 4. 1. Gaz parfait dans des conditions isothermes
L’équation fondamentale de la statique des fluides : dp
= − g
dz
L’équation d’état d’un gaz parfait :

nRT RT pM 1
pV = nRT  p = = = rT  = = p
V M RT rT

dp M M  M 
=−p g  Ln p = − g z + C te  p (z) = C te
exp −
 RT g z 
dz RT RT  

Si p = p 0 et z = z 0 :

 M 
p (z) = p 0 exp  − g (z − z0 )
 RT 

130
Cours MDF Aroua Aouadi Enit 2022-2023
 IV. 5. Forces hydrostatiques
IV. 5. 1. Force de pression exercée sur une plaque plane
IV. 5. 1. 1. Cas général
La résultante des pressions élémentaires est appelée “poussée totale” et notée F[N]. Elle s'applique sur chaque face
d'une surface S [m²] quelconque placée verticalement, horizontalement ou inclinée.
Dans le cas général, cette résultante s’exprime comme suit :

dF = − p n dS et p = p 0 +  g h
 F =  − p n dS  F = − n  ( p 0 +  g h ) dS
s s
 
 F = − n  p 0 S +  g  h dS 
 
 s 
Le terme  h dS représente le moment statique de la surface AB par rapport à ox, qui est défini comme :  h dS = h G S
s s

h G : Profondeur du barycentre de la surface (profondeur du centre de gravité)

131
Cours MDF Aroua Aouadi Enit 2022-2023
 IV. 5. Forces hydrostatiques
IV. 5. 1. Force de pression exercée sur une plaque plane
IV. 5. 1. 2. Détermination de point d’application d’une force hydrostatique
Pour déterminer le point d’application noté A (selon la figure ci-dessous) d’une force hydrostatique F s’exerçant sur
une surface S, Il faut calculer le moment de cette force par rapport à un point O quelconque (de préférence le point
O appartient à la surface S), puis identifier ce moment à la résultante des moments élémentaires par rapport à ce
même point O :

OA  F = s OM  dF

132
Cours MDF Aroua Aouadi Enit 2022-2023
 IV. 5. Forces hydrostatiques
IV. 5. 1. Force de pression exercée sur une plaque plane
IV. 5. 1. 2. Détermination de point d’application d’une force hydrostatique
Cas d’une paroi verticale
OA  F = s OM  dF  OA . F e1  e 2 = s OM . dF e1  e 2  OA . F = s OM . dF
e3 e3

F = Fhydro − Fatm  F = (p0 +  g hG ) S − p 0 S

h
h hG =
 F =  g hG S =  g S 2
2

Au fait :
h h2

h
dF =  g h dS =  g h dh L  F = gL h dh =  g L = gS
0 2 2 Pression variable
h

OA =
 0
h  g h dS
 OA =
2
h = hc
h 3
gS
2 133
Cours MDF Aroua Aouadi Enit 2022-2023
 IV. 5. Forces hydrostatiques
IV. 5. 1. Force de pression exercée sur une plaque plane
IV. 5. 1. 2. Détermination de point d’application d’une force hydrostatique
Cas d’une paroi inclinée
G : xG et A : x c

OA . F = s OM . dF
 g x G sin  S x c = 
s
 g x 2 sin  dS

xc =
 s
x 2 sin  dS
=
 s
x 2 dS
x G sin  S xG S


s
x 2 dS = I : Le moment d'inertie de la paroi par rapport à la droite horizontale de trace O

IG
Le théorème d’Huygens permet d’écrire : I = I G + S x G2 xc = x G +
S xG

134
Cours MDF Aroua Aouadi Enit 2022-2023
 IV. 5. Forces hydrostatiques
IV. 5. 1. Force de pression exercée sur une plaque plane
IV. 5. 1. 2. Détermination de point d’application d’une force hydrostatique
Retour sur la paroi verticale
h
 
2
x dS x 2 L dx I
xc = s
= 0
=
xG S xG S xG S

xc =
 s
x 2 dS
=

s
x 2 L dx
xG S xG S

h h3 h3 2
I= 0 x dx = xc = =
2
3
h
3 h 3
3
2

135
Cours MDF Aroua Aouadi Enit 2022-2023
 V. DYNAMIQUE DES
FLUIDES PARFAITS

136

Cours MDF Aroua Aouadi Enit 2022-2023

 V. 1. Forces agissant sur une particule fluide
La loi fondamentale de la dynamique appliquée à un élément de fluide :
Dv Dv
dF = dFs + dFV = dm =  dV
Dt Dt

Avec dFV =  g dV et dFs = −  p dV

Ainsi la loi fondamentale de la dynamique appliquée à un élément de fluide s’écrit :


Dv
= −  p + g
D 
Avec = + v .
Dt Dt t

v
 +  (v .) v = −  p +  g Equation d’Euler qui décrit le mouvement d’un fluide parfait
t

137
Cours MDF Aroua Aouadi Enit 2022-2023
 V. 2. Equation d’Euler
v v
 +  (v .) v = −  p +  g 1
+ (v .) v = −  p + g
t t 
On note (u, v, w) les composantes du vecteur vitesse dans un repère cartésien, (x, y, z) où z est l’axe vertical orienté
vers le haut de sorte que l’accélération de la pesanteur s’écrive :

g = − g iz

La projection de l’équation d’Euler :

  u u u u  p
 p 
  t + u + v + w  = −

 xa = −    x
x   x y z 
    v
 p v v v  p

 ya = − 
  + u + v + w  = −
y    t  x  y  z  y

 p   w w w w  p

 za = − − g   + u + v + w  = − − g
 z    t x y z  z

138
Cours MDF Aroua Aouadi Enit 2022-2023
 V. 2. Equation d’Euler

Exercices et applications de l’équation d’Euler : chapitre 5 du cours

(diapos 12 à 17)

139
Cours MDF Aroua Aouadi Enit 2022-2023
 V. 3. Equation de Bernoulli
V. 3. 1. De l’équation d’Euler à l’équation de Bernoulli
L’équation de Bernoulli est une forme d’intégration scalaire de l’équation d’Euler en écoulement de fluide parfait. En
considérons les hypothèses suivantes :

• Un écoulement permanent : =0
t
• Un fluide parfait : τ = 0

• Un fluide incompressible
v
L’équation d’Euler s’écrit dans ce cas :  +  (v .) v = −  p +  g  (v .) v = −  p +  g
t

Soit un repère cartésien (O, x y, z) tel que l’axe Oz est dirigé selon la verticale vers le haut. Le champ de pesanteur
g dérive d’un potentiel scalaire de gravité (énergie potentielle). Mathématiquement, on écrit :
g = − g iz

 g = −  ( g z)

140
Cours MDF Aroua Aouadi Enit 2022-2023
 V. 3. Equation de Bernoulli
V. 3. 1. De l’équation d’Euler à l’équation de Bernoulli

 (v .) v = −  p +  g  g = −  ( g z)
1
D’autre part le terme convectif ( v .  ) v peut s’écrire sous la forme : (v .) v =  ( v v ) − v  (  v )
2

1 rot v  v
 ( v2 )
2
1
L’équation d’Euler s’écrit donc comme :  (  v 2 ) +  p +  ( g z) =  v  (  v) = 2  v  
2
2

Si l’écoulement est irrotationnel :  = 0

1 1 2
 (  v 2 + p +  g z) = 0  v + p +  g z = C te
2 2
Equation de Bernoulli pour un écoulement permanent
irrotationnel d’un fluide parfait incompressible
141
Cours MDF Aroua Aouadi Enit 2022-2023
 V. 3. Equation de Bernoulli
V. 3. 2. Différentes écritures de l’équation de Bernoulli
1 2 Les terme ont la dimension d’une pression (Pa)
 v + p +  g z = C te
2

P v2
Equation généralisée de Bernoulli +z+ = C te Les terme ont la dimension d’une longueur (m.ce)
g 2g

1 2 p
v + + g z = C te Les terme ont la dimension d’une énergie par unité de masse
2 

142
Cours MDF Aroua Aouadi Enit 2022-2023
 V. 3. Equation de Bernoulli
V. 3. 3. Représentation graphique de l’équation de Bernoulli

Ligne de charge totale horizontale

143
Cours MDF Aroua Aouadi Enit 2022-2023
 V. 3. Equation de Bernoulli
V. 3. 3. Représentation graphique de l’équation de Bernoulli

Ligne de charge totale horizontale

144
Cours MDF Aroua Aouadi Enit 2022-2023
 V. 3. Equation de Bernoulli
V. 3. 3. Représentation graphique de l’équation de Bernoulli
Différence entre tube piézométrique et tube de Pitot

v2
2g
p H
h = ρg + z

Tube piézométrique et tube de Pitot

v2 p v2
+ +z=H +h=H
2g g 2g

145
Cours MDF Aroua Aouadi Enit 2022-2023
 V. 3. Equation de Bernoulli
V. 3. 4. application du théorème de Bernoulli
Ecoulement à travers un orifice (exemple : vidange d’un réservoir – théorème de Toricelli
On applique le théorème Bernoulli entre les deux point A et B de même ligne de
courant : La charge H est constante suivant une même ligne de courant : HA = HB
𝐐 = 𝐪𝐯
v A2 v B2
pA + ρ g z A + ρ = pB + ρ g z B + ρ
2 2

p A = p B = P atm vA ≃ 0

S >>> ac La vitesse au point B : la vitesse de vidange (la vitesse théorique)

vB = 2 g ZA − ZB = 2 g H0 Théorème de Toricelli

Q réel = C q Q théorique = C c S 2 g H

146
Cours MDF Aroua Aouadi Enit 2022-2023
VI. BILAN DE QUANTITÉ DE
MOUVEMENT - THÉORÈME
D’EULER

147

Cours MDF Aroua Aouadi Enit 2022-2023

 VI. 1. Conservation de la quantité de mouvement
La quantité de mouvement contenue dans un volume de contrôle Vs s’écrit : M =
  v dV
Vs

En appliquant la relation fondamentale de la dynamique, la variation de quantité de mouvement est égale à la

somme des forces extérieures (forces volumiques et surfaciques) agissant sur le volume de fluide.

  v dV =  F
dM d
= ext = P + R
dt dt
Vs

P =   g dV
Vs
Forces volumiques (Poids)

R =  Σ n dS
Sc
Forces surfaciques

148
Cours MDF Aroua Aouadi Enit 2022-2023
 VI. 1. Conservation de la quantité de mouvement
dF f
En utilisant le théorème de transport de Reynolds :
dt
= 
Vc
t
dV +  f v . n dS
Sc

L’équation de la variation de quantité de mouvement devient :

 ( v )
  v dV = P + R   ( v) v . n dS
dM d
= P + R = dV +
dt dt
Vs Vc
t Sc

L’équation ci-dessus indique que la somme des forces agissant sur le volume de contrôle est égale à somme de la
variation de QDM dans le volume de contrôle et du flux net de QDM sur la surface qui limite le volume de contrôle.

Pour les écoulements permanents, la variation de la quantité de mouvement se réduit à :

P + R =  ( v) v . n dS
Sc

 ( v) v . n dS = −  p n dS +  τ . n dS +   g dV
Sc Sc Sc Vc
149
Cours MDF Aroua Aouadi Enit 2022-2023
 VI. 2. Théorème d’Euler
Si le volume de contrôle est un tube de courant caractérisé par deux
sections S1 et S2. S est la surface latérale. Le tube de courant correspond
à plusieurs cas d’exemples : un écoulement en charge dans un tronçon
d’une conduite ; un jet d’eau dans l’air atmosphérique...


S1 + S2 + S
( v ) v . n dS = P + R

 ( v ) v− .vn
S1
1 1 1 dS +  ( v ) v
S2
2 2 . n 2 dS +  ( v ) v
S
. n dS = P + R
1
v2 0

  v 2 v 2 S2 −  v1 v1 S1 = P + R
D’autre part, la conservation de la masse se traduit par une conservation de débit :

Q m1 = Q m 2 = Q m    v
S2
2 . n 2 dS = −   v
S1
1 . n1 dS   v 2 S2 =  v1 S1

150
Cours MDF Aroua Aouadi Enit 2022-2023
 VI. 2. Théorème d’Euler
 v 2 v 2 S2 −  v1 v1 S1 = P + R

Q m1 = Q m 2 = Q m    v
S2
2 . n 2 dS = −   v
S1
1 . n1 dS   v 2 S2 =  v1 S1

Théorème d’Euler Q m ( v 2 − v1 )=P+R

151
Cours MDF Aroua Aouadi Enit 2022-2023