CH 8
CH 8
CH 8
VIII.1 Introduction
Nous venons d’étudier dans le chapitre précédent l’écoulement visqueux dans les
conduits et canaux, i.e. l’écoulement interne. Ces écoulements ont des caractéristi-
ques de l’écoulement couche limite suivant :
L’écoulement visqueux de couche limite dans les écoulements externes est par
contre complètement différent comme nous pouvons voir pour l’écoulement sur
une plaque plane à la figure 8.1 :
VIII-1
Le paramètre le plus important dans un écoulement externe est le nombre de Re-
ynolds défini comme :
ρU ∞ x U ∞ x
Re = =
µ ν
où ρ = densité de fluide, µ = viscosité dynamique de fluide, ν = viscosité cinémati-
que de fluide, U ∞ = vitesse caractéristique de l’écoulement, x = dimension caracté-
ristique de l’écoulement.
500,000 <Recr<3,000,000
Pour des écoulements incompressibles sur une plaque plane le nombre critique cor-
respondant à xcr est :
ρU ∞ xcr
Re = = 500,000
µ
Donc, pour x<xcr l’écoulement est laminaire et pour x>xcr est turbulent.
Écoulement Laminaire
La solution de problèmes de couche limite est obtenue des équations de couche li-
mite, dite les équations de Prandtl. Elles sont obtenues des équations de Navier-
Stokes, en appliquant les hypothèses et approximations de couche limite :
VIII-2
i.e. la vitesse dans la direction y est négligeable par rapport à celle dans la direction
x et la variation dans la direction x est plus faible que celle dans la direction y.
Donc, les composants en x et y des équations de Navier-Stokes donnent :
⎛ ∂u ∂u ⎞ ∂p ⎛ ∂ 2u ⎞
ρ ⎜⎜ u + v ⎟⎟ = − + µ ⎜⎜ 2 ⎟⎟
⎝ ∂x ∂y ⎠ ∂x ⎝ ∂y ⎠
∂p
≈ 0 i.e. p = p(x)
∂y
dp dU dp dU
+ ρU = 0 ou = − ρU
dx dx dx dx
i.e. le gradient de pression est donné et le profil de vitesse U = U(x) est connu de
l’analyse de l’écoulement potentiel.
∂u ∂v
+ =0
∂x ∂y
⎛ ∂u ∂u ⎞ dU ⎛ ∂ 2u ⎞
ρ ⎜⎜ u + v ⎟⎟ = U + ν ⎜⎜ 2 ⎟⎟
⎝ ∂x ∂y ⎠ dx ⎝ ∂y ⎠
η = 0 : f = 0, f ′(0) = 0
η = ∞: f ′ =1
VIII-3
La solution est obtenue par intégration numérique. Le profile de vitesse u/U = f’(η)
ainsi obtenu concorde très bien avec des résultats expérimentaux. La contrainte de
cisaillement τw sur la paroi est donnée par
∂u ⎡ ∂η ⎤
1/ 2
⎛U ⎞
τ w = µ | y =0 = ⎢ µUf '' ⎥ = µUf '' (0)⎜ ⎟
∂y ⎣ ∂y ⎦ y =0 ⎝ νx ⎠
Ce résultat peut être exprimé à l’aide d’un coefficient de frottement cf défini selon :
τw
cf =
1
ρU 2
2
On obtient :
0.664
cf =
Re x
La force de traînée par unité de largeur de la plaque et uniquement pour une face
ainsi que le coefficient de traînée sont :
x
Fx 1.328
Fx = ∫ τ w ( x)dx et C x = =
1 Re x
0 ρU 2 x
2
VIII-4
Point Important : Donc l’écoulement est laminaire le long de la plaque et nous
pouvons utiliser les équations pour l’écoulement laminaire à une position x de 0 à 1
m.
τ w = 0.0030 N / m2 [ Pa ]
FD = 0.0259 N
VIII-5
Points importants :
1. Chacun de ces calculs peut être fait à une position x soit x ≤ L .
Problème 3 : L’air fait un mouvement sur une plaque plane de L = 60 cm telle que
ReL = 2E5. Tracez l’épaisseur de couche limite pour 0 < x <L.
Écoulement Turbulent
Dans le cas de l’écoulement turbulent c’est difficile de faire une analyse analytique
semblable à celle pour laminaire en utilisant les équations de mouvement. Pourtant,
nous avons la possibilité de faire une analyse en utilisant la méthode intégrale et
d’obtenir des équations adéquates et pratiques en utilisant des résultats expérimen-
taux et les paramètres pertinents.
VIII-6
0.031
CD =
Re1/7
L
Ces relations sont valables pour une plaque plane et lisse. Nous avons des corréla-
tions pour les plaques avec surface rugueuse. Comme dans le cas de l’écoulement
interne, la rugosité est exprimée comme le rapport de la rugosité absolue à la lon-
gueur caractéristique. Dans ce cas on utilisera l’inverse de ce rapport comme : x/ε
ou L/ε. Dans le régime complètement rugueux, CD est indépendant du nombre de
Reynolds. En conséquence, le frottement varie directement avec U2. Les relations
empiriques pour le coefficient local et le coefficient de frottement sont :
−2.5
⎛ x⎞
c f = ⎜ 2.87 + 1.58 log ⎟
⎝ ε⎠
− 2 .5
⎛ L⎞
C D = ⎜1.89 + 1.62 log ⎟
⎝ ε⎠
La figure 8.2 montre la relation CD = f(ReL, L/ε) où les relations sont tracées pour
l’écoulement laminaire, pour l’écoulement turbulent sur une plaque plane lisse et
l’écoulement turbulent dans le régime complètement rugueux. Fig. 8.2 montres
aussi la région de transition qui est 500x103<ReL<80x106. Notons que la transition
de la laminaire au turbulent commence à 500x103 comme limite inférieure et à
3x106 comme limite supérieure. Schlichting a proposé les corrélations suivantes
pour la région de transition avec ces deux limites :
0.031 1440
Re trans = 500 x10 3 : C D = −
Re1L/ 7 Re L
0.031 8700
Re trans = 3 x10 6 : C D = −
Re1L/ 7 Re L
Ces deux équations peuvent être utilisées pour calculer le coefficient de frottement
dans le cas où l’écoulement serait laminaire entre 0<x<xcr et devient turbulent pour
x>xcr.
VIII-7
Fig. 8.2 Coefficient de frottement pour les plaques lisses et rugueuses
L’écoulement de couche limite sur une plaque plane consiste en un écoulement la-
minaire pour x<xcr et un écoulement turbulent pour x>xcr, comme montré à la fi-
gure 8.1. Donc, normalement nous devons calculer l’épaisseur de couche limite et
le cisaillement en utilisant les expressions appropriées pour ces deux modes
d’écoulement :
5x 1.328
δ ( x) = et C D = pour x < xcr
Re x Re L
0.16 x et 0.031
δ (x ) =
pour x>xcr
1/ 7 C D = 1/7
Re x Re L
Si la longueur de la plaque est très longue ou bien un fil de déclenchement est utili-
sé pour faire l’écoulement turbulent à partir du bord d’attaque, nous pouvons cal-
culer le frottement moyen en utilisant seulement le coefficient de frottement pour
l’écoulement turbulent :
0.031
CD =
Re1/7
L
VIII-8
Il faut noter que si nous calculons de cette façon, le frottement total moyen sera lé-
gèrement surestimé puisque le frottement de l’écoulement turbulent est plus grand
que celui laminaire.
Exemple 2 : L’eau s’écoule sur une plaque plane lisse d’une dimension de L =
2.55 m et w = 1 m. U = 2 m/s, ρ = 1000 kg/m3, ν = 1.02E-6 m2/s. Estimez l’erreur
en FD si l’on suppose que l’écoulement est turbulent sur la plaque entière.
U∞L (2 m / s )(2.55 m )
Re L = = = 5E 6 > 500, 000
υ 1.02 E − 6 m 2 / s
Pour Recr = 500x103 on trouve xcr = 0.255 m; i.e. sur 10% de la plaque
l’écoulement est laminaire.
0.031 0.031
CD ≈ 1/7 = = 0.003423
Re L (
5x106 1/ 7
)
1 1 kg 2 m 2
FD = ρ U ∞ C D A = ⋅ 1000 3 ⋅ 2 2 ⋅ 0.003423 ⋅ 2.55m 2
2
2 2 m s
FD =17.46 N
La valeur trouvée doit être élevée puisqu’on a supposé que la plaque entière soit
turbulente tandis que jusqu’à xcr = 0.255 m l’écoulement est laminaire.
1 1 kg 2 m 2
FD = ρ U ∞ C D A = ⋅ 1000 3 ⋅ 2 2 ⋅ 0.003135 ⋅ 2.55m 2
2
2 2 m s
VIII-9
FD =15.99 N
17.46 − 15.99
Erreur = ⋅ 100 = 9.2% plus élevée
15.99
Donc, la traînée est surestimée par 9.2% quand on suppose l’écoulement est turbu-
lent sur la plaque entière.
(c) Calculons maintenant séparément pour les parties de couche limite laminaire et
turbulente :
La partie laminaire :
1.328
CD = = 0.001878
Re L
FD = (0.5)(1000)(22)(0.255x1)(0.001878) = 0.96 N
La partie turbulente :
0.031
CD = = 0.00342
Re1/7
L
2
FD = (0.5)(1000)(2 )((2.55 - 0.255)x1)(0.00342) = 15.70 N
Problème 5 : Un chaland avec une carène plate est tiré avec une vitesse de 1.543
m/s. Il a une longueur de 40 m et largeur de 14 m. On suppose que les côtés du
chaland ont un frottement négligeable. Calculez la puissance nécessaire en kW
pour le tirer : (a) pour la surface lisse, (b) pour la surface rugueuse avec une rugosi-
té absolue de 3 mm. Les données : ρ = 1025 kg/m3, µ = 1.07E-3 kg/m.s. Rép. : (a)
2.535 kW, (b) 4.911 kW
VIII-10
Méthode Intégrale
L’épaisseur de Déplacement δ*
VIII-11
Fig. 8.3 Effet de d/placement dans la couche limite
Nous pouvons voir que l’épaisseur de déplacement varie seulement avec le profil
de vitesse sans dimension u/ U∞ et elle est δ*= f(δ).
Exemple 3 :
u ⎛ y ⎞ ⎛ y⎞ 2
=2 −
U∞ ⎝δ ⎠ ⎝δ⎠
Solution : Nous notons que ce profil quadratique satisfait les conditions aux fron-
tières, notamment :
(a) pour y = 0 : u = 0
(b) pour y = δ : u = U∞
VIII-12
Nous pouvons voir qu’en utilisant le profil quadratique, nous obtenons une épais-
seur de déplacement égale à un tiers de l’épaisseur de couche limite. Cette valeur
donne une approximation satisfaisant en comparaison de la valeur obtenue par la
méthode analytique de 0.344 δ.
Point important : en choisissant le profil de vitesse de couche limite, faites sûre
que les conditions aux frontières sont satisfaites et le profil u(y) soit rapproché as-
sez bien.
δ u ⎛ u⎞
θ=∫ ⎜1 − ⎟dy
0U∞ ⎝ U∞ ⎠
VIII-13
D/A
CD = 1
2
ρ U 2∞
Nous pouvons montrer que :
2 θ (L )
CD =
L
Donc nous voyons qu’en connaissant le profil de vitesse de couche limite u = f(y),
nous pouvons déterminer aussi le coefficient de frottement.
τw
ρ
= δ * U∞
d U∞
dx
+
d
dx
(U∞2 θ )
d U∞ δ* d P
δ U∞
*
=−
dx ρ dx
dθ τ 1 dU ∞ Cf dU ∞
= w2 − (2θ + δ * ) = − ( H + 2)
dx ρU ∞ U ∞ dx 2 dx
où H = δ∗/θ définit le paramètre de forme qui joue un rôle important dans la théorie.
Pour un écoulement non accéléré sur une plaque plane, on peut montrer que
d U∞
P = const., U∞ = const., =0
dx
L’équation intégrale de von Karman devient :
VIII-14
τw d dθ
=
ρ dx
(U∞2 θ )= U∞2
dx
Si l’on suppose le profile de vitesse de Blasius pour la couche limite :
u ⎛ y ⎞ ⎛ y⎞ 2
=2 − = 2η − η 2
U∞ ⎝δ ⎠ ⎝δ⎠
d u⎞ ⎧ 2 2 y⎫ 2 µ U∞
τw = µ ⎟ = 2U∞ ⎨ − 2 ⎬ =
d y⎠ w ⎩δ δ ⎭ y=0 δ
δ u ⎛ u⎞
θ=∫ ⎜1 − ⎟dy
0U∞ ⎝ U∞ ⎠
δ 2δ
θ = ∫ (2η − η 2 )(1 − 2η + η 2 )δ d η =
0 15
dθ 2 dδ
τ w = ρ U∞2 = ρ U∞2
d x 15 dx
2 d δ 2 µ U∞
τw = ρ U∞2 =
15 dx δ
15 µ
δ dδ = dx
ρ U∞
On trouve :
1/ 2
⎧ 30 µ x ⎫ δ 5.5
δ =⎨ ⎬ ou =
⎩ ρ U∞ ⎭ x Re x
VIII-15
Nous voyons que l’épaisseur de couche limite trouvée est satisfaisante par rapport
à δ obtenu en utilisant la méthode analytique par Blasius.
Pour une plaque plane nous devons considérer seulement le frottement de friction
et pas celui de pression; donc, nous pouvons écrire :
τw (x) 2 µU ∞ 1
cf = =
1
2
ρ U ∞2 δ 1
2
ρ U ∞2
En substituant pour δ on obtient :
2 µU ∞ Re 0.73
cf = =
5.48 1
2
ρ U 2
∞ Re x
u y
= 2η − 2η 3 + η 4 avec η =
U∞ δ
Ces forces externes appliquées et les moments sont généralement une fonction de
paramètres suivants :
VIII-16
Les forces et les moments sont généralement exprimés dans une forme adimen-
sionnelle, comme par exemple, la traînée,
FD /A
CD = .
2 ρ U∞
1 2
(a) L’aire frontale, i.e. la projection de corps : Elle est généralement utilisée
pour les corps non aérodynamiques, comme les édifices, véhicules ;
(b) L’aire vue de haut, i.e. la projection sur le plane horizontal : Elle est utilisée
pour les corps minces, e.g. les ailes, hydrofoil ;
(c) L’aire mouillée : L’aire totale en contact avec le fluide. Elle est utilisée pour
les bateaux, péniches.
Autre que la traînée de frottement due à la friction, il y a aussi la traînée de pres-
sion due à la pression ou la forme qui peut être plus important. La traînée de pres-
sion est une conséquence de la distribution de pression sur la surface d’un corps.
Le coefficient combiné de la traînée peut être exprimé comme :
C D = CD,press + C D,friction
ou
VIII-17
La contribution majeure de chacun dépend de la forme aérodynamique de corps et
aussi bien que des conditions de l’écoulement.
Considérons les deux corps montrés à la figure 8.6. L’écoulement autour de profile
aérodynamique reste toujours attaché à la surface, sans séparation de l’écoulement.
En conséquence, la traînée de pression est faible. Par contre, l’écoulement autour
de cylindre, une forme non aérodynamique, sépare de sa surface ainsi créant une
région de sillage à basse pression; comme résultat, la traînée de pression devient
importante.
(b)
Fig. 8.6 (a) corps non aérodynamique, (b) corps aérodynamique sans séparation de
couche limite
Ce phénomène illustré à la figure 8.7. Pour un profil mince, i.e. t/c = 0, toute la
traînée est due à friction. Pour un cylindre, i.e. t/c = 1, la traînée de friction est seu-
lement 3% tandis que la traînée de pression est de 97%.
VIII-18
Fig. 8.7 Effet de rapport d’épaisseur sur la traînée de pression
Dans le cas de corps non aérodynamiques la traînée de pression peut être changée
en réduisant la grandeur de sillage à basse pression comme illustré à la figure 8.8.
Fig. 8.8 Écoulement autour d’un cylindre circulaire (a) séparation laminaire, (b)
séparation turbulente
VIII-19
Fig. 8.9 Distribution de pression théorique et expérimentale dans un écoulement
autour d’un cylindre circulaire
Finalement, l’effet de forme aérodynamique sur la traînée totale est illustré à la fi-
gure 8.10. On put observer que
(1) Comme la forme du corps change d’une forme contondante à une forme plus
aérodynamique l’effet de la traînée de pression diminuera.
(2) La transformation de (a) à (b) et puis de (b) à (c) peut causer une augmenta-
tion de traînée de friction mais puisque la traînée dominante est due à la
pression, le résultat net est une traînée réduite.
Le tableau 8.2 donne les CD pour les corps de deux dimensions à Re > 104. Dans ce
cas, la troisième dimension, perpendiculaire à la page est très longue.
Remarque : les corps non aérodynamiques avec le point de séparation fixe ont un
seul coefficient de CD sans tenir compte de nombre de Reynolds, e.g. carré, cylin-
dre, le demi-tube, etc. Les corps aérodynamiques ont usuellement réduction de CD
pour changement de l’écoulement laminaire à turbulent; ceci est à cause de chan-
gement favorable de point de séparation de couche limite, e.g. cylindres elliptiques.
VIII-20
Tableau 8.2 Coefficient de traînée pour les corps de deux dimensions
Le tableau 8.3 montre les CD pour les corps de trois dimensions à Re > 104. La troi-
sième dimension perpendiculaire à la page n’est pas longue. La même remarque
faite pour les corps de deux dimensions s’applique aussi pour ce cas.
VIII-21
Tableau 8.3 Coefficient de traînée pour les corps de trois dimensions
Exemple 4 :
Un pilotis d’une section carrée est 0.15 mx0.15 m et 6 m longue. Il situe dans l’eau
comme montrée à la figure 8.11. La vitesse de l’eau est de 1.5 m/s. Estimez le mo-
ment maximum appliqué au fond du pilotis. La densité de l’eau est de 1000 kg/m3,
sa viscosité cinématique est de 1.307x10-6 m2/s.
VIII-22
Fig. 8.11 Schéma de l’exemple de pilotis
Solution : Supposons que le pilotis peut être traité comme deux dimensionnel et
l’effet de bout soit négligeable. On calcule le nombre de Reynolds comme
Re = (1.5m/s)(0.15m)/1.307x10-6m2/s = 2.1x105.
Le tableau 8.2 donne CD = 2.1. L’aire frontale est de (6m)(0.15m) = 0.9 m2. On ob-
tient pour la traînée :
VIII-23