Mécanique Des Fluides 2 Eme Année
Mécanique Des Fluides 2 Eme Année
Mécanique Des Fluides 2 Eme Année
2e année
Mesures Physiques
Mécanique des Fluides
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 .
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 .
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
r r
∫ ρ(M ) F (M )dω + ∫ − p(M *)n(M *)dA = 0
M ∈D M *∈S
(I-4)
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
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.
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
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
∂p ( M ) ∂p( M ) r ∂p( M ) r
grad p( M ) = x+ y+ z (I-8)
∂x ∂y ∂z
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)
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.
r p( M )
F ( M ) = grad (I-12)
ρ( M )
p( M )
∫
M ∈D
ρ( M ) grad
ρ( M )
dω + ∫ − p(M *) n(M *)dA =0
M *∈S
(I-13)
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
U ( M ) = gz ( M ) + U 0 (II-3)
Puisque la masse volumique est constante en tout point de D, ρ(M ) = ρ , l’équation fondamentale SDF-3
s’écrit :
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)
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
V= ∫ dω(M )
M ∈K
(II-6)
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.
Mécanique_des_fluides_2A_MP.doc 7 21/11/2007
dp ( M )
= −ρ( M ) g (III-1)
dz
ρ( M ) = ρ( M 0 )[1 + χ( p ( M ) − p( M 0 ) )]
dp ( M )
ρ( M 0 )[1 + χ( p ( M ) − p( M 0 ) )] g = − (III-2)
dz
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
p(M ) T (M 0 )
ρ( M ) = ρ ( M 0 ) (III-5)
p(M 0 ) T (M )
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 )
Mécanique_des_fluides_2A_MP.doc 9 21/11/2007
DYNAMIQUE DES FLUIDES
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
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
Si l’écoulement est stationnaire les trajectoires et les lignes de courant sont confondues.
Il s’écrit symboliquement :
r r r
rot V ( M , t ) = ∇ ∧ V ( M , t ) (VI-1)
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.
∫
∂ρ( 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
∂ρ( 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-
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.
S1 n*
SL surface latérale
n1
q
V1
S2
q
V2
n2
on a donc :
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.
On applique ci-dessous le PFD à un domaine D de fluide que l’on suit dans son mouvement.
M est un point 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)
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
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).
∂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 )
∂V ( M , t ) r r
+ (rot V ( M , t ) ∧ V ( M , t ) = −grad H ( M , t ) (IX-10)
∂t
∂V ( M , t )
=0 (IX-11)
∂t
De plus comme
r
dH = grad H ( M , t ) ⋅ d l (IX-12)
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.
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 :
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)
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.
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
p g ( M ) = p ( M ) + ρg Z ( M ) (XII-1)
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
,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…
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.
U fonction du temps
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
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)
µ
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 :
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
1 ∆p g 2 r2
u (r ) = R 1 − 2 (XIII-17)
4µ L R
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)
8µ
et que l’on note
∆ pg
α= (XIII-21)
L
Mécanique_des_fluides_2A_MP.doc 25 21/11/2007
qV
qV = ασ ou : α= (XIII-22)
σ
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.
1 ∆p g 2 r2 u r2
= 1− 2
u (r ) = R 1 − 2 ou : (XIV-1)
4µ L Vq R
R
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)
Mécanique_des_fluides_2A_MP.doc 26 21/11/2007
Conduite
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
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
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)
1
e = u + V 2 + gZ (XIV-9)
2
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ù :
δQ
∆Pt1−2 = ρ(u1 − u 2 ) − ρ (XIV-11)
δm
∆Pt1−2
∆H1−2 = (XIV-12)
ρg
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.
L 1 2
∆ Pt = Λ ρV (XIV-15)
Dh 2 q
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.
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
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
Mécanique_des_fluides_2A_MP.doc 30 21/11/2007
1
Vq = ∫
u ( M )dA
A M ∈S
(XV-2)
k = −1 (XV-3)
∆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µ
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 .
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)
l’équation du mouvement permanent du fluide suivant l’axe du tube se ramène à une équation de Poisson :
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
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.
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 :
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 = µ
On représente le développement de la couche limite sur une plaque plane dans la figure ci-dessous :
y Zone de Transition
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
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)
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.
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
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
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.
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
Mécanique_des_fluides_2A_MP.doc 37 21/11/2007
0,730
Cf = (XVI-9)
Re( x)
ρU 0 x
Re( x) = (XVI-10)
µ
δ 5,48
= (XVI-11)
x Re( x)
δ* 1,835
= (XVI-12)
x Re( x)
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
δ 0,37
= (XVI-15)
x Re( x) 1 / 5
0,0577
Cf = 1/ 5
, lorsque : 5 × 10 5 < Re( x ) < 10 7 (XVI-16)
Re( x)
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
où
ρ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