Cinematique Et Dynamique Des Fluides Uee PDF

UEE PA6 – IFMAT et MASTER IMCE
UNITE D’ENSEIGNEMENT
AERO-HYDRODYNAMIQUE INTERNE DES MACHINES
TOME I
CINEMATIQUE ET DYNAMIQUE DES FLUIDES
Septembre 2008 R. REY – Professeur

Robert REY - Professeur
-Chapitre
A- A
PROPRIETES
PROPRIETES CINEMATIQUES
CINEMATIQUES DES
ECOULEMENTS
DES ECOULEMENTS

Ch.A – Propriétés cinématiques des écoulements
I) GENERALITES p. 5
1.1. Définitions - Propriétés des lignes de courant
1.2. Accélération locale - Rotationnel
1.3. Equation de continuité
1.4. Ecoulement irrotationnel - Fonction potentiel
1.5. Circulation du vecteur vitesse
II) ECOULEMENTS BIDIMENSIONNELS p. 15
2.1. Ecoulements plans en coordonnées cartésiennes
2.2. Ecoulements bidimensionnels de révolution
2.3. Méthode simplifiée de tracé des lignes de courant
2.4. Solutions analytiques de l’équation de Laplace
2.5. Application de la superposition des écoulements
CINEMATIQUE ET DYNAMIQUE DES FLUIDES 2

- PROPRIETES CINEMATIQUES DES ECOULEMENTS -
La mécanique des fluides constitue depuis plus d'un siècle le champ d'application des
plus grands noms des mathématiques et de la physique, parmi lesquels :
. au 18ème siècle : d'Alembert, Euler, Lagrange, Laplace ;
. au 19ème siècle : Cauchy, Dirichlet, Poiseuille, Rankine, Reynolds, Riemann, Stokes ;
. et au siècle dernier : Blasius, Joukowski, Karman et Prandtl.
Tour à tour domaine d'excellence du fluide parfait, du fluide réel puis de la physique
expérimentale, la mécanique des fluides subit aujourd'hui une grande mutation par
l'introduction des méthodes numériques qui permettent, grâce au développement des
moyens de calcul, de traiter conjointement les équations cinématiques, dynamiques et
énergétiques.
Seuls quelques aspects intrinsèques tels que les couches limites, la turbulence, les
écoulements secondaires, les régimes transitoires, limitent, à des degrés divers, la fiabilité
des modélisations.
Malgré ces progrès très importants, la détermination des vitesses autour d'un obstacle,
d'un profil portant, est toujours une question d'actualité.
Pour le concepteur, les problèmes posés demeurent inchangés et se formulent en termes

identiques. Ils se ramènent pour la plupart à la résolution des deux problèmes types : le
problème direct et le problème inverse.
Problème direct :
Pour un profil ou une succession de profils donnés (grille d'aubes annulaire, cylindrique,
etc...) immergé dans un écoulement défini en amont de l'obstacle, quelle est la
distribution de vitesses au voisinage des parois ?
L'application a posteriori des équations dynamiques permettra d'évaluer les pressions

locales puis les efforts mis en jeu.
C'est le problème qui se traite dans le cadre de l'évaluation des performances d'une
machine existante.

Problème inverse :
C'est la démarche adoptée lors du projet d'une nouvelle machine. Pour une cinématique
d'écoulement définie à l'avance, à savoir : distribution de vitesse à l'amont et à l'aval du
profil ou de la grille, champ de vitesses imposé à la paroi, quelle est la géométrie du ou
des profils correspondants ?
La définition a priori du champ de vitesses est souvent inspirée des considérations

suivantes :
- obtention des performances attendues,
- maîtrise spatiale des échanges énergétiques en vue de l'optimisation du rendement, du
bruit, de la tenue à la cavitation, de la stabilité en débit partiel.
Pour simplifier le traitement du problème direct et du problème inverse, l'écoulement

principal est souvent ramené à une juxtaposition d'écoulements bidimensionnels
élémentaires compatibles entre eux. La résolution de l'aspect purement cinématique n'est
rendue possible que par l'introduction d'un certain nombre d'hypothèses simplificatrices :
- fluide non visqueux,
- écoulement irrotationnel.

I) GENERALITES
1.1. Définitions - Propriétés des lignes de courant
Dans le domaine des turbomachines, la plupart des écoulements et leurs équations

associées, sont traités en coordonnées cartésiennes (x, y, z) ou en coordonnées
cylindriques (r, θ, z).
x r direction radiale
repères y ou θ direction tan gentielle
z z direction axiale
→
Dans ces repères la vitesse C en un point quelconque M aura pour composantes :
dx dr
u= Cr =
dt →
dt
→ dy r dθ
C= v= ou C = Cθ ou Cu =
dt dt
dz
dz Ca =
w= dt
dt
Relations entre coordonnées :
 x = r cos θ

 y = r sin θ (1)
z = z

et
 r = x ² + y²

 y
 tg θ =
 x (2)
 z = z

Figure 1 : Définition du tube de courant
La description lagrangienne de l'écoulement consiste à observer, en fonction du temps, le

déplacement d'une particule donnée et identifiée ; son déplacement s'effectue sur une
trajectoire. Inversement, la description eulérienne se place en un point de l'écoulement et
observe les variations de la vitesse en fonction du temps (plusieurs particules se
succèdent en ce même point).
Si l'on observe au même instant l'ensemble des points de l'écoulement pendant une durée
assez courte (photographie) on peut définir les lignes de courant.
→
La courbe (ψ) est une ligne de courant si en tout point M appartenant à (ψ), la vitesse C
est tangente à (ψ). On observera donc, pour l'ensemble des points :
dx dy dz
= = (3)
u v w
Les composantes du vecteur vitesse dépendent à la fois des coordonnées du point M et du

temps.
Lorsque l'écoulement est permanent (ou stationnaire), les composantes du vecteur vitesse
sont indépendantes du temps. Les trajectoires et les lignes de courant sont alors
confondues. La majorité des écoulements industriels sont instationnaires : l'écoulement
dans une machine, l'écoulement turbulent en conduite. Cependant, il est souvent aisé et
suffisamment précis de considérer qu'ils sont permanents en moyenne.
L'enveloppe des lignes de courant s'appuyant sur un contour fermé forme un tube de
courant. Aucun débit ne traversant l'enveloppe du tube de courant, le débit massique Q
circulant dans le tube reste constant, il est calculable par :
→ →
q m = ∫∫ ρ C. n dS (4)
s
avec : dS : élément de surface

→
n : vecteur unitaire normal à dS
ρ : masse volumique du fluide
q m = ∫∫ ρ ( u dy dz + v dx dz + w dx dy) (5)
s
On définit également le débit volume comme le flux à travers la surface S.
qv = ∫∫s u dy dz + v dx dz + w dx dy (6)
Dans une section normale à la ligne de courant (S1 et S2 de la figure 1), il est commode
d'introduire la vitesse moyenne spatiale (ou vitesse débitante) :
q m = ρ1 S1 C1 = ρ2 S2 C 2 (7)

ρ1 et ρ2 étant les masses volumiques moyennes dans les sections S1 et S2 (calculables à

partir de l'équation d'état du fluide).
1.2. Accélération locale - Rotationnel
→ →
Le vecteur accélération Γ est calculable localement à partir de la vitesse C :
→ d →
Γ = (C) (8)
dt
et ses trois composantes en coordonnées cartésiennes :
d
( u) = Γx
dt
→ d
Γ= ( v) = Γy
dt
d
( w ) = Γz
dt (9)
Sachant que dans le cas général u, v, w sont des fonctions de x, y, z et t on obtiendra la

première composante à partir de la différentielle totale d'une fonction de quatre variables
(annexe R8) :
∂u ∂u ∂u ∂u
du = . dx + . dy + . dz + .dt
∂x ∂y ∂z ∂t
du ∂u ∂u ∂u ∂u
et Γx = = .u + .v + .w +
dt ∂x ∂y ∂z ∂t
dx dy dz →
où l'on retrouve : u= ,v = , w= les composantes de C .
dt dt dt
→
Soit finalement les composantes de Γ :
∂u ∂u ∂u ∂ u
Γx = u +v +w +
∂x ∂y ∂z ∂t
→ ∂v ∂v ∂v ∂ v
Γ= Γy = u +v +w + (10)
∂x ∂y ∂z ∂t
∂w ∂w ∂w ∂w
Γz = u +v +w +
∂x ∂y ∂z ∂t

Pour chaque composante, le dernier terme est l'accélération locale traduisant la non
permanence de l'écoulement, les autres termes correspondent à l'accélération convective
associée aux variations de vitesse dues aux variations de section de passage.
En coordonnées cylindriques, on peut calculer de la même façon les coordonnées Γr, Γθ,
Γz :
∂Cr ∂Cr ∂Cr Cu ² ∂Cr
Γr = Cr+ Cu + Ca − +
∂r r∂θ ∂z r ∂t
→ ∂Cu ∂Cu ∂Cu CrCu ∂Cu
Γ = Γθ = Cr + Cu + Ca + + (11)
∂r r∂θ ∂z r ∂t
∂Ca ∂Ca ∂Ca ∂Ca
Γz = Cr + Cu + Ca +
∂r r∂θ ∂z ∂t
On pourra retenir d'une façon générale la relation vectorielle (forme de Lamb) :
→
→ ∂C → C ² → → →
Γ = + grad − C Λ rot C
∂t 2 (12)
Nous verrons par la suite qu'il est important d'utiliser le vecteur rotationnel de la vitesse
→
C.
En coordonnées cartésiennes, celui-ci sera défini par :
∂w ∂v
− = 2Ω x
∂ y ∂z
→ → ∂u ∂w
rot C = − = 2Ω y (13)
∂x ∂y
∂ v ∂u
− = 2Ω z
∂ x ∂y
et en coordonnées cylindriques :
1 ∂Ca ∂Cu
− = 2 Ωr
r ∂θ ∂z
→ → ∂Cr ∂Ca
rot C = − = 2 Ωθ (14)
∂z ∂r
1 ∂ 1 ∂Cr
( r Cu ) − = 2 Ωz
r ∂r r ∂θ

→
On rappellera enfin la définition du vecteur tourbillon Ω :
→ 1 → →
Ω = rot C (15)
2
1.3. - Equation de continuité
Cette équation exprime la conservation de la masse contenue dans un volume élémentaire

dV (Lagrange) ou encore pour un élément fixe (Euler), que la diminution de masse
contenue dans le volume dV est égale à la masse ayant traversé la surface extérieure de
l'élément (débit sortant) :
En utilisant cette seconde configuration, on aura (figure 2) :
Masse totale entrant dans un parallélépipède ayant pour côtés dx, dy, dz :
ρ u. dy dz + v dx dz + w dx dy dt
Masse totale sortant au même instant :
 ∂( ρu )   ∂( ρv )   ∂( ρw ) 
 ρu + ∂x dx  dy dz dt +  ρv + ∂y dy dx dz dt +  ρw + ∂z dz dx dy dt
z
∂ (ρw )
ρw + . dz
∂z
ρu
∂ ( ρv )
ρv + dy
∂y
ρv
y
ρw
∂ (ρu )
ρu + dx
∂x
x
Figure 2 : Flux à la surface de l'élément fluide

La masse sortante est donc égale à :
 ∂( ρu ) ∂( ρv ) ∂( ρw ) 
dm =  + +
∂z 
. dx dy dz dt
 ∂x ∂y
La diminution de masse à l'intérieur de l'élément pendant la même durée dt:

∂ρ
dm = − . dt . dx dy dz
∂t
En égalant les deux valeurs de dm on obtient l'équation de continuité :
∂ρ ∂ ( ρ u ) ∂ ( ρv ) ∂ ( ρ w )
+ + + = 0 (16)
∂t ∂x ∂y ∂z
Ou encore sous forme plus compacte (annexe A4):
∂ρ →
+ div ( ρ C ) = 0 (17)
∂t
En coordonnées cylindriques, l'équation de continuité s'écrira :
∂ρ 1 ∂ ∂( ρCu) ∂( ρCa )
+
∂t r ∂r
( r ρ Cr ) +
r∂θ
+
∂z
= (18)
Cas particuliers :
Le premier terme s'annule dans les cas particuliers suivants:
- Ecoulement permanent:
ρ en un point est constant même si sa valeur peut varier d'un point à un autre.
- Ecoulement d'un fluide incompressible:

ρ reste constant quelles que soient les valeurs de la pression et de la température au sein
du fluide.

Dans ce cas, la relation (17) se ramène à :

→
div C = 0 (19)
1 ∂ ∂(Cu) ∂(Ca)
ou r Cr + + = 0 (20)
r ∂r r ∂θ ∂z
N.B. :
Quelle que soit la forme géométrique du volume V, l'expression de l'équation de
continuité est accessible en utilisant le théorème de Green de la divergence permettant
de remplacer la relation (4) par une intégrale de volume (annexe R7).
1.4. Ecoulement irrotationnel - Fonction potentiel
"L'écoulement est à potentiel des vitesses" ou encore le "champ des vitesses dérive d'un
potentiel" est la terminologie employée dans le cas où la relation vectorielle suivante est
vérifiée :
→ 
→
C = grad φ (21)
avec φ fonction potentiel dépendant des coordonnées au même titre que les composantes
→
du vecteur vitesse C .
On l'écrit également (annexe R3), en coordonnées cartésiennes :
∂φ
u =
∂x
→ ∂φ
C = v =
∂y
∂φ
w =
∂z (22)
∂φ
Cr =
∂r
→ ∂φ
C = Cu =
r∂θ
∂φ
Ca =
∂z (23)

Propriétés de la fonction potentiel φ :

L'écoulement à potentiel des vitesses présente la propriété remarquable de vérifier la
relation vectorielle :
→ → →
rot C = 0 (24)
En remplaçant dans le système (13) les termes u, v et w par leurs valeurs issues de (22),
on obtiendra :
∂² φ ∂² φ
− = 0
∂z ∂y ∂y ∂z
→ → ∂² φ ∂² φ
rot C = − = 0
∂x ∂z ∂z ∂x
∂² φ ∂² φ
− = 0
∂y ∂x ∂x ∂z
→ →
Les trois composantes du vecteur rot C sont nulles, on qualifie alors l'écoulement
d'irrotationnel.
Propriétés de l'écoulement irrotationnel :
D'après la relation (15), l'écoulement irrotationnel n'est le siège d'aucun vecteur

tourbillon. Il ne peut donc être utilisé que pour la description des écoulements en fluide
parfait ou par extension des écoulements en fluide réel exempts de décollement. Cette
propriété n'est pas applicable aux écoulements réels présentant :
• une forte incidence sur des obstacles minces (désadaptation)
• de grandes variations de direction ou de section de passage
On démontre qu'en tout point d'un écoulement irrotationnel, l'énergie totale est constante.
Lorsqu'il existe un échange d'énergie entre les particules fluides et le milieu environnant,
l'écoulement est rotationnel et l'énergie totale des particules varie d'une ligne de courant à
l'autre. C'est le cas de l'écoulement relatif dans les machines tournantes.
Réciproquement, l'écoulement irrotationnel est vérifié par l'introduction de la fonction
potentiel φ sous la forme d'une différentielle totale :
∂φ ∂φ ∂φ → →
dφ = dx + dy + dz = u dx + v dy + w dz = C . ds (25)
∂x ∂y ∂z

ou en coordonnées cylindriques :
∂φ ∂φ ∂φ → →
dφ = dr + rd θ + dz = Cr dr + Cu rd θ + Ca dz = C . ds (26)
∂r r∂θ ∂z
Les relations (25) et (26) seront développées plus loin.
Si l'on utilise la même fonction définie en (22) dans l'équation de continuité (16) en fluide
incompressible (ρ = cte), on obtiendra l'équation de Laplace :
∂²φ ∂²φ ∂²φ
∇ ²φ = + + = 0 (27)
∂x ² ∂y ² ∂z²
1.5. Circulation du vecteur vitesse
→
Par définition, la circulation Γ du vecteur vitesse C de A et B le long de la courbe (C) est
donnée par :
→ →
Γ = ∫ C . ds = ∫ C . cos α . ds (28)
AB AB
→→
dx
→
avec ds dS vecteur tangent en M à (C) : dS = dy
dz
Figure 3 : Définitions
Le produit scalaire entre les deux vecteurs permet de calculer Γ (annexe R2) :
Γ = ∫ ( u dx + v dy + w dz) (29)
AB

Dans le cas général, Γ dépend des points A et B et du chemin suivi pour aller de A à B.
Par contre, dans le cas où le terme (udx + vdy + wdz) est une différentielle totale exacte,
la circulation ne dépend que des bornes A et B. La fonction φ satisfait à cette condition
(annexe R8) :
Γ = ∫ dφ = φB − φA (30)
AB
Ainsi, la circulation est constante de A à B si l'écoulement est irrotationnel (ou

identiquement à potentiel de vitesses).
L'équation (30) montre que la circulation est nulle le long d'une courbe fermée.

II) ECOULEMENTS BIDIMENSIONNELS -
Une grande partie des écoulements réels, qu’ils soient internes ou externes, présentent un
caractère tridimensionnel. Au moins au stade de l’avant-projet, un certain nombre de
méthodes, dont la transformation conforme, permettent de décomposer ces écoulements
complexes en une juxtaposition d’écoulements bidimensionnels compatibles entre eux.
Toujours en fluide parfait, nous introduisons ici les propriétés de l'écoulement

bidimensionnel, propriétés particulièrement intéressantes dans le cadre de l'étude du
problème direct :
Pour une géométrie donnée le problème direct consiste à définir le champ de vitesses
correspondant. Cette étude est grandement facilitée par l'introduction de la fonction
potentiel φ et de la fonction de courant ψ, toutes deux harmoniques et satisfaisant
l'équation de Laplace. La fonction potentiel satisfaisant l'irrotationnalité de l'écoulement
et la fonction de courant l'équation de continuité. Cette représentation de l'écoulement en
coordonnées "naturelles" constitue le point de départ des maillages les plus sophistiqués.
Le maillage (φ,ψ) étant réalisé, le module et la direction de la vitesse sont accessibles en
tout point de l’écoulement.
On introduira ensuite les fonctions analytiques, à l'origine de la méthode des singularités.

Deux méthodes qui, grâce au développement des moyens de calcul, bénéficient
aujourd'hui d'un regain d'intérêt. On préférera l'une à l'autre des méthodes en fonction du
type de problème à traiter.
Figure 4 : Représentation des lignes de courant et des équipotentielles.

2.1. Ecoulements plans en coordonnées cartésiennes
On se place dans l'hypothèse d'un écoulement plan en fluide incompressible (ρ= cte) pour
lequel la vitesse en un point reste constante en module et direction en fonction de
l'épaisseur.
Les lignes de courant sont définies par l'équation (3) :
u v
= (40)
dx dy
→
u et v étant les composantes du vecteur vitesse C :
→
→ dS
C = (41)
dt
→
où d S est l'abscisse curviligne sur la ligne de courant :
→ dx
dS =
dy
→ →
D'où les composantes du vecteur vitesse C colinéaire à dS :
dx
→
u =
dt
C =
dy
v =
dt
Dans ces conditions, l'équation de continuité (16) se réduit à l'expression :
∂u ∂v
+ = 0
∂x ∂y
(42)
On suppose l'écoulement irrotationnel et on introduit la fonction potentiel φ à partir des

→
composantes de la vitesse C :
∂φ
u =
→ ∂x
C = (43)
∂φ
v =
∂y

On peut calculer φ à partir de sa différentielle totale (annexe A8) :
∂φ ∂φ
dφ = dx + dy = u dx + v dy (44)
∂x ∂y
dφ = C . ds (45)
∆φ
ou encore : C = (46)
∆s
En tout point, le module de la vitesse est donné par le taux de variation de la fonction φ
en fonction de l'abscisse curviligne (figure 4).
Fonction de courant Ψ :
Au même titre que la fonction potentiel, on introduit la fonction de courant ψ à partir des
composantes de la vitesse. Pour cela on posera :
∂ψ ∂φ
→ u= =
C = ∂y ∂x (47)
∂ψ ∂φ
v=− =
∂x ∂y
La propriété essentielle de cette fonction est de vérifier l'équation de continuité. En effet,

en remplaçant u et v par leurs valeurs dans l'équation (42), on obtient :
∂²ψ ∂²ψ
− = 0
∂x . ∂y ∂y . ∂x
Dans le cas d'un écoulement irrotationnel, la troisième composante de (13) donnera :
→ → ∂²ψ ∂²ψ
rot C = + = 0 (48)
∂x ² ∂y ²
ψ est donc une fonction harmonique.
Si l'on reprend l'équation (40) définissant les lignes de courant :
u dy - v dx = 0

et si on introduit ψ par l'intermédiaire des composantes u et v (47), on aura :
∂ψ ∂ψ
dy + dx = 0 (49)
∂y ∂x
On reconnaîtra pour le premier terme la différentielle totale de la fonction ψ :
dψ = 0 (50)
Montrant que sur une ligne de courant la fonction ψ est constante.
Propriétés importantes :
L'écoulement irrotationnel est décrit par la fonction φ et les lignes de courant par la
fonction ψ.
φ et ψ sont des fonctions des coordonnées du point considéré et sont, dans tout l’espace
fluide, orthogonales (figure 4).
On remarquera que les frontières de l'écoulement (parois) sont également des lignes de
courant.
→
On démontre enfin que le module de la vitesse C est égal en tout point à la variation
transversale de la fonction ψ :
∆ψ
C = (51)
∆n
En reprenant les relations (46) et (51) on observe qu'en tout point du maillage, on aura :
∆φ ∆ψ
C = = (52)
∆S ∆n
Dans les applications on admettra que C est la vitesse au centre de la maille (intersection
des pseudo diagonales).

Condition de glissement à la paroi :
Certaines lignes de courant (ψ1 et ψ5 de la figure 4) forment les frontières de

l'écoulement. Pour ces lignes de courant particulières d'équation S(x, y) = 0, la condition
→
de glissement sera observée lorsque le produit scalaire entre la vitesse C et la normale à
la paroi est nul :
∂S ∂S
u + v = 0 (53)
∂x ∂y
2.2. Ecoulements bidimensionnels de révolution
L'étude de l'écoulement de révolution peut se ramener à celle des écoulements plans. Il

est représenté dans le système de coordonnées (r, z) où l'axe z est confondu avec l'axe de
symétrie.
Figure 5 : Ecoulement à symétrie de révolution

Les lignes de courant sont des surfaces de révolution, dans le plan d'étude désigné par
plan méridien, elles prennent les cotes ψ1, ψ2, ψ3, ψ4. En un point M quelconque du
→
plan, la vitesse méridienne C est tangente à la ligne de courant et a pour composantes :
∂φ
→
Cr =
∂r
C = (54)
∂φ
Ca =
∂z

→
La vitesse C dérive du potentiel φ satisfaisant la propriété d'irrotationnalité.
Comme pour l'écoulement plan, on introduit la fonction de courant ψ satisfaisant

l'équation de continuité (relation 18) :
1 ∂ ∂ Ca
r Cr + = 0 (55)
r ∂r ∂z
Cette équation est satisfaite en posant :
∂φ 1 ∂ψ
Cr = =
∂r r ∂z (56)
∂φ 1 ∂ψ
Cz = = −
∂z r ∂r
Les fonctions φ et ψ ont les mêmes propriétés que celles définies au chapitre II.1
précédent :
- elles forment deux réseaux orthogonaux
- elles satisfont à l'équation de Laplace.
En un point du maillage, le module de la vitesse C est calculable à partir des relations :
∆φ ∆ψ
C = = (57)
∆m 2π r. ∆n
2.3. Méthode simplifiée de tracé des lignes de courant et des équipotentielles
Le tracé des lignes de courant et des équipotentielles s'inspire de la méthode graphique de

PRASIL. Il s'obtient par approximations successives respectant à la fois la propriété
d'orthogonalité et l'équation (52).
Seul un maillage très raffiné permet de connaître en détail le champ des vitesses mais la
méthode proposée ici donne une première approximation relativement précise. On se
place dans le cadre du problème direct et l'on se propose de définir le champ de vitesses
dans le canal formé entre deux aubages de distributeur. L'écoulement à l'infini amont est
supposé uniforme en grandeur et direction.
Dans les conditions d'adaptation, celui-ci aura la même orientation que la direction
d'entrée d'aubage ; dans le cas présent, il s'agira d'un écoulement axial.

En supposant également que la direction de l'écoulement en sortie est la même que celle
de la ligne moyenne, on peut définir sensiblement l'enveloppe de l'écoulement. Elle est
formée, partiellement, dans l'espace interaubages, des lignes intrados et extrados du profil
(figure 6a).
Condition de périodicité :
La condition de périodicité implique l'unicité du vecteur vitesse en tout point de

l'écoulement. Suivant ce principe, à la même abscisse x, les lignes i et e de la figure 6a
doivent présenter une vitesse égale : par exemple, le point 20 de la ligne i et le point 28
de la ligne e. Le respect de cette condition permet de définir la valeur exacte de la vitesse
de sortie par retouches successives de la direction issue du bord de fuite.
A) Maillage simplifié d'un espace inter-aubages
On trace avec un pas sensiblement uniforme l'ensemble des cercles ayant pour
enveloppes les deux frontières de l'écoulement. Ces lignes enveloppes (ψi et ψe) sont
considérées comme des lignes de courant au même titre que le lieu des centres de cercles
que l'on désignera par ligne moyenne (ψ0) ou ligne de partage.
Si l'on joint par une même ligne les trois points qui constituent les points de tangence à
Ψi et Ψe et le centre d'un même cercle, on obtient l'allure des équipotentielles. On
veillera au moment de ce tracé, à l'orthogonalité avec les lignes de courant.
Champ de vitesses :
On se place sur la ligne moyenne et l'on suppose que la vitesse moyenne (ou vitesse
débitante) est obtenue sur cette ligne. On calculera alors successivement :
qv
* la vitesse moyenne : Cj = (60)
h. b j
avec h : épaisseur de la veine fluide
bj : diamètre du cercle correspondant à l'équipotentielle j (de 1 à 28)
* la cote des équipotentielles par : φj + 1 = φj + ∆φj (61)
avec, d'après (52) : ∆φj = C j . ∆sj (62)

∆sj représente la distance curviligne séparant, sur ψo deux équipotentielles j et j+1.
* les vitesses Ci et Ce sur les lignes de courant ψi et ψe : on réutilise pour cela la relation
(52) :
∆φ j
Ci, j = (63)
∆si, j
∆φ j
Ce, j = (64)
∆se, j
∆si,j = distance sur ψi séparant deux équipotentielles.
∆se,j = distance sur ψe séparant les mêmes équipotentielles.
On remarquera que sur l'extrados du profil (ligne ψe) les équipotentielles sont plus
serrées que sur l'intrados. On obtiendra ainsi des vitesses importantes sur ψe au voisinage
des points j = 9, 10, 11, 12 et inversement sur ψi, de faibles vitesses pour les mêmes
valeurs de j.
Nous verrons plus loin qu'à ce champ de vitesses est associé un champ de pression d'où
→
prendra naissance l'effort aérodynamique global F , représentant l'action du fluide sur
l'aubage.
Application :
Les différentes grandeurs sont calculées à partir des relations (60) à (64). Les valeurs de
∆s sont mesurées directement sur la figure 6a à une échelle convenable. La figure 6a
donne également l'évolution de Ce et Ci le long de l'abscisse curviligne de la ligne
moyenne ainsi qu'une comparaison avec les résultats obtenus par Bonataki par une
méthode numérique inspirée de la transformation conforme (référence 201). La figure 6b
montre une comparaison des isovitesses obtenues par les deux méthodes : en haut la
méthode graphique simplifiée, en bas la méthode numérique.
On trouve enfin en figure 6c le maillage simplifié issu de la méthode de l’hodographe
décrite au chapitre G. Ce maillage, réalisé en φ - ψ, donne accès directement au champ de
vitesses.

Figure 6a : Etude de l'espace inter-aubages d'un distributeur de turbine

Figure 6b : Comparaison des isovitesses

Figure 6c : Maillage interaubages issu de la méthode de l’hodographe.
B) Maillage simplifié d'une section méridienne
On se propose de définir les équipotentielles puis la distribution de vitesses dans la

section méridienne correspondant à une pompe centrifuge. S'agissant d'un écoulement
bidimensionnel, la section est dépourvue d'aubages à l'origine des composantes giratoires
Cu normales au plan de l'étude.
Pour cet écoulement à symétrie de révolution, la méthode utilisée est identique à celle
décrite en A mais la relation générique est ici la relation (57) au lieu de la relation (52).
La figure 7a ci-dessous donne le maillage (φ-ψ) de la section étudiée. Obtenu par

itérations successives, il a fait apparaître les équipotentielles et les lignes de courant
satisfaisant à l'équation (57), il peut être aisément automatisé (réf. 202) pour donner la
vitesse de l'écoulement en tout point.
La figure 7b présente le maillage simplifié. Il permet de calculer les vitesses de

l'écoulement à la ceinture et au plafond (frontière extérieure et intérieure de la section
méridienne).

A titre de comparaison, la figure 7c montre l'évolution du rapport de vitesses Cc/Cp (Cc

et Cp étant définies sur la même équipotentielle) en fonction de l'abscisse curviligne de la
ligne moyenne. Ce rapport est défini à partir du maillage simplifié et du maillage
complet.
Sur cette même figure 7c, on a porté quelques points issus de la modélisation de Flügel
décrite dans l'ouvrage de Karassik (réf. 104) et faisant appel aux rayons de courbure des
enveloppes. On peut noter dans tous les cas une forte accélération à la ceinture dans la
zone des points 18-19-20.
Remarque :
Contrairement à l'écoulement plan, le centre des cercles osculateurs d'un écoulement à
symétrie de révolution n'est pas la ligne de courant moyenne. Dans ce cas précis, le
respect de la relation (57) implique le tracé de deux séries de cercles respectant la
propriété de division du débit en deux parties égales. La démarche constructive sera donc
la suivante.
Suite à la première construction, on trace le segment joignant deux points homologues de
la ceinture et du plafond. Deux cercles sont centrés sur ce segment, ils ont la propriété
d'être à la fois tangents entre eux et aux enveloppes de la section d'écoulement. D'autre
part, ils respectent la relation d'égalité de débit :
ri.bi = re.be
La ligne de courant moyenne est le lieu des tangentes communes aux différents couples
de cercles (figure 7d). Cette seconde construction donne de meilleurs résultats dans les
zones à forte courbure.
Figure 7a : Equipotentielles et lignes de courant

Figure 7b : Maillage simplifié
Figure 7c : Comparaison des résultats
Figure 7d : Construction de la ligne moyenne

2.4. - Solutions analytiques de l'équation de Laplace
L'étude des écoulements plans à potentiel des vitesses est facilitée par l'introduction des
fonctions complexes.
On peut définir dans le plan d'écoulement un certain nombre de propriétés permettant

d'accéder à la définition précise des champs de vitesses y compris dans les zones à forte
courbure comme les bords d'attaque.
Dans ce paragraphe, on étudiera en détail quelques uns des écoulements élémentaires les
plus souvent utilisés dans les méthodes de superpositions (méthode des singularités).
Au point M de coordonnées (x,y), représenté sur la figure 8, on fait correspondre une

affixe complexe z telle que :
z = x + iy = r (cos θ + i sin θ) = r exp (i θ) (65)
avec r = x² + y² (module de z) (66)
y
et θ = arc tg (argument de z) (67)
x
Figure 8 : Représentation dans le plan complexe.

→
De même, la vitesse C au point M de composantes (u,v) peut être écrite sous forme
vectorielle complexe :
C = u + i v = C (cos α + i sin α) = C exp (iα) (68)
Nous avons vu également que la fonction de courant ψ et la fonction potentiel φ peuvent

→
être introduites à partir des composantes de la vitesse C , l'une satisfaisant l'équation de
continuité, l'autre l'irrotationnalité de l'écoulement (relation 47) :
∂ψ ∂φ
u = =
→ ∂y ∂x
C = (69)
∂ψ ∂φ
v =− =
∂x ∂y
Dans le cadre de l'étude analytique, on introduit ψ et φ par l'intermédiaire de la fonction
potentiel complexe F(z) :
F( z ) = φ + i ψ (70)
La fonction F(z) pour laquelle φ et ψ satisfont aux relations (69) (conditions de

CAUCHY-RIEMANN) est une fonction analytique dont les parties réelle et imaginaire
satisfont à l'équation de LAPLACE.
Propriétés complémentaires de la fonction potentiel complexe :
dF ∂F ∂φ ∂ψ
= = +i = u − iv (71)
dz ∂x ∂x ∂x
dF ∂F 1  ∂φ ∂ψ 
et = =  +i  = − iv + u (72)
dz i ∂y i  ∂y ∂y 
La fonction potentiel complexe possède une dérivée unique en un point M quelconque.

On retiendra ensuite que la dérivée de F par rapport à z :
dF
= u − iv (73)
dz

→ →
est une vitesse complexe C' appelée vitesse conjuguée de la vitesse C .
A) Ecoulement uniforme : F( z) = Co. z (74)
Le potentiel complexe peut être développé à partir de la définition de l'affixe complexe

(65) :
F(z) = Co.z = Co (x + i y)
et identifié à la définition (70) :
F(z) = φ + i ψ
Identification d'où l'on tire l'équation des équipotentielles et des lignes de courant :
φ = Co.x et ψ = Co.y (75)
Equations d'un réseau maillé que l'on peut tracer (figure 9).
Figure 9 : Ecoulement uniforme parallèle à x

La vitesse locale est calculable à partir de :
∂φ ∂ψ
→ u = = = Co
C = ∂x ∂y
∂φ ∂ψ
v = = − = O
∂y ∂x
La vitesse d'écoulement est donc uniforme en grandeur et direction en tout point du plan.
Dans le cas où l'écoulement n'est pas parallèle à l'axe Ox mais forme un angle α avec
celui-ci, les équations seront les suivantes :
F ( z ) = Co z . exp ( − iα ) = Co r exp ( iθ ) exp ( − iα )

F( z ) = Co r exp i ( θ − α ) = Co r cos ( θ − α ) + i sin ( θ − α )
[
F( z ) = Co ( x cos α + y sin α) + i ( − x sin α + y cos α) ]
soit : φ = Co x cos α + y sin α
ψ = Co − x sin α + y cos α
∂φ ∂ψ
→ u= = = Co cos α
C = ∂ x ∂ y
∂φ ∂ψ
v= =− = Co sin α
∂y ∂x
L'opérateur exp (- iα) correspond à une rotation d'angle α dans le sens trigonométrique.
Remarque :
On accède plus rapidement à ce résultat en utilisant la relation (73) :
F(z) = Co . z . exp (- iα)
= Co . exp ( − iα) = Co [ cos α − i sin α] = u − iv

dF
dz
D’où les composantes de la vitesse :
u = Co cos α
→
C=
v = Co sin α

B) Ecoulement de type source ou puits : F( z) = ± Co ln z (76)
On traite ce cas plus facilement en coordonnées polaires (deuxième partie de la relation

65) :
Cas de la source (signe +)
F(z) = Co ln [r exp (iθ)] = Co ln r + i Co θ (77)
et F(z) = φ + i ψ
donnant après identification :
φ = Co ln r = Co ln x² + y² (78)
y
ψ = Co.θ = Co arc tg (79)
x
Les équipotentielles sont des cercles centrés sur la source (r = cte) et les lignes de courant
des droites passant par l'origine (y/x = cte).
Champ des vitesses :

→
La vitesse locale C est également calculable à partir des relations de base (23) en
coordonnées polaires :
∂φ Co
→
Cr = =
∂r r (80)
C
∂φ
Cu = = O
r∂θ
→
En tout point de l'écoulement, la vitesse C est radiale et varie en raison inverse du rayon.
Débit :
Le débit volume est calculé sur un rayon quelconque :
qv = 2πr h . Cr = 2π Co (81)
avec h : épaisseur de l'écoulement prise égale à l'unité.

qv
soit finalement : Co = (82)
2π
qv
et C = Cr = (83)
2 πr
Dans le cas d'un puits, la vitesse d'écoulement est de sens opposé.
source puits
Figure 10 : Source et puits de même débit.
Si la source ou le puits sont décentrés de l'affixe zo, ils seront définis par la relation :
F(z) = ± Co ln (z-zo) (84)

C) Ecoulement tourbillonnaire F( z) = ± i Co ln z (85)
Traitons le cas positif, en développant la fonction F(z) :
F(z) = i Co ln [r exp (iθ)]
F(z) = i Co ln r - Co θ (86)
en identifiant avec :
F(z) = φ + i ψ
On obtient :
y
φ = − Co . θ = − Co arc tg (87)
x
ψ = Co ln r = Co ln x ² + y ² (88)
Les lignes de courant sont cette fois des cercles concentriques tandis que les
équipotentielles forment un faisceau de droites centré sur l'origine.
Vitesse locale :
∂φ
→
Cr = = O
∂r (89)
C
∂φ Co
Cu = = −
r∂θ r
La vitesse radiale est nulle, la vitesse circonférentielle est inversement proportionnelle à

r. Il s'agit de l'écoulement tourbillonnaire centré en O et désigné aussi par vortex libre.
Le signe + de la relation (85) correspond à une rotation inverse du sens trigonométrique,
le signe - de la même relation correspond au sens trigonométrique.
L'écoulement tourbillonnaire est caractérisé par la circulation du vecteur vitesse (relation

28) :
→ →
Γ = ∫ C . dS (90)
( ψ)
avec dS abscisse curviligne le long du cercle de rayon r formant une ligne de courant :

dS = − r. dθ
Co
C = Cu = −
r
Co
Γ = ∫ . r . dθ = 2π Co (91)
( ψ) r
Quel que soit le cercle retenu, la circulation est une constante. En utilisant (89) et (91) on
peut donner une autre expression à la vitesse C :
Γ
C = Cu = −
2πr
Figure 11 : Ecoulement tourbillonnaire (ou vortex libre)
Remarque :
Compte tenu de la relation (91), l'écoulement tourbillonnaire sera représenté par :
Γ
F( z) = i . ln z (93)
2π

Γ
φ =− θ (94)
2π
avec :
Γ
ψ =+ ln r (95)
2π
Γ
et : C = Cu = − (96)
2πr
Pour un tourbillon, décentré par rapport au centre du repère de l'affixe zo, on écrira :
Γ
F( z) = i ln ( z − z o ) (97)
2π
2.5. Application de la superposition d'écoulements élémentaires
Le principe de superposition des écoulements se traduit par l'addition des fonctions

potentiel complexes et par conséquent par l'addition des fonctions de courant ψ et des
fonctions potentiel φ.
Si l'on a, pour deux écoulements distincts :
F1(z) = φ1 + i ψ1
F2(z) = φ2 + i ψ2
La superposition de ces deux écoulements aura comme potentiel complexe :
F(z) = F1(z) + F2(z) = φ + i ψ (100)
avec φ = φ1 + φ2 (101)
et ψ = ψ1 + ψ2 (102)
Concernant les vitesses, celles-ci s'additionnent vectoriellement.

On démontre ci-dessous que les composantes exprimées dans le même repère, possèdent
la propriété d'additivité :
∂φ ∂ ( φ1 + φ2 ) ∂φ1 ∂φ2
u = = = + = u1 + u 2
∂x ∂x ∂x ∂x
∂φ ∂ ( φ1 + φ2 ) ∂φ1 ∂φ2
v = = = + = v1 + v 2
∂y ∂y ∂y ∂y

Nous proposons dans les pages suivantes quelques unes des superpositions les plus
significatives et des plus utiles en mécanique des fluides.
A) Superposition d'un écoulement uniforme et d'une source
Rappelons les caractéristiques de ces fonctions complexes, avec Co vitesse de

l'écoulement uniforme et qv le débit de la source :
F1 (z) = Co. z avec φ1 = Co x et ψ 1 = Co y
qv qv qv
F2 (z) = ln z avec φ 2 = ln r et ψ 2 = θ
2π 2π 2π
qv qv y
ou φ2 = ln x ² + y ² et ψ2 = arc tg
2π 2π x
La superposition de ces deux écoulements donnera :
qv
F ( z ) = Co z + ln z = φ + i ψ
2π
D'où l'on tire :

qv
φ = Co x + ln x ² + y ²
2π (103)
qv y
et ψ = Co y + arc tg
2π x
qv
La ligne de courant ψ = 0 est obtenue pour y = 0 et yo = − θ correspondant à la
2π Co
qv
courbe en trait fort (figure 12). C'est une courbe à deux asymptotes distantes de et
Co
obtenues pour θ = - π et + π. Le calcul des différentes lignes de courant se fait suivant
l’organigramme de l’annexe R19.
Les fonctions φ et ψ permettent de calculer les composantes de la vitesse locale :
∂ψ ∂φ qv x
u = = = Co +
∂y ∂x 2π x ² + y ²
(104)
∂ψ ∂φ qv y
v =− = =
∂x ∂y 2π x ² + y ²
Le point A est le point d'arrêt de l'écoulement correspondant à l'annulation conjointe de u

et de v. Les coordonnées de ce point sont les suivantes :

qv
xA = −
A 2π Co
yA = O
Figure 12 : Ecoulement autour d'un obstacle (solide de Rankine).
Lorsque x tend vers l'infini, v tend vers zéro et u tend vers Co, vitesse de l'écoulement
uniforme.
Ce principe de superposition permet donc de représenter l'écoulement autour d'un

obstacle (correspondant notamment à ψ = o).
On peut en principe obtenir n'importe quelle forme d'obstacle en disposant dans le plan
un nombre suffisant de sources et de puits.

B) Superposition d'une source et d'un puits de mêmes débits
La source centrée à l'affixe zs (O,a) est caractérisée par les fonctions suivantes (R9) :
qv
F1 ( z) = ln ( z − zS )
2π
qv
φ1 = ln x ² + ( y − a )²
2π
qv y−a
ψ1 = arc tg
2π x
Le puits centré à l'affixe zp (O, - a) est caractérisé par les fonctions :
qv
F2 ( z) = − ln ( z − z p )
2π
qv
φ2 = − ln x ² + ( y + a )²
2π
qv y+a
ψ2 = − arc tg
2π x
L'écoulement correspondant à la superposition des deux écoulements précédents sera lui

caractérisé par (Annexe R10) :
qv z − zS
F( z) = F1 ( z) + F2 ( z) = ln
2π z − zp
qv x ² + ( y − a )²
φ = φ1 + φ2 = ln
2π x ² + ( y + a )²
qv 2 ax
ψ = ψ1 + ψ 2 = arc tg
2π a ² − x² − y²
(105)
Lignes de courant :
La ligne de courant principale ψ = o, obtenue pour x = o, est constituée de l'axe y. Les
autres lignes de courant sont constituées de la famille de cercles centrés sur l'axe x et
passant par les points S et P. Leur rayon Rm est obtenu par la relation géométrique issue
de la figure 13 :
Rm² = a² + (Rm-x1)²
a ² + x12
soit après développement : Rm = (106)
2 x1

Figure 13 : Superposition d'une source et d'un puits.
où x1 représente l'abscisse de l'intersection du cercle avec l'axe Ox. L'abscisse m1 du

centre du cercle est donnée par m1 = - (Rm-x1).
A cette ligne de courant est associée la fonction ψ1 ayant pour valeur :
qv 2 ax1
ψ = ψ1 = arc tg (107)
2π a ² − x12
Les cotes des autres lignes de courant sont obtenues grâce à l'abscisse xi de leur inter-
section avec l'axe Ox.
Equipotentielles :
Les équipotentielles sont constituées de la familles des cercles conjugués centrés sur l'axe
Oy. Si y1 représente l'ordonnée de l'intersection avec l'axe Oy, on observera les relations
géométriques :
Rayon : Rn = n² − a ²
Ordonnée du centre : n = Rn + y1

soit finalement :
a 2 − y12
Rn = (108)
2 y1
et n = Rn + y1 (109)
Comme précédemment, l'équipotentielle est caractérisée par son intersection avec l'axe y
et prend sur ce cercle la valeur :
qv y −a
φ −1 = ln 1 (110)
2π y1 + a
Pour y1 ∈ {o, a}, φ est négatif d'où la présence de l'indice -1.
Représentation générale d'une association source-puits :

Le diagramme suivant donne l'allure générale des lignes de courant et des
équipotentielles d'une association source-puits. Elle est représentée sur la figure 14 pour
une valeur de a égale à 100mm. L’association source-puits sert de base à la méthode de
l’hodographe donnée en détail au chapitre G.
Figure 14 :
Lignes de courant et équipotentielles correspondant
à une association source-puits (valeurs obtenues pour a = 100 mm)

Cas particulier du doublet :
Un doublet est constitué d'une association source-puits centrés au même point D.
Si l'on reprend la relation générale (105) :
qv z − zS
F( z ) = ln (111)
2π z − z P
Le segment PS formant l'angle α avec l'axe Ox et 2a étant sa longueur, on peut écrire :
zS − z P = 2 a exp ( iα )
Le potentiel complexe aura alors pour valeur :
F( z) =
qv LM
ln 1 −
2 a exp ( iα ) OP
2π N z − zP Q
On constitue le doublet en adoptant les limites suivantes :
a  → o , zP  → zD , ln (1 − ε ) → − ε
qv
et − 2a 
→ − m ² (valeur arbitraire notée intensité du doublet)
2π
Le potentiel complexe sera alors donné par :
m² . exp (iα)
F( z) = − (112)
z − zD
Le doublet est représenté sur la figure 15.
Figure 15 : Doublet centré au point D

Cas particulier :
Lorsque D est centré en 0 et α = 0, le doublet est décrit par l'équation :
m²
F( z) = − (113)
z
En posant z = x + i y, l'équation précédente se transforme en :
m² m ² ( x − iy )
F( z ) = − = −
( x + iy ) x² + y²
en posant d'autre part :
F( z) = φ + i ψ
On obtient les expressions de φ et ψ :
m² x
φ = −
x² + y² (114)
m² y
ψ =
x² + y²
C) Superposition d'un écoulement uniforme et d'un doublet
On procède comme pour les cas précédents :
. écoulement uniforme : F1 ( z ) = Co z
m²
. doublet centré en O : F2 ( z ) = Co
z
D'où le potentiel complexe de l'écoulement :
 m² 
F( z) = Co  z +  = φ+iψ (115)
 z 
On accède à φ et ψ en remplaçant z par x + iy et en identifiant les parties réelles et

imaginaires :

 m² x 
φ = Co  x + 
 x² + y² 
(116)
 m² y 
ψ = Co  y − 
 x² + y² 
Définissons la ligne de courant principale (ψ = o), elle est constituée par :
y = o ( axe Ox )
x ² + y ² = m² ( cercle de centre O et de rayon m)
Le champ de vitesses est donné par les relations classiques :
∂φ ∂ψ  ( x² + y² )m² − 2x² m² 
u = = = Co 1 + 
∂x ∂y  ( x² + y² )² 
(117)
∂φ ∂ψ 2x y m ²
v = = − = − Co
∂y ∂x ( x² + y² )²
Les autres lignes de courant seront obtenues en suivant l’organigramme de l’annexe R20.
Champ de vitesses à la surface du cylindre :

On remplace (x²+y²) par m² (équation du cercle) dans les équations définissant u et v, on
obtient les valeurs particulières uo et vo :
 x² 
u o = 2 Co 1 −
 m ² 
xy
v o = − 2 Co
m²
y x
En remarquant que sur le cercle sin θ = et cos θ = , on aura :
m m
u o = 2 Co sin ² θ
v o = − 2 Co sin θ cos θ
→
D'où l'on peut déduire le module de la vitesse C :
C = 2 Co sin θ (118)

Les points A (θ = π) et A' (θ = o) sont les points d'arrêt de l'écoulement.
La superposition d'un écoulement uniforme et d'un doublet se présente comme

l'écoulement irrotationnel autour d'un cylindre de révolution de rayon m (figure 16) :
Figure 16 : Ecoulement autour d'un cylindre de rayon m
D) - Superposition d'une source et d'un vortex centrés en O :
On étudie ici une source de débit qv centrée en O associée à un vortex centré également
en 0 et d'intensité Γ.
Le potentiel complexe correspondant est donné par :
qv Γ
F( z ) = ln z + i ln z (119)
2π 2π
En retenant les coordonnées polaires et en remplaçant z par r exp (iθ), on obtient :
qv Γ
F( z ) = (ln r + i θ) + i (ln r + i θ)
2π 2π
donnant :

 qv Γ
 φ = 2π ln r − 2π θ
 (120)
 ψ = qv θ + Γ ln r
 2π 2π
ainsi que le champ de vitesses :
∂φ qv
Cr = =
∂r 2 πr (121)
∂φ Γ
Cu = =−
r∂θ 2 πr
L'angle β formé entre la direction radiale et la tangente à la ligne de courant est donnée
par (figure 17) :
Cu Γ
tg β = =− (122)
Cr qv
L'angle β est donc constant quel que soit le rayon proposé. Les lignes de courant sont des
spirales logarithmiques d'équation :
 2 π ψ  θ 
r = exp  . exp − 
 Γ   tg β 
L'écoulement est axisymétrique puisque l'on observe simultanément :
∂ Cr ∂ Cu
= = O
∂θ ∂θ
Figure 17 : Superposition d'une source et d'un vortex

Il s'agit de l'écoulement relatif simplifié produit par une roue centrifuge équipée
d'aubages en forme de spirales logarithmiques. On démontre que l'écoulement n'est
véritablement axisymétrique que dans le cas d'un nombre infini d'aubages (hypothèse
d'Euler).

Ch.B – Dynamique des fluides incompressibles
Chapitre
- B -B
DYNAMIQUE DES FLUIDES
DYNAMIQUE DES FLUIDES
INCOMPRESSIBLES
INCOMPRESSIBLES

I - EQUATIONS DE LA DYNAMIQUE DES FLUIDES PARFAITS p. 4

1.1. Intégration le long d’une ligne de courant
1.2. Théorème des quantités de mouvement
II - PROPRIETES DU FLUIDE REEL p. 11

2.1. Viscosité
2.2. Description du frottement visqueux
2.3. Généralisation
2.4. Tenseur des contraintes
2.5. Liaison entre contraintes et taux de déformation
III - EQUATION DE NAVIER-STOKES p. 23
IV - CARACTERISTIQUES DE L’ECOULEMENT VISQUEUX p. 29

4.1. Régimes d’écoulement
4.2. Propriétés des écoulements laminaires
4.3. Propriétés des écoulements turbulents
4.4. Modélisation de la turbulence
V - COUCHE LIMITE BIDIMENSIONNELLE p. 43

5.1. Equations régissant la couche limite laminaire
5.2. Evolution de la couche limite laminaire à la surface d’une paroi de forme quelconque
5.3. Autres paramètres caractéristiques de la couche limite
5.4. Couche limite turbulente
5.5. Couche limite turbulente sur une surface quelconque
5.6. Ordre de grandeur des caractéristiques de la couche limite sur une plaque plane
VI - ECOULEMENTS EN CONDUITE p. 76
6.1. Régime laminaire
6.2. Régime turbulent lisse
6.3. Influence de la rugosité

DYNAMIQUE DES FLUIDES INCOMPRESSIBLES
Au même titre que le chapitre A, ce second chapitre constitue la partie classique et

académique de la mécanique des fluides. En général, l’étude des écoulements visqueux
nécessite une étude préalable en fluide parfait qui sert de référence aux calculs mettant en
évidence le profil de la couche limite, son éventuel décollement et les pertes par frottement
associées.
L’équation de la dynamique des fluides réels est donnée dans son expression la plus générale
sous la forme de l’équation de Navier-Stokes.
Elle est obtenue par le développement des principes fondamentaux de la mécanique :

r
( F = mΓ)
r
l’équation de Newton ou le théorème des quantités de mouvement. Son
établissement dans sa forme actuelle est basé sur les travaux de Stokes, Saint-Venant, Navier
et Poisson.
Exceptés quelques cas très limités où il existe une solution exacte, l’équation de Navier-
Stokes ne possède pas de solution analytique générale et constitue ainsi le champ
d’application privilégié des numériciens.
Les diverses formes de l’équation de l’énergie seront examinées dans le cadre de la

dynamique des écoulements compressibles. Dans une large mesure, les écoulements de gaz
peuvent être considérés comme incompressibles tant que le nombre de Mach reste inférieur à
0,4.

I - EQUATIONS DE LA DYNAMIQUE DES FLUIDES PARFAITS
L’équation dynamique est l’équation de Newton traduisant, pour un élément fluide, l’égalité
des efforts extérieurs et des efforts d’inertie :
→ →
F = Γ . dm (1)
avec dm : masse de l’élément fluide dm = ρ . dV

→
Γ : accélération dans le repère fixe (galiléen)
→
F : résultante des efforts extérieurs appliqués
Les efforts agissant sur la particule sont classés en deux groupes principaux :
• Forces agissant à la surface de l’élément fluide :

- efforts de pression agissant perpendiculairement aux parois : Fp ;
- efforts de frottement agissant tangentiellement et se manifestant en présence d’un fluide
visqueux : Ff ;
- efforts matériels tels que l’effort d’un aubage sur l’élément de fluide : Fa.
• Forces de volume : Fv
Efforts à distance tels que la pesanteur, les effets magnétiques ou tout autre effort extérieur
agissant sans contact. Plaçons nous en coordonnées cartésiennes dans l’hypothèse d’un fluide
parfait (Ff = 0). Considérons un volume de fluide élémentaire éloigné d’une paroi (Fa = 0) sur
→
lequel agit un effort de volume FV de composantes (Fx, Fy, Fz).
→
Examinons les efforts en présence suivant la direction Ox (Figure 1) :
Cliquez ici pour taper du texte.
Figure 1 : Disposition des efforts de pression suivant Ox

 ∂p  ∂p
Forces de pression : p dy dz −  p + dx dy dz = − dx dy dz
 ∂x  ∂x
Forces de volume : Fx
Masse de l’élément : ρ dV = ρ dx dy dz
Accélération : Γx
Le bilan établit l’équation d’équilibre issue de l’équation (1) :
∂p
− dx dy dz + Fx = ρΓx dx dy dz (2)
∂x
→
En extrapolant aux autres composantes et en retenant f comme effort par unité de masse :
→
→ FV
f = (3)
ρ dx dy dz
On obtient l’équation vectorielle :
1 → → →
− grad p + f = Γ (4)
ρ
→
L’accélération Γ de la particule fluide est exprimée par ses composantes cartésiennes au
chapitre A (relation 10). En reportant celles-ci dans l’équation (4), on obtient le système
d’équations différentielles suivant, désigné par équations d’Euler :
1 ∂p ∂u ∂u ∂u ∂u
− + fx = u +v +w +
ρ ∂x ∂x ∂y ∂z ∂t
1 ∂p ∂v ∂v ∂v ∂v
− + fy = u +v +w + (5)
ρ ∂y ∂x ∂y ∂z ∂t
1 ∂p ∂w ∂w ∂w ∂w
− + fz = u +v +w +
ρ ∂z ∂x ∂y ∂z ∂t
qui s’écrivent aussi en coordonnées cylindriques :

1 ∂p ∂Cr ∂Cr ∂Cr Cu ² ∂Cr

− + fr = Cr + Cu + Ca − +
ρ ∂r ∂r r∂θ ∂z r ∂t
1 ∂p ∂Cu ∂Cu ∂Cu CrCu ∂Cu
− + fθ = Cr + Cu + Ca + + (6)
ρ r∂θ ∂r r∂θ ∂z r ∂t
1 ∂p ∂Ca ∂Ca ∂Ca ∂Ca
− + fz = Cr + Cu + Ca +
ρ ∂z ∂r r∂θ ∂z ∂t
En général, les forces de volume dérivent d’un potentiel et, plus particulièrement, de celui de
la pesanteur. On notera dans ce cas :
→ →
f V = − grad (gh ) (7)
avec g : accélération de la pesanteur

h : cote ou altitude du point considéré
Ainsi, l’expression la plus générale de l’équation dynamique s’écrira en fluide parfait :
1 → → →
− grad p − grad ( gh ) = Γ (8)
ρ
En remplaçant l’accélération de la particule par sa valeur issue de la relation vectorielle (12)

du chapitre A, l’équation dynamique s’écrira sous la forme de Lamb :
→
1 → → ∂C → C² → → →
− grad p − grad ( gh ) = + grad − C ∧ rot C
ρ ∂t 2
Pour un fluide baratrope pour lequel la masse volumique ne dépend que de la pression
(ρ = ρ(p)), le premier terme se transformera en :
1 → → dp
−
ρ
grad p = − grad
ρ ∫
→
Enfin, le groupement des termes grad donnera :

→
→  dp C²  ∂ C → → →
− grad 
∫
 ρ
+ gh +  =
2  ∂t
− C ∧ rot C (9)
1.1. Intégration le long d’une ligne de courant
La relation générale (9) peut être intégrée le long d’une ligne de courant d’abscisse curviligne
ds.
Figure 2 : Ligne de courant, ligne tourbillonnaire, surface de Helmotz
On démontre que le dernier terme de la relation (9) peut s’exprimer sous la forme :
→ → → →
C ∧ rot C = grad H
où H est une constante sur la surface de Helmotz contenant à la fois la ligne de courant
et la ligne tourbillonnaire portant le vecteur tourbillon :
→ 1 → →
Ω= rot C
2
L’intégration donne le résultat suivant :
∫ dpρ + gh + C2² + ∫ ∂∂Ct ds = H ( t ) (10)

où H(t) est une constante variant d’une ligne de courant à l’autre ainsi qu’en fonction
du temps. L’écoulement rotationnel est donc le siège d’un échange d’énergie avec le
milieu environnant. Cependant, entre deux points 1 et 2 de la même ligne de courant,
on pourra écrire à chaque instant :
  dp ∂C  C12   dp ∂C  C22
 
 ρ
∫ + ds
∂t 
 + gh1 + = 
2  ρ
 + ds∫

∂t 
 + gh2 +
2
(11)
1 2
Cas de l’écoulement irrotationnel :
Dans cette hypothèse, le dernier terme du second membre de (9) s’annule et la vitesse dérive
d’un potentiel φ :
→ →
rot C = 0
→ →
C = grad φ
→
∂C ∂  →  →  ∂φ 
= grad φ = grad  
∂t ∂t  
  ∂t 
L’équation dynamique prend la forme :
→  dp C² ∂φ 
grad 
 ρ
∫
+ gh +
2
+ =0
∂t 
soit, après intégration :
∫ dpρ + gh + C2² + ∂φ∂t = Cte(t) (12)
C’est la forme la plus générale du théorème de Bernoulli où le terme Cte(t) est une constante
d’intégration dépendant du temps mais identique pour tous les points de l’écoulement.
Au même instant, en deux points quelconques 1 et 2, situés ou non sur la même ligne de
courant, on observera l’égalité pratique :
 dp ∂φ  C1 ²  dp ∂φ 

 ρ
∫
+  + gh1 +
∂t 1 2
= 
 ρ
+ ∫
 + gh 2 +
∂t  2
C2 ²
2
(13)

Autres cas particuliers :
• Ecoulement permanent : On observe simultanément que ∂φ/∂t = 0 et que la constante

d’intégration est indépendante du temps.
• Fluide incompressible : Dans le cas où ρ = cte, le terme de pression se réduit à

l’expression :
∫ dρp = 1ρ ∫ d p = pρ
Ainsi, la forme la plus simplifiée de l’équation de Bernoulli sera donnée par :
p C²
+ gh + = cte (14)
ρ 2
Elle correspondra, sous cette forme, à l’écoulement permanent et irrotationnel d’un fluide
parfait incompressible.
1.2. Théorème des quantités de mouvement

Il s’agit d’une autre expression de l’équation dynamique que nous établirons sous forme
intégrale dans le cas particulier d’un tube de courant en régime stationnaire.
L’énoncé de ce théorème est le suivant :
La quantité de mouvement sortant d’un volume fluide est égale à la somme des forces
extérieures appliquées à ce volume.
S’agissant d’un tube de courant (Figure 3), aucun fluide ne traverse l’enveloppe extérieure et
seules les sections terminales du tube débitent.
La quantité de mouvement est un effort obtenu par le produit du débit massique Q par la
vitesse locale.
Figure 3 : Equation des quantités de mouvement appliquée à un tube de courant

En considérant les vitesses uniformes à l’entrée et la sortie du tube de courant, la quantité de

mouvement sortante a pour expression :
 → → → →
Q  C2 − C1  = ρ2 qv 2 C2 − ρ1 qv1 C1 (15)
 
→ →
Les modules de C1 et C 2 sont calculables à partir du débit volume et des sections :
qv1
C1 =
S1
qv2
C2 =
S2
Au titre des forces extérieures, on retiendra comme précédemment les forces de volume Fv (la
pesanteur essentiellement), les forces de contact Fa correspondant à l’action des parois ou du
fluide environnant, les forces de frottement Ff agissant tangentiellement à la surface du tube et
enfin les efforts de pression agissant sur les surfaces S1 et S2.
→
En regroupant tous les efforts sous la forme résultante F hormis les efforts de pression aux
extrémités, le théorème des quantités de mouvement s’écrira :
→ → → → →
Q  C 2 − C1  = F + p 1S1 − p 2S 2 (16)
 
L’avantage de cette forme de l’équation dynamique est qu’elle permet d’évaluer l’effort
extérieur global sans connaître par le détail la distribution de pression à la surface de contact
avec le milieu extérieur : fluide environnant ou paroi matérielle (aubages, conduite, etc...).

II - PROPRIETES DU FLUIDE REEL
Le fluide réel se caractérise, en opposition au fluide parfait, par deux propriétés importantes :
- La viscosité, qui est une caractéristique physique du fluide se manifestant par une résistance
de celui-ci aux déformations et plus particulièrement aux vitesses de déformation.
Elle est due à la combinaison des efforts de cohésion et d’agitation moléculaire s’opposant au
déplacement relatif des couches liquides les unes par rapport aux autres.
La viscosité dépend de la nature du fluide et varie considérablement d’un fluide à un autre.

Pour un même fluide, elle dépend de la pression et de la température, pour les liquides
incompressibles, elle est pratiquement invariable avec la pression.
- L’adhérence du fluide aux parois solides : contrairement au fluide parfait pour lequel la
vitesse d’écoulement possède une valeur non nulle à la paroi, le fluide réel adhère
parfaitement à celle-ci et sa vitesse est nulle à cet endroit.
Dans cette zone de fort gradient de vitesse désignée par couche limite, les effets du frottement
interne sont importants (Figure 4).
Figure 4 : Ecoulement au voisinage de la paroi
Nous verrons que l’épaisseur δ de la couche limite dépend à la fois de la viscosité du fluide et
du régime d’écoulement.

A l’extérieur de la couche limite, la cinématique du fluide est pratiquement identique à celle

décrite en fluide parfait. Ainsi, dans certains modèles d’écoulement, les conditions de
glissement ne s’écriront pas à la paroi elle-même mais à la frontière de la couche limite.
2.1. Viscosité
L’expérience de Couette est souvent utilisée pour mettre en évidence l’existence de la

viscosité.
On dispose d’une plaque de surface S infiniment mince flottant à la surface libre d’un liquide
au repos. La hauteur de liquide est égale à h. Pour atteindre une vitesse de déplacement égale
à U, il faut exercer sur la plaque, en régime permanent, un effort F (Figure 5).
La contrainte de cisaillement à la surface de la plaque est égale à :
F
τ= (17)
S
Figure 5 : Expérience de Couette
Pour un fluide de viscosité donnée, on observe que τ est proportionnelle à la vitesse et

inversement proportionnelle à la hauteur de fluide :
U
τ = µ. (18)
h
µ : coefficient de proportionnalité qualifié de coefficient de viscosité ou encore de viscosité

dynamique

Pour les conditions d’adhérence déjà présentées, la vitesse du fluide est nulle au fond et égale
à U à la surface de la plaque. La vitesse en tout point du fluide compris entre le fond et la
plaque étant supposée linéaire, celle-ci est donnée par :
U
u= y (19)
h
Le gradient transversal s’écrit donc d’après l’équation (19) :
∂u U
= (20)
∂y h
L’hypothèse de Newton consiste à généraliser cette propriété en écrivant :
∂u
τ =µ. (21)
∂y
La viscosité dynamique µ s’exprime dans le système SI en N.s/m² ou en kg/m.s ou encore

Poiseuille (Pl). On la trouve très souvent exprimée à partir des unités du système CGS en
g/cm.s ou Poise (Po).
D’autres lois de comportement peuvent être observées, dépendant de la nature du fluide :
∂u
τ = τ0 + µ (22)
∂y
La présence de la contrainte critique ou contrainte de Bingham τ0 montre la nécessité de

dépasser un certain seuil de contrainte pour engendrer l’écoulement du fluide. Ces fluides,
qualifiés généralement de non-newtoniens, se rencontrent industriellement aussi couramment
que les fluides newtoniens.
A l’inverse des fluides newtoniens qualifiés de liquides purs, les fluides non newtoniens sont
généralement constitués de mélanges liquide/solide : sang, pâte à papier, lait, etc... ou encore
d’émulsions diverses liquide/liquide ou gaz/liquide.
2.2. Description du frottement visqueux

Examinons dans le plan (x, y) la déformation de la face du parallélépipède rectangle
élémentaire représentant une particule fluide en mouvement. Les gradients de vitesse seront
supposés positifs (Figure 6).

Figure 6 : Déformations d’une face de la particule
A l’instant t = t0, la face considérée est supposée rectangulaire de côtés ∆x, ∆y. A l’instant
t = t0 + ∆t, le point 0 s’est déplacé jusqu’au point A de composantes :
u∆t
A = v∆t
u et v étant les composantes du vecteur vitesse au point 0.
Pendant le même instant, le point O’ s’est déplacé en B. Le théorème des accroissements finis
au premier ordre indique pour le point O’ une vitesse ayant pour composantes :
∂u
u+ ∆x
∂x
C0'
∂v
v+ ∆x
∂x
D’où les coordonnées du point B :
 ∂u 
∆x +  u + ∆x ∆t
 ∂x 
B
 ∂v 
v + ∆x ∆t
 ∂x 

Le segment OO’ a subi à la fois une translation, un allongement et une rotation qui sont tous
trois calculables. S’agissant de faibles déformations, l’angle ∆α est assimilable à sa tangente :
y − yA
∆α = B (23)
xB − xA
∆α est calculable à partir des coordonnées de A et B :
 ∂v 
v + ∆x ∆t − v∆t
 ∂x 
∆α =
 ∂u 
∆x +  u + ∆x ∆t − u∆t
 ∂x 
qui, après simplifications, s’écrit :
∂v
∆α = ∆t (24)
∂x
De la même façon, on peut calculer les coordonnées du point C transformé du point O’’ :
 ∂u 
u + ∆y ∆t
 ∂y 
C
 ∂v 
∆y +  v + ∆y ∆t
 ∂y 
et la rotation ∆β :
 ∂u 
u + ∆y ∆t − u∆t
xC − xA  ∂y 
∆β = =
yC − yA  ∂v 
∆y + v + ∆y ∆t − v∆t
 ∂y 
∂u
∆β = . ∆t (25)
∂y
La rotation de la diagonale AD notée ∆ε sera déduite des relations (24) et (25) :
∆α − ∆β
∆ε =
2

1  ∂v ∂u 
∆ε =  ∆t − ∆t 
2  ∂x ∂y 
D’où l’on peut calculer la vitesse angulaire de la particule :
∆ε 1  ∂v ∂u 
Ωz = lim =  −  (26)
∆t → 0 ∆t 2  ∂x ∂y 
Le dernier membre de cette équation représente la composante sur z du vecteur tourbillon

r 1 → r
Ω = rot C .
2
Dans le cas d’un écoulement irrotationnel, la particule se déformera sans changement de

direction de sa diagonale et subira donc une déformation sans rotation (∆ε = 0).
Taux de déformation :
Suivant la direction Ox, l’élongation du segment AB s’écrit :
∆X = (xB - xA) - ∆x
  ∂u  
∆X =  ∆x +  u + ∆x ∆t − u ∆t  − ∆x
  ∂x  
∂u
∆X = ∆x ∆t
∂x
La vitesse d’élongation sera quant à elle définie par :
∆X ∂u
V= lim = . ∆x
∆t → 0 ∆t ∂x
et le taux de déformation par :
V ∂u
=
∆ x ∂x
∂u
Suivant la direction Oy, le taux de déformation du même segment est égal à .
∂y

Concernant le segment AC, les taux de déformation suivant Ox et Oy sont respectivement :
∂v ∂v
et
∂x ∂y
On peut ainsi construire le tenseur des taux de déformation dans le plan (x, y) :
∂u ∂u
∂x ∂y
T= (27)
∂v ∂v
∂x ∂y
que l’on peut décomposer en deux tenseurs : un tenseur d’élongation symétrique et un tenseur
de rotation antisymétrique :
T=D+Ω (28)
∂u 1  ∂u ∂v  1  ∂u ∂v 
 +  0  − 
∂x 2  ∂y ∂x  2  ∂y ∂x 
T= +
1  ∂v ∂u  ∂v 1  ∂u ∂v 
 +  −  −  0
2  ∂x ∂y  ∂y 2  ∂y ∂x 
Concernant le tenseur D , les termes diagonaux représentent les taux d’élongation des côtés de
la particule dans les directions x et y tandis que les autres termes non diagonaux représentent
la déformation par cisaillement. Les deux côtés de la particule aboutissant au point A
subissent une rotation respective ∆α et ∆β de sens inverse. Le taux de cisaillement dans le
plan (x, y) est donc égal à :
 ∆α + ∆β 
lim  
∆t → 0  ∆t 
soit d’après les relations (24) et (25) :
 ∂v ∂u 
 +  (29)
 ∂x ∂y 
Cette valeur est proportionnelle à l’allongement et au raccourcissement des deux diagonales

du parallélogramme.

2.3. Généralisation
Si l’on considère une particule de forme parallélépipédique ayant pour dimensions respectives
∆x, ∆y, ∆z, le principe développé précédemment peut être aisément généralisé.
Dans un champ de vitesses quelconque, le tenseur des taux de déformation s’écrira comme la
superposition d’un tenseur d’élongation symétrique et d’un tenseur de rotation antisymétrique.
∂u ∂u ∂u
∂x ∂y ∂z
∂v ∂v ∂v
T= ∂x ∂y ∂z =D+Ω (30)
∂w ∂w ∂w
∂x ∂y ∂z
∂u 1  ∂u ∂ v  1  ∂w ∂u 
 +   + 
∂x 2  ∂y ∂x  2  ∂x ∂ z 
1  ∂u ∂ v  ∂v 1  ∂ v ∂w 
T=  +   +  (31)
2  ∂y ∂x  ∂y 2  ∂z ∂y 
1  ∂w ∂u  1  ∂v ∂ w  ∂w
 +   + 
2  ∂x ∂z  2  ∂z ∂ y  ∂z
1  ∂u ∂v  1  ∂w ∂u 
0  −  −  − 
2  ∂y ∂x  2  ∂x ∂z 
1  ∂u ∂v  1  ∂v ∂w 
+ −  −  0  −  (32)
2  ∂y ∂x  2  ∂z ∂y 
1  ∂w ∂u  1  ∂v ∂w 
 −  −  −  0
2  ∂x ∂z  2  ∂z ∂y 
Remarque :
A l’instant t = t0 + ∆t les dimensions principales de la particule deviennent (Figure 6) :
 ∂u 
AB ∆x 1+ ∆t 
 ∂x 

 ∂v 
AC ∆y 1+ ∆t 
 ∂y 
 ∂w 
sur Oz ∆z 1 + ∆t 
 ∂z 
Le volume de celle-ci devient au premier ordre :
  ∂u ∂v ∂w  
∆x ∆y ∆z 1 +  + +  ∆t 
  ∂x ∂y ∂z  
soit encore :
r
∆x ∆y ∆z [1+ div C . ∆t] (33)
r
Pour un fluide incompressible (ρ = cte), le terme div C est nul, satisfaisant ainsi l’équation de
continuité (relation A17). Au premier ordre, en fluide incompressible, la particule tourne et se
déforme sans changer de volume.
2.4. Tenseur des contraintes

Le tenseur C des contraintes appliquées à la surface d’une particule de forme parallélépipé-
dique de dimensions élémentaires ∆x, ∆y, ∆z est composé des contraintes normales et des
contraintes tangentielles. Examinons sur l’une des faces supposée contenue dans le plan yOz,
les contraintes en présence (Figure 7).
→
L’effort F exercé par le fluide environnant donne lieu à deux contraintes :
σxx contrainte normale

τxy  composantes de la contrainte de cisaillement τx
τxz 

Le double indice associé à ces contraintes doit être interprété de la manière suivante :
Le premier (x) correspond à la face orientée par l’axe correspondant, le second (x, y ou z)
correspond à la direction de projection.
Pour une face orientée de manière quelconque dans l’espace, le tenseur des contraintes
s’écrira :
σxx τxy τxz

C = τyx σyy τyz (34)
τzx τzy σzz
On démontre que les termes symétriques sont égaux entre eux :
τxy = τyx
τxz = τzx
τyz = τzy
Pour un fluide au repos ou un fluide de viscosité nulle, les contraintes tangentielles sont
nulles, les contraintes normales sont alors égales et opposées à la pression hydrostatique p.
On utilisera cette propriété pour écrire le tenseur C sous la forme de la superposition de deux
états :
−p 0 0 σxx + p τxy τxz

C= 0 −p 0 + τyx σyy + p τyz (35)
0 0 −p τzx τzy σzz + p
soit C = − pI + τ (36)
1 0 0
avec I= 0 1 0 tenseur unitaire
0 0 1
σxx + p τxy τxz

et τ = τyx σyy + p τyz
τzx τzy σzz + p
tenseur déviateur représentant uniquement les effets du frottement.

2.5. Liaison entre contraintes et taux de déformation
Trois principes fondamentaux sont utilisés dans cette partie :La rotation sans déformation de
la particule fluide n’engendre aucune contrainte de cisaillement.
• Les contraintes sont proportionnelles aux taux de déformation (généralisation de l’équation
de Newton).
• Le fluide est isotrope.
On observera donc entre les tenseurs la relation générale :
r
τ = ε (div C ) . I + 2µ D (37)
avec ε : viscosité de dilatation

µ : viscosité dynamique
Relation permettant d’accéder aux différentes composantes des contraintes :
→ ∂u
σxx = ε div C + 2 µ −p
∂x
Contraintes normales σyy = ε div C + 2 µ ∂v − p

→
(38)
∂y
→ ∂w
σzz = ε div C + 2 µ −p
∂z
 ∂u ∂v 
τxy = τyx = µ  + 
 ∂y ∂x 
 ∂w ∂u 
Contraintes tangentielles τxz = τzx = µ  +  (39)
 ∂x ∂z 
 ∂v ∂w 
τyz = τyz = µ  + 
 ∂z ∂y 
Pour un fluide compressible ou non, l’addition des trois contraintes normales conduit à :
→
σxx + σyy + σzz = (3ε + 2µ) div C − 3 p
permettant d’exprimer la pression hydrostatique locale :
→
1  2 
p = − (σxx + σyy + σzz) +  ε + µ  div C (40)
3  3 

→
En fluide incompressible, l’équation de continuité div C = 0 montre que la pression est égale
et opposée à la moyenne des contraintes normales.
En fluide compressible, on admet qu’il en va de même, ce qui implique que le terme

2
ε + µ = 0.
3
C’est une hypothèse que l’on adopte couramment sous forme de l’équation de Stokes :
2
ε=− µ (41)
3
Dans l’écoulement plan pour lequel v = w = 0, les relations (39) se ramènent à :
∂u
τxy = µ .
∂y
τxz = 0 (42)
τyz = 0
correspondant au cas de l’expérience de Couette décrite au § 2.1.

III - EQUATION DE NAVIER-STOKES
C’est l’équation de la dynamique des fluides établie en fluide réel pour lequel la viscosité
engendre des contraintes de cisaillement.
Reprenons les concepts et les efforts développés au § I concernant une particule fluide et
→
examinons les efforts en présence suivant la seule direction Ox . A l’inverse de la Figure 1, on
observe dans le cadre de l’écoulement d’un fluide réel non seulement les contraintes normales
mais également les contraintes tangentielles (Figure 8).
→
Figure 8 : Actions à la surface de l’élément fluide suivant Ox
 ∂σxx 
Forces de surface :  dx dydz
 ∂x 
 ∂τyx 
+ dy dxdz
 ∂y 
 ∂τzx 
+ dz dxdy
 ∂z 
Forces de volume : Fx
Accélération : Γx

L’équation de la dynamique (1) appliquée à l’élément fluide donne l’équation suivante :
 ∂σxx ∂τyx ∂τzx 

Γxdm = Fx +  + +  dxdydz
 ∂x ∂y ∂z 
se réduisant après simplification à :
Fx ∂σxx ∂τyx ∂τzx

ρ Γx = + + + (43)
dx dy dz ∂x ∂y ∂z
En remplaçant les contraintes σxx, τyx, τzx par leurs valeurs issues de (38) et (39), on aura :
Fx ∂  → ∂u  ∂   ∂u ∂v  ∂   ∂w ∂u 
ρ Γx = +  ε div C + 2µ − p + µ  +  + µ  + 
dx dy dz ∂x  ∂x  ∂y   ∂y ∂x  ∂z   ∂x ∂z 
En supposant µ constant, elle s’écrit également :
Fx ∂p ∂  →  ∂u ∂v ∂w   ∂² u ∂² u ∂² u 
ρ Γx = − + ε div C + µ  + +  + µ  + +  (44)
dx dy dz ∂x ∂x   ∂x ∂y ∂z   ∂ x ² ∂ y² ∂ z ² 
r
soit encore en identifiant div C et le laplacien de u.
Fx ∂p ∂  →

ρ Γx = − + ( ε + µ ) div C  + µ ∇² u (45)
dx dy dz ∂x ∂x  
r
En retenant le cas des fluides incompressibles pour lesquels le terme div C est nul (équation
de continuité), en posant d’autre part :
→
f = effort à distance par unité de masse
→
→ F
f =
ρ dx dy dz
µ
et ν = viscosité cinématique,
ρ
l’équation de Navier-Stokes se réduit à :
→ → 1 → →
Γ= f − grad p + ν ∇ ² C (46)
ρ

En coordonnées cartésiennes (46) s’écrit comme (5) mais figurent en plus les termes
visqueux :
∂u ∂u ∂u ∂u 1 ∂p  ∂² u ∂² u ∂² u 
+u +v +w =− + fx + ν  + + 
∂t ∂x ∂y ∂z ρ ∂x  ∂x ² ∂y ² ∂z ² 
∂v ∂v ∂v ∂v 1 ∂p  ∂² v ∂² v ∂² v 
+u +v +w =− + fy + ν  + +  (47)
∂t ∂x ∂y ∂z ρ ∂y  ∂x ² ∂y ² ∂z ² 
∂w ∂w ∂w ∂w 1 ∂p  ∂² w ∂² w ∂² w 
+u +v +w =− + fz + ν  + +
∂t ∂x ∂y ∂z ρ ∂z  ∂x ² ∂y ² ∂z ² 
et, en coordonnées cylindriques, la forme plus complexe (48) :
∂Cr ∂Cr Cu ∂Cr ∂Cr Cu ² 1 ∂p

+ Cr + + Ca − =− + f r + ...
∂t ∂r r ∂θ ∂z r ρ ∂r
∂Cu ∂Cu Cu ∂Cu ∂Cu CrCu 1 ∂p
+ Cr + + Ca + =− + f + ...
∂t ∂r r ∂θ ∂z r ρ r∂θ θ
∂Ca ∂Ca Cu ∂Ca ∂Ca 1 ∂p
+ Cr + + Ca =− + f z + ...
∂t ∂r r ∂θ ∂z ρ ∂z
 ∂²Cr 1 ∂Cr 1 ∂²Cr ∂²Cr Cr 2 ∂Cu 

... + ν  + + + − −
 ∂r ² r ∂r r ² ∂θ² ∂z ² r ² r ² ∂θ 
 ∂²Cu 1 ∂Cu 1 ∂²Cu ∂²Cu Cu 2 ∂Cr 
... + ν  + + + − + (48)
 ∂r ² r ∂r r ² ∂θ² ∂z ² r ² r ² ∂θ 
 ∂²Ca 1 ∂Ca 1 ∂²Ca ∂²Ca 
... + ν  + + +
 ∂r ² r ∂r r ² ∂θ² ∂z ² 
Lorsque la viscosité du fluide est nulle, on retrouve naturellement l’équation (4) d’Euler.
→
L’accélération Γ de la particule fluide est exprimée par ses composantes cartésiennes et
cylindriques au chapitre A. Les propriétés du laplacien sont rappelées en annexe R.6.
La viscosité cinématique ν s’exprime dans le système d’unité SI en m²/s et dans le système

CGS en cm²/s ou Stoke (St). L’unité pratique est le plus souvent le centi-Stoke (cSt) égal à
1 . 10-6 m²/s.
On retiendra quelques valeurs importantes :
• pour l’eau à 20°C sous 1 atm (760 mm Hg) : ν = 1 cSt = 1 . 10-6 m²/s
• pour l’air dans les mêmes conditions : ν = 15,1 cSt = 15,1 . 10-6 m²/s

Cas du fluide compressible :
Si l’on reprend l’expression (45) et l’hypothèse formulée précédemment indiquant que la

pression hydrostatique est la moyenne des contraintes normales (41) et conduisant à
ε = - 2/3 µ, l’équation de Navier-Stokes s’écrit :
→ → 1 →  1 →
 →
Γ= f − grad p − µ div C + ν∇² C (49)
ρ  3 
Cas du fluide non-newtonien :
Pour les fluides non-newtoniens, la loi de comportement est donnée par exemple par la forme
de Bingham (22) :
∂u
τ xy = τ 0 + µ
∂y
où τ0 est une contrainte de cisaillement critique dépendant de plusieurs facteurs.
• d’une part, de sa composition, le fluide non-newtonien étant le plus souvent constitué d’un
mélange diphasique ou de fluides de même phase mais non miscibles, par exemple : eau
boueuse, pâte à papier, eau-pétrole, gaz-pétrole, etc... ;
• d’autre part, de la forme physique du mélange : concentration des divers constituants,

forme du mélange par micro-bulles en suspension ou par macro-inclusions (poches), etc... ;
toutes ces propriétés étant non-isotropes et non permanentes en un point donné.
Par analogie avec les écoulements turbulents, on peut écrire que τ0 est constitué d’une valeur
moyenne et d’une composante fluctuante dépendant de l’hétérogénéité du mélange :
τ 0 = τ 0 + τ'0
τ 0 est la valeur moyenne spatiale =

1
L
∫ τ0 . ds
0

La composante fluctuante étant caractérisée par une moyenne spatiale nulle à condition que
celle-ci soit effectuée sur une échelle importante :
τ'0 =
1
L
∫ τ'0 ds = 0
0
L’hétérogénéité du mélange est caractérisée par le terme H :
( τ' ox 4 + τ 'oy 4 + τ 'oz 4 )

H=
τ0
r
valeur correspondant au module du vecteur H de composantes :
τ 'ox ²
τ0
→ τ 'oy ²
H
τ0
τ 'oz ²
τ0
Ainsi, l’équation de Navier-Stokes s’écrira, en s’inspirant de la forme (46) :
→ → 1 → → 1 →
Γ= f − grad p + ν∇ ² C + grad H (49b)
ρ ρ
Autre formulation :
L’expression (46) correspondant au fluide incompressible peut être transformée en
introduisant l’expression A-12 de l’accélération :
→
→ ∂C → C² → → →
Γ= + grad − C ∧ rot C
∂t 2
Dans le cas où les forces de volume dérivent d’un potentiel tel que la pesanteur :
→ →
f = − grad (gh )

La relation de Navier-Stokes s’écrit :
→
∂C → → → →  C² dp  →
∂t
− C ∧ rot C = − grad  + gh +
2 ∫ ρ 
+ ν . ∇ ² C
r
Introduisons le vecteur tourbillon Ω défini par :
→ 1 → →
Ω= rot C
2
→ →
et prenons le rotationnel des deux membres de cette équation. En retenant que rot (grad) = 0 et
→ → → → → →
que ∇² C = rot ( rot C) = rot ( 2 Ω) , on obtient :
∂ → → → → → →
→ →
( rot C) − rot ( C∧ 2Ω) = ν rot rot 2Ω
∂t  
→
∂Ω → → → →
+ rot (Ω ∧ C ) = ν∇ ² Ω (50)
∂t
C’est une des formes alternatives de l’équation de Navier-Stokes. C’est la forme

tourbillonnaire établie par Helmotz.
Elle indique que l’écoulement irrotationnel satisfait l’équation de Navier-Stokes mais ne peut
satisfaire les conditions aux frontières. A la surface d’un corps solide, il se forme une couche
tourbillonnaire constituant la couche limite qui diffuse ensuite dans l’écoulement et se
transforme en chaleur.

IV - CARACTERISTIQUES DE L’ECOULEMENT VISQUEUX

L’observation de l’écoulement des fluides est aujourd’hui facilitée par des moyens de mesures
de plus en plus performants tels que le fil chaud et, plus récemment, l’anémomètre laser.
La connaissance de la structure fine des écoulements et en particulier des phénomènes

turbulents, des couches limites, des sillages a permis de mettre en place des modélisations de
plus en plus performantes.
Sur la base des expériences de Reynolds, un pas important a été franchi à la fin du dernier
siècle, concernant la structuration des régimes d’écoulement. Prandtl et Von Karman ont
présenté au début du vingtième siècle les premiers développements théoriques concernant la
couche limite. Celle-ci fait encore l’objet de nombreuses études, tant théoriques qu’expéri-
mentales, portant essentiellement sur la transition et le décollement en présence de gradients
de pression positifs ou négatifs et du taux de turbulence.
4.1. Régimes d’écoulement
L’expérience de Reynolds consiste à observer l’écoulement de l’eau dans un tube transparent,

écoulement au centre duquel est introduit un fin filet de liquide coloré.
On observe principalement deux régimes distincts :
- un régime laminaire montrant un filet coloré parfaitement rectiligne ne se mélangeant pas à

l’écoulement principal. Ce régime est observé aux faibles vitesses, il est parfaitement stable
et pratiquement imperturbable.
- un régime turbulent se manifestant à partir d’une certaine distance de l’entrée par une
oscillation du filet coloré puis à un mélange intime avec le fluide principal.
En utilisant plusieurs diamètres de tube et différentes viscosités de fluide (en faisant varier par
exemple la température), Reynolds a mis en évidence que le régime laminaire était toujours
observé lorsque le produit C.D/ν était inférieur à 2300. Ce produit adimensionnel, appelé
nombre de Reynolds, quantifie l’ensemble des caractéristiques de l’écoulement : la vitesse
moyenne du fluide, le diamètre de la conduite, la viscosité cinématique du fluide.
Au-delà de cette valeur, suivant la qualité de l’entrée du tube, la régularité de la surface

intérieure, la présence de vibrations, etc..., l’écoulement laminaire peut ou non se transformer
en régime turbulent. Dans des conditions idéales, l’écoulement laminaire peut se maintenir
jusqu’à un nombre de Reynolds de 4.104. Au-delà, l’écoulement organisé se dégrade
systématiquement en écoulement turbulent.
4.2. Propriétés des écoulements laminaires
L’écoulement laminaire se rapproche de l’idée que l’on se fait de l’écoulement d’un fluide
parfait pour lequel on imagine un glissement des diverses couches de fluide les unes sur les
autres sans échange de particules entre les différentes couches.
Dans ce régime d’écoulement, les effets visqueux sont pourtant prépondérants et on est très
loin de l’écoulement potentiel, en particulier au voisinage des parois.

Cependant, par opposition à l’écoulement turbulent, c’est l’organisation et la stabilité qui

dominent parmi les propriétés principales de l’écoulement laminaire.
Les applications industrielles les plus importantes correspondent aux paliers et patins
hydrodynamiques et à tous les domaines de la lubrification.
Ils correspondent également à l’écoulement en milieu poreux : filtration, applications
pétrolières, etc...
L’écoulement laminaire unidimensionnel voire bidimensionnel fournit quelques unes des

solutions exactes de l’équation de Navier-Stokes.
Ecoulement entre deux plaques parallèles
L’écoulement s’effectue comme indiqué sur la Figure 9 entre deux plaques planes parallèles
de largeur très importante par rapport à la longueur l. La distance séparant les plaques est
égale à 2h. Suivant les cas, il peut exister ou non un gradient de pression et la plaque
inférieure peut être fixe ou se déplacer parallèlement à elle-même (vitesse de translation U0).
Figure 9 : Ecoulement de Poiseuille
Reprenons l’équation de continuité et l’équation de Navier-Stokes (47) exprimée en

coordonnées cartésiennes.
Hypothèses simplificatrices :
• Les composantes v et w de l’écoulement suivant y et z sont nulles.

∂
• L’écoulement est permanent, ce qui implique que les termes en sont nuls.
∂t
• On suppose les forces de pesanteur négligeables (écoulement horizontal par exemple).
Compte
r tenu du caractère unidimensionnel de l’écoulement, l’équation de continuité div
C = 0 se réduit à :
∂u
=0
∂x

indiquant que la vitesse u n’est fonction que de y :
u = u(y) (51)
Dans ces conditions, les projections de l’équation de Navier-Stokes se réduisent aux trois
équations suivantes :
1 ∂p ∂² u
0=− +ν (52)
ρ ∂x ∂y ²
1 ∂p
0=− (53)
ρ ∂y
1 ∂p
0=− (54)
ρ ∂z
Les équations (53) et (54) indiquent que la pression ne dépend que de x, l’équation du
mouvement se réduit alors à :
1 dp d ²u
=ν (55)
ρ dx dy ²
La composante u ne dépendant que de y, l’équation (55) montre que le premier membre ne

peut être fonction de x et, par conséquent, doit être nul ou constant :
dp ∆p
En posant : =− (56)
dx l
avec ∆p = p(x = 0) - p (x = l) (57)
l’équation (55) s’intègre en :
∆p y ²
u=− + C1 y + C 2 (58)
µl 2
où C1 et C2 sont des constantes d’intégration dépendant des conditions aux limites.
Conditions aux limites :
Les conditions aux limites sont données sur la Figure 9 et traduisent les conditions
d’adhérence du fluide à la paroi :
y = - h, u (- h) = U0
y = + h, u (+ h) = 0

On obtient ainsi, d’après (58), deux équations permettant de calculer C1 et C2 :
∆p h ²
U0 = − − C1h + C 2
µl 2
∆p h ²
0=− + C1h + C 2
µl 2
donnant pour C1 et C2 :
U0 U0 ∆ p h ²
C1 = − et C2 = +
2h 2 µl 2
d’où l’équation générale du mouvement :
U0  y  ∆p h ²   y 2 
u= 1− + 1 −    (59)
2  h  µl 2   h  
Un certain nombre de situations peuvent se rencontrer :
A) Ecoulement de Poiseuille
La plaque inférieure est immobile (U0 = 0) et l’écoulement s’effectue sous l’action du

gradient de pression ∆p/l. La distribution de vitesse est parabolique et donnée par l’équation
(60) :
∆p h ²   y  2 
u= 1 −    (60)
µl 2   h  
Figure 10 : Distribution de vitesse en écoulement laminaire

La vitesse maximale est obtenue au centre du canal et vaut :
∆p h²
u max = (61)
µl 2
La vitesse moyenne spatiale ou vitesse débitante est obtenue à partir du calcul du débit sur la
largeur L :
qv = ∫ u L dy = 2h L u (62)
−h
∆p h² 2
u= = u (63)
µl 3 3 max
Pour une largeur L, le débit est calculable par la relation (64) montrant que celui-ci est
proportionnel à la pression motrice ∆p :
2 ∆p 3
q v = 2h L u = h L (64)
3 µl
Le nombre de Reynolds de l’écoulement est calculable à partir du diamètre hydraulique DH de

la section. Il est donné par la relation :
4S
DH = (65)
P
S est la section de passage : 2hL

P est le périmètre : 4h + 2L
Soit DH pour un écoulement tel que h << L :
8hL
DH = = 4h
2L
d’où le nombre de Reynolds :
R = 4 hν u (66)
L’écoulement reste laminaire et l’on observe une distribution de vitesse parabolique tant que
le nombre de Reynolds reste inférieur à la valeur critique de 2300 à 40000 suivant les
conditions d’établissement.

Contraintes au sein du fluide :
L’état de contrainte au sein de l’écoulement est donné par l’équation de Newton (21) :
du
τ=µ
dy
Relation dans laquelle la direction y désigne la direction allant de la paroi vers le centre de
l’écoulement.
Cette contrainte s’annule au centre de l’écoulement où le terme du/dy s’annule (sommet de la

parabole) et se trouve maximale à la paroi (Figure 11) :
du du
τp = µ . y = − h = − µ. y=h
dy dy
∆p
τp = .h (67)
l
La contrainte interne est linéaire du centre à la paroi.
Figure 11 : Distribution des contraintes dans l’écoulement
Coefficient de perte de charge :
La pression motrice à mettre en jeu pour assurer l’écoulement est souvent ramenée à la
pression dynamique engendrée par l’intermédiaire d’un coefficient λ appelé coefficient de
perte de charge :
u²
∆p = λ . ρ (68)
2
La pression motrice est égale à la perte de charge par frottement sur les parois et se calcule
par l’intermédiaire de la contrainte de paroi τp.

A l’équilibre, le bilan des forces exercées sur le volume fluide compris entre l’entrée et la
sortie s’écrit :
∆p . 2h L = 2 τp l . L
∆p
τp = h
l
équation évidemment identique à celle établie en (67).

On introduit également le coefficient de perte de charge linéique Λ défini par la relation :
l l
λ=Λ. =Λ. (69)
DH 4h
La relation (68) permet de calculer Λ :
4h ∆p
Λ=
l u²
ρ
2
En exprimant u à partir de (63) et en introduisant le nombre de Reynolds défini en (66), on
obtient finalement :
96
Λ= (70)
R
B) Ecoulement de Couette généralisé
C’est l’écoulement général décrit par l’équation (59) dans laquelle la vitesse de déplacement
de la paroi inférieure est égale à U0 :

U0 y  ∆p h 2   y 2 
u= 1 − + 1 −    (71)
2 h  µl 2   h  
Suivant la valeur de ∆p, on observe les allures suivantes pour la distribution de vitesse (Figure
12) :
Figure 12 : Différentes distributions de vitesse

On peut noter que le présent écoulement 12-a est susceptible de constituer une pompe capable
de débiter en présence d’une contre-pression (∆p < 0).
4.3. Propriétés des écoulements turbulents
Le passage du régime laminaire au régime turbulent s’effectue dans une conduite donnée par
augmentation du débit.
La connaissance des propriétés des écoulements turbulents est essentiellement d’origine

expérimentale. En régime turbulent, l’écoulement unidimensionnel est instationnaire et
présente des fluctuations dans la direction d’écoulement mais également sur les directions
normales. Ces fluctuations de caractère aléatoire (régies par les lois du hasard) prennent
naissance de l’existence au sein du fluide de tourbillons de taille variable en interaction
permanente.
Une étude théorique basée sur l’intégration de perturbations dans l’équation de Navier-Stokes,
montre que dans la couche limite, celles-ci sont soit amorties pour de faibles nombres de
Reynolds, soit amplifiées pour donner à l’aval un régime turbulent.
En un point quelconque d’un écoulement unidimensionnel, l’enregistrement des valeurs

instantanées de la composante u présente l’allure de la Figure 13.
Figure 13 : Fluctuations de la vitesse en fonction du temps
On présente généralement la vitesse comme la somme de la composante moyenne u et d’une

composante fluctuante u’ (pouvant atteindre localement 20 à 30% de u ) :
u = u + u’ (72)
La composante moyenne u représente la moyenne temporelle établie par :

t0 + T
u=
1
T ∫ u . dt (73)
t0
Pour un écoulement permanent en moyenne, sa valeur reste constante si le temps T
d’observation est suffisamment important.

Par contre, la composante fluctuante possède une moyenne temporelle nulle :
t0 + T
u' =
1
T ∫ u' . dt = 0 (74)
t0
Les fluctuations de vitesses se produisent également sur les directions normales à

l’écoulement principal, les composantes v et w ne sont plus nulles mais égales à leurs
composantes fluctuantes v’ et w’.
v = v’
(75)
w = w’
Egalement de moyennes temporelles nulles, v’ et w’ satisfont à la relation (74).
En réalisant les moyennes quadratiques, on observe que celles-ci sont non nulles :
t0 + T
u' ² =
1
T ∫ u' ² . dt > 0 (76)
t0
On définit ainsi l’intensité de turbulence I par le rapport :
1
3
( u' ² + v' ² + w' ² )
I= (77)
u
Dans le cas d’un écoulement tridimensionnel, le dénominateur sera constitué de la vitesse

moyenne C .
Coefficients de corrélation :
L’analyse des variations du coefficient de corrélation Ru entre les fluctuations de deux

composantes présente un grand intérêt :
u'1 u'2
Ru = (78)
u'12 u'22
Lorsque u’1 et u’2 sont mesurées en deux points différents et que l’on fait varier la distance
séparant ces deux points, on peut définir la corrélation spatiale et la longueur de corrélation
par l’intégrale :
L= ∫ R u (l) dl (79)
0
représentant l’aire sous la courbe de la Figure 14.

En portant l’évolution de Ru en fonction de la distance l séparant les points de mesure, on

observe l’évolution suivante : Ru = 1 pour l = 0 et Ru = 0 pour l → ∞.
Figure 14 : Corrélation spatiale pour les fluctuations de même nature
Pour une corrélation proche de 1, la distance l0 définit la taille des petits tourbillons, L définit
l’échelle des grands tourbillons.
Dans une conduite de diamètre D, l’ordre de grandeur de l0 est de 0,03 D.

La corrélation peut porter également sur deux composantes orthogonales de la vitesse :
u'1 v '2
R uv = (80)
u'1 ² v '2 ²
Il existe en tout 9 coefficients de corrélation.
Analyse spectrale :
S’agissant de composantes fluctuantes, l’analyse spectrale permet de décomposer le signal

sous la forme d’une série de Fourier et de reporter l’amplitude des différentes sinusoïdes en
fonction de leur fréquence f.
On observe en général un spectre monotone et continu sans fréquence caractéristique.
Si on analyse le terme u' 2 , on obtiendra l’énergie cinétique des fluctuations longitudinales de

la vitesse (analyse de Taylor).
Les formes les plus classiques sont données par la Figure 15.

Figure 15 : Spectre énergétique de la turbulence
Interprétation :
La turbulence s’entretient par prélèvement d’énergie sur l’écoulement principal, assurant la

création de grands tourbillons qui se disloquent en formant des tourbillons plus petits jusqu’à
une échelle telle qu’ils sont dissipés par viscosité.
4.4. Modélisation de la turbulence
En écoulement turbulent, l’influence des termes fluctuants devient plus importante que celle
des termes visqueux déjà examinés.
Leur introduction dans les équations nécessite une modélisation qui reste encore aujourd’hui
de caractère semi-empirique.
Equation de continuité :
En se référant au fluide incompressible, l’équation de continuité instantanée s’écrit :
∂u ∂v ∂w
+ + =0
∂x ∂y ∂z
Comme chacune des composantes est constituée de la somme de deux termes, on calculera les
moyennes des différents termes de la manière suivante :
t + ∆t t + ∆t
0 0
∂u 1 ∂u 1 ∂
=
∂x T ∫
t0
∂x
dt =
T ∂x ∫ u dt
t0

Equation dans laquelle on reconnaît la définition (73) de u :
∂u ∂ u
=
∂x ∂ x
L’équation de continuité s’écrit alors :
∂ u ∂ v ∂w
+ + =0 (81)
∂x ∂y ∂z
En remplaçant chaque terme par sa valeur issue de (72), on observera également :
∂u' ∂v' ∂w '

+ + =0 (82)
∂x ∂y ∂z
montrant que les fluctuations respectent également l’équation de continuité.
Equation dynamique :
Toujours en fluide incompressible, considérons la projection sur x de l’équation de Navier-

Stokes sous la forme donnée par (43).
∂p ∂τxx ∂τxy ∂τxz

ρΓx = ρfx − + + + (83)
∂x ∂x ∂y ∂z
∂u ∂u ∂u ∂u
avec Γx = +u +v +w
∂t ∂x ∂y ∂z
∂u
Le terme instationnaire a une moyenne nulle :
∂t
t0 + T
∂u 1 ∂
=
∂t T ∫
t0
∂t
u dt = [u ( t 0 + T) − u ( t 0 )] = 0
1
T
(84)
En effet, lorsque l’intervalle de temps T est suffisamment grand, la moyenne du terme

instationnaire est nulle pour un écoulement permanent en moyenne.
On démontre que les autres termes se mettent sous la forme :
∂u ∂u' ² ∂u ∂u' v' ∂u ∂u' w '

Γx = u − +v − +w −
∂x ∂x ∂y ∂y ∂z ∂z
∂u' ² ∂u' v' ∂u' w'

Γx = Γx − − −
∂x ∂y ∂z

qui donne finalement comme projection :
∂p ∂ ∂ ∂
ρΓx = ρfx − +
∂x ∂x
( τxx − ρu' ²) +
∂y
( τxy − ρu' v') +
∂z
( τxz − ρu' w ') (85)
Le tenseur des contraintes est donc modifié puisqu’au tenseur des contraintes visqueuses est
soustrait le tenseur des contraintes turbulentes appelé tenseur de Reynolds.
τ = 2 µ D − ρ τt (86)
Rappelons la forme de ces tenseurs :
D : tenseur des taux de déformation donné en (31)
u' ² u' v' u' w'
et τ t = v' u' v' ² v' w' (87)

w' u' w' v' w' ²
Chacune des contraintes du tenseur correspond à un coefficient de corrélation entre les

composantes fluctuantes de la vitesse en un même point de l’écoulement.
Sauf au voisinage immédiat des parois où les termes ∂ u /∂y sont importants, les contraintes
turbulentes sont généralement supérieures aux tensions visqueuses. On s’autorise ainsi à
concevoir l’écoulement turbulent d’un fluide comme étant constitué de zones distinctes : les
zones visqueuses constituées des couches limites où l’écoulement subit à la fois l’influence
des tensions visqueuses et des tensions turbulentes et enfin la zone dite « saine » ou ne
subsistent que les tensions turbulentes.
La modélisation des fluctuations de vitesses est donné par une relation semi-empirique établie
par Prandtl montrant que la fluctuation locale est directement proportionnelle au gradient de
vitesse moyenne locale.
On posera donc :
du
u' = l . (88)
dy
où l est la longueur de mélange de Prandtl.
On admet également que v’ est du même ordre de grandeur que u’ mais de signe opposé, cette
estimation étant justifiée par l’équation de continuité des composantes fluctuantes (82).
Ainsi, la tension turbulente s’écrira :
du du
τ t xy = − ρ u' v' = ρl ² (89)
dy dy

du
Le terme ρl² est donc assimilable à une viscosité que l’on qualifie de viscosité turbulente.
dy
Toujours au titre des modèles empiriques, Von Karman propose comme longueur de mélange
une fonction de la distance y à la paroi :
l = ky (90)
La contrainte turbulente devient alors :
du du
τ t xy = ρk ² y ² (91)
dy dy
L’expérience montre que k est voisin de 0,41.

V - COUCHE LIMITE BIDIMENSIONNELLE

L’une des propriétés des fluides visqueux newtoniens ou non est d’imposer une condition
d’adhérence à la surface d’une paroi solide.
r
r r le fluide visqueux, lar condition à la paroi s’écrira donc C = 0 alors qu’elle s’écrit
Pour
C . n = 0 en fluide parfait, n représentant la normale locale à la paroi.
Ce ralentissement important n’occupe en réalité qu’une épaisseur réduite par rapport aux
dimensions du corps en contact avec le fluide et, très rapidement, en s’éloignant de la paroi,
on atteint les valeurs de vitesses égales ou comparables à celles obtenues en fluide parfait
(Figure 16).
Cette zone à fort gradient de vitesse est qualifiée de couche limite. Son épaisseur δ varie
continûment depuis le point d’arrêt et son régime d’écoulement passe du laminaire au
turbulent en présentant entre les deux un régime de transition.
Figure 16 : Développement de la couche limite à la surface d’un corps solide
La couche limite se sépare ensuite du corps solide vers le bord de fuite en formant, en aval de
celui-ci, un sillage généralement très turbulent qui se dissipe dans l’écoulement principal.
En dehors de la couche limite, l’écoulement est pratiquement non visqueux puisque les
gradients de vitesses sont relativement modestes.

5.1. Equations régissant la couche limite laminaire
En s’imposant un repère local dont l’axe Ox est tangent à la paroi et l’axe Oy est normal et
dirigé vers l’extérieur, les équations bidimensionnelles d’un écoulement incompressible et
stationnaire s’écrivent en négligeant les forces de volume :
∂ u ∂v
+ =0 (92)
∂ x ∂y
∂u ∂u 1 ∂p  ∂² u ∂² u 
u +v =− + ν +  (93)
∂x ∂y ρ ∂x  ∂x² ∂y² 
∂v ∂v 1 ∂p  ∂² v ∂² v 
u +v =− + ν +  (94)
∂x ∂y ρ ∂y  ∂x² ∂y² 
Un certain nombre de simplifications ont été apportées, permettant d’accéder aux équations
caractéristiques dites équations de Prandtl.
On définit pour cela un certain nombre de références permettant d’adimensionnaliser les

variables.
Ue : vitesse de référence (par exemple la vitesse à l’extérieur de la couche limite dans la

direction x).
L: longueur caractéristique dans la direction x
δ: longueur caractéristique dans la direction y
δ
On retiendra : << 1
L
δ
Ue : vitesse de référence dans la direction y
L
ρUe² : la pression de référence
L
: le temps de référence
Ue
Il est alors possible de déterminer l’ordre de grandeur de chacun des membres des équations
du mouvement.

∂u ∂v
+ =0 (92)
∂x ∂y
Ue Ueδ Ue
=
L Lδ L
∂u ∂u 1 ∂p ∂²u ∂² u
u +v =− +ν +ν (93)
∂x ∂y ρ ∂x ∂x² ∂y ²
U U ² δ U e Ue ² Ue ² U U
Ue . e = e Ue . = ν. e ν. e
L L L δ L L L² δ²
∂v ∂v 1 ∂p ∂²v ∂² v
u +v =− +ν +ν (94)
∂x ∂y ρ ∂y ∂x² ∂y ²
Ueδ δ U eδ 1 Ue ² Uδ 1 Uδ 1
Ue Ue . . ν. e . ν. e
L² L L δ δ L L² L δ²
Pour l’équation (92), les deux termes subsisteront tandis que dans les termes visqueux de (93),
seul le second terme sera conservé (fort gradient transversal et faible gradient longitudinal).
Dans l’équation (94), seul le terme de pression et le second terme de viscosité subsistent.
Si l’on effectue pour l’équation (93) le rapport entre le terme de pression et le terme visqueux,
on obtient :
Ue ² δ² Ueδ² Ue L  δ 2
. = = . 
L νUe L. ν ν L
En effectuant la même opération pour l’équation (94), on aura :
Ue ² δL UeL
. =
δ νUe ν
montrant que le terme visqueux du second membre est négligeable devant le terme de
pression pour l’équation (94).
Compte tenu de ces simplifications, les équations de Prandtl de la couche limite laminaire
sont données ci-dessous :
∂u ∂v
+ =0
∂x ∂y
∂u ∂u 1 ∂p ∂² u
u +v =− +ν (95)
∂x ∂y ρ ∂x ∂y ²
∂p
=0
∂y

La dernière équation indique que la couche limite transmet à la paroi la pression existant à la
périphérie de la couche à la même abscisse :
p (y = 0) = p (y = δ)
C’est un point très important concernant les prises de pression effectuées à la paroi.
En présence d’une forte courbure de rayon RC, on utilisera :
1 ∂p u²
= (96)
ρ ∂y R C
Autre expression :
A la périphérie de la couche limite, l’écoulement est pratiquement non visqueux et peut se

calculer selon la théorie potentielle.
Selon le théorème de Bernoulli (14), la vitesse à la limite Ue(x) est liée à la pression pe(x) par
la relation :
Ue² (x)
pe( x) + ρ = cte
2
Sachant que pe ne dépend pas de y, on peut écrire :
1 dpe dUe( x)
− = Ue (97)
ρ dx dx
d’où l’expression finale :
∂u ∂v
+ =0
∂x ∂y
∂u ∂u dUe ∂² u (98)
u +v = Ue . +ν
∂x ∂y dx ∂y ²
où Ue est une fonction de x décrite par la théorie potentielle.
• A la paroi : y = 0 et u = v = 0
• A la périphérie de la couche limite : y = δ(x) et u = Ue(x)

A) Couche limite laminaire à la surface d’une plaque plane
C’est un cas très intéressant pour lequel une solution a été proposée par Blasius au début du
vingtième siècle.
On imagine une plaque plane lisse et infiniment mince disposée parallèlement à la direction
de l’écoulement.
dUe
L’écoulement potentiel est de la forme Ue(x) = cte, ainsi le terme est nul (Figure 17).
dx
En introduisant la fonction de courant ψ, les composantes de la vitesse prennent la forme :
∂ψ
u=
∂y
(99)
∂ψ
v=−
∂x
La fonction de courant satisfaisant l’équation de continuité, l’écoulement dans la couche

limite est alors régi par une seule équation issue de (98) :
∂ψ ∂2ψ ∂ψ ∂2ψ ∂3ψ

. − . 2 =ν 3 (100)
∂y ∂y∂x ∂x ∂y ∂y
Pour simplifier cette équation, on introduit la variable adimensionnelle s et la fonction f(s) :
Ue y
s= y. = Rx (101)
νx x
ψ = νx Ue . f (s) (102)
Rx est un nombre Reynolds construit avec la distance parcourue à la surface de la plaque

plane à partir du point d’arrêt :
Rx = Ue x (103)
ν

Figure 17 : Evolution de la vitesse à la surface de la plaque plane
Les relations (99) et (100) s’écrivent alors :
Ue df df
u = νxUe . . = Ue (104)
νx ds ds
1 νUe  df 
v=  s . ds − f 
(105)
2 x
d ²f d 3f
f. +2 3 =0 (106)
ds² ds
df
A la paroi : s=0 =0 f=0 (d’après 104 et 105)
ds
Sachant que dans ce cas précis, la vitesse ne varie pas de y = δ(x) jusqu’à l’infini, on écrira :
df
s→∞ u = Ue et d’après (104) : =1
ds
L’intégration numérique de la relation (106) permet d’accéder aux caractéristiques locales de

la couche limite laminaire. La fonction f est tabulée en Annexe R17.

Profil de vitesse :
u df
Le profil de vitesses est obtenu en traçant l’évolution de = en fonction de s (Figure
Ue ds
18).
Figure 18 : Profil de vitesse de la couche limite laminaire

à la surface d’une plaque plane
Le profil théorique de Blasius est comparé avec les résultats d’essais de Dhawan (NACA TN
2567-1952).
Contrainte à la paroi :
Celle-ci est calculée à partir de la relation de Newton :
∂u
τp = µ
∂y y = 0
que l’on peut calculer à partir des relations (104) et précédentes :
Ue d ² f
τ p = µ . Ue . . (0) (107)
νx ds²
On définit le coefficient de frottement à la paroi par le rapport adimensionnel ci-dessous :

τp
Cf = (108)
1
ρ Ue²
2
En remplaçant τp par sa valeur issue de (107) et en tirant de l’annexe R17 la valeur de

d ²f
(0) = 0,33206 , on obtient finalement :
ds²
0, 664
Cf = (109)
(Rx )
où le nombre de Reynolds Rx est défini en (103).
Epaisseur de la couche limite :
Par convention, on admet que la limite de la couche frontière est atteinte lorsque la vitesse
u(y) atteint 99% de la valeur atteinte à la limite, à savoir :
u
= 0,99 (110)
Ue
L’Annexe 17 et la Figure 18 indiquent que cette valeur est obtenue pour s = 5.

On aura ainsi, d’après (101), en désignant par δ cette épaisseur :
νx
δ=5 (111)
Ue
ou encore :
5x
δ= (112)
Rx
Utilisant l’équation (112), on obtiendrait les valeurs de δ suivantes :
Pour de l’air et pour x = 1 m :
δ = 19 mm pour Ue = 1 m/s
Pour de l’eau et pour x = 1 m :
δ = 2,9 mm pour Ue = 3 m/s

Comme le montre la relation (112), lorsque x augmente, l’épaisseur de couche limite

augmente aussi. Cependant, à partir d’un certain nombre de Reynolds Rx, l’écoulement dans
la couche limite devient transitoire puis turbulent et les relations proposées ne sont plus
applicables.
B) Cas général
La méthode décrite précédemment nécessite la connaissance de l’écoulement potentiel Ue(x)

dUe
au voisinage d’une paroi pour donner un ordre de grandeur au terme Ue figurant dans
dx
l’expression générale (98).
Falkner et Skan ont étendu la méthode précédente au cas général :
Ue(x) = Co xm (113)
Relation dans laquelle Co est une constante et où m peut prendre des valeurs réelles positive,
négative ou nulle (ramenant au cas de la plaque plane étudié précédemment).
La relation (113) correspond à l’écoulement potentiel sur un dièdre pour lequel l’axe Ox est
confondu avec la face étudiée (Figure 19).
Figure 19 : Ecoulement sur un dièdre

Rappel : (chapitre A - §2.4.)
La fonction potentiel F(z) correspondant à l’équation :
zm+1
F(z) = Co (114)
m+1
s’écrit en coordonnées cylindriques :

r [ cos ( m + 1) θ + i sin ( m + 1) θ ]
Co m+1
F( z) =
m +1
En l’identifiant à la forme générale φ + iψ, on peut déduire la fonction de courant ψ :
Co m+1
ψ= r sin ( m + 1) θ (115)
m +1
La forme du profil correspondant à ψ = 0 est obtenue pour :
mπ
sin (m + 1) θ0 = 0 ou θ0 = avec m entier.
m+1
β
En remarquant que θ0 est identique à π , on retient pour toutes les valeurs de m :
2
2m
β= (116)
m+1
• m = 0 et β = 0 correspondant à une plaque plane disposée dans la direction de

l’écoulement.
• m = 1 et β = 1 correspondant à une plaque plane perpendiculaire à l’écoulement.
Les composantes de la vitesse sont données par :
 ∂ψ
 Cr = = Co r m cos ( m + 1) θ
r∂θ

 ∂ψ
 Cu = − ∂r = Co r sin ( m + 1) θ
m
Si enfin on effectue un changement de repère correspondant à une rotation d’angle (m + 1)θ,

les deux composantes deviennent :

 Cr = Co r m

 Cu = 0
Soit avec les notations retenues jusqu’ici :
Ue = Co . xm (117)
identique à la relation (113) proposée.
• β > 0 et m > 0 correspondent à une accélération de l’écoulement en fonction de x et à un

gradient de pression négatif.
• β < 0 et m < 0 correspondent à une décélération et à un gradient de pression positif.

Lorsque ce gradient atteint une valeur suffisante, on obtient un décollement de la couche
limite de la paroi.
En introduisant s et f(s) sous les formes suivantes, similaires à (101) et (102) :
y
s = m +1 Rx (118)
x
νx Ue . f (s )
1
ψ= (119)
m +1
On obtient alors les relations :
∂ψ df
u= = Ue (120)
∂y ds
1 ν Ue  m − 1 df 
v=− m + 1 f + s  (121)
2 x  m + 1 ds 
L’équation de Navier-Stokes s’écrit alors :
2 m   df  
2
d ²f d 3f
f +2 3 + 1 −    = 0 (122)
ds² ds m + 1   ds  
La résolution numérique proposée par Hartree, ainsi que les conditions aux limites, sont
identiques à celles retenues précédemment, le terme m est un paramètre.
La Figure 20 donne les valeurs du profil de vitesse pour diverses valeurs de m.

Figure 20 : couche limite laminaire sur différents dièdres
Parmi ces caractéristiques, la plus importante est celle correspondant à m = - 0,0904.
 ∂u 
Pour cette valeur limite, on observe   = 0 , conduisant à l’annulation de la contrainte
 ∂y  y = 0
de cisaillement à la paroi. Cette condition est celle qui sera retenue comme critère de
décollement de la couche limite laminaire.
Le demi-angle correspondant est égal d’après (116) à - 0,3121 radian soit encore à - 17.88°.
u
Pour m = 1 (Figure 21), l’épaisseur de couche limite est donnée pour = 0,99 , soit pour
Ue
s = 3,4 correspondant d’après (118) à :
δ Co x²
2 = 3,4
x ν
ν
soit : δ = 2,4 (123)
Co
Valeur constante en fonction de x et de y.

Figure 21 : Couche limite à la surface d’un plan perpendiculaire à l’écoulement
5.2. Evolution de la couche limite laminaire à la surface d’une paroi de forme

quelconque
Lorsque l’écoulement potentiel n’est pas de la forme Ue(x) = Co xm, on qualifie la couche
limite de non-autosimilaire. C’est le cas que l’on doit traiter lors de l’étude de l’écoulement
autour de formes comme celles des profils aérodynamiques.
La méthode la plus sommaire consiste à remplacer les résultats obtenus aux différentes
abscisses curvilignes par ceux obtenus pour la plaque plane aux mêmes abscisses linéaires
pour les mêmes valeurs de m.
L’abscisse x est définie à partir du point d’arrêt sur l’obstacle.
Figure 22 : Définition locale de la couche limite

A l’abscisse x, on définit la tangente à la paroi par rapport à la vitesse à l’infini Co. Cet angle
π
est égal à β (Figure 19), d’où l’on calcule le coefficient β.
2
La relation (116) permet d’en déduire m :
β
m= (124)
2−β
L’épaisseur de la couche limite δ(x) sera alors définie à partir des valeurs de s correspondant à
u
= 0,99 illustrées en Figure 20 :
Ue
0,437
s99 = 3,028 + (125)
0,1987 + m
Partant de cette valeur, on peut déduire l’épaisseur de couche limite δ en se référant à la

définition de s donnée en (118) :
s99 1
δ= x. (126)
m+1 Ue x
ν
Cette application nécessite évidemment la connaissance de Ue(x) qui aura été déterminée au
préalable grâce à l’écoulement potentiel en fluide parfait. Le chapitre C permet de résoudre ce
type de problème lorsque le profil est défini géométriquement.
Si la couche limite reste laminaire, le décollement s’effectuera à partir du point D, formant

l’angle de - 17.88° avec la direction de Co (Figure 22).
Cette méthode présente l’avantage d’être rapide et permet d’obtenir avec une précision
acceptable l’épaisseur de couche limite δ. Aucune indication n’est fournie concernant le profil
de vitesse. Celui-ci n’est accessible que par intégration numérique de l’équation de Prandtl
fournie en (98).
Contrainte à la paroi :
La contrainte à la paroi τp est donnée par la relation de Newton :
∂u
τp = µ
∂y y = 0

Compte tenu de la définition de s donnée en (118) et de la définition de u donnée en (120) :
y Ue x
s= m+1
x ν
df
u = Ue
ds
∂u
Il est possible de déterminer τp à partir du calcul de :
∂y
∂u d ² f ∂s 1 Uex d ² f
= Ue . = Ue m + 1 . (127)
∂y ds² ∂y x ν ds²
et d’en déduire Cf d’après (108) :
τp 2 m + 1 d2f
Cf = 1 = 2 (0)
(128)
ρUe2 Ue x ds
2 ν
Les valeurs de la pente à l’origine sont disponibles en fonction de m (Figure 23).
Figure 23 : Evolution du gradient de vitesse à la paroi en fonction de m

Cette fonction peut être modélisée par la corrélation suivante :
d ²f
(0) = 0,841 ( m + 0,0904)0,387
ds²
conduisant à la relation générale :
m+1
Cf = 1,682 ( m + 0,0904)0,387 (129)
Ue x
ν
On retrouve pour la plaque plane (m = 0) la valeur classique donnée par (109) :
0,664
Cf =
Ue x
ν
Pour une définition plus précise, on se reportera aux travaux de Von Karman, Polhausen et
Thwaites.
Remarque :
Pour améliorer la précision, et selon le développement de la théorie de la couche limite auto-

similaire, divers auteurs [105], ont proposé de remplacer l’angle β défini géométriquement
par un développement de la fonction potentielle Ue(x) :
β=
2 dUe
Ue ² dx
∫ Ue(x) dx (130)
0
En adoptant pour Ue la forme Co xm, on retrouve l’expression déjà établie :
2m
β=
m+1
Cette démarche donne des résultats tout à fait acceptables et plus précis que ceux obtenus en
remplaçant β par sa valeur géométrique.

5.3. Autres paramètres caractéristiques de la couche limite
D’autres paramètres caractéristiques sont utilisés dans l’étude et la modélisation de la couche

limite. En plus de l’épaisseur δ, on introduit δ1 l’épaisseur de déplacement et δ2 l’épaisseur de
quantité de mouvement, grandeurs utilisées dans les équations dynamiques intégrées.
Epaisseur de déplacement δ1 :
Elle correspond au déficit de débit lié à la présence de la couche limite. On écrit, par
définition :
δ1 Ue = ∫ ( Ue − u) dy (131)
0
Figure 24 : Epaisseur de déplacement de la couche limite
Comme u ne varie pratiquement plus à partir de δ, la relation (131) permet de calculer δ1

suivant la forme :
δ1 = ∫ 
1 −

u 
 dy
Ue 
(132)
0
Epaisseur de quantité de mouvement δ2 :
On définit de la même façon l’épaisseur δ2 correspondant du bord d’attaque jusqu’à l’abscisse

x, la variation de quantité de mouvement :
ρ U e2 δ2 = ∫ ρu ( Ue − u) dy (133)
0
d’où la valeur de δ2 :
δ2 = ∫ u 
1 −
Ue 
u 
 dy
Ue 
(134)
0
qui représente la perte de quantité de mouvement associée à la couche limite.
On observe toujours l’inégalité :
δ > δ1 > δ2
Equation intégrale de Von Karman :
Von Karman a intégré l’équation de Navier-Stokes de la paroi y = 0 jusqu’à l’épaisseur y = δ,

on obtient ainsi une équation différentielle ordinaire en δ2 [réf. 105-1985].
d dUe τp
( U e2 . δ2 ) + Ue . δ1 = (135)
dx dx ρ
En divisant par Ue2, on fait apparaître le coefficient de frottement Cf et on obtient :
d δ2 Cf δ2 dUe
= − (H + 2) (136)
dx 2 Ue dx
où H représente le paramètre de forme :
δ
H= 1 (137)
δ2
Polhausen a introduit un deuxième paramètre de forme K :
δ2 dUe
K= 2
ν dx
qui permet d’écrire l’équation intégrale de la couche limite sous la forme (138) :
d  δ2 
Ue  2  = F( K) (138)
dx ν

Walz et Thwaites ont enfin proposé une forme très simple pour F(K) permettant l’intégration
de (138) :
F(K) = 0,45 - 6K (139)
donnant l’expression générale :
x
δ22
ν
= 0,45 . Ue −6 ∫ Ue5 . dx (140)
0
Thwaites a proposé un tableau numérique permettant de traiter avec une excellente précision
les principaux cas lorsque la fonction Ue(x) est connue.
Ce tableau numérique reporté en Annexe R18 permet, comme au chapitre précédent, de

calculer les différentes épaisseurs de couche limite et les valeurs de la contrainte de paroi sans
pour autant définir le profil de vitesse. Le décollement de la couche limite s’effectue pour
K = - 0,09.
dUe
Dans le cas de la plaque plane où = 0 , l’équation de Von Karman s’écrit d’après (136) :
dx
x
δ2 = ∫ Cf
2
dx
0
En retenant l’expression (109) pour Cf, on obtient finalement :
δ2 0, 664
= (141)
x Rx

5.4. Couche limite turbulente
La couche limite laminaire est caractérisée par un écoulement sans fluctuations de vitesse
locale. Lorsque le nombre de Reynolds Rx prend une valeur suffisante, on observe
l’apparition de bouffées turbulentes de faible durée et de faible amplitude longitudinale.
En introduisant un facteur d’intermittence γ qui est la fraction du temps correspondant au

régime turbulent, on peut définir la zone de transition par le critère suivant :
Elle débute lorsque γ = 2% et se termine lorsque γ vaut 90%. La couche limite est alors
considérée comme turbulente. Le début et la fin de cette zone de transition dépendent
beaucoup de la turbulence de l’écoulement hors couche limite mais également de la courbure
et de la rugosité de la paroi.
Dans de nombreux cas, on néglige la longueur de la zone de transition en supposant que la

couche limite change brutalement de régime en passant du laminaire au turbulent.
Plusieurs corrélations de prévision de transition sont disponibles, la plus simple concerne la

plaque plane pour laquelle on observe :
Rx = 3 . 105 à 3 . 106 (142)
Du début à la fin de la transition, une forte variation du gradient de vitesse de paroi est
constaté. Il se traduit par une chute du facteur de forme et un accroissement tout aussi brutal
du coefficient de frottement Cf (figure 25).
Figure 25 : Evolution de H et Cf au voisinage dans la zone de transition

Les travaux de Smith concluent sur une zone de transition définie par la démarche suivante :
A chaque étape du calcul de la couche limite, δ1 et δ2 sont accessibles ainsi que le nombre de
Reynolds associé à δ2 :
Rδ 2 = Ue δ 2 (143)
ν
La zone stable (laminaire) est définie par :
R δ2 < ( R δ2) tr (144)
 22400 
avec (Rδ 2 ) tr = 1,174 1 +  Rx 0.46 (145)
 Rx 
La transition xT apparaît à l’intersection des deux nombres de Reynolds en fonction de

l’abscisse x :
Figure 26 : Définition de la transition selon Smith
Equation de la couche limite en régime turbulent
Les équations de la couche limite turbulente bidimensionnelle incompressible se réduisent à

des expressions comparables à (98) dans lesquelles apparaissent les tensions turbulentes
jusqu’ici absentes. Elles s’écrivent à partir des moyennes temporelles des composantes de la
vitesse :

∂u ∂ v
+ =0
∂x ∂y
∂u ∂u ∂² u ∂ (146)
u
∂x
+v
∂y
=−
1 dp
ρ dx
+ν −
∂y ² ∂y
( u' v')
Les tensions turbulentes difficilement modélisables sont maintenant prépondérantes sauf au

voisinage immédiat de la paroi.
La contrainte de cisaillement au sein du fluide est, d’après (146), formée de deux termes :
∂u
τ = µ. − ρ u' v' (147)
∂y
Profils de vitesses associés à la plaque plane
Sous-couche laminaire :
A la paroi, le terme ρu' v' est nul et l’on retrouve :
du
τp = µ (148)
dy y = 0
L’influence de la viscosité ne se fait sentir que dans une zone appelée sous-couche laminaire
inférieure à 5% de l’épaisseur de la couche limite δ.
Dans cette zone, on définit la vitesse de paroi u* ou vitesse de frottement :
τp
u*= (149)
ρ
La combinaison de (148) et (149) permet d’écrire :
u u*y
= (150)
u* ν
On effectue souvent un changement de variables en introduisant les paramètres

adimensionnels suivants :
u yu *
u+ = , y+ = (151)
u* ν
L’équation du profil de vitesse s’écrira alors :

u+ = y+ (152)
.
Ce profil de vitesse est valable pour la sous-couche laminaire délimitée par :
y+ < 10 (153)
.
δ δ
Région proche de la paroi : <y<
98 10
Par extension, on adoptera dans cette zone une autre loi de paroi de la forme :
u+ = f(y+) (154)
.
La fonction f sera établie en partant des ordres de grandeur des contraintes visqueuses et
turbulentes dans la couche limite (figure 27).
Figure 27 : Evolution de la contrainte de cisaillement dans la couche limite
Dans la zone de paroi, la valeur de τ ne diffère pratiquement pas de τp.
En introduisant la longueur du mélange de Prandtl, on obtient la forme établie en (91) :
 du 
2
τ p = ρk² y²   (155)
 dy 

u
Soit en introduisant la vitesse de paroi u* et les termes adimensionnels u + = et
u*
u*y
y+ = :
ν
1 dy+
du+ = (156)
k y+
Soit, après intégration, la première loi logarithmique :
u + = A log10 y + + B (157)
A et B sont des constantes introduites d’après l’expérience. La forme finale de la loi de paroi
s’écrira finalement :
u + = 5,85 log y + + 5,56 (158)
Elle est utilisable dans l’intervalle :
30 < y + < 500 (159)
On utilise aussi la forme (160), issue également de l’expérimentation :
u+ = 8,74 . y1+/ 7 (160)
Reste de la couche limite :
Dans la zone restante comprise entre δ/10 et δ, on introduit un profil de vitesse appelé loi
déficitaire et dépendant cette fois du gradient de pression extérieur dp/dx.
Ue − u y
= A log + C (161)
u* δ
pour 30 < y+ < 1000 (162)
La relation (161) se superpose pratiquement à (158).
Au-delà, on utilise :
Ue − u  y b
= a 1 −  (163)
u*  δ

où a et b sont des constantes dépendant du gradient de pression extérieure (b est voisin de 2).
Pour y = δ, la forme (163) satisfait la condition à la limite u = Ue .

On formule souvent, sous la forme simplifiée d’une loi en puissance, l’ensemble de la
distribution de vitesse dans la couche limite turbulente :
1/ n
u  y
=   avec n ≈ 7 (164)
Ue  δ 
Relations relatives à la plaque plane :
Les relations établies dans la zone de paroi permettent d’introduire un certain nombre de
relations importantes.
La définition de la vitesse de frottement permet de faire apparaître Cf :
τp
u*=
ρ
τp
et Cf = 1
ρ U2e
2
Soit en éliminant τp :
u* Cf
= (165)
Ue 2
Introduisons cette expression dans les différentes lois de paroi en puissance (160) et (164) :
u+ = 8,74 . y1+/ 7
u  y 1/ n
= 
Ue  δ 
Développons (160) :
u Ue  y . u *1/7
. = 8,74 . 
Ue u *  ν 
En remplaçant u* par sa valeur issue de (165), on peut introduire le coefficient Cf :
1/7
u  Cf +1/2 y  Cf 1/2 
=  . 8,74  Ue .   
Ue  2  ν  2  

En identifiant cette valeur avec celle issue de (164), pour n = 7, on obtient :
2 ν 1/ 4
Cf =  
(8,74) 7 / 4  Ue . δ 
 ν 1/ 4
soit Cf = 0,045 .   (166)
 Ue . δ 
Epaisseur de quantité de mouvement δ2 :
Comme en laminaire, l’expression de δ2 se calcule à partir de (134) :

δ
δ2 = ∫ u 
1 −
u 
Ue  Ue 
 dy
0
qui donne, d’après (164) :
δ2 n
= (167)
δ ( n + 1) ( n + 2)
et pour n = 7 :
δ2
= 0,0972 (168)
δ
On démontre que l’équation intégrale de Von Karman conserve son expression dans le cadre
de la couche limite turbulente, les termes de fluctuation de vitesse apparaissant sous la forme
d’une différence pratiquement négligeable :
 δ 
∫( )
2
( ) U 
d δ2 Ue2 + δ1 Ue dUe = Cf e dx + d 
2

u'2 − v'2 dy
 0 
Ainsi, pour la plaque plane où dUe = 0 :
d δ2 τp Cf
= =
dx ρ Ue² 2
En remplaçant δ2 par sa valeur issue de (168) et Cf issu de (166), on obtient :
dδ  Ue δ −1/ 4
= 0,231 .  
dx  ν 

En intégrant de 0 à x :
δ  Ue x −1/ 5
= 0,370 .   (169)
x  ν 
soit finalement, pour l’écoulement turbulent :
0, 058
Cf = (170)
(Rx )1/ 5

5.5. Couche limite turbulente sur une surface quelconque
A) Méthode générale
La définition complète des paramètres de la couche limite dépend comme en régime

laminaire du gradient de pression. Cette influence s’exerce sur le profil de vitesse et les
conditions de décollement.
Le point de départ est constitué de l’équation intégrale de quantité de mouvement qui s’écrit
en régime turbulent :
 δ 
 
∫
d ( δ2 Ue ²) + δ1 Ue dUe = Cf Ue² ds + d  ( u' ² − v' ²) dy  (171)
 0 
Les termes u' ²et v' ² étant du même ordre de grandeur, le dernier terme est négligeable par
rapport aux autres, réduisant ainsi l’équation intégrale à la même expression qu’en laminaire :
Cf
2 δ2 Ue dUe + Ue² d δ2 + δ1 Ue dUe = Ue² ds
2
dδ2 δ dUe Cf
+ ( H + 2) 2 = (172)
ds Ue ds 2
Cette équation contient trois inconnues : δ2, H et Cf, le terme Ue(s) étant supposé connu par
une étude en écoulement potentiel où s représente l’abscisse curviligne. La résolution sera
possible en introduisant deux relations complémentaires, l’une évaluant Cf et l’autre la forme
du profil de vitesse définissant H (cette forme dépend ici du gradient de pression).
Coefficient de frottement Cf :
Sous une forme simplifiée, on peut retenir la relation de Prandtl établie à partir de (166)
et(168) pour la plaque plane, elle est indépendante de dp/ds :
1
Cf  Ueδ2 − 4
= 0,0128   (173)
2  ν 
Pour plus de précision, on adoptera la relation de Ludwieg-Tillmann basée sur des résultats
expérimentaux en présence d’un gradient de pression non nul :
 Ueδ 2 −0,268
Cf = 0,246 . 10 −0,678H   (174)
 ν 
D’autres méthodes couplent à la fois le profil de vitesse et le coefficient de frottement par le

paramètre G de Rotta-Clauser.

Profil de vitesse :
L’influence du gradient de pression sur le profil de vitesse est très important comme le
montrent les mesures de Clauser (1954).
Figure 28 : Influence du gradient de pression sur le profil de vitesse
Le gradient de pression est décrit par le paramètre D dont la valeur nulle correspond à la
plaque plane :
δ dp
D= 1 (175)
τ p ds
Deux formes sont utilisées pour décrire l’influence de D :
• la forme de Le Foll [réf. 209] retenant la loi déficitaire suivante :
y y
g1   + λ g2  
Ue − u δ δ
= (176)
u* k
 y  1  y 2  y  y  y 3  y 2
avec : g1   =   − 1 − ln   ; g2   =2   − 3  + 1 ; k = 0,41
 δ  2  δ   δ δ δ δ
et λ paramètre dépendant implicitement de D par les relations de Huo.

• La forme en puissance beaucoup plus simple :
1
u  y n
=  (177)
Ue  δ 
pour laquelle le gradient de pression est pris en compte par le facteur n.
On a déjà examiné que n = 7 pour D = 0 (plaque plane). En reprenant les définitions

classiques de δ1 et δ2, on obtient :
δ
δ1 1 H −1
=
δ δ
∫ 
1 −
u 
 Ue 
 dy =
1
=
n +1 H +1
(178)
0
δ
δ2 1 H −1
=
δ δ
∫ 
1 −
u  u

 Ue  Ue
dy =
n
=
( n + 1) ( n + 2) H ( H + 1)
(179)
0
δ n+2 2
H= 1 = ou n= (180)
δ2 n H −1
Equations donnant deux formes à la loi en puissance :
H −1
u y 2
=  (181)
Ue  δ 
H −1
u  y H −1  2
et =  (182)
Ue  δ2 H ( H + 1) 
Equation définissant le paramètre de forme :
Le paramètre G de Rotta-Clauser est défini par le rapport suivant :
δ
 Ue − u 2
∫ 
 u* 
 dy
G= 0δ (183)
 Ue − u 
∫ 
 u* 
 dy
0
Quel que soit le profil de vitesse retenu, on démontre qu’il peut se mettre sous la forme :

H −1 1
G= . (184)
H Cf
2
Suivant la forme de profil retenu, il prendra aussi l’expression :
• Pour la forme de Le Foll :
1,244 + 115
, λ + 0,3714 λ2
G= (185)
2 λ
k + 
3 2
• Pour la forme en puissance :
1 2
G= (186)
Cf n + 2
2
Huo a montré que G et D étaient reliés par la relation :
G = 7,55 D + 18
, − 355
, pour G ≤ 13 (187)
G = 6,10 D + 1,8 - 0,40 pour G > 13 (188)
Mode opératoire :
A partir du point de transition, on connaît les valeurs de xT, δ2, δ2, δ1 et H. Partant de ces
valeurs, on peut définir successivement :
Cf. par la relation (174)
D par la relation (175)
G par la relation (187) ou (188)
H par la relation (184)
λ ou n par les relations (185) ou (186)
Le profil de vitesse est alors parfaitement défini par (176) ou (177). En se fixant un pas ∆s
suffisamment faible, on déduit la nouvelle valeur de δ2 en utilisant l’équation (172) :
∆δ2 δ dUe Cf
+ ( H + 2) 2 . =
∆s Ue ds 2
soit δ2 = δ2 + ∆δ2
et ainsi de suite jusqu’à l’extrémité de la paroi ou jusqu’au point de décollement.

B) Méthode de Head
Head a intégré l’équation de continuité dans la direction y de 0 à δ :
δ δ
∂u
∫ ∂x
dy + ∫ ∂∂vy dy = 0 (189)
0 0
On développe d’abord le premier terme :
δ δ
d
ds
∫ u dy =
d
ds
∫ Ueu . Ue dy = Ue dsd (δ − δ1) + dUe
ds
( δ − δ1) = [ Ue ( δ − δ1)]
d
ds
0 0
 δ − δ1 
Par hypothèse, le second terme est défini comme une fonction de   , autre facteur de
 δ2 
forme de la couche limite noté H1 :
δ − δ1
H1 = (190)
δ2
L’équation de continuité s’écrit alors :
d
ds
( Ue δ2 H1) = Ue . F ( H1) (191)
On lui associe l’équation de Von Karman :
dδ2 Cf δ2 dUe
= − (H + 2) (192)
ds 2 Ue dx
En précisant la forme des fonctions F, H1 et Cf en fonction de H, le système (191) - (192) se

ramène à deux équations différentielles en δ2 et H qui peut s’intégrer numériquement avec
pour conditions initiales les valeurs données par la transition.
Les formes retenues sont les suivantes :
H1 = 0,8234 ( H − 11, )−1,287 + 3,3 ; H < 16,

−3,064
H1 = 15501
, ( H − 0,6778) + 3,3 ; H > 16
,
−0,6169 (193)
F = 0,0306 ( H1 − 3)
 Ueδ2 −0,268
Cf = 0,246   . 10−0,678.H
 ν 

C) Décollement de la couche limite turbulente
La valeur moyenne correspondant au décollement de la couche est donnée par H = 3 au lieu

de 3,55 pour la couche limite laminaire. Le décollement ne se produit pas précisément pour
Cf = 0 mais pour un Cf décroissant très rapidement.
La couche limite turbulente décolle plus difficilement de la paroi que la couche limite
laminaire à cause du fort échange énergétique existant avec l’écoulement extérieur.
5.6. Ordre de grandeur des caractéristiques de la couche limite sur une plaque plane
A titre d’application et d’illustration, on retient une plaque plane de largeur infinie et de

longueur égale à 200 mm, soumise à un écoulement d’air ( ν = 15 E-6 m²/s et
µ = 18 E-6 kg/ m.s ) de vitesse Ue = 150 m/s.
Les résultats des différents modèles proposés en laminaire par Thwaites ou déduits des
travaux de Falkner et Skan, et en turbulent par Head sont présentés ci-dessous.
Arbitrairement le critère de transition retenu est ici : R x = 1,72 106
Comme on l’a déjà noté en figure 25, le passage du régime laminaire au régime turbulent
s’effectue avec un accroissement important du coefficient de frottement Cf.
Caractéristiques de la couche limite sur une plaque plane.

VI) ECOULEMENTS EN CONDUITE
Il s’agit d’un écoulement interne par excellence où la couche limite sur les parois occupe
pratiquement toute la section d’écoulement.
L’écoulement est régi par les équations de Navier-Stokes et sont conditionnées, en fluide
incompressible newtonien, par deux nombres caractéristiques : le nombre de Reynolds et le
nombre de Froude.
Etudions le cas général en disposant l’axe z vertical, les forces du volume étant supposées
réduites à l’action de la pesanteur. Exprimées par unité de masse, elles se réduisent à leur
seule composante verticale (équation 7) :
0
→
f = 0 (200)
−g
Pour éviter d’alourdir la démonstration, considérons dans un premier temps les termes
convectifs négligeables.
Le système d’équation de Navier-Stokes (47) se ramène alors à :
du 1 ∂p
=− + ν ∇²u
dt ρ ∂x
1 ∂p
0= − + ν ∇²v (201)
ρ ∂y
1 ∂p
0= − − g + ν∇ ² w
ρ ∂z
On peut adminensionnaliser les variables en utilisant comme précédemment les grandeurs de

référence suivantes.
C: La vitesse de référence correspondant par exemple à la vitesse débitante dans la

section de conduite.
D: La longueur de référence correspondant à une dimension géométrique caractéristiques,
soit le diamètre de la conduite ou son diamètre hydraulique (relation 65).
D/C : le temps de référence
ρC² : la pression de référence
Les variables réduites s’écriront alors :
u v w
u+ = v+ = w+=
C c C
x y z
x+ = y+ = z+ =
D D D

t C. t p
t+= = , p+ =
D D ρC²
C
Les différents termes aisément calculables :
du C² du +
=
dt D dt +
1 ∂p C ² ∂p +
− =−
ρ ∂x D ∂x +
C
ν ∇² u = ν . ∇² u +
D²
D
En multipliant tous les termes par , on obtiendra le système en variables réduites :
C²
du + ∂p + ν
=− + ∇²u +
dt + ∂x + CD
dv + ∂p + ν
=− + ∇²v + (202)
dt + ∂y + CD
dw + ∂p + gD ν
=− + + ∇²w +
dt + ∂z + C² CD
système dans lequel apparaît le nombre de Reynolds de l’écoulement :
R = C.D (203)
ν
qui représente le rapport des forces d’inertie aux forces de viscosité.
et le nombre de Froude :
F = C² (204)
gD
qui représente le rapport des forces d’inertie aux forces de pesanteur.

Pour les fluides non newtoniens dont la loi de comportement est de la forme :
∂u
τ = τ0 + µ
∂y

où τ0 correspond aux tensions interfaciales prenant naissance par exemple dans un mélange de
fluides non miscibles ou d’une émulsion, un terme complémentaire figurera dans l’équation
de Navier-Stokes (relation 49b) :
1 →
grad H
ρ
En introduisant la tension interfaciale σ comme contrainte de référence, l’équation de Navier-

Stokes contiendra en plus le nombre de Weber :
W = ρ C²
σ (205)
Les nombres de Reynolds et de Froude servent de base aux règles de similitude des
écoulements. L’influence du nombre de Froude n’est sensible que dans les écoulements à
surface libre. Par contre, pour les écoulements en charge, c’est le nombre de Reynolds qui est
prépondérant puisqu’il fixe à la fois le régime d’écoulement mais également les gradients de
vitesse à la paroi, la forme de la distribution de vitesse, etc...

6.1. Régime laminaire
Nous étudions dans ce chapitre l’écoulement laminaire dans un tube de section circulaire en
régime incompressible et permanent (écoulement de Poiseuille).
Reprenons l’équation de continuité et l’équation de Navier-Stokes en coordonnées

cylindriques (A-18) et (B-48) en retenant les points suivants :
• les forces de pesanteur sont supposées négligeables,
• les composantes Cr et Cu de l’écoulement sont nulles,
• l’écoulement est axisymétrique et les termes de vitesse en ∂/∂θ sont nuls,
• le régime d’écoulement est établi, ce qui signifie que si l’on se trouve loin de l’entrée du
tube, ou des phénomènes de mise en vitesse se produisent [réf. 101].
L’équation de continuité s’écrit :
∂Ca
=0 soit Ca = Ca(r)
∂z
L’équation de Navier-Stokes s’écrit :
1 ∂p
0=−
ρ ∂r
1 ∂p
0=−
ρ r∂θ
1 ∂p  ∂²Ca 1 ∂Ca 
0=− + ν +
ρ ∂z  ∂r ² r ∂r 
Les deux premières équations indiquant que la pression ne dépend que de z, d’où :
∂p dp
=
∂z dz
L’équation dynamique se réduit à :
1 d  dCa  1 dp
ν r =
r dr  dr  ρ dz
(206)
La composante Ca ne dépendant que de r et la pression ne dépendant que de z, l’équation

(206) indique que dp/dz ne peut être que constant (fluide en mouvement avec Ca ≠ 0) ou nul
(fluide au repos avec Ca = 0).
dp ∆p
En posant =− (207)
dz l
avec l : longueur de la conduite et ∆p la chute de pression ou perte de charge sur la longueur l.

L’équation (206) s’intègre en donnant :
d  dCa  ∆p
 r  =− r
dr  dr  µl
dCa ∆p
r =− . r ² + C1
dr 2 µl
dCa ∆p C1
=− .r +
dr 2 µl r
∆p r ²
Ca ( r ) = − + C1 . ln r + C2 (208)
2 µl 2
C1 et C2 sont des constantes d’intégration permettant de satisfaire les conditions aux limites
qui sont ici ramenées à une seule condition : la vitesse est nulle à la paroi du tube.
Pour obtenir une valeur finie sur l’axe du tube, il faut nécessairement C1 = 0.
La condition de paroi s’écrit alors pour un tube de rayon R :
∆p
Ca (r = R) = R² + C2 = 0
4 µl
D’où l’équation finale :
∆p
Ca ( r ) =
4 µl
( R² − r ²) (209)
Le profil de vitesse est donc parabolique et symétrique par rapport à l’axe de la canalisation.
Par la suite, pour simplifier, nous noterons la vitesse axiale par C au lieu de Ca.
Figure 29 : Distribution de vitesse en régime laminaire

Propriétés de la distribution de vitesse :
La vitesse maximale est obtenue sur l’axe de la canalisation :
∆p . R²
Cmax = (210)
4 µl
Le débit volume est calculé par l’intégrale suivante :
R R
∆p
qv = ∫ C . dS = ∫ 4 µl
( R² − r ²) 2π r dr (211)
0 0
∆p R 4
soit qv = π (212)
4 µl 2
On observe donc que le débit dans la canalisation est proportionnel à la différence de pression
s’exerçant entre l’entrée et la sortie.
Le débit s’écrivant aussi sous la forme du produit de la vitesse moyenne (ou vitesse débitante)
par la section de passage, nous aurons par définition :
qv = C . πR² (213)
donnant comme vitesse moyenne par identification de (212) et (213) :
1 ∆p . R ²
C= (214)
2 4 µl
La vitesse moyenne est donc égale à la moitié de la vitesse au centre de la conduite.
Contrainte au sein du fluide :
Pour un fluide Newtonien, l’état de contrainte au sein du fluide est donné par l’équation :
dC dC
τ=µ =−µ (215)
dy dr
Equation dans laquelle la direction y est dirigée de la paroi vers le centre de l’écoulement.
∆p r
Soit : τ= (216)
2l

La contrainte est donc nulle au centre de la conduite et varie linéairement jusqu’à la paroi où
elle atteint la valeur maximale τp :
∆p . R
τp = (217)
2l
Figure 30 : Distribution de contrainte au sein du fluide
La pression motrice totale à mettre en jeu pour assurer l’écoulement est ramenée à la pression
dynamique engendrée pour définir le coefficient de perte de charge.
C2
∆pt = λ . ρ (218)
2
∆pt représente la différence de pression totale entre l’entrée et la sortie de la conduite :
 Ce  
2
Cs 
2
∆pt = pe + ρ  −  ps + ρ  (219)
 2   2 
Dans le cas particulier qui nous intéresse ici, la section est constante et la pression motrice
totale se ramène à la différence de pression.
On introduit également le coefficient de perte de charge unitaire Λ sous la forme :
l C²
∆p =Λ ρ
D 2
où D est le diamètre de la canalisation.

En remplaçant C par sa valeur issue de (214), on obtient pour Λ :
64
Λ= (220)
R
où R représente le nombre de Reynolds calculé comme en (203) par la relation :
R =

(221)
ν
Le coefficient de frottement à la paroi Cf est introduit à partir de la contrainte de paroi et de la

pression dynamique (108) :
τp
Cf = (222)
C²
ρ
2
En remplaçant τp et C par leurs valeurs issues de (217) et (214), on obtient :
Λ
Cf = (223)
4
Comme nous l’avons examiné au § 4.1., l’écoulement laminaire dans une conduite de section
circulaire se maintient jusqu’à un nombre de Reynolds de 2300. Au-delà s’établit, suivant les
conditions d’écoulement, le régime turbulent.

6.2. Régime turbulent lisse
La surface du tube est lisse et on se référera dans cette partie à un certain nombre de résultats
obtenus pour la plaque plane (§ 5.4.).
On démontre que la distribution de contrainte radiale reste proportionnelle au rayon comme

indiqué sur la Figure 30. Compte tenu du modèle proposé pour décrire l’écoulement turbulent,
la contrainte en tout point s’exprime à partir de la longueur de mélange de Prandtl :
dC dC dC
τ=µ + ρl² (224)
dr dr dr
Suivant que l’on est proche ou éloigné de la paroi, c’est le premier terme ou le second qui est
prépondérant.
Dans la sous-couche laminaire, on obtiendra en référence à l’équation (155) une équation de

la a forme :
C+ = y+ (225)
C y . u*
avec C+ = ; y+ =
u* ν
et où u* est la vitesse de frottement donnée par :
τp
u* = (226)
ρ
La valeur de τp est toujours calculable par (217). L’équation (225) est valable pour y+ < 10.
Hors de cette zone et toujours au voisinage de la paroi, on observe la loi logarithmique donnée
pour (157) où les constantes A et B sont légèrement différentes :
C+ = 5,65 log y + + 5,00 (227)
y
Cette relation est valable pour < 0,2 .
R
Dans la zone centrale de la conduite, on fait appel à une loi de type déficitaire :
C max − C  y 2
= 7,2 1 −  (228)
u*  R
Cmax représentant la vitesse au centre de la conduite.

D’autres auteurs utilisent les lois en puissance déjà examinées pour la plaque plane (164) :
1
y n
C
=  (229)
Cmax  R 
Elle est en principe utilisable sur toute la section de la conduite mais n varie avec le nombre
de Reynolds moyen :
C .D 4 . 103 105 2 . 106

R =
ν
n 6 7 10
Calcul du débit et de la vitesse moyenne
Le débit se calcule par la relation d’intégration classique :
qv = C . πR² = ∫ C . 2 πr dr
0
Les zones de paroi étant relativement réduites, on retient pour C la relation (228) sur toute la
section et on remplace y par (R - r).
On obtient la relation :
= 5, 65 log R* + 2

(230)

avec R* = u* . R (231)
ν
Coefficient de perte de charge
Le coefficient Λ peut se déterminer en utilisant les relations précédentes :
Λ = 4 Cf
τp Cf Λ
u* = =C =C
ρ 2 8
1 1 C
soit =
Λ 8 u*

C
En remplaçant par sa valeur issue de (230), on obtient la relation :
u*
1 2
= 2 log R* +
Λ 8
Le nombre de Reynolds R* peut être développé :
Λ
R* = u * R = C . D . u * = R
ν ν 2C 4 2
soit finalement :
1
Λ
= 2 log (R )
Λ − 0,8 (232)
C’est l’équation de Von Karman, forme à résoudre par itération en partant par exemple de
0,02.
Blasius a établi une forme explicite en calculant l’intégrale définissant le débit à partir de la
forme (160) donnée pour la plaque plane :
0,316
Λ= (233)
(R )1/ 4
Elle n’est valable que pour des nombres de Reynolds inférieurs à 105.
Forme de la distribution de vitesse :
Comme l’indique l’ensemble des relations, le profil de vitesse est continûment variable avec
le nombre de Reynolds. Il est de forme beaucoup plus aplatie que dans le cas de l’écoulement
laminaire.
Figure 31 : Distribution de vitesse en régime turbulent

Le rapport entre C et Cmax dépend également du nombre de Reynolds, on le calcule par la

relation de Coantic :
C  R 
= 0, 71 + 0,1 log   (234)
Cmax  4200 

6.3. Influence de la rugosité
L’influence de la rugosité a été étudiée expérimentalement par Nikuradse qui a disposé à la

surface d’un tube de paroi lisse une rugosité standard constituée de grains de sable de
granulométrie déterminée (rugosité uniformément répartie de hauteur ε petite par rapport au
diamètre de la conduite).
Partant de cette rugosité, on peut construire deux termes :

ε
. la rugosité relative
D
ε .u *
. le nombre de Reynolds associé ε + =
ν
Lorsque ε+ est inférieur à 5, les distributions de vitesses ne sont pas perturbées par la rugosité.
Au-delà, la distribution proche de la paroi est affectée et la relation (227) devient :
∆u
C+ = 5,65 log y + + 5 − (235)
u*
∆u
où est une fonction de la rugosité.
u*
Pour ε+ > 70, on vérifie les relations :
∆u
= 5,65 log ε+ − 3,2
u*
D’où la relation en régime hydrauliquement rugueux :
y
C+ = 5,65 log + 8,2 (236)
ε
Quelle que soit la rugosité, la zone centrale n’est pas affectée par la rugosité.
En régime laminaire, Λ est inchangé. En régime turbulent, on peut utiliser la relation de

Colebrook :
1  2, 51 ε 
= − 2 log  +  (237)
Λ  R Λ 3, 71 D 

Pour ε = 0, on retrouve l’équation de Prandtl (232). L’évolution de Λ est donnée sur la Figure
32.
Figure 32 : Evolution de Λ en fonction du nombre de Reynolds

et de la rugosité de la paroi
Dissipation énergétique :
La puissance hydraulique Phy dégradée sous forme de perte de charge se transforme

intégralement en chaleur et correspond à une élévation de température ∆T calculable par :
Phy = ∆p . qv = ρ Cp qv ∆T (238)

Chapitre C
-C-
DEVELOPPEMENTS
DEVELOPPEMENTS
COMPLEMENTAIRES
COMPLEMENTAIRES

Ch.C – Développements complémentaires
I) GRILLE PLANE CIRCULAIRE EN ROTATION
L'écoulement inter-aubages, déjà abordé, est repris ici dans le cadre de la méthode des
singularités.
Dans cette modélisation, on superposera une distribution de tourbillons à la surface des
aubages avec une source centrée en 0 et assurant le débit.
Examinons d'abord la fonction d'influence obtenue par l'existence de trois vortex de
2π
même intensité γ centrés au rayon ro et décalés angulairement de (figure 1).
3
Figure 1 : Définitions
La vitesse complexe C' sera donnée par la sommation suivante (relation A-73) :
iγ 3
1
C' = u − iv =
2π
∑ z − zn
(1)
n =1
avec z affixe d'un point quelconque M

z1 = ro exp (iθ)
 2π 
z 2 = ro exp i θ + 
 3 
 4π 
z 3 = ro exp i θ + 
 3 
L'expression de C' passe par le développement de la somme algébrique :
1 1 1
A= + + (2)
z − z1 z − z2 z − z3
En réduisant au même dénominateur et en groupant les termes, on obtient la fonction

complexe :
3z2 − 2 z ( z1 + z2 + z3 ) + z1z2 + z2 z3 + z1z3

A= 3
z − z2 ( z1 + z2 + z3 ) + z ( z1z2 + z2 z3 + z1z3 ) − z1z2 z3
On démontre en annexe R.15 que la fonction A se réduit pour trois vortex à la

relation :
3z 2
A= (3)
z3 − z13
qui peut être étendue aisément au cas de Na aubages :
Na zNa −1
A = Na (4)
z − z1Na
Si l'on dispose le long de chaque aubage, supposé infiniment mince, la distribution de

tourbillon γ(s), la vitesse conjuguée sera donnée par :
Na z Na−1
S2
i
C' = u − iv = ∫
2π S
γ(s) Na
z − z sNa
ds (5)
1

où s est l'abscisse curviligne le long de l'aubage et zs l'affixe du centre de chacun des

vortex γ(s).
On superposera à cette distribution de singularité une source de débit qv centrée en 0,
caractérisée par le potentiel complexe :
qv
F( z ) = ln z
2π
et de vitesse conjuguée :
dF qv
=
dz 2πz
Cette superposition conduit à l'expression :
Na z Na −1
S2
qv i
C' = u − iv = +
2πz 2π ∫
S
z Na − z sNa
γ(s ) ds (6)
1
La vitesse d'entraînement sera représentée par un tourbillon forcé d'intensité ω centré

en 0 :
U = iωz (7)
avec ω : vitesse angulaire de la roue
On accédera à la vitesse relative en appliquant les règles de composition des vitesses :
→ → →
W=C−U
r
d'où les composantes sur x et y de W :
Re( C' ) + ω r sin θ

→
W (8)
− Im( C' ) − ω r cos θ
Condition de glissement :
La distribution γ(s) sera définie en écrivant la condition de glissement à la surface des
aubes (figure 1) :
− Im( C' ) − ωr cos θ

tg ( θ − β) = (9)
Re( C' ) + ωr sin θ

avec β : angle d'aubage
L'annexe 15 donne les valeurs de Re(C') et Im(C') et la condition de glissement (9)

s'écrit :
 R Na 
S2  r cos θ − Na −1
cos ζ 
qv  r 
2π r
sin θ − ωr cos θ − ∫ R 2 Na
R Na
γ(s ) ds
S 1 r + 2( Na −1) − 2 Na −2 cos ( θ + ζ )
2
tg( θ − β ) = r r (10)
 R Na 
S2  r sin θ − Na −1 cos ζ 
qv  r 
2π r
cos θ + ωr sin θ + ∫ R 2 Na
R Na
γ(s ) ds
S 1 r + 2( Na −1) − 2 Na −2 cos ( θ + ζ )
2
r r
z = r exp (iθ)
avec zs = R exp (iθ')
ζ = Naθ' − ( Na − 1) θ

II) METHODE DE L'HODOGRAPHE
Cette méthode tombée en désuétude mérite une attention particulière dans le cadre de
l'avant-projet. C'est une des rares méthodes inverses qui contrôle, au stade de la
conception, le champ de vitesses en tout point d'un écoulement bidimensionnel
incompressible : écoulement inter-aubages ou entre parois solides.
Son handicap majeur est qu'elle est inopérante dans les zones à fort gradient proches
du bord d'attaque et du bord de fuite, particulièrement importantes au demeurant.
La méthode de l'hodographe est une transformation conforme qui repose sur

l'écoulement potentiel source-puits reprenant pour l'écoulement en projet, la même
organisation, à savoir les mêmes fonctions de courant et potentiel. [réf.227]
a - plan d'étude b - plan réel
Figure 2 : Principe de la méthode

A la source est associé l'écoulement à l'infini amont et au puits l'écoulement à l'infini

aval.
La figure 2 illustre cette démarche.
Comme pour la transformation conforme, on définit le plan d'étude où l'on trace
d'abord les lignes de courant et les équipotentielles de l'écoulement général source-
puits (figures 13 et 14 du chapitre A).
Dans le cadre du problème inverse, les directions des vitesses d'entrée C1 et de sortie
C2 d'un écoulement à définir sont des données. On connaît donc les angles α1 d'entrée
à l'infini amont et α2 de sortie à l'infini aval. On notera (sur la figure 2) la présence
d'un pôle PO et le respect de l'équation de continuité s'écrivant en fluide
incompressible :
C1 cos α1 = C2 cos α2
→ →
La vitesse C 1 est associée à la source S et la vitesse C 2 au puits P.
→
En un point quelconque M de l'écoulement, la vitesse C sera définie par son module C
et sa direction α ; d'autre part, ce même point M est associé à une ligne de courant
ψ = cte et à un potentiel φ = cte.
En se déplaçant sur une même ligne de courant, on pourra écrire pour le plan réel
(relation A - 52) :
∆x2 ∆y2
∆φ = C . ∆s = C = C. (11)
cos α sin α
De même, en se déplaçant sur une même équipotentielle, on pourra écrire également :
∆x2 ∆y2
∆ψ = C . ∆n = C . = −C. (12)
sin α cos α
La construction des lignes de courant et des équipotentielles du plan réel s'effectue

dans l'ordre suivant :

- On se déplace sur la ligne de courant principale ψ = 0 en partant du point O, centre

du repère qui est un point double de la transformation :
point 0 x1 = 0 et x2 = 0
y1 = 0 y2 = 0
- Suivant le maillage retenu pour le plan d'étude, on lit directement le ∆φ

correspondant à un déplacement vers S par exemple.
Pour ce même point O, on relève également les valeurs du module C et de la

direction α de la vitesse.
- On calcule grâce à la relation (11) les valeurs de ∆x2 et ∆y2 :
∆φ cos α
∆x 2 =
C
(13)
∆φ sin α
∆y 2 =
C
Les coordonnées du point 1 de la ligne de courant ψ = 0 sont alors définie.
- On se déplace d'un pas sur la même ligne de courant et on procède de même en

relevant ∆φ, C et α. Les relations (13) permettent de définir les coordonnées du point
2 et ainsi de suite jusqu'à la source S.
- On procède de même en se déplaçant du point O vers le point P et l'on obtient dans le

plan réel les coordonnées des points allant du point O vers l'infini aval.
- La ligne de courant ψ = 0 étant construite, on décrit, toujours à partir de O, les points

correspondant à l'équipotentielle φ = 0.
En se déplaçant d'abord vers les ψ croissants, on relève la valeur de ∆ψ et la relation

(12) permet de calculer les coordonnées du premier point de l'équipotentielle :
sin α
∆x 2 = ∆ψ .
C
(14)
cos α
∆y 2 = − ∆ψ .
C

On observe facilement que, dans le déplacement en ψ croissants, la vitesse C

associée est continuellement croissante. En limitant cette vitesse à une valeur
maximale, on obtiendra la valeur maximale admissible pour ψ, constituant ainsi l'une
des limites de l'écoulement (extrados de l'aubage dans le plan réel).
Inversement, le déplacement en ψ décroissant fera apparaître des vitesses de plus en

plus faibles. En s'imposant une vitesse minimale, on définira ainsi la seconde
frontière de l'écoulement constituant l'intrados de l'aubage dans le plan réel.
- On construit l'ensemble des points correspondant à une autre ligne de courant ou à

une autre équipotentielle en ajoutant les valeurs correspondantes de ∆x2 et ∆y2 aux
coordonnées des points précédemment définis.
En chacun des points du maillage ψ - φ, le module et la direction des vitesses sont

parfaitement définis. On peut alors accéder au champ de pression et aux efforts sur les
aubages.
Suivant les valeurs maximales et minimales admissibles pour la vitesse, on peut noter
que la forme de l'extrados et de l'intrados de l'aubage sont légèrement modifiés et que
les aubages se rapprochent ou s'écartent.
On remarquera enfin que la forme du bord d'attaque et du bord de fuite ne sont pas
définis et restent à l'initiative du concepteur.
La condition de périodicité n'étant pas respectée au bord de fuite, ceci impliquera un

recalage de la direction α2 imposée en sortie (écart angulaire réduisant la déflexion).
Cet écart angulaire, pouvant atteindre plusieurs degrés, rend par trop imprécise cette
méthode en cas de faible déflexion.


Chapitre
- D D-
EXERCICES D’APPLICATION
EXERCICES D’APPLICATION
ET
ET PROBLEMES
PROBLEMES

Ch.D – Exercices d’application et problèmes
EXERCICES SE RAPPORTANT AU CHAPITRE A
Exercice 1 :
On donne le champ de vitesse suivant :
→ u = x 2y + y2
C
v = x 2 − y2x
1 - Vérifier l’équation de continuité.

2 - L’écoulement est-il irrotationnel ?
3 - Déterminer la fonction de courant.
Exercice 2 :
→
Le vecteur C décrit un écoulement cylindrique tournant (sortie d’une grille d’aubes
tournante) :
Cr = 0
→ K3
C Cu = K 2 +
r
Ca = g ( r )
K 2 et K 3 sont des constantes et g(r) est une fonction à définir.
1 - Examiner les conditions nécessaires pour que l’écoulement soit irrotationnel.

2 - Exprimer le potentiel φ.
Exercice 3 :
Un écoulement plan incompressible et irrotationnel est engendré par une source de

débit Q située en (-a, 0), d’un puit de même débit disposé en (+a, 0) et d’un
écoulement uniforme de vitesse C0.
1 - Définir le potentiel complexe F(z).

2 - Déterminer le potentiel φ et la fonction de courant ψ.
3 - Quelle est la forme de l’obstacle ainsi constitué.
4 - Calculer et localiser la vitesse maximale sur l’obstacle.

Exercice 4 :
Le champ de vitesse d’un écoulement hélicoïdal est donné par le vecteur :
Cr = 0
→ K
C Cu = 1
r
Ca = K 2
où K 1 et K 2 sont des constantes.
1 - Parmi les propositions suivantes, encadrer celles qui vous paraissent

justifiées :
Ecoulement permanent
Axisymétrique
Incompressible
Irrotationnel
→
2 - Définir le vecteur accélération Γ
K 12
Γr = − 3
→
r
Γ Γθ = 0
ΓZ = 0
3 - En déduire la distribution de pression radiale en supposant les forces de

volume.
K 12
p (r) = − ρ + cte pour ρ = cte
2 r2

Exercice 5 :
1 - Ecrire le potentiel complexe correspondant à l’association d’un écoulement

→
plan infini Co parallèle à Ox et d’un vortex d’intensité Γ centré en (0, 0)
F (z) =
→
2 - Calculer la vitesse conjuguée C′ et en déduire u et v du champ de vitesses :
u=
v=
3 - Définir les coordonnées du point d’arrêt A :
xA =
yA =

EXERCICES SE RAPPORTANT AU CHAPITRE B
Problème 1 : Caractéristiques de la couche limite à la surface de la plaque plane
• Caractéristiques du fluide : ν = 15.10-6 m² / s

µ = 18.10-6 kg / m.s
• Vitesse de l’écoulement : Ue = 150 m / s
• Longueur de la plaque plane : l = 200 mm
1 - Définir l’abscisse de la transition en retenant pour critère : Rxcr = 1.72 10-6
2 - Dresser le tableau suivant, sur la base de travaux de Blasius.
x (mm) Rx δ (mm) τp (N/mm²) Cf
50
100
150
3 - Comparer ces valeurs avec celles données par la méthode de Thwaites (R17) :
x (mm) δ2 (mm) δ1 (mm) τp (N/mm²) Cf
50
100
150
4 - Calculer, dans la zone turbulente, en adoptant la distribution de vitesse en

puissance (n=7) les valeurs de δ , δ2 et Cf pour x = l.

Problème 2 : Définition de la couche limite à la surface d’un cylindre
Le calcul proposé correspond à l’écoulement bidimensionnel d’un fluide autour d’un

cylindre de rayon R. la couche limite sera supposée entièrement laminaire du point
d’arrêt jusqu’au décollement. On retiendra pour ce calcul les méthodes simplifiées
consistant à utiliser les relations de la plaque plane étendues au dièdre.
1- Dans le cadre de l’écoulement potentiel, définir la vitesse Ue à la surface du

cylindre en fonction de Co la vitesse à l’infini amont et de l’angle φ défini sur la
figure ci-dessous.
2- L’abscisse X, abscisse curviligne à la surface du cylindre, est une fonction de R et

de φ. Montrer que l’expression du coefficient β de l’angle de dièdre, en fonction de
l’angle φ, est donnée par la relation
3- Dresser le tableau ci-dessous dans les conditions suivantes :Co = 10 m/s

R=100 mm
V=15 /s (air)
φ (°) β m Ue m/s X mm Rx ∆ mm Cf
0
20
40
60
80
φ max=
4- Préciser la valeur de correspondant au décollement de la couche limite

( Rappeler la définition de ce phénomène ). Vérifier l’hypothèse selon laquelle la
couche limite est entièrement laminaire.

Problème N°3 : Etude d’un ballon géostationnaire
L’étude proposée concerne le projet d’un ballon géostationnaire de moyenne

altitude destiné au marché des télécommunications locales.
Placé dans la stratosphère, le ballon sera soumis aux vents locaux et devra
garder une position pratiquement stationnaire pour remplir correctement sa mission.
Un système de propulsion, à définir, permettra de réaliser cet objectif.
Nous nous proposons de modéliser l’action du vent sur un tel ballon.
Définition géométrique :
Le ballon gonflé à l’hélium, est un corps ovoïde de révolution présentant une

longueur L de 71 m et un diamètre maximal de 25 m pour un volume total de 43 500
. Les deux extrémités I et II sont raccordées, respectivement, sur des sphères de
diamètre 8 m et 22 m ( figure 1 ).
Conditions de vol :
L’altitude nominale est égale à 20 000 m, altitude où la pression absolue est

égale à 0,05466 bars et où les caractéristiques physiques de l’air sont :
/s
La vitesse maximale du vent local est égal à 25m/s ( quelques heures par an).
C’est cette vitesse, supposée dirigée suivant l’axe Ox, que l’on retiendra comme valeur
dimensionnante pour le propulseur.
I – Etude en fluide parfait
Dans un premier temps, au stade de l’avant-projet, on retiendra les hypothèses

suivantes :
- Le fluide est supposé non visqueux.

- L’étude de l’ écoulement tridimensionnel axisymétrique autour du ballon est
ramenée à celle d’un écoulement bidimensionnel autour d’un profil cylindrique de
section droite identique à celle de la figure 1.
L’écoulement au voisinage de l’extrémité I du ballon est assimilé à la

superposition d’un écoulement uniforme de vitesse Co et d’une source unitaire
centrée à l’origine du repère (figure 2).

1-Rappeler la forme analytique des potentiels complexes de ces deux écoulements.
2 – Calculer la fonction de courant. Montrer que l’équation du ballon dans cette partie
peut s’écrire ( avec θ en radian ) :
3- Montrer que l’expression littérale des composantes de la vitesse à la surface du

ballon est donnée par :
En déduire les composantes du point d’arrêt.
4 – En identifiant ce point d’arrêt avec celui de la figure 2, calculer la valeur
numérique du rapport .
II – Calcul de la couche limite :

La couche limite sera calculée à partir des données de la figure 2 où l’on
relèvera directement, pour les divers points allant de A à 4, les valeurs de l’angle du
dièdre équivalent ainsi que l’abscisse curviligne depuis le point d’arrêt.
5- Déterminer l’ensemble des éléments du tableau suivant en précisant le principe

général du calcul et les hypothèses retenues.
pts Θ Ue Rx β m Cf
rd m/s
A O

6- Commenter l’allure générale de la figure 3 issue d’un calcul complet. Préciser la

position de la transition dans le tableau précédent et en déduire la valeur de Reynolds
critique retenu ici.
Déduire de la même figure le point de décollement, quelle est sa position dans la
région II du ballon.
7- Le ballon est placé en position inverse de manière à présenter l’extrémité II au vent.

La figure 4 donne l’évolution du coefficient de frottement. Quels sont vos
commentaires sur cette figure et quelle conclusion en tirez-vous.
Figure 1 : Définition géométrique du ballon.



Problème N° 4 - Pompe périphérique à viscosité :
Considérons une pompe de modestes caractéristiques dont les propriétés seront

définies à partir des équations de Navier Stokes en écoulement unidimensionnel
( figure 1).
La machine est constituée d’un rotor lisse de rayon R, tournant à la vitesse angulaire
ω. Le canal d’écoulement est constitué d’une section annulaire de hauteur j et de
largeur l. un secteur de largeur négligeable assure l’étanchéïté entre l’entrée à la
pression et la sortie à la pression .
Figure 1 : pompe périphérique à viscosité
Le fluide compris entre le rotor et l’enveloppe extérieure est régi par l’équation de
continuité et l’équation de Navier Stokes écrite en coordonnées cylindriques. Le fluide
est supposé incompressible, l’écoulement permanent, les composantes axiales et
radiale de l’écoulement nulles et les forces de pesanteur négligeables.
Questions :
1- Ecrire l’équation de continuité.

2- Ecrire les projections de l’équation de Navier Stockes en r, θ, z.
3- En admettant, dans l’hypothèse où j est petit devant R, que l’équation du
mouvement se ramène à :
avec

Montrer que Cu peut s’écrire sous la forme :
4- En considérant les conditions aux limites en r = R et r = R+j, déterminer les

constances et
5- Montrer que, par intégration, le débit peut s’écrire sous la forme :
6- Le fluide traversant la pompe est considéré comme fluide Newtonien, la

contrainte visqueuse, s’exerçant sur le rotor peut alors s’écrire
Montrer que la contrainte à la paroi du rotor peut s’écrire de la manière suivante :
7- En déduire le couple résistant s’exerçant sur l’arbre moteur et la puissance

correspondante.
8- Donner l’expression du rendement en fonction du coefficient de débit :
Problème 5 - Comportement d’un hydrofoil orientable
On se propose de réaliser la première étape du dimensionnement d’un ensemble

portant et propulsif destiné à l’équipement d’une vedette de transport rapide.
Deux dispositifs sont mis en place : le dispositif avant I est constitué d’une aile
hydrodynamique de corde I et d’envergure L permettant de supporter sensiblement la
moitié du poids P du bateau et le dispositif arrière II, intégrant le même dispositif
qu’à l’avant avec le propulseur (figure 1). L’objectif poursuivi est d’obtenir, pour la
vitesse de croisière, le déjaugeage (*) du bateau. La résistance à l’avancement dans
l’eau ( majoritaire ), est alors réduite à la seule traînée des dispositifs I et II.

Figure 1 : Dispositif portant et système de propulsion.
(*) Elévation du bateau au-dessus de la surface de l’eau.

I-Dimensionnement global du dispositif portant.
L’aile hydrodynamique orientable est constituée d’un profil Joukowski dont les
caractéristiques géométriques sont les suivantes :
Angle de Cambrure β= - 10°
Coefficient d’épaisseur K =
On néglige les effets d’extrémité et l’on suppose que l’envergure de chaque aile est 4
fois supérieure à la corde ( L=4l).
Les caractéristiques physiques adoptées pour l’eau sont données par :
/s
1- En fluide parfait, définir la corde «l » du profil pour les conditions nominales

suivantes :
Incidence nulle
P = 280 000 N
Co=8 m/s
2-Définir l’angle d’incidence à adopter pour supporter un poids total égal à à la

vitesse d’avance Co=3,56 m/s ( fonctionnement à vide en vitesse réduite ).
II – Approche en fluide réel.
La figure 2 donne le champ de vitesse à la surface du profil en écoulement potentiel

pour l’incidence nulle et l’incidence de 16°. La figure 3 détaille, pour une échelle ½, la
forme de la première moitié du profil.
On se propose de calculer la couche limite à la surface du profil pour deux
configurations d’écoulement. On retiendra comme critère de transition la valeur
supérieure généralement admise.
3- Pour l’incidence nulle et la vitesse d’avance de 8 m/s, établir le tableau suivant pour
les points 1, 2 et 3 définis sur les figures 2 et 3.
pt X Ue Rx Régime δ Cf
m m/s

Pour simplifier, on adoptera pour ce calcul les relations établies pour le dièdre tangent
et on relèvera directement sur les figures 2 et 3 les éléments nécessaires. Commenter
les résultats numériques obtenus.
4- Montrer que sous certaines hypothèses à définir, le décollement de la couche limite

s’effectue au point 4 pour .
III – Cavitation du profil.
Le profil est immergé à 2 m de profondeur. On se place dans l’hypothèse de

l’écoulement potentiel.
5- Expliciter la condition à respecter pour le rapport afin d’éviter la cavitation au

point de pression minimale de la surface au profil. Localiser ce point pour
6-Vérifier la condition pour .
7-A partir de la figure 2, localiser la poche de cavitation lorsque Co. atteint 8 m/s pour

Figure 2 : Champs de vitesse à la surface du profil

Figure 2 : Champs de vitesse à la surface du profil
Figure 3 : Allure du profil à l’échelle 1/2

Problème 6 – Mise en équilibre d’un cylindre
Un cylindre de masse M, de rayon dont la base est parfaitement plane, repose sur
un plan horizontal par l’intermédiaire d’une couche fluide de viscosité dynamique
et d’épaisseur h. Sur ce plan horizontal est situé une cavité circulaire de rayon , un
orifice situé au fond de cette cavité permettra l’alimentation en fluide à une pression po
( voir figure ).
Nous vous proposons d’étudier m’équilibre dynamique en supposant que la base du

cylindre reste constamment parallèle au plan horizontal.
Pour les applications numériques on prendra les valeurs suivantes :
M=1000 kg
1) Dans un premier temps le cyclindre repose sur le plan horizontal sans qu’il ait
écoulement. Calculer la pression nécessaire pour soulever celui-ci.
2) On admet que l’écoulement reste partout en parallèle aux faces, qu’il est
axisymétrique, stationnaire et sans rotation.
a) Quelles sont les conséquences de ces hypothèses sur les composantes des
vitesses et leurs dérivées.
b) Ecrire l’équation de continuité ainsi simplifiée.
c) Montrer qu’elle se réduit à

3) En tenant compte de ces hypothèses,

a) Ecrire les équations de Navier-Stokes.
b) Montrer que la pression qui règne dans la couche fluide ne dépend que du
rayon.
c) Montrer que l’équation dynamique se réduit à la forme :
(1)
4) Dans cette partie, le fluide est supposé non visqueux.
a) Donner l’expression de p(r).

b) Donner l’expression de la vitesse débitante moyenne
en fonction du débit de fuite à travers la levée h.
c) En utilisant les conditions aux limites suivantes ; )=
Ecrire l’équation qui régit la répartition radiale de la pression.

d) En écrivant l’équilibre des forces de poids et de pression
dans la cavité et dans le canal, donner la relation entre la pression dans la
cavité, le débit de fuite et la levée h.
5) Etude de l’influence de la viscosité :
Dans l’équation (1) on pose : (2)

Et on considère les conditions aux limites suivantes :
Conditions d’adhérence sur les parois et
a) A partir de l’équation (2) déterminer l’expression de Cr(z)

b) Donner la nouvelle expression de p(r).

Chapitre E
-E-
ANNEXESET
ANNEXES ET RAPPELS
RAPPELS

ANNEXES
______
R.1. - Produit vectoriel
→ →
Soient deux vecteurs V1 et V2 définis par leurs composantes respectives (coordonnées
cartésiennes) :
u1 u2
→ →
V1 = v1 et V2 = v2
w1 w2
On définit le produit vectoriel par :
→ → v1w 2 − w1v2
V1 Λ V2 = w1u2 − u1w 2
u1v2 − v1u2
→ → → → →
Si V1 et V2 sont colinéaires, V1 Λ V2 = O
R.2. - Produit scalaire
→ →
Le produit scalaire de V1 par V2 est défini par :
→ →
V1. V2 = u1u2 + v1v2 + w1w2 = V1V2 cos α
→ → → → →
Si V1 et V2 sont perpendiculaires, V1 . V2 = O
R.3. - Gradient d'une fonction scalaire ϕ
La fonction ϕ (x, y, z) a pour gradient le vecteur suivant :

Ch.E – Annexes et rappels
∂ϕ ∂ϕ
→
∂x ∂r
∂ϕ → ∂ϕ
grad ϕ = grad ϕ =
∂y r∂θ
∂ϕ ∂ϕ
∂y ∂z
coordonnées cartésiennes coordonnées cylindriques
R.4. - Divergence
u
→
Soit le vecteur V v , on définit sa divergence par le scalaire
w
→ ∂u ∂v ∂w
div V = + + (coordonnées cartésiennes)
∂x ∂y ∂z
R.5. - Rotationnel
→
Le vecteur V a pour rotationnel le vecteur :
∂w ∂v
−
∂y ∂z

→
r ∂u ∂w
rot V = −
∂z ∂x
∂v ∂u
−
∂x ∂y
R.6. - Laplacien
Laplacien d'une fonction ϕ(x, y, z) :
∂²ϕ ∂²ϕ ∂²ϕ

∇ ²ϕ = + +
∂ x ² ∂y ² ∂z ²
∂²ϕ 1 ∂ϕ 1 ∂²ϕ ∂²ϕ

∇ ²ϕ = − + +
∂r ² r ∂r r ² ∂θ² ∂z²

R.7. - Théorème de Green-Ostrogradski

r r
∫∫S nr . AdS = ∫∫∫V div A . dV
R.8. - Différentielle totale
Soit φ une fonction des coordonnées et du temps :
φ = φ ( x, y, z, t )
La différentielle totale de φ sera donnée par :
∂φ ∂φ ∂φ ∂φ
dφ = dx + dy + dz + dt
∂x ∂y ∂z ∂t
Inversement, la fonction K = K1 dx + K2 dy + K3 dz + K4 dt est une différentielle totale

exacte égale à dU si l'on a simultanément :
∂U ∂U ∂U ∂U
K1 = , K2 = , K3 = , K4 =
∂x ∂y ∂z ∂t
R.9. - Changement de repère
r ' = x2 + ( y − a ) 2
y−a
θ' = arc tg
x
R.10. - Rappel
a−b
arc tg a − arc tg b = arc tg
1 + ab
R.11. - Rappel
du
f = arc tg u df =
1 + u2
R.12. – Autres rappels
z (x dx− a) n
=−
1
( n − 1)( x − a ) n − 1
z x dx+ a
2 2
=
1
a
arc tg
x
a
z x xdx
2
+a 2
=
1
2
ln ( x2 + a 2 )
R.13. - En utilisant les relations suivantes :
cos (α + 2β) = cos α cos 2β − sin α sin 2β
sin (α + 2β) = sin α cos 2β + cos α sin 2β
sin 2β = 2 sin β cos β
cos 2β = 1 − 2 sin2 β = 2 cos2 β − 1
On obtient :
Γ
Co 2 cos α sin2 β + 2 sin α sin β cos β = sin β
2 πa
Γ
Co 2 sin α cos β + 2 cos α sin β cos β = cos β
2
2 πa
conduisant à la même relation :
Γ
2 Co sin ( α + β) =
2 πa
soit : Γ = 4πa Co sin (α + β)
R.14. - Mouvement relatif
Ecrivons l'irrotationnalité de l'écoulement absolu :

→ r
rot C = 0

En coordonnées cylindriques, cette relation vectorielle s'écrit :

(1)
1 ∂Ca ∂Cu
− =0
r ∂θ ∂z
→ r ∂Cr ∂Ca
rot C = − =0
∂z ∂r
1 ∂( rCu) 1 ∂Cr
− =0
r ∂r r ∂θ
En remplaçant chaque composante de la vitesse absolue par se valeur en fonction des

composantes relatives et d'entraînement :
0 Wr Cr = Wr
r r r
C = U + W = ωr + Wu = Cu = ωr + Wu
0 Wa Ca = Wa
Le système (1) s'écrit :
(2)
1 ∂Wa ∂Wu
− =0
r ∂r ∂z
→ r ∂Wr ∂Wa
rot C = − =0
∂z ∂r
1 ∂( rWu) 1 ∂Wr
2ω + − =0
r ∂r r ∂θ
r
En observant que le rotationnel de U s'écrit :
→ 0
r
rot U = 0
2ω
Le système (2) peut s'écrire :

→ r →r r
rot C = rot W + 2 ω
ou encore :
(3)
1 ∂Wa ∂Wu
− =0
r ∂r ∂z
→r ∂Wr ∂Wa
rot W = − =0
∂z ∂r
1 ∂( rWu) 1 ∂Wr
− = − 2ω
r ∂r r ∂θ
Le mouvement relatif est donc rotationnel.
uez ici pour taper du texte.



R.21. - Euler 3 D
1 ∂p  ∂Wm ∂δ   Wm  ∂Wm
− = Wm sin δ  sin δ +Wm cos δ  +  tg β + ω sin δ
ρ ∂r  ∂r ∂r   r  ∂θ
 ∂Wm ∂δ  Cu2 ∂Wm
+ Wm cos δ  sin δ +Wm cos δ  − + sin δ − fD cos β sin δ
 ∂z ∂z  r ∂t
1 ∂p  ∂Wm ∂β   Wm  ∂Wm
−
ρ r∂θ
= Wm sin δ 
 ∂r
(
tg β + Wm 1 + tg2β )
+
∂r   r
tg β + ω
 ∂θ
tg β
 ∂Wm ∂β  Cu ∂Wm dω
+ Wm cos δ 
 ∂z
(
tg β + Wm 1 + tg2β
∂z 
)
+
r
Wm sin δ +
∂t
tg β + r
dt
− fD sin β
1 ∂p  ∂Wm ∂δ   Wm  ∂Wm
− = Wm sin δ  cos δ − Wm sin δ  +  tg β + ω cos δ
ρ ∂z  ∂r ∂r   r  ∂θ
 ∂Wm ∂δ  ∂Wm
+ Wm cos δ  cos δ − Wm sin δ  + cos δ − fD cos β cos δ
 ∂z ∂z  ∂t
En utilisant les propriétés géométriques de (28) et en remarquant que :
∂Wm ∂Wm dm dt ∂Wm dt

= . = Wm ,
∂t ∂t dt dm ∂t dm
les composantes s’écrivent :
1 ∂p Wm sin δ  ∂Wm ∂Wm ∂Wm ∂Wm 

− =  dr + rdθ + dz + dt 
ρ ∂r dm  ∂r r∂θ ∂z ∂t 
Wm 2
 ∂δ ∂δ  Cu 2
∂Wm
+ cos δ  dr + dz − +ω sin δ − fD cos β sin δ
dm  ∂r ∂z  r ∂θ
1 ∂p Wm tg β  ∂Wm ∂Wm ∂Wm ∂Wm 

− = dr + rdθ + dz + dt 
ρ ∂r dm  ∂r r∂θ ∂z ∂t 
Wm2  ∂β ∂β  Cu dω ∂Wm
+
dm
( )
1 + tg2β  dr +
 ∂r
dz +
∂z  r
Wm sin δ + r
dt
+ω
∂θ
tg β − fD sin β
1 ∂p Wm cos δ  ∂Wm ∂Wm ∂Wm ∂Wm 

− =  dr + rdθ + dz + dt 
ρ ∂z dm  ∂r r∂θ ∂z ∂t 
Wm 2
 ∂δ ∂δ  ∂Wm
− sin δ  dr + dz + ω cos δ
dm  ∂r ∂z  ∂θ

Les termes entre crochets sont des différentielles totales.
∂Wm ∂Wm ∂Wm ∂Wm

dWm = dr + rdθ + dz + dt
∂r r∂θ ∂z ∂t
∂δ ∂δ
dδ = dr + dz
∂r ∂z
∂β ∂β
dβ = dr + dz
∂r ∂z

Chapitre F
-F-
REFERENCES BIBLIOGRAPHIQUES
REFERENCES BIBLIOGRAPHIQUES

Ch.F- Références bibliographiques
1) OUVRAGES GENERAUX
101 - R. COMOLET : « Mécanique expérimentale des fluides »

1982 - Ed. Masson - Paris
102 - R. OUZIAUX et J. PERRIER : « Mécanique des fluides appliquée »

1972 - Ed. Dunod - Paris
103 - A. HOUBERECHT : « Les turbines »

1972 - Ed. Vander - Louvain
104 - KARASSIK : « Pump Handbook »

1976 - Ed. Mc Graw Hill Co - New York
105 - I. RYHMING : « Dynamique des fluides »

1985 - Ed. Presses Polytechniques Romandes
106 - G.T. CSANADY : « Theory of turbomachinery »

1964 - Ed. Mc Graw Hill Co - New York
107 - G. MEAUZE et A. CARRERE : « Etude aérodynamique de l’écoulement dans

les aubages de machines tournantes »
1984 - Cours à l’Ecole Nationale Supérieure de l’Aéronautique et de l’Espace
108 - F. LEBOEUF, R. VOUILLARMET, P. FERRAND : « Aérothermodynamique

des turbomachines »
Cours de l’Ecole Centrale de Lyon
109 - B. LAKSHINARAYANA : « An assessment of computational fluid dynamic

techniques in the analysis and design of turbomachinery »
1991 - Journal oof Fluid Engineering - ASME - Vol. 113
110 - D. ALTHAUS, F.X. WORTMANN : « Stutgarter profilkatalog I »

1981 - Ed. VIEWEG - Wiesbaden
111 - F.W. RIEGELS : « Results from wind-tunnel investigations »

1961 - Ed. Butterworths - London
Ch.F – Références bibliographiques
112 - F. GIRERT : « Aérodynamique B2 »

1977 - Cours du CNAM
113 - D. ALTHAUS : « Profilpolaren für den modellflug »

1980 - Ed Neckar - Verlag - Villingen
2) PUBLICATIONS ET COMMUNICATIONS
201 - E. BONATAKI, P. CHAVIAROPOULOS, K.D. PAPAILIOU : « An inverse

inviscid method for the design of quasi-three-dimensional turbomachinery
cascade »
ASME - journal of Fluids Engineering - March 1993 - Vol. 115
202 - Document CETIM - HYD/CD - 85 Ni049

Mai 1985
203 - M.H. FAGARD : « Etude des écoulements décollés en grilles d’aubes en

régime sous-sonique compressible »
Thèse de Doctorat - LEMFI - Paris VI - 1993
204 - T.S. LUU, B. VINEY, L. BENCHERIF : « Inverse problem using S2-S1

approach for the design of the turbomachine with splitter blades »
Revue Française de Mécanique - p. 209-224 - n° 1992.3 - FRANCE
205 - A.A. FAY : « Simple performance model for pumps »

Pumping Machinery - ASME - FED. Vol. 154 - p. 119-125 - 1993
206 - M. JAREMCZAK : « Analyse théorique et expérimentale du bruit hydraulique

des pompes centrifuges »
Thèse de Doctorat - Université de Lille - FRANCE - 1992
207 - F.C. VISSER, J.J.H. BROUWERS : « Review of status of potential flow

calculations in radially bladed pump impellers »
Pumping machinery - ASME FED. Vol. 154 - p. 219-226 - 1993

208 - J. CHAUVIN, R.A. VAN DEN BRAEMBUSSCHE, P. LUBUNGA,

C. LOUDET : « Calcul de l’écoulement dans le plan circonférentiel d’une roue
mixte » Institut Von Karman - Rhodes - Belgique - 1973
209 - G. ASSASSA et K.D. PAPAILIOU : « An integral method for calculating
turbulent boundary layers with separations »
1979 - ASME - Journal of fluid engineering - Vol. 100 - p. 110-116
210 - WU C.H. : « A general theory of three dimensional flow in subsonic and

supersonic turbomachines of axial, radial and mixed-flow types »
NACA - TN 2604 - 1952
211 - MARSH M. : « A digital computer program for the through flow fluid
mechanics in an arbitrary turbomachine using a matrix method »
Reports et Memoranda - N° 3509 - 1966
212 - BOSMAN C. et MARSH M. : « Improved method for calculating the flow in

turbomachine, including a consistent loss model »
JSME - Vol. 16 - N° 1 - 1974
213 - BOSMAN C. et EL SHAARAWI . : « Quasi three dimensional numerical

solutions of flow in turbomachines »
ASME - 1977
214 - DAVIS W.R. et MILLAR D. : « A comparison of the matrix and steamline

curvative methods of axial flow turbomachinery analysis from a user’s point of
view »
ASME - 1975
215 - RIBAUD Y. et SOVRANO R. : « Calculs théoriques d’écoulements dans les

turbomachines radiales. Calcul de l’écoulement moyen »
ONERA - NT N° 1/3086 - 1973
216 - NOVAK R.A. : « Streamline curvative computing procedures for fluid-flow

problems »
ASME - 1967
217 - SENOO Y. et NAKASE Y. : « An analysis of flow through a mixed flow

impeller »

ASME - 1972
218 - FROST D. : « A steamline curvature through-flow computer program for

analysing the flow through axial-flow turbomachines »
Reports and Memoranda - N° 3687 - 1970
219 - WILKINSON D.H.
220 - NOVAK R.A. et HEARSEY R.H. : « A nearly three-dimensional intrablade

computing system for turbomachinery »
ASME - 1977
221 - HIRSCH C.H. et WARZEE G. : « A finite element method for the axisymetric
flow computation in a turbomachine »
INJME - Vol. 10 - 1976
222 - HONG TUAN HA : « Etablissement d’un programme de calcul de l’écoulement

méridien d’une turbomachine radiale par la méthode des éléments finis »
Journal de Mécanique Appliquée - N° 4 - 1979
223 - DENTON J.D. et SINGH U.K. : « Time marching methods for turbomachinery
flows calculation »
VKI - Lectures series - 1979
224 - VEUILLOT J.P. : « Calculation of the quasi-three dimensional flow in a

turbomachine blade row »
ASME - N° 76- GT-56 - 1976
225 - CROCCO L. : « A suggestion for the numerical solution of the steady Navier-
Stokes equations »
AIAA - Vol. 3 - N° 10
226 - KATSANIS T. : « Computer program for calculating velocities and streamlines

on a blade-to-blade surface of a turbomachine »
NASA - TN 4525 - 1968
NASA - TN 8430 - 1977

227 - WEISS P. : « Etude théorique et expérimentale de la compression diphasique.

Application au pompage de liquides aérés »
1996 - Thèse de Doctorat - LEMFI - ENSAM - Paris
228 - BELAMRI T. : « Modélisation des écoulements par la méthode des singulatités »

1994 - DEA Conversion de l’Energie - Filière I - LEMFI - ENSAM - Paris
229 - ROUIBI M. : « Mise en oeuvre numérique et analyse critique de différents

modèles de couche limite »
1995 - DEA Conversion de l’Energie - Filière I - LEMFI - ENSAM - Paris
230 - KUS U. : « Contribution à la prédiction des performances globales et locales

des grilles d’aubes »
231 - HERRIG L., EMERY C., ERWIN R. : « Systematic two dimensionnal cascade
tests of NACA 65 - Series compressor blades at low speeds »
1957 - NACA TN 3916
232 - REY R. : « Methode générale de détermination d’un étage de turbomachine

axiale de compression »
1981 - Thèse de Doctorat d’Etat - Paris VI
233 - NOGUERA R. : « Contribution à la maîtrise du dimensionnement et du

fonctionnement des turbomachines axiales - Etude des débits partiels et de la
cavitation »
1987 - Thèse de Doctorat d’Etat - Paris VI
234 - BAKIR F. : « Méthode de dimensionnement et d’analyse des machines de

compression hélico-centrifuges en régime incompressible »
235 - KOUIDRI S. : « Etude théorique et expérimentale de l’influence des paramètres

constructifs sur les caractéristiques aéroacoustiques des ventilateurs axiaux »


Cinematique Et Dynamique Des Fluides Uee PDF

Transféré par

Droits d'auteur :

Formats disponibles

Cinematique Et Dynamique Des Fluides Uee PDF

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Cinematique Et Dynamique Des Fluides Uee PDF

Transféré par

Droits d'auteur :

Formats disponibles

UEE PA6 – IFMAT et MASTER IMCE

CINEMATIQUE ET DYNAMIQUE DES FLUIDES

Septembre 2008 R. REY – Professeur

CINEMATIQUE ET DYNAMIQUE DES FLUIDES

CINEMATIQUE ET DYNAMIQUE DES FLUIDES

1.1. Définitions - Propriétés des lignes de courant

1.2. Accélération locale - Rotationnel

1.3. Equation de continuité

1.4. Ecoulement irrotationnel - Fonction potentiel

1.5. Circulation du vecteur vitesse

II) ECOULEMENTS BIDIMENSIONNELS p. 15

2.1. Ecoulements plans en coordonnées cartésiennes

2.2. Ecoulements bidimensionnels de révolution

2.3. Méthode simplifiée de tracé des lignes de courant

2.4. Solutions analytiques de l’équation de Laplace

2.5. Application de la superposition des écoulements

CINEMATIQUE ET DYNAMIQUE DES FLUIDES 2

- PROPRIETES CINEMATIQUES DES ECOULEMENTS -

Pour le concepteur, les problèmes posés demeurent inchangés et se formulent en termes

L'application a posteriori des équations dynamiques permettra d'évaluer les pressions

CINEMATIQUE ET DYNAMIQUE DES FLUIDES 3

La définition a priori du champ de vitesses est souvent inspirée des considérations

Pour simplifier le traitement du problème direct et du problème inverse, l'écoulement

CINEMATIQUE ET DYNAMIQUE DES FLUIDES 4

1.1. Définitions - Propriétés des lignes de courant

Dans le domaine des turbomachines, la plupart des écoulements et leurs équations

Figure 1 : Définition du tube de courant

La description lagrangienne de l'écoulement consiste à observer, en fonction du temps, le

Les composantes du vecteur vitesse dépendent à la fois des coordonnées du point M et du

avec : dS : élément de surface

On définit également le débit volume comme le flux à travers la surface S.

CINEMATIQUE ET DYNAMIQUE DES FLUIDES 6

ρ1 et ρ2 étant les masses volumiques moyennes dans les sections S1 et S2 (calculables à

1.2. Accélération locale - Rotationnel

et ses trois composantes en coordonnées cartésiennes :

Sachant que dans le cas général u, v, w sont des fonctions de x, y, z et t on obtiendra la

CINEMATIQUE ET DYNAMIQUE DES FLUIDES 7

On pourra retenir d'une façon générale la relation vectorielle (forme de Lamb) :

CINEMATIQUE ET DYNAMIQUE DES FLUIDES 8

1.3. - Equation de continuité

Cette équation exprime la conservation de la masse contenue dans un volume élémentaire

En utilisant cette seconde configuration, on aura (figure 2) :

Masse totale sortant au même instant :

Figure 2 : Flux à la surface de l'élément fluide

CINEMATIQUE ET DYNAMIQUE DES FLUIDES 9

La masse sortante est donc égale à :

La diminution de masse à l'intérieur de l'élément pendant la même durée dt:

En égalant les deux valeurs de dm on obtient l'équation de continuité :

Ou encore sous forme plus compacte (annexe A4):

En coordonnées cylindriques, l'équation de continuité s'écrira :

- Ecoulement d'un fluide incompressible:

CINEMATIQUE ET DYNAMIQUE DES FLUIDES 10

Dans ce cas, la relation (17) se ramène à :

1.4. Ecoulement irrotationnel - Fonction potentiel

On l'écrit également (annexe R3), en coordonnées cartésiennes :

CINEMATIQUE ET DYNAMIQUE DES FLUIDES 11

Propriétés de la fonction potentiel φ :

Propriétés de l'écoulement irrotationnel :

D'après la relation (15), l'écoulement irrotationnel n'est le siège d'aucun vecteur