Chapitre 1

Notions fondamentales de la mécanique des Fluides

et la rhéologie
______________________________________________

Si le modèle de fluide newtonien décrit bien la très grande majorité des fluides composés de
molécules simples, il existe un bon nombre de fluides, dont certains sont d’usage très courant, qui
ont un comportement sous écoulement plus complexe. La définition d’un fluide newtonien est
assez restrictive : les contraintes de cisaillement sont proportionnelles au gradient de vitesse, ce
qui implique que :

– Dans un écoulement de cisaillement simple, les seules contraintes crées par l’écoulement sont
des contraintes de cisaillement.
– La viscosité est indépendante de la vitesse de cisaillement.
– La viscosité est indépendante du temps et les contraintes s’annulent immédiatement lorsque
l’écoulement est arrêté.
Toute déviation de ces règles est le signe d’un comportement non-newtonien. La description
des ces comportements et leur interprétation en relation avec la structure microscopique du fluide
constitue la discipline appelée rhéologie. Cette discipline est assez récente; elle a connu un
développement considérable avec l’apparition des polymères synthétiques.

1.1- Viscosité de cisaillement

Les expériences de Robert Hooke en 1678 ont mené à la loi qui porte son nom et qui constitue
une équation constitutive pour un solide élastique, pour de petites déformations. Cette loi exprime
que la force de traction est proportionnelle à la déformation, soit :

σ=Eε (1.1)

où σ est la contrainte de traction appliquée mesurée en Pascal (force par unité de surface), ε est
la déformation relative, définie comme étant l’allongement du solide divisée par sa longueur
initiale (donc sans dimension), et E le module d’Young, ou module élastique longitudinal,
mesurée en Pascal (abréviation : Pa), et qui est intrinsèque au solide.
En 1687, quelques années après les expériences de Hooke, Issac Newton exprime intui-
tivement dans ses Philosophiae Naturalis Principia Mathematica1 l’idée de viscosité pour un
liquide visqueux: la contrainte appliquée est proportionnelle au gradient de vitesse, et le
coefficient de proportionnalité est la viscosité2. L’équation correspondante s’écrit :

    (1.2)

où σ est la contrainte de cisaillement (mesurée en Pa dans le Système International (SI) de

mesure),  le gradient de vitesse dans l’épaisseur du fluide (mesuré en s-1), et η la viscosité
dynamique (mesurée en Pa.s). L’unité de la viscosité dynamique est parfois appelée Poiseuille (1
Poiseuille = 1 Pa.s = 1 kg.m-1.s-1). Dans le système CGS (centimètre, gramme, seconde), la
viscosité dynamique s’exprime en Poise (1 Poise = 0,1 Poiseuille = 0,1 Pa.s).

On définit également la viscosité cinématique υ qui est le rapport entre la viscosité dynamique
η et la masse volumique ρ du liquide. Elle est mesurée en m2 s-1 (dans le système CGS, l’unité de
mesure est le Stokes, St en abrégé, avec la correspondance suivante: 1 St = 100 cSt = 10- 4 m2 s-1 =
1 cm2 s-1). Il joue le rôle d’un coefficient de diffusion de la quantité de mouvement à travers les
différentes couches de fluides qui glissent imparfaitement les unes par rapport aux autres lors de
l’écoulement. Dans la suite de l’exposé, nous n’utiliserons que la viscosité dynamique η.
Les vérifications expérimentales n’arrivent que bien plus tard, avec les travaux de Poiseuille en
1856 et de Couette en 1890. L’équation de proportionnalité    est ainsi vérifiée pour de
nombreux liquides, que l’on qualifie alors de «newtonien». Historiquement, un liquide newtonien
est donc un liquide non élastique dont la viscosité ne varie pas suivant le gradient de vitesse.
Précisons la signification physique de l’équation    . Soit un liquide entre deux plans
parallèles (voir fig. 1.1). Le plan inférieur est fixe, et on exerce sur le plan supérieur une contrainte
σ : le plan supérieur se met alors à une vitesse v. On parle, pour une telle configuration, de
géométrie plan-plan ou de Couette plan.

Fig. 1.1a- Profil de vitesse pour un écoulement de cisaillement simple d’un liquide newtonien entre deux plans

En régime laminaire, le liquide peut être modélisé par une succession de couches infiniment minces
parallèles aux plaques et parallèles les unes aux autres. La déformation du fluide se fait par glissement
successif des couches les unes sur les autres par frottement, sans transport de matière d’une couche à l’autre
mais avec transport de la quantité de mouvement. Les conditions aux limites sont la vitesse nulle sur le plan
inférieur immobile, et la vitesse v sur le plan supérieur.
La déformation γ est le rapport entre la longueur parcourue l(t) par un élément de liquide du plan
supérieur et l’épaisseur e entre les deux plaques (épaisseur appelée entrefer) :

l (t )
 (1.3)
e
---------------------------------------------------------------------------------------------------------------------------------
1
Principes mathématiques de la philosophie naturelle.
2
Newton formule l’’hypothèse suivante [Newton, 1687] : « La résistance qui vient du défaut de lubricité des parties
d’un fluide doit être, toutes choses égales, proportionnelle à la vitesse avec laquelle les parties de ce fluide
peuvent être séparées les unes des autres ». La « résistance» est la viscosité, c’est-à-dire la réaction à la contrainte
appliquée, et la « vitesse avec laquelle les parties de ce fluide peuvent être séparées les unes des autres» est le gradient de
vitesse. Le « défaut de lubricité» des couches de fluides, c’est-à-dire le glissement imparfait des couches entres elles,
provient des échanges de quantité de mouvement entre les différentes couches du fluide. Le mot « lubricité» est à
prendre ici dans le sens de « lubrification »...
 La déformation γ dépend du temps t puisque le liquide se déforme continûment au cours du
temps : plus l(t) est grand pour e fixé, plus γ est grand; de même plus e est petit pour un même
l(t), plus γ est grand. Le liquide n’arrête cependant pas de se déformer avec le temps. À la place
de la déformation, on s’intéresse donc à une grandeur physique plus pertinente : la déformation par
unité de temps dγ/dt =  , que l’on appelle également taux de déformation ou encore taux de
cisaillement5. L’unité SI est le s-1. Le taux de cisaillement  n’est autre que la variation de la
vitesse du plan supérieur en fonction de l’épaisseur e: il s’agit donc également du gradient de
vitesse, qui s’écrit v/e. On a donc:

v
  (1.4)
e

Le concept est identique pour une couche infiniment mince qui compose le liquide, et on définit
donc un taux de cisaillement local  ( y )  v x ( y ) / y . La contrainte σ xy est la contrainte
appliquée par une couche de fluide sur la couche de liquide inférieure. En géométrie plan-plan,
σ xy est constante dans l’épaisseur e par équilibre mécanique entre les différentes couches de
fluide, et par conséquent, à l’état stationnaire, le gradient de vitesse  est également constant
dans l’épaisseur e. D’où la relation pour un liquide quelconque en géométrie plan-plan, qui définit
la viscosité de cisaillement η (  ) pour le liquide (constante dans l’épaisseur e en géométrie plan-
plan à l’état stationnaire) :

σ xy = η (  ) .  (1.5)
où encore plus explicitement :

 xy   ( ) v x / y (1.6)

Pour un liquide newtonien, la

viscosité η ne dépend pas du taux
de cisaillement  : la contrainte
σ xy est proportionnelle au taux de
cisaillement.
Fig. 1.1b – Contraintes engendrées par un écoulement de cisaillement..

1.2- Contraintes et déformations dans un fluide

1.2.1- Tenseur des contraintes

La relation σ xy = η (  )  relie la contrainte de cisaillement à la vitesse de déformation, ou taux

de cisaillement, avec comme coefficient de proportionnalité la viscosité. Pour un liquide non-
newtonien, il est nécessaire de tenir compte des contraintes perpendiculaires à l’écoulement,
plus communément appelées contraintes normales. Nous introduirons le tenseur des contraintes
à cet effet.
 Le tenseur des contraintes fait apparaître toutes les contraintes qui peuvent s’exercer sur un
matériau7, c’est-à-dire non seulement les contraintes tangentielles à l’origine des cisaillements,
mais également les contraintes à l’origine des compressions et des élongations du matériau.
Considérons un élément de surface d S placé à l’horizontale : la normale n à la surface d S
est suivant l’axe y (voir fig. 1.2a). Une contrainte σ appliquée à cette surface peut être
décomposée en trois composantes suivant les axes x, y et z: ce sont respectivement σ xy, σ yy et
σ zy. Plus généralement, dans la notation σij, i est la direction de la contrainte et j la direction de
la normale n à la surface dS. σ xy et σ zy sont des contraintes tangentielles et σ yy est une
contrainte normale.

Fig. 1.2– a) Composantes σ xy, σ yy et σ zy. de la contrainte s’exerçant sur une surface dS de normale orientée suivant
l’’axe y. b) Force de contrainte dF s’exerçant sur une surface dS dont la normale n est orientée de façon quelconque.

Pour une force dF de direction quelconque appliquée à la surface dS de normale n (voir fig.
1.2b), la contrainte σn appliquée à dS s’écrit :

σn = σ . n = dF/dS (1.7)

Soit, en notation matricielle,

  xn    xx  xy  xz   nx   dF 
    1  x
  yn     yx  yy  yz   n y    dFy  (1.8)
      dS  dF 
 zn   zx  zy  zz   nz   z

σ est le tenseur total des contraintes; il peut se décomposer sous la forme :

σ=τ–pδ (1.9)

p étant la pression hydrostatique, δ le delta de Kronecker (δ = δij avec δij = 1 si i = j, 0 si i ≠

j) et τ le tenseur des contraintes de viscosité qui apparaît lorsque le fluide est en mouvement.
Par convention, le signe négatif de la pression p traduit le fait que le liquide est en
compression : la pression est de sens opposé à la normale n de l’élément de surface dS
considéré (un signe positif traduirait une traction). Pour un fluide au repos, les seules
contraintes non nulles du tenseur des contraintes sont les termes diagonaux, et on a σ xx, σ yy et
σ zz = -p, le fluide étant uniquement soumis à la pression hydrostatique. Notons que pour les
contraintes de cisaillement, on a : σij = τi j , i≠j. Enfin, seules les composantes du tenseur total
des contraintes sont mesurables expérimentalement.
En rhéologie, le matériau étudié est souvent considéré comme étant incompressible. Or, pour
de tels fluides, on ne parvient pas à différencier la pression p des contraintes normales de
déformation τi j , lors de la mesure de σij = τi j - p [Bird e t a l . ( 1987), Morrison (2001)]. On
cherche alors à s’affranchir de la contribution de la pression isotrope p. Ceci est possible en
considérant les différences de contraintes normales, qui s’écrivent pour un écoulement de
direction x (y étant la direction normale aux surfaces entre lesquelles le fluide est cisaillé et z
une direction neutre, comme pour le cisaillement plan) :

N1 = σx x - σyy = τxx - τyy première différence de contraintes normales,

N2 = σyy - σzz = τyy – τzz deuxième différence de contraintes normales.

Pour l’eau, les différences de contraintes normales N1 et N2 sont nulles. L’absence de

différences de contraintes normales est donc prise comme une autre caractéristique des fluides
Newtoniens. Pour des fluides non Newtoniens présentant des différences de contraintes normales,
généralement, N1 est positif tandis que N2 est négatif. N1 est au moins une fois plus grand en ordre
de grandeur que N2 [Larson (1999), Bird et al. (1987)].

L’existence de différences de contraintes normales

dans certains liquides non Newtoniens est à l’origine
de phénomènes spectaculaires tels que l’effet
Weissenberg ou le gonflement de jets.
Considérons un liquide qui présente une différence
importante de contraintes normales lors d’un
écoulement de cisaillement (par exemple une solution
de polyacrylamide dans du glycérol [Bird et al.
(1987)]). Pour observer l’effet Weissenberg, on fait
tourner une tige cylindrique plongée dans ce liquide :
on constate alors que le liquide remonte le long de la
tige. Au contraire, dans un liquide Newtonien tel que
l’eau, le liquide se creuse autour de la tige par inertie Fig. 1.3- Mise en évidence de l’effet Weissenberg
[Barnes et al. (1989)].

1.2.2- Équation constitutive

L’étude d’un écoulement d’un fluide nécessite de déterminer la relation qu’il existe entre les
contraintes, explicitées par le tenseur τ, et les taux de déformation qui en résultent. Les
déformations sont données par les coefficients du tenseur G des gradients de vitesse :

vi
G  Gij  (1.10)
x j

Le tenseur G peut être décomposé en deux parties : une partie symétrique Dij, et une partie
antisymétrique ωi j [Guyon e t a l . ( 2001)]. On a alors :
Gij = Dij + ωi j (1.11)
 1  v v j  1  v v j 
avec Dij   i  et wij   i 
2  x j xi  2  x j xi 

La partie symétrique Dij est la somme de deux tenseurs: le tenseur diagonal tij,
représentant la variation de volume des éléments du fluide (variation nulle pour un fluide
incompressible), et le tenseur symétrique de trace nulle dij, traduisant les déformations des
éléments de fluide, sans variation de volume. La partie antisymétrique ωij traduit la rotation en
bloc des éléments du fluide (rotation dite pure ou solide). La rotation solide ne faisant pas
intervenir d’interactions entre éléments voisins, elle n’est pas liée aux contraintes appliquées
au fluide. On ne s’intéresse donc qu’à la partie symétrique de G: le tenseur des déformations
D= Di j , qui caractérise la déformation du fluide.
Une équation constitutive d’un fluide est une équation qui relie les tenseurs τ et D, c’est-
à-dire une relation qui relie les contraintes (en dehors des contraintes liées à la pression du
fluide au repos) et les déformations du fluide. L’équation constitutive pour un fluide
newtonien est [Larson (1999)] :

τ=2ηD (1.12)

Cette équation se réduit, pour le cisaillement plan de la fig. 1.1, à :

 xy   v x / y (1.13)

L’équation constitutive, ou loi de comportement, permet de caractériser l’écoulement d’un

fluide en indiquant quelle est la réponse d’un fluide à une contrainte imposée, ou quelle contrainte
doit être appliquée afin d’obtenir une déformation donnée. Pour un fluide réel, il est difficile
d’expliciter une équation constitutive pour un écoulement quelconque qui serait à la fois de
cisaillement et élongationnel, et pour lequel par ailleurs les coefficients du tenseur des contraintes
dépendraient du temps. C’est pourquoi on se contente de déterminer une équation constitutive
pour des écoulements spécifiques: de tels écoulements, plus simples à étudier, mettront en
évidence une ou plusieurs caractéristiques, non Newtoniennes ou non, du fluide.
L’équation constitutive est ainsi un des éléments permettant d’étudier un écoulement. Pour
résoudre complètement un problème d’écoulement, il est nécessaire de connaître :

 l’équation constitutive du fluide: τ = f (D), ce qui permet d’expliciter le tenseur des

contraintes :

σ=τ-pδ (1.14)
 l’équation du mouvement (pour un fluide compressible ou non) [Guyon e t a l .
( 2001)]:
 v 
  v grad v   div  F (1.15)
 t 
Or div   div (  p )  div   grad p (1.16)
 v 
  v grad v    grad p  div   F (1.17)
 t 

ρ étant la masse volumique du fluide et F les forces extérieures (gravitation, etc.)

  l’équation de conservation de masse (ou de continuité) :


 div ( v )  0 (1.18)
t
Soit div v = 0 pour un fluide incompressible;
 les conditions aux limites.

Dans le cas d’un fluide newtonien, div τ devient [Guyon et al. (2001), Morrison (2001)] :


div    v  (   ) grad (div v)
3 (1.19)

où ς est la « viscosité de volume », dont les contraintes associées (termes diagonaux du tenseur
τ) correspondent aux variations de volume par compression. Si de plus le fluide est
incompressible (div v = 0), on obtient : div τ = ηΔv. L’équation du mouvement pour un fluide
newtonien incompressible, appelée équation de Navier-Stokes, est donc:

v
(  v.grad v )   grad p   v  F (1.20)
t

Lorsque le nombre de Reynolds Re est petit, on peut négliger le terme non linéaire. En se
plaçant par ailleurs en régime permanent, l’équation de Navier-Stokes se réduit à l’équation de
Stokes:

grad p   v  F (1.21)

1.2.3- Viscosité élongationnelle d’un liquide newtonien

Dans le paragraphe précèdent, nous avons défini la viscosité de cisaillement η résultante des
contraintes de cisaillement. Nous montrons à titre d’information un autre type de défor-
mation: l’élongation, qui nous amène à définir la viscosité élongationnelle. Un écoulement
élongationnel est un écoulement dans lequel les éléments de fluide sont étirés sans rotation ni
cisaillement. La fig. 1.4 montre l’élongation simple d’un matériau dans la direction x.

Fig. 1.4– Élongation simple dans la direction x.

 Pour un tel écoulement, l’élongation ε(t) dans la direction x est la déformation relative à
l’instant t dans cette direction. L’élongation ε(t) et le taux d’élongation έ s’écrivent:

x(t  dt )  x(t ) dx 1  dx 
 (t )   et    
x(t ) x(t ) dt  x(t )  (1.22)

d’où, en posant la vitesse vx=dx/dt , l’obtention de la relation: vx = έ x. Que valent vy et vz ?

Par symétrie, vy = vz. En tenant compte de l’équation de conservation pour un fluide
incompressible (div v = 0), on obtient finalement :

1 1
v x   x , v y    y et v z    z
2 2

Le tenseur des déformations s’écrit donc :

 
 
  0 0 
1  v v  1 (1.23)
D  Dij   i  j    0   0 
 
2  x j xi   2 
 1 
0 0   
 2 

Le tenseur total des contraintes τ s’écrit donc, pour un fluide newtonien (pour lequel l’équation
constitutive est τ = 2 η D) :
 2 0 0 
 
   0   0  (1.24)
 0 0   


On définit la viscosité élongationnelle comme étant :

 xx   zz
e  (1.25)


Pour un fluide newtonien, d’après les coefficients τij du tenseur des contraintes de
déformation τ, la viscosité élongationnelle ηe vaut:

ηe = 3 η (1.26)

La viscosité élongationnelle ηe est donc une constante pour un fluide newtonien. Le

coefficient de proportionnalité est appelé rapport de Trouton Tr: Tr ≡ ηe /η, et vaut donc 3 pour
un fluide newtonien. Cela signifie en particulier que, pour un fluide newtonien, la
connaissance de la viscosité de cisaillement η permet de déterminer la viscosité
élongationnelle ηe. La viscosité élongationnelle η (έ) d’un fluide non-newtonien dépend, elle, du
taux d’élongation έ.
1.2.4- Paramètre influençant la viscosité la viscosité
a- Influence de la température sur la viscosité
On observe que la viscosité des gaz croit avec la température. On pourra interpoler cet effet
à partir de la relation de Arnold [Arnold (1933)]:
C1T 3 / 2
 (1.27)
T  1,47 Teb

où C1 est une constante pour un gaz donné et Teb sa température d’ébullition

Par contre la viscosité des liquides newtoniens décroît avec la température. On pourra
utiliser entre la température de fusion et le point d’ébullition, la formule de Guzman-
Andrade [Andrade (1930)] :
 C2 
  C1 exp  (1.28)
T 
où C1 et C2 sont des constantes.

b- Influence de la pression sur la viscosité

L’influence de la pression sur la viscosité dynamique des gaz reste faible jusqu’à 20 bars
environ. Au-delà, la technique des états correspondants permet une bonne évaluation de η en
fonction de P.
On admet fréquemment pour les liquides, la relation empirique :
 P 
 ( p) 
P
1 
 a 0 
(1.29)
 ( p0 )
1.2.5- Valeurs de la viscosité
a- Viscosité dans le cas des mélanges
La viscosité d’un mélange n’est pas fournie par la règle d’additivité qui ne s’applique même
pas aux mélanges de gaz permanents.
- Pour les gaz, Wilke [Wilke (1950)] propose :
2
   1 / 2  M 1 / 4 
1   i   j  
    M  
n
i  j  i 
m   avec ij   (1.30)
n  yj   M 
1/ 2


i 1
 ij  2 2 1  i 

j 1  yi
  M 
  j 

Où ηi , yi et Mi sont la viscosité, le titre et la masse molaire du iè composant.

- Pour les mélanges binaires de liquide si  2  1  0,015 Pa.s , on pourra utiliser la
corrélation de Kendall-Monroe [Kendall et Monroe(1921)] :

 m1 / 3  x1 11 / 3  x2  21 / 3 (1.31)

Où xi est le titre molaire de l’espèce de viscosité ηi.

b- Viscosités dans le cas des solutions ou dispersions

Depuis la relation de Stokes-Einstein, qui donne la viscosité d’une suspension de faible

fraction volumique en fonction de la phase continue, et de la fraction volumique [Einstein (1906)],
de nombreux travaux ont tenté d’établir des relations empiriques ou des équations
phénoménologiques pour modéliser la viscosité d’un milieu concentré.
En 1906, Einstein donne une expression de la viscosité en fonction de la concentration pour
des suspensions diluées de sphères dures (φ < 2%). Il ne prend pas en compte les effets
d’interactions entre particules. La relation qu’il propose est de la forme :


   1  K 1
r
s (I.32)

Vs
φ étant la fraction volumique du solide définie par :  
V  Vs

Où K1 = 2.5, ηr est la viscosité relative de la suspension, η est la viscosité de la suspension, ηs est

la viscosité de la phase suspendante (solvant), Vs volume de solide et V volume de solvant
(liquide). Des valeurs de K1 ont aussi été calculées par différents auteurs et pour différentes formes
de particules. Généralement K1 > 2,5 pour des particules allongées et K1 < 2,5 pour des particules
déformables.

Pour étendre la validité de l’équation d’Einstein à des concentrations plus élevées, il est
nécessaire de prendre en compte les effets d’interactions entre les sphères. Cela introduit dans
l’équation d’Einstein des termes complémentaires en φ2.
En 1942, De Bruijn (1942) donne une expression de la viscosité relative en prenant compte la
perturbation hydrodynamique due à la présence d’autres particules de la forme :

 r  1  2,5  1,552 2 (I.33)

Ainsi, lorsqu’on introduit un soluté dans un solvant, les interactions moléculaires

augmentent et la viscosité également.
Afin de différencier les influences du soluté (ou de la phase dispersée) et du solvant, on
définit quelques termes de viscosité:

- La viscosité relative : ηr = η / ηsolvant

C’est le rapport des viscosités dynamiques de la solution et du solvant.

 - La viscosité spécifique : ηsp = ( η - ηsolvant ) / ηsolvant = ηr – 1
C’est l’augmentation relative de la viscosité à cause du soluté, elle dépend de la
concentration du soluté.

- La viscosité réduite : ηred = ηsp /C où C est la concentration en soluté

Unité usuelle : C en g/cm3 donc ηred en g-1.cm3.

- La viscosité intrinsèque : [red ]  lim
red
c 0

C’est une caractéristique de la nature de la solution

Intérêt :
 Pour les suspensions et les émulsions, la viscosité intrinsèque peut nous renseigner
sur la forme des particules. Par exemple, pour des particules sphériques, on trouve
[ηRed] = 2,5/ ( masse volumique), la relation est autre pour une suspension de
particules ellipsoïdales ou en bâtonnets.
 pour les solutions de macromolécules, la viscosité intrinsèque permet de déterminer
la masse molaire moyenne selon la loi de Mark Houwink :

[ red ]  K M 

Où K est une constante pour un couple polymère /solvant donné.