Université Moulay Ismail
Ecole Nationale Supérieure d’Arts et Métiers ENSAM-Meknès
Département d’Energétique
Cycle Ingénieur - 3ème année
Pr. Mohammed SRITI
1 M. Sriti Cours Mécanique des fluides - ENSAM
Sommaire
Chap1: Définitions et généralités
Chap2: Statique des Fluides
Chap3: Cinématique des Fluides
Chap4: Dynamique des Fluides Parfaits
Chap5: Dynamique des Fluides Réels
Chap6: Pertes de Charges
2 M. Sriti Cours Mécanique des fluides - ENSAM
Chapitre 3
Cinématique des Fluides
Sommaire
1. Définitions
2. Description Lagrangienne
3. Description Eulérienne
4. Dérivée matérielle ou particulaire
5. Flux et divergence du vecteur vitesse
6. Circulation et rotationnel du vecteur vitesse
7. Loi de conservation de la masse
4 M. Sriti Cours Mécanique des fluides - ENSAM
1. Définitions :
La cinématique des fluides est l’étude du mouvement d’un fluide sans tenir
compte des forces qui le produisent.
On s’intéresse, au mouvement du fluide par rapport à un solide fixe ou en
mouvement.
En mécanique, il existe deux modes de description du mouvement:
Point de vue Euler,
Point de vue Lagrange.
Pour étudier le mouvement, on choisit généralement un référentiel
cartésien, de base
i , j,k d’origine O, fixe par rapport à l’un des
solides étudiés.
5 M. Sriti Cours Mécanique des fluides - ENSAM
2. Description Lagrangienne :
2.1 Définition:
L’observateur suit chaque particule fluide dans son mouvement à partir
d’une position initiale à un instant donné. x3
Une particule fluide qui,
M o, to
M ,t
à l’instant t=to se trouve en Mo de
a x
coordonnées a a1 , a2 , a3 OM 0 O
x2
à l’instant t>to se trouve en M de
coordonnées x x1,x2 ,x3 OM x1
a1 , a 2 , a3 , t : sont les variables de Lagrange nécessaires à la description,
M. Sriti elles sont
Cours Mécanique indépendantes
des fluides - ENSAM entre elles.
6
Pour connaître l’évolution du fluide, il faut déterminer les fonctions:
x x a, t x i f i a j ,t
x 1 f 1 a 1 , a 2 , a 3 , t
Soit:
x 2 f 2 a 1 , a 2 , a 3 , t
x f a , a , a , t
3 3 1 2 3
2.2 Propriété de base:
x x a , t permet d’associer à chaque particule fluide à l’instant initial to une
et une seule particule fluide à l’instant t
La transformation est une fonction bijective, pour chaque particule à l’instant
t lui correspond un et un seul point d’origine .
x x a , t a a x , t
7 M. Sriti Cours Mécanique des fluides - ENSAM
2.3 Trajectoire des particules fluides:
Pour connaître la trajectoire, il suffit de fixer le point d’origine et suivre
l’évolution d’une particule fluide au fil du temps, cela revient à fixer (a1, a2,
a3) et à faire varier le temps t.
La trajectoire est le lieu des positions successives d’une particule fluide au
cours du temps. Elles sont directement fournies par
xi fi a j , t à aj fixés
2.4 Vitesse et accélération:
8 M. Sriti Cours Mécanique des fluides - ENSAM
Vitesse: z
M,t
dx
Soit une particule fluide entre l’instant t au M’
x
t t
point M et l’instant (t+Δt ) au point M’:
x dx
o
x OM y
x dx OM MM ' x
dx MM ' : déplacement de la particule de M à M’ pendant Δt
La vitesse est définie par: dx
v (a ,t)
dt
a j
dx i
Soit en indices: vi (a j , t)
dt a j
9 M. Sriti Cours Mécanique des fluides - ENSAM
Accélération:
L’accélération d’ une particule fluide est:
dv
(a ,t)
dt
a j
dv i d 2 xi
Soit en indices: i ( a j ,t )
dt
a j
dt 2
a j
10 M. Sriti Cours Mécanique des fluides - ENSAM
3. Description Eulérienne :
3.1 Définition:
L’observateur est placé en un point M fixe de l’espace de coordonnées xi
dans le référentiel attaché au solide de référence, et regarde passer les
particules fluides devant lui.
Pour des instants différents, ce n’est pas la même particule qui occupe la
position x (M )
On notera la propriété F au point M à
l’instant t, telle que:
F x,t
11 M. Sriti Cours Mécanique des fluides - ENSAM
(x1, x2, x3, t) sont appelées les variables d’Euler :
v1 x1 , x2 , x3 ,t
La vitesse d’une particule fluide s’écrit: v x ,t v2 x1 , x2 , x3 ,t
v x , x , x , t
3 1 2 3
3.2 Lien entre les deux descriptions :
Soit une propriété physique F du fluide (scalaire, vecteur, tenseur) :
En Lagrangienne F a,t
En Eulérienne
F x,t
Par égalité des deux propriétés, on trouve:
F a , t F x a , t , t
Ou F x , t F a x , t , t
12 M. Sriti Cours Mécanique des fluides - ENSAM
La connaissance d’une des deux descriptions, conduit à la connaissance de
l’autre.
Dans chaque problème, on utilise la description la plus adaptée.
3.3 Trajectoires:
L’ Ensemble des positions occupées par une particule fluide donnée, est une solution:
dx
v ( x ,t ) dxi vi ( x ,t ) dt
dt
L’ équation de la Trajectoire est:
dx1 dx 2 dx 3
dt
v1 ( x , t ) v 2 ( x , t ) v3 ( x , t )
On aboutit à trois équations du premier ordre et trois constantes d’intégration
arbitraires (aj ).
13 M. Sriti Cours Mécanique des fluides - ENSAM
3.4 Lignes de courant:
Ligne dont chacun de ses points le vecteur vitesse de la particule fluide est tangent
à un instant t0 fixé.
MM ' d x // v ( x , t o )
L’équation des lignes de courant est:
dx1 dx2 dx3
v1 ( x , to ) v2 ( x , to ) v3 ( x , to )
Cas particulier : Ecoulement stationnaire
Les trajectoires et les lignes de courant sont confondues.
Dans la description Eulérienne, les grandeurs sont indépendantes du temps.
Dans la description Lagrangienne, les grandeurs dépendent en général du
temps.
14 M. Sriti Cours Mécanique des fluides - ENSAM
f ( x , t )
t
4. Dérivée matérielle ou particulaire
4.1 Dérivée matérielle d’une grandeur scalaire : f ( x , t )
M, t
x
df ( x , t ) f ( x x , t t ) f ( x , t ) M’, t+t
lim
dt t 0 t x
df f ( x x, t t ) f ( x, t t ) f ( x, t t ) f ( x, t ) x x
lim
dt t 0 t t
O
f ( x x , t t ) f ( x , t t ) grad f ( x , t t ) x 0(x 2 )
grad f ( x , t t ). v ( x , t ) t 0(t 2 ) 0(v 2 t 2 )
grad f ( x , t t ).v ( x , t ) t 0(t 2 )
Et f ( x , t t ) f ( x , t ) f ( x , t )
lim
t 0
t
t
df f
Soit enfin : v grad f
15
dt t
4.2 Dérivée matérielle d’une grandeur vectorielle : Accélération ( x , t )
En description Eulérienne , l’accélération est définie par :
dv ( x , t ) v ( x x , t t ) v ( x , t )
(x, t) lim
dt t 0 t
dv v v2
Soit : rot v v grad
dt t 2
Translation Rotation Déformation
En général, le mouvement d’une particule fluide est la superposition d’une
translation, d’une rotation et d’une déformation .
16 M. Sriti Cours Mécanique des fluides - ENSAM
Soit en indice muet, l’accélération s’écrit :
dv i v i v i
i vj avec i=1, 2, 3
dt t x j
dv 1 v1 v1 v1 v1
1 v v v
x3
1 2 3
dt t x 1
x 2
dv 2 v2 v2 v2 v2
Soit en composantes : 2 v v v
x3
1 2 3
dt t x 1
x 2
dv 3 v3 v3 v3 v 31
3 v v v
x3
1 2 3
dt t x 1
x 2
17 M. Sriti Cours Mécanique des fluides - ENSAM
4.3 Dérivée matérielle d’une intégrale de volume :
On considère un domaine D occupé par un fluide
de volume V et de surface S, à l’instant t. Soit f une •M’
n
fonction scalaire et soit l’intégrale:
v
I
V
f dV M
•
dS V
dI
dt
1
t 0 t V V
lim f (
x , t t ) dV V f (
x, t ) dV V
S
n étant le vecteur unitaire normal à S, orienté vers l’extérieur de manière
que le flux fluide quitte le volume V à l’instant t vers l’extérieur avec une
vitesse v , c. à d. que dV soit positif et négatif dans le cas contraire.
V+V étant le volume occupé à l’instant t+ t par les particules contenues dans
V à l’instant t.
18 M. Sriti Cours Mécanique des fluides - ENSAM
La dérivée matériel de I peut s’écrire :
dI 1
lim f ( x , t t ) dV f ( x , t t ) dV
dt t 0 t V V V
f ( x , t t ) dV f ( x , t ) dV
V V
L’élément de volume peut s’écrire :
dV MM ' n dS v ( x , t ) t n dS
Après développement, on obtient :
dI
dt t f
V
dV
S
f v n dS
Premier terme : Partie instationnaire ou dérivée matérielle;
Second terme : Partie convective ou flux de ( f v ) à travers la surface S.
19 M. Sriti Cours Mécanique des fluides - ENSAM
Remarque: Si V ne renferme pas de surface de discontinuité pour f et v
f
L’équation en forme locale s’écrit :
t
V
f dV
V t
dV
Par théorème d’Ostrogradski, on a :
S
f v n dS div ( f v ) dV
V
dI f
d’où :
dt
V
t div ( f v ) dV
Par analogie, pour une intégrale de fonction vectorielle , telle que :
Soit : J A dV
V
dJ
On a:
dt t AdV A (v n ) dS
V
S
20 Partie instationnaire Partie convective
5. Flux et divergence du vecteur vitesse
5.1 Flux du vecteur vitesse::
Soit S une surface quelconque,
Pour un élément de surface dS de normale n
Le volume élémentaire sortant de dS pendant le
temps dt est : dV v n dt dS
v : Le vecteur vitesse
dV
Le volume élémentaire par unité de temps est : dqV v n dS
dt
Le flux du vecteur vitesse est le volume total sortant de S par unité de
temps, c’est aussi le débit volumique à travers cette surface:
qv v n dS vn dS
S S
21 M. Sriti Cours Mécanique des fluides - ENSAM
Remarque1: Le flux du vecteur v est égal à la masse de fluide sortant de (S)
par unité de temps est le débit masse tel que:
qm
S
v n dS
Remarque2: Si S est une surface de courant, donc le débit volumique qv ou le
débit massique qm est nul, car :
v .n 0
Remarque3:
Le flux est une grandeur scalaire qui dépend de l’orientation de S.
D’une manière générale, le flux d’un vecteur quelconque A à travers une
surface S s’écrit:
A n dS
S
22 M. Sriti Cours Mécanique des fluides - ENSAM
5.2 Divergence du vecteur vitesse::
Le flux du vecteur vitesse à travers une surface
fermée S contenant le volume V à un instant t
donné est :
qv n dS
S
v
Le flux a la propriété d’additivité.
Remarque:
Si qv > 0, on dit que D contient un puits de débit qv
Si qv < 0, D contient une source de débit qv.
Si qv = 0, D ne contient ni puits ni sources, ou ils se compensent.
23 M. Sriti Cours Mécanique des fluides - ENSAM
Si on admet que les sources et les puits sont répartis dans le domaine D .
La densité volumique de sources ou de puits pour un élément de volume
autour d’un point M est :
dq v
div v Soit dqv div v . dV
dV
Donc le débit volumique total pour le volume V du domaine D est:
qv div v .dV v n dS
V S
C’est la formule d’Ostrogradski pour la divergence définissant
le passage d’une intégrale de volume à celle de surface.
24 M. Sriti Cours Mécanique des fluides - ENSAM
Ax A y Az
div A
x y z
25 M. Sriti Cours Mécanique des fluides - ENSAM
6. Circulation et rotationnel du vecteur vitesse
6.1 Circulation du vecteur vitesse:
L’intégrale curviligne d’un champ vectoriel le
long d’un arc de courbe orienté (C ) par la n
relation est :
v (r ) . dr v ( r ). .dl
t t
C C
Si la courbe est fermée, l’intégrale s’appelle la circulation du vecteur
A
le long de (C ), notée:
v (r ) dr v ( r ) dl
t t
C C
26 M. Sriti Cours Mécanique des fluides - ENSAM
Soit une surface (S) limitée par le contour (C),
(C)
si on découpe S en trois surfaces S1, S2 et S3, (C3, S3)
limitées par les contours C1, C2 et C3. (S) (C1,S1)
La circulation possède la propriété d’additivité. (C2,S2)
Γ Γ1 Γ 2 Γ 3
6.2 Rotationnel du vecteur vitesse:
En utilisant un raisonnement analogue à celui appliqué pour le flux, on peut
définir la densité de circulation par unité de surface élémentaire dS par:
d
dS
t
rot v . n
n est la normale à S orientée en concordance avec le sens choisi pour (C ).
27 M. Sriti Cours Mécanique des fluides - ENSAM
Soit :
t
d rot v . n . dS
La circulation du vecteur vitesse s’écrit :
t t
v (r ) . dr rotv . n.dS
C S
d
rotv Indique la direction suivant laquelle dS
est maximum.
La circulation s’applique pour tout vecteur quelconque A
t
t
A(r ) . dr rot A . n.dS
C C
28 M. Sriti Cours Mécanique des fluides - ENSAM
7. Loi de conservation de la masse
7.1 Formulation intégrale : dS
dm dV
On considère un domaine fluide D de volume (S)
V limité par une surface S et caractérisé par (V)
une densité ρ.
Une particule fluide de masse dm, de volume
dV et de surface dS
La masse totale contenue dans le domaine D à l’instant t est:
m (t ) V
( x , t ) dV
Par analogie à la dérivée matérielle d’une intégrale de volume (f=), la
dérivée de la masse totale par unité de temps dans le volume V s’écrit :
dm(t ) d
dt
dt
(x, t) dV t dV v.n dS
V V S
29
Hypothèse:
En description Eulérienne, la masse d’un système matériel étant constante
au cours du temps, le domaine matériel D ne renferme ni source ni puits la
dérivée particulaire de la masse est nulle.
La variation de la masse totale est nulle :
dm(t )
dV v.n dS 0
dt V
t S
D’après le théorème de la divergence, on a:
dm ( t )
dt
V t div ( v ) dV 0
C’est la loi de la conservation de la masse en forme intégrale
30 M. Sriti Cours Mécanique des fluides - ENSAM
7.2 Formulation Locale :
pour un volume arbitraire qui ne renferme pas de discontinuité, l’expression
locale de la loi de conservation de la masse, s’écrit:
div ( v ) 0
t
Ou encore: d
div ( v ) 0
dt
Soit en forme d’indices : d ( v j )
0 et i 1 , 2 , 3
t dx j
Pour un fluide incompressible: cte et div(v ) 0
vi
Soit : 0 ; i 1, 2, 3
31 xi
7.3 Opérateurs vectoriels:
-Soit un repère dans une base orthonormée:
-Soit un vecteur V et une fonction scalaire U :
-Soit l’opérateur vectoriel nabla:
Gradient:
Divergence :
32 M. Sriti Cours Mécanique des fluides - ENSAM
Rotationnel:
33
Exemple1:
Un réservoir de volume V = 1 m3 est rempli d’air à la température 20°C et à
pression p0 = 100 kPa. A l’instant initial, une pompe à vide permet d’évacuer
l’air avec un débit volumique Qv = 80L/min. On suppose que l’air est un gaz
parfait et que le système étudié est isotherme. Quel est le temps nécessaire pour
vider l’air jusqu’à une pression de 20 kPa.
Réponse:
V p
t ln 1,2 10 3 s 20 min
Q po
34 M. Sriti Cours Mécanique des fluides - ENSAM