pour le groupe : transport et distribution des hydrocarbures

METICHE MEHDI
Matre Assistant
Charg de cours
guembour rabah

MECANIQUE DES
FLUIDES APPROFONDIE
ECOULEMENTS VISQUEUX, TURBULENCE, COUCHE
LIMITE ET ECOULEMENTS TRANSITOIRES

EXERCICES RESOLUS

U
B

D
dx

Centre Universitaire de Bchar

Octobre 2004

MECANIQUE DES
FLUIDES APPROFONDIE
ECOULEMENTS VISQUEUX, TURBULENCE, COUCHE
LIMITE ET ECOULEMENTS TRANSITOIRES

EXERCICES RESOLUS
transport et distribution des hydrocarbures ::::: universite de boumerdes

METICHE MEHDI
Matre Assistant
Charg de cours
Centre Universitaire de Bchar (Algrie)

A
Amira ..
Signe de reconnaissance ;
Une autre fois encore.

transport et distribution des hydrocarbures ::::: universite de boumerdes

transport et distribution des hydrocarbures

SOMMAIRE

Avant propos---------------------------------------------------------------------------- 1
Ecoulements Visqueux----------------------------------------------------------------

La Turbulence--------------------------------------------------------------------------

19

La Couche Limite----------------------------------------------------------------------

24

Ecoulements Transitoires------------------------------------------------------------

47

Rfrences bibliographiques--------------------------------------------------------- 72

Avant propos

transport et distribution des hydrocarbures ::::: universite de boumerdes

AVANT PROPOS

Le prsent manuel intitul Mcanique des Fluides Approfondie : Exercices Rsolus ;

est ralis suivant le programme du module TEC058, portant le titre de lHydraulique
Avance appartenant au cursus de formation de lIngnieur Hydraulicien.

Le travail est le fruit denseignement du module pendant quatre ans au niveau de

lInstitut dArchitecture, de Gnie Civil et dHydraulique, au niveau du Centre Universitaire
de Bchar. Il comporte alors aussi, les examens rsolus du module.

Son utilit apparat sur son contexte de prsenter des exercices rsolus pour des sujets
dont de tels sont un peu rares. En plus des ingnieurs hydrauliciens, il peut servir aussi la
formation des ingnieurs mcaniciens, chimistes et physiciens, comme il forme un outil
intressant pour les tudiants de poste graduation des diffrentes spcialits prcites, vu le
manque flagrant des ouvrages de ce genre.

Quatre sujets fondamentaux y sont traits : les coulements visqueux, la turbulence, la

couche limite, et les coulements transitoires, dont plus de soixante (60) exercices sont
prsents, exposs, traits et rsolus dune manire dtaille et exhaustive.

transport et distribution des hydrocarbures

1
Mcanique des fluides approfondie : Exercices rsolus

Ecoulements visqueux

ECOULEMENTS VISQUEUX

Exercice 1 :
Donner lquation de mouvement dun coulement permanent, irrotationnel, pour un fluide
incompressible pour :
- un fluide rel ?
- un fluide parfait ?
Solution :
Pour un fluide rel, lquation dEuler scrit :
V

1
+ rotV V + .grad p = grad p + g + V
t

V
Pour un coulement permanent, nous avons :
=0
t

Pour un coulement irrotationnel, nous avons : rotV =0

Avec rot (rotV ) = grad (divV ) div( grad V )
On admet que la force de pesanteur drive dun potentiel U : g = gradU

1
( gradV 2 ) = grad p + gradU + V
2
( gradU +

.grad (U +

grad p +

1
gradV 2 ) = .V
2

1
p + V 2 ) = .V

On a : U = -g.z .grad ( gz +

1
p + V 2 ) = .V

Pour un fluide parfait, nous avons : = 0

1
1
.grad ( gz + p + V 2 ) = 0

2
2
Mcanique des fluides approfondie : Exercices rsolus
transport et distribution des hydrocarbures

Ecoulements visqueux
guembour rabah

1
gz + p + V 2 = Cons tan te

2
Divisant les deux termes de lquation par g ; nous aurons :
z+

1 2
V = Cte
2g

Cette dernire nest autre que lquation de Bernoulli.

Exercice 2 :

Etablir lquation de mouvement dun coulement permanent dun fluide rel incompressible
se produit entre deux plaques planes lisses parallles ?
(Ecoulement monodimensionnel suivant x)
y
b

u
x

Solution :

Lquation dEuler scrit :

u
u
u
u
p
Suivant x : .( + u
+v
+w )=
+ u
t
x
y
z
x
v
v
v
v
p
Suivant y : .( + u + v + w ) =
+ v
t
x
y
z
y
w
w
w
w
p
Suivant z : .(
+u
+v
+w )=
+ w
t
x
y
z
z
On a aussi coulement monodimensionnel :
v = w = 0

2
2
v = v = w = w = 0
y y 2 z z 2

v v 2 v 2 v
= 2 = 2 =0
=

x
z
x
z

2
w w w 2 w
x = y = x 2 = y 2 = 0

Et puisque cest incompressible et permanent, on peut crire de lquation de continuit :

3
Mcanique des fluides approfondie : Exercices rsolus

Ecoulements visqueux
u u u 2 u 2 u
=
=
=
=0
=
x y z x 2 z 2

u = v = w = 0
t t t
Les trois quations dEuler deviennent :

p
p
2u
p
+ 2 = 0,
= 0 , et
=0
x
y
z
y

p
2u
= 2
x
y
p
=0
y
p
=0
z

Qui dmontre que cest un coulement de Couette.

Exercice 3 :

Un coulement dun liquide de viscosit dynamique = 2.10-2 N.s/m2, sur une plaque plane
fixe, est caractris par le profil donn par le schma ci-dessous :

u = 2 m/s

si lpaisseur de lcoulement est e = 5 cm, dterminer la valeur de la contrainte de

cisaillement :
a- la paroi ?
b- une distance de 2 cm de la paroi ?
c- une distance e de la paroi ?
Solution :

Nous avons :
du
=
dy
Et daprs le profil, on a : tg =

u max u
u
2
= u = max . y =
. y = 40 y
e
y
e
0,05

4
Mcanique des fluides approfondie : Exercices rsolus

Ecoulements visqueux
u = 2 m/s

du
= 40
dy
du
= .40 = 2.10-2.40 = 0,8 N/m2
dy
du
est constante, la valeur de la tension de cisaillement est constante nimporte
Puisque
dy
quel point de lcoulement est gale 0,8 N/m2.

Exercice 4 :

La distribution de vitesses pour un coulement permanent incompressible bidimensionnel est

donne par :
x

u = x 2 + y 2

v = y

x2 + y2

a- montrer que cette distribution satisfait lquation de continuit ?

b- montrer que lcoulement vrifie lquation de Laplace si le champ de vitesse drive dun
potentiel ?
Solution :

a- lquation de continuit pour un coulement incompressible bidimensionnel est donne

par :
u v
+
=0
divV = 0
x y
Et nous avons :
2x 2
1
u
= 2
+
x x + y 2 ( x 2 + y 2 ) 2
2y2
1
v
= 2
+
y x + y 2 ( x 2 + y 2 ) 2

5
Mcanique des fluides approfondie : Exercices rsolus
transport et distribution des hydrocarbures ::::: universite de boumerdes

Ecoulements visqueux

2x 2
2y2
2( x 2 + y 2 )
1
1
2
u v
= 2
+
=
+
+
+
+
x y x + y 2 ( x 2 + y 2 ) 2 x 2 + y 2 ( x 2 + y 2 ) 2
x 2 + y 2 (x 2 + y 2 )2

divV =

divV =

2
+2
+ 2
=0
2
x +y
x + y2
2

Donc divV = 0 : lcoulement est continu.

b- nous avons : V = u + v
Si le champ de vitesse drive dun potentiel, on peut crire :
u=

et v =
x
y

V =

+
= grad
x y

Nous avons aussi : divV = 0 div( grad )= 0 = 0

2 2
+
= 0 ; cette dernire est bien vident, lquation de Laplace
x 2 y 2

Exercice 5 :

Le profil de vitesse pour un coulement plan dun liquide, est donn par :
V(y) = 3.y3 + 2.y2.
Si la viscosit dynamique du liquide est = 3,5.10-2 N.s/m2.
Calculer la valeur de la tension de cisaillement la paroi et 30 cm de celle-ci ?
Solution :

On a :
du
dy
Avec V = 3.y3 +2.y2.

dV
= 9.y2 + 4.y
dy

A la paroi : 0 =

du
dy

y =0

= 3,5.10-2.(9.0 + 4.0) = 0

6
Mcanique des fluides approfondie : Exercices rsolus

Ecoulements visqueux
A 30 cm de la paroi : 0,3 =

du
dy

y =0,3

= 3,5.10-2.(9.0,3 + 4.0,3) = 7,035.10-2 N/m2

Exercice 6 :

Soit un coulement plan dun liquide de viscosit cinmatique = 5.10-4 m2/s et de masse
volumique = 103 kg/m3 sur une plaque plane.
Le profil de vitesse est donn par :

1 3
y .
2
Dterminer la valeur de la tension de cisaillement :
- la paroi ?
- 7 cm de la paroi ?
V(y) =

Solution :

On a :

dV
dy

Avec = . = 5.10 4.10 3 = 5.10 1 kg / m = 5 N .s / m 2

dV 3 2
= y
dy 2

y =0 =

dV
dy

y =0, 7 =

dV
dy

y =0

=0

y =0,7

= 5.(3/2.0,072) = 3,67.10-2 N.m2

Exercice 7 :

Soit un coulement dont le potentiel des vitesses est donn par :

= x2 -2.y y2 avec v = grad
Dmontrer que lcoulement est :
- bidimensionnel
- permanent
- continu (vrifiant lquation de continuit) ?
Solution :

- Nous avons : V(u,v,w)

7
Mcanique des fluides approfondie : Exercices rsolus
transport et distribution des hydrocarbures ::::: universite de boumerdes

Ecoulements visqueux

Avec :

= 2.x
u=
x

v=
= -2.y-2
y

w=
=0
z
Il est vidant, que la vitesse ne dpend pas de z : lcoulement est donc bidimensionnel.

V grad
=
=0
t
t
Il est vident, que la vitesse ne dpend pas du temps : lcoulement est donc permanent.
- Dautre part, nous avons :

- on a aussi : divV =

u v w
+
+
= 2 2 +0 = 0
x y z

Lquation de continuit est vrifie, donc lcoulement est continu.

Exercice 8 :

Supposant que les composantes du vecteur vitesse sont :

u = x2 + z2
v = x2 + y2
Dterminer les composantes du vecteur vitesse suivant la direction z, qui satisfont lquation
de continuit ?
Solution :

Pour satisfaire lquation de continuit, divV = 0.

u v w
w
u v
+
+
=0
=-[ +
]
x y z
z
x y

Nous avons :
u = x2 + z2
v = x2 + y2

u
v
= 2.x et
= 2.y
x
y

w
= - (x + y) w = -2(x + y).z
z
8
Mcanique des fluides approfondie : Exercices rsolus
transport et distribution des hydrocarbures ::::: universite de boumerdes

transport et distribution des hydrocarbures ::::: universite de boumerdes

Ecoulements visqueux
Exercice 9 :

Le profil de vitesses pour un coulement plan sur une plaque plane fixe, dun liquide de
viscosit dynamique = 10-2 N.S/m2, est donn par :
V(y) = y3 + 2.y2 + 5.y
Dterminer la valeur de la tension de cisaillement 10 cm de la paroi ?
Solution :

On a :
dV
dy
Et nous avons : V(y) = y3 + 2.y2 + 5.y

dV
= 3. y 2 + 4. y + 5
dy

dV
= 2.10-2.[3.0,12 + 4.0,1 + 5] = 2.10-2.[0,03 + 0,4 + 5] = 10,86. 10-2 N/m2.
dy

= 10,86. 10-2 N/m2.

Exercice 10 :

Sur une plaque plane lisse, faisant avec lhorizontal un angle , en mouvement permanent
bidimensionnel tabli et sous une paisseur a, coule sous leffet de la pesanteur, un liquide de
viscosit .
Dterminer le dbit par unit de largeur et la vitesse dbitante (moyenne) ?
a
z

w
u

Solution :

Les quations de Navier-Stokes scrivent 2 dimensions :

9
Mcanique des fluides approfondie : Exercices rsolus
transport et distribution des hydrocarbures ::::: universite de boumerdes

transport et distribution des hydrocarbures ::::: universite de boumerdes

Ecoulements visqueux

p g
2u 2u
u
u
u
+ .( 2 + 2 )
+u
+w )=
x
t
x
z
x
z
2
p g
w 2w
w
w
w
Suivant z : .(
+ .( 2 + 2 )
+u
+w )=
z
t
x
z
x
z

Suivant x : .(

De lquation de continuit, on crit :

u w
+
=0
x z
Lcoulement est permanent :

Lcoulement est tabli :

=0
t

=0
=0
x
z

w est nul et u ne dpend que de z, ce qui nous permet dcrire :

p g

2u
= 2
x
z
p g
= 0 avec Pg = p + .g.h
z

p
h
p
+ .g . =0 soit :
+ .[Link] =0
z
z
z

Aprs intgration de cette quation, nous obtenons :

p = .[Link] + f(x)
Puisque la pression est constante sur la surface libre : z = a f(x) = 0.

2 u p
h
= + .g. = 0 - .[Link]
2
x
x
z

En intgrant, nous trouvons :

u
= .[Link] + A
z
z2
.u = - .g. .Sin + A.z + B
2
Sur la paroi nous avons : z = 0, u = 0 B = 0

A la surface libre, la contrainte est nulle : z = a :

u
= 0 A = .[Link]
z

10
Mcanique des fluides approfondie : Exercices rsolus
transport et distribution des hydrocarbures ::::: universite de boumerdes

transport et distribution des hydrocarbures

Ecoulements visqueux
u = .g .

Sin
z2
.(a.z )
u
2
a

Le dbit par unit de largeur : qv = [Link] = .g .

0

a3
a 3 .g
.Sin
.Sin =
3
3

Et la vitesse dbitante (moyenne) sera :

Um = qv / a =

a 2 .g
.Sin
3

Exercice 11 :

1 dP 2
. ( R r 2 ) , dune huile dans un
4 dx
tube de 12 mm de diamtre, de longueur L, inclin dun angle dont Sin = 0,0445.
Si la pression statique lintrieure est constante tout le long du tube, et le dbit mesur gal
20 l/h, dterminer la viscosit cinmatique de cette huile ?
Un coulement de Poiseuille (rgime laminaire, u =

u
P2
P1
R

Solution :

Si P2 est rfrence des pressions, on crit :

P2 = p2 , et : P1 = p1 + .g .[Link]
P = P1 - P2 = p1 + .g .[Link] - p2
A pression statique constante : p1 = p2
P = .g .[Link]

P
= .[Link]
L

Puisque le rgime est laminaire : dQv = u. 2. .[Link]

11
Mcanique des fluides approfondie : Exercices rsolus

Ecoulements visqueux
1 dp
1

Qv = . .( R 2 r 2 ).2 .[Link] =
. .g .Sin .( R 2 r 2 ).2 .rdr =
. .g .Sin .D 4
4u dx
4
128
0
0
R

On a aussi : =

[Link]

avec =
. .g .Sin .D 4
[Link]

.[Link] .D 4 = 4.10-5 m2/s .

Exercice 12 :

Un rcipient cylindrique de 5 cm de diamtre, laisse couler un liquide visqueux quil contient

par un tube horizontal de diamtre 1 mm et de longueur 40 cm. Sachant quil a fallu 75
minutes pour que la charge en amont du tube passe de 5 2,5 cm, dterminer la viscosit
cinmatique du liquide (on ngligera les effets dus aux extrmits du tube ; lcoulement dans
P 4
.
d ?
le tube est permanent de type Poiseuille : qv =
128 L.
5cm
2,5 cm

5 cm

Solution :

Pour un coulement permanent, on crit :

P = .g.h
qv =

g .h 4
d
128 L.
.

De lquation de continuit, on crit :

qv = S

.d 2 h
.5 2 h
g .h 4
h
. =
. =
=
.
d
4 t
4 t 128 L.
t

h g 4 4
h d 4 g
.
.h = 0
=
0

+
=
.
d .
.
h
t 128 L.
t 800 L.
.5 2

d4 g
h
.
.dt
=
Ou bien :
800 L.
h
12
Mcanique des fluides approfondie : Exercices rsolus

transport et distribution des hydrocarbures

Ecoulements visqueux
Lintgration de cette quation, donne :
h1 d 4 g
ln =
.
.T
h2 800 L.

d4 g
T
(10 3 ) 4
10
75.60
. .
=
.
.
= 2.10-6 m2/s
2
800 L h1
800 40.10
5
ln
ln
2,5
h2

= 2.10-6 m2/s = 2.10-2 Stokes

Exercice 13 :

Soit un tube cylindrique de 3 km de long, de 10 cm de diamtre, parcouru par un liquide de

viscosit dynamique = 0,4 Poise.
On suppose que la distribution des vitesses dans la section droite du tube est donne par
lquation parabolique u = 10.( y y 2 ) en unit CGS.
u tant la vitesse la distance y de la paroi ?
Calculer :
- la force de frottement visqueux par unit de surface contre la paroi ?
- la force de frottement visqueux par unit de surface 2 cm de la paroi ?
- la force totale de frottement sexerant sur le tube ?

v
y

Solution :
1- =

du
= .(10 20.y)
dy

A la paroi, nous avons :

y=0;u=0
y = 10 cm ; u = 0.
= .10 = 0,4.10 = 4 dynes/cm2.

2- pour y = 2 cm :
du
= 10 20.0,2 = 10 4 = 6 .
dy

13
Mcanique des fluides approfondie : Exercices rsolus

Ecoulements visqueux

du
= 0,4.6 = 2,4 dynes/cm2.
dy

3- la surface totale du tube :

S = .D.L = 3. .106
La force de frottement sera :
F = 0 .S =

du
. .D.L = 4.3 .106 = 3,77. 107 dunes = 377 Newtons
dy y =0

Exercice 14 :

Un rfrigrant est compos dun groupe de 100 tubes cylindriques parallles, de diamtre D =
1 cm et de longueur L = 4 m. A la vitesse V = 2 m/s, on y fait circuler de lhuile ( =900
Kg/cm3, 1 = 3.10-2 Poiseuille lentre et varie jusqu la sortie, cause de refroidissement,
pour atteindre 2 = 10-1 Pl ).
1- montrer que lcoulement est laminaire ?
2- donner -pour un tube- la relation ( x) : f ( 1 , 2 , L, x) avec x la variable longueur entre
lentre et la sortie du tube ?
dP 32.V ( 1 + 2 )
3- lcoulement est de type Poiseuille dont
= 2 .
.L , calculer la variation de
dx
2
D
pression entre lentre et la sortie dun tube ?
4- en dduire la puissance ncessaire W, pour faire couler lhuile dans un tube, ensuite pour
lensemble du groupe ?
5- donner la fonction W de n tubes sans tenir compte du diamtre D ?
Solution :

1
1

2
x

x
L
1- on a Dmax =

.V .D 9.10 .2.10
=

3.10 2
2

= 600 < 2000

le rgime est laminaire.

14
Mcanique des fluides approfondie : Exercices rsolus

Ecoulements visqueux
2- (x) =

x
.( 1 - 2) + 1.
L

dP 32.V ( 1 + 2 )
32.2 (3.10 2.10 5 + 10 1.10 5 )
= 2 .
3- P =
.L =
.4 =1,66 bars
.
dx
2
2
D
0,012

4- W = Qv. P =

.(10-2) .2.1,66.105 = 2610 Wat.

5- W= n. P. Qv = n.32.

V ( 1 + 2 ) D 2
.
.L.
.V = 4 .L.V 2 .( 1 + 2 ).n
2
4
D2

W= 4 .L.V 2 .( 1 + 2 ).n

Exercice 15 :

Soit la figure ci-aprs, reprsentant un cylindre plein de rayon R1 lintrieur dun cylindre
creux de rayon intrieur R2 et entre lesquels existe un coulement laminaire.
1- lcoulement est de type Couette dont la distribution de vitesse est donne par
u=

P 2
.r + C1. ln(r ) + C 2
4 .L

u = 0, pour : r = R1

Avec u = 0, pou : r = R2 ,
et : P = P P
1
2

R2

r
R1

Etablir la relation donnant le rayon pour lequel la vitesse est maximale ?

2- trouver la loi du dbit volumique de cet coulement annulaire et en dduire la vitesse
moyenne de dbit ?

Solution :

1- Nous posons : K =

P P2
P
= 1
4.L
4 .L
15

Mcanique des fluides approfondie : Exercices rsolus

transport et distribution des hydrocarbures ::::: universite de boumerdes

Ecoulements visqueux
2

C1.lnR1 + C2 K.R1 = 0
2

Et C1.lnR2 + C2 K.R2 = 0
2

Ce qui donne : C1 = K .

u=

P1 P2
4. .L

R1 R2
R
ln 1
R2

et C2 = K.R1 -

( R1 R2 ). ln R1
R
ln 1
R2

2
ln(r / R2 )
2
2
2
( R2 r ) ( R2 R1 ).

ln( R1 / R2 )

Le rayon pour lequel, la vitesse est maximale :

2

R R2
-[Link] - 1
/ Rmax = 0 Rmax =
R1
ln
R2

du
=0
dr
2

R2 R1
2. ln( R1 / R2 )

2- le dbit volumique scrit comme :

Qv =

R2

R2

R1

R1

[Link] = u.2 .rdr

R2

Qv =

R1

Qv =

P1 P2
4. .L

2
ln(r / R2 )
2
2
2
( R2 r ) ( R2 R1 ).
.2 .rdr
ln( R1 / R2 )

2
2
( R R1 )
.( P1 P2 ) 2
2
2
2
.( R2 R1 ).[ R2 R1 2
]
8. .L
ln( R2 / R1 )

La vitesse moyenne sera :

2
2
QV
QV
( P1 P2 ) 2
( R2 R1 )
2
. R2 R1
=
V=
=

8. .L
ln( R2 / R1 )
A .( R2 2 R1 2 )

Exercice 16 :

Un viscosimtre cylindres coaxiaux a les dimensions reprsentes sur la figure en cidessous.

On se propose de calculer la viscosit dynamique du liquida existant entre les deux
cylindres, sachant que le couple mesur C vaut 9.103 dynes = 9.10-4 N.
On donne :
16
Mcanique des fluides approfondie : Exercices rsolus

transport et distribution des hydrocarbures

Ecoulements visqueux

R1 = 77 mm, R2 = 80 mm, e = 3 mm.

1
2
L = 200 mm, a = 4 mm, f =
Tours / sec onde =
(1 / s )
30
30
La vitesse de Couette est donne par :
V =

1
2

R2 R1

.[r.( f .R2 )
2

R1 .R2
.f ]
r

R2
L

r
R1

e
a

Solution :

Le couple C1 cre par les surfaces latrales est donn par :

V
2
C1 = F.R2 = .S. R2 = .2 . R2. R2 = .
r = R2 .2 .R2
r
Avec V = R2.f
C1 = 4 . .L. f .

R2 .R1
2

R2 R1

Pour le couple C2 du au fond du cylindre, il est donn par :

dC2 = [Link] = r. .Ds
o : contrainte tangentielle au fond du cylindre et dS = 2 .rdr
dC2 = r. . 2 .rdr = r. .

avec = .

V
.2 .[Link]
r

V
r. f
V
= . = .
a
a
r

dC2 = r. .

car : V = r.f

r. f
2. . . f 3
.2 .[Link] =
.r dr
a
a
17
Mcanique des fluides approfondie : Exercices rsolus

Ecoulements visqueux
R1

C2 =

2. . . f R1
2. . . f 3
.
.r dr =
a
4
a

Le couple total mesur sera :

C = C1 + C2
C = 4 . .L. f .

= C / 2 . f .[

R2 .R1
2

R2 R1

2.L.R2 .R1
2

R2 R1

R
2.L.R2 .R1
2. . . f R1
.
= 2 . . f .[
+ 1 ]
2
2
a
4
4a
R2 R1
4

R
+ 1 ] = 0,199 Poise.
4a

18
Mcanique des fluides approfondie : Exercices rsolus

La Turbulence

LA TURBULENCE

Exercice 1 :
Dmontrer, pour un coulement turbulent que :
a- U .V = U .V + U ' .V '
b- U (t ) = U

c- U .V .W = U . V . W + U .V ' .W ' + V .U ' .W ' + W .U ' .V ' + U ' .V ' .W '

Solution :

a- U .V = (U + U ' ).(V + V ' ) = U .V + U .V '+U '.V + U '.V ' = U .V + U .V ' + U '.V + U '.V '
= U .V + U .V ' + U '.V + U '.V '
Et puisque U .V ' = U '.V = 0
U .V = U .V + U ' .V '

b- U (t ) = U + U ' = U + U ' = U + U '

Et puisque U ' = 0 U (t ) = U

c- U .V .W = (U + U ' ).(V + V ' ).(W + W ' )

= U .V .W + U .V '.W + U '.V .W + U '.V '.W + U .V .W '+U .V '.W '+U '.V .W '+U '.V '.W '
= U .V .W + U .V '.W + U '.V .W + U '.V '.W + U .V .W ' + U .V '.W ' + U '.V .W ' + U '.V '.W '
Et puisque U .V '.W = U '.V .W = U .V .W ' = 0
U .V .W = U . V . W + U .V ' .W ' + V .U ' .W ' + W .U ' .V ' + U ' .V ' .W '
19
Mcanique des fluides approfondie : Exercices rsolus

La Turbulence
Exercice 2 :

Montrer pour un coulement de leau 15 C dans une conduite circulaire de diamtre D, que
la condition de la turbulence du rgime dcoulement est conditionne par :
0,0026
V>
, tout en considrant que le nombre de Reynolds critique est de 2300 ?
D
Montrer que pour un diamtre de lordre de 1 cm, le rgime turbulent est assur pour des
vitesses V > 0,3 m/s ?
Solution :

On a : =

0,0178
1 + 0,0337.t + 0,000221.t 2

A 15 C : = 1,14.10-6 m2/s
Et on a : =

V .D

> 2300 : condition de la turbulence.

V>

2300. 2300.1,14.10 6 0,0026

=
=
D
D
D

V>

0,0026
D

Pour D = 1 cm V > 0,26 m/s

Pour V > 0,3 m/s : le rgime turbulent est suffisamment assur.

Exercice 3 :

Dmontrer pour un coulement turbulent dun fluide incompressible, que lquation de

continuit scrit :
u ' v ' w '
+
+
=0?
x
y
z
Solution :

Pour un fluide incompressible : divV = 0,

u v w
(u + u ' ) (v + v' ) ( w + w' )
+
+
=0
=0
+
+
x y z
z
y
x

u v w u ' v' w'

+
+
+
+
+
=0
x y z
x y z

20
Mcanique des fluides approfondie : Exercices rsolus

La Turbulence

Et puisque

u v w
=
=
=0
x y z

u ' v' w'

+
+
=0
x y z

Exercice 4 :

Dmontrer que pour un coulement turbulent :

1-

u
=0
t

2-

u
vu
u
= v. + v .
y
y
y
2

3- u 2 = u + u 2
Solution :

u u u ' u u ' u u
=
+
=
+
=
+
1t t t
t t
t
t

'

'

u
u
= 0 et
=0
Et puisque
t
t

u
=0
t

2-

(u + u ' )
u
u
u '
vu
u '
= (v + v' ).
= v. + v'. + v.
+ v'.
y
y
y
y
y
y

= v.

u
u
u '
u '
+ v'. + v.
+ v'.
y
y
y
y

Et puisque v' =

u '
=0
y

u
u '
vu
= v. + v'.
y
y
y
2

3- u.v = (u + u ' ).(u + u ' ) = u + u ' 2 + 2.u.u '

21
Mcanique des fluides approfondie : Exercices rsolus

La Turbulence
Et puisque .u.u ' = u.u ' = 0
2

u 2 = u + u2

Exercice 5 :

Comparer entre le tenseur des tensions visqueuses dun coulement laminaire, et un

coulement turbulent pour un fluide incompressible ?
Solution :

Pour un fluide incompressible, nous avons : div(V) = 0.

Le tenseur des tensions sexprime comme :

xx

xy

xz

yx

yy

yz

zx

zy

zz

Pour un coulement laminaire, les lments du tenseur scrivent comme :

2
u
u
= p x + 2
x
3
x
2
v
v
yy = p y . .divV + 2 . = p y + 2
y
3
y
2
w
w
zz = p z . .divV + 2 .
= p z + 2
z
3
z
u v
xy = yx = .( + )
y x
u w
xz = zx = .( + )
z x
v w
yz = zy = .( + )
z y

xx = p x . .divV + 2 .

Pour un coulement turbulent, les lments du tenseur scrivent comme :

u
+ .u ' 2
x
v
yy = p y + 2 . + .v' 2
y
w
zz = p z + 2 . + .w' 2
z
u v
xy = yx = .( + ) + - .u '.v'
y x

xx = p x + 2 .

22
Mcanique des fluides approfondie : Exercices rsolus
transport et distribution des hydrocarbures ::::: universite de boumerdes

La Turbulence
u w
+
) + - .u '.w'
z x
v w
= .( +
) + - .v'.w'
z y

xz = zx = .(
yz = zy

Le tenseur des tensions visqueuses dun coulement turbulent comme est prsent, est la
somme du tenseur dun coulement laminaire et le tenseur de Reynolds.
Car le tenseur de Reynolds scrit comme :

xx = .u ' 2

yy = .v' 2
zz = .w' 2

xy = yx = - .u '.v'
xz = zx = - .u '.w'
yz = zy = - .v'.w'

transport et distribution des hydrocarbures ::::: universite de boumerdes

23
Mcanique des fluides approfondie : Exercices rsolus

La Couche limite

LA COUCHE LIMITE

Exercice 1 :
Une plaque plane de 3 cm de longueur, et de 30 cm de long est remorque paralllement
elle-mme dans le sens de sa longueur dans leau la vitesse de 6 m/s.
Dterminer la force de frottement sexerant sur lune des faces de la plaque, et la force
sexerant sur les trois premiers mtres de la plaque ?
On donne = 10 6 m 2 / s.
Solution :

6.30
= 1,8.108
6

10
Le coefficient de frottement moyen est :

Nous avons : L =

U .L

Pour > 10 7 on admet que pour une plaque plane lisse parallle U de longueur L :
Cx = 0,455. (log10 L ) 2,58

0,455
= 1,965.10 3
231,8
Et nous avons lair dune face de la plaque : S = 3.30 = 90 m2.

Cx =

La force sexerant sur une face :

U2
= 1,965.10-3.90.1000.36/2 = 3183 N.
2
Si le nombre critique de transition correspond = 5.10 5
La position de la ligne de transition est situe la distance x du bord dattaque :

F = Cx. .S .

U.x

6.x
x = 0,083 m.

10 6
En ngligeant linfluence du frottement dans la couche limite laminaire (x = 0,083 m)

5.10 5 =

Cx sur les 3 m de la plaque : Cx = 0,455. (log10 L ) 2,58

U .L 6.30
= 6 = 1,8.108
Avec L =

10
-3
Cx = 2,74.10 .
La force de frottement sur S :

24
Mcanique des fluides approfondie : Exercices rsolus

La Couche limite
S = 3.3 = 9 m2
Fx = Cx. .S .

U2
= 2,74.10-3.9.1000.36/2 = 444 N.
2

Exercice 2 :

Dans lcoulement laminaire dun fluide sur une plaque mince et plate, on admet que la
distribution des vitesses dans la couche limite rpond lquation :
u
y
= Sin( . )
2
U
Avec :
U : vitesse du fluide libre (coulement)
u : vitesse la distance y de la paroi
: paisseur de la couche limite.
1- calculer latralement : x , x , lx , H , C fx , C F , et local x ?
2- on donne :
U = 15 m/s
L = 4.10-2 m.
= 10-6 m2 /s
= 1,2 kg/m3.
1
2
Sin xdx = 2.( x [Link])
Solution :
1du
a- =
dy

* = (1
0

= .U .

u
).dy =
U

* = +

b- =

y =0

y
.Cos ( . )
2
2

= .U .

y =0

1 Sin( 2 . ).dy = +
0

.(Cos

Cos 0) =

y
.Cos( . )

= 0,363.

u
u
y
y
0 U (1 U ).dy = 0 Sin( 2 . ).[1 Sin( 2 . )].dy =

y
2 y
0 Sin( 2 . ).dy - 0 Sin ( 2 . ).dy

1
On a aussi : Sin 2 xdx = .( x [Link])
2
1
Sin 2 udu = .(u [Link] )
2

25
Mcanique des fluides approfondie : Exercices rsolus

La Couche limite
= -

2
2
.Cos Cos 0

2

y
y
1 y
. .( . Sin . .Cos . ) =
2
2 0
2 2

2 1
2
.( 0) (0 0) =
= .( ) = 0,137.
2
2
2

= 0,137.

c- local = .U 2 .

d
d
= .U 2 .0,137.
= .U.
2.
dx
dx

d
d ( 2 )
= .

= 22,931.
dx U 0,137
dx
U
.x
2 = 22,931.
, multipliant le deuxime terme de cette quation par x/x, nous obtenons :
U
.x 2 22,931 2
x
=
2 = 22,931.
.x =
.4,789
x
U .x
x
local = 2 .

= 0,137.

x
x

.4,789 = 0,656.

Et * devient : * = 1,738.
d- H =

x
x

*
= 2,649

e- Cf = 2.

f- CF =

d 0,656
=
dx
x

2. 1,312
=
L
L

Et local devient :

local = .U .

. .U .
=
= 0,328. .U . x / x
2.
2

2- application numrique :
*
L = 6.10 5 ; L = 0,247 mm ; L = 0,0897mm
L = 0,0338mm ; CF = 1,694.10-3 ; Cf = 8,468.10-4.

26
Mcanique des fluides approfondie : Exercices rsolus

La Couche limite
Exercice 3 :

Un filtre nid dabeilles, plac en avant du collecteur dune soufflerie, est constitu de
mailles carres dont le cot est de 3 m, la profondeur des lamelles dans le sens de la vitesse
est de 4 cm. On suppose que chaque lamelle se comporte comme une plaque plane la vitesse
U = 15 m/s.
Les caractristiques de lair sont : = 1,2 kg/m3 et = 15.10-2 St.
a- montrer que la couche limite nest pas en rgime totalement turbulent ( eL < 5.105) ?
b- calculer lpaisseur de la couche limite ?
c- en dduire les valeurs de l , , Cf, CF, local ?
d- calculer la somme des forces de frottement sur les quatre faces dune maille ?
e- quelle est la chute de pression, la traverse du filtre, pour une maille ?
f- en dduire le coefficient de perte de charge singulire pour ce filtre ?
Solution :

a- eL =
b- =

U .L

4,646.L
eL

15.4.10 2
= 4.10 4 > 5.10 5
15.10 6
=

4,646.4.10 2
= 0,928mm
2.10 2

c- 1 = 0,375. = 0,348mm

= 0,139. = 0,129mm
Cf =

CF =

0,646
eL
1,292

local =

eL

0,646
= 3,23.10 3
2
2.10

= 6,46.10 3 = 2.C f

3.U 3.15.10 6.15.1,2

=
= 0,436 Pa
2.
2.0,928.10 3

d- si la surface de frottement la paroi est S :

1
1
Fv = p .S = .C F . .U 2 .S = 6,46.10 3.1,2.15 2.12.10 4 = 1,046.10 3 N
2
2
4. Fv = 4.1,046.10-3 = 4,186.10-3 N

e- si A est la section de passage : A = 32.10-4.

[Link] 4,186.10 3
P =
= 2 4 = 4,65 Pa
A
3 .10
27
Mcanique des fluides approfondie : Exercices rsolus

transport et distribution des hydrocarbures ::::: universite de boumerdes

La Couche limite

f- hp = .

U 2 P
=
2 g g

2.P 2.4,65
=
= 0,034
.U 2 1,2.15 2

Exercice 4 :

Dans un coulement laminaire dun fluide sur une plaque mince et plate, on admet que le
profil des vitesses dans la couche limite rpond lquation :
u 3 y 1 y
= . .
U 2 2

Avec :
U : vitesse dcoulement libre.
u : vitesse la distance y de la paroi.
: paisseur de la couche limite.
a- pour un calcul latral appliqu une unit de largeur, montrer que lpaisseur de la couche
limite la sortie de la plaque est donne par :
L
= 4,646.

eI

1/ 2

Sachant que :

P = local = .

du
et o L : est la largeur de la plaque.
dy y =0

b- en dduire les paisseurs de quantit de mouvement , de dplacement * , les coefficients

de frottement local Cf et moyen CF ainsi que le facteur de forme H ?
Solution :

a- on a local

u 3 y 1 y
du
= .
avec :
= . .
U 2 2
dy y =0

local =

3 U
3 U y
3 U
.u ( ) .u. . = . .
2
2
2

Dautre part : local = .U 2 .

d
dx

28
Mcanique des fluides approfondie : Exercices rsolus
transport et distribution des hydrocarbures ::::: universite de boumerdes

La Couche limite

u
u
3 y 1 y
3 y 1 y
Avec = (1 ).dy = [ . ( ) 3 ].[1 . + .( ) 3 ]dy
U
U
2 2
2 2
0
0

3

.
= 0,139.
10
8 28

local = .U 2 .0,139.
2. .

d 3 U
= ..
dx 2

d d ( 2 )
=
= 21,583.
dx
dx
U

On a : x = 0 : = 0 2 = 21,583.

Pour la longueur L : =

21,583.

.x

x2
= 21,583.
U
ex

L2
el

b- lpaisseur de quantit de mouvement sera :

= 0,139. = 0,646.

eL
- lpaisseur de dplacement sera :
u
3 y 1 y
= l = (1 ).dy = [1 . + ( ) 3 ].dy
U
2 2
0
0
*

Par intgration, on obtient :

3
8

* = . = 0,375.
Pour la longueur L : * = l = 1,742.

L
eL

1/ 2

d
d
d
d (4,646.x)
1 4,646

= .
avec
- Cf = 2.
= 2.0,139.
=
U
2 ex
dx dx
dx
dx
( )1 / 2

1 4,646 0,646
=
et Cf = 2.0,139. .
2 ex
ex

29
Mcanique des fluides approfondie : Exercices rsolus

transport et distribution des hydrocarbures ::::: universite de boumerdes

La Couche limite

- CF = 2.

-H=

= 2.

0,646
eL

1,292
eL

l 1,742
=
= 2,7
0,646

Exercice 5 :

Si lexpression de la distribution de la vitesse dun coulement laminaire sur une plaque plane
est donne par :
u
u
u
=
=
= 2 2 3 + 4
U U E Ul
Avec =

1- donner lexpression du frottement parital ?

2- calculer les paisseurs de dplacement * , et de quantit de mouvement ?
3- en dduire lpaisseur de la couche limite, pour une plaque de longueur L, en fonction du
nombre de Reynolds local L ?
Solution :

D
L

30
Mcanique des fluides approfondie : Exercices rsolus
transport et distribution des hydrocarbures ::::: universite de boumerdes

transport et distribution des hydrocarbures ::::: universite de boumerdes

La Couche limite
Section
AD

Bilan de masse
0

Bilan quantit de mouvement

0

. U .dy

AB

.L. U 2 .dy

. [Link]

CD

.L. u 2 .dy

BC

. (U u ).dy

.L. U .(U u ).dy

.L u.(u U ).dy

p .dx = . u.(U u ).dy

0

p ( x) = .

d (U 2 . )
u
u
avec = (1 ).dy
dx
U
U
0

- lpaisseur de dplacement est : = (1

*

- lpaisseur de dplacement est : =

u
3
).dy =
.
10
U

U (1 U ).dy =
0

37
.
315

En la comparant avec lexpression de p :

p ( x) = .

d (U 2 . )
dx

et on a : p = .

2
du
= . .U
dy y =0

2
d (U 2 . )
. .U = .

dx
Et dautre part, on a : .d =

Et comme x =

U .x

dx

2
2

.x =

.x 2
U .x

x2
=

2 x

31
Mcanique des fluides approfondie : Exercices rsolus

La Couche limite
2 .x

- le coefficient de frottement local :

2

. .U
4.
=
=
=
Cf =
1
1

.U
2
2
. .U
. .U
p

Cf =

Cf =

2. 2 . . x
x.U

4.
2 .x
x

.U

2. 2 . x
x

2. 2
x

- le coefficient de frottement global :

L

CF = C f .dx =
0

2. 2 .L
L

Exercice 6 :

Un sous marin a une longueur L = 84 m et une surface totale de coque de 1800 m2.
Calculer la rsistance due aux forces de frottement visqueux sexerant sur la coque, quand en
plonge, le sous marin a une vitesse V = 5 m/s ? Pour faire ce calcul, on admettra qu la
valeur Cx du coefficient moyen de frottement tabli pour une plaque plane est applicable ici.
On a dailleurs :
-2,58
V .L
> 107.
pour L =
Cx = 0,455.(lg10 L )

Et on a :

= 1025 kg/m3, = 1,2.10-6 m2/s.

Solution :

On a
L =

V .L

5.84
= 3,5.10 8 > 10 7
6
1,2.10

Cx = 0,455.[log10 L ]

2 , 58

= 1,79.10 3

La force de frottement sera :

32
Mcanique des fluides approfondie : Exercices rsolus

La Couche limite
Fx = Cx . .

V2
52
.S = 1,79.10 3.1025. .1800 41280 N.
2
2

Exercice 7 :

Soit un coulement laminaire permanent, bidimensionnel dun fluide incompressible le long

dune paroi quelconque dont le rayon de courbure est trs suprieur lpaisseur de la
couche limite, de sorte que les quations de Prandtl soient valables :
u

u
u
1 dp
2u
+v
=
+ 2
x
y
dx
y

u v
+
=0
x y
U2
Et p +
= Cte pour y =
2
u
Avec U
v
1- calculer le dbit en masse travers la couche limite en fonction de et * ?
2- calculer le dbit de quantit de mouvement en fonction de , * et ?
3- utiliser ces relations et le thorme des quantits de mouvement appliqu un lment
ABCD, dpaisseur dx, pour tablir lquation :
U d ( 2 )
2 dU
*
.
= .
.(2 + ) + 0 .
2 dx
dx

(A)

4- si le profil des vitesses est donn par :

2

u
y m y
= L. + . : 0 y
U
2

(B)

u
=1:y>
U
Avec L et m : paramtres ne dpendant que de x.
a- calculer et montrer que L et m sont lis par une relation de type : L = (m) ?
b- calculer

*
*
et montrer quon peut avoir une relation de la forme :
= (m) ?

c- en utilisant les quations de Prandtl et les conditions la paroi, montrer que :

33
Mcanique des fluides approfondie : Exercices rsolus

La Couche limite

2 dU
.
?
dx

m=-

d- vrifier que le 2eme membre de lquation (A) ne dpend que de m ?

0,45 + 6m
. Donner dans ces conditions, une
e- ce 2eme membre est approximativement gal
2
formule permettant de calculer en fonction de U ?
Calculer , et m lorsque U = [Link](kx) o U0 = constante et 0 k.x ?

U
B

C u

A
x

D
dx

Solution :

1- nous avons :

* = (1
0

u
u
u
)dy = dy - dy = - dy
U
U
U
0
0
0

Par consquent : q =

.[Link] = .U.(

-)

2- =

u
u
u
u2
u2
*

=
=

(
1

).
dy
.dy
.dy
0 U U
0 U
0 U 2
0 U 2 .dy

Par consquent : le dbit de quantit de mouvement : M =

.u

.dy = .U 2 . ( - * - )

3- considrant le volume lmentaire ABCD, dpaisseur unitaire :

Le dbit entrant par AB est : q
q
dx
Le dbit sortant CD est : q +
x
q
dx la vitesse C
Le dbit entrant par BC est :
x
Le dbit sortant par AD est : 0

34
Mcanique des fluides approfondie : Exercices rsolus

transport et distribution des hydrocarbures ::::: universite de boumerdes

La Couche limite

q
dx
x

q+

q
A

q
dx
x

D
0

En projection parallle la paroi ; les dbits de quantit de mouvement seront :

Par AB entre le dbit : M,

M
dx
x
q
Par BC entre le dbit : U. dx
x
Par AD cest 0
Par CD sort le dbit : M +

U.
B

q
dx
x

M+

M
A

M
dx
x

D
0

Lapplication du thorme des quantit de mouvement donne :

- M + (M +

Soit :

M
q
dx ) - U. dx = p. - (p+ dx ).( +
.dx ) + p. dx 0 .dx
x
x
x
x
x

M
q
p
- U. = . 0
x
x
x

En utilisant les valeurs de M et q calcules prcdemment :

.U 2 .

U
p

+ .U .
.( * 2 ) = . 0
x
x
x

Et puisque p + .

U
p
U2
= Cte
= .U .
2
x
x
35
Mcanique des fluides approfondie : Exercices rsolus

La Couche limite

U 2.

dU
+ U.
.( * + 2. ) = 0

dx
dx

Multipliant les deux termes maintenant par

; on obtient :
.U

U 2
2 dU
*
.
= .
.(2 + ) + 0 .
2 dx
dx

4a- nous avons :

u
u
= (1 ).dy = .
U
U
0
/

1=

u
u
y
(1 ).d ( )
U
U

y m y 2 y
y m y 2
L. + 2 .( ) .1 L. 2 .( ) .d ( )

L 2 1 m
3 L.m 4 m 3 5
2
1= .( ) + .( L ).( )
.( )
.( )

2
3 2
4
20

(C)

En faisant y = , le profil des vitesses nous donne :

1 = L.

m 2
+ .( )
2

(D)

m 2
L
L.m 4 m 3 5

1 m
+ .( ) = 1= .( ) 2 + .( L2 ).( ) 3
.( )
.( )
2

2
3 2
4
20

Posant maintenant : =

1 = L.

L.m 3 m 3 4
1
1 m
L + .( L m). + ( L2 ). 2
.
. = 0
2
3 2
4
20
Ce qui montre que L ne dpend que de m : L = (m)
b- on a :

* = (1
0

u
u
y
y m y
y
)dy = . (1 )d ( ) = . (1 L. .[ ] 2 )d ( )
U
U

0
0

* L 2 m 3
= .( ) .( )
2
6

(E)

36
Mcanique des fluides approfondie : Exercices rsolus

La Couche limite

Entre les quations (C), (D) et (E) ; on peut liminer

*
*
: ce qui implique que
ne dpend

*
que de m
= (m)

c- sur la paroi : y = 0, on a u = v = 0,
Les quations de Prandtl donneront donc :
1 dp
2u

.
+
.
=0

dx
y 2

dp = .U . dU
dx
dx
Lquation (B) donnera :

2 u m.U
= 2
y 2

dU
U
U.
+ .m. 2 = 0
dx

2 dU
m= .
dx

d- nous avons :

0 = .

.
du
U.L
= .
0 = L = ( m)
dy y =0

.U

Lquation (A) peut scrire donc :

U d ( 2 )
.
= m.[2 + (m)] + (m)
2 dx
Il est vident que le deuxime terme ne dpend que de m.
e- lquation prcdente scrit :
U d ( 2 ) 0,45 + 6m
2 dU
U d ( 2 )
=
.
.
= 0,45 + 6m = 0,45 - 6. .
dx
2 dx
2
dx

2 dU
U d ( 2 )
.
+ 6. .
= 0,45.
dx
dx

En multipliant les deux termes par U5 ; nous obtenons :

d
(U 6 . 2 ) = 0,45. .U 5
dx

37
Mcanique des fluides approfondie : Exercices rsolus
transport et distribution des hydrocarbures ::::: universite de boumerdes

transport et distribution des hydrocarbures ::::: universite de boumerdes

La Couche limite
0,45.
=
. U 5 dx
6
U
0
Posons z = k.x
x

Si U = [Link](kx) U 0 .Sin 5 (kx)dx =

5

U 0 kx 5
U z
. Sin (kx)d (kx) = - 0 . Sin 4 z.d (Cosz )
k 0
k 0

U
Cos 5 z 2
= 0 .(Cosz +
.Cos 3 z )
k
5
3
2 =

Cos 5 (kx)
0,45.
2
3
Cos
kx

Cos
kx

.(
.
(
)
(
)
5
k .U 0 .Sin 6 (kx) 3

De mme, on a : m =

m=

dU
2 dU
.
avec
= U 0 .k .Cos (kx)
dx
dx

0,45.
Cos 6 (kx) 2
2
.[
Cos
(
kx
)
+
.Cos 4 (kx)]
5
3
Sin 6 (kx)

Exercice 8 :

Exprimer le coefficient de frottement en fonction de * , et U et cela laide de lquation

de Von Karman ?
Solution :

Lquation de Von-Karman scrit :

( * + 2 ).

p
1 dU d
.
+
=
U dx dx .U 2

Et puisque Cf =

p
1
. .U 2
2

p
1
= Cf
2
2
.U

1 dU d

+
Cf = 2.( * + 2. ). .
U dx dx

Exercice 9 :

Une couche limite se dveloppant sur une plaque plane en coulement uniforme de vitesse u
dont le profil de vitesse est :

38
Mcanique des fluides approfondie : Exercices rsolus

La Couche limite
u ( x, y )
= F( ) = 2 2 3 + 4
U
o : constante.

*
1- montrer que les quantits
,
et p . sont des constantes dont on donnera les valeurs ?

2- dduire de lintgration de lquation de Von Karman, la loi de variation de lpaisseur de
la couche limite (x) ?
3- donner les expressions de (x)/x, * (x)/x et (x)/x ainsi que le coefficient de frottement
U.x
?
p (x)/ U 2 en fonction du seul nombre de Reynolds local x =

Solution :

1-

* ( x) 1
3
A
= [1 f ( )]d =

( x) 0
10 120

2 ( x) dU
Avec : A =
.

dx

* 3
=
10

1
37
37
A
A2
= [1 f ( )]. f ( ).d =

=
0
315 945 9072 315

p . = .

U ( x)
A
.(2 + ). = 2 .U
( x)
6

Qui sont bien des constantes.

p
d
=
dx .U 2
Aprs substitution des expressions prcdentes, on trouve :
2- lquation de Von-Karman se rduit :

.d =

630
. .dx
37 U

En imposant la condition (0) = 0;

Lintgration de cette quation diffrentielle conduit :

( x) =

1260 .x
.x
5,84.
.
37
U
U

En divisant membre membre lexpression ci-dessus par x, il vient immdiatement :

39
Mcanique des fluides approfondie : Exercices rsolus

La Couche limite

( x)
x

Et

5,84
x

* ( x)
x

( x)
x

1,75
x

0,69
x

p
2
U
0,34
=
. .

2
5,84 U .x
.U
x

Exercice 10 :

Soit un coulement caractris par le profil de vitesses suivant :

u ( x, y )
= [Link] ( ) + [Link]( )
U
o =

y
et a, b, et sont quatre paramtres priori fonctions de la seule variable x.
(x)

1- donner les valeurs de a, b, et compatibles avec les conditions aux limites de

lcoulement ?
2- exprimer les valeurs de

*
,
et p . ? vrifier quil sagit bien toujours de constantes

pour cet coulement ?

3- calculer la loi dpaisseur de la couche limite (x) ?
4- comparer les valeurs adimensionnelles (x)/x, * (x)/x, (x)/x et p (x)/ U 2 avec celles
obtenues par le profil de vitesses suivant (exercice prcdant) :
u
u
u
=
=
= 2 2 3 + 4 ?
U U E Ul

5- calculer les paisseurs en bout dune plaque de longueur L = 2 m place dans un

coulement de vitesse U = 20 m/s dun fluide de viscosit cinmatique = 1,4.10-5 m2/s ?
Solution :

1- on imposera :
u(x,0) = 0, V(x,0), u(x, ) = U,

u
y

=0
y =

40
Mcanique des fluides approfondie : Exercices rsolus

La Couche limite

La premire condition amne a = 0

u ( x, y )
= [Link]( . )
U

Les deux dernires conditions donnent alors : [Link] = 1 et b. .Cos = 0

La seule solution non triviale (banale) correspond :

et b = 1

le profil de vitesses de la forme prescrite compatible avec trois conditions aux limites, scrit :

u ( x, y )
= f ( ) = Sin( )
2
U
2- avec cette expression, on obtient :
* 1

2
= 1 Sin( . ).d =

2
0

4
= 1 Sin( . ).Sin( . ).d =
0
2
2
2.

Et

Le domaine de variation physiquement admissible de cette fonction est de 0 1 : [0,1] au lieu

de la variation de 0 .
Et p = .

U
U
du
= . . f ' (0) = . .
2

dy y =0

3- en substituant maintenant dans lquation intgrale, on aura :

.d =

2
. .dx
4 U

Qui avec la mme condition dintgration : (0) = 0 , donne :

( x)
x

= .

2
1
4,8

.
4 x
x

4- la comparaison des rsultats obtenus avec les deux formulations est prsente dans le
tableau suivant, en rfrence aux valeurs de la solution exacte de Blasius :

41
Mcanique des fluides approfondie : Exercices rsolus

La Couche limite

Paramtres ( x)
. x
Profil
x

* ( x)

( x)

. x

p ( x)
. x
.U 2

0,66

0,33

2,64
2,54

. x

Sin( )
2
2. 2. 3 + 4

4,8

x
1,74

5,84

1,75

0,69

0,33

Blasius

4,92

1,72

0,66

0,33

5- avec les valeurs indiques, on trouve :

( L) = 5,7mm
* ( L) = 2,1mm
( L) = 0,8mm

Exercice 11 :

Calculer les valeurs des paisseurs adimensionnelles

* **
,
ainsi que le rapport de forme H

pour le profil de vitesse linaire :

y
u ( x, y )
=
?
U ( x, y )
Solution :

On a :
y
u
=
U

(1 ).dy
0

1 y2
(y .
2
=

1
2

1
2

**
:

Nous avons :

u2
y
y2
y y3
.u
** =
.(1 2 ).dy = .(1 2 ).dy = ( 3 ).dy
.U

0
0
0
=(

2 1 4 1
y2 1 y4
.
=
. 3) =
2 4 0 2 4 3 4

42
Mcanique des fluides approfondie : Exercices rsolus

La Couche limite

** 1
1
= / =

4
4

-H=

Avec * =

*
1
. = .

1
u
y
y
y y2
1 y2 1 y3
.u
=
=
=
=

(
1

).
dy
.(
1

).
dy
(

).
dy
(
.
.
)
2
0 .U U
0
0 2
2 3 0 6

Et =

* 2
H=
=
=3
1

Exercice 12 :

Soit un coulement caractris par le profil de vitesses :

u ( x, y )

= Sin( . )
2
U
O =

y
(x)

- exprimer les valeurs de

*
, / et p . ? Vrifier quils sagit bien toujours de constantes

pour cet coulement ?

- calculer la loi dpaisseur de la couche limite (x) ?
Solution :

a- on a :
u

y
= Sin( . ) avec =
2
U

u
y
= Sin( . )
U
2

u
y

- = (1 )dy = 1 Sin( . ).dy =

U
2
0
0
2
2.
2
+ .Cos
=

2
*

2
y

y + .Cos ( 2 . ) =

43
Mcanique des fluides approfondie : Exercices rsolus

La Couche limite
* / = 1

= 0,363 : constante.

u
u
- = (1 ).dy =
U
U
0

y
y
y
2 y
0 Sin( 2 . ).[1 Sin( 2 . )].dy = 0 Sin( 2 . ).dy - 0 Sin ( 2 . ).dy

1
Et on a : Sin 2 xdx = .( x [Link])
2
1
Sin 2 udu = .(u [Link] )
2
= -

2
.Cos Cos 0
2

y
y
1 y
. .( . Sin . .Cos . ) =
2
2 0
2 2

2 1

2
.( 0) (0 0) =
= .( )
2
2
2

/ = .(

1
2 1
) / = 0,137 : constante
2
2

y
(U .Sin( . )
du
2
=
- P = .
y
dy y =0

P = .U .

= .U .
y=0

.Cos ( . ) y =0 = .U . .Cos 0
2
2
2

2.

P . = . .U : constante.
2

Daprs lquation de Von-Karman :

( * + 2 ).

p
1 dU d
+
=
.
U dx dx .U 2

Et puisque U : constante

p
d
=
dx .U 2

=(

1
).
2

dU
=0
dx

avec / =

1
2

p
d =
.dx
.U 2
44

Mcanique des fluides approfondie : Exercices rsolus

guembour rabah

La Couche limite

Avec P = .U .

.U
2 .dx
2

2 1
d ( ). =
2.
.U
2

.
2 1
dx
=
.dx
( ).d = . .
2
2 U . 2.U .
.
2 1
dx
( ). .d = . .
=
.dx
2
2 U . 2.U
.d =

2.(

1 U
)
2

.dx =

2
. .dx
4 U

1 2
2. 2
2
. =
. .x 2 =
. .x
2
4 U
4 U

( x) =

2. 2
.
.x ( x) = 4,79.
.x
U
4 U
1
.x
x

Ou bien ( x) = 4,79.

Exercice 13 :

Calculer les valeurs de

u y
=
U

* **
,
, H et
pour un profil de vitesses linaire :

1/ 7

Solution :

u y
On a :
=
U

a- * = (1
0

1/ 7

u
y
1
7

)dy = (1 ( )1 / 7 ).dy = y ( )1 / 7 . . y 8 / 7
U

1
7
= =
8
8

* 1

=
8

45
Mcanique des fluides approfondie : Exercices rsolus

La Couche limite

u
u2
u
u
y
y
b- = .(1 2 ).dy = ( ( ) 3 ).dy = [( )1 / 7 ( ) 3 / 7 ].dy
U U
U

U
0
0
0
**

7
1
7
7
1
7
1
1

= ( )1 / 7 . . y 8 / 7 ( ) 3 / 7 . . y 10 / 7 = ( )1 / 7 . . 8 / 7 ( ) 3 / 7 . . 10 / 7
8
10
8

10

7
7
7
. . = .
40
10
8

** 7

40
*
c- H =

Avec =

u
u
u
u 2
y 1/ 7
y 2/7
0 U (1 U ).dy = 0 [U (U ) ].dy = 0 [( ) ( ) ].dy

7
7
1
7
7

7 7
1

= ( )1 / 7 . . y 8 / 7 ( ) 2 / 7 . . y 9 / 7 = . .( )1 / 7 .( ) 2 / 7 = .
9
8
9

8 9

= [

7
63 56
]=
72
72
1

72 36 9
*
H =
= 8 =
=
=
7

56 28 7

72

d- = 72

7
72

Exercice 14 :

Une plaque plane de longueur L est place paralllement un coulement uniforme de vitesse
U linfini.
1- en utilisant les rsultats de la solution de Blasius, tablir lexpression de la trane de la
plaque par unit denvergure en fonction du nombre de Reynolds U.L/ ?
2- en application du thorme dEuler un domaine de fluide que lon prcisera, tablir
lexpression de cette mme trane en fonction de lpaisseur de quantit de mouvement.
Retrouver alors le rsultat de la premire question ?
3- donner lexpression du champ de vitesse longitudinale dans la zone de sillage des couches
limites an aval de la plaque ?
46

Mcanique des fluides approfondie : Exercices rsolus

Ecoulements transitoires

transport et distribution des hydrocarbures ::::: universite de boumerdes

ECOULEMENTS TRANSITOIRES

Exercice 1 :
Donner la clrit de londe de pression provoque par la fermeture rapide dune vanne dans
une conduite, tout en supposant que la paroi de la conduite est rigide ?
Solution :
Lnergie cintique de leau qui va tre transforme en nergie lastique est :
2

M.

V1
1 S .L V1
= ..
.
2
2
g 2

O L : longueur de la conduite
M : masse de leau
V1 : vitesse initiale
= .g
Le module dlasticit cubique de leau est :

p
p
=
volume [Link] V / V0

Avec V =

V .p

S .L. .g.h

Le travail de compression est :

1
1
2
. .S .L.V1 = . .g .h.S .L. .g.h /
2
2
2

h =

V1 .
g 2 .
2

(1)

En ngligeant les frottements, la conservation de la quantit de mouvement donne :

M .V1 ( FX .dt ) = M .V2

Avec Fx : la force qui produit la variation de quantit de mouvement.

.g .h.S = .Q.(0 V1 ) .g .h.S = .S .a.V1
47
Mcanique des fluides approfondie : Exercices rsolus

Ecoulements transitoires

h = a.V1/g

En remplaant dans lquation (1), nous trouvons :

a 2 .V1
V .
= 12
2
g
g .
2

a=

Exercice 2 :

Calculer laugmentation de pression produite par la fermeture instantane dune vanne sur une
conduite ?
Solution :

Soit P est la variation de pression due la fermeture de la vanne.

Daprs lquation de conservation de la quantit de mouvement :
La variation de la quantit de mouvement dans la direction x est :

.Q

.(V2 V1 ) avec = .g
g
Si on nglige les frottements : la force qui produit la variation de la quantit de mouvement
est : P.S
Fx =

- P.S =

.Q
g

.(0 V1 ) = .S .a.(0 V1 )

P = .a.V1 et puisque P = .g.h'

Avec h : la hauteur de pression

P = .g.h' = .a.V1 h =

a.V1
g

Exercice 3 :

Comparer les vitesses des ondes de pression circulant le long dun tuyau rigide contenant :
- de leau 15 C de module cubique B = 2,16.109 Pa. ?
- de lhuile de densit 0,8 et de module cubique 1,38.109 Pa ?
Solution :

Pour leau : a =

2,16.10 9
= 1470 m/s
1000
48

Mcanique des fluides approfondie : Exercices rsolus

Ecoulements transitoires

Pour lhuile : a =

1,38.10 9
= 1310 m/s
0,8.1000

Exercice 4 :

Calculer la clrit de propagation de londe de pression, dans une conduite en acier normal
(E = 2,1.1011 N/m2), o scoule de leau 20 C. ( = 998 kg/m3, = 21,39.108 N/m2) ?
Le diamtre de la conduite est de 0,5 m, et lpaisseur de 6 mm.
Solution :

Nous avons :
a=

1 D
.
E e

a = 998.

1
1
500
+
.
8
11
6
21,39.10
2,1.10

= 1077 m/s

Exercice 5 :

Mme exercice prcdant pour une conduite rigide en PVC (E = 2,5.1010 Pa), de diamtre D =
40 cm, et dpaisseur 2 mm ?
Exercice 6 :

Quelle est laugmentation de pression produite par larrt brutal de la circulation de lhuile, de
densit 0,85 et de module dlasticit = 1,7.109 Pa, raison de 0,57 m3/s dans un tuyau
dacier de 60 cm de diamtre suppos rigide ?
Solution :

Nous avons :
P = .a. u = .a.u0

avec a =

1,7.10 9
= 1414,2 m/s
0,85.1000

u0 = Q/S = Q/[ .D2/4] = 0,57/ [ .0,62/4] = 2,02 m/s

P = .a.u0 = 0,85.1000.1414,2.2,02 = 2428,2.103 Pa

49
Mcanique des fluides approfondie : Exercices rsolus
transport et distribution des hydrocarbures

Ecoulements transitoires
Exercice 7 :

On ferme une vanne brusquement dans un tuyau de 75 mm de diamtre transportant de lhuile

de densit 0,8 et de module dlasticit = 1,38.109 Pa.
Si laugmentation de la pression est de 600 kPa, quel est le dbit probable ? Le tuyau est
suppos rigide.
Solution :

On a :
a=

1,38.10 9
= 1310 m/s
0,8.1000

P = .a. u = .a.u0 = 0,8.1000.1310.u0 .

u0 =

600.10 3
= 0,572 m/s.
0,8.1000.1310

Q = u0. .D 2 / 4 = 5,22.10-3 m3/s.

Exercice 8 :

Si un tuyau dacier de 60 cm de diamtre et de 2440 m de long, a t conu pour rsister une

pression de 103 MPa sous une charge statique maximale de 330 m deau, de combien la
pression dur les parois du tuyau va-t-elle augmenter quand une vanne fermeture rapide
arrte le dbit de 0,85 m3/s ?
On donne le module dlasticit de lacier E = 2.109 Pa.
En supposant que le tuyau est rigide.
Solution :

Nous avons :

P = .a. u

Avec a =
=

2.10 9
= 1,41.103 m/s
1000

Et u = u0 u1 = u0 0 = u0 = Q/S = Q/[ .D2/4] = 0,85/ [ .0,62/4] = 3 m/s

P = 1000.1,41.103.3 = 4,23.106 Pa = 4,23 MPa

Et P = P0 + P
Avec P0 = .g.h = 1000.9,81.330 = 3,24 MPa
50
Mcanique des fluides approfondie : Exercices rsolus

Ecoulements transitoires

P = 4,23 + 3,24 = 7,47 MPa

Exercice 9 :

Un tuyau dacier de 1,2 m de diamtre et de 10 mm dpaisseur, transporte de leau 16 C

( = 2,16.109 Pa) une vitesse de 1,8 m/s. si le tuyau a une longueur de 3000 m, et si la vanne
cot sortie est ferme en 2,50 s, quelle augmentation de pression dans les parois du tuyau
doit-on sattendre ?
On donne pour lacier E = 207.109 Pa.
Solution :

La clrit de londe est :

=
D

a=

2,16.10 9

1 + .
E e

2,16.10 9 1200
10 1 +
.

9
10
207.10

= 980 m/s

Le temps daller et retour de londe est :

T0 = 2L/a = 2.3000/980 = 6,1 s
Puisque T = 2,5 s < T0 = 6,1 s : il sagit donc dun coup de blier direct
P = .a. u = .a.u0 = 1000.980.1,8 = 1,764.106 Pa = 1764 kPa

Exercice 10 :

Dterminer la vitesse de propagation de londe de coup de blier et la surlvation de la

pression en cas dune fermeture instantane dune conduite dacier de diamtre D = 450 mm,
lpaisseur des parois e = 8 mm et la vitesse initiale de leau V0 = 1,8 m/s ?
On donne pour lacier

= 0,01, et E = 207.108 kgf.

Solution :

On a :

a=
1+

avec a0 =
.D

2,07.10 8
= 1425 m/s
102

E.e

51
Mcanique des fluides approfondie : Exercices rsolus

Ecoulements transitoires

a=

1425
450
1+
.0,01
8

= 1137 m/s

Et nous avons : P = .a. h = .g.

a
a
.u0 car h = .(u u 0 )
g
g

P = 1000.1137.1,8 = 2,05.106 Pa = 2050 kPa

Exercice 11 :

Dterminer llvation de la pression une fermeture instantane dune tuyauterie dacier de

diamtre D = 500 mm, parois dune paisseur e = 9 mm ?
La vitesse initiale de leau tant V0 = 1,47 m/s.
Solution :

= 1425 m/s

Nous avons : a0 =

Et a =
1+

=
.D
E.e

1425
500
1+
.0,01
9

= 1142,74 m/s

P = .a. u0 = 1000.1142,74.1,47 = 1,68.103 kPa

Exercice 12 :

Dterminer llvation de la pression lors de fermeture instantane dune conduite de

diamtre D = 600 mm, et dpaisseur des parois e = 5 mm, vhiculant un dbit de 200 m3/h ?
On donne :
E = 18.1010 N/m2.
= 21.108 N/m2.
Solution :

Nous avons :
P = .a. u
Avec :

52
Mcanique des fluides approfondie : Exercices rsolus

Ecoulements transitoires

a=

=
D

1 + .
E e

21.10 8
21.10 8 0,6
10 1 +
.

10
18.10 0,005

= 935,41 m/s

Dautre part, nous avons : Q = 200 m3/h

v = Q/[ .D2/4] = (200/3600)/ [ .0,62/4] = 0,196 m/s
P = 1000.935,41. 0,196 = 183,34 kPa.

Ou bien : P =

u.a 0,196.935,41
=
= 18,69 m.
g
9,81

Exercice 13 :

De lhuile circule dans un tuyau en acier de 2440 m de long raison de 0,57 m3/s, combien
faut-il de temps pour fermer la vanne pour viter la valeur de coup de blier maximale ?
On donne :
La densit de lhuile est 0,85 et son module dlasticit = 1,7.109 Pa.
Le tuyau est suppos rigide de 60 cm de diamtre.
Solution :

Nous avons : T0 = 2L/a

Avec a =

1,7.10 9
= 1414,2 m/s
0,85.1000

Pour viter le coup de blier maximum T > T0

T0 = 2L/a = 2.2440/1414,2 = 3,45 s
Pour viter le coup de blier maximum, le temps de fermeture de la vanne T doit tre
suprieure 3,45 s
T > 3,45 s.

Exercice 14 :

La longueur dune tuyauterie dacier entre le rservoir et la vanne est L = 1800 m, le diamtre
D = 450 mm, lpaisseur e = 6 mm. Le dbit circulant est Q = 127 l/s.
Dterminer llvation maximale de la pression Pmax prs de la vanne sa fermeture
progressive durant T = 3 s et la loi linaire de la variation de la vitesse ?

transport et distribution des hydrocarbures

53
Mcanique des fluides approfondie : Exercices rsolus

Ecoulements transitoires
Solution :

On a : T0 = 2L/a
Avec a0 =

a=

1+

= 1425 m/s

=
.D
E.e

1425
450
1+
.0,01
6

= 1077,18 m/s

T0 = 2L/a = 2.1800/1077,18 = 3,34 s > T = 3 s

Il sagit dun coup de blier direct.

Pmax = .a. u0

Avec : u0 = Q/S = Q/[ .D2/4]

Pmax = .a.Q/[ .D2/4] = 1000.1077,18.0,127/( .0,452/4) = 860,6 kPa

Exercice 15 :

Une conduite de 460 m de long et de 1,2 m de diamtre, vhiculant de leau une vitesse de
0,6 m/s, et la clrit de londe qui la traverse est de 1143 m/s.
Calculer la valeur maximale de coup de blier sur la vanne et 152 m et 304 m de la vanne :
- si la vanne est brusquement et compltement ferme ?
- si la vanne est brusquement et partiellement ferme rduisant la vitesse 0,18 m/s ?

Exercice 16 :

Une conduite de 300 m de long, o leau scoule 0,8 m/s dont la clrit de londe est de
1000 m/s.
Calculer la valeur maximale du coup de blier sur la vanne, et 100 m et 200 m de la
vanne :
- si la vanne est ferme compltement et ?
- si la vanne est brusquement et partiellement ferme permettant la rduction de la vitesse
0,2 m/s ?

Exercice 17 :

De leau scoule une vitesse de 3,048 m/s dans une conduite de 122 m de long, de 0,2 m de
diamtre et de 6 mm dpaisseur.
transport et distribution des hydrocarbures

Mcanique des fluides approfondie : Exercices rsolus

54

Ecoulements transitoires
Calculer le temps dun aller et retour de londe de pression et la valeur maximale de la
surpression produite par la fermeture brusque de la vanne pour :
- une fermeture complte ?
- une rduction de la vitesse 1,83 m/s ?

Exercice 18 :

Nous supposons que la clrit de londe dans une conduite est de 975 m/s. si cette conduite
est de 610 m de long, et de 1,2 m de diamtre, calculer la valeur maximale du coup de blier
sur la vanne, et 910 m du rservoir pour :
- la vanne est partiellement ferme en 4 s permettant le dbit de passer de 0,85 m3/s 0,28
m3/s ?
- le dbit initial est de 0,42 m3/s et la vanne est compltement ferme en 1 s ?

Exercice 19 :

Une conduite en acier de diamtre D = 300 mm, dpaisseur e = 5 mm, de longueur L = 3 km

et comportant une vanne son extrmit, et vhiculant de leau un dbit de 200 l/s.
Si on ferme la vanne partiellement en 3 s permettant de rduire la vitesse 50% de sa valeur
initiale, dterminer laugmentation maximale de la pression dans la conduite ?
On donne :
Le module dlasticit volumique de leau = 21.108 N/m2.
Le module dlasticit longitudinale de lacier E = 2.1011 N/m2.
Solution :

Nous avons :
a=

=
D

1 + .
E e

21.10 8
21.10 8 0,3
10 1 +
.

2.1011 0,005

= 1135,05 m/s

Le temps daller et retour de londe est :

T0 = 2L/a = 2.3000/1135,05 = 5,286 s
Donc T0 = 5,286 s > T = 3 s
Il sagit dun coup de blier direct.
Pmax = .a. u = .a. (v0 v) = 1000.1135,05.(v0 v0/2) = 1000.1135,05.v0/2

Avec v0 = Q0/[ .D2/4]

Pmax = 1000.1135,05.(Q0/2)/[ .D2/4] = 1000.1135,05.(0,2/2)/[ .0,32/4] = 1606,58 kPa

Pmax = 1606,58 kPa

55
Mcanique des fluides approfondie : Exercices rsolus
transport et distribution des hydrocarbures

Ecoulements transitoires
1
a. Q0 / S
a.Q0
1135,05.0,2
u.a
=
=
= 163,77 m
= 2
Ou bien hmax =
g
g
2 g.S
0,3 2
2.9,81. .
4

Exercice 20 :

Une conduite en acier de longueur L = 2 km, et de diamtre D = 300 mm, reliant un rservoir
et une vanne, et vhiculant un dbit Q = 100 l/s.
Dterminer llvation maximale de la pression sur la vanne, si elle se ferme progressivement
pendant 4 s ?
Sachant que la loi de variation de la vitesse est linaire.
On prend :
e = 5 mm.
= 21.108 N/m2.
E = 2.1011 N/m2.
Solution :

Nous avons :
a=

=
D

1 + .
E e

21.10 8
21.10 8 0,3
10 1 +
.

2.1011 0,005

= 1135,05 m/s

T0 = 2.2000/1135,05 = 3,52 s
Donc T0 = 3,52 s < T = 4 s.
Il sagit donc dun coup de blier indirect.
Pmax =2. .L.u0 /T = 2. .(L/T).(Q0 /S) =

2.10 3.2.10 3
.
4

0,1
= 1415,43 kPa
0,3 2
.
4

Pmax = 1415,43 kPa

Ou bien :
2 L.u 0 2 L.Q0 / S
hmax =
=
=
g.T
g.T

2.2.10 3.0,1
= 144,28 m
0,3 2
9,81.4. .
4

Exercice 21 :

Calculer laugmentation de pression produite par le fermeture partielle dune vanne monte
sur une conduite en acier de module dlasticit E = 205.109 Pa, dpaisseur e = 2 mm, de
diamtre D = 70 cm et de longueur L = 2 km, vhiculant de leau 20 C de module
dlasticit = 21,39.108 Pa un dbit Q0 = 0,8 m3/s, sachant que la fermeture partielle de la
vanne conduit une rduction de 50% de la vitesse pendant 4 s ; et cela pour les deux cas :

56
Mcanique des fluides approfondie : Exercices rsolus

Ecoulements transitoires
a- la conduite est rigide ?
b- la conduite est lastique ?
Solution :

a- conduite rgide :
La clrit de londe est :
a=

21,39.10 8
= 1,46.103 = 1462,5 m/s
1000

Le temps de laller et le retour de londe sera :

T0 = 2.L/a = 2.2.103 /1462,5 = 2,73 s
Le temps de fermeture est T = 4 s
Donc : T0 = 2,73 s < T = 4 s
Il sagit dun coup de blier indirect
Pmax =

2 .L
2L 1
2 L 1 Q0
=
. .u 0 =
. .
.(u u 0 ) =
g .T
g .T 2
g .T 2 S

Pmax = 106 m 1,06 MPa 10,6 bars

2.10 3.0,8
= 106 m
.0,7 2
9,81.4.
4

b- conduite lastique :
La clrit de londe :
a=

=
D

1 + .
E e

21,39.10 8
21,39.10 8 700
10 1 +
.

2
205.10 9

= 678,2 m/s

Le temps de laller et le retour de londe sera :

T0 = 2L/a = 2.2.103/278,2 = 5,9 s
Nous avons T0 = 5,9 s > T = 4 s
Il sagit dun coup de blier direct :
1 Q
a. . 0
Pmax = a.(u - u0)/g = 2 S =
g

a.Q0
678,2.0,8
=
= 71,89 m
2
D
0,7 2
2 g. .
2.9,81. .
4
4

Pmax = 71,89 m 0,72 MPa 7,19 bars

57

Mcanique des fluides approfondie : Exercices rsolus

Ecoulements transitoires
Exercice 22 :

Une conduite lvatoire de longueur L = 500 m et de section s = 1 m2, o scoule un dbit Q0

= 1 m3/s, est protge dans sa partie initiale par une chemine de section constante = 5 m2.
Calculer loscillation dans la chemine, dans le cas dun arrt brusque des pompes, et compte
tenu de la perte de charge dans la conduite quest de H = 4,5 m ?

Exercice 23 :

Une galerie damnagement hydrolectrique de longueur L et de section s, relie une retenue

de superficie assez grande dont les variations de niveau sont ngligeables ; une chemine
dquilibre verticale cylindrique de diamtre D.
En ngligeant les pertes de charge, trouver lexpression de la priode des petites oscillations
de lensemble ?
D (S)
h

d (s)

L = 16km

On donne s = 16m 2
D = 16m

Solution :

La masse de leau dans la galerie est : m = .s.L

(1)

Les forces dinertie aux quelles soumise cette masse : m.

Les forces de pression : S. .g.h
Lquation de continuit : s.V = S.

dV
dV
= .s.L.
dt
dt

(2)

dh
dt

Avec :
S : section de la chemine
dV S d 2 h

= . 2
dt
s dt

(3)

58

Mcanique des fluides approfondie : Exercices rsolus

Ecoulements transitoires
En remplaant (1) dans (2), et faisant galit des forces de pression aux forces dinertie, on
trouve :
.s.L.

dV
= S . .g .h
dt

(4)

Remplaant maintenant (3) dans (4) pour retrouver :

L.S d 2 h
.
+h=0
g .s dt 2

La solution oscillante non amortie de cette quation diffrentielle est de type :

h = [Link]( 2.

t
)
T

16.10 3. .16 2
L.S
= 2. .
= 891 s
O la priode est : T = 2. .
g .s
10.4.16
h0 = m Q0 .

L
L.s
= m v0 .
s.g.S
g.S

Exercice 24 :

On suppose que la galerie de lexercice prcdant dbite Q = 48 m3/s.

1- si lon ferme rapidement la conduite immdiatement en aval de la chemine dquilibre,
quelle hauteur leau montrait-elle dans la chemine dquilibre ? (en ngligeant les
frottements).
2- si la hauteur de la chemine dquilibre est gale 20 m, quel sera le volume deau
dvers ?
Solution :

1- Daprs le mme raisonnement prcdant, nous avons :

L.S d 2 h
.
+h=0
(1)
g .s dt 2
Dont lintgration de cette quation, donne :
h = [Link]( 2.

Avec T = 2. .

t
)
T
L.S
g.s

59
Mcanique des fluides approfondie : Exercices rsolus
guembour rabah

Ecoulements transitoires
La valeur de h0 qui correspond au maximum de hauteur dans la chemine, sera obtenue en
considrant les conditions initiales.
On a : s.V0 = S.

dh
s
dh 2
t
.V0 =
=
.h0 .Cos (2 . )
dt
S
dt
T
T

Pour t0 = 0 Cos(

2
s
2
.0 ) = 1 V0 =
h0
T
S
T

L.s
16.10 3.16
h0 = V0 .
= 3.
= 33,9m
g .S
.16 2
10.
4

2- lquation prcdente (n1) nest valable que pour des niveaux ne dpassant pas la hauteur
de la chemine (20 m), si non ; on considre que la force de pression reste constante :

.S .L.

dV
L dV
+ .g.S .h1 = 0 .
+ h1 = 0
dt
g dt

Avec h1 = 20 m.
V = V1 -

g .h1 .t
L

(2)

O V1 : vitesse dans la conduite linstant o h atteint 20 m dans la chemine, le tant tant

compt maintenant partir de cet instant.
T1

Le volume dvers U =

s.V .dt = s.V .T

1

s.g .h1 2
.T1
2L

O T1 est le temps ncessaire pour annuler la vitesse dans la galerie.

Daprs (2) : T1 =
2

V1 .L
g .h1
2

s.V1 .L s.g .h1 V1 .L

s.V .L
U=
. 2 2 = 1

g .h1
2 L g .h1
2 g .h1

La vitesse dans la premire phase de lcoulement est :

V=

2 .t1
S dh S
2
2 .t
.
= .h0 . Cos (
) = V0 .Cos (
)
s dt s
T
T
T

O t1 : temps coul entre la fermeture et le dbut de dversement :

60
Mcanique des fluides approfondie : Exercices rsolus

Ecoulements transitoires
on a : Sin( 2.

t1
) = h1/h0 = 20/33,9 = 0,59
T

t1 = 98,48 s

V1 = v0.0,807 = 2,42 m/s

U = 3740 m3.
T1 = 194 s.
V = 3m / s
Avec bien sure : 0
T = 981s

Exercice 25 :

On dsire fermer un orifice de section s0, situ lextrmit avale dune conduite longue, de
section S, de telle sorte que le dbit dcroisse linairement en fonction du temps.
1- dterminer la loi de pression devant lorifice pendant la fermeture, en supposant que le
coefficient de dbit m est constant ?
s
2- quelle doit tre la loi de fermeture de lorifice en fonction du temps :
?
s0
On admettra que la dure T de fermeture est infrieure au temps que mettent les ondes de
pression mises depuis lorifice, pour y revenir aprs rflexion sur lextrmit amont de la
conduite ?
Solution :

Le dbit traversant lorifice chaque instant :

qv = m.s. 2 g.h
Avec h : charge devant lorifice
La variation temporelle de ce dbit est :
qv = qv0 qv0.t/T dqv = - [Link]/T
La surpression devant lorifice est:
dh =

a
a
du =
dq v
g
g .S

dh =

a
.q v 0 .dt
g .S .T

h0

dh =

a qv 0
. .dt
g.S T

61
Mcanique des fluides approfondie : Exercices rsolus
guembour rabah

Ecoulements transitoires

h=

a.q v 0
.t + h0
g .S .T

O h0 : charge initiale devant lorifice.

En fin de fermeture : t = T h1 =

a.u 0
a.q v 0
+ h0 =
+ h0
g
g.S

O u0 : vitesse initiale dans la conduite.

La surpression totale produite par la fermeture est :

h = h1 h0 =

a.u 0 a.q v 0
=
g
g.S

2- Nous avons :
qv = m.s. 2 g.h
2

h=

qv
qv 0
t
=
.(1 ) 2
2 2
2 2
T
m .s .2 g 2 g.m .s

Avec h = h0 +

a.q v 0 t
.
g.S T
2

qv 0
a.q v 0 t
t
.(1 ) 2
h0 +
. =
2 2
T
g.S T 2 g.m .s
2

A t = 0, nous avons : qv0 = m.s0. 2 g .h0

qv0
et h0 =
2 g.m 2 .s 2

(1 ) 2
2

a.q v 0 .t
qv 0
1
T

=
.
2
2
2
g .S .T 2 g .m s
s0

t
s
T
=
La loi de fermeture est donne par :
2
s0
2.a.m 2 .s 0 t
1+
.
q v 0 .S
T
On peut crire aussi :
t
t
1
1
s
T
T
=
=
s0
a.u 0 t
h t
1+
.
1+
.
h0 T
g.h0 T
1

62
Mcanique des fluides approfondie : Exercices rsolus

Ecoulements transitoires
Exercice 26 :

Un manomtre est constitu par deux tubes verticaux de diamtre D, relis par un tube plus
fin de longueur L et de diamtre d.
En admettant que les pertes de charge dans ce dernier tube sont chaque instant, gales
celles dun coulement laminaire permanent, calculer la condition pour laquelle, la viscosit
du liquide ralise lamortissement critique ?
D

L
Solution :

La masse du fluide contenue dans le tube horizontal est :

m = .L.

.d 2
4

Sa force dinertie chaque instant t est :

dV
d 2 dV
m.
= .L. . .
dt
4 dt
Cette force est quilibre par des forces de pression et de frottement.
.d 2
- les forces de pression sont :
. .g .h
4
.d 2
L V2
avec P = f . . .
- les forces de frottement sont : P.
4
2
d
64 64.
En coulement laminaire f =
=
V .d
P.

.d 2
4

= 8 . .L. .V

En quilibre ; la somme des forces est nulle :

63
Mcanique des fluides approfondie : Exercices rsolus

Ecoulements transitoires
.L.

.d 2 dV
4

dt

+ .g .

.d 2
4

.h + 8 . .L. .V = 0

(1)

Dautre part, lquation de continuit des volumes fluides est :

V.

.d 2
4

1 D 2 dh
D 2 dh
. V =
.
= . .
2
4 dt
2.d 2 dt

dV
D2 d 2h
=
.
dt
2d 2 dt 2

(2)

On remplace (2) dans (1), et on simplifie, nous trouvons :

1 D2 d 2h
L. .D 2 dh
. 2 . 2 + 16.
. +h=0
2 g d dt
g.d 4 dt
Cette dernire, est une quation diffrentielle linaire du second ordre coefficients constants,
de type :
A.

d 2h
dh
+ B. + C.h = 0
2
dt
dt
2

O lamortissement critique est obtenu avec : B = 4.A.C

L. 2 .D 2
128.
=1
g.d 6
L. 2 .D 2
Et le systme sera oscillatoire si 128.
<1
g.d 6
Application numrique : avec g = 9,81 m/s2.
= 10-6 m2/s , d = D = 0,01 m L = 780 m.
= 10-6 m2/s , L = 10 m d = D = 3,4 m.
= 10-6 m2/s , D = 0,01 m , d = 0,032 m L = 0,84 m.

Exercice 27 :

On opte pour une chemine dquilibre pour protger une conduite en fonte de 800 m de long,
et de 1,5 m de diamtre, contre le coup de blier. La chemine est place sur la conduite 35
m au dessous de la surface libre du chteau deau.
Pour un dbit de 12 m3/s dans la conduite, quel est le diamtre minimal donner au chemine
pour viter de dpasser 18,5 m au dessus de la surface libre du chteau deau ?
On donne le coefficient de perte de charge dans la conduite C = 11,1 ; et on nglige les effets
de linertie et les pertes dans la chemine.

64
Mcanique des fluides approfondie : Exercices rsolus

Ecoulements transitoires
zmax
D (A)
z=0
z0
d (S)

Solution :

Normalement le niveau dans la chemine se stabilise z = 0 (mme niveau que celui dans le
chteau deau) mais cause des pertes de charge dans la conduite, le niveau se stabilise z0
(infrieur z).
Nous avons :
2

V
-z0 = C. 0 avec C : coefficient de perte de charge.
2g
Pour une vitesse V donne (pendant une fermeture dune vanne) :
La cote deau dans la chemine sera :
Faisant forces dinertie gales aux forces de pression :
m.

dV
= S . .g.z + S . .g.z
dt
2

V
V
Et puisque m = .S .L et z = z z 0 = 0 C. 0 = C. 0
2g
2g
Nous remplaons et nous simplifions, nous trouvons :
.S .L.

dV
V2
= S . .g .z + S . .g.C.
dt
2g
2

z + C.

V0
L dV
= .
2g
g dt

Si on donne
S : section de la conduite
A : section de la chemine
De lquation de continuit, on crit :

65
Mcanique des fluides approfondie : Exercices rsolus

Ecoulements transitoires

S.V = A.

dz
A dz
V= .
dt
S dt
2

V
L S
dV
-z = C. 0 + . .V .
2g
g A
dz
2

z + C.

V0
L A d 2z
=- . . 2
g S dt
2g
2

V0
C A 2 dz
L A d 2z
. .
. . 2 + z = - C.
=2g
2 g S 2 dt
g S dt

V
L S
dV
Ou bien z = C. 0 + . .V .
2g
g A
dz
Lintgration de cette quation donne :
2

V =

2 g L.S

.
z C1 .e
C C. A

C . A
.z
L. S

(1)

Avec C1 : coefficient dintgration.

Pour les conditions aux limites, nous avons :
t = t0 : V = V0 et z = z0.
2 g L.S

V =
.
z 0 C1 .e
C C. A

C2.A

.V0
2 g.L.S
2
= C1 .e 2 g .L.S
C .A

C . A
. z0
L.S

C.A
2
.C . 0
V0
2 g L.S
L. S
2g
=
.
+ C.

C
.
e

1
C C. A
2 g
V

(2)

Et pour z = zmax : V = V0.

En remplaant dans lquation (1) :
0=

2 g L.S

.
z max C1 .e
C C. A

2 g L.S

.
z max = C1 .e

C C. A

C . A
. z max
L. S

C . A
. z max
L.S

2 g L.S 2 g
.

.z max = C1 .e
C
C2 A

C . A
. z max
L.S

(3)

66
Mcanique des fluides approfondie : Exercices rsolus

Ecoulements transitoires
De lquation (2), on tire :
C2.A

2 g L.S 2 g .L.S .V0

.
.e
C1 =
C2 A
Nous remplaons maintenant dans (3) :
2

C2.A

2 g L.S 2 g
2 g L.S 2 g .L.S .V0
2.

.z max = 2 .
.e
A
A
C
C
C
C. A

C. A
. z max
L.S

V0 2

.( C .
+ z max )
C. A
1.z max = e L.S 2 g
L.S

(4)

Et puisque le 2eme terme est toujours positif :

C. A
.z max < 1
L.S

Application numrique :
12
V0 = Q0/S =
= 6,79 m/s
.1,5 2 / 4
V2
6,79 2
= -11,1.
= -26,09 m : 26,09 < 35
z0 = -C.
2g
2.9,81

zmax = 18,5 m
L.S
800. .1,5 2 / 4
C. A
= 6,88
.z max < 1 A <
=
On a
C.z max
11,1.18,5
L.S
Rsolvant maintenant lquation (4) par itrations successives :
A
6,88
6,80
6,00
6,06

C. A
1.z max
L.S
0,00064
0,01226
0,12847
0,11975

V 2
C.A
.( C . 0 + z max )
L.S
2g

e
0,08994
0,09249
0,12239
0,11985

A 6,06 m2 d = 2,78 m.

Exercice 28 :

Une conduite dadduction deau de longueur L = 2 km et de diamtre D = 2 m, est protge

contre les phnomnes transitoires par une chemine dquilibre, place 30 m sous le niveau
de la surface libre du rservoir (voir le schma en dessous).
67
Mcanique des fluides approfondie : Exercices rsolus

Ecoulements transitoires
Trouver le diamtre minimal donner au chemine pour viter que leau ne dpasse pas les
20 m au dessus de la surface libre du rservoir ?
On donne le coefficient des pertes de charge dans la conduite C = 10,2 ; et en ngligeant les
effets de linertie, et les pertes de charge dans la chemine.
Le dbit vhicul par la conduite est Q = 15 m3/s.
zmax = 20 m
d
z=0
z0
30 m

D=2 m

L=2m
Solution :

Nous avons :
C. A
1
.z max = e
L.S
1>

v 2
C . A
.( C . 0 + z max )
L.S
2g

C. A
.z max
L.s

Avec v0 = Q0/s =

15
2
.2

= 4,78 m/s
4

zo = -C.v2/2g = -10,2.4,782/2.9,81 = -11,88 m

O 11,88 < 20 m

Et 1 >

C. A
L.S
.z max A < 1 >
C.z max
L.s

22
4 = 30,78 m2.
=
10,2.20
2000. .

1-0,0325.A = e-0,0518.A

Cest une quation transcendante, pour la rsoudre, on utilise les itrations successives :

68
Mcanique des fluides approfondie : Exercices rsolus

Ecoulements transitoires

A
30,78
30
20
19,8
19,7
19,6
19,65

1-0,0325.A
-0,00035
0,025
0,35
0,3565
0,35975
0,363
0,36137

e-0,0518.A
0,203
0,2114
0,3548
0,3585
0,3604
0,3623
0,36136

Donc A = 19,65 m2 = .d 2 / 4
d=

4.19,65

= 5,0032 5 m.

Exercice 29 :

Une chemine dquilibre de forme cylindrique, de diamtre D = 5 m, est choisie comme

moyen de protection anti-blier dune conduite dadduction relie un chteau deau.
1- si la conduite est de diamtre d = 1,5 m et de longueur L = 12 km, vhiculant un dbit Q0 =
1 m3/s, en ngligeant les pertes de charge, calculer laugmentation du niveau deau dans la
chemine aprs une fermeture rapide de la conduite ?
2- si le niveau deau dans le rservoir est H = 10 m, donner la hauteur totale de la chemine ?
D
Rservoir
Chemine
H
d

L
Solution :

Il sagit dune oscillation en masse de la masse deau existante dans la conduite, entre le
rservoir et la chemine.
- La masse deau dans la conduite est : m = .s.L
- les forces dinertie aux quelles est soumise cette masse :
dV
dV
m.
= .s.L.
dt
dt
- les forces de pression exerces sur cette masse :
69

Mcanique des fluides approfondie : Exercices rsolus

Ecoulements transitoires

s. .g.h

dV
(1)
dt
Daprs lquation de continuit :
s. .g.h = .s.L.

dh
dV S d 2 h

= . 2
s.V = S.
dt
dt
s dt

L.S d 2 h
.
+h=0
g .s dt 2

Et lquation (1) devient :

S d 2h
s. .g.h = .s.L. . 2
s dt
Pour que la masse oscille entre la chemine et le rservoir, la solution de la dernire quation
est de type :
h = [Link](

2t
)
T

O la priode des oscillations est T = 2 .

Et leur amplitude est h0 = m Q0 .

h0 = m 1.

12.10 3

L
L
= m Q0 .
2
s.g.S

.g.d 2 .D 2
2
4

= 5,94 m
2

.9,81.1,5 .5

42

L S
.
g s

h = 5,94 m.

La hauteur totale de la chemine est h = H + h = 10 + 5,94 = 15,94 16 m.

Exercice 30 :

Une conduite de diamtre D = 2 m et de longueur L = 20 km, vhiculant un dbit Q = 3 m3/s,

reliant un rservoir de grande section et une chemine dquilibre verticale de diamtre d =
0,5 m.

70

Mcanique des fluides approfondie : Exercices rsolus

transport et distribution des hydrocarbures ::::: universite de boumerdes

transport et distribution des hydrocarbures ::::: universite de boumerdes

Ecoulements transitoires

Si on ferme rapidement la vanne qui se trouve juste aprs la chemine, en ngligeant les
pertes de charge :
- calculer la priode des petites oscillations de lensemble ?
- calculer llvation maximale de leau dans la chemine ?
D
Hmax

L
Solution :

Nous avons :
T = 2. .

H0 = m Q0

L.S
20.10 3. .0,5 2 / 4
= 2. .
= 70,89 s
g .s
9,81. .2 2 / 4
L
= m3
s.g .S

Ou bien h0 = m v0.

20.10 3
22
0,5 2
. .9,81. .
4
4

= m 172,56 m

L.s
= m 172,56 m
g.S

transport et distribution des hydrocarbures

Mcanique des fluides approfondie : Exercices rsolus

71
guembour rabah

Rfrences bibliographiques

Rfrences bibliographiques
1 [Link], Mcanique des fluides : lments dun premier parcours, Edition Cpadus,
Toulouse 1997.
2 [Link], Mcanique des milieux continus : concepts de base, Edition Dunod, Paris 1997.
3 [Link], Problmes rsolus de mcanique des fluides, Edition Dunod, Paris 1995.
4 [Link], Mcanique des fluides, Edition Marketing, Paris 1988.
5 [Link] et [Link], Physique thorique : mcanique des fluides, Edition Mir, Moscou,
1994.
6 [Link], Mcanique des fluides, Edition OPU, Alger 1994.
7 [Link], Mcanique des fluides, Edition Ellipses, Paris 1991.
8 [Link] et [Link], Exercices de mcanique des fluides : tome 1 et 2, Edition
Eyrolles, Paris 1994.
9 [Link], Mcanique exprimentale des fluides: tome 1,2,et 3, Edition Masson, Paris
1994.
10 Researche & education association, Problem solvers: fluid mecanics, USA 1996.
11 [Link], Dynamique des fluides, Presses polytechniques modernes, Lausanne 1985.
12 [Link], Fluides en coulement, Edition Masson, Paris 1991.
13 [Link], Modlisation asymptotique en mcanique des fluides Newtoniens, Edition
Springer-Verlag, Villeneuve dAscq (France) 1994.
14 [Link] et [Link], Nonlinear instability, chaos and turbulence, Edition WIT press,
Canada 1998.
15 [Link], Exercices resoles de mcanique des fluides, Edition El-Oumma, Alger 1999.
16 [Link], Mcanique et rhologie des fluides, Edition Lavoisier, Paris 1985.
17 [Link] et [Link], Fluid mecanics, Edition McGraw-Hill companies, USA
1997.
18 [Link] and al., Experimentation, modelling and computation in flow,
turbulence and combustion: tome 2, Edition John wiley & sons Lrd. Chichester (England)
1996.
19 [Link] et [Link]-zade, Hydraulique gnrale, Edition OPU, Alger 1993.
20 [Link], Hydraulic design handbook, Edition McGraw-Hill companies, USA 1999.
21 [Link], Mcanique des fluides, Edition Dunod, Paris 1996.
22 Edmond [Link], Andr Martinot-Lagarde et Jean Mathieu, Mcanique des fluides, Edition
Dunod, Paris 1970.

transport et distribution des hydrocarbures ::::: universite de boumerdes

transport et distribution des hydrocarbures ::::: universite de boumerdes
rabah

guembour

transport et distribution des hydrocarbures ::::: universite de boumerdes

transport et distribution des hydrocarbures

72
Mcanique des fluides approfondie : Exercices rsolus

guembour rabah
transport et distribution des hydrocarbures ::::: universite de boumerdes