Ecoulement de Stokes 2020/2021

ECOULEMENT A TRES FAIBLE NOMBRE DE REYNOLDS

ECOULEMENT DE STOKES

Introduction
Le rapport entre forces d'inertie et de viscosité peut varier d'un écoulement à un autre et, pour
un écoulement donné, d'un point à l'autre du champ. La prépondérance relative de l'une ou
l'autre de ces forces dépend de la valeur globale ou locale du nombre de Reynolds. Dans le
cas où ce paramètre est globalement petit, l'écoulement est de type rampant et relève d'un
modèle spécifique. Lorsque le comportement visqueux du fluide est de type Newtonien, ce
modèle correspond à une forme simplifiée linéaire du modèle de Navier-Stokes, connue sous
le nom de modèle de Stokes. Objectif de cette partie.

Hypothèses
Fluide visqueux Newtonien
1. le milieu est incompressible ;
2. les forces d'inertie sont en tout point et à tout instant, négligeables vis-à-vis des forces de
viscosité.
U0 L
Re  1

Uo, est une échelle représentative de la vitesse du mouvement du fluide,
L une échelle de longueur
ν la viscosité cinématique du milieu

Le modèle gouvernant la dynamique de l'écoulement se réduit aux deux équations suivantes :

  C ste  divV  0

V 1
 F  gradP   V
t 

Pr E.G. FILALI Page 1

 Ecoulement de Stokes 2020/2021
L'équation de la dynamique ainsi simplifiée fait intervenir quatre échelles de temps
caractérisant respectivement l'instationnarité (Ti) ainsi que les forces extérieures de volume
(TF), de pression (TP) et de viscosité, (Tν).

L2
Rapportée à l'échelle de longueur L, cette dernière peut être estimée par : Tv 

Nous supposons en plus que les conditions suivantes sont vérifiées :

1 . Ti >> Tν, de sorte que l’écoulement peut être considéré comme quasi-stationnaire ;
2. TP ~ TF ~ Tν, ce qui signifie que les trois forces extérieures sont du même ordre de

grandeur
3. les forces extérieures de volume se réduisent aux seules forces de gravité conduisant à un
champ de pression hydrostatique Po(r) = pgz en l’absence de mouvement, conformément à
l'équation :
1
F grad  P0   g

Equation de Stokes
Si le nombre de Re est très faible (effets important de la viscosité) le terme inertiel peut de
plus être négligé même si l’écoulement n’est pas du tout parallèle.
1
0g  grad  P    V divV  0


Equations en formulation vitesse-pression

Les équations qui régissent le champ de vitesse et de pression :
grad  P  P 0   V

L’équation de Stokes étant linéaire il y a unicité de la solution une fois connues les conditions
aux limites.

Equations en formulation pression-rotationnel

Il existe d’autres formes équivalentes de l’équation de Stokes :
 Première variante :

 
grad  P  P 0     rot 

Pr E.G. FILALI Page 2

 Ecoulement de Stokes 2020/2021


En effet V   divV  rot rotV    mais divV  0 pour un fluide incompressible et la

vorticité est définie par   rot V  

 Deuxième variante en prenant la divergence de l’équation précédente :

2  P  P0   0 ou encore   P  P0   0

En effet on a 
  P  P0   div   P  P0   
et div rot    0
La pression satisfait donc l’équation de Laplace à très faible nombre de Reynolds

 Troisième variante :

 
2
   0 ou encore   0

En effet,
2
 
   
    div rotV   rot rot      1 rot   P  P   0
0

Ceci montre qu’à très faible nombre de Reynolds, il n’y a plus transport de la vorticité. La
vorticité est dans un état d’équilibre et ne diffuse plus.

Validité des écoulements rampants

1. déplacements à très faible vitesse ; déplacements de géomatériaux. Malgré
l’importance des échelles de longueur mises en jeu, la lenteur des déplacements est
telle que le nombre de Reynolds global reste très faible. A titre d'exemple, on peut
citer le déplacement des glaciers ou celui de l'écorce terrestre.

2. mouvements à très petite échelle ; les mouvements à échelle "microscopique". Il peut

s'agir d'écoulements internes dans des géométries de taille réduite telles que les canaux
et cavités d'un milieu poreux. Il peut s'agir encore du déplacement d'objets de très
faibles dimensions tels que bactéries ou suspensions diverses.

3. écoulements de fluides très visqueux : écoulements de fluides tels que certaines

peintures, les huiles, le bitume, le verre fondu ou encore le miel.

Pr E.G. FILALI Page 3

 Ecoulement de Stokes 2020/2021
Ordres de grandeurs de configurations pratiques.
 En tectonique, la vitesse de déplacement des plaques est de l'ordre de quelques
-8
centimètres par an, soit pour une amplitude de 30 cm, une vitesse d'environ 10
m/s. Avec une viscosité  ≈ 1020 m2/s cela conduit à un nombre de Reynolds de
-23
l'ordre de l0 pour une échelle de longueur de 100 km.
 Pour un objet dont la taille est de l'ordre d’un micromètre est qui se déplace dans
de l'eau  ≈ 10-6 m2/s à une vitesse de 10 μm/s, le nombre de Reynolds
-5
correspondant est de l'ordre de 10 .
 Enfin, pour du verre fondu  ≈ 10-2 m2/s soufflé par un orifice de 1 cm de
-2
diamètre à une vitesse de 1 cm/s, le nombre de Reynolds est de l'ordre de 10 .

Ecoulement rampant de Hele-Shaw

En écoulement de Hele Shaw (1898), un fluide très visqueux se déplace à faible vitesse Ua
entre deux plans parallèles distants d'une faible épaisseur (2h). Si l'on dispose
transversalement un obstacle de dimension caractéristique L, nous allons voir que sous
l'hypothèse de régime rampant, cette configuration permet de simuler, à distance de
l'obstacle, le champ de l'écoulement à potentiel de fluide parfait incompressible autour de
celui-ci.

Pr E.G. FILALI Page 4

 Ecoulement de Stokes 2020/2021
h Ua L
On pose :   et Re 
L 
On néglige les forces de gravité (pesanteur), et on suppose que les conditions suivantes sont
vérifiées :  << 1 et Re << 1
On désigne par U, V, W et P les composantes de la vitesse et la pression en tout point (x,y,z)
de l’écoulement dans la cellule.

1) estimation des ordres de grandeurs :

U V W
Equation de continuité   0
x y z

x  L

Ordres de grandeurs y  L avec  =h/L << 1
z  h


U  U a

V  U a

W  U a
'

U 0 U 0 U 0' U 0 U 0 U 0'
   0 ou encore   0
L L h L L L

U V W
Les termes et étant du même ordre de grandeur que  U a'  U a
x y z

On peut répéter cette analyse dimensionnelle pour les équations de Navier-Stokes (en
particulier les termes d’inertie et de viscosité), il vient en projection sur x, y et z :

Pr E.G. FILALI Page 5

 Ecoulement de Stokes 2020/2021

• les contributions du Laplacien sont telles que :

2 2 2

x 2 y 2 z 2

• les termes d’inertie et de viscosité vérifient :

 Inertie   2 Re Vis cos ité 

• les termes de viscosité vérifient :

Vis cos ité x  Vis cos ité y et Vis cos ité x   2 Vis cos ité z

Ne retenant que les termes d’ordre supérieurs, les gradients de pression sont équilibrés par les
forces de viscosité en z :
 P 1  2U
0   x   z 2

 P 1  2V
0   
 y  z 2
 P
0  
 z

La pression est donc indépendante de z, et P(x,y,z) = P(x,y).

Pr E.G. FILALI Page 6

 Ecoulement de Stokes 2020/2021

L’intégration par rapport à z donne :

 1  z2  P
U  x, y, z      z   
  2  x

V x, y, z  1  z   ' z   '  P
2

    
  2  y
Les constantes d’intégration sont déterminées par les conditions d’adhérence sur les plans
z   h où U  x, y  h   0 et V  x, y  h   0

1 2 2 P 1 2 2 P
U  x, y, z  
2
 z h 
x
et V  x, y, z  
2
 z  h  y
Dans le plan médian z = 0
h2 P h2 P
U  x, y , 0   U 0  x, y    et V  x, y, 0   V0  x, y   
2 x 2 y

Et ce vecteur vitesse médian V0 ( x, y ) , vérifie les deux équations :

divV0  x, y   0 P  0
 
 h 2  P P  soit  h2
 0       0     grad  P  x, y  
2  x y 
V x , y i j V x , y
  2

P  0
 écoulement de Stokes
 

 
h2
2
 
rot V0  x, y     rot grad  P  x, y    0 écoulement à potentiel des vitesses

rot (grad()) = 0

Ces relations indiquent que la pression est solution de l’équation de Laplace et le champ de
vitesse est irrotationnel, comme l'est celui de fluide à potentiel. quelque soit la forme de
l’obstacle.
Les lignes de courant en cellule Hele-Shaw ( Re ~ 0) fournissent donc une analogie du champ
de l'écoulement à potentiel de fluide parfait (Re→∞).

APPROXIMATION DE LA LUBRIFICATION HYDRODYNAMIQUE

C’est une approximation très importante et est abondamment utilisée pour ramener un
problème d’écoulement 2D ou 3D à un problème 1D. Ecoulements quasiparallèles quasi-
stationnaires
Qu’appelle-t-on lubrification

Pr E.G. FILALI Page 7

 Ecoulement de Stokes 2020/2021
- Ecoulement entre surfaces proches et presque parallèles ou de films fluides minces à
surfaces libres.
- Applications : réduire les frottements entre pièces proches

Une des circonstances importantes dans laquelle l’inertie du fluide peut être négligée
concerne les écoulements dits de lubrification. Ce terme recouvre les écoulements de
fluides visqueux confinés entre deux parois solides très proches en mouvement relatif. Les
deux parois délimitent un espace de très grand rapport d’aspect : l’épaisseur moyenne < h >
est très petite devant la longueur L.
Si l’épaisseur h varie d’un point à un autre de l’écoulement, on se limite au cas où cette
variation est très lente, c’est-`a-dire : dh/dx très inférieur à 1.

Cette géométrie particulière a deux conséquences :

 la composante longitudinale de vitesse u est beaucoup plus grande que la
composante transverse v.
u v u v L
    0 avec u  v v
x y L h h

 La force qui s’exerce sur les parois est beaucoup plus grande dans la direction
normale à l’écoulement que dans la direction parallèle à l’écoulement.

Pr E.G. FILALI Page 8

 Ecoulement de Stokes 2020/2021

La force normale Fn est l’intégrale de la pression sur une des parois solides :

Fn    P  P0  dx   P  P0  L
L

Toujours en tenant compte du rapport d’aspect L/h >> 1, à partir de l’équation de Stokes :
  2u  2u   2u U P  P0
  2 
   2 
 x y  y
2 2
h L

P  P0 UL UL2
  2 et Fn   2
L h h
La force Ft dans la direction parallèle à l’écoulement est l’intégrale de la contrainte
tangentielle sur la paroi solide :
L L u UL
Ft    xx dx    dx  
0 0 y h
h
Ft  Fn 1
L

Ceci permet de supporter des charges importantes tout en conservant la possibilité d’un
mouvement relatif avec une résistance faible.
C’est le principe utilisé dans la lubrification des pièces mécaniques en rotation, dans les
articulations entre les os, dans la sustentation des têtes de lecture de disques magnétiques ...

Pr E.G. FILALI Page 9