Mécanique Des Fluides 2 Eme Année

Télécharger au format pdf ou txt
Télécharger au format pdf ou txt
Vous êtes sur la page 1sur 39

Mécanique des Fluides

2e année

Mesures Physiques
Mécanique des Fluides

STATIQUE DES FLUIDES


I Equations
I1 Equation intégrale
I2 Equation locale
II Fluide homogène incompressible dans un repère terrestre (galliléen approché)
II1 Equation locale
II2 Equation des isobares et de la surface libre
II3 Equation intégrale – Théorème d’ARCHIMEDE
III Fluide homogène compressible dans un repère terrestre (galliléen approché)
III1 Equation locale : cas d’un liquide, cas d’un gaz

DYNAMIQUE DES FLUIDES


I Coordonnées de LAGRANGE et coordonnées d’EULER
II Vitesse
III Dérivée particulaire
IV Accélération
V Champ des vitesses dans un écoulement permanent (ou stationnaire)
VI Ecoulement irrotationnel
VII Equation locale de conservation de la masse
VIII Equation intégrale de conservation de la masse. Débit de liquide à travers une section de passage
VIII1 Equation intégrale
VIII2 Débit à travers la section droite d’un tuyau
IX Principe fondamental de la dynamique appliqué à un fluide parfait – Equation d’EULER
IX1 Théorème de la résultante en fluide parfait
IX2 Intégrale première du mouvement stationnaire. Equation de Daniel BERNOULLI (1700-1782)
IX3 Equation de BERNOULLI
X Aspect énergétique de l’équation de BERNOULLI. Energie mécanique totale d’une particule fluide
XI L’équation de BERNOULLI en pratique
XII Pression motrice dans un écoulement unidirectionnel stationnaire

XIII Ecoulement d’un fluide visqueux newtonien incompressible


XIII1 Ecoulement interne et écoulement externe
XIII2 Ecoulement laminaire et écoulement turbulent
XIII3 Ecoulement laminaire établi entre deux plaques parallèles infinies
XIII4 Ecoulement laminaire développé établi dans une conduite de section circulaire
XIV Ecoulement dans les conduites en charge
XIV1 Profils des vitesses dans les tuyaux
XIV2 Répartition des contraintes de cisaillement (ou frottement) dans les écoulements établis (cas des conduites)
XIV3 Bilan d’énergie dans un tube : PERTE DE PRESSION TOTALE et PERTE DE CHARGE
XIV4 Perte de charge régulière et perte de charge singulière
XV Calcul des pertes de charge
XV1 Coefficient de perte de charge régulière
XV1 1 Cas de l’écoulement laminaire dans une conduite de section circulaire : Loi de POISEUILLE
XV1 2 Cas des conduites cylindriques de section non circulaire
XV1 3 Cas de l’écoulement turbulent établi
XV2 Coefficient de perte de charge singulière

XVI Quelques éléments sur la théorie de la couche limite


XVI1 Introduction
XVI2 Equation intégrale
XVI3 Ecoulement dans la couche limite d’une plaque plane
XVI4 Equation intégrale pour un écoulement sans gradient de pression
XVI5 Coefficients de frottement des plaques lisses en écoulement turbulents

XVII MESURES DANS LES FLUIDES


XVII 1Mesure de pression
XVII 2 Mesure de vitesses
XVII 3 Mesure de température
XVII 4 Mesure de débit

Mécanique_des_fluides_2A_MP.doc 2 21/11/2007
STATIQUE DES FLUIDES

I- Equations
I-1 Equation intégrale
Soit un domaine D de fluide, de surface S, et au repos dans le repère d’étude R (x,y,z). En statique des fluides,
r
on choisit en général la direction z, (de vecteur unitaire z ) confondue avec la verticale du lieu. Dans ce cas,
r
l’orientation de z sera toujours « vers le haut » (ascendante).

Repère de référence R
et domaine fluide D n(M*)
Direction z M*
M(x,y,z)

z
D
y Direction y

x
O
n
ct i o
Dir
e x

r
M désignera un point dans D : M ∈ D ; n ( M ) est le vecteur unitaire normal à S au point M.
r
M* désignera un point de S : M * ∈ S ; n (M *) est le vecteur unitaire normal à S au point M*.
Le fluide extérieur à D est appelé complémentaire de D, et noté D .

Ce domaine fluide est soumis à :

Des forces à distance (comme la gravité)


r r
df g ( M ) = ρ( M ) F ( M )dω (I-1)

r
de densité massique F (M ) exercée sur chaque particule élémentaire (au sens de la MMC) de masse ρ( M )dω .
Ces forces sont aussi appelées «forces de volume », c-à-d réparties sur tout le volume de fluide contenu dans D.

Des forces de pression exercées par le fluide extérieur à D. Ces «forces de surface » sont exercées par le com-
plémentaire D de D en chaque point M* de la surface S et s’ écrivent, en chaque point M* :
r r
df ( M *) = − p ( M *) n ( M *) dA (I-2)

r
où on rappelle que n (M *) est la normale extérieure à D en M* c-à-d dirigée vers le fluide de D .

La condition d’équilibre de D s’exprime à l’aide des deux torseurs d’actions extérieures :

T ( F / D) + T ( D / D) = 0 (I-3)

où T ( F / D) est le torseur des forces à distance, et, T ( D / D) est le torseur des forces de pression.

Mécanique_des_fluides_2A_MP.doc 3 21/11/2007
Les éléments de réduction de ces deux torseurs (résultante et moment) sont définis par :

r r


T ( F / D)=  r
R ( F / D ) = ∫ ρ
M ∈D
r
( M ) F ( M ) dω
r
O

M (O, F / D) = OM ∧ ρ( M ) F ( M )dω
M ∈D
r r



R( D / D) = ∫ − p( M *) n ( M *) dA
M *∈S
T ( D / D) =  r r r
 ∫
M (O, D / D) = OM ∧ − p( M *) n ( M *) dA
M *∈S

L’équation d’équilibre entre les résultantes s’écrit :

r r
∫ ρ(M ) F (M )dω + ∫ − p(M *)n(M *)dA = 0
M ∈D M *∈S
(I-4)

I-2 Equation locale

Décomposons le domaine D en « particules » élémentaires parallélépipédiques (« briques ») dont chaque face est
perpendiculaire à l’une des directions du repère R, et dont le centre est le point M. Le volume d’une « particule »
s’écrit :

dω = dx dy dz (I-5)
r
Les forces extérieures appliquées à une « particule » sont les forces à distance, ρ( M ) F ( M ) dx dy dz , et les for-
ces de pression sur chaque face du parallélépipède dont nous cherchons les expressions ci-dessous.
r
Si l’on appelle p(M ) la pression en M, on a pour la face orientée par x et située à l’abscisse x + dx :

 1 ∂p ( M )  r
−  p( M ) + dx  x dydz
 2 ∂x 

et pour la face opposée située à l’abscisse x :

 1 ∂p ( M )  r
 p( M ) − 2 ∂x dx  x dydz .
 
r
La force de pression suivant la direction x est la somme des deux expressions précédentes, soit :

∂p( M ) r
- dx dy dz x (I-6)
∂x

On obtient le même type d’expressions pour les quatre autres faces en permutant les vecteurs unitaires et les
coordonnées.

La condition d’équilibre s’écrit donc, pour la particule considérée :

 r ∂p( M ) r ∂p ( M ) r ∂p ( M ) r
ρ( M ) F ( M ) − x− y− z  dx dy dz = 0 (I-7)
 ∂x ∂y ∂x 

dω = dxdydz ne pouvant être physiquement nul, on doit avoir en tout point M de D:

Mécanique_des_fluides_2A_MP.doc 4 21/11/2007
r  ∂p( M ) r ∂p( M ) v ∂p( M ) r
ρ( M ) F ( M ) −  x+ y+ z = 0
 ∂x ∂y ∂z 

On appelle gradient de pression en M, le vecteur défini par :

∂p ( M ) ∂p( M ) r ∂p( M ) r
grad p( M ) = x+ y+ z (I-8)
∂x ∂y ∂z

D’où l’équation locale (c-à-d au point M) :


r
ρ( M ) F ( M ) − grad p ( M ) = 0 (I-9)

Si ρ(M ) est constante en tout point M de D, ou si ρ(M ) ne dépend que de la pression, l’équilibre du fluide n’est
r
possible que si F (M ) dérive d’un potentiel scalaire U tel que :

r
F ( M ) = −grad U ( M ) (I-10)

(U sera précisé ultérieurement) d’où l’équation fondamentale de la statique des fluides:

− ρ( M )grad U ( M ) − grad p ( M ) = grad[ p( M ) + ρ( M )U ( M )] = 0 ⇒ p( M ) + ρ( M )U ( M ) = Constante, ∀M ∈ D


(I-11)

Les surfaces équipotentielles ( U = Cte ) et isobares ( p = Cte) sont confondues.


r
En chaque point M, la normale commune a la direction de F (M ) .

Remarque 1 : Résultante des forces de pression appliquées sur les faces d’un élément de volume. Quelle
que soit la situation (fluide au repos ou en mouvement) la somme vectorielle des forces de pression appliquées
sur les 6 faces de l’élément s’exprime avec le vecteur gradient conformément à l’équation I-8:
r
[
dF(M) = − grad p(M) dω ]
Le gradient de pression est déterminé physiquement par l’équilibre des actions mécaniques existant dans le
r
milieu fluide. Ainsi, en statique, il est égal au vecteur −ρ g dω qui n’est autre que le poids changé de signe de
r
l’élément de fluide. En dynamique, le vecteur grad p(M) est couplé à : (i) l’accélération ρ dV / dt , (ii) la
variation grad [ ρ U(M) ] de l’énergie potentielle du fluide en mouvement, (iii) aux forces de surfaces.

Remarque 2: Equivalence entre l’intégrale de surface et l’intégrale de volume. Sachant que

r p( M )
F ( M ) = grad (I-12)
ρ( M )

d’après (I-9) (fluide incompressible), on a en reportant dans (I-4),

p( M )

M ∈D
ρ( M ) grad
ρ( M )
dω + ∫ − p(M *) n(M *)dA =0
M *∈S
(I-13)

soit la relation remarquable :

∫ p(M *) n(M *) dA = ∫ grad p(M )dω


M *∈S M ∈D
(I-14)

Mécanique_des_fluides_2A_MP.doc 5 21/11/2007
II- Fluide homogène incompressible dans un repère terrestre (galiléen ap-
proché).
II-1 Equation locale
r
La densité massique de force à distance est alors l’accélération locale de la pesanteur g et dérive du potentiel U
tel que :

r  ∂U r ∂U r ∂U r
g = −grad U = −  x+ y+ z (II-1)
 ∂x ∂y ∂z 
r
si z est la verticale ascendante du lieu, les composantes du vecteur g dans R sont alors :

∂U ∂U ∂U
0= , 0= , g= (II-2)
∂x ∂y ∂z

On en déduit que U ne dépend que de z et a pour expression :

U ( M ) = gz ( M ) + U 0 (II-3)

où U 0 est une constante d’intégration.

Puisque la masse volumique est constante en tout point de D, ρ(M ) = ρ , l’équation fondamentale SDF-3
s’écrit :

p( M ) + ρgz ( M ) = Cte, ∀M ∈ D (II-4)

On appelle pression motrice le terme p( M ) + ρgz ( M ) et on le note p g (M ) . L’équation SDF-5 devient donc :

p g ( M ) = Cte, ∀M ∈ D (II-5)

Dans un fluide au repos, la pression motrice est constante.

II-2 Equation des isobares et de la surface libre

Sur une isobare, p( M ) = C M = Cte , on peut donc écrire : ρgz ( M ) = Cte , soit : z ( M ) = H M = cte , c’est
l’équation d’un plan horizontal de cote H M .

La surface libre est une isobare particulière, sur laquelle p(M *) est une constante par exemple la pression at-
mosphérique, et d’équation : z ( M *) = H 0 .

Mécanique_des_fluides_2A_MP.doc 6 21/11/2007
II-3 Equation fondamentale intégrale – Théorème d’ARCHIMEDE

n(M*)
M*
K D

Soit un corps K, de frontière (surface) S et de volume

V= ∫ dω(M )
M ∈K
(II-6)

immergé dans un liquide.

D’après le principe fondamental appliqué au domaine D de fluide qui remplace fictivement l’espace occupé par
ce corps, on a :

r r
ρg ∫
M ∈D
dω + ∫ − p(M *) n(M *)dA = 0
M *∈S
(II-7)

D’autre part les forces de pression exercée par le fluide extérieur à D, (c’est-à-dire D ) sur la surface S
s’écrivent :
r
R( D / K ) = ∫ − p(M *) n(M *) dA
M *∈S
(II-8)

On en déduit que :
r r
R ( D / K ) = −ρ g V (II-9)

C’est la formule donnée par le théorème d’Archimède : tout corps K immergé dans un liquide est soumis de la
r
part de celui-ci à une poussée R ( D / K ) verticale dirigée vers « le haut », de module égal au poids du volume V
de liquide déplacé par le corps.

Interprétation.

En statique, la présence d’un corps immergé ne modifie pas le champ de pression dans le liquide. Le champ de
pression p(z) est identique, qu’il y ait ou qu’il n’y ait pas de corps. L’intégrale II-8 est donc identique en
l’absence du corps, si on la calcule sur une surface (S) identique à celle que le corps occupera quand on va
r
l’immerger. En l’absence du corps, la résultante R(D / K) est égale à l’opposé du poids du volume de liquide
enfermé dans (S) et en équilibre. D’où II-9 et la phrase consacrée qui suit.

III-Fluide homogène compressible dans un repère terrestre (galiléen appro-


ché).
III-1 Equation locale
On a toujours l’équation différentielle au point M :

Mécanique_des_fluides_2A_MP.doc 7 21/11/2007
dp ( M )
= −ρ( M ) g (III-1)
dz

A/ Cas d’un LIQUIDE compressible

Soient : ρ(M ) la masse volumique du liquide à la pression en un point M,


ρ 0 = ρ( M 0 ) et p( M 0 )
la masse volumique et la pression en un point M 0 de référence dans le fluide supposé à température constante.
L’équation d’état du liquide s’écrit alors en fonction du coefficient de compressibilité isotherme du fluide, χ :

ρ( M ) = ρ( M 0 )[1 + χ( p ( M ) − p( M 0 ) )]

L’équation différentielle s’écrit alors :

dp ( M )
ρ( M 0 )[1 + χ( p ( M ) − p( M 0 ) )] g = − (III-2)
dz

elle s’intègre par séparation des variables :

[1 + χ( p( M ) − p(M 0 ) )] g = e χρ(M 0 ) g [z (M 0 ) − z (M )] (III-3)

où z ( M 0 ) et z (M ) sont respectivement les cotes de M 0 et de M.

En général, l’exposant est très petit devant 1, on peut alors remplacer l’exponentielle par son approximation au
1er ordre, d’où :

 χ ρ( M 0 ) g z ( M 0 ) − z ( M ) 
p( M ) − p ( M 0 ) = ρ( M 0 ) g [z ( M 0 ) − z ( M )] 1 + 
 2 

Si on choisit le point de référence sur la surface libre de cote z ( M 0 ) = 0 , et que M est à la profondeur
h = − z (M ) , on a :

 χ ρ( M 0 ) g h 
p( M ) − p ( M 0 ) = ρ( M 0 ) g h 1 +  (III-4)
 2 

Le terme de compressibilité ne devient significatif qu’aux profondeurs importantes. Pour l’eau, χ = 5 E − 10 Pa −1


dans les conditions normales, le terme χρgh vaut 0,01 à partir de 2000 mètres.

B/ Cas d’un gaz.

Soient encore : ρ( M 0 ) , p( M 0 ) , T ( M 0 ) , la masse volumique, la pression, la température, au point de référence


M0.
En considérant que le gaz est parfait, on a :

p(M ) T (M 0 )
ρ( M ) = ρ ( M 0 ) (III-5)
p(M 0 ) T (M )

Dans un gaz à température constante, on a :

Mécanique_des_fluides_2A_MP.doc 8 21/11/2007
p(M ) dp ( M )
ρ( M 0 ) g=−
p(M 0 ) dz

L’intégration conduit à :

ρ(M 0 ) g
p(M )
[z ( M 0 ) − z ( M )]
p( M 0 )
g =e (III-6)
p(M 0 )

Cette formule permet d’estimer l’évolution (décroissante) de la pression avec l’altitude.

Mécanique_des_fluides_2A_MP.doc 9 21/11/2007
DYNAMIQUE DES FLUIDES

I- Coordonnées de Lagrange et coordonnées d’Euler


Les coordonnées de LAGRANGE sont un ensemble de trois paramètres a1 , a 2 , a 3 , qui caractérisent la posi-
tion de la particule à l’instant initial t = 0. Il s’agit des coordonnées dans le repère fixe de référence du point
avec lequel coïncide la particule à l’instant t = 0.

Les coordonnées d’EULER sont un ensemble de trois paramètres x, y, z , caractérisant la position de la particule
à l’instant t considéré. Il s’agit des coordonnées dans le même repère fixe que précédemment, du point avec le-
quel coïncide la particule à l’instant t.

On notera bien que la direction z du repère fixe n’est pas nécessairement confondue avec la verticale du
lieu.

Pour une particule donnée, donc que l’on « suit dans son mouvement » , on a pour t > 0 :

x = f 1 ( a1 , a 2 , a 3 , t ) a1 = g 1 ( x , y , z , t )
y = f 2 ( a1 , a 2 , a 3 , t ) a 2 = g 2 ( x, y , z , t ) (I-1)
z = f 3 ( a1 , a 2 , a 3 , t ) a 3 = g 3 ( x, y , z , t )

A l’instant initial, on a :

f 1 (a1 , a 2 , a 3 ,0) = a1
f 2 (a1 , a 2 , a 3 ,0) = a 2 (I-2)
f 3 (a1 , a 2 , a 3 ,0) = a 3

II- Vitesse
r r r
La vitesse d’une particule est décrite par ses trois composantes u, v, w, sur la base orthonormée ( x , y, z ) du re-
père de référence, c’est-à-dire :
r r r r
V = u ( x, y, z ) x + v( x, y, z ) y + w( x, y , z ) z (II-1)

∂x
u ( x, y , z , t ) =
∂t
(II-2)

= f 1 ( a1 , a 2 , a 3 , t )
∂t
∂y
v ( x, y , z , t ) =
∂t
(II-3)

= f 2 ( a1 , a 2 , a 3 , t )
∂t
∂z
w( x, y, z , t ) =
∂t
(II-4)

= f 3 ( a1 , a 2 , a 3 , t )
∂t

Mécanique_des_fluides_2A_MP.doc 10 21/11/2007
III- Dérivée « particulaire »
Soit une grandeur physique scalaire φ liée à la particule de fluide en mouvement (concentration, masse volumi-
que, composante de vitesse,….).

Si cette grandeur est exprimée en fonction des coordonnées d’EULER de la particule (x,y,z), on peut écrire :

φ ( M , t ) = φ ( x, y , z , t ) (III-1)

On appelle dérivée particulaire de φ , sa dérivée totale par rapport au temps pour a1 , a 2 , a 3 , donnés, soit :

dφ ∂φ ∂φ dx ∂φ dy ∂φ dz
= + + +
dt ∂t ∂x dt ∂y dt ∂z dt
(III-2)
∂φ r
= + (V ⋅ grad) φ
∂t
r
V est le vecteur vitesse de la particule située en M à l’instant t.
r
(V ⋅ grad) désigne l’opérateur de convection appliqué ici à la grandeur φ .
r r r
Dans le cas d’une grandeur vectorielle quelconque, A( M , t ) = A( x, y, z , t ) , le vecteur dA / dt a pour composantes
les dérivées particulaires de chaque composante :

 dAx ∂Ax r
 dt = ∂t + (V ⋅ grad) Ax
r 
dA  dAy ∂Ay r
= = + (V ⋅ grad) Ay (III-3)
dt  dt ∂t
 dA ∂A r
 z = z + (V ⋅ grad) Az
 dt ∂t

IV- Accélération
r
Lorsque la grandeur vectorielle A est la vitesse de la particule au point M à l’instant t, on obtient le vecteur ac-
célération et ses trois composantes cartésiennes :

 du ∂u r
 dt = ∂t + (V ⋅ grad) u
r r 
r dV ∂V r r  dv ∂v r
Γ= = + (V ⋅ grad) V =  = + (V ⋅ grad) v (IV-4)
dt ∂t  dt ∂t
 dw ∂w r
 = + (V ⋅ grad) w
 dt ∂t

r
∂V
Le terme représente l’accélération locale ou temporelle.
∂t
r r
Le terme (V ⋅ grad) V représente l’accélération convective.

Le vecteur accélération peut se mettre sous la forme qui fait apparaître la vorticité

r 1 r
ω = rot V (IV-5)
2

Mécanique_des_fluides_2A_MP.doc 11 21/11/2007
r r
dV ∂V V2 r r
= + grad + (rot V ) ∧ V
dt ∂t 2
r (IV-6)
∂V V2 r r
= + grad + 2ω ∧ V
∂t 2

V- Champ des vitesses dans un écoulement permanent (ou stationnaire).


On dit que l’écoulement est permanent ou stationnaire si ses propriétés ne dépendent pas explicitement du temps.

En particulier, la vitesse ne dépend pas explicitement du temps, c’est-à-dire que les composantes sont des fonc-
tions uniquement des coordonnées d’EULER x, y, z.
On a donc :

∂u ∂v ∂w
= = =0 (V-1)
∂t ∂t ∂t

Les composantes du vecteur accélération se réduisent à :


r r r
Γx = (V ⋅ grad) u Γy = (V ⋅ grad) v Γz = (V ⋅ grad) w (V-2)

L’expression du vecteur accélération ne contient alors que l’accélération convective :


r
dV V2 r r
= grad + (rot V ) ∧ V (V-3)
dt 2

Si l’écoulement est stationnaire les trajectoires et les lignes de courant sont confondues.

VI. Ecoulement irrotationnel.


On dit qu’à un instant t donné, l’écoulement est irrotationnel dans le domaine D si le rotationnel de la vitesse
r
V ( M , t ) est nul en chacun des points M de D. Le rotationnel est un vecteur (défini en cinématique des fluides)
qui caractérise en un point M, la rotation dans l’espace de la particule fluide.

Il s’écrit symboliquement :
r r r
rot V ( M , t ) = ∇ ∧ V ( M , t ) (VI-1)

Il s’écrit par exemple en coordonnées cartésiennes :

r  ∂w( M , t ) ∂v( M , t )  r  ∂u ( M , t ) ∂w( M , t )  r  ∂v( M , t ) ∂u ( M , t )  r


rot V ( M , t ) =  − x+ − y+ − z (VI-2)
 ∂y ∂z   ∂z ∂x   ∂x ∂y 

VII. Equation locale de conservation de la masse.


Soit un domaine D de fluide que l’ « on suit dans son mouvement ».
La masse m de liquide contenue dans D est constante dans le temps et dans l’espace (c-à-d au cours de
l’écoulement).

Mécanique_des_fluides_2A_MP.doc 12 21/11/2007
dm
La dérivée particulaire de m est donc nulle : =0
dt
On écrit m sous forme de l’intégale de volume :

m= ∫ ρ(M , t ) dω
M ∈D

dm d
=
dt dt ∫ ρ( M , t )dω = ∫
 ∂ρ( M , t )
 ∂t [ r
] 
+ div ρ( M , t ) V ( M , t )  dω = 0 (VII-1)
M ∈D M ∈D
 

Cette relation1 s’applique à tout élément de volume, et donc en particulier au volume élémentaire dω ; ainsi, en
supprimant le signe d’intégration, on obtient l’équation dite équation de continuité ou de conservation de la
masse :

∂ρ( M , t )
∂t
[ r
]
+ div ρ( M , t ) V ( M , t ) = 0 (VII-2)

Comme

r
[ ] r r
(
div ρV = ρ divV + V ⋅ grad ρ )
on peut faire apparaître la dérivée particulaire de la masse volumique, fonction scalaire du point M et du temps t :

dρ ( M , t ) r
+ ρ( M , t ) ε = 0 ε = divV ( M , t ) (VII-3)
dt

r
Lorsque la dérivée particulaire de la masse volumique dρ / dt est nulle en tout point, on a divV ( M , t ) = 0 et
l’écoulement est isovolume.

On montre en effet que ε représente le taux de déformation cubique (ou vitesse de dilatation) de la particule si-
tuée au point M.
r
Il est important de se souvenir dans quel cas dρ / dt = ∂ρ / ∂t + (V ⋅ grad)ρ = 0 .
r
Statique des Fluides : ∂ρ / ∂t = 0 et V = 0 ⇒ dρ / dt = 0 ;
Fluides homogènes, incompressibles et indilatables : ∂ρ / ∂t = 0 et grad ρ = 0 ⇒ dρ / dt = 0 ;
Lignes de courant confondues avec les lignes isochores (où ρ = Cte ) et masse volumique constante dans le
r
temps (fluide inhomogène) : ∂ρ / ∂t = 0 et V est perpendiculaire à grad ρ ⇒ dρ / dt = 0

Lorsque l’écoulement est permanent ou stationnaire (le mouvement ne dépend pas explicitement du temps), la
dérivée partielle de la masse volumique par rapport au temps est nulle :

∂ρ( M , t )
∂t
[ r
] [ r
+ div ρ( M , t ) V ( M , t ) = div ρ( M , t ) V ( M , t ) = 0 ] (VII-4)

1
r r r r r ∂u ∂v ∂w
On appelle « divergence » d’un vecteur V = ux + vy + wz la quantité (scalaire) : div V = + + .
∂x ∂y ∂z

Mécanique_des_fluides_2A_MP.doc 13 21/11/2007
r
On dit alors que le champ de vecteurs ρV est conservatif.

VIII. Equation intégrale de conservation de la masse. Débit de liquide à


travers une section de passage.

VIII-1 Equation intégrale


Reprenons l’équation de conservation :


 ∂ρ( M , t )
 ∂t [ r 
+ div ρ( M , t ) V ( M , t )  dω = ]
∂ρ( M , t )
∂t ∫ ∫ [
r
]
dω + div ρ( M , t ) V ( M , t ) dω (VIII-1)
M ∈D
  M ∈D M ∈D

et transformons la deuxième intégrale de volume en une intégrale de surface en exprimant que la divergence du
r
vecteur ρ( M , t ) V ( M , t ) dans le domaine D est égal au flux de ce vecteur à travers S:

∫ [ ]
r r r
M ∈D
div ρ( M , t ) V ( M , t ) dω = ∫
M *∈S
ρ( M *, t ) V ( M *, t ) n ( M *)dA (VIII-2)

Or la projection
r r r r
V ( M *, t ) n ( M *)dA = V ( M *, t )dA = dqV (VIII-3)
r
est le débit élémentaire volumique dqV à travers l’élément de surface dA orienté par n (M *) . Ainsi on a :

∂ρ( M , t )

M ∈D
∂t
dω + ∫ ρ(M *, t ) dq
M *∈S
V ( M *) = 0

Lorsque le mouvement du fluide (écoulement) ne dépend pas du temps, on a ∂ρ / ∂t = 0 et donc :

∂ρ( M , t )

M ∈D
∂t
dω = 0 ⇒ ∫ ρ(M *, t )dq
M *∈S
V ( M *) = 0 (VIII-4)

Dans les applications où le fluide est incompressible, ρ( M *, t ) = ρ = Cte peut être « sortie » de l’intégrale :

∫ ρ(M *, t )dq
M *∈S
V ( M *) = ρ ∫ dq
M *∈S
V ( M *) = 0 ⇒ qV = 0 (VIII-5)

L’intégrale de surface ∫ dq
M *∈S
V ( M *) est tout simplement le débit en volume total à travers l’enveloppe S du do-

maine fluide D. On voit que ce flux total est nul.

Si la masse volumique n’est pas constante, il faut raisonner sur le débit en masse, dq m = ρdq V , d’où :

∫ ρ(M *, t )dq
M *∈S
V ( M *) = ∫ dq
M *∈S
m ( M *) = 0 ⇒ qm = 0 (VIII-6)

Mécanique_des_fluides_2A_MP.doc 14 21/11/2007
VIII-2 Débit à travers la section droite d’un tuyau.
Ici le domaine D de fluide est un tronçon de tube constitué par une surface amont S1, une surface aval S2, et une
surface latérale SL. Pour le débit en volume, en tout point de SL la projection
r r r r
V ( M *, t ) n ( M *)dA = V ( M *, t )dA = dqV est nulle, puisque le liquide ne peut traverser la paroi (la vitesse est per-
pendiculaire à la normale). Il en est de même pour le débit en masse.

Le débit en volume est utilisé en pratique pour les fluides incompressibles.


Le débit en masse est utilisé en pratique pour les fluides compressibles (gaz).

S1 n*
SL surface latérale
n1
q
V1
S2
q
V2
n2
on a donc :

qV = qV 1 + q V 2 = 0 ⇒ q V 1 = − qV 2 ou : qV 1 = qV 2 = CTE Le débit en volume est constant à travers


les sections droites

q m = q m1 + q m 2 = 0 ⇒ q m1 = − q m 2 ou : q m1 = q m 2 = CTE Le débit en masse est constant à travers


les sections droites
r
Remarque : le signe du débit dépend de l’orientation choisie pour la normale à la section droite. Si n est dans le
sens de l’écoulement, le débit est positif, négatif dans le cas contraire.

Si on raisonne (cas fréquent) sur les valeurs absolues, on utilise la valeur absolue qV pour le débit en volume,
q m pour le débit en masse.

IX. Principe fondamental de la dynamique appliqué à un fluide parfait–


Equation d’EULER
On cherche la relation entre le mouvement des particules fluides et les forces qui s’exercent sur ces particules.

On applique ci-dessous le PFD à un domaine D de fluide que l’on suit dans son mouvement.

M est un point de D.

S désignera l’enveloppe (ou surface) de D.

M* est un point de S
r r r
Le mouvement est étudié dans un repère absolu R (O, x , y, z )

Mécanique_des_fluides_2A_MP.doc 15 21/11/2007
IX-1 Théorème de la résultante en fluide parfait
r
A chaque instant, la quantité d’accélération A/ R ( D, t ) du domaine D est égale à la somme des forces extérieures
qui s’appliquent sur D (on reprend les notations vues en statique des fluides):
r v r
A/ R ( D, t ) = R( F / D) + R ( D / D) (IX-1)
v r
Ces forces comportent les forces de gravité R ( F / D) de densité massique g ,

r r
R( F / D) = ∫ ρ ( M , t ) g dω
M ∈D
(IX-2)
r
et la résultante R ( D / D) des forces de surface exercées par le complémentaire D de D (fluide extérieur) sur
l’enveloppe S ,
r r
R( D / D) = ∫ − p(M *, t ) n(M *) d A
M *∈S
(IX-3)
r
où n (M *) est la normale extérieure (unitaire) au point M* de S.
r
R ( D / D) représente la résultante des forces de pression qui s’exercent sur S à l’instant considéré.

r
dV ( M , t )
La quantité d’accélération d’une particule étant égale à ρ( M , t ) dω , la quantité d’accélération de D (ou
dt
résultante des quantités d’accélération du domaine D) est donnée par l’intégrale :
r
r dV ( M , t )
A/ R ρ( M , t ) = ∫
M ∈D
ρ( M , t )
dt
dω (IX-4)

Finalement, le PFD stipule que :


r
dV ( M , t ) r r

M ∈D
ρ( M , t )
dt
dω =
M ∈D
∫ ρ(M , t ) g dω + ∫ − p(M *, t ) n(M *)dA
M *∈S
(IX-4)

L’intégrale des pressions sur la surface peut s’écrire aussi sous la forme d’une intégrale de volume (ceci se tra-
duit physiquement par la poussée d’Archimède en statique des fluides) :

r
∫ − p(M *, t ) n(M *)dA = ∫ − grad [ p(M *, t ) ] dA
M *∈S M ∈D
(IX-5)

Il s’ensuit :
r
 dV ( M , t ) r r 
∫ 

M ∈D 
ρ ( M , t )
dt
− ρ( M , t ) g + grad p ( M *, t ) n ( M *)dω = 0


Localement, pour une particule fluide de centre M, l’équation du mouvement s’écrit donc :
r
dV ( M , t ) r
ρ( M , t ) = ρ( M , t ) g − grad p( M , t ) (IX-6)
dt

Cette relation est appelée équation d’EULER.

L’accélération de la pesanteur s’écrit en introduisant le potentiel scalaire U(M):

Mécanique_des_fluides_2A_MP.doc 16 21/11/2007
r
g = −grad[U ( M )]
(IX-7)
U ( M ) = g Z( M )

où Z(M) est la cote du point M sur l’axe vertical ascendant du lieu (et non la coordonnée z).

Si la masse volumique est constante, on a :


r
dV ( M , t )  p( M , t ) 
= −grad[ U ( M )] − grad  
dt  ρ( M , t ) 
 p(M , t ) 
= −grad  U ( M ) +
 ρ( M , t ) 

Utilisons l’expression de l’accélération :


r r
dV ∂V V2 r r
= + grad + (rot V ) ∧ V (IX-8)
dt ∂t 2

∂V ( M , t ) V 2 (M , t ) r r  p(M , t ) 
+ grad + (rot V ( M , t ) ∧ V ( M , t ) = −grad  U ( M ) +
∂t 2  ρ ( M , t ) 

∂V ( M , t ) r r V 2 ( M , t ) p( M , t ) 
+ (rot V ( M , t ) ∧ V ( M , t ) = −grad  + U (M ) +
∂t  2 ρ( M , t ) 

Posons :

V 2 (M , t ) p(M , t )
H (M , t ) = + U (M ) + (IX-9)
2 ρ( M , t )

le premier membre de l’équation précédente vérifie donc :

∂V ( M , t ) r r
+ (rot V ( M , t ) ∧ V ( M , t ) = −grad H ( M , t ) (IX-10)
∂t

IX-2 Intégrale première de l’équation du mouvement stationnaire


Equation de Daniel BERNOULLI (1700-1782)..
En mouvement stationnaire, on a en tout point M :

∂V ( M , t )
=0 (IX-11)
∂t

De plus comme
r
dH = grad H ( M , t ) ⋅ d l (IX-12)

la fonction H doit vérifier :


r r r
− (rot V ( M , t ) ∧ V ( M , t ) ⋅ d l = dH (IX-13)

Mécanique_des_fluides_2A_MP.doc 17 21/11/2007
r r
Ainsi, selon que la direction du vecteur déplacement élémentaire d l est confondue avec la vitesse V ( M , t ) ou
r
le rotationnel rot V ( M , t ) , on obtient respectivement, pour les intégrales H, les propriétés suivantes :

r r
d l colinéaire à V(M,t) ⇒ H ( M , t ) = Constante sur une ligne de courant ;
r r
d l colinéaire à rot V(M,t) ⇒ H ( M , t ) = Constante sur une ligne tourbillon

Une ligne tourbillon est tangente en chacun de ses points au vecteur tourbillon.
Dans un écoulement rotationnel, H est diffère d’une ligne tourbillon à une autre.

IX-3 Equation de BERNOULLI


On appelle équation de BERNOULLI l’intégrale constante de l’équation du mouvement sur une ligne de cou-
rant. Elle s’applique au cas d’un fluide parfait, en écoulement permanent et irrotationnel.

On a donc pour une ligne de courant donnée Li:

V 2 (M , t ) p( M , t )
∀M ∈ Li, H (M , t ) = + U (M ) + = Constante (IX-14)
2 ρ( M , t )

Si, dans l’écoulement, le rotationnel est nul en tout point , alors la fonction H est la même pour toutes les lignes
de courant, c’est-à-dire dans tout le fluide :

X- Aspect énergétique de l’équation de BERNOULLI


Raisonnons sur le travail des forces subies par la particule de volume dω lorsqu’elle effectue un déplacement
donné par le vecteur déplacement
r r r r
d l = dx x + dy y + dz z (X-1)

Ces forces sont :


r
[
-Les forces de pression au point M, d’intensité dF p ( M ) = − grad p ( M ) dω , (X-2)]
r r r
-Les forces de pesanteur : dFg ( M ) = −ρdω g z = ρdω g = −ρ grad U ( M ) dω . (X-3)

r
Lorsque la particule effectue un déplacement d l , le travail des forces de pression et de pesanteur est égal à la
variation d’énergie cinétique :

δ W p + δW g = δ E c (X-4)

[
δE c = d ρdω V 2 ( M ) / 2 ]
r r
[ ] r
δW p = dF p ( M ) ⋅ d l = − grad p ( M ) dω ⋅ d l
r r
[ ] r
δW g = dFg ( M ) ⋅ d l = ρ − grad U ( M ) dω ⋅ d l (X-5)

On se rappellera ici que la différentielle d’une fonction scalaire f(M) peut s’écrire :
r
df ( M ) = grad f ( M ) ⋅ d l (X-6)
Le travail de la pression et de la pesanteur s’écrivent donc :

Mécanique_des_fluides_2A_MP.doc 18 21/11/2007
δW p = dp ( M ) dω δW g = ρ dU ( M ) dω = d [ρU ( M )] dω

Soit :

[ ]
d ρdω V 2 ( M ) / 2 = dp ( M )dω + d [ρU ( M )]dω (X-7)
[ ]
0 = d p ( M ) + ρU ( M ) + ρ V 2 ( M ) / 2 dω

Comme le volume dω de la particule ne peut être nul, on en déduit que c’est la différentielle qui est nulle :

[ ]
d p( M ) + ρU ( M ) + ρ V 2 ( M ) / 2 = dPt = 0 Pt ( M ) = p ( M ) + ρU ( M ) +
1
2
ρ V 2 (M )

La différentielle d’une fonction constante étant nulle, on en déduit que la fonction Pt (M ) , que l’on appelle
pression totale, est constante :

Pt ( M ) = p ( M ) + ρU ( M ) + ρ V 2 ( M ) / 2 = Constante (X-8)

C’est l’expression mathématique du théorème de Daniel BERNOULLI : dans un écoulement irrotationnel de


fluide parfait incompressible, dans le champ de pesanteur, la pression totale est constante sur une ligne de cou-
rant.

Dans cette démonstration, on s’est attaché à une particule fluide donnée. Le résultat concerne donc implicitement
la même ligne de courant.

Si l’écoulement est de plus irrotationnel en tout point, la pression totale est identique en tout point.

XI- L’équation de BERNOULLI, en pratique


D’abord :

On remarquera que les trois termes de l’équation de Bernoulli doivent se compenser en tout
point M de la ligne de courant considérée, de manière à ce que la pression totale soit CONS-
TANTE.

Ensuite :
La pression totale a la dimension d’une pression ( ML−1 T −2 ), son unité SI est le Pa.
La pression totale étant homogène à une pression, elle peut être exprimée par une hauteur de fluide équivalente
qui a la dimension d’une longueur (L).

Selon l’unité de travail choisie, les différents termes de l’équation de Bernoulli sont dénommés comme suit :

Formulation homogène à une pression ( ML−1 T −2 ) Formulation homogène à une longueur (L) ou
hauteur équivalente
p(M ) V 2 (M )
Pt ( M ) = p ( M ) + ρg Z( M ) + ρV 2 ( M ) / 2 H (M ) = + Z( M ) +
ρg 2g

Mécanique_des_fluides_2A_MP.doc 19 21/11/2007
H ( M ) = Charge, ou hauteur totale
Pt ( M ) = Pression totale
p(M )
p( M ) + ρgZ( M ) = Pression motrice + Z( M ) = Hauteur piézométrique
ρg
ρV 2 ( M ) / 2 = Pression dynamique
V 2 (M )
= Hauteur cinétique
2g

XII- Pression motrice dans un écoulement unidirectionnel stationnaire


On rappelle que la pression motrice est définie en un point M à l’instant t par :

p g ( M ) = p ( M ) + ρg Z ( M ) (XII-1)

et que Z(M) est la cote du point M sur la verticale ascendante du lieu.

Il n’y pas forcément une relation simple entre z(M) coordonnées de M sur l’axe z du repère d’étude et Z(M),
contrairement à la statique des fluides. Dans la figure ci-dessous on a par exemple représenté un écoulement des-
cendant. ∂Z / ∂x représentera la pente de l’écoulement.

Axe vertical Z
g

Plan y-z où la pression motrice est constante

,z)
M(x,y
y

V= u
x

Considérons le cas simple mais instructif d’un écoulement stationnaire de fluide incompressible de la forme :
r r
V (M ) = u(M ) x (XII-2)

r r
Les vitesses sont toujours colinéaires à x . L’équation locale de continuité se réduit ici à divV = 0 , d’où :

r ∂u
div V = = 0 ⇒ u est indépendante de x (XII-3)
∂x

Ainsi dans un tel mouvement, la vitesse ne peut être que constante ; c’est un écoulement uniforme.

Les composantes Γ y et Γz de l’accélération sont nulles ; seules Γx peut être non nulle :

Mécanique_des_fluides_2A_MP.doc 20 21/11/2007
dv 1 ∂
Γy = =− [ p + ρg Z ] = 0 ⇒ p + ρg Z est indépendante de y
dt ρ ∂y
(XII-4)
dw 1 ∂
Γz = =− [ p + ρg Z ] = 0 ⇒ p + ρg Z est indépendante de z
dt ρ ∂z

Ainsi la pression motrice p g = p + ρgZ est constante perpendiculairement au plan de l’écoulement uniforme. On
dit alors que la pression p a une répartition hydrostatique car cette répartition est identique à celle donnée par la
statique des fluides.

Mécanique_des_fluides_2A_MP.doc 21 21/11/2007
XIII. Ecoulement d’un fluide visqueux newtonien incompressible.
XIII-1 Ecoulement interne et écoulement externe
Les écoulements de fluide visqueux peuvent être divisés en deux catégories :

1/ Les écoulements internes qui sont limités par des parois (écoulements dans les tuyaux, canaux, tuyères…) ;

2/ Les écoulements externes, qui sont des mouvements de fluides autour de corps tels que des avions, fusées,
voitures, etc…

XIII-2 Ecoulement laminaire et écoulements turbulent


Un écoulement turbulent peut être caractérisé par une orientation aléatoire (ou fluctuation) des vecteurs vitesses
r
en chaque point ; autrement dit chaque composate u(M,t), v(M,t), w(M,t), de V ( M , t ) obéissent à des lois de dis-
tributions (au sens probabiliste) en fonction de l’espace et du temps. On dit que le régime de l’écoulement est
r
turbulent. Un écoulement turbulent peut posséder une composante moyenne U (le mouvement global se fait
d’ouest en en est par exemple) pouvant dépendre du temps mais de manière plus « régulière »; les fluctuations
ont alors lieu autour de cette composante moyenne, de sorte que la vitesse s’écrit :
r r
V ( M , t ) = U + v' (XIII-1)
C’est alors sur v' que l’on fait apparaître des lois statistiques.

Dans un écoulement laminaire, le mouvement du fluide s’effectue en « couches » parallèle entre elles, chaque
couche possédant sa propre vitesse ; le profil de vitesses dans le fluide en mouvement est bien ordonné. Les
vecteurs vitesse conservent une orientation stable au cours du temps. On dit que le régime de l’écoulement est
laminaire.

Lorsque le mouvement laminaire du fluide dégénère en un écoulement turbulent , il perd son caractère ordonné
et stable ; on dit qu’il y a transition du régime laminaire vers le régime turbulent, ou plus simplement transition
laminaire turbulent.

La différence entre le régime turbulent et le régime laminaire peut être représentée par les courbes donnant la
vitesse en fonction du temps en un point M fixe.

Evolution des vitesses avec le temps en un point fixe

Ecoulement turbulent instationnaire


Vitesse

U fonction du temps

Ecoulement turbulent stationnaire

Ecoulement laminaire stationnaire


U constante

temps

XIII-3 Ecoulement laminaire établi entre deux plaques parallèles infinies


Le fait de considérer un écoulement stationnaire permet de se ramener à un problème indépendant du temps.

Mécanique_des_fluides_2A_MP.doc 22 21/11/2007
Le fait de considérer un écoulement établi permet de se ramener à un problème indépendant de l’abscisse x : le
mouvement n’est ni accéléré ni ralenti selon la position sur l’axe x ; u ne dépend pas de x.
Le fait de considérer des plaques infinies permet de se ramener à un problème plan r (à deux dimensions) ; autre-
r
ment dit, la répartition des vitesses est invariante par translation suivant z, soit : V ( M ) = u ( y ) x
Comme le fluide est visqueux, les parois des plaques fixes imposent une condition limite à la vitesse : la vitesse
du fluide qui touche une paroi doit être égale à celle de la paroi : c’est la condition d’adhérence.

Lorsque le fluide est visqueux, il faut rajouter dans le 2e membre de l’équation du mouvement (c’est l’éq.
d’Euler si le fluide est parfait) les forces de frottement d’origine visqueuse. Si le fluide est newtonien, sa visco-
sité est constante pour une température donnée. Les efforts visqueux sont alors proportionnels aux gradients de
vitesses. La formulation de ces effortsrvisqueux est simple dans le cas de l’écoulement entre deux plaques car le
r
champ de vitesse est unidirectionnel, V ( M ) = u ( M ) x .

Volume de contrôle
y

Ecoulement dy a
dx
x
Cas où les deux plaques sont fixes.
On a comme conditions limites :

y=0 ⇒ u=0
(XIII-2)
y=a ⇒ u=0

r
on cherche une vitesse V dont les composantes sont telles que :

u = u( y) v = w = 0 (XIII-3)
Après calculs, on trouve pour la composante u :

a 2 dp g  y 2 y 
u( y) = − (XIII-4)
2µ dx  a 2 a 
dp g / dx est appelé gradient de pression motrice de l’écoulement ; c’est une chute de pression motrice propor-
tionnelle à la longueur parcourue : la chute de pression divisée par la longueur parcourue est constante car
l’écoulement est établi (la pression motrice décroît linéairement, voir la figure ci-dessous). Ainsi, entre deux
sections distantes d’une longueur L on a :

dp g pg2 − pg2 ∆p g
= =− (XIII-5)
dx x 2 − x1 L
Pour déterminer la dérivée de la pression motrice par rapport à x il n’est donc pas obligatoire de se ramener au
voisinage d’un point. La distinction des deux cas où dp g / dx = Cte et dp g / dx = f ( x) sont illustrés ci-dessous.

Mécanique_des_fluides_2A_MP.doc 23 21/11/2007
Pression motrice pg
L
∂pg ∆pg
= Cte =
∂x L

∆pg
∂pg
= f(x)
∂x

Distance parcourue

Entre les deux plaques, on a une distribution de contraintes visqueuses donnée par :

1 dp g 3
σ yx = − a (XIII-6)
12µ dx
Le débit volumique est l’intégrale de u sur une section droite ;

a a 3 ∂p g

qV = u ( y )dy dz = −
0 ∫
z
12µ ∂x ∫ dz
z
(XIII-7)

On choisit d’exprimer ici le débit par unité de largeur (la largeur est perpendiculaire à la figure, soit z) ; ce débit
unitaire est donné par :

qV
qV ' =
∫ dzz (XIII-8)
a a ∂p g
3
a ∆p g 3


qV ' = u ( y )dy = −
0
=
12µ ∂x 12µ L

La vitesse moyenne de l’écoulement est obtenue à partir du débit et de l’aire A d’une section de passage de réfé-
rence S perpendiculaire à l’écoulement:

∫ ∫
A = dy dz
0 z
(XIII-9)
q a 2 ∆p g
Vq = V =
A 12µ L

Le point où la vitesse est maximale est situé à :

a 3
y= ⇒ u max = V (XIII-10)
2 2 q

La vitesse maximale est ici 1,5 fois plus grande que la vitesse moyenne.

XIII-4 Ecoulement laminaire développé établi dans une conduite de section circulaire

La « longueur d’entrée », L E , est la longueur de tuyau entre l’entrée du tube et la position où l’écoulement lami-
naire est complètement développé et établi. Cette longueur est donnée par une formule empirique :

L E = 0,058 Re D (XIII-11)

Mécanique_des_fluides_2A_MP.doc 24 21/11/2007
où D = 2 R est le diamètre intérieur du tube et re, le nombre de REYNOLDS de l’écoulement, défini à partir de
la vitesse moyenne V q du diamètre D, de la masse volumique ρ et de la viscosité dynamique µ du fluide :

ρV q D
Re = (XIII-12)
µ

La vitesse de l’écoulement établi est ici aussi de la forme :


r r
V = ux (XIII-13)

où l’on peut montrer que u ne dépend que du rayon r dans une section droite ( 0 ≤ r ≤ R ).
La vitesse vérifie la condition d’adhérence sur la paroi :

r=R ⇒ u=0 (XIII-14)

Le profil des vitesses dans la section droite est :

1 ∂p g 2
u (r ) =
4µ ∂x
(
r − R2 ) (XIII-15)

Le gradient de pression motrice est ici aussi une constante pour un écoulement donné ; on a donc, entre deux
sections distantes de L :

∂p g ∆p g
=− (XIII-16)
∂x L

Le profil de vitesse peut donc s’écrire :

1 ∆p g 2  r2 
u (r ) = R 1 − 2  (XIII-17)
4µ L  R 

la distribution de contrainte est donnée par :

du r ∂p g
σ rx = µ = (XIII-18)
dr 2 ∂x

Le débit en volume vérifie la formule célèbre ci-dessous, connue sous le nom de « loi de POISEUILLE » :

πR 4 ∆ p g
qV = (XIII-19)
8µ L

On remarque que la loi de Poiseuille donne une relation de proportionnalité entre : le débit et le gradient de pres-
sion motrice. Si on appelle «conductance» le terme :

πR 4
σ= (XIII-20)

et que l’on note
∆ pg
α= (XIII-21)
L

le gradient de pression motrice de l’écoulement, on a la relation :

Mécanique_des_fluides_2A_MP.doc 25 21/11/2007
qV
qV = ασ ou : α= (XIII-22)
σ

On appelle conductance le terme σ .

La vitesse moyenne de l’écoulement est :

1 R R2 ∆ pg
Vq =
A 0 ∫
u (r ) 2πr dr =
8µ L
(XIII-23)

Le point où la vitesse est maximale se situe sur l’axe du tube, soit en r = 0 , lorsque du / dr = 0 :

du R 2 ∆p g
r=0 ⇒ = 0 ⇒ u = u max = = 2 Vq (XIII-24)
dr 4µ L

Dans un écoulement laminaire établi dans une conduite de section circulaire, la vitesse maximale est deux fois
plus grande que la vitesse moyenne.

XIV- Ecoulement dans les conduites en charge.


XIV-1 Profils des vitesses dans les tuyaux
Dans un écoulement laminaire complètement développé, le profil des vitesses vérifie la loi de Poiseuille :

1 ∆p g 2  r2  u r2
= 1− 2
u (r ) = R 1 − 2  ou : (XIV-1)
4µ L Vq R
 R 

la vitesse est maximale au centre et nulle sur la paroi.

Dans un écoulement turbulent stationnaire, le profil des vitesses peut être représenté par la formule empirique :

1/ n
u  r
= 1 −  (XIV-2)
Vq  R 

où l’exposant n est fonction du nombre de REYNOLDS que l’on a déjà défini plus haut.

Cette formule empirique est appelée « loi puissance ». On peut préciser quelques valeurs de l’exposant n :

Re = 4 × 103 , n=6
Re = 1,1 × 105 , n=7 (XIV-3)
Re = 3,2 × 10 6 , n = 10

Le rapport entre la vitesse moyenne et la vitesse au centre du tube est donné par :

u 2 n2
= (XIV-4)
Vq (n + 1)(2n + 1)

On voit que si n augmente, ce rapport augmente.


Dans la figure qui suit, on a représenté différents profils de vitesses pour l’écoulement laminaire et un écoule-
ment turbulent dont le nombre de Reynolds est de 4000.

Mécanique_des_fluides_2A_MP.doc 26 21/11/2007
Conduite

r/R Laminaire (Poiseuille)


umax
Vq
Turbulent (loi puissance)

XIV-2 Répartition des contraintes de cisaillement (ou frottement) dans les écoulements
établis (cas des conduites).
On fait les hypothèses suivantes :

-Ecoulement horizontal. La pente du tube est nulle, donc : dZ / dx = 0 . La pesanteur n’intervient pas dans le
mouvement du liquide : dp g / dx = d ( p + ρgZ ) / dx = dp / dx + ρgdZ / dx = dp / dx .

-Ecoulement stationnaire

-Fluide incompressible

-Ecoulement complètement développé (établi)

L’écoulement étant stationnaire la projection des forces extérieures au volume de contrôle suivant l’axe x est
nulle :

Conduite

τr y
p M p
1 2
x

dx
r r
∑F ⋅x = 0
r r (XIV-5)
 ∂p dx  2  ∂p dx  2
∑ F ⋅ x = ( p1 − p2 )πr 2 + τ r 2πrdx =  p ( M ) −
 ∂x 2 
πr −  p ( M ) +
 ∂x 2 
πr + τ rx 2πr dx = 0

D’où :

r ∂p
τ rx (r ) = − (XIV-6)
2 ∂x

En particulier, sur la paroi du tube ( r = R ) ,

R ∂p
τ rx = − (XIV-7)
2 ∂x

Mécanique_des_fluides_2A_MP.doc 27 21/11/2007
XIV-3 Bilan d’énergie dans un tube : PERTE DE PRESSION TOTALE et PERTE DE
CHARGE
Considérons une conduite horizontale comportant un changement de direction ou « coude ».
Si on applique le premier principe de la thermodynamique au volume de contrôle D limité par la surface S, on a :

P + PS + Pfrottement + Pautre =

∂t
{∫ ρ(M ) e(M ) dω +∫ [e(M *) + p(M *)v(M *)] ρVr(M *) ⋅ dA}
M ∈D M *∈S
(XIV-8)

où le premier membre représente les puissances « produites » (au sens algébrique), et :

1
e = u + V 2 + gZ (XIV-9)
2

u est l’énergie mécanique par unité de masse

Si on fait les hypothèses suivantes :

(i) PS = Pfrottement = Pautre = 0


(ii) écoulement stationnaire de fluide incompressible
(iii) pression et énergie interne uniforme dans le volume de contrôle,

le bilan d’énergie précédent se réduit à :

 1 2   1 2  δQ
 p1 + ρgZ1 + 2 ρ α1 V q 1  −  p2 + ρgZ 2 + 2 ρ α 2 V q 2  = ρ(u1 − u 2 ) − ρ δm (XIV-10)
   

où :

α est le coefficient d’énergie cinétique


Vq est la vitesse moyenne
 1 2 
 p1 + ρgZ1 + 2 ρ α1 V q 1  est l’énergie mécanique totale par unité de volume dans la section S1
 
 1 2 
 p2 + ρgZ 2 + 2 ρ α 2 V q 2  est l’énergie mécanique totale par unité de volume dans la section S2
 
δQ
ρ(u1 − u 2 ) − ρ est la différence d’énergie mécanique totale entre les sections S1 et S2.
δm

On appelle perte de pression totale le terme :

δQ
∆Pt1−2 = ρ(u1 − u 2 ) − ρ (XIV-11)
δm

et perte de charge le terme équivalent exprimé en hauteur de fluide :

∆Pt1−2
∆H1−2 = (XIV-12)
ρg

On peut donc écrire :

 1 2   1 2 
 p1 + ρgZ1 + 2 ρ α1 V q 1  −  p2 + ρgZ 2 + 2 ρ α 2 V q 2  = ∆ Pt1−2 (XIV-13)
   

Mécanique_des_fluides_2A_MP.doc 28 21/11/2007
On remarquera que l’énergie mécanique totale dans une section est écrite à l’aide de la vitesse moyenne V q dans
cette section alors même que le profil des vitesses n’est pas « plat » (voir les formules donnant les vitesses). Ceci
n’est possible que grâce à l’introduction du coefficient d’énergie cinétique α qui « corrige » la valeur moyenne
de l’énergie cinétique. Il n’y a pas lieu, par contre, d’introduire de coefficient de pression pour le terme p + ρgZ
puisque celui est constant dans le plan de la section droite.

XIV-4 Perte de charge régulière et perte de charge singulière


La perte d’énergie mécanique totale dans un tuyau se décompose en deux parties : la perte de charge régulière
qui se produit tout le long de l’écoulement et les pertes de charge singulières qui sont localisées en des endroits
précis du conduit (coude, valve, raccords, etc..)
Pour tenir compte de ceci, on écrit :

∆ Pt1−2 = (∆ Pt1−2 )rég + (∆ Pt1−2 )sing (XIV-14)

XV – Calcul des pertes de charge.


XV-1 Coefficient de perte de charge régulière
Quel que soit le régime de l’écoulement (laminaire ou turbulent développé), la perte de charge peut être mise
sous la forme suivante :

L 1 2
∆ Pt = Λ ρV (XIV-15)
Dh 2 q

C’est la formule dite de DARCY-WEISBACH.


Cette formule est obtenue par analyse dimensionnelle, en appliquant le théorème de VASHY-BUCKINGHAM.
Elle donne la chute de pression totale se produisant entre deux sections distantes de la longueur L, et en écoule-
ment complètement développé.

Dh est le diamètre hydraulique ; il est défini à partir de l’aire A de la section droite et du périmètre mouillé
χ par :

4A
Dh = (XIV-16)
χ

On appelle coefficient de perte de charge régulière le coefficient sans dimension Λ ; c’est une fonction du
nombre de Reynolds et de la rugosité ε intérieure du tube (si l’écoulement est turbulent):

ρV q D h
Λ = f ( Re, ε ) = f × Re k , Re = (XIV-17)
µ

Cette fonction de Re est notamment donnée dans le diagramme de MOODY. On peut citer aussi les travaux de
NIKURADSE concernant la détermination de cette fonction.

La détermination de cette fonction consiste à trouver le coefficient f et l’exposant k du nombre de Reynolds.

On appelle coefficient de forme le coefficient sans dimension f. Il ne dépend, pour un régime d’écoulement
donné, que de la forme de la section droite du tube.

Mécanique_des_fluides_2A_MP.doc 29 21/11/2007
Exercice : Etablir la formule de Darcy-Weibach pour un écoulement établi en supposant que la chute de pression totale ∆ Pt
peut-être écrite comme une fonction de : D h , µ , L , V q , pour un tube à paroi lisse, fonction telle que :
∆Pt = f ( Dh , µ, ρ, L, Vq ) = C Dha µ b ρ c Ld V qe , (1)
où C est une constante sans dimension.

Il s’agit donc de trouver les exposants a,b, c, d, e, de telle sorte que la dimension de C Dh µ ρ L V q soit identique à
a b c d e

celle de ∆Pt , soit : [∆ Pt ] = ML−1T −2 .


Les dimensions de chacun des termes influençant la perte de charge sont :
[Dh ] = L ; [µ ] = ML−1 T −1 ; [ρ] = ML−3 ; [Vq ] = LT −1 ; [µ ] = ML−1 T −1 .
En exprimant que la formule (1) est homogène du point de vue dimensionnel :
(
ML−1T −2 = (L ) a ML−1T −2 ) (ML ) (L) (LT )
b −3 c
.d −1 e

Expérimentalement on montre que dans un écoulement établi, la perte de charge est proportionnelle à la longueur : on peut
donc écrire : d = 1 .
En développant puis en identifiant respectivement les exposants de chaque grandeur (Masse, Longueur, Temps), on a :
 b+c =1 a = e − 3
 
a − b − 3c − e = −2 ⇒  b = −e − 2
 − b − e = −2 c = e −1
 
a, b, c, sont paramétrés ici par e. D’où :
∆Pt = C D he−3 µ − e+ 2 ρ e−2 L V qe .
Faisons apparaître dans la perte de charge le nombre de REYNOLDS de l’écoulement ainsi que l’énergie cinétique
moyenne ρVq2 / 2 :
∆Pt = C D eh − 3 µ − e + 2 ρ e − 2 L V qe
 ρ e − 2Vqe − 2 D eh − 2  −1 1 L 1
= 2C  e−2  L Dh ρVq2 = 2C Re e − 2 ρVq2
 µ  2 Dh 2

En posant :

f = 2C Coefficient de forme
k = e−2 Exposant dépendant du régime d’écoulement
Λ = f ⋅ Re k
Coefficient de perte de charge

on obtient :
L 1
∆ Pt = Λ ρVq2
Dh 2

XV-1-1 Cas de l’écoulement laminaire dans une conduite de section circulaire cons-
tante : loi de Poiseuille
L’écoulement de Poiseuille possède une solution exacte ; il est donc possible d’obtenir théoriquement le coeffi-
cient de forme directement à partir du profil des vitesses u.

En effet,

∆Pt 1
Λ= D (XV-1)
L h ρVq2

L’intégrale de u figure dans la vitesse moyenne, puisque :

Mécanique_des_fluides_2A_MP.doc 30 21/11/2007
1
Vq = ∫
u ( M )dA
A M ∈S
(XV-2)

Les courbes expérimentales montrent que si l’écoulement est laminaire développé,

k = −1 (XV-3)

On a donc, en se rappelant que Vq = qV / A , avec A = πR 2 et χ = 2πR :

∆Pt 2 ρVq Dh ∆Pt 2 ∆Pt 32 A3 ∆Pt 32π 3 R 6


f = Re Dh = Dh = = (XV-4)
L ρVq2
µ L ρVq2
L µχ qV
2
L µχ 2 qV

On exprime le débit grâce la loi de Poiseuille en utilisant :

∆Pt πR 4
qV = ασ = (XV-5)
L 8µ

D’où :

∆Pt 32π 3 R 6
f = (XV-6)
L ∆Pt πR 4
µ 4π 2 R 2
L 8µ

Après simplifications, on obtient la valeur entière :

f = 64 (XV-7)

On notera bien que c’est une valeur exacte (et non déduites de mesures).
L’écoulement entre deux plaques fixes parallèles possède aussi une solution exacte (écoulement de Couette) ;
son coefficient de forme vaut : f = 96 .

XV-1-2 Cas des conduites cylindriques de section non circulaires


En écoulement laminaire, et lorsque la section n’est pas circulaire, (cas des échangeurs de chaleur…), on peut
calculer le coefficient de forme à condition que l’on ait au préalable déterminé le profil des vitesses u dans la
section droite.

Si on défini les variables sans dimension suivantes : ( µ viscosité dynamique du fluide, α est le gradient de
pression motrice, a est une longueur caractéristique de la section droite, par exemple la longueur du côté si la
section est rectangulaire)

µ
u* = u (vitesse normalisée)
αa 2
x* = x / a , y* = y / a , z* = z / a (coordonnées normalisées)

S devient S* (section en coordonnées normalisées)

l’équation du mouvement permanent du fluide suivant l’axe du tube se ramène à une équation de Poisson :

∆u* = −1 dans la section de passage,


u* = 0 , sur le bord de la section de passage S* (condition d’adhérence sur la paroi)

Mécanique_des_fluides_2A_MP.doc 31 21/11/2007
La résolution de cette équation de Poisson est possible par des méthodes numériques (différences finies ou élé-
ments finis si le contour est «complexe»). Le profil des vitesses u* étant connu, on peut trouver par intégra-
tion la conductance (en valeur adimensionnelle) :

µ
σ* =
a4 ∫
σ = u * ( M ) dy * dz *
M ∈S *
((XV-8)

D’où la valeur du coefficient de forme (en notant A * l’aire de S* et χ* = χ / a le périmètre normalisé de S*)

32 A *3
f = (XV-9)
σ * χ *2

XV-1-3 Cas de l’écoulement turbulent établi.


Pour les écoulements turbulents on doit utiliser des résultats expérimentaux pour déterminer le coefficient de
perte de charge Λ en fonction du nombre de REYNOLDS Re. Ces résultats sont présentés dans un diagramme
dit de MOODY, tel que celui de la figure ci-dessous qui s’applique aux conduites circulaires∗.

0,12
0,10 D
0,08 ε ε=
0,06 1 30 t

Λ
0,05 61
0,04 Zone de transition 4 120
0,03 222

0,02
2 504
1014

0,01
3
3 4 5 6
10 10 10
4 6 8 10 2 4 6 8
Vq Dh
Re =
ν
La « rugosité » est notée ε, c’est le rapport entre le diamètre intérieur et la taille moyenne, t, des aspérités :
ε = D / t . C’est un paramètre important pour l’écoulement turbulent.

Les différentes courbes correspondent aux formules suivantes :

64
1→ Loi de POISEUILLE : Λ= ;
Re
0,316
2→ Loi de BLASIUS, conduites lisses : Λ= ;
Re1/4
3→ Formule de KARMAN-PRANDTL :
1
Λ
[
= 2 log10 Re Λ − 0,8]
1 D
4→ Formule de KARMAN-PRANDTL : = 2 log10 + 1,14
Λ t


Les résultats tracés ici sont ceux de NIKURADSE.

Mécanique_des_fluides_2A_MP.doc 32 21/11/2007
En résumé, si on veut calculer les pertes de charge régulières d’un écoulement en régime turbulent, on doit sui-
vre la démarche suivante :

(a) Evaluer le nombre de Reynolds par la formule : Re = ρVq Dh / µ ; en régime turbulent, Re >> 2000 .
(b) Evaluer la rugosité intérieure ε de la conduite
(c) En déduire le coefficient de perte de charge Λ grâce au diagramme de MOODY ou à des formules empiri-
ques fournies dans les ouvrages d’hydraulique.
(d) Calculer le perte de charge, soit en unité SI de pression (Pa), soit en hauteur équivalent de liquide, en appli-
L 1
quant : ∆Pt = Λ ρV 2
Dh 2 q
Si la conduite est lisse, (rugosité nulle) on se rappellera que la loi de BLASIUS est applicable ( f = 0,316 et
k = −1 / 4 ) pour le régime turbulent, soit :

Λ = 0,316 × Re −1/ 4 (XV-10)

XV-2 Coefficient de perte de charge singulière


La perte de charge singulière est exprimée en fraction de l’énergie cinétique du fluide à l’entrée de la singularité.
On utilise la formule suivante si cette perte singulière est exprimée en unité de pression :

1
∆ Ptsing = K ρV 2 Homogène à une pression ( ML−1T −2 , Pa en unité SI) (XV-11)
2 q

ou la formule suivante si cette perte singulière est exprimée en hauteur équivalent de fluide :

1 2
∆ H sing = K V Homogène à une longueur (L, mètre en unité SI) (XV-12)
2g q

Le coefficient K est appelé coefficient de perte charge singulière ; il est sans dimension. Ce coefficient est par-
fois noté ξ .
Le problème pour l’ingénieur est alors de savoir déterminer le coefficient K.

De nombreuses études expérimentales ont été réalisées et leurs résultats sont publiés dans des bases de données
pour l’hydraulique. On peut ainsi obtenir pour la plupart des singularités usuelles (coudes, réductions, vannes,
raccords, etc…) la valeurs du coefficient K en fonction des caractéristiques géométriques des singularités et du
nombre de Reynolds Re .

Il est toujours possible de déterminer soi-même K. Pour cela il faut mesurer simultanément la perte de charge
entre l’entrée et la sortie de la singularité et la vitesse moyenne en entrée (ou le débit qV ), puis tracer le graphe
de ∆Ptsing en fonction de ρVq2 / 2 ; le coefficient cherché est alors la pente du nuage de points.
( Pa )
sing
∆ Pt

Droite de pente K

1 ρ V 2 ( Pa )
q
2

Mécanique_des_fluides_2A_MP.doc 33 21/11/2007
Il faut noter que pour le même nombre de Reynolds, des singularités « homothétiques » auront le même coeffi-
cient K.
Par exemple, pour le même nombre de Reynolds Re, deux coudes « droits » (angle de 90°) ayant le même rap-
R
port rayon de courbure/diamètre (soit 0 ) auront le même coefficient K.
D

R
0

ρ Vq D
Vq Re = µ

XVI – Quelques éléments sur la théorie de la couche limite


XVI-1 Introduction
Les écoulements de couche limite sont caractéristiques des fluides visqueux en mouvement relatifs près d’une
paroi solide.
La couche limite est la région du fluide en mouvement proche d’une paroi solide. Cette zone est particulièrement
sensible aux effets de la viscosité.
Une zone laminaire débute au bord d’attaque ( x = 0 ) et son épaisseur augmente à mesure que l’écoulement
avance le long de la paroi.
Il se crée alors une zone de transition laminaire-turbulent après laquelle l’épaisseur de la couche limite aug-
mente.
Il existe dans la zone turbulente une sous couche visqueuse dans laquelle les effets visqueux (liés aux
frottements entre couches de fluides animées de vitesses différentes) sont prédominants devant les effets
d’inertie (liés à la mise en mouvement des masses fluides).

On représente le développement de la couche limite sur une plaque plane dans la figure ci-dessous :

δ(x) épaisseur de la couche limite

y Zone de Transition

x Sous couche visqueuse

Zone laminaire Zone turbulente

On appelle épaisseur de la couche limite la distance δ entre le fluide et paroi à partir de laquelle la vitesse du
fluide est égale à 99% la vitesse U 0 du fluide libe (non perturbé par la paroi).

Mécanique_des_fluides_2A_MP.doc 34 21/11/2007
y
Fluide libre U0

δ(x) épaisseur de la couche limite


y

u( xM ,y) Couche limite


x
xM

On appelle épaisseur de déplacement δ * la longueur du déplacement qu’il faudrait appliquer à la paroi solide
pour conserver le même débit en masse à un écoulement de fluide parfait hypothétique.

On a, en fluide incompressible :


 u ( xM , y ) 
δ* = ∫
0
1 − U
 0
 dy

(XVI-1)

XVI-2 Equation intégrale


Admettant que la pression motrice s’identifie à la pression statique p dans la section droite et que de plus p est
constante dans cette section, la résultante des forces suivant la direction x s’écrit :

dF X = σ( x) − σ( x + dx) + σ(δ) − τ p dx (XVI-2)

Frontière de la couche limite Volume de contrôle


Couche limite

y
δ(x) δ(x)+dδ
x

Plaque

Mécanique_des_fluides_2A_MP.doc 35 21/11/2007
σ(δ)
y
σ(x) σ(x + dx)

x τp dx

x dx x+dx
Où τ p est la contrainte de frottement sur la paroi.

En considérant une unité de longueur perpendiculairement à la figure, on a comme force :

 dp 
σ( x) = p × δ , σ( x + dx) =  p( x) + dx × (δ + dδ )
 dx 
 1 dp 
σ (δ ) =  p ( x ) + dx × dδ
 2 dx 

 1 dp   1 dp 
dFx = p × δ −  p( x) + dx  × (δ + dδ ) +  p( x) + dx dδ − τ p dx (XVI-3)
 2 dx   2 dx 

Si l’on combine le flux de quantité de mouvement suivant x avec l’équation de conservation de la masse pour le
volume de contrôle choisi, on obtient la forme générale de l’équation intégrale de quantité de mouvement :

dp d  δ 2  d  δ 
−δ
dx
− τ p dx = 
dx  0 ∫
ρ u dy  − U 0
 ∫
dx  0
ρ u dy 

(XVI-4)

Si l’on veut évaluer l’épaisseur de couche limite pour un écoulement laminaire de fluide incompressible au des-
sus d’une plaque plane, il faut :

(i) Déterminer le gradient de pression dp / dx en dehors de la couche limite (par une étude en écoulement
potentiel) ;
(ii) Se donner une expression analytique plausible pour le profil de vitesses u ( y ) dans la couche limite ;
(iii) Utiliser comme formule pour la contrainte visqueuse (frottement), la loi de Newton :

∂u
τp = µ
∂y

XVI-3 Ecoulement dans la couche limite d’une plaque plane


On prend les hypothèses suivantes :

L’écoulement extérieur possède une vitesse constante : U 0 = Cte ;


dp
La pression est constante : = 0.
dx
Le fluide est incompressible : ρ = Cte .

Mécanique_des_fluides_2A_MP.doc 36 21/11/2007
L’équation intégrale de quantité de mouvement se « simplifie » et conduit à :

d  δ  u2 u   2 d 
δ
u  u  
∫ ∫
2
τ p = −ρU 0   2 −  dy  = ρU  1 − U  dy 
dx  0
U 0 U 0   dx  0 U0  0  
 

Soit la variable sans dimension définie en x donné :

y
η= ⇒ dy = δ dη
δ

On obtient :

dδ 1
u  u 

2
τ p = ρU 0 1 − U  dy (XVI-5)
dx 0U 0  0 

C’est l’équation intégrale des quantités de mouvement pour un écoulement sans gradient de pression ;
cette expression de la contrainte de frottement est restreinte au cas suivant, rappelons-le :

(i) Ecoulement stationnaire sans forces de pesanteur (le fluide est dit « non-pesant ») ;
(ii) Ecoulement incompressible bidimensionnel ;
(iii) Ecoulement le long d’une plaque plane.

XVI-4 Equation intégrale pour un écoulement sans gradient de pression ( dp / dx = 0 )

A/ Ecoulement laminaire le long d’une plaque plane

On se donne la forme du profil de vitesses :

u ( y ) = A + By + Cy 2

c’est un polynôme de degré 2 en y. Il doit vérifier les conditions limites : adhérence sur la plaque (1ère condition),
raccordement avec l’écoulement extérieur de vitesse U 0 (2ème et 3ème condition) :

y=0 ⇒ u=0
y = δ ⇒ u =U0 (XVI-6)
du
y=δ ⇒ =0
dy

En injectant l’expression polynomiale ci-dessus dans la formule de τ p vue plus haut et en appliquant la loi de
Newton pour le membre de gauche, on obtient :

2µU 0
τp = (XVI-7)
δ

Par définition, le coefficient de frottement est donné par le rapport entre la contrainte de frottement et l’énergie
cinétique (ou pression dynamique) de l’écoulement extérieur :

τp
Cf = (XVI-8)
1
ρU 02
2

Dans le cas présent, on peut montrer que :

Mécanique_des_fluides_2A_MP.doc 37 21/11/2007
0,730
Cf = (XVI-9)
Re( x)

Où Re(x) est le nombre de REYNOLDS à l’abscisse x :

ρU 0 x
Re( x) = (XVI-10)
µ

L’épaisseur de la couche limite est donnée par :

δ 5,48
= (XVI-11)
x Re( x)

L’épaisseur de déplacement est donnée par :

δ* 1,835
= (XVI-12)
x Re( x)

B/ Ecoulement turbulent le long d’une plaque plane

On se donne la forme du profil de vitesse :

1/ 7
u  y
=  (XVI-13)
U0  δ 
la contrainte de frottement, déterminée expérimentalement (PRANDTL), vérifie :

1/ 4
 µ 
τ p = 0,0225 ρU  2
 (XVI-14)
 ρU 0 δ 
0

L’épaisseur de la couche limite est donnée par :

δ 0,37
= (XVI-15)
x Re( x) 1 / 5

Le coefficient de frottement sur la paroi est donné par :

0,0577
Cf = 1/ 5
, lorsque : 5 × 10 5 < Re( x ) < 10 7 (XVI-16)
Re( x)

XVI-5 Coefficients de frottement des plaques lisses en écoulement turbulent

Pour des nombres de Reynolds supérieurs à 10 5 , le coefficient de frottement C f est exprimé par une formule
déduite de résultats expérimentaux :

−1 / 5
C f = 0,074 Re L , lorsque : 5 × 10 5 < Re L < 10 7 (XVI-17)

0,455
Cf = , lorsque : Re L > 10 7 (d’après H. SCHLICHTING) (XVI-18)
(log Re L )2,58

Mécanique_des_fluides_2A_MP.doc 38 21/11/2007

ρU 0 L
Re L =
µ
est le nombre de Reynolds défini à partir de la longueur de la plaque.

La figure ci-dessous représente le coefficient de frottement C f du fluide léger sur la paroi d’une plaque plane
lisse en fonction du nombre de Reynolds Re L .

6 E-3

4 E-3
Cf

1 2
2 E-3 3

1 E-3
1 E5 2 5 1 E6 1 E7 1 E8 1 E9

Re L
Courbe 1 : Ecoulement laminaire.
Courbe 2 : Ecoulement de transition laminaire-turbulent
Courbe 3 : Ecoulement turbulent.

Mécanique_des_fluides_2A_MP.doc 39 21/11/2007

Vous aimerez peut-être aussi