Université d’Angers - LPA

CH III Dynamique des Fluides

Il s’agira
’ i iicii d
de ttenir
i compte
t des
d différentes
diffé t forces
f agissant
i t sur les
l
particules fluides en mouvement. Une description quantitative de
l’écoulement pourra ainsi être déduite d’équations fondamentales
locales.
locales

1 - Forces de surface - Tenseur des contraintes

Pour un fluide réel (visqueux) en mouvement, les forces de surface
ne sont p
plus seulement normales à la surface : il existe des
contraintes tangentielles dues à la viscosité (frottements).
r
En un point M d’une
d une surface dS
dS, la force de n
surface s’exprime comme : r r r
dF = Tn dS M Tn
r
ù Tn est la
où l contrainte
t i t quii s’exerce
’ sur la
l
r
surface de normale n. dS
Université d’Angers - LPA CH III - Dynamique

Considérons une surface ⊥ à ll’axe

axer x. r
y
La normale à cette surface est : n = ex
σ yx r
Tx La contrainterexercée sur cette surface est
alors notée Tx et peut se décomposer
r comme : r
n σ xx r r r
Tx = σ xx ex + σ yx ey + σ zx ez
σ zx x
On remarque que les composantes σ yx et σ zx
z sont des composantes tangentielles : on les
notera plutôt τ yx et τ zx pour les distinguer
de la composante normale σ xx .
y
On peut de même considérer la surface σ yy r
⊥ à l’axe y. On a ainsi la contrainte : Ty
r r r r r
Ty = τ xy ex + σ yy ey + τ zy ez n
τ xy
Et pour la surface ⊥ à l’axe z la contrainte
s’exprime
s exprime : τ zy x
r r r r
Tz = τ xz ex + τ yz ey + σ zz ez z
Université d’Angers - LPA CH III - Dynamique

Considérons maintenant une surface dont l’orientation est

quelconque. Dans le repère cartésien, sa normale peut se
décomposer
p en : r r r r
n = nx ex + ny ey + nz ez
Dans ce cas, la contrainte s’exerçant sur cette surface s’exprime
comme : r r r r
Tn = nx Tx + ny Ty + nz Tz
En développant, on obtient :
r r r r
nx (σ xx ex + τ yx ey + τ zx ez ) ⎧(nxσ xx + nyτ xy + nzτ xz ) ex
r r r r ⎪ r
Tn = + ny (τ xyy ex + σ yy ey + τ zyy ez ) = ⎨(nxτ yyx + nyσ yy + nzτ yyz ) ey
r r r ⎪(n τ + n τ + n σ ) er
+ nz (τ xz ex + τ yz ey + σ zz ez ) ⎩ x zx y zy z zz z

Ce qui revient à effectuer le produit d’une

d une matrice par la normale :

⎛ σ xx τ xy τ xz ⎞⎛ nx ⎞ Tenseur
r ⎜ ⎟⎜ ⎟ r r
Tn = ⎜ τ yx σ yy τ yz ⎟⎜ ny ⎟ ⇒ Tn = T n des
⎜τ Contraintes
⎝ zx τ zy σ zz ⎟⎠⎜⎝ nz ⎟⎠
 Université d’Angers - LPA CH III - Dynamique

Ce tenseur des contraintes peut alors se décomposer en la somme

d’un tenseur sphérique et d’un tenseur de trace nulle :

T = α I + T'
⎛1 0 0⎞ ⎛ σ 'xx τ xy τ xz ⎞ Tenseur des
⎜ ⎟ ⎜ ⎟
α ⎜0 1 0⎟ ⎜ τ yx σ 'yy τ yz ⎟ contraintes
t i t
⎜0 0 1⎟ ⎜τ τ zy σ 'zz ⎟⎠ de viscosité
⎝ ⎠ ⎝ zx
⎧σ xx = α + σ 'xx
⎪ σ 'xx +σ 'yy +σ 'zz = 0 = −3α + σ xx + σ yy + σ zz
Ainsi : ⎨σ yy = α + σ 'yy et
⎪σ = α + σ '
⎩ zz zz
⇒ α = 13 (σ xx + σ yy + σ zz ) = 1
3
Tr(T )
r r r r r r
Tn = T n = α I n + T ' n = αn + T ' n
r r r r r r r
dF = TndS = αndS + T ' ndS ⇒ d F = − pn d S + T ' n d S
force normale force de pression
⇒
à la surface hydrostatique où p = − 13 Tr(T )
α = −p et Tr(T ') = 0
Université d’Angers - LPA CH III - Dynamique

2 - Application du Principe Fondamental de

la Dynamique
De la même façon que l’on a établi l’équation fondamentale de la
statique des fluides, nous allons établir l’équation fondamentale de
la Dynamique en appliquant le principe fondamental de la
dynamique à un élément de volume fluide en mouvement. Il nous
faut donc faire le bilan des forces s’exerçant en surface et en
volume.
r r r r
dv
dF = dFS + dFV = ρ dV
dt
forces de forces de
surface volume
r
⇓ ª poids : ρ dV
dVg
r r r
dFS x ex + dFS y ey + dFS z ez

Intéressons-nous à la composante suivant y

 Université d’Angers - LPA CH III - Dynamique
r r
n = ez
r r
[τ yz ]z + dz dxdy
dF = T n dS
r
− [τ yx ]x dydz n = −e
r
− [σ yy ]y dxdz x [σ yy ]y + dy dxdz

z r r [τ yyx ]x + dx dydz r r
n = −ey n = ey
r r
n = ex
⎛σ τ xy τ xz ⎞
− [τ yz ]z dxdy r ⎜ xx ⎟r
y dF = ⎜ τ yx σ yy τ yz ⎟ n dS
r r ⎜τ τ zy σ zz ⎟⎠
x ⎝ zx
n = −ez
⎧ ∂τ yx
[τ yx ]x + dx = [τ yx ]x + dx
dFS y = ([τ yx ]x + dx − [τ yx ]x ) dydz ⎪ ∂x
⎪ ∂σ yy
⎪
+ ([σ yy ]y + dy − [σ yy ]y ) dxdz où ⎨[σ yy ]y + dy = [σ yy ]y + dy
⎪ ∂y
+ ([τ yz ]z + dz − [τ yz ]z ) dxdy ⎪ ∂τ yz
⎪⎩ [τ yz ]z + dz = [τ yz ]z + ∂z
∂z
dz
∂τ yx ∂σ yy ∂τ yz
⇒ dFS y = dxdydz + dydxdz + dzdxdy
∂x ∂y ∂z
Université d’Angers - LPA CH III - Dynamique

∂τ yx ∂σ yy ∂τ yz
dFS y = dxdydz + dydxdz + dzdxdy
∂x ∂y ∂z
dV dV dV
⎛ ∂τ yx ∂σ yy ∂τ yz ⎞
dFS y = ⎜⎜ + + ⎟⎟ dV
⎝ ∂x ∂y ∂z ⎠

Et par analogie, suivant les autres axes :

⎛ ∂σ ∂τ ∂τ ⎞ ⎛ ∂τ ∂τ ∂σ zz ⎞
dFS x = ⎜⎜ xx + xy + xz ⎟⎟ dV dFS z = ⎜⎜ zx + zyy + ⎟⎟ dV
⎝ ∂x ∂y ∂z ⎠ ⎝ ∂x ∂y ∂z ⎠

Ce qui
q i pe
peutt s’écrire
s’éc i e comme :

⎛ σ xx τ xy τ xz ⎞
r r⎜ ⎟ r
dFS = ∇ ⎜ τ yx σ yy τ yz ⎟ dV = ∇T dV
⎜τ τ zy σ zz ⎟⎠
⎝ zx
Université d’Angers - LPA CH III - Dynamique

r r r r
dv
Par conséquent : dF = dFS + dFV = ρ dV
dt
r
r r dv
⇒ ∇T dV + ρ g dV = ρ dV
dt
r
r r dv
⇒ ∇T + ρ g = ρ
dt
On peut alors faire apparaître les deux
parties du tenseurs des contraintes : T = − pI + T '
⎛ p 0 0⎞
r r⎜ ⎟ r r r
∇T = −∇⎜ 0 p 0 ⎟ + ∇T ' = −∇p + ∇T '
⎜ 0 0 p⎟
⎝ ⎠
r
dv r r r équation fondamentale
Soit finalement : ρ = −∇p + ∇T ' + ρ g de la dynamique
dt (équation locale)
Université d’Angers - LPA CH III - Dynamique

3 - Fluide newtonien
et équation de Navier-
Navier-Stokes

Par définition, les fluides « newtoniens » sont ceux pour lesquels les
composantes du tenseur des contraintes de viscosité T ' dépendent
linéairement des composantes du tenseur des taux de déformation
pure e.

Remarque :
Une rotation pure n’engendre aucune déformation : par conséquent
il n’y a pas de contrainte. C’est pourquoi T ' et ω ne sont pas liés.
é

Remarque :
Tous les fluides que l’on étudiera pourront être considérés newtoniens.
Université d’Angers - LPA CH III - Dynamique

Considérons les éléments tensoriels de e :

1⎜
⎛ ∂v i
∂v j ⎞ on rappelle que ce tenseur
eij = + ⎟
2⎜ ⎟ est symétrique, car eij = e ji
⎝ ∂x j ∂xi ⎠
On admettra alors que pour un fluide isotrope, les éléments
éé tensoriels
de T ' et e sont liés par la relation suivante :

σ 'ij = 2μ eij + μ' (exx + eyy + ezz )δ ij

viscosité viscosité de symbole de

dilatation Kronecker
∂v x ∂v y ∂v z rr
On remarque que : exx + eyy + ezz = + + = ∇v
∂x ∂y ∂z
rr
Donc, si le fluide est incompressible, on a ∇v = 0 et dans ce cas :

σ 'ij = 2μ eij ⇒ T ' = 2μ e

Université d’Angers - LPA CH III - Dynamique

Reprenons l’équation fondamentale de la dynamique :

r
dv r r r
ρ = −∇p + ∇T ' + ρ g
dt
Pour un fluide incompressible
p newtonien,, cette équation
q devient donc :
r
dv r r r
ρ = −∇p + 2μ ∇e + ρ g
dt
r
Explicitons le terme ∇e :

r ⎛ ∂eij ⎞ r ⎛ ∂v i ∂v j ⎞
∇e = ∑i ⎜⎜ ∑j ∂x ⎟⎟ ei où eij = 1
2
⎜
⎜ ∂x
+
∂x i
⎟
⎟
⎝ j ⎠ ⎝ j ⎠

r ⎛ 1 ∂ 2v i ∂ 2v j ⎞ r
⇒ ∇e = ∑i ∑j ⎜⎜ 2 ∂x 2 + 1
2
⎟ ei
∂x j ∂x i ⎟⎠
⎝ j
Université d’Angers - LPA CH III - Dynamique

r ⎛ 1 ∂ 2v i ∂ 2v j ⎞ r
∇e = ∑i ∑j ⎜⎜ 2 ∂x 2 + 1
2
⎟ ei
∂x j ∂x i ⎟⎠
⎝ j

⎛ ∂ 2
vi ⎞r ∂ ⎛ ∂v j ⎞r r 1 r rr
= 2∑ ∑
1 ⎜ ⎟ ei + 1
∑i ∂x ⎜∑ ⎟ ei = 1
Δv + 2 ∇(∇v )
⎜
i ⎝ j ∂x j
2 ⎟ 2 ⎜ j ∂x ⎟ 2
⎠ i ⎝ j ⎠
car fluide
rr incompressible
Δv i ∇v
r r rr
Δv ∇(∇v )
r r
Il reste
t alors
l : ∇e = 1
2
Δv
Ce qui conduit à :
r r
dv r r r dv r r r
ρ = −∇p + 2μ ∇e + ρ g ⇒ ρ = −∇p + μ Δv + ρ g
dt dt
équation fondamentale de la
dynamique pour un fluide équation de Navier-Stokes
newtonien incompressible (équation locale)
 Université d’Angers - LPA CH III - Dynamique

r dérivée instantanée
dv r r r
ρ = −∇p + μ Δv + ρ g dérivée convective
dt r r
rr r
dv ∂v
ª dérivée particulaire : = + (v ∇) v
dt ∂t
r
⎛ ∂v r r r⎞ r r r
⇒ ρ⎜ + (v ∇) v ⎟ = −∇p + μ Δv + ρ g
⎝ ∂t ⎠
r r
En posant g = −gez , la projection de l'équation de Navier-Stokes sur
les 3 axes du repère cartésien donne :

⎧ ⎛ ∂v x ∂v x ∂v x ∂v x ⎞ ∂p ⎛ ∂ 2v x ∂ 2v x ∂ 2v x ⎞
⎪ ρ ⎜⎜ + vx + vy + vz ⎟⎟ = − + μ ⎜⎜ 2
+ 2
+ 2 ⎟
⎟
⎪ ⎝ ∂t ∂x ∂y ∂z ⎠ ∂x ⎝ ∂x ∂y ∂z ⎠
⎪⎪ ⎛ ∂v y ∂v y ∂v y ∂v y ⎞ ∂p ⎛ ∂ 2v y ∂ 2v y ∂ 2v y ⎞
⎨ ρ ⎜⎜ + vx + vy + vz ⎟⎟ = − + μ⎜ + + ⎟
⎜ 2 2 2 ⎟
⎪ ⎝ ∂t ∂x ∂y ∂z ⎠ ∂y ⎝ ∂x ∂y ∂z ⎠
⎪ ⎛ ∂v ∂v z ∂v z ∂v z ⎞ ∂p ⎛ ∂ 2v z ∂ 2v z ∂ 2v z ⎞
⎪ ρ ⎜⎜ z
+ vx + vy + vz ⎟⎟ = − + μ ⎜⎜ 2
+ 2
+ 2 ⎟
⎟ − ρg
⎪⎩ ⎝ ∂t ∂x ∂y ∂z ⎠ ∂z ⎝ ∂x ∂y ∂z ⎠
Université d’Angers - LPA CH III - Dynamique

4 - Equation fondamentale en volume

et théorème d’Euler
d’Euler
a) Théorème de transport de Reynolds

Considérons
C idé une grandeur
d
r scalaire
l i fonction
f ti des
d coordonnées
d é ded
l’espace et du temps : f (r , t )

Sur ler volume Vs d d’un

un système de particules fluides,
fluides ll’intégrale
intégrale
de f (r , t ) s’écrit : r
F =
∫∫∫VS
f (r , t ) dV

Si l’on souhaite évaluer les variations de F dans le temps, il nous

faut
aut calculer
ca cu e :
dF
dF d r
dt
=
dt ∫∫∫ VS (t )
f (r , t ) dV

Le problème est qu’ici Vs est une fonction du temps :

en effet le système de particules fluides est en mouvement.
Université d’Angers - LPA CH III - Dynamique

La parade consiste à utiliser un volume fixe Vc (volume de contrôle),

délimité par une surface Sc (surface de contrôle) à travers laquelle
on pourra comptabiliser le flux de f :

∫∫∫
S dF d r
= f (r , t ) dV
c dt dt VS (t )

∫∫∫ ∫∫
dF ∂f rr
= dV + f v.n dS
dt Vc ∂t Sc

VS(t+dt) dérivée
VS(t) = Vc dérivée
locale convective

Variations instantanées
é de f
dans le volume de contrôle

Flux de f à travers la
surface de contrôle
Université d’Angers - LPA CH III - Dynamique

r r
Prenons un exemple concret : f (r , t ) ≡ ρ (r , t ) masse volumique
r r
Dans ce cas : F =
∫∫∫ VS
f (r , t ) dV ≡
∫∫∫ VS
ρ (r , t ) dV = M masse de VS

La variation de masse s’exprime donc comme :

dM ∂ρ rr
dt
=
∫∫∫ Vc ∂t
dV +
∫∫
Sc
ρ v.n dS

variation de masse due variation de masse due au

aux variations flux massique à travers la
instantanées de ρ surface de contrôle

Utilisons la formule d’Ostrogradski pour transformer l’intégrale de

surface en une intégrale de volume :

r r rr r r r r
∫∫ S
A.n dS =
∫∫∫ V
∇A dV ⇒
∫∫ Sc
( ρv )).n dS =
∫∫∫
Vc
∇( ρv ) dV
Université d’Angers - LPA CH III - Dynamique

Il s’en suit :

dM ∂ρ rr ∂ρ r r
dt
=
∫∫∫ Vc ∂t
dV +
∫∫Sc
ρ v.n dS =
∫∫∫
Vc ∂t
dV +
∫∫∫
Vc
∇( ρv ) dV

dM ⎛ ∂ρ r r ⎞
⇒
dt
=
∫∫∫ ⎜
Vc ⎝ ∂t
+ ∇( ρv ) ⎟dV
⎠
équation de continuité

Ainsi, si l’écoulement est conservatif, on retrouve bien :

dM
dM
=0 autrement dit, il n’y a pas de variation de masse
dt
Université d’Angers - LPA CH III - Dynamique

b) Théorème
é è de transport appliqué
é à la quantité
é de mouvement

La q
quantité de mouvement d’un système
y fluide de volume VS
s’écrit : r
∫∫∫ VS
ρv dV

Or, le principe fondamental de la dynamique nous dit que la dérivée

par rapport au temps de la quantité de mouvement doit être égale
à la somme des forces agissant sur le système :

d r r r
dt ∫∫∫ VS
ρv dV = R + P
forces de volume (poids)
r
∫∫∫
forces de surface
ρ g dV
r
∫∫
VS
T n dS
SS

calcul trivial
calcul d’accès
difficile ⇒ il est généralement plus simple de calculer
la variation de quantité de mouvement
Université d’Angers - LPA CH III - Dynamique

r r d ⎡d
∑ ⎤r
r
Ainsi : R + P =
dt ∫∫∫ VS
ρv dV =
i
⎢⎣ dt ∫∫∫ VS
ρvi dV ⎥ ei
⎦

∂( ρvi )
∑ ∫∫∫
⎡ ⎤r
rr
=
i
⎢⎣ Vc ∂t
dV +
∫∫Sc
( ρv i ) vn dS ⎥ ei
⎦

⎡ ∂( ρv i ) r ⎤ ⎡ r ⎤ rr
=
∫∫∫ ⎢ ∑
Vc ⎣ i ∂t
ei ⎥ dV +
⎦ ∫∫
⎢ ∑
Sc ⎣ i
( ρv i )ei ⎥ vn dS
⎦
r
∂( ρv ) r rr
=
∫∫∫ Vc ∂t
dV +
∫∫
Sc
( ρv ) vn dS

dérivée instantanée débit de quantité de

de la quantité de mouvement à travers
mouvement la surface de contrôle

Pour un écoulement stationnaire,

stationnaire r r r rr
la dérivée instantanée est nulle, et donc : R+P =
∫∫Sc
( ρv ) vn dS
Université d’Angers - LPA CH III - Dynamique

c) Théorème
é è d’Euler

Appliquons
pp q le résultat p
précédent au cas d’un tube de courant :
S2
r
v2
S1 r r
r v1 n2
n1
r r
nl vl Sl

On supposera la vitesse constante en tout point d’une même

section ; l’écoulement
é sera supposé
é stationnaire.
r r r rr
R+P =
∫∫ S1 + S2 + Sl
( ρv ) vn dS
r r r r r r r r r r r
R+P =
∫∫S1
( ρv1 ) v1n1 dS +
∫∫ S2
( ρv2 ) v2n2 dS +
∫∫
Sl
( ρv l ) v l nl dS
− v1 v2 0
r r
− ρv1v1S1 ρv2v2S2
Université d’Angers - LPA CH III - Dynamique

r r r r
Il reste : R + P = − ρv1v1S1 + ρv2v2S2

Or, on sait que dans un tube de courant le débit massique est

Or
conservé :
qm = ρv1S1 = ρv2S2
r r r r
D’où
D où le résultat simple suivant : R + P = qm(v2 − v1 )
Théorème d’Euler

Remarque : ce théorème permet d’obtenir simplement la

résultante des forces de surface (notamment
les forces de frottement) sans avoir à recourir
au tenseur des contraintes.