PC - Lycée Dumont D’Urville

TD 3 mécanique des fluides

I. Ecoulement de Poiseuille

En 1835 un médecin français, Poiseuille fit une série

d’expériences pour déterminer comment un fluide
visqueux s’écoule dans un tuyau droit. Son but était
de comprendre la dynamique de la circulation san-
guine chez l’homme sachant que le plasma sanguin
se comporte comme un fluide newtonien. Le sang
s’écoule dans le corps humain sous l’action du coeur
(pompe) et d’un gradient de pression ( pression de
l’ordre de 27 hP a aux pieds et pression de l’ordre de
8 hP a à la tête).

On considère l’écoulement d’un fluide visqueux dans

un tube cylindrique de rayon R, parallèlement à l’axe R
r
V(r)

Oz du tube, de sorte que la vitesse est de la forme

→
−
V = V (r)− →
ez . L’écoulement est du à un gradient z=0 z=L Oz

de pression ente les deux extrémités du tube dis-

tants de L, de sorte que la pression est de la forme
P (z) = P0 − Kz avec K constante positive. On
néglige les forces de pesanteur.

→ 1 d dV →
−
1. Ecrire l’équation de Navier-Stokes, la simplifier et la projeter sur Oz. On donne ∆ V = (r )−
ez .
r dr dr
→
−
2. En déduire le champ des vitesses V .
3. Calculer le débit volumique à travers une section du tube.
4. Représenter sur une section du tube le vecteur vitesse et en déduire le sens de la force exercée par le
fluide sur le tube.
K 2 KR4 π
Réponses: V (r) = (R − r2 ), Dv = .
4η 8η
II. Ecoulement de Couette

Oy
On considère l’écoulement unidimensionnel perma-
nent d’un liquide de masse volumique ρ et de vis-
cosité η entre deux plans parallèles de surface S L
très grande et distants de H. Le fluide sera sup- fluide
H
posé soumis aux seules forces de pression et viscosité. O visqueux Ox

L’écoulement, appelé écoulement de Couette, est de

→
−
la forme V = V (x, y)− →
ex . On suppose la pression au
sein du fluide uniforme et égale à P0 . L’écoulement Oz
est du au déplacement du plan supérieur à la vitesse
→
−
V = V0 − →
ex , le plan inférieur restant fixe.

1. Que peut-on déduire de l’équation locale de conservation de la masse sur l’expression V (x, y)?
2. Déduire de l’équation de Navier-stokes, l’expression de V (y).
3. Représenter le profil des vitesses sur un section perpendiculaire à l’écoulement et calculer le débit volu-
mique à travers cette section de hauteur H et de largeur L.

→
− dV −
→
4. Calculer la force exercée par le fluide sur chacun des deux plans. On donne F = ±ηS ex . Justifier le
dy
→
−
signe dans l’expression de F par un raisonnement physique.
→ V0 y −
− → V0 HL2
Réponses : 1- V ne dépend pas de x 2- V = ex 3- Dv =
H 2

1
III. Ecoulement de Poiseuille entre deux plans

Oy
On considère l’écoulement unidimensionnel perma-
nent d’un liquide de masse volumique ρ et de viscosité
η entre deux plans parallèles de surface S très grande L
et distants de H. Le fluide sera supposé soumis aux fluide
H
seules forces de pression et viscosité. L’écoulement O visqueux Ox

→
−
est de la forme V = V (y)− →
ex .

Oz

1. Ce modèle d’écoulement est-il compatible avec un écoulement incompressible?

On suppose ici que les deux plans sont fixes et qu’ il existe un gradient de pression uniforme au sein du fluide,
−−→
de la forme gradP = −K − →
ex où K est une constante positive. L’écoulement, du à ce gradient de pression est
appelé écoulement de Poiseuille.
2. On donne P (x = 0) = P0 . Exprimer P (x).
3. Déduire de l’équation de Navier-stokes, l’expression de V (y).
4. Représenter le profil des vitesses sur un section perpendiculaire à l’écoulement et calculer le débit volu-
mique à travers cette section de hauteur H et de largeur L.

→
− −
→ K KH 3 L
Réponses : 1- div( V ) = 0 2- P (x) = P0 − Kx 3- V = − (y 2 − Hy)−
→
ex 4- Dv =
2η 12η

IV. Ecoulement sur un plan incliné

Oy Un liquide de viscosité η, de masse volumique ρ

s’écoule sur un plan incliné faisant un angle α par rap-
port à l’horizontale sur une hauteur e. L’écoulement
P0 g est permanent et incompressible. La vitesse s’écrit
e →
−
V = V (x, y)− →
ex et la pression ne dépend que de y, on
fluide
la note P (y).
α

Ox

1. Que peut-on déduire de l’équation locale de la conservation de la masse? En déduire que l’accélération
des particules fluide est nulle.
2. Projeter sur Ox et Oy l’équation de Navier-Stokes.
3. L’atmosphère impose la pression P0 uniforme à l’interface avec le fluide. Montrer que P (y) = P0 −
ρg cos α(y − e).
4. En déduire V (y) en fonction de deux constantes d’intégration. Que vaut V (y = 0)? justifier votre
dV
réponse. On a (y = e) = 0 par continuité de la force de viscosité en y = e. En déduire les deux
dy
constantes d’intégration.
5. Calculer le débit volumique à travers une section de largeur L selon Oz.

ρg sin α g sin αe3 L

Réponses: 4- V (y) = 2η (2ey − y 2 ) 5- Dv = 3ν

2
V. Equation de propagation du son

On étudie la propagation du son selon Ox dans un fluide supposé parfait. Le mouvement du fluide est
décrit sans prendre en compte les effets de viscosité et les échanges thermiques à l’intérieur du fluide. Les
détentes et compressions locales du fluide sont isentropiques, le coefficient de compressibilité isentropique
−1 ∂V
s’écrit: χS = |S . La propagation des ondes sonores est associée à un écoulement irrotationnel. Les
V ∂P
effets de pesanteur ne sont pas pris en compte.
On pose:
→
−
v (x, t) = v(x, t)−
→
ex : vitesse des particules en x à t
P (x, t) = P0 + p(x, t) : pression en x à t
ρ(x, t) = ρ0 + µ(x, t) : masse volumique en x à t
−−−−→
1. On suppose que v(x, t), p(x, t) et µ(x, t) sont des infiniment petits d’ordre 1. Quel nom porte cette
approximation?
1 ∂ρ
2. Expliquer le signe − dans l’expression de χS donnée dans l’énoncé. Montrer que χS = |S et
ρ ∂P
exprimer χS en fonction de p, ρ0 et µ.
3. Ecrire l’équation de conservation de la masse et la simplifier en ne gardant que les termes d’ordre 1.
4. Ecrire l’équation de Navier-Stokes et la simplifier en ne gardant que les termes d’ordre 1.
5. En déduire l’équation de propagation vérifiée par p(x, t). En déduire la vitesse de propagation C du son.
6. Le fluide est de l’air assimilé à un gaz parfait à la température. Après avoir déterminé en fonction du
rapport γ des capacités thermiques molaires du gaz et de la pression P0 , établir l’expression de C en fonction
de T0 , du rapport γ, de la masse molaire M du fluide et de la constante R des gaz parfaits. AN pour l’air à
T0 = 293 K avec γ = 1, 4 et M = 29 g/mol.

∂v ∂µ ∂v ∂p ∂2p ∂2p
Réponses : ρ0 = − , ρ0 =− et µ = ρ0 χS p, équation de d’Alembert : 2
= ρ0 χ S 2
∂x ∂t ∂t ∂x ∂x ∂t

VI. Forces de viscosité sur une paroi

Pe Ps<Pe
On considère un écoulement d’un fluide de viscosité
η dans une canalisation de rayon R, d’axe Oz et de
V(r)
longueur L, sous l’action d’un gradient de pression r
−−→
uniforme gradP = K − →
ez . On note P (z = 0) = Pe
z=0 z=L Oz
et P (z = L) = Ps . Le champ des vitesses s’exprime
→
− K 2
V = (R − r2 )−
→
ez .
4η

1. Exprimer K.
2. Représenter le profil des vitesses sur une section de la canalisation et en déduire qualitativement le sens
et la direction des forces de viscosité sur les parois.
→
− dV
3. Dans ce type de géométrie, la force surfacique de viscosité s’écrit fs = −η (r)−
→
ez . En déduire la
dr
force totale exercée par le fluide sur le tube qui le contient. Vérifier le sens de la force et l’homogénéité de
l’expression.

→
− πR2 −
→
Réponse: F = KL ez
2

3
VII. Temps de convection et temps de diffusion
−
→
∂V
Dans l’équation de Navier-Stokes, lorsqu’on ne conserve que le terme convectif, on obtient: ρ =
− −−→ −
→ → ∂t
−ρ( V .grad) V .
→
−
∂V →
−
Dans l’équation de Navier-Stokes, lorsqu’on ne conserve que le terme diffusif, on obtient: ρ = η∆ V .
∂t
Déduire de ces deux équations, une expression donnant l’ordre de grandeur du temps caractéristique de
convection τconv et du temps caractéristique de diffusion τdif f .
Exprimer le nombre de Reynolds en fonction de ces deux temps et commenter les cas particuliers Re << 1
et Re >> 1.

VIII. Force de trainée

On considère une sphère de rayon R, de vitesse V , en mouvement uniforme dans un fluide de viscosité η et
de masse volumique ρ.
π
1. La force de trainée exercée par le fluide sur la sphère a pour expression F = Cx (Re )Rα V γ ρλ où Cx
2
est un coefficient sans dimension et α, γ et λ sont des entiers naturels. Déterminer α, γ et λ par une analyse
dimensionnelle.
2. Pour un écoulement laminaire (Re < 1), la loi de Stokes donne l’expression de la force de trainée soit
→
− →
−
F = −6πRη V . Exprimer dans ce cas, Cx en fonction de Re .
3. Que devient cette force pour un fluide parfait?

24
Réponses : Cx =
Re

IX. Viscosimètre à chute de bille

Une bille solide sphérique de rayon R et de masse volumique ρs tombe dans un fluide visqueux de masse
→
− →
−
volumique ρl et de viscosité η. On suppose que la force de traı̂née s’écrit −6πηR V où V est la vitesse de la
bille par rapport au fluide.
1. Exprimer la vitesse limite U0 atteinte très rapidement par la bille. En déduire η en fonction des données
et de U0 .
2. AN : calculer η et ν pour ρs = 300 kg/m3 , ρl = 126 kg/m3 , U0 = 8 mm/s et R = 4, 5 mm.
3. Calculer le nombre de Reynolds et justifier a postériori l’expression de la force de traı̂née.

4
X. Correction TD 3 exo IX: Viscosimètre à bille

La bille, dans le référentiel terrestre supposé galiléen, est soumise à:

4
- son poids ρs πR3 − →g
3
4
- la poussée d’Archimède −ρl πR3 − →g
3
→
−
- la force de trainée −6πRη V
Au sujet de la poussée d’Archimède : c’est le nom que l’on donne à la résultante des forces de pression
appliquées sur un solide dans un fluide. C’est une force verticale ascendante (donc de meme direction que
le poids mais de sens opposé). La norme de cette force est égale au poids du volume de fluide déplacé : on
la trouve en faisant le produit : masse volumique du fluide fois volume du solide immergé dans le fluide fois
4
g. Ici le volume de fluide déplacé est le volume de la sphère πR3 donc la poussée d’Archimède s’écrit :
3
4
−ρl πR3 − →g (le moins signifie que la force est vers le haut, contrairement au poids).
3
Lorsqu’on plonge la bille dans l’eau, elle se met en mouvement et accélère sous l’action de son poids. Sa
vitesse augmente et donc la force de trainee augmente. Au bout d’un certain temps, la somme des trois
forces est nulle et le solide a un mouvement rectiligne uniforme, sa vitesse est constante et elle est appelée
vitesse limite.

en régime transitoire, en régime permanent,

la force de trainee est de plus la somme des trois
en plus grande car la vitesse de forces est nulle, la bille
la bille augmente a une vitesse constante
appelée vitesse limite
Farchi
Farchi

Ft
Ft

P
P

On écrit donc la RFD dans le régime permanent où la vitesse de la bille est constante:
→
− 4 4 →
−
0 = ρs πR3 − →
g − ρl πR3 − →
g − 6πRη Vl
3 3
→
−
On en déduit la vitesse limite Vl :
→ 2(ρs − ρl )R2 −
− →
V = g.
9η
Vérification de l’unité: η est en P a.s = N.m−2 .s = kg.m.s−2 .m−2 .s = kg.m−1 .s−1 , on vérifie alors facilement
que Vl est en m.s−1 .
2(ρs − ρl )R2 g 2(300 − 126).(4, 5.10−3)2 .10
AN : η = = = 0, 98 P l.
9U0 9.8.10−3
ρLV 126.4, 5.10−3.8.10−3
On calcule le nombre de Reynolds : Re = = = 5.10−3 << 1 : à faible nombre de
η 0, 98
Reynolds, on sait que la force de trainée est bien proportionnelle à V (remarque : à fort nombre de Reynolds,
on rappelle que la trainée est proportionnelle à V 2 , ce n’est pas le cas ici). Donc le modèle de la trainée
linéaire donnée dans l’énoncé convient pour décrire cette bille.

5
XI. Correction TD 3 exo IV : Ecoulement sur un plan incliné
→
−
−
→ ∂V
1. L’équation locale de conservation de la masse s’écrit div(ρ V ) + = 0. Ici le fluide est un liquide donc
→
−∂t
il est en écoulement incompressible, et cette équation conduit à div V = 0.
→
− −
→
On a V = V (x, y)− →
ex donc Vx = V (x, y), Vy = 0 et Vz = 0, on calcule la divergence par div V =
∂Vx ∂Vy ∂Vz ∂V (x, y)
+ + = = 0 : on en déduit donc que V ne peut pas dépendre de x, la vitesse
∂x ∂y ∂z ∂x
→
−
s’écrit donc V = V (y)−→
ex .

→
−
−
→ ∂V → −−→−
− → − → −
→
On calcule l’accélération par a = + ( V .grad V ). V . Faites le calcul, on trouve bien 0 .
∂t
→
− −−→ →
−
2. L’équation de Navier-Stokes s’écrit donc 0 = ρ− →g − gradP + η∆ V (il n’y a pas d’autres forces, le
référentiel d’étude est galiléen).
−−→ dP −→
Ici P ne dépend que de y donc gradP = ey .
dy
→
− ∂2 −→ ∂ 2 V (y) −
→
Ici la force de viscosité s’écrit η∆ V = η 2 ( V ) = η ex .
∂y ∂y 2
∂ 2 V (y)
On projette sur Ox : 0 = ρg sin α + η
∂y 2
dP
On projette sur Oy : 0 = −ρg cos α − .
dy
dP dP
3. On intègre l’équation 0 = −ρg cos α − ou = −ρg cos α soit P (y) = −ρg cos αy + A avec
dy dy
P (y = e) = P0 : surface libre soit A = ρg cos αe + P0 et P (y) = −ρg cos α(y − e) + P0 .

∂ 2 V (y) ∂ 2 V (y) ρg sin α ρg sin α y 2

4. On intègre l’équation 0 = ρg sin α+η ou = − soit V (y) = − +B.y +C.
∂y 2 ∂y 2 η η 2
Conditions aux limites :
- V (y = 0) = 0 = C (en effet le fluide visqueux s’accroche aux parois, ici les parois sont immobiles donc la
vitesse du fluide est nulle en y = 0).
dV (y)
- La force de viscosité dépend de (voir cours), sur la surface lire, cette force est nulle car l’air n’est
dy
dV ρg sin αe ρg sin αe
pas visqueux donc on a (y = e) = 0 = − + B soit B = .
dy η η
ρg sin α
On en déduit V (y) = 2η (2ey − y 2 ).

5. La vitesse n’est pas uniforme donc on calcule le débit par une intégrale:
ZZ Z y=e Z z=L
→
− →
−
Dv = V .dS. n = V (y)dy.dz en effet −
→
n =− →
ex (sens et direction de l’écoulement) et un point
y=0 z=0
M de la section perpendiculaire à l’écoulement est repéré par y et z qui varient. Donc dS = dy.dz. On
g sin αe3 L
calcule l’intégrale, on trouve Dv = .
3ν