Expérience de cisaillement: Cas d’un solide

Contrainte de cisaillement, t = F/s

U0

Par cette inclinaison, le corps enregistre une déformation de cisaillement constante

Déformation de cisaillement = g = tan(q) = U0/h

Rigidité: G = t
g
[ g ] = 1 ou 100 %
 Expérience de cisaillement: cas d’un liquide
t = F/s
U(t4) U(t5) U(t6) U(t7)
U(t2)U(t3)
U(t1)

Succession d’inclinaison: q(t1); q(t2); q(t3); ……

g(t)= tanq(t)  la déformation de cisaillement dépend du temps 

.g
.g L’unité de est s-1
d’où l’ de la vitesse de déformation , = dg/dt
 Condition d’adhérence ou
condition de non glissement

• La couche du fluide qui est

en contact à la paroi va
l’adhèrer. Elle aura le
même comportement que
la paroi.
• Si la paroi est immobile
(ex: une table) , la
vitesse de cette couche
=0
• Si la paroi est mobile de
vitesse (une spatule),
alors la couche fluide
aura la même vitesse .
 C’est quoi un étalement ?

tgq(t1)
tgq(t2)

L’étalement est le résultat de la création successive de plusieurs couches fluides qui

adhérent à la paroi inférieure. Ainsi, l’épaisseur du corps diminue jusqu’à obtention d’un
étalement complet. Le glissement des couches fluides les uns sur les autres est
« l’écoulement »
 y Profil de vitesse
dans un
´
𝜏 écoulement de
V0 Couette Plan
Plan A mobile
V(y+dy)
dy
V(y)

𝜃
Plan B fixe

est la vitesse déplacement d’une couche fluide à la hauteur y.

est le déplacement d’une couche fluide à la hauteur y

Normalement, le profil de vitesse est non linéaire, toutefois, en considérant la distance

entre les deux parois, h, très petite, on peut alors supposer que le profil est linéaire.
La vitesse de déformation s’écrit :

La déformation est donnée par :

On remplace dans l’expression de la vitesse de déformation:

La vitesse de déformation représente alors un gradient de vitesse

Vitesse de déformation de cisaillement
.
g (s-1)
 Ecoulement de poiseuille: la paroi cylindrique est immobile

Ecoulement de couette rotatif Ecoulement de Couette Plan

 Dispersion & division par cisaillement:
principe du rotor-stator ?

Stator, partie fixe

contenant des
ouvertures

Rotor, partie motrice

générant un écoulement
soit:
1. radial,
2. axial
3. radial et axial
L’homogénéisateur rotor-stator : Quel est
son principe ?

Stator, partie fixe

contenant des
ouvertures

Rotor, partie motrice

générant un
écoulement soit:
1. radial,
2. axial
3. radial et axial
Dispersion d’un agrégat par cisaillement
Division d’une gouttelette par cisaillement
 Domaine de la vitesse de cisaillement en fonction de
l’application pratique
Phénomène Vitesse de cisaillement Application
typique (s−1)
Sédimentation – crémage 10−3 − 10−1 Supensions - crèmes
Agitation, pompage – re- 101 − 103
dispersion des suspensions

Pulvérisation & injection 103 – 104 Suspensions pulvérisables (nasales),

Enrobage par pulvérisation
Suspensions injectables
Etalement par friction 104 − 105 cosmétique, lotion et crème sur la peau

Attention!
Les fluides se rappellent du traitement antérieur!
(Ce sont des fluides à mémoire !)
Prélever un échantillon avec une pipette,
génère une vitesse de cisaillement de :
1000 s-1 pour la glycérine
3000s-1 pour l’eau
Systèmes de mesure du Rhéomètre générant
des cisaillements simples
Paramètre \ Type Couette Plateaux parallèles Cône-plan

Diamètre ⌀ (mm) Cylindre : de 20 mm à 80mm Plan : 8mm – 60 mm Cône : 20 mm – 60mm

Entrefer h (mm) De 0,1 à quelques mm 0,2mm – 3mm Dépend de θ et du ⌀ du cône
Volume d’échantillon (ml) important ⇒ temps de mise en Faible Faible (~ 0,5ml – 2ml)
température assez long
Inerties mécanique et Élevées Très faibles Très faibles
thermique
Vitesse de cisaillement dans Quasi constant si h < R ⇒ mesure Variable ⇒ mesure de η relative Quasi constant si θ ne dépasse pas 4° ⇒
l’entrefer de η absolue possible mesure de η absolue possible
Forts taux de cisaillement Non (à cause de l'inertie) Oui si h faible Oui si θ faible
Gamme de fréquence en Basses fréquences uniquement (à Large ⇒ rhéomètre dynamique Large ⇒ rhéomètre dynamique
oscillation cause de l'inertie) (mesures avec h < 1 mm)
Mesures de faibles η Oui (surface de contact élevée ⇒ Possible avec plan ⌀ 60 Possible avec cône ⌀ 60
bonne sensibilité)
Mesures de fortes η Non Avec plan ⌀ 20 Avec cône ⌀ 20
Echantillons contenant de Oui si h élevé Oui si h élevé Non
grosses particules
Facilité d'emploi Remplissage très délicat, Nettoyage facile Assez sensible à un mauvais
nettoyage parfois difficile remplissage ; nettoyage facile
Rhéomètre MCR 302 Anton-Paar

Moteur

Plan supérieur ,
rotor

Plan inférieur
thermo régulé ,
stator
 Viscosité
(shear viscosity)
C’est le rapport de la sollicitation mécanique au
cisaillement résultant

(Pa.s) Loi de Newton

t (Pa)  (Pa.s)

( s-1)

 ( s-1)
Le Brookfiled et ses
spindles…
 Les viscosimètres pour les liquides de
viscosité constante

• Le viscosimètre brookfield est utilisé

normalement seulement pour les fluide
newtoniens. Leur viscosité est constante,
donc ne dépend pas de la vitesse de
déformation.

• Le viscosimètre à chute de bille qui est

plus précis, plus robuste, sauf nécessité
l’utilisation de liquides transparents,
pour pouvoir voir la bille.

• Le viscosimètre à capillaire est intéressant

pour les solutions de polymère. Ca permet
de déterminer le poids moléculaire,
notamment. Ce type de viscosimètre mesure
une viscosité cinématique. Pour obtenir une
viscosité dynamique il faut connaitre au
préalable la masse volumique.
 Fluides newtonien et non
Newtonien
(Pa.s) Loi de Newton

h ne dépend pas du h dépend du

cisaillement cisaillement

h absolue ou h Apparente
dynamique Fluide Fluide non
Newtonien Newtonien
 Courbe d’écoulement de viscosité:
𝟎 → ∞
10000
Produit plastique (
= Cte 1000 Produit peudoplastique ()
Viscosité apparente Pa.s

100

10
Rég
io ns r
1 héo
fl uidi
fian
te  ∞=𝑪𝒕𝒆
0,1
s

0,01
10-4 10-3 10-2 10-1 100 101 102 103 104 105

Vitesse de cisaillement s-1

 Dynamique des fluides réels
Les forces de Dans le cas d’un
frottement jouent un Objectif: fluide parfait:
rôle significatif  (rappel)

déterminer les Les équations du

La viscosité du équations du
fluide ne pourra mouvement des
mouvement des fluides sont les
plus être négligée. fluides équations d’Euler.

Si Newtoniens 
« fluide parfait »  Les équations de
« fluide réel ».
Navier-Stokes
 Les équations du mouvement d’un
fluide parfait: les équations d’Euler.
Dans le cas des
fluides • Le poids: (extérieur ou intérieur ?): dû à
incompressibles la présence de la pesanteur de la terre.
 = Cte, on • Les forces de pression: agissent
détermine les normalement et vers l’intérieur sur les six
forces faces de l’élément dV
extérieures qui • La force d’inertie proportionnelle à
agissent sur un l’accélération: : représente la résultante
élément de des forces extérieures.
volume
dV = dxdydz:
 Les forces de pression: Agissent normalement et
vers l’intérieur sur les six faces de l’élément dV
Z
⃗
𝐹𝑝 𝑧 (𝑧+𝑑𝑧)
z+dz

⃗
𝐹𝑝 𝑥 (𝑥)
⃗ ⃗
𝐹𝑝 𝑦 ( 𝑦 +𝑑𝑦)
𝐹𝑝 𝑦 ( 𝑦 )
⃗
𝐹𝑝 𝑥 (𝑥+𝑑𝑥)
Y
x, y ,z
y+dy

x+dx ⃗
X 𝐹𝑝 𝑧 (𝑧)
 Les forces de pression: Agissent normalement et
vers l’intérieur sur les six faces de l’élément dV
On détermine la résultante de toutes
ces forces de pression:

⃗ ⃗ ⃗
∑ 𝑭𝒑=−( 𝑭𝒑𝒙(𝒙+𝒅𝒙)−𝑭𝒑𝒙(𝒙))𝒊−(𝑭𝒑𝒚 (𝒚 +𝒅𝒚 )−𝑭𝒑𝒚 (𝒚 ) 𝒋−(𝑭𝒑𝒛 (𝒛+𝒅𝒛)−𝑭𝒑𝒛(𝒛)) 𝒌
⃗
 Les forces de pression: Sur les six faces de l’élément dV

⃗ ⃗ ⃗
∑ 𝑭𝒑=−( 𝑭𝒑𝒙(𝒙+𝒅𝒙)−𝑭𝒑𝒙(𝒙))𝒊−(𝑭𝒑𝒚 (𝒚 +𝒅𝒚 )−𝑭𝒑𝒚 (𝒚 ) 𝒋−(𝑭𝒑𝒛 (𝒛+𝒅𝒛)−𝑭𝒑𝒛(𝒛)) 𝒌
⃗
On applique le théorème des accroissement finis:

𝒇 ( 𝒙 +𝒅𝒙 ) − 𝒇 ( 𝒙 ) 𝝏𝒇 𝝏𝒇
= 𝒇 ’ ( 𝒙 )=  𝒇 ( 𝒙+ 𝒅𝒙 ) − 𝒇 ( 𝒙 )= 𝒅𝒙
𝒅𝒙 𝝏𝒙 𝝏𝒙

on applique ce théorème à la résultante de Force de pression

∑⃗
𝑭𝒑=− ( 𝝏 𝑭𝒑 𝒙
𝝏𝒙
⃗
) (
𝒅𝒙 𝒊 −
𝝏 𝑭𝒑 𝒚
𝝏𝒚
⃗
𝒅𝒚 𝒋 − ) (
𝝏 𝑭𝒑 𝒛
𝝏𝒛
𝒅𝒛 ⃗
𝒌 )
 Les forces de pression: Sur les six faces de l’élément dV

dxdy

dxdz
dydz

On remplace l’expression des forces

de pression par la pression: P=Fp/S

Py est la pression qui règne sur la surface dxdz surface

dxdz orientée suivant l’axe y

Px est la pression qui règne sur la surface dydz surface

dydz orientée suivant l’axe x

Pz est la pression qui règne sur la surface dxdy

dxdy surface orientée suivant l’axe Z
 Les forces de pression: Sur les six faces de l’élément dV

La pression est une grandeur isotrope: Sa valeur ne dépend pas de l’orientation

de la surface sur laquelle règne.

𝑷𝒙= 𝑷𝒚= 𝑷𝒛= 𝑷

dydz dxdz dxdy

On remplace ces trois expressions dans l’équation de la résultante de pression:

∑⃗
𝑭𝒑=− ( ) (
𝝏 𝑭𝒑 𝒙
𝝏𝒙
⃗
𝒅𝒙 𝒊 −
𝝏 𝑭𝒑 𝒚
𝝏𝒚 ) (⃗
𝒅𝒚 𝒋 −
𝝏 𝑭𝒑 𝒛
𝝏𝒛
𝒅𝒛 ⃗
𝒌 )
∑ 𝑭𝒑=−( 𝝏 𝒙 𝒅𝒙 ) 𝒊 −( 𝝏 𝒚 𝒅𝒚 ) 𝒋 − ( 𝝏 𝒛 𝒅𝒛 ) ⃗𝒌
⃗ 𝝏( 𝑷𝒅𝒚𝒅𝒛) ⃗ 𝝏(𝑷𝒅𝒙 𝒅𝒛) ⃗ 𝝏( 𝑷𝒅𝒙 𝒅 𝒚 )

∑ 𝑭𝒑=−( 𝝏 𝒙 𝒅𝒚𝒅𝒛 𝒅𝒙 ) 𝒊 − ( 𝝏 𝒚 𝒅𝒙𝒅𝒛 𝒅𝒚 ) 𝒋 − ( 𝝏 𝒛 𝒅𝒙𝒅𝒚 𝒅𝒛 ) ⃗𝒌

⃗ 𝝏( 𝑷 ) ⃗ 𝝏( 𝑷) ⃗ 𝝏(𝑷 )
 Les forces de pression: Sur les six faces de l’élément dV

⃗
(
𝝏( 𝑷 ) ⃗ 𝝏( 𝑷)
) (
⃗ 𝝏(𝑷 )
) (
∑ 𝑭𝒑=− 𝝏 𝒙 𝒅𝒚𝒅𝒛 𝒅𝒙 𝒊 − 𝝏 𝒚 𝒅𝒙𝒅𝒛 𝒅𝒚 𝒋 − 𝝏 𝒛 𝒅𝒙𝒅𝒚 𝒅𝒛 ⃗𝒌 )
Or:
 ∑ ⃗
𝑭𝒑=− (𝝏𝒙 ) (𝝏𝒚 ) (𝝏𝒛 )
𝝏𝑷
𝒅 𝑽 ⃗
𝒊 −
𝝏𝑷
𝒅𝑽 ⃗
𝒋 −
𝝏𝑷
𝒅𝑽 ⃗
𝒌

(∑ ⃗
𝑭𝒑 )
( 𝝏𝑷 ⃗
) (
𝝏𝑷 ⃗ 𝝏𝑷 ⃗
) ( )
On divise les deux
membres par dV  =− 𝒊− 𝒋− 𝒌
𝐝𝐕 𝝏𝒙 𝝏𝒚 𝝏𝒛

Ceci est l’expression du ? Gradient ⃗

(∑ ) 𝒈𝒓𝒂𝒅 𝑷
𝑭𝒑 / 𝐝𝐕 =−⃗
Conclusion
La résultante des forces de pression statique par unité de volume qui s’exercent
sur un élément de volume d’un fluide se résume en un gradient de la pression
qui règne sur ses facettes.
 Les équations du mouvement d’un
fluide parfait: les équations d’Euler.
Dans le cas des
fluides • Le poids: (extérieur ou intérieur ?): dû à
incompressibles la présence de la pesanteur de la terre.
 = Cte, on • Les forces de pression: agissent
détermine les normalement et vers l’intérieur sur les six
forces faces de l’élément dV
extérieures qui • La force d’inertie proportionnelle à
agissent sur un l’accélération: : représente la résultante
élément de des forces extérieures.
volume
dV = dxdydz:
 Les équations du mouvement d’un
fluide parfait: les équations d’Euler.
Dans le cas d’un fluide incompressible,  = Cte, les forces
extérieures (par unité de volume, dV = dxdydz) sont:

Les forces de La force

Le poids/dV:
pression/dV: d’inertie/dV:

: Les équations d’Euler.

C’est quoi son
expression ?
 Détermination de l’expression de
l’accélération
L’accélération, dans la cadre de la physique du point matériel ou des
solides, est le résultat de la variation de la vitesse par rapport temps:
Les distances inter-particulaires sont considérées constantes dans le
temps.
Dans le cas du liquide, l’accélération tient compte de la
variation de la vitesse :
1. Par rapport temps
2. Par rapport à l’espace: Les distances inter-particulaires
ne sont plus considérées constantes en mouvement.
 Détermination de l’expression de
l’accélération dans le cas des liquides

La vitesse varie alors en fonction du temps, t, et de

la position (x, y et z) 
•
• est la composante de la vitesse suivant l’axe
• est la composante de la vitesse suivant l’axe
• est la composante de la vitesse suivant l’axe
 Détermination de l’expression de
l’accélération dans le cas des liquides
• Ecrivant la différentielle totale d’une des trois composantes: exemple:

En divisant les deux membres par dt:

𝑑𝑡
= ( ) ( ) ( ) ( )
𝑑𝑢 ( 𝑥 , 𝑦 , 𝑧 ,𝑡 ) 𝜕 𝑢 𝑑𝑥 𝜕𝑢 𝑑𝑦 𝜕 𝑢 𝑑𝑧 𝜕 𝑢 𝑑𝑡
+ + +
𝜕 𝑥 𝑑𝑡 𝜕 𝑦 𝑑𝑡 𝜕 𝑧 𝑑𝑡 𝜕 𝑡 𝑑𝑡

Or: , et

Ce qui 𝑑𝑢 ( 𝑥 , 𝑦 , 𝑧 ,𝑡 ) 𝜕 𝑢 𝜕𝑢 𝜕𝑢 𝜕𝑢
implique: = ( 𝑢) + (𝑣 )+ (𝑤)+ ( 1)
𝑑𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑡
 Détermination de l’expression de l’accélération dans
le cas des liquides
Nous avons trouvé l’expression suivante
𝑑𝑢 ( 𝑥 , 𝑦 , 𝑧 ,𝑡 ) 𝜕 𝑢 𝜕𝑢 𝜕𝑢 𝜕𝑢
= ( 𝑢) + (𝑣 )+ (𝑤)+ ( 1)
𝑑𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑡
Qui peut s’écrire sous la forme d’un produit scalaire:

𝑑𝑡 (
𝑑𝑢 ( 𝑥 , 𝑦 , 𝑧 ,𝑡 ) 𝜕𝑢 ⃗ 𝜕 𝑢 ⃗ 𝜕 𝑢 ⃗
=
𝜕𝑥
𝑖+
𝜕𝑦
𝑗 +
𝜕𝑧 ) ⃗ ⃗ ⃗
𝑘 ∙ ( 𝑢 𝑖 +𝑣 𝑗 +𝑤 𝑘 ) +
𝜕𝑢
𝜕𝑡

𝑑𝑢 𝜕 𝑢 ⃗ ⃗
= + 𝑉 ∙ 𝑔𝑟𝑎𝑑 𝑢
𝑑𝑡 𝜕 𝑡
La même chose, on peut écrire pour les deux autres composantes:
𝑑 𝑣 𝜕𝑣 ⃗ ⃗ 𝑑𝑤 𝜕𝑤 ⃗ ⃗
= + 𝑉 ∙ 𝑔𝑟𝑎𝑑 𝑣 et = + 𝑉 ∙ 𝑔𝑟𝑎𝑑 𝑤
𝑑𝑡 𝜕𝑡 𝑑𝑡 𝜕𝑡
 Détermination de l’expression de
l’accélération dans le cas des liquides
l’expression de l’accélération sera alors:

{
𝑑𝑢 𝜕𝑢 ⃗ ⃗
= + 𝑉 ∙ 𝑔𝑟𝑎𝑑 𝑢
𝑑𝑡 𝜕𝑡
⃗
𝛾= 𝑑𝑣 𝜕𝑣 ⃗ ⃗
= + 𝑉 ∙ 𝑔𝑟𝑎𝑑 𝑣
𝑑𝑡 𝜕𝑡
𝑑𝑤 𝜕𝑤 ⃗ ⃗
= + 𝑉 ∙ 𝑔𝑟𝑎𝑑 𝑤
𝑑𝑡 𝜕𝑡

L’accélération L’accélération
locale convective
Les équations du mouvement d’un fluide Réel: Représentation
des contraintes de frottements visqueux.
Z ⃗ 𝒛𝒛 (𝒛 +𝒅𝒛 )
𝝉
Les contraintes de frottements z+dz
visqueux naissent en ⃗ 𝒛𝒚 (𝒛 +𝒅𝒛)
𝝉
écoulement. ⃗ 𝒛𝒙 ( 𝒛 +𝒅𝒛)
𝝉
⃗ 𝒚𝒛 ( 𝒚 +𝒅𝒚 )
𝝉
En statique, elles sont nulles
⃗ 𝒙𝒛 (𝒙+𝒅𝒙 )
𝝉
Sur chaque surface, il existe ⃗ 𝒚𝒚 (𝒚 +𝒅𝒚 )
𝝉
trois composantes:
• 02 de cisaillement ⃗ 𝒙𝒚 (𝒙 +𝒅𝒙)
𝝉
(tangentielles, , i=1,3 ⃗ 𝒚𝒙 (𝒚 +𝒅𝒚 )
𝝉 Y
et j=1,3) x, y ,z ⃗ 𝒙𝒙 (𝒙+𝒅𝒙 )
𝝉
• 01 normale
y+dy
(Elongationnelle, ,
i=1,3)

x+dx
X
Les contraintes qui sont de 09, se transmettent d’une couche à un autre (couches adjacentes:
voisines par frottement visqueux. Ceci est responsable du transfert de quantité de
mouvement (moteur ou résistant).
Les déformations qui résultent des trois contraintes:
02 de cisaillement: Deux inclinaisons (écoulements de cisaillement)
01 normale: Ecoulement Elongationnel

Z z+dz

⃗ 𝒚𝒛 ( 𝒚 +𝒅𝒚 )
𝝉

⃗ 𝒚𝒚 (𝒚 +𝒅𝒚 )
𝝉

⃗ 𝒚𝒙 (𝒚 +𝒅𝒚 )
𝝉 Y
x, y ,z
y+dy

x+dx
X
 Contraintes de cisaillement et
contraintes normales
Un fluide purement visqueux développe des contraintes de
cisaillement. lien exemple visqueux
En écoulement tridimensionnel (3D), il aura 06 composantes de
cisaillement txy, txz, tyx, tyz, tyz et tzy
Un fluide viscoélastique développe des contraintes de cisaillement et des
contraintes normales.
En écoulement tridimensionnel il aura 09 composantes:
• 06 de cisaillement txy, txz, tyx, tyz, tyz et tzy
• et
• 03 contraintes normales txx, tyy, tzz

Lien exemple viscoélastique

 Le tenseur des contraintes
Le tenseur des contraintes de 09 composantes s’écrit sous forme d’une
matrice 3x3.
C’est un tenseur d’ordre 2.
Note:
Un tenseur d’ordre 1: est un vecteur
Un tenseur d’ordre 0 est un scalaire.

[ ]
𝜏 𝑥𝑥 𝜏 𝑥𝑦 𝜏 𝑥𝑧
τ́ = 𝜏 𝑦𝑥 𝜏 𝑦𝑦 𝜏 𝑦𝑧
𝜏 𝑧𝑥 𝜏 𝑧𝑦 𝜏 𝑧𝑧
Les termes de la diagonale:
les composantes normales
Les équations du mouvement dans le cas d’un fluide parfait:
: Les équations d’Euler.

Les équations du mouvement dans le cas d’un fluide réel:

: Les équations du mouvement.

Les équations du mouvement dans le cas d’un fluide réel:

: Les équations du mouvement.
 Calcul de la divergence du tenseur
de contrainte

[ ] [ ]
𝜏 𝑥𝑥 𝜏 𝑥𝑦 𝜏 𝑥𝑧 𝜏 𝑥𝑥 𝜏 𝑥𝑦 𝜏 𝑥𝑧
´ = 𝑑𝑖𝑣 𝜏 𝑦𝑥
𝒅𝒊𝒗 𝝉 𝜏 𝑦𝑦 𝜏 𝑦𝑧 ¿⃗
𝛻 ∙ 𝜏 𝑦𝑥 𝜏 𝑦𝑦 𝜏 𝑦𝑧
𝜏 𝑧𝑥 𝜏 𝑧𝑦 𝜏 𝑧𝑧 𝜏 𝑧𝑥 𝜏 𝑧𝑦 𝜏 𝑧𝑧

[ ]
𝜏 𝑥𝑥 𝜏 𝑥𝑦 𝜏 𝑥𝑧
¿( 𝜕 ⃗
𝜕𝑥
𝑖+
𝜕 ⃗
𝜕𝑦
𝑗+
𝜕 ⃗
𝜕𝑧 )
𝑘 ∙ 𝜏 𝑦𝑥
𝜏 𝑧𝑥
𝜏 𝑦𝑦
𝜏 𝑧𝑦
𝜏 𝑦𝑧
𝜏 𝑧𝑧

() () ()
𝜏 𝑥𝑥 𝜏𝑥 𝑦 𝜏 𝑥𝑧
( 𝜕⃗ 𝜕⃗ 𝜕⃗
𝜕𝑥 𝜕𝑦 𝜕𝑧
𝜏 𝑧𝑥
) (
⃗ 𝜕⃗ 𝜕⃗ 𝜕⃗
𝜕𝑥 𝜕𝑦 𝜕𝑧
𝜏 𝑧𝑦
⃗
) (
𝜕⃗ 𝜕⃗ 𝜕⃗
¿ 𝑖+ 𝑗+ 𝑘 ∙ 𝜏 𝑦𝑥 𝒊+ 𝑖+ 𝑗+ 𝑘 ∙ 𝜏 𝑦𝑦 𝒋+ 𝑖+ 𝑗+ 𝑘 ∙ 𝜏 𝑦𝑧 ⃗𝒌
𝜕𝑥 𝜕 𝑦 𝜕𝑧
𝜏 𝑧𝑧
)
 () () ()
𝜏 𝑥𝑥 𝜏𝑥 𝑦 𝜏 𝑥𝑧
(𝜕⃗ 𝜕⃗ 𝜕⃗
𝜕𝑥 𝜕𝑦 𝜕𝑧 ) (
𝜏 𝑧𝑥
⃗ 𝜕⃗ 𝜕⃗ 𝜕⃗
𝜕𝑥 𝜕𝑦 𝜕𝑧
𝜏 𝑧𝑦
⃗
) (
𝜕⃗ 𝜕⃗ 𝜕⃗
¿ 𝑖+ 𝑗+ 𝑘 ∙ 𝜏 𝑦𝑥 𝒊+ 𝑖+ 𝑗+ 𝑘 ∙ 𝜏 𝑦𝑦 𝒋+ 𝑖+ 𝑗+ 𝑘 ∙ 𝜏 𝑦𝑧 ⃗𝒌
𝜕𝑥 𝜕 𝑦 𝜕𝑧
𝜏 𝑧𝑧
)

( )( )(
𝜕𝜏 𝑥𝑥 𝜕𝜏 𝑦𝑥 𝜕𝜏 𝑧𝑥 ⃗ 𝜕𝜏 𝑥𝑦 𝜕𝜏 𝑦𝑦 𝜕𝜏 𝑧𝑦 ⃗ 𝜕𝜏 𝑥𝑧 𝜕𝜏 𝑦𝑧 𝜕𝜏 𝑧𝑧 ⃗
¿ + + 𝒊+ + + 𝒋+ + + 𝒌
𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑥 𝜕𝑦 𝜕𝑧 )
Finalement la divergence d’un tenseur d’ordre 2 donne un vecteur.
Les équations du mouvement dans le cas d’un fluide réel:
: Les équations du mouvement.

On remplace les expressions des trois composantes de dans les

équations du mouvement
et aussi on remplace les expressions de l’accélération par ses trois
composantes.
Les contraintes de cisaillement  des écoulements de cisaillement
=> vitesses de déformation de cisaillement
Les contraintes normales  des écoulements élongationnels.
=> Vitesses de déformation normale
On doit remplacer les 09 contraintes par leur expression en vitesse
de déformation de cisaillement
On a besoin de la relation entre:
• le tenseur des contraintes
• et le tenseur taux de déformation

´
´ =η 𝐷
𝜏
 Pour un fluide Newtonien: ´
´ =η 𝐷
𝜏
 = Cte, est la viscosité dynamique

Pour un fluide non Newtonien:

est la viscosité apparente qui dépend des 𝜏 ´
´ = η𝑎𝑝𝑝 𝐷
vitesses de déformation de cisaillement:

Dans le cadre de cette matière, TQM, on considérera que le

fluide a un comportement newtonien

[ ] [ ]
𝜏 𝑥𝑥 𝜏 𝑥𝑦 𝜏 𝑥𝑧 𝛾
˙ 𝑥𝑥 𝛾
˙ 𝑥𝑦 𝛾
˙ 𝑥𝑧
𝜏 𝑦𝑥 𝜏 𝑦𝑦 𝜏 𝑦𝑧 =η 𝛾
˙ 𝑦𝑥 𝛾
˙ 𝑦𝑦 𝛾
˙ 𝑦𝑧
𝜏 𝑧𝑥 𝜏 𝑧𝑦 𝜏 𝑧𝑧 𝛾
˙ 𝑧𝑥 𝛾
˙ 𝑧𝑦 𝛾
˙ 𝑧𝑧
´
´ =η 𝐷
𝜏
Les expressions des vitesses de déformation
Exemple de lecture: : Variation par rapport à de
la vitesse dirigée suivant

Finalement, on déduit qu’il existe

des termes de composantes du
tenseur taux de déformation qui
sont les mêmes.

On dit que le tenseur taux de déformation est symétrique

Dans ce qui suit, on remplace toutes les expressions de

dans les équations du mouvement du fluide réel.
 On fera le traitement mathématique sur une seule équation, et les deux autres
seront déduites.

[ ]
𝝏𝒖 2 ⃗ 𝜕𝑣 𝜕𝑢
𝜕(2 − 𝑑𝑖𝑣 𝑉 ) 𝜕 ( + )
𝜕𝑝 𝝏𝒙 3 𝜕𝑥 𝜕𝑦 𝜕𝑤 𝜕𝑢 ¿
=− + + +𝜕 + ¿ ¿𝜕 𝑧
𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑥 𝜕𝑧
 [ ]
𝝏𝒖 2 ⃗ 𝜕𝑣 𝜕𝑢
𝜕(2 − 𝑑𝑖𝑣 𝑉 ) 𝜕 ( + )
𝜕𝑝 𝝏𝒙 3 𝜕𝑥 𝜕𝑦 𝜕𝑤 𝜕𝑢 ¿
=− + + +𝜕 + ¿ ¿𝜕 𝑧
𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑥 𝜕𝑧
Dans le cas d’un fluide incompressible, =Cte.
Dans le cadre de la cinématique: L’équation de continuité décrivant la
conservation de masse s’écrit:

=0
et
=>

[ ( )+
]
𝝏𝒖 𝜕𝑣 𝜕𝑢 𝜕𝑤 𝜕𝑢
𝜕 2 𝜕( + ) 𝜕( + )
𝜕𝑝 𝝏𝒙 𝜕𝑥 𝜕𝑦 𝜕𝑥 𝜕𝑧
=− + +
𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑧

[ ]
𝟐 2 2 2 2
𝜕𝑝 𝝏 𝒖 𝜕 𝑣 𝜕 𝑢 𝜕 𝑤 𝜕 𝑢
=− + 2 𝟐
+ + + +
𝜕𝑥 𝝏𝒙 𝜕 𝑥𝜕 𝑦 𝜕𝑦
2
𝜕 𝑥 𝜕 𝑧 𝜕 𝑧2
 [ ]
𝟐 2 2 2 2
𝜕𝑝 𝝏 𝒖 𝜕 𝑣 𝜕 𝑢 𝜕 𝑤 𝜕 𝑢
=− + 2 𝟐
+ + + +
𝜕𝑥 𝝏𝒙 𝜕 𝑥𝜕 𝑦 𝜕𝑦
2
𝜕 𝑥 𝜕 𝑧 𝜕 𝑧2

[ ]
𝟐 𝟐 2 2 2 2
𝜕𝑝 𝝏 𝒖 𝝏 𝒖 𝜕 𝑣 𝜕 𝑢 𝜕 𝑤 𝜕 𝑢
=− + 𝟐
+ 𝟐
+ + + +
𝜕𝑥 𝝏𝒙 𝝏𝒙 𝜕 𝑥 𝜕 𝑦 𝜕 𝑦2 𝜕 𝑥 𝜕 𝑧 𝜕 𝑧 2

[ ]
𝟐 2 2 𝟐 2 2
𝜕𝑝 𝝏 𝒖 𝜕 𝑢 𝜕 𝑢 𝝏 𝒖 𝜕 𝑣 𝜕 𝑤
=− + 𝟐
+ + + 𝟐
+ +
𝜕𝑥 𝝏𝒙 𝜕𝑦
2
𝜕𝑧
2
𝝏𝒙 𝜕 𝑥𝜕 𝑦 𝜕𝑥 𝜕𝑧

[ ( )]
𝟐 2 2
𝜕𝑝 𝝏 𝒖 𝜕 𝑢 𝜕 𝑢 𝜕 𝝏𝒖 𝜕𝑣 𝜕𝑤
=− + + + + + +
𝜕𝑥 𝝏𝒙
𝟐
𝜕𝑦
2
𝜕𝑧
2
𝜕𝑥 𝝏𝒙 𝜕𝑦 𝜕𝑧

∆𝑢:𝑙𝑒 𝐿𝑎𝑝𝑙𝑎𝑐𝑖𝑒𝑛 ⃗ =0
𝑑𝑖𝑣 𝑉
𝜕𝑢 ⃗ ⃗ 𝜕𝑝
𝜌 + 𝜌 𝑉 ∙ 𝑔𝑟𝑎𝑑𝑢 = − + ∆ 𝑢
𝜕𝑡 𝜕𝑥
 suivant OX De le même
suivant OY manière on écrit
les deux autres
suivant OZ équations suivant
OY et OZ.
Ce sont les équations du mouvement de Navier-Stokes

Claude Louis Marie Henri Navier, né en

George Gabriel Stokes (13 août 1819 – 1er février
1785 et mort en1836 est un français. Il est
1903); Il publie, en 1845, le résultat de ses
spécialiste de mécanique newtonienne, il
travaux sur les mouvements des fluides. Son
établit en 1822 les
approche mathématique décrivant l'écoulement
équations décrivant le mouvement des fluides
d'un fluide newtonien incompressible dans un
, que nous connaissons sous une forme
espace tridimensionnel, en ajoutant une force de
élaborée avec le double nom de Navier-
viscosité à partir des équations d'Euler est à
Stokes, déterminante pour la
l'origine des équations de Navier-Stokes.
 Le régime d’écoulement
Dépend de la valeur d’un certain nombre adimensionnel (sans
dimension), c’est le Nombre de Reynolds, .
𝑙𝑎 𝑓𝑜𝑟𝑐𝑒 𝑑 ’ 𝑖𝑛𝑒𝑟𝑡𝑖𝑒 𝜌 𝑉𝐿
𝑅𝑒= ⇒ 𝑅𝑒=
𝑙𝑎 𝑓𝑜𝑟𝑐𝑒 𝑑𝑒 𝑓𝑟𝑜𝑡𝑡𝑒𝑚𝑒𝑛𝑡 𝑣𝑖𝑠𝑞𝑢𝑒𝑢𝑥 η
La masse volumique du fluide
: Viscosité du fluide
Vitesse moyenne d’écoulement (vitesse débitante)
Une longueur caractéristique de l’écoulement (exemple: : dans le cas d’un écoulement
dans une conduite cylindrique de diamètre D)

Stokes 𝑅𝑒 <1
Régimes
d’Ecoulement
Laminaire 𝑅𝑒≤2000
le nombre de
Reynolds,
Transitoire
𝑅𝑒≈ 4000
5
Turbulent 𝑅𝑒≥ 10
 Régimes d’Ecoulement et
le nombre de Reynolds,

Transitoire
Stokes Laminaire Turbulent

Osborne Reynolds
𝑹𝒆<𝟏 𝑹𝒆≤𝟐𝟎𝟎𝟎 𝑹𝒆≈𝟒𝟎𝟎𝟎 𝑹𝒆≥𝟏𝟎 𝟓 1842 -1912
(Angleterre)

Ingénieur et Physicien

Irlandais
hydrodynamique
dynamique des fluides

Introduction du
nombre de Reynolds
en 1883.
 Ecoulement de Stokes (écoulement rampant)

Lorsqu'un fluide visqueux s'écoule lentement en un lieu étroit ou

autour d'un petit objet, les effets visqueux dominent les effets inertiels
qui sont quasiment nuls.
L'équation de Stokes permet en particulier de décrire les écoulements
de liquide dans les systèmes microfluidiques.

Les écoulements de Couette

et de Poiseuille sont aussi décrits par cette équation.

 Ecoulement Laminaire

L’écoulement en régime laminaire (d'un fluide):

• se manifeste par ses différentes couches qui glissent les
unes sur les autres sans se mélanger.
• Les lignes de courant ne s’interpénètrent pas.
• Les effets visqueux dominent les effets inertiels.
• Il est régulier (il ne présente pas trop de variations spatiales ou
temporelles), bien souvent stationnaire.
Exemples: Fontaines laminaires lumineuses
 Application des équations de Navier-Stokes:
Ecoulement plan de Couette
Un écoulement de Couette désigne l'écoulement
Laminaire d’un fluide visqueux entre deux surfaces dont 𝒚
l'une est en mouvement par rapport à l'autre.

L'écoulement est dû à la contrainte de 𝒙

frottement visqueux de cisaillement qui agit sur le fluide.

Le nom de cet écoulement se réfère au physicien français

Maurice Couette, qui a inventé le premier viscosimètre𝒛 à
rotation.
Enoncé:
On considère l ’écoulement d’un fluide confiné
entre deux plaques infinies planes parallèles et
horizontales.
Le liquide est newtonien et incompressible .
Le champs de vitesse est orienté suivant .
L’écoulement est supposé stationnaire (pas de
1858 - 1943 dépendance au temps).
Application des équations du mouvement de L’équation de continuité
Navier-Stokes dans le cas de l’écoulement de dans le cas d’un fluide
Couette. incompressible:

/ OX

/ OY 𝑑𝑖𝑣 ⃗
𝑉 =0
/ OZ

Première étape: Détermination des équations simplifiées de Navier-Stokes:

=>

⃗ 𝝏𝒖 𝝏𝒗 𝝏𝒘 𝝏𝒖
𝒅𝒊𝒗 𝑽 = + + =𝟎 ⇒ =𝟎
𝝏𝒙 𝝏𝒚 𝝏𝒛 𝝏𝒙
Ceci veut dire que la vitesse dirigée suivant On dit aussi que la vitesse est
l’axe ne dépend pas de l’abscisse x. constante le long de la ligne de
On dit que l’écoulement est uniforme: courant
 , Ecoulement uniforme suivant la direction

𝒓
O 𝒙 𝒙 +𝒅𝒙

Mais,

Attention, ne pas confondre avec le profil uniforme dont l’écoulement peut ne pas être
uniforme

𝒓
O 𝒙 𝒙 +𝒅𝒙

Mais:
 Première étape: Simplification des équations de Navier-Stokes:

L’écoulement étant permanent / OX

(stationnaire):
/ OY

/ OZ

/ OX

/ OY

/ OZ Ecoulement Plaques infinies

uniforme (situation répétitive)

𝒚
𝝏𝒖 ⃗
𝒌
𝝏𝒛
𝒙

𝒛
 / OX

Produit scalaire = 0
/ OY

/ OZ Car

On remplace l’expression du Laplacien dans l’équation / :

OX: OX:
OY: OY:
OZ:  OZ:

OX:
OY:
On peut négliger la gravitation, n’étant pas
responsable de l’écoulement (contraintes
 OZ:
visqueuses)
 OX:
OY:
Les équations simplifiées de
Navier-Stokes OZ:

Deuxième étape: OX:

Détermination du OY:
profil de vitesse: OZ: 𝜕𝑝 𝜕𝑝
= =0 ⇒ 𝑷= 𝑷 ( 𝒙)
𝜕𝑦 𝜕𝑧
Dans cette équation,, on remarque que donc
Aussi, qui peut être égale à
Donc: cette équation impose la relation =.
On a égalité entre deux fonctions qui dépendent de deux variables indépendantes.
La seule solution possible à cette équation est: =.
la variation de la pression le long de l’écoulement est constante, .
On a aussi: qui par une première intégration donne:
Une deuxième intégration donne:
 Deuxième étape: 1 𝝏𝒑 𝟐
Détermination du profil 𝑢 ( 𝑦 )= 𝑦 + 𝐵𝑦 +𝐶
de vitesse: 𝟐 𝝏 𝒙
Dans cette équation représentant le profil de vitesse, nous avons deux constantes et à
déterminer
Pour cela on utilise les conditions aux limites qui traduisent la condition d’adhérence au
niveau des deux parois.

Exemple:
On considère la paroi supérieure animée d’un mouvement
avec une vitesse constante ,
La paroi inférieure est fixée dans un bâtit

On remplace dans l’expression de par les deux conditions: et

=> =>
Deuxième étape: 1 𝝏𝒑 𝟐
Détermination du profil 𝑢 ( 𝑦 )= 𝑦 +𝐵𝑦 +𝐶
de vitesse: 𝟐 𝝏 𝒙

On remplace dans l’équation de les expressions des deux constantes et .

On fait sortir en en facteur commun

Attention:

(( ) )
𝟐
1 𝝏𝒑 𝟐 𝒚 𝒚 𝑈 Il faut toujours vérifier la validité des
𝑢 ( 𝑦 )= 𝒉 − + 𝑦
𝟐 𝝏 𝒙 𝒉 𝒉 h deux conditions aux limites au niveau
de ce résultat du profil de vitesse.
Deuxième étape:

(( ) )
Détermination du 𝟐
1 𝝏𝒑 𝟐 𝒚 𝒚 𝑈
profil de vitesse: 𝑢 ( 𝑦 )= 𝒉 − + 𝑦
𝟐 𝝏 𝒙 𝒉 𝒉 h

Le profil de vitesse dépend de la valeur de la variation de la pression (=Cte)

Cas où (la pression est constante le Profil de vitesse

long de l’écoulement)
linéaire

La vitesse de déformation de cisaillement dans ce cas est constante sur n’importe quelle
couche fluide (couches horizontales).

La contrainte de cisaillement est également constante sur toutes les couches,

Loi de Newton: :
Deuxième étape:

(( ) )
Détermination du 𝟐
1 𝝏𝒑 𝟐 𝒚 𝒚 𝑈
profil de vitesse: 𝑢 ( 𝑦 )= 𝒉 − + 𝑦
𝟐 𝝏 𝒙 𝒉 𝒉 h
Le profil de vitesse dépend de la valeur de la variation de la
pression (=Cte)

On fait sortir en facteur commun de l’équation de la vitesse , La vitesse .

{ (( ) ) }
𝟐 𝟐
𝟏 𝝏𝒑 𝒉 𝒚 𝒚 𝑦
⇒ 𝑢 ( 𝑦 )=𝑼 − +
𝟐 𝝏𝒙 𝑼 𝒉 𝒉 h
Quelle est la dimension de ce terme ? Ce terme est adimensionnel

On va le nommer

{( ( )) }
EN fonction de la valeur que 𝟐
𝒚 𝒚 𝑦
prend P, on aura un profil de 𝑢 ( 𝑦 ) =𝑼 𝑃 − +
vitesse différent 𝒉 𝒉 h
 { ( ( )) }
EN fonction de la valeur que 𝟐
𝒚 𝑦
prend P, on aura un profil de 𝑢 ( 𝑦 ) =𝑼 𝑃 𝟏 − +
vitesse différent
𝒉 h
𝟐
𝟏 𝝏𝒑 𝒉 Couches plus rapides
P=−
𝟐 𝝏 𝒙 𝑼 que la paroi mobile!

3
2
P=
P=
- 1 0
= -3 -2 P= P=
P P=

1
P=
Couches se
déplaçant dans le
sens contraire de la
paroi mobile! 0

Pour , le mouvement du fluide créé par la plaque supérieure n'est pas assez fort pour
surmonter le gradient de pression positif défavorable, par conséquent on aura un reflux
(c'est-à-dire un négatif) qui se produit dans une région inférieure.

Pour , le mouvement du fluide créé par la plaque supérieure est amplifié par le gradient de
pression négatif favorable, donc est positif et peux même dépasser la vitesse de la paroi
au niveau de certaines couches.
 Deuxième étape: 1 𝝏𝒑 𝟐
Détermination du profil 𝑢 ( 𝑦 )= 𝑦 + 𝐵𝑦 +𝐶
de vitesse: 𝟐 𝝏 𝒙
Dans cette équation représentant le profil de vitesse, nous avons deux constantes et à
déterminer
Pour cela on utilise les conditions aux limites qui traduisent la condition d’adhérence au
niveau des deux parois.

Exemple:
On considère la paroi supérieure animée d’un mouvement
avec une vitesse constante ,
La paroi inférieure est fixée dans un bâtit

On remplace dans l’expression de par les deux conditions: et

et ,
Deuxième étape: 1 𝝏𝒑 𝟐
Détermination du profil 𝑢 ( 𝑦 )= 𝑦 +𝐵𝑦 +𝐶
de vitesse: 𝟐 𝝏 𝒙

On remplace dans l’équation de les expressions des deux constantes et .

On fait sortir en en facteur commun

Attention:

(( ) )
𝟐
1 𝝏𝒑 𝟐 𝒚 𝒚 𝑈 Il faut toujours vérifier la validité des
𝑢 ( 𝑦 )= 𝒉 − + 𝑦
𝟐 𝝏 𝒙 𝒉 𝒉 h deux conditions aux limites au niveau
de ce résultat du profil de vitesse.
Deuxième étape:

(( ) )
Détermination du 𝟐
1 𝝏𝒑 𝟐 𝒚 𝒚 𝑈
profil de vitesse: 𝑢 ( 𝑦 )= 𝒉 − + 𝑦
𝟐 𝝏 𝒙 𝒉 𝒉 h

Le profil de vitesse dépend de la valeur de la variation de la pression (=Cte)

Cas où (la pression est constante le Profil de vitesse

long de l’écoulement)
linéaire

La vitesse de déformation de cisaillement dans ce cas est constante sur n’importe quelle
couche fluide (couches horizontales).

La contrainte de cisaillement est également constante sur toutes les couches,

Loi de Newton: :