Introduction

à l’hydrodynamique

Cours et exercices  Licence 2

S.E. Belhenniche
Préface

Ce Cours est destiné aux étudiants de la deuxième année universitaire filière génie
maritime. Parmi les objectifs du cours est de permettre à l’étudiant d’acquérir les
notions de base de la dynamique des fluides qui seront utiles à résoudre des
problèmes liés à l’hydrodynamique navale au cours de leurs cursus.
En commençant avec des définitions et propriétés des fluides et allant jusqu’à la
démonstration des équations de d’Euler et de Navier stokes instantanées et
moyennées. Ce cours est consacré aussi à l’utilité de l’analyse dimensionnelle, la
tension superficielle et acquérir des notions fondamentales sur la couche limite et
ces grandeurs caractéristiques, ainsi que les écoulements autours des obstacles. Les
equations moyennées sont développées dans le chapitre turbulence ; les méthodes
de modélisations des tenseurs de Reynolds sont mentionnées.
A la fin de chaque chapitre, des exercices sont proposés avec des réponses
permettant de tester les connaissances des étudiants et de se préparer aux examens.

I
Sommaire

I DEFINITIONS ET PROPRIETES DES FLUIDES ………..…. 01

I.1 Introduction .………………………………………………….……….………………. 02
I.2 Analyse dimensionnelle ………………………………………………….………… 02
I.2.1 Equation aux dimensions …………………………………..…………….. 03
I.2.2 Système international (SI) ………………………………………………… 03
I.2.3 Système CGS …………………………………………………………………… 04
I.3 Définition d’un fluide ………………………………………………………………. 05
I.4 Fluide parfait et fluide réel ………………………………………………….…… 06
I.5 Fluide Compressible et incompressible ………………………………….….. 06
I.6 Les caractéristiques d’un fluide …………………………………………………. 06
I.6.1 La viscosité …………………………………..………………………………... 06
I.6.2 La masse volumique ………………………………………………………… 07
I.7 Calcul du couple des disques en fonction de la contrainte de
cisaillement ……………………………………………………………………………. 08
I.8 Ecoulements laminaires et turbulents …………………………..…………… 09
I.9 Ecoulements internes et externes ………………………………….………….. 10
I.9.1 Ecoulements internes ………………………………….………………….. 10
I.9.2 Ecoulements externes ……………………………………………………… 11
I.10 Exemples évolution de la forme géométrique (aérodynamique) ………. 11
I.10.1 Forme géométrique d’une voiture …………………………………….. 11
I.10.2 Forme géométrique d’une balle ………………………………………… 12
I.11 Etude des écoulements externes ……………………………………………….. 13
I.11.1 Etude expérimentale et visuel : Moyens expérimentaux ………. 13
I.11.2 Etude théorique ……………………………………………………………….. 13
I.12 Exercices avec réponses …………………………………………………………... 14
II TENSION SUPERFICIELLE …………………………………………………. 16
II.1 Introduction .………………………………………………….………………….……. 17
II.2 Forces de cohésion et d’adhésion ……………………………………………… 17
II.3 La loi de Young Laplace …………………………………………………………… 20
II.4 Angles de contact et capillarité ……………………………………..........…… 21
II.5 Phénomène de mouillabilité ……………………………………………..……… 23
II.6 Exercices avec réponses …………………………………………………….……… 25
III LES EQUATIONS DE CONSERVATION …………………………. 26
III.1 Conservation de la masse (Equation de continuité) ……………………. 27
III.1.1 Le Théorème de Transport de Reynolds ……………………………. 27
III.1.2 Forme intégrale de l’équation de continuité ……………………….. 27
III.1.3 Forme différentielle de l’équation de continuité …………………. 28
III.1.4 Forme vectorielle de l’équation de continuité ……………………… 30

II
 III.2 Conservation de la quantité de mouvement (Equations d’Euler) …… 31
III.2.1 La 2 ème loi de Newton ……………………………………………………. 32
III.2.2 Forme vectorielle de l’équation d’Euler ……………………………… 35
III.3 Conservation de la quantité de mouvement (Equations de Navier
Stokes) ……………………………………………………………………………………. 36
III.3.1 Le tenseur de contraintes total …………………………………………. 36
III.3.2 Les contraintes normales et les contraintes tangentielles …….. 36
III.3.3 Forme vectorielle de l’équation de Navier Stokes ………………… 38
III.4 Exercices avec réponses …………………………………………………………… 40
IV LA COUCHE LIMITE ……………………………………………………………. 42
IV.1 Introduction .………………………………………………….………………….……. 43
IV.2 Grandeurs caractéristiques de la couche limite ………………………….. 44
IV.2.1 Epaisseur ………………………………………………………………………… 44
IV.2.2 Epaisseur de déplacement ………………………………………………… 45
IV.2.3 Epaisseur de la quantité de mouvement …………………………….. 46
IV.2.4 Autres facteurs et coefficients …………………………………………… 46
IV.3 Evolution de la couche limite ……………………………………….…………… 47
IV.3.1 Plaque plane …………………………………………………………………… 47
IV.3.2 Une forme arrondie (profil, cylindre….) …………………………….. 48
IV.4 Formules utilisées pour le calcul des paramètres de
la la couche limite (cas d’une plaque plane) ………………………………….. 50
IV.5 Exercices avec réponses …………………………………………………………… 51
V ECOULEMENTS AUTOUR DES OBSTACLES ………………. 53
V.1 Introduction .………………………………………………….………………….……. 54
V.2 Evolution expérimentale de 𝐶𝑥 (𝑅𝑒 ) et description de l’écoulement
autour d’une sphère lisse ………………………………………………………….. 56
V.2.1 Evolution expérimentale de 𝐶𝑥 (𝑅𝑒 ) ……………………………………. 56
V.2.2 Description de l’écoulement ……………………………………………… 57
V.3 Evolution expérimentale de 𝐶𝑥 (𝑅𝑒 ) et description de l’écoulement
autour d’un cylindre ………………………………………..……………………….. 59
V.3.1 Evolution expérimentale de 𝐶𝑥 (𝑅𝑒 ) …………………………………… 59
V.3.2 Description de l’écoulement ……………………………………………… 60
V.4 Ecoulement autour d’un hydrofoil ………………..…………………………… 61
V.5 Ecoulement autour d’une hélice ………………………………………………... 62
V.6 Ecoulement autour d’une plaque plane (calcul du 𝐶𝑥 ) ………………….. 63
V.7 Quelques paramètres influents sur le coefficient de
Trainée 𝐶𝑥 ……………………………………………………………………………….. 65
V.7.1 Position de la plaque ………………………………………………………... 65
V.7.2 La surface …………………………………………………………………....... 65

III
 V.8 Optimisation des formes géométriques ……………………………………… 65
V.9 Exercices avec réponses ……………………………………………………………. 67
VI TURBULENCE ……………………………………………………………………….. 69
VI.1 Introduction .………………………………………………….………………….……. 70
VI.2 Les equations moyennées …………………………………………………………. 70
VI.2.1 Principe ………………………………………………………………………….. 70
VI.3 Rappel des equations instantanées ……………………………………………. 71
VI.3.1 Equations principales ………………………………………………………. 71
VI.4 Critère de choix et classification ………………………………………………… 74
ANNEXES …………..………………………………………………………….……….. 75
REFERENCES …………..………………………………………….……………….. 78

IV
Liste des figures

Fig (I.1) Branches de la mécanique des fluides ………………………………………… 02

Fig (I.2) Etapes à suivre pour le calcul de la puissance du moteur à bord ……. 04
Fig (I.3) Profil de vitesse (écoulement autour d’une plaque plane) …………….. 06
Fig (I.4) Cylindre mobile autour d’un cylindre fixe …………………………………… 08
Fig (I.5) (a) écoulement laminaire (b) écoulement turbulent …………………….. 09
Fig (I.6) Ecoulement dans une conduite (écoulement interne) ………………….. 10
Fig (I.7) Expérience de Reynolds (différents régimes d’écoulements) ………… 11
Fig (I.8) Ecoulement autour d’un cylindre (écoulement externe) ………………. 11
Fig (I.9) Evolution des formes géométriques des voitures …………………………. 12
Fig (I.10) (a) écoulement autour d’une balle lisse (b) écoulement autour
d’une balle de golf …………………………………………………………………….. 12
Fig (II.1) Efforts intermoléculaires des différentes molécules dans un liquide
avec une surface libre ……………………………………………………………….. 17
Fig (II.2) Exemples de la tension superficielle …………………………………………… 18
Fig (II.3) Experience d’un film de savon entre 02 tiges ………………………………. 18
Fig (II.4) L’expérience de l’anneau …………………………………………………………… 20
Fig (II.5) Bulle sphérique d’eau ……………………………………………………………….. 20
Fig (II.6) Bulle d’eau savonneuse …………………………………………………………….. 21
Fig (II.7) Formation d’un ménisque dans un tube de verre plongé dans un
liquide mouillant ……………………………………………………………………… 22
Fig (II.8) Ascension et dépression capillaires ……………………………………………. 23
Fig (II.9) Mouillabilté d’une goutte d’eau sur différents matériaux ……………… 24
Fig (II.10) Goutte d’eau sur une feuille de cactus ………………………………………… 24
Fig (III.1) Vecteur normal et vitesse sur une surface élémentaire 𝑑𝐴 …………… 27
Fig (III.2) Débit entrant et sortant pour un petit élément 𝑑𝑣 ……………………….. 29
Fig (III.3) Vecteur normal et vitesse sur une surface élémentaire 𝑑𝐴 …………… 31
Fig (III.4) La Contrainte normale et les contraintes tangentielles sur une
surface élémentaire 𝑑𝐴 ……………………………………………………………… 33
Fig (III.5) La Contrainte normale et les contraintes tangentielles sur le plan X 34
Fig (IV.1) Schéma descriptif de la couche limite ………………………………………… 43
Fig (IV.2) Formation d’une couche limite sur une plaque plane …………………… 43
Fig (IV.3) Décollement de la couche limite ………………………………………………… 44
Fig (IV.4) Epaisseur de la couche limite …………………………………………………….. 44
Fig (IV.5) Epaisseur de déplacement ………………………………………………………… 45
Fig (IV.6) Epaisseur de quantité de mouvement ………………………………………… 46
Fig (IV.7) Structure de l’écoulement dans la couche limite ………………………….. 47
V
Fig (IV.8) Evolution du profil de vitesse ……………………………………………………. 48
Fig (IV.9) Décollement de la couche limite (écoulement autour d’un cylindre) 49
Fig (V.1) Les forces de pression et les forces de frottement exercées sur
un profil…………………………………………………………………………………… 54
Fig (V.2) La force de trainée et la force de portance exercées sur un profil …… 55
Fig (V.3) Influence de la forme géométrique sur le coefficient de trainée …….. 55
Fig (V.4) L’évolution du coefficient de trainée en fonction du nombre de
Reynolds 𝑅𝑒 (cas d’une sphère lisse) ………………………………………….. 56
Fig (V.5) Description de l’écoulement en fonction du nombre de Reynolds 𝑅𝑒
(cas d’une sphère lisse) …………………………………………………………….. 59
Fig (V.6) L’évolution du coefficient de trainée en fonction du nombre de
Reynolds 𝑅𝑒 (cas d’un cylindre) …………………………………………………. 59
Fig (V.7) Description de l’écoulement en fonction du nombre de Reynolds 𝑅𝑒
(cas d’un cylindre) ……………………………………………………………………. 61
Fig (V.8) Coefficient de portance 𝐶𝑧 en fonction de l’angle d’incidence
(cas d’un hydrofoil) ………………………………………………………………….. 62
Fig (V.9) Cas de l’écoulement autour d’une hélice marine ………………………….. 63
Fig (V.10) Coefficient de trainée 𝐶𝑥 en fonction de nombre de Reynolds 𝑅𝑒
(cas d’une plaque plane à différents coefficient de rugosité) …………. 64
Fig (V.11) Ecoulement autour d’une plaque plane. 65
Fig (V.12) Choix de la surface pour le calcul du Coefficient de trainée 𝐶𝑥
(différentes formes géométriques) …………………………………………….. 65
Fig (V.13) Coefficnet de trainée 𝐶𝑥 en fonction du nombre de Froud 𝐹𝑟
(une coque avec et sans un bulbe) ……………………………………………… 66
Fig (VI.1) Champs de vitesse en un point dans un régime turbulent …………….. 70
Fig (VI.2) Structure de l’écoulement laminaire et turbulent (cas d’une fumée
de cigarette) …………………………………………………………………………….. 70

VI
Liste des tableaux

Tab (I.1) Unités fondamentales du système international (SI) ……………….. 03

Tab (I.2) Dimensions et unités en système international (SI) et anglo-
saxon …………………………………………………………………………………… 04
Tab (II.1) Angle de contact et capilarité entre 02 matériaux différents ……… 23
Tab (IV.1) Les formules de Blasius et Prandtl pour le calcul des paramètres
de la couche limite (cas d’une plaque plane) ……………………………. 50
Tab (V.1) Valeurs numériques du coefficient de trainée en fonction du
nombre de Reynolds 𝑅𝑒 (cas d’une sphère lisse) ……..……………….. 56
Tab (V.2) Description des zones caractéristiques de la courbe du
coefficient de trainée (cas d’une sphère lisse) ………………………….. 56
Tab (V.3) Valeurs numériques du coefficient de trainée en fonction du
nombre de Reynolds 𝑅𝑒 (cas d’un cylindre) …………………………….. 60
Tab (VI.1) Inconnus à calculer pour le cas des équations instantanées et
moyennées …………………………………………………………………………… 74

VII
Liste des symboles

SYMBOLES LATINS

𝐶𝐿 , 𝐶𝑧 Coefficient de portance [−]

𝐶𝐷 , 𝐶𝑥 Coefficient de trainée [−]
𝐶𝑓 Coefficient de frottement local [−]
𝐶𝐹 Coefficient de frottement global [−]
𝑑𝑇 Poussée élémentaire [𝑁]
𝑑𝑄 Couple élémentaire [𝑁. 𝑚]
𝑑, 𝐷 Diamètre [𝑚]
𝐷⃗ , 𝐹𝑥 La force de trainée [𝑁]
𝐸 Energie [𝑁. 𝑚] ou Joule
𝐹 Force [𝑘𝑔. 𝑚⁄𝑠 2 ]
𝐹𝑎 Force d’adhésion [𝑘𝑔. 𝑚⁄𝑠 2 ]
𝐹𝑐 Force de cohésion [𝑘𝑔. 𝑚⁄𝑠 2 ]
𝑓 Fréquence [𝑠 −1 ]
𝐹𝑟 Nombre de Froude [−]
𝑔 Accélération de pesanteur [𝑚⁄𝑠 2 ]
ℎ Hauteur de capilarité 𝑚
𝐻 Facteur de forme [−]
𝐼 Intensité du courant électrique [𝐴]
𝑖 Angle d’incidence [°]
𝐿 Longueur [𝑚]
⃗
𝐿 La force de portance [𝑁]
𝑚 Masse [𝑘𝑔]
Pression 2
[𝑁⁄𝑚 ] ou Pascal
𝑃
𝑃𝑖 Pression intérieur [𝑁⁄𝑚2 ] ou Pascal
𝑃0 , 𝑃𝑜𝑢𝑡 Pression extérieur [𝑁⁄𝑚2 ] ou Pascal
𝑃̅ Pression moyenne [𝑁⁄𝑚2 ] ou Pascal
𝑃′ Fluctuation de la pression [𝑁⁄𝑚2 ] ou Pascal
𝑃𝐸 Puissance de remorquage [𝑁]
𝑞𝑚 Débit massique [𝑘𝑔/𝑠]
𝑞𝑣 Débit volumique [𝑚3 /𝑠]
𝑟, 𝑅 Rayon [𝑚]
𝑅𝑒 Nombre de Reynolds [−]
𝑅𝑒𝑥 Nombre de Reynolds local [−]
𝑅𝑒𝑐 Nombre de Reynolds critique [−]
𝑅𝑇 Résistance totale [𝑁]
𝑆 Surface [𝑚2 ]
𝑇 Couple [𝑁. 𝑚]
𝑈, 𝑉 ⃗ Vitesse [𝑚⁄𝑠]
𝑢 Vitesse suivant X [𝑚⁄𝑠]
𝑢̅ Vitesse moyenne [𝑚⁄𝑠]
𝑢́ Fluctuation de la vitesse [𝑚⁄𝑠]

VIII
𝑣 Vitesse suivant Y [𝑚⁄𝑠]
𝑤 Vitesse suivant Z [𝑚⁄𝑠]
𝑊 Travail [𝑁. 𝑚 ] ou Joule

SYMBOLES GRECS

𝜏 Contrainte de cisaillement [𝑁⁄𝑚2 ] ou Pascal

𝜎 Tension superficielle [𝑁/𝑚]
𝜎𝑠𝑙 Tension superficielle solide/liquide [𝑁/𝑚]
𝜎𝑙𝑣 Tension superficielle liquide/vapeur [𝑁/𝑚]
𝜎𝑠𝑣 Tension superficielle solide/vapeur [𝑁/𝑚]
𝜌 Densité [𝑘𝑔⁄𝑚3 ]
𝛾 Poids spécifique [𝑁⁄𝑚3 ]
𝜇 Viscosité dynamique [𝐾𝑔/𝑚. 𝑠]
𝛿𝑖𝑗 Symbole de Kronecker [−]
𝛿 Epaisseur de la couche limite [𝑚]
𝛿∗ Epaisseur de déplacement [𝑚]
𝜃 Epaisseur de la quantité de mouvement [𝑚]
∇ Volume [𝑚3 ]
𝛺 Vitesse angulaire [𝑠 −1 ]
𝜀 Rugosité [𝑚]

IX
X
Chapitre I

DEFINITIONS ET PROPRIETES DES

FLUIDES

I DEFINITIONS ET PROPRIETES DES FLUIDES ………………………... 01

I.1 Introduction .……………………………….…………………….……………………. 02
I.2 Analyse dimensionnelle …………………….……………………………………… 02
I.2.1 Equation aux dimensions …………………………………..…………….. 02
I.2.2 Système international (SI) ………………………………………………… 03
I.2.3 Système CGS …………………………………………………………………… 04
I.3 Définition d’un fluide ………………………………………………………………. 05
I.4 Fluide parfait et fluide réel …………………………………………….………… 06
I.5 Fluide Compressible et incompressible ……………………………….…….. 06
I.6 Les caractéristiques d’un fluide …………………………………………………. 06
I.2.6 La viscosité …………………………………..………………………………... 06
I.2.6 La masse volumique ………………………………………………………… 07
I.7 Calcul du couple des disques en fonction de la contrainte de
cisaillement ……………………………………………………………………………. 08
I.8 Ecoulements laminaires et turbulents ……………………………..………… 09
I.9 Ecoulements internes et externes …………………………………….……….. 10
I.9.1 Ecoulements internes ………………………………….………………….. 10
I.9.2 Ecoulements externes ……………………………………………………… 11
I.10 Exemples évolution de la forme géométrique (aérodynamique) …… 11
I.10.1 Forme géométrique d’une voiture …………………………………….. 11
I.10.2 Forme géométrique d’une balle ………………………………………… 12
I.11 Etude des écoulements externes ……………………………………………….. 13
I.11.1 Etude expérimentale et visuel : Moyens expérimentaux ………. 13
I.11.2 Etude théorique ……………………………………………………………….. 13
I.12 Exercices avec réponses …………………………………………………………... 14

1
 Définitions et propriétés des fluides

I.1 Introduction :
La mécanique des fluides est la science des lois de I ‘écoulement des fluides. Elle est
la base du dimensionnement des conduites de fluides et des mécanismes de transfert
des fluides. C’est une branche de la physique qui étudie les écoulements de fluides
c'est-à-dire des liquides et des gaz lorsque ceux-ci subissent des forces ou des
contraintes. Elle comprend deux grandes sous branches:

 La statique des fluides, ou hydrostatique qui étudie les fluides au repos. C'est
historiquement le début de la mécanique des fluides, avec la poussée d'Archimède
et l'étude de la pression.

 La dynamique des fluides qui étudie les fluides en mouvement. Comme autres
branches de la mécanique des fluides.

On distingue également d’autres branches liées à la mécanique des fluides :

L’hydraulique, l'hydrodynamique, l'aérodynamique, …Une nouvelle approche a vu le
jour depuis quelques décennies: la mécanique des fluides numérique (CFD ou
Computational Fluid Dynamics en anglais), qui simule l'écoulement des fluides en
résolvant les équations qui les régissent à l'aide d'ordinateurs.

Fluide Mécanique des fluides

(Liquide, Gaz) (Hydro, Aéro)

Compressibilité
Statique des fluides Dynamique des fluides
Densité

Viscosité

Pression Fluide au repos Fluide en mouvement

Tension Statique des fluides ∑ 𝑭𝒊 = 𝟎 Dynamique des fluides

superficielle ∑ 𝑭𝒊 > 𝟎

Fig (I.1) Branches de la mécanique des fluides.

I.2 Analyse dimensionnelle :

L'utilisation des dimensions permet de savoir si une expression littérale est
homogène ou non. Cela permet de rechercher d’éventuelles erreurs. Mais l'analyse
des dimensions ou analyse dimensionnelle permet de retrouver, de deviner des lois
physiques lorsque la résolution théorique est trop complexe.

2
 Définitions et propriétés des fluides

I.2.1 Equation aux dimensions :

On appelle équation aux dimensions une équation reliant la dimension d'une
grandeur G a celles des grandeurs de base.

[𝐺] = 𝑇 𝛼 . 𝐿𝛽 . 𝑀𝛾 . 𝐼 𝛿 . Θ𝜖 . 𝐽𝜁 . 𝑁 𝜂

Ces grandeurs sont respectivement : temps, longueur est le mètre, masse, intensité de
courant, température, intensité lumineuse, quantité de matière.
Exemple : Conversion d’unité :
Pression : 𝑃 = 𝐹 ⁄𝑆 = 𝑀. 𝐿. 𝑇 −2 . 𝐿−2 = 𝑀. 𝐿−1 . 𝑇 −2
Pour les unités mécaniques le choix le plus courant est de prendre la longueur, la
masse et le temps. Il existe 2 systèmes d’unités :

I.2.2 Système international (SI) :

Le SI est fondé sur un choix de sept unités de base bien définies et considérées par
convention comme indépendantes du point de vue dimensionnel : le mètre, le
kilogramme, la seconde, l'ampère, le kelvin, la mole et la candela. Les unités dérivées
sont formées en combinant les unités de base d'après les relations algébriques qui
lient les grandeurs correspondantes

Tab (I.1) Unités fondamentales du système international (SI).

Grandeur Nom Symbole Dimension
Longueur Mètre 𝑚 𝐿
Masse Kilogramme 𝑘𝑔 𝑀
Temps Seconde 𝑠 𝑠
Intensité du courant électrique Ampère 𝐴 𝐼
Température thermodynamique Kelvin 𝐾 Θ
Quantité de matière Mole 𝑚𝑜𝑙 𝑁
Intensité lumineuse Candela 𝑐𝑑 𝐽

Le tableau I.2 montre les dimensions et les unités en système internationale et anglo-
saxon de quelques grenadeurs physiques.

3
 Définitions et propriétés des fluides

Tab (I.2) Dimensions et unités en système international (SI) et anglo-saxon.

Grandeur Symbole Dimension S.I Anglo-Saxon
Surface 𝐴 𝐿2 𝑚2 𝑓𝑡 2
Volume ∇ 𝐿3 𝑚3 𝑓𝑡 3
Vitesse 𝜈 𝐿⁄𝑇 𝑚⁄𝑠 𝑓𝑡⁄𝑠𝑒𝑐
Accélération 𝑎 𝐿 ⁄𝑇 2 𝑚⁄𝑠 2 𝑓𝑡⁄𝑠𝑒𝑐 2
Vitesse angulaire Ω 𝑇 −1 𝑠 −1 𝑠𝑒𝑐 −1
Force 𝐹 𝑀𝐿⁄𝑇 2 𝑘𝑔. 𝑚⁄𝑠 2 𝑠𝑙𝑢𝑔. 𝑓𝑡⁄𝑠𝑒𝑐 2 ou 𝑙𝑏
Densité 𝜌 𝑀⁄𝐿3 𝑘𝑔⁄𝑚3 𝑠𝑙𝑢𝑔⁄𝑓𝑡 3
Poids spécifique 𝛾 𝑀⁄𝐿2 𝑇 2 𝑁⁄𝑚3 𝑙𝑏⁄𝑓𝑡 3
Fréquence 𝑓 𝑇 −1 𝑠 −1 𝑠𝑒𝑐 −1
Pression 𝑝 𝑀⁄𝐿 𝑇 2 𝑁⁄𝑚2 ou Pascal 𝑙𝑏⁄𝑓𝑡 2
Contrainte 𝜏 𝑀⁄𝐿 𝑇 2 𝑁⁄𝑚2 ou Pascal 𝑙𝑏⁄𝑓𝑡 2
Tension superficielle 𝜎 𝑀⁄𝑇 2 𝑁⁄𝑚 𝑙𝑏⁄𝑓𝑡
Travail 𝑊 𝑀𝐿2 ⁄𝑇 2 𝑁. 𝑚 ou Joule 𝑓𝑡. 𝑙𝑏
Energie 𝐸 𝑀𝐿2 ⁄𝑇 2 𝑁. 𝑚 ou Joule 𝑓𝑡. 𝑙𝑏
Taux de chaleur 𝑄̇ 𝑀𝐿2 ⁄𝑇 3 𝐽/𝑠 𝐵𝑡𝑢⁄𝑠𝑒𝑐

I.2.3 Système CGS :

Système CGS Le système centimètre-gramme-seconde (CGS abrégé ou CGS) est un
système métrique d'unités physiques basées sur le centimètre comme unité de
longueur, le gramme comme unité de masse, et la seconde comme unité de temps. Il
permet d’exprimer toutes les grandeurs de la mécanique.

Exemple : Résistance à l’avancement d’un navire : 𝑹𝑻

Les Concepts résistance et propulsion sont indissociable et étroitement liés entre eux
car pour déterminer la puissance, il faut en premier connaitre la résistance à
l’avancement du navire qui est un sujet très important de l’hydrodynamique navale
car il s’agit de déterminer l’effort à vaincre pour déplacer un engin maritime, et par
conséquent la motorisation requise pour propulser le navire à une vitesse déterminée.

Résistance La puissance La puissance La puissance

totale RT de remorquage de l’hélice du moteur
PE

Fig (I.2) Etapes à suivre pour le calcul de la puissance du moteur à bord.

D’après les hypothèses faites sur les paramètres les plus influents pour décrire un
phénomène physique, l’analyse dimensionnelle, permet de déterminer les nombre
sans dimensions indépendants.

4
 Définitions et propriétés des fluides

Donc, la première étape pour étudier un phénomène physique est l´identification des
variables importantes. La relation mathématique entre ces variables constitue une loi
de la physique.
Pour des systèmes complexes ceci peut être difficile. Dans ce cas une méthode de
modélisation telle que l´analyse dimensionnelle est indispensable.

Les propriétés physiques (les paramètres) qui influent sur la résistance à

l’avancement 𝑅𝑇 doivent être au minimum:

𝜌 Masse volumique de l’eau (𝑘𝑔⁄𝑚3 )

𝑈 Vitesse d’avance (𝑚⁄𝑠)
𝐿 Longueur caractéristique du navire (𝑚)
𝜇 Viscosité de l’eau (𝑘𝑔⁄𝑚. 𝑠)
𝑔 Accélération de la pesanteur (𝑚⁄𝑠 2 )
𝑃 Pression (𝑘𝑔⁄𝑚. 𝑠 −2 )

On peut donc supposer que la résistance à l’avancement 𝑅𝑇 est en fonction des

paramètres cités ci-dessus :
𝑅𝑇 = 𝑓( 𝜌, 𝑈, 𝐿, 𝜇, 𝑔, 𝑃)

L’équation peut être réécrite comme une relation entre nombre adimensionnels à
partir des paramètres intervenants dans cette équation (6 paramètres), cette
homogénéité est connue sous le nom de Vaschy-Buckingham ou théorème Pi, est un
des théorèmes de base de l'analyse dimensionnelle, qui suppose que si une équation
physique met en jeu n variables physiques, celles-ci dépendant de (𝑘) unités
fondamentales alors il existe une équation équivalente mettant en jeu (𝑛 −
𝑘) variables sans dimensions construites à partir des variables originelles.
Comme le nombre de grandeurs fondamentales est égal à 3, les nombres
adimensionnelles entrant dans la relation finale doivent être égal à 6-3=3.
Application :
La formule de Stokes-Oseen pour le calcul de la force de trainée 𝐹 d’une sphère de
diamètre d et pour un écoulement de faible vitesse 𝑣, est donnée par 𝐹 = 3𝜋𝜇 𝑑 𝑣 +
(9𝜋/16)(𝜌𝑣 2 𝑑 2 ). La formule est-elle valide dimensionnellement ?

Réponse :
1. [𝑀𝐿. 𝑇 −2 ] = [𝑀𝐿. 𝑇 −2 ] la formule est valide dimensionnellement.

I.3 Définition d’un fluide :

Un fluide peut être considéré comme étant une substance formé d'un grand nombre
de particules matérielles, très petites et libres de se déplacer les unes par rapport aux
autres. C’est donc un milieu matériel continu, déformable, sans rigidité et qui peut
s'écouler. Les forces de cohésion entres particules élémentaires sont très faibles de
sorte que le fluide est un corps sans forme propre qui prend la forme du récipient qui

5
 Définitions et propriétés des fluides

le contient, par exemple: les métaux en fusion sont des fluides qui permettent par
moulage d'obtenir des pièces brutes de formes complexes.

I.4 Fluide parfait et fluide réel :

Un fluide est dit parfait s'il est possible de décrire son mouvement sans prendre en
compte les effets de frottement. Contrairement à un fluide parfait, qui n’est qu’un
modèle pour simplifier les calculs, pratiquement inexistant dans la nature, dans un
fluide réel les forces tangentielles de frottement interne qui s’opposent au glissement
relatif des couches fluides sont prises en considération.

I.5 Fluide Compressible et incompressible:

Un fluide est dit incompressible lorsque le volume occupé par une masse donné ne
varie pas en fonction de la pression extérieure. Les liquides peuvent être considérés
comme des fluides incompressibles (eau, huile, etc.). Un fluide est dit compressible
lorsque le volume occupé par une masse donnée varie en fonction de la pression
extérieure. Les gaz sont des fluides compressibles. Par exemple, l’air, l’hydrogène, le
méthane à l’état gazeux, sont considérés comme des fluides compressibles.

I.6 Les caractéristiques d’un fluide :

I.6.1 La viscosité :

𝑢 = 𝑢( 𝑦)
Couche fluide supérieur𝑒
A A'
𝐹
𝜏 = 𝜏( 𝑦)
Couche fluide inférieure

Fig (I.3) Profil de vitesse (écoulement autour d’une plaque plane).

La viscosité est une propriété du fluide, elle se manifeste quand le fluide est en
mouvement. On considère un écoulement avec un profil de vitesse montré ci-dessus.
Prenant une section adimensionnelle AA’, Les particules de la couche fluide avec une
vitesse supérieure se glissent sur les particules fluides de la couche inférieure, ces
dernières et avec le frottement freinent les particules de la couche supérieurs (Force
𝐹 exercée entre les 2 couches due aux frottements).
La force de frottement 𝐹 qui s'exerce à la surface de séparation de ces deux couches
s'oppose au glissement d'une couche sur l'autre. Chaque couche fluide ne s'écoule pas
à la même vitesse.
6
 Définitions et propriétés des fluides

𝑑𝑢
𝐹 = −𝐴𝜇
𝑑𝑦
Relation proportionnelle :
𝑑𝑢
𝜏∝
𝑑𝑦
La loi de Newton sur la viscosité :
𝑑𝑢
𝜏=𝜇
𝑑𝑦

𝜇: Coefficient de viscosité ou viscosité dynamique.

𝑑𝑢⁄𝑑𝑦 : Taux de cisaillement.
𝜏: Contrainte de cisaillement.

Unité : Dans le système international (𝑆𝐼), l'unité de viscosité est le Poiseuille (𝑃𝑙) :
1 𝑃𝑙 = 1 𝑘𝑔⁄𝑚. 𝑠
On trouve encore les tables de valeurs numériques le coefficient de viscosité dans un
ancien système d'unités (CGS) :
L'unité est le Poise (𝑃𝑜); 1 𝑃𝑙 = 10 𝑃𝑜
Dans de nombreuses formules apparaît le rapport de la viscosité dynamique 𝜇 et de
la masse volumique 𝜌, ce rapport est appelé viscosité cinématique.
𝜇
𝜈=
𝜌
Unité SI : m2/s
Système CGS : le Stoke (St) 1m2/s = 106 cSt.

La viscosité cinématique caractérise le temps d'écoulement d’un liquide. Par contre,

la viscosité dynamique correspond à la réalité physique du comportement d’un fluide
soumis à une sollicitation (effort). En d’autre terme, cette dernière exprime la «
rigidité » d’un fluide à une vitesse de déformation en cisaillement.

I.6.2 La masse volumique :

C’est le rapport de la masse (𝑚) sur le volume (∇) du corps considéré :

𝑚
𝜌=
∇

Unité SI : 𝑘𝑔⁄𝑚3

Comme les solides et liquides sont incompressibles, le volume peu dépendre de la

température et pas du tout de la pression. 𝜌 diminue un peu si la température
augmente.

7
 Définitions et propriétés des fluides

Les gaz sont compressibles : ∇ augmente si 𝑇 augmente et ∇ diminue si 𝑃 augmente.

𝜌 varie d’un fluide à un autre.

Exemple : 𝜌𝐻2 𝑂 =1000 𝑘𝑔⁄𝑚3

𝜌𝐴𝑖𝑟 = 1.225 𝑘𝑔⁄𝑚3

I.7 Calcul du couple des disques en fonction de la contrainte de

cisaillement :
Pour mesurer la viscosité, on considère un cylindre long tournant à l'intérieur d'un
deuxième cylindre, comme le montre la figure I.4. Afin de faire tourner le cylindre
intérieur avec la vitesse de rotation Ω, un couple 𝑇 doit être appliqué.
La vitesse du cylindre interne est 𝑅Ω et la vitesse du cylindre extérieur est nulle. La
distribution de vitesse dans l'intervalle h entre les cylindres est une distribution
linéaire telle que représentée, de sorte que :

Fig (I.4) Cylindre mobile autour d’un cylindre fixe.

𝑑𝑢 𝑅𝛺
𝜏=𝜇 =𝜇
𝑑𝑦 ℎ

Le couple 𝑇 est calculé comme suit :

𝑇 = 𝑟. 𝜏. 𝐴
𝑇 = 𝑅. 𝜏. 2𝜋𝑅𝐿

𝑅𝛺 𝑅3 𝛺𝐿𝜇
𝑇 = 𝑅. 𝜇 . 2𝜋𝑅𝐿 = 2𝜋
ℎ ℎ

I.8 Ecoulements laminaires et turbulents :

La résistance interne dans un fluide à s’écouler est appelée viscosité cinématique (qui
fait référence à la vitesse et au mouvement). Cette viscosité est la cause de la

8
 Définitions et propriétés des fluides

résistance hydrodynamique de l’eau sur la surface d’un bateau. En 1883, le physicien

Reynolds a réalisé des expériences sur l’écoulement d’un fluide dans une conduite
cylindrique, en augmentant progressivement la vitesse du fluide et en mettant
l’évidence de 02 régimes d’écoulement laminaire et turbulent.
Reynolds a montré que le paramètre qui permettait de déterminer si l’écoulement est
laminaire ou turbulent est un nombre sans dimension appelé nombre de Reynolds.
𝜌𝑉𝐿
𝑅𝑒 =
𝜇

𝜌: La masse volumique.
𝑈: La vitesse du fluide.
𝐿: Longueur du navire.
𝜇: La viscosité dynamique.

Le nombre de Reynolds mesure l’importance relative des forces d’inerties et les forces
visqueuses. Pour n’importe quel écoulement laminaire, il existe une valeur bien
déterminée du nombre de Reynolds pour laquelle l’existence de ce régime
d’écoulement est menacée, et pour des valeurs inférieures à ce nombre critique de
Reynolds l’écoulement est stable et il est décrit comme un ensemble de lames
(couches) glissant les unes sur les autres d’une manière ordonnée.
Les particules suivent des trajectoires parallèles et le fluide se comporte comme s'il
était constitué de lames superposées glissant parfaitement les unes sur les autres. Au
voisinage d'une surface solide, la vitesse de l'écoulement ralentit au fur et à mesure
que l'on s'en approche pour finalement s'annuler au contact de celle-ci. Cette couche
d'air freinée s'appelle " la couche limite ". Pour un nombre de Reynolds élevé,
l’écoulement est toujours turbulent c.-à-d aléatoire et désordonné, et impossible
d’analyser mais aménageables pour étudier des valeurs moyennes. La figure ci-dessus
montre les lignes de courant dans une conduite pour 02 régimes différents.

(a)

(b)

Fig (I.5) (a) écoulement laminaire (b) écoulement turbulent.

I.9 Ecoulements internes et externes :

I.9.1 Ecoulements internes : Exemple : écoulement dans une conduite.

9
 Définitions et propriétés des fluides

Fig (I.6) Ecoulement dans une conduite (écoulement interne).

 La force de viscosité est dominante sur tout le champ de l’écoulement.

 L’influence de la viscosité (force de viscosité) s’étend du premier point jusqu’au
dernier point de l’écoulement.
 Gradient de vitesse sur tout le champ de l’écoulement.

L’expérience de Reynolds a montré que si le nombre de Reynolds 𝑅𝑒 < 2300,

l’écoulement est considéré laminaire, et 𝑅𝑒 > 2300 il devient facilement turbulent.

Laminaire : 𝑅𝑒 < 2300

Transitoire: 2300 < 𝑅𝑒 < 4200

10
 Définitions et propriétés des fluides

Turbulent : 𝑅𝑒 < 4200

Fig (I.7) Expérience de Reynolds (différents régimes d’écoulements).

I.9.2 Ecoulements externes :

Pour les écoulements externes, la zone proche des parois, la force de viscosité est
dominante, loin des parois, les forces de viscosité sont négligeables.

Fig (I.8) Ecoulement autour d’un cylindre (écoulement externe).

On peut diviser l’écoulement externe en zones :

Zone 1 : Couche limite (influence de la viscosité).
Zone 2 : Région invicide (les forces de viscosité sont négligeables).
Pour les écoulements externes si : 𝑅𝑒 < 5. 105 l’écoulement est considéré laminaire.
Il est dit turbulent lorsque le nombre de Reynolds 𝑅𝑒 > 5. 105 .

I.10 Exemples évolution de la forme géométrique (aérodynamique)

I.10.1 Forme géométrique d’une voiture :
Le but d’étudier l’écoulement externe, est de l’utiliser en design. La forme s’est
développée entre 1920 jusqu’à maintenant, on remarque la réduction du coefficient
de trainée (𝐶𝑑 ), synonyme d’une économie de fuel, réduire les émissions…

11
 Définitions et propriétés des fluides

Fig (I.9) Evolution des formes géométriques des voitures.

I.10.2 Forme géométrique d’une balle :

La balle de golf généré une trainée moins importante par rapport à une balle lisse
grâce à sa forme géométrique constituée des alvéoles sur toute la surface de la balle,
lorsque l’écoulement est générée autour de la balle de golf, les alvéoles se remplissent
d’air ce qui permet à retarder le point de décollement, synonyme d’une réduction de
trainée.

(a) Vortex géants générant une

trainée supérieur.

(b) Vortex petits générant une

trainée inférieure.

Fig (I.10) (a) écoulement autour d’une balle lisse (b) écoulement autour d’une balle de golf.

La même chose pour les carènes des navires, on étude l’écoulement externes pour
réduire le coefficient de résistance total, en jouant sur le conception géométrique.

12
 Définitions et propriétés des fluides

I.11 Etude des écoulements externes :

Il existe 2 approches :
I.11.1 Etude expérimentale et visuel : Moyens expérimentaux
Exemples :
 Bassin d’essais des carènes pour les maquettes navires.
 Tunnel hydrodynamique pour les hélices.
 Soufflerie.
But : Les mesures des champs de vitesse, la pression, le coefficient de la résistante
totale, la poussée et le couple des hélices…
I.11.2 Etude théorique :
Résoudre les équations gouvernantes l’écoulement (purement théorique).
 Conservation de masse.
 Conservation de quantité de mouvement.
 Conservation d’énergie (Première loi de la thermodynamique).
 2 ème loi de la thermodynamique.

Numérique : CFD (Computational Fluid Dynamics).

13
 Définitions et propriétés des fluides

I.12 Exercices avec réponses :

Exercice 01 :
Un fluide Newtonien ayant une densité 𝑑 = 0.92 et une viscosité cinématique
𝜈 = 4. 10−4 𝑚2 /𝑠 s’écoule sur une plaque plane en régime permanant. Le profil de
vitesse près de la surface (voir annexe I.1), et suit la loi 𝑢(𝑦)⁄𝑈 = 𝑠𝑖𝑛(𝜋𝑦⁄2𝛿).
Déterminer l’amplitude et la direction de la contrainte de cisaillement développée
sur la plaque en fonction de 𝑈 et 𝛿.

Réponse :
𝑈
1. 𝜏 = 0,578 𝑁. 𝑚−2
𝛿

Exercice 02 :
Un fluide à une viscosité absolue de 0.048 𝑃𝑎. 𝑠 et une densité de 0,913 (voir
annexe I.2). Calculer le gradient des vitesses et l’intensité de la contrainte
tangentielle à la paroi et aux points situés à 25 𝑚𝑚, 50 𝑚𝑚, et 75 𝑚𝑚 et de celle-ci,
en admettant (𝑎) une distribution de vitesse linéaire, (𝑏) une distribution de
vitesse parabolique. La parabole de la figure a son sommet en (𝐴). L’origine est en
(𝐵).

Réponse :
𝑑𝑈
1. = 15
𝑑𝑦
2. 𝜏 = 0,72 𝑃𝑎

Exercice 03 :
Une plaque plane de dimensions 0,5 𝑚 x2 𝑚 est remorquée horizontalement sur
une paroi à une vitesse de 5 𝑚⁄𝑠. Sachant que l’huile 𝑆𝐴𝐸 − 30 à 38°𝐶 sépare la
plaque plane et la paroi d’une couche de 2 𝑚𝑚, et supposant que le profil de
vitesse est linéaire, déterminer la force de frottement F. (F=𝜏. 𝐴 ) (voir annexe I.3).

Réponse :
1. 𝐹 = 250 𝑁

Exercice 04 :
De l’huile ayant une viscosité remplit l'espace 𝑌 situé entre un cône rotatif et la
base stationnaire (voir annexe I.4), le cône tourne à une vitesse angulaire 𝛺. On
négligeant les contraintes sur le fond du cône et on suppose que la distribution de
la vitesse est linéaire. Donner l’expression de : la contrainte de cisaillement 𝜏 en
fonction de 𝛼 et 𝑌 ; la surface élémentaire 𝑑𝐴 en fonction de 𝛼 et 𝑦, la force
élémentaire 𝑑𝐹 en fonction de 𝛼, 𝜇 et 𝑌 ; le couple élémentaire 𝑑𝑇 en fonction de
𝛼, 𝜇 et 𝑌 et le couple 𝑇 en fonction de 𝑎 𝑒𝑡 𝑏. On donne : 𝑑𝑦 = 𝑑𝑠 𝑐𝑜𝑠 𝛼.

Réponse :
𝑑𝑣 𝜇.𝑦.𝑡𝑔𝛼.𝛺 𝑡𝑔𝛼 2𝜋.𝜇.𝛺.𝑡𝑔2 𝛼
1. = 2. 𝑑𝐴 = 2𝜋𝑦. (𝑐𝑜𝑠𝛼) 𝑑𝑦 3. 𝑑𝐹 = ( ) 𝑦 2 𝑑𝑦
𝑑𝑦 𝑌 𝑌.𝑐𝑜𝑠𝛼
2𝜋.𝜇.𝛺.𝑡𝑔3 𝛼 2𝜋.𝜇.𝛺.𝑡𝑔3 𝛼
4. 𝑑𝑇 = ( ) 𝑦 3 𝑑𝑦 5. 𝑇 = [(𝑎 + 𝑏)4 − 𝑎4 ]
𝑌.𝑐𝑜𝑠𝛼 𝑌.𝑐𝑜𝑠𝛼

14
 Définitions et propriétés des fluides

Exercice 05 :
Un pivot de palier utilisé sur l'arbre d'un instrument électrique est indiqué sur la
figure en annexe (voir anx I.5). De l’huile ayant une viscosité remplit l'espace 𝑏
situé entre l’arbre rotatif et la base stationnaire. Donner l’expression du couple
visqueux sur l'arbre en fonction de 𝜇, 𝛺, 𝜃, 𝑅 𝑒𝑡 𝑏 (On suppose que la distribution de
la vitesse est linéaire).

Réponse :
𝜋𝜇Ω
1. 𝑇 = 2𝑏𝑠𝑖𝑛𝜃 𝑅 4

15
Chapitre II

TENSION SUPERFICIELLE

II TENSION SUPERFICIELLE ………………………………………..……………….. 16

II.1 Introduction .………………………………………………….………………….……. 17
II.2 Forces de cohésion et d’adhésion ……………………………………………… 17
II.3 La loi de Young Laplace …………………………………………………………… 20
II.4 Angles de contact et capillarité ……………………………………..........…… 21
II.5 Phénomène de mouillabilité ……………………………………………..……… 23
II.6 Exercices avec réponses …………………………………………………….……… 25

16
 Tension superficielle

II.1 Introduction:
L’interface entre le liquide et le gaz est caractérisé par un paramètre qui dépend de la
nature physico-chimique du fluide et du gaz ainsi que la température et de la pression
dans le gaz. Ce paramètre est appelé tension superficielle. Ce paramètre est lié aux
interactions moléculaires d’un fluide (IMA- Intermolecular attraction). La molécule
d’eau a des charges (+) et (-), les charges différentes crées une attraction
intermoléculaire.
La figure ci-dessous représente qualitativement les efforts intermoléculaires des
différentes molécules dans un liquide avec une surface libre. On peut remarquer que
lorsque la molécule est totalement immergée dans le liquide, cette dernière est en
équilibre puisque les efforts s’annulent deux à deux dans les deux directions du
repère cartésien (attraction est dans tous les directions). En revanche lorsque la
molécule s’approche de plus en plus de la surface libre, il y a alors un déséquilibre des
forces dû à la tension de surface (les deux forces horizontales s’annulent) et par
conséquent la molécule est tirée vers le bas.

Fig (II.1) Efforts intermoléculaires des différentes molécules dans un liquide avec une surface libre.

II.2 Forces de cohésion et d’adhésion :

Le phénomène de tension de surface provient de deux forces intermoléculaires :
 Forces de cohésion : ce sont les forces d’attraction entre les molécules du
liquide ;

17
 Tension superficielle

 Forces d’adhésion : ce sont les forces entre les molécules de différents liquides
ou entre un liquide et un solide.

Exemples :

(a.1) (a.2) (b)

Fig (II.2) Exemples de la tension superficielle.
(a.1) et (a.2) la pièce d’alluminium et l’aiguille ne coulent pas au fond du verre.
(b) l’insecte peut se déplacer.

Pour en comprendre la signification physique considérons l’expérience représentée

sur la figure (II.3). Elle consiste à étirer un film de savon entre deux tiges, distantes
de 𝑏, en venant exercer un effort 𝐹 dans la direction horizontale. Notons bien que
cette force ne sert pas seulement à étirer le film mais aussi à maintenir le film avec
une aire constante. En effet sans l’action de cette force, le film tendra à se rétracter
pour retrouver une surface minimale.

Fig (II.3) Experience d’un film de savon entre 02 tiges.

Pour déplacer le coté 𝐴𝐵 d’une distance 𝛿𝑥 , il faut fournir un travail :

𝛿𝑊 = 𝐹 𝛿𝑥

18
 Tension superficielle

Dans ce cas l’aire de la surface du fil va être augmentée de la quantité (en tenant
compte les 2 faces du film).
𝛿𝐴 = 2𝐿 𝛿𝑥

La tension superficielle est définie comme le travail produit par unité de surface (c’est
une force de traction par unité de longueur agissant sur un élément de surface) :
𝛿𝑊
𝜎=
𝛿𝐴

On voit que pour l’étirage du film de savon on a alors :

𝐹
𝜎=
2𝐿

L’unité de mesure de la tension superficielle est 𝑁/𝑚.

On notera que la tension superficielle diminue lorsque la température augmente :
cela s’explique par le fait que les forces de cohésion intermoléculaires sont réduites
par l’agitation thermique.
Exemple :
(a) Avec la fibre de verre fabriquant un anneau sur lequel on fixe un fil de coton, le
fil ne doit pas être tendu mais au contraire, il tendra assez librement.

(b) Plongeant l’anneau dans de l’eau additionnée de produit vaisselle, lorsqu’on le

ressort, un film de savon apparait de chaque côté du fil qui reste non tendu.

(c) Crevant alors l’une des surfaces savonneuse cette fois le fil se tend en formant
un arc de cercle, en effet le film savonneux restant tire sur le fil pour pouvoir
diminuer sa surface.

(d) Si on crève ce film, le fil de coton le plus soumis à la tension superficielle créé
par l’eau savonneuse épand de nouveau librement.

(a) (b)

19
 Tension superficielle

(c) (d)
Fig (II.4) L’expérience de l’anneau.

II.3 La loi de Young Laplace :

On considère le cas d’une bulle sphérique de diamètre 𝐷 , les forces de tension
superficielle, qui sont dirigées vers l’intérieur de la goutte, exercent une compression
à l’intérieur de celle-ci. La pression 𝑃𝑖𝑛 dans la goutte est donc supérieure à celle du
milieu extérieur 𝑃0 .

𝐹𝑐 𝐹𝑎
𝐹𝑎 < 𝐹𝑐

𝑃𝑖𝑛 > 𝑃𝑜𝑢𝑡

Fig (II.5) Bulle sphérique d’eau.

Les forces de tension superficielle :

𝐹𝑐 = 𝐹 = 𝜎. 𝐿 = 𝜎. (2𝜋𝑅)

Les forces de pression :

∆𝑃(𝜋𝑅 2 ) = (𝑃𝑖𝑛 − 𝑃𝑜𝑢𝑡 )(𝜋𝑅 2 )
2𝜎
𝑃𝑖𝑛 − 𝑃𝑜𝑢𝑡 =
𝑅

Cette dernière relation est appelé la loi de Young-Laplace et elle montre que la
pression à l’intérieur d’une goutte liquide est donc d’autant plus grande que la goutte
est petite.

20
 Tension superficielle

Remarque :
Une bulle est formée d’une membrane comportant deux surfaces (interne et externe)
supposées de même rayon 𝑅, chacune d’elles étant le siège d’une tension superficielle.
Les forces de pression qui ont globalement tendance à faire dilater la bulle, doivent
donc compenser les forces de tension superficielle sur les deux interfaces.

𝐹𝑐 𝐹𝑎
𝐹𝑎 < 𝐹𝑐

𝑃𝑖𝑛 > 𝑃𝑜𝑢𝑡

Fig (II.6) Bulle d’eau savonneuse.

Les forces de tension superficielle :

𝐹𝑐 = 𝐹 = 𝜎. 𝐿 = 2𝜎. (2𝜋𝑅)

Les forces de pression :

∆𝑃(𝜋𝑅 2 ) = (𝑃𝑖𝑛 − 𝑃𝑜𝑢𝑡 )(𝜋𝑅 2 )
4𝜎
𝑃𝑖𝑛 − 𝑃𝑜𝑢𝑡 =
𝑅

II.4 Angles de contact et capillarité :

Un tube de verre de faible diamètre est plongé dans un liquide mouillant, de l’eau par
exemple. Dans le tube, le niveau du liquide est supérieur au niveau de la surface libre
du récipient. Le ménisque concave fait un angle 𝜃 avec la surface du tube.
L’ascension capillaire est due aux forces superficielles appliquées en tout point du
contour du ménisque. La résultante 𝐹 de ces forces équilibre le poids 𝑃 du liquide
soulevé. L’élévation du liquide dans le tube compense la différence de pression entre
les deux côtés de la paroi (loi de Laplace).

21
 Tension superficielle

Fig (II.7) Formation d’un ménisque dans un tube de verre plongé dans un liquide mouillant.

Le poids de la colonne de liquide dans le tube :

𝑃 = 𝑚𝑔 = 𝜋𝑅 2 ℎ𝜌𝑔

Equilibré par la force de tension superficielle :

𝐹 = 2𝜋𝑅 𝜎 𝑐𝑜𝑠𝜃

S’exerçant sur la ligne de raccordement entre le liquide et la paroi du tube.

Et on obtient ainsi la relation que l’on appelle Loi de Jurin :

2 𝜎 𝑐𝑜𝑠𝜃
ℎ=
𝑅𝜌𝑔

𝑅: Rayon intérieur du tube.

𝜌: Masse volumique du fluide.
𝑔: Accélération de la pesanteur.
𝜃: Angle de raccordement liquide/solide.

Dans le cas du mouillage parfait 𝑐𝑜𝑠𝜃 = 1

Si l’angle 𝜃 dépasse 90°, la loi de Jurin donne ℎ négatif. On parle alors de dépression
capillaire. C’est le cas du mercure au contact du verre. Cette fois les forces de
cohésion sont supérieures aux forces d’adhésion, le liquide ne mouille pas les parois
du tube. Le niveau du liquide s’abaisse dans le tube au-dessous du niveau de la
surface libre du récipient. Le ménisque est convexe et forme l’angle 𝜃 > 90° avec la
paroi du tube.

22
 Tension superficielle

Fig (II.8) Ascension et dépression capillaires.

Tab (II.1) Angle de contact et capilarité entre 02 matériaux différents.

Interface dans l’air Angle de contact Ascension capillaire
Eau-verre 0
Alcool-verre 0
Mercure-verre 140
Kérosène-verre 26
Eau- acier 90 Nulle
Eau-bois 0

Application:
Trouver la hauteur de capillarité d’un tube de verre de 1 𝑚𝑚 rayon à une température
de 20°𝐶, le tube est immergé dans un récipient d’eau. L’angle d’interface formé entre
le liquide/air/verre est 0°, calculer la hauteur ℎ pour le cas du mercure si l’angle
d’interface est 130°.
σeau = 0.0728 𝑁/𝑚
σmercure = 0.514 𝑁/𝑚
Réponse :
1. ℎ𝑒𝑎𝑢 = 14,84 𝑚𝑚
2. ℎ𝑚𝑒𝑟𝑐𝑢𝑟𝑒 = −5 𝑚𝑚

II.5 Phénomène de mouillabilité :

Lorsque l’on dépose une goutte d’eau sur une plaque en téflon puis sur une plaque en
verre, on constate que la goutte ne s’étale pas du tout de la même façon. Dans le cas
du téflon, la goutte à tendance à avoir une conformation sphérique pour limiter au
maximum la surface de contact entre le domaine fluide et le domaine solide. Dans le
cas du verre la goutte a tendance à s’étaler. On appelle ce phénomène la mouillabilité
d’une surface solide par un fluide. On dira donc que le verre est mouillable à l’eau,
alors que le téflon ne l’est pas. On se reportant sur la configuration représentée sur la

23
 Tension superficielle

figure, on nommera 𝜃 l’angle de contact (ou angle de mouillabilité). On distingue

alors les configurations suivantes :
 Soit s’étaler, on dit que le liquide mouille parfaitement le solide 𝜃 = 0°
 Soit forme une lentille avec deux cas de figure :
1) Le liquide mouille imparfaitement le solide 𝜃 < 90°
2) Le liquide ne mouille pas le solide 𝜃 > 90°

(a) 𝜃 = 0° (b) 𝜃 < 90° (c) 𝜃 > 90°

Fig (II.9) Mouillabilté d’une goutte d’eau sur différents matériaux.

La ligne de raccordement caractérise le périmètre de la surface de contact entre fluide

et solide. Sur cette ligne nous avons coexistence des trois phases : liquide (𝑙), solide
(𝑠) et vapeur (𝑣).
Trois paramètres sont donc à prendre en compte :
 La tension superficielle 𝜎𝑠𝑙 entre le solide et le liquide ;
 La tension superficielle 𝜎𝑙𝑣 entre le liquide et sa phase vapeur ;
 La tension superficielle 𝜎𝑠𝑣 entre le solide et la vapeur.
L’équilibre de la goutte se traduit par :
𝜎𝑠𝑣 = 𝜎𝑠𝑙 + 𝜎𝑙𝑣 𝑐𝑜𝑠𝜃

Soit :
𝜎𝑠𝑣 − 𝜎𝑠𝑙
𝑐𝑜𝑠𝜃 =
𝜎𝑙𝑣

Fig (II.10) Goutte d’eau sur une feuille de cactus.

24
 Tension superficielle

II.6 Exercices avec réponses:

Exercice 01:
Un tube de verre immergé dans un récipient contenant du mercure (voir annexe
II.1). Calculer la force de tension superficielle?

On donne:
𝜎 = 0.514 𝑁/𝑚

Réponse :
1. 𝐹 = 0,0623 𝑁

Exercice 02:
Quelle est la différence de pression ∆𝑃 entre l’intérieur et l’extérieur d’une
gouttelette d’eau ayant une tension superficielle 𝜎𝑒𝑎𝑢 = 0.0728 𝑁/𝑚 et un diamètre
𝑑 = 7.6 10−3 𝑐𝑚.
Calculer la tension superficielle pour le cas d’une bulle de savon sachant que son
rayon est de 𝑑 = 4 𝑐𝑚 et la différence de pression ∆𝑃 = 2,5 𝑃𝑎.

Réponse :
1. ∆𝑃 = 3831,6 𝑃𝑎
2. 𝜎 = 0.025 𝑁/𝑚

Exercice 03:
Développer une formule de remontée capillaire entre deux tubes de verre
concentriques de rayon 𝑟0 et 𝑟1en function de l’angle de capilarité et le poids
spécifique 𝛾 (voir annexe II.2).

Réponse :
2𝜎.𝑐𝑜𝑠𝜃
1. ℎ = 𝛾(𝑟0 −𝑟1 )

Exercice 04:
Quelle est la force nécessaire pour soulever un mince anneau de fil de 6 𝑐𝑚 de
diamètre à partir d'une surface d'eau à 20°𝐶.

Réponse :
1. 𝐹 = 0,0274 𝑁

Exercice 05:
Un petit jet circulaire de mercure de 200 𝜇𝑚 de diamètre sort d'une ouverture
(voir annexe II.3). Quelle est la différence de pression entre l'intérieur et
l'extérieur du jet à 20°𝐶.

Réponse :
1. ∆𝑃 = 5140 𝑃𝑎

25
Chapitre III

LES EQUATIONS DE CONSERVATION

III LES EQUATIONS DE CONSERVATION ………………………….…………. 26

III.1 Conservation de la masse (Equation de continuité) ……………………. 27
III.1.1 Le Théorème de Transport de Reynolds ……………………………. 27
III.1.2 Forme intégrale de l’équation de continuité ……………………….. 27
III.1.3 Forme différentielle de l’équation de continuité …………………. 28
III.1.4 Forme vectorielle de l’équation de continuité ……………………… 30
III.2 Conservation de la quantité de mouvement (Equations d’Euler) … 31
III.2.1 La 2 ème loi de Newton ……………………………..………………….…. 32
III.2.2 Forme vectorielle de l’équation d’Euler ……………………………… 35
III.3 Conservation de la quantité de mouvement (Equations de Navier
Stokes) …………………………………………………………………..………………. 36
III.3.1 Le tenseur de contraintes total …………….………………………...…. 36
III.3.2 Les contraintes normales et les contraintes tangentielles ... 36
III.3.3 Forme vectorielle de l’équation de Navier Stokes ………………… 38
III.4 Exercices avec réponses ………………………………………………………….… 40

26
 Les equations de conservation

III.1 Conservation de la masse (Equation de continuité) :

III.1.1 Le Théorème de Transport de Reynolds :

𝑁: Propriété
𝑛: 𝑁 par unité de masse

𝑑𝑁 𝜕𝑁
| = ⃗⃗⃗ 𝜂̂ ) 𝑑𝐴
| + ∫ 𝑛. 𝜌(𝑉.
𝑑𝑡 𝑆𝑦𝑠 𝜕𝑡 𝑐𝑣 𝑐𝑠

𝜂: Vecteur normal sur la surface

⃗⃗⃗ 𝜂̂ . 𝑑𝐴: Débit volumique
𝑉.
𝜌. ⃗⃗⃗
𝑉. 𝜂̂ . 𝑑𝐴: Débit massique
⃗⃗⃗ 𝜂̂ . 𝑑𝐴: Transport total de la propriété
𝑛. 𝜌. 𝑉.

Pour un système :(conservation de la masse)

𝑑𝑀𝑠𝑦𝑠
=0 Fig (III.1) Vecteur normal et vitesse
𝑑𝑡
sur une surface élémentaire 𝑑𝐴.

𝑁=𝑚
𝑛=1

Forme générale :

𝑑𝑚 𝜕𝑚
| = ⃗⃗⃗⃗⃗𝑅 . 𝜂̂ ) 𝑑𝐴 = 0
| + ∫ 𝜌 (𝑉
𝑑𝑡 𝑆𝑦𝑠 𝜕𝑡 𝑐𝑣 𝑐𝑠

𝜕𝑚 𝜕
| = ∫ 𝜌 𝑑𝑉
𝜕𝑡 𝑐𝑣 𝜕𝑡 𝑐𝑣

III.1.2 Forme intégrale de l’équation de continuité

𝑑𝑚 𝜕
| = ∫ 𝜌 𝑑𝑉 + ∫ 𝜌(⃗⃗⃗⃗⃗
𝑉𝑅 . 𝜂̂ ) 𝑑𝐴 = 0
𝑑𝑡 𝑆𝑦𝑠 𝜕𝑡 𝑐𝑣 𝑐𝑠

27
 Les equations de conservation

1. Volume de contrôle non déformable :

𝜕 𝜕
∫ 𝜌 𝑑𝑉 = ∫ (𝜌) 𝑑𝑉
𝜕𝑡 𝑐𝑣 𝑐𝑣 𝜕𝑡

2. Volume de contrôle stationnaire :

⃗⃗⃗⃗
𝑉𝑅 =⃗⃗⃗𝑉

𝜕
∫ ⃗ . 𝜂̂ ) 𝑑𝐴 = 0
𝜌 𝑑𝑉 + ∫ (𝜌. 𝑉
𝑐𝑣 𝜕𝑡 𝑐𝑠

On applique le théorème de divergence de Gauss :

⃗ . 𝜂̂ ) 𝑑𝐴 = ∫ ∇(𝜌. 𝑉
∫ (𝜌. 𝑉 ⃗ ) 𝑑𝑉
𝑐𝑠 𝑐𝑣

𝜕𝜌
[ ⃗ )] 𝑑𝑉 = 0
+ ∇(𝜌. 𝑉
𝜕𝑡

III.1.3 Forme différentielle de l’équation de continuité :

𝜕𝜌
⃗)=0
+ ∇(𝜌. 𝑉
𝜕𝑡

Pour un volume de control :

𝜕
⃗ 𝑛⃗ 𝑑𝐴 = 0
∫ 𝜌 𝑑𝑉 + ∫ 𝜌 𝑉
𝜕𝑡 𝑐𝑣 𝑐𝑠

On considère un petit élément : 𝛿𝑣 = 𝛿𝑥𝛿𝑦𝛿𝑧

𝜕 𝜕𝜌
∫ 𝜌 𝑑𝑉 = 𝛿𝑥 𝛿𝑦 𝛿𝑧
𝜕𝑡 𝑐𝑣 𝜕𝑡

𝜌: Uniforme en 𝛿𝑣 (élément petit)

28
 Les equations de conservation

𝑧
𝑥

𝑦
𝛿𝑧

𝜕 𝛿𝑥 𝜕 𝛿𝑥
𝝆 [𝜌𝑢 + 𝜌𝑢 ] 𝛿𝑦 𝛿𝑧
[𝜌𝑢 − 𝜌𝑢 ] 𝛿𝑦 𝛿𝑧 𝜕𝑥 2
𝜕𝑥 2

𝛿𝑦

𝛿𝑥

Fig (III.2) Débit entrant et sortant pour un petit élément 𝑑𝑣 .

Débit massique :
𝜌𝑉
𝑞𝑚 = [𝑘𝑔/𝑠]
𝑡

Débit volumique :
𝑉
𝑞𝑣 = [𝑚3 /𝑠]
𝑡
𝑞𝑚 = 𝜌 𝑞𝑣

Débit massique net :

Sur le plan X :
𝜕(𝜌𝑢)
𝛿𝑥𝛿𝑦𝛿𝑧
𝜕𝑥
Sur le plan Y :
𝜕(𝜌𝑣)
𝛿𝑥𝛿𝑦𝛿𝑧
𝜕𝑦

Sur le plan Z :
𝜕(𝜌𝑤)
𝛿𝑥𝛿𝑦𝛿𝑧
𝜕𝑧

29
 Les equations de conservation

𝜕𝜌 𝜕(𝜌𝑢) 𝜕(𝜌𝑣) 𝜕(𝜌𝑤)

+ + + =0
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧

III.1.4 Forme vectorielle de l’équation de continuité :

𝜕𝜌
⃗ =0
+ ∇. 𝜌𝑉
𝜕𝑡

Stationnaire :
𝜕𝜌
=0
𝜕𝑡

Incompressible :
𝜕(𝜌𝑢) 𝜕𝑢
→𝜌
𝜕𝑥 𝜕𝑥
𝜕𝑢 𝜕𝑣 𝜕𝑤
+ + =0
𝜕𝑥 𝜕𝑦 𝜕𝑧

Application :
Vérifier l’équation de continuité pour les cas suivants :
𝑢 = 2𝑥 2 − 𝑥𝑦 + 𝑧 2 , 𝑣 = 𝑥 2 − 4𝑥𝑦 + 𝑦 2 , 𝑤 = −2𝑥𝑦 − 𝑦𝑧 + 𝑦 2
𝑢 = 4𝑥𝑦 + 𝑦 2 , 𝑣 = 4𝑥𝑦 + 3𝑥

Réponse :
𝜕𝑢 𝜕𝑣 𝜕𝑤
1. + + =0 1. L’équation de continuité est vérifiée.
𝜕𝑥 𝜕𝑦 𝜕𝑧
𝜕𝑢 𝜕𝑣 𝜕𝑤
2. + 𝜕𝑦 + = 4𝑦 + 4𝑥 2. L’équation de continuité n’est pas vérifiée.
𝜕𝑥 𝜕𝑧

30
 Les equations de conservation

III.2 Conservation de la quantité de mouvement (Equations

d’Euler) :

Le Théorème de Transport de Reynolds :

𝑁: Propriété.
𝑛: 𝑁 par unité de masse.

𝑑𝑁 𝜕𝑁
| = ⃗⃗⃗ 𝜂̂ ) 𝑑𝐴
| + ∫ 𝑛. 𝜌(𝑉.
𝑑𝑡 𝑆𝑦𝑠 𝜕𝑡 𝑐𝑣 𝑐𝑠

𝜂: Vecteur normal sur la surface

⃗⃗⃗
𝑉. 𝜂̂ . 𝑑𝐴: Débit volumique
⃗⃗⃗ 𝜂̂ . 𝑑𝐴: Débit massique
𝜌. 𝑉.
𝑛. 𝜌. ⃗⃗⃗
𝑉. 𝜂̂ . 𝑑𝐴: Transport total de la propriété

Fig (III.3) Vecteur normal et vitesse sur

une surface élémentaire 𝑑𝐴.

Pour un système : (conservation de la quantité de mouvement) :

𝑑
∑ ⃗⃗⃗⃗
𝐹𝑁 = ⃗ 𝑑𝑚
∫𝑉
𝑑𝑡

𝑁 = 𝑚𝑉⃗
𝑛=𝑉 ⃗
𝑑(𝑚𝑉⃗) 𝜕
| = ∫ 𝜌𝑉⃗ 𝑑𝑉 + ∫ (𝜌𝑉
⃗ )(𝑉
⃗⃗⃗⃗⃗𝑅 . 𝜂̂ ) 𝑑𝐴
𝑑𝑡 𝑆𝑦𝑠 𝜕𝑡 𝑐𝑣 𝑐𝑠

1. Volume de contrôle non déformable :

𝜕 𝜕
∫ 𝜌𝑉⃗⃗⃗⃗𝑅 𝑑𝑉 = ∫ ⃗⃗⃗⃗𝑅 ) 𝑑𝑉
(𝜌𝑉
𝜕𝑡 𝑐𝑣 𝑐𝑣 𝜕𝑡

2. Volume de contrôle stationnaire :

⃗⃗⃗⃗
𝑉𝑅 =⃗⃗⃗𝑉

31
 Les equations de conservation

𝑑(𝑚𝑉⃗) 𝜕
| =∫ ⃗ ) 𝑑𝑉 + ∫ (𝜌𝑉
(𝜌𝑉 ⃗ )(⃗⃗⃗𝑉 . 𝜂̂ ) 𝑑𝐴
𝑑𝑡 𝑆𝑦𝑠 𝑐𝑣 𝜕𝑡 𝑐𝑠

On applique le théorème de divergence de Gauss :

𝜕
∑ ⃗⃗⃗⃗
𝐹𝑁 = ∫ ⃗ ) 𝑑𝑉 + ∫ (𝜌𝑉
(𝜌𝑉 ⃗ )(⃗⃗⃗𝑉 . 𝜂̂ ) 𝑑𝐴
𝑐𝑣 𝜕𝑡 𝑐𝑠

⃗ )(⃗⃗⃗𝑉 . 𝜂̂ ) 𝑑𝐴 = ∫ ∇(𝜌. 𝑉
∫ (𝜌𝑉 ⃗)𝑉
⃗ 𝑑𝑉
𝑐𝑠 𝑐𝑣

𝜕
∑ ⃗⃗⃗⃗
𝐹𝑁 = ∫ ⃗ ) 𝑑𝑉 + ∫ ∇(𝜌. 𝑉
(𝜌𝑉 ⃗ )⃗⃗⃗𝑉 𝑑𝑉
𝑐𝑣 𝜕𝑡 𝑐𝑣

On considère un petit élément : 𝛿𝑣 = 𝛿𝑥𝛿𝑦𝛿𝑧

𝑑 ⃗
𝑑𝑉
⃗⃗⃗⃗𝑁 =
𝛿𝐹 ⃗ 𝛿𝑚) = 𝛿𝑚
(𝑉
𝑑𝑡 𝑑𝑡
⃗⃗⃗⃗𝑁 = 𝛿𝑚 𝑎
𝛿𝐹

III.2.1 La 2 ème loi de Newton :

∑ ⃗⃗⃗⃗
𝐹𝑁 =Masse x Accélération

Forces extérieures = Forces d’inertie

∑ ⃗⃗⃗⃗
𝐹𝑁 = 𝐹𝑆𝑢𝑟𝑓𝑎𝑐𝑒 + 𝐹𝑣𝑜𝑙𝑢𝑚𝑒

Forces extérieures :
Poids (Body forces) :
Ex : gravité, champs magnétique, champs électrique…
Forces de surfaces (Surface forces), elles influent sur la surface:
Ex : Force normale : 𝑃 (pressions) et les forces tangentielles (forces de viscosité) 𝜏.

32
 Les equations de conservation

𝝈𝒙𝒙 : 𝐶𝑜𝑛𝑡𝑟𝑎𝑖𝑛𝑡𝑒 𝑛𝑜𝑟𝑚𝑎𝑙𝑒

𝑧
𝑥 𝝉𝒙𝒚 , 𝝉𝒙𝒛 : 𝐶𝑜𝑛𝑡𝑟𝑎𝑖𝑛𝑡𝑒𝑠 𝑑𝑒 𝑐𝑖𝑠𝑎𝑖𝑖𝑙𝑒𝑚𝑒𝑛𝑡

𝝉𝒙𝒛
𝑦 𝝈𝒙𝒙
𝛿𝑧 𝝉𝒙𝒚

𝛿𝑦
𝛿𝑥

Fig (III.4) La Contrainte normale et les contraintes tangentielles sur une surface élémentaire 𝑑𝐴.

Sur le plan X :

𝑧
𝑥

𝜕 𝛿𝑥 𝜕 𝛿𝑥
[𝝈𝒙𝒙 − 𝝈𝒙𝒙 ] 𝛿𝑦 𝛿𝑧 [𝝈𝒙𝒙 + 𝝈𝒙𝒙 ] 𝛿𝑦 𝛿𝑧
𝜕𝑥 2 𝜕𝑥 2

𝛿𝑥

𝑧
𝑥

𝜕 𝛿𝑦
[𝝉𝒚𝒙 − 𝝉 ] 𝛿𝑥 𝛿𝑧
𝑦 𝜕𝑥 𝒚𝒙 2

𝜕 𝛿𝑦
[𝝉𝒚𝒙 + 𝝉𝒚𝒙 ] 𝛿𝑥 𝛿𝑧
𝜕𝑥 2

𝛿𝑥

33
 Les equations de conservation

𝑧
𝑥

𝜕 𝛿𝑧
[𝝉𝒛𝒙 + 𝝉 ] 𝛿𝑥 𝛿𝑦
𝑦 𝜕𝑥 𝒛𝒙 2

𝜕 𝛿𝑧
[𝝉𝒛𝒙 − 𝝉 ] 𝛿𝑥 𝛿𝑦
𝜕𝑥 𝒛𝒙 2

𝛿𝑥

Fig (III.5) La Contrainte normale et les contraintes tangentielles sur le plan X.

Forces de surfaces :

Sur le plan X :

𝜕𝜎𝑥𝑥 𝜕𝜏𝑦𝑥 𝜕𝜏𝑧𝑥

𝛿𝐹𝑠𝑥 = ( + + ) 𝛿𝑥 𝛿𝑦𝛿𝑧
𝜕𝑥 𝜕𝑦 𝜕𝑧

Sur le plan Y :
𝜕𝜏𝑥𝑦 𝜕𝜎𝑦𝑦 𝜕𝜏𝑧𝑦
𝛿𝐹𝑠𝑦 = ( + + ) 𝛿𝑥 𝛿𝑦𝛿𝑧
𝜕𝑥 𝜕𝑦 𝜕𝑧

Sur le plan Z :

𝜕𝜏𝑥𝑧 𝜕𝜏𝑦𝑧 𝜕𝜎𝑧𝑧

𝛿𝐹𝑠𝑧 = ( + + ) 𝛿𝑥 𝛿𝑦𝛿𝑧
𝜕𝑥 𝜕𝑦 𝜕𝑧

𝛿𝐹 = 𝛿𝑚 𝑎

𝛿𝑚 = 𝜌 𝛿𝑥 𝛿𝑦𝛿𝑧

⃗
𝑑𝑉
∑ 𝐹 = 𝑚𝑎𝑠𝑠𝑒 . 𝑎𝑐𝑐𝑒𝑙𝑒𝑟𝑎𝑡𝑖𝑜𝑛 = 𝑚
𝑑𝑡

34
 Les equations de conservation

Les équations de mouvement s’écrivent :

𝜕𝜎𝑥𝑥 𝜕𝜏𝑦𝑥 𝜕𝜏𝑧𝑥 𝜕𝑢 𝜕𝑢 𝜕𝑢 𝜕𝑢
𝜌𝑔𝑥 + + + = 𝜌( +𝑢 +𝑣 +𝑤 )
𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧
𝜕𝜏𝑥𝑦 𝜕𝜎𝑦𝑦 𝜕𝜏𝑧𝑦 𝜕𝑣 𝜕𝑣 𝜕𝑣 𝜕𝑣
𝜌𝑔𝑦 + + + = 𝜌( +𝑢 +𝑣 +𝑤 )
𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧
𝜕𝜏𝑥𝑧 𝜕𝜏𝑦𝑧 𝜕𝜎𝑧𝑧 𝜕𝑤 𝜕𝑤 𝜕𝑤 𝜕𝑤
𝜌𝑔𝑧 + + + = 𝜌( +𝑢 +𝑣 +𝑤 )
{ 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧

Bilan :
Nombre d’équations : 03 équations+ 01 equation de continuité= 04 équations
Inconnus : 𝑢, 𝑣, 𝑤, + les contraintes.

On a :

𝜎𝑖𝑗 = −𝑃𝛿𝑖𝑗 + 𝜏𝑖𝑗

Ecoulements parfaits : Contraintes de cisaillement sont nuls.

𝜏𝑖𝑗= 0

𝜎𝑥𝑥 = 𝜎𝑦𝑦 = 𝜎𝑧𝑧 = −𝑃

𝜕𝑃 𝜕𝑢 𝜕𝑢 𝜕𝑢 𝜕𝑢
𝜌𝑔𝑥 − = 𝜌( +𝑢 +𝑣 +𝑤 )
𝜕𝑥 𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧
𝜕𝑃 𝜕𝑣 𝜕𝑣 𝜕𝑣 𝜕𝑣
𝜌𝑔𝑦 − = 𝜌( +𝑢 +𝑣 +𝑤 )
𝜕𝑦 𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧
𝜕𝑃 𝜕𝑤 𝜕𝑤 𝜕𝑤 𝜕𝑤
𝜌𝑔𝑧 − = 𝜌( +𝑢 +𝑣 +𝑤 )
{ 𝜕𝑧 𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧

Les équations ci-dessus sont appelées les équations d’Euler.

III.2.2 Forme vectorielle de l’équation d’Euler:

⃗
𝜕𝑉
⃗ ∇⃗⃗⃗𝑉 )
𝜌𝑔 + ∇𝑃 = 𝜌 ( + 𝑉
𝜕𝑡

35
 Les equations de conservation

III.3 conservation de la quantité de mouvement (Equations de

Navier Stokes) :
III.3.1 Le tenseur de contraintes total :
𝜎𝑖𝑗 = −𝑃𝛿𝑖𝑗 + 𝜏𝑖𝑗

𝑃: La contrainte normale (pression hydrostatique).

𝜏𝑖𝑗 : Le tenseur de contraintes visqueuses (contraintes de cisaillements).
𝛿𝑖𝑗 : Le symbole de Kronecker.

𝛿𝑖𝑗 = 1 𝑠𝑖 𝑖 = 𝑗
{
𝛿𝑖𝑗 = 0 𝑠𝑖 𝑖 ≠ 𝑗

𝜕𝑢𝑖 𝜕𝑢𝑗 𝜕𝑢𝑘

𝜏𝑖𝑗 = 2𝜇𝜀𝑖𝑗 + 𝜂𝛿𝑖𝑗 𝑒 = 𝜇 ( + ) + 𝜂𝛿𝑖𝑗
𝜕𝑥𝑗 𝜕𝑥𝑖 𝜕𝑥𝑘

𝜇 et 𝜂 sont deux paramètres caractérisant les propriétés visqueuses (le frottement

interne) des fluides.

𝜇: Viscosité de cisaillement
𝜂: Viscosité de dilatation
𝜀𝑖𝑗 : Le tenseur de taux de déformation
𝑒: Le taux de dilatation

𝜕𝑢 𝜕𝑣 𝜕𝑤
⃗ =
𝑒 = 𝑑𝑖𝑣 𝑉 + +
𝜕𝑥 𝜕𝑦 𝜕𝑧

III.3.2 Les contraintes normales et les contraintes tangentielles:

Les contraintes normales :

𝜕𝑢 𝜕𝑢 𝜕𝑣 𝜕𝑤
𝜎𝑥𝑥 = −𝑃 + 2𝜇 +𝜂( + + )
𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑧
𝜕𝑣 𝜕𝑢 𝜕𝑣 𝜕𝑤
𝜎𝑖𝑗 = −𝑃𝛿𝑖𝑗 + 𝜏𝑖𝑗 avec 𝑖 = 𝑗 𝜎𝑦𝑦 = −𝑃 + 2𝜇 +𝜂( + + )
𝜕𝑦 𝜕𝑥 𝜕𝑦 𝜕𝑧
𝜕𝑤 𝜕𝑢 𝜕𝑣 𝜕𝑤
𝜎𝑧𝑧 = −𝑃 + 2𝜇 +𝜂( + + )
{ 𝜕𝑧 𝜕𝑥 𝜕𝑦 𝜕𝑧

36
 Les equations de conservation

Les contraintes tangentielles :

𝜕𝑢 𝜕𝑣
𝜎𝑥𝑦 = 𝜎𝑦𝑥 = 𝜏𝑥𝑦 = 𝜏𝑦𝑥 = 𝜇 ( + )
𝜕𝑦 𝜕𝑥
𝜕𝑣 𝜕𝑤
𝜎𝑖𝑗 = 𝜏𝑖𝑗 avec 𝑖 ≠ 𝑗 𝜎𝑦𝑧 = 𝜎𝑧𝑦 = 𝜏𝑦𝑧 = 𝜏𝑧𝑦 = 𝜇 ( + )
𝜕𝑧 𝜕𝑦
𝜕𝑤 𝜕𝑢
{ 𝜎𝑧𝑥 = 𝜎𝑥𝑧 = 𝜏𝑧𝑥 = 𝜏𝑥𝑧 = 𝜇 ( 𝜕𝑥 + 𝜕𝑧 )

Les forces de surface :

Sur le plan X :

𝜕𝜎𝑥𝑥 𝜕𝜏𝑦𝑥 𝜕𝜏𝑧𝑥

𝛿𝐹𝑠𝑥 = ( + + ) 𝛿𝑥 𝛿𝑦𝛿𝑧
𝜕𝑥 𝜕𝑦 𝜕𝑧
𝜕𝑃 𝜕 2𝑢 𝜕 2𝑢 𝜕 2𝑣 𝜕 2𝑤 𝜕 2𝑢 𝜕 2𝑣
= [(− + 2𝜇 2 + 𝜂 ( 2 + + )+𝜇( 2 + )
𝜕𝑥 𝜕𝑥 𝜕𝑥 𝜕𝑥𝜕𝑦 𝜕𝑥𝜕𝑧 𝜕𝑦 𝜕𝑦𝜕𝑥
𝜕 2𝑤 𝜕 2𝑢
+𝜇( + ))] 𝛿𝑥 𝛿𝑦𝛿𝑧
𝜕𝑧𝜕𝑥 𝜕𝑧 2

𝜕𝑃 𝜕 2𝑢 𝜕 2𝑢 𝜕 2𝑢 𝜕 2𝑢 𝜕 2𝑣 𝜕 2𝑤
𝛿𝐹𝑠𝑥 = [− + 𝜇 ( 2 + 2 + 2) + 𝜇 ( 2 + + )
𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑥 𝜕𝑥𝜕𝑦 𝜕𝑧𝜕𝑥
𝜕 2𝑢 𝜕 2𝑣 𝜕 2𝑤
+𝜂( 2 + + )] 𝛿𝑥 𝛿𝑦𝛿𝑧
𝜕𝑥 𝜕𝑥𝜕𝑦 𝜕𝑥𝜕𝑧

𝜕𝑃 𝜕 2𝑢 𝜕 2𝑢 𝜕 2𝑢 𝜕 𝜕𝑢 𝜕𝑣 𝜕𝑤
𝛿𝐹𝑠𝑥 = [− + 𝜇 ( 2 + 2 + 2) + 𝜇 ( + + )
𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑧
𝜕 𝜕𝑢 𝜕𝑣 𝜕𝑤
+𝜂 ( + + )] 𝛿𝑥 𝛿𝑦𝛿𝑧
𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑧

𝜕𝑃 𝜕 2𝑢 𝜕 2𝑢 𝜕 2𝑢 𝜕 𝜕𝑢 𝜕𝑣 𝜕𝑤
𝛿𝐹𝑠𝑥 = [− + 𝜇 ( 2 + 2 + 2 ) + (𝜇 + 𝜂) ( + + )] 𝛿𝑥 𝛿𝑦𝛿𝑧
𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑧

𝜕𝑃 𝜕 2𝑢 𝜕 2𝑢 𝜕 2𝑢 𝜕𝑒
𝛿𝐹𝑠𝑥 = [− + 𝜇 ( 2 + 2 + 2 ) + (𝜇 + 𝜂) ] 𝛿𝑥 𝛿𝑦𝛿𝑧
𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑥

37
 Les equations de conservation

Incompressibilité :

⃗ =0
𝑒 = 𝑑𝑖𝑣 𝑉

𝜕𝑢 𝜕𝑣 𝜕𝑤
+ + =0
𝜕𝑥 𝜕𝑦 𝜕𝑧

Les équations de mouvement s’écrivent :

𝜕𝑢 𝜕𝑢 𝜕𝑢 𝜕𝑢 𝜕𝑃 𝜕 2𝑢 𝜕 2𝑢 𝜕 2𝑢
𝜌( +𝑢 +𝑣 + 𝑤 ) = 𝜌𝑔𝑥 − + 𝜇 ( 2 + 2 + 2)
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑧
𝜕𝑣 𝜕𝑣 𝜕𝑣 𝜕𝑣 𝜕𝑃 𝜕 𝑣 𝜕 𝑣 𝜕 2𝑣
2 2
𝜌( +𝑢 +𝑣 + 𝑤 ) = 𝜌𝑔𝑦 − + 𝜇 ( 2 + 2 + 2)
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑦 𝜕𝑥 𝜕𝑦 𝜕𝑧
𝜕𝑤 𝜕𝑤 𝜕𝑤 𝜕𝑤 𝜕𝑃 𝜕 𝑤 𝜕 𝑤 𝜕 2𝑤
2 2
𝜌( +𝑢 +𝑣 +𝑤 ) = 𝜌𝑔𝑧 − +𝜇( 2 + + )
{ 𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑧 𝜕𝑥 𝜕𝑦 2 𝜕𝑧 2

Les équations ci-dessus sont appelées les équations de Navier Stokes pour les fluides
incompressibles.

III.3.3 Forme vectorielle de l’équation de Navier Stokes:

⃗
𝜕𝑉
𝜌( ⃗ ∇⃗⃗⃗𝑉) = 𝜌𝑔 + ∇𝑃 + 𝜇 ∇2 ⃗⃗⃗𝑉
+𝑉
𝜕𝑡

⃗ ∇⃗⃗⃗𝑉 : Terme convectif.

𝑉
𝜕𝑉⃗
: Terme instationnaire.
𝜕𝑡
𝜌𝑔: Gravité (force de volume)
∇𝑃: Pression (force de surface).
𝜇 ∇2 ⃗⃗⃗𝑉 : Force de viscosité (force de surface).

Avec :

2
𝜕2 𝜕2 𝜕2
∇ = 2+ 2+ 2
𝜕𝑥 𝜕𝑦 𝜕𝑧

Bilan :
Nombre d’équations : 03 équations+ 01 equation de continuité= 04 équations.
Inconnus : 𝑢, 𝑣, 𝑤, +𝑃 𝑒𝑡 𝜌 = 05 inconnus.

38
 Les equations de conservation

On ajoute l’équation d’état.

Problème :

 Le terme convectif est non linéaire.

 Equation différentielle de 2 ème ordre.
 Très peu de solutions (cas d’un écoulement stationnaire incompressible
laminaire d’un fluide visqueux newtonien).

La seule approche utilisée pour résoudre ces équations c’est la CFD (Computational
Fluid Dynamics).

39
 Les equations de conservation

III.4 Exercices avec réponses:

Exercice 01:
Pour un écoulement incompressible plan, lequel de ces champs de vitesse suivants
satisfait l’équation de conservation de masse :

(a) 𝑢 = −𝑥 , 𝑣 = 𝑦 (b) 𝑢 = 3𝑦 , 𝑣 = 3𝑥 (c) 𝑢 = 4𝑥 , 𝑣 = −4𝑦

(d) 𝑢 = 3𝑥𝑡 , 𝑣 = 3𝑦𝑡 (e) 𝑢 = 𝑥𝑦 + 𝑦 2 𝑡 , 𝑣 = 𝑥𝑦 + 𝑥 4 𝑡 (f) 𝑢 = 4𝑥 2 𝑦 3 , 𝑣 = −2𝑥𝑦 4

Les composantes de vitesse dans les directions 𝑥 et 𝑧 d’un écoulement

incompressible en régime stationnaire sont données comme suit :

𝑢 = 𝑥 4 +3𝑧 4 , 𝑤 = 𝑦 4 − 3𝑦𝑧

Donner la forme générale de la deuxième composante de vitesse qui satisfait

l’équation de continuité.
Le champ de vitesse d’un écoulement incompressible est donné par :

𝑉 = 𝑥 2 𝑖 − 𝑧 2 𝑗 − 3𝑥𝑧𝑘.

Si la viscosité du fluide est de 0,04 𝑃𝑎. 𝑠, évaluer le tenseur des contraintes

visqueuses 𝜏𝑖𝑗 au point (𝑥, 𝑦, 𝑧) = (3,2,1).

Réponse :
1. L’équation de continuité est satisfaite pour (a), (b), (c), et (f), et elle n’est pas
satisfaite pour (d) et (e).
3
2. 𝑣 = 2 𝑦 2 − 4𝑦𝑥 3 + 𝑓(𝑥, 𝑦)
0,48 0 −0,12
3. 𝜏𝑖𝑗 = | 0 0 −0,08|
−0,12 −0,08 −0,72

Exercice 02:
On appliquant l’équation de conservation de masse et la quantité du mouvement
sur un écoulement 2𝐷 suivant 𝑋 entre deux plaques parallèles inclinées (voir
annexe III.1)

Calculer :
La vitesse 𝑢 ;
La contrainte de Cisaillement 𝜏 ;
Le débit 𝑄.
Sachant que le fluide est incompressible, newtonien et laminaire et pour un
régime stationnaire.

40
 Les equations de conservation

Réponse :
1 𝑃2 −𝑃1
1. 𝑢 = 2𝜇 . [− ( ) + 𝜌𝑔. 𝑠𝑖𝑛𝜃] (ℎ2 − 𝑦 2 )
𝐿
𝜕𝑃
2. 𝜏(𝑦=0) = 0 , 𝜏(𝑦=ℎ) = − [− (𝜕𝑥 ) + 𝜌𝑔𝑥 ] ℎ
3 ∆𝑃
2ℎ [− +𝜌𝑔.𝑠𝑖𝑛𝜃]
𝐿
3. 𝑄 = 3𝜇

Exercice 03:
Un champ de vitesse à deux dimensions est donné par :

𝐾𝑦
𝑢=
𝑥2 + 𝑦2

𝐾𝑥
𝑣=
𝑥2 + 𝑦2

𝐾 est une constante. Le champ satisfait-il l’équation de continuité? Transformer

ces vitesses en composantes polaires 𝑣𝑟 et 𝑣𝜃 sachant que :

𝑥2 + 𝑦2 = 𝑟2
𝑥
cos 𝜃 =
𝑟
𝑦
sin 𝜃 =
𝑟

Réponse :
𝜕𝑢 𝜕𝑣 𝜕𝑤
1. + 𝜕𝑦 + =0 L’équation de continuité est satisfaite.
𝜕𝑥 𝜕𝑧
𝐾
2. 𝑣𝑟 = 0 , 𝑣𝜃 = 𝑟

41
Chapitre IV

LA COUCHE LIMITE

IV LA COUCHE LIMITE ……………………………………………………………….…. 42

IV.1 Introduction .………………………………………………….………………….……. 43
IV.2 Grandeurs caractéristiques de la couche limite ………………………….. 44
IV.2.1 Epaisseur ………………………………………………………………………… 44
IV.2.2 Epaisseur de déplacement ………………………………………………… 45
IV.2.3 Epaisseur de la quantité de mouvement …………………………….. 46
IV.2.4 Autres facteurs et coefficients …………………………………………… 46
IV.3 Evolution de la couche limite ……………………………………….…………… 47
IV.3.1 Plaque plane …………………………………………………………………… 47
IV.3.2 Une forme arrondie (profil, cylindre….) …………………………….. 48
IV.4 Formules utilisées pour le calcul des paramètres de la couche
limite (cas d’une plaque plane) ……………………..………………………….. 50
IV.5 Exercices avec réponses ……………………………………………….…………… 45

42
 La couche limite

IV.1 Introduction :
Considérons l’écoulement d’un fluide autour d’un obstacle, les effets de viscosité sont
négligeables loin, de la paroi, le fluide peut être considéré parfait.
Lorsqu’on se rapproche da la paroi, les effets de la viscosité sont considérables et la
vitesse du fluide devient progressivement nulle.
Le domaine de transition où la vitesse devient progressivement nulle est appelé
couche limite comme présente la figure IV.1.

Fluide
parfait

Fluide
visqueux

Fig (IV.1) Schéma descriptif de la couche limite.

L’épaisseur de la couche limite dépend de Reynolds 𝑅𝑒 . Au sein de la couche limite,

l’écoulement peut être soit laminaire soit turbulent (cela dépend de 𝑅𝑒 aussi).

Fig (IV.2) Formation d’une couche limite sur une plaque plane.
(a) écoulement laminaire (b) écoulement transitoire (c) écoulement turbulent

43
 La couche limite

Au contact de la paroi, quand le profil de vitesse présente une pente infinie, on dit
qu’il y a un décollement de la couche limite. Après le point de décollement 𝐷, la
couche limite devient turbulente. On cherche toujours à optimiser le profil de façon à
minimiser le sillage. Le sillage est responsable d’une dissipation d’énergie.

Fig (IV.3) Décollement de la couche limite.

IV.2 Grandeurs caractéristiques de la couche limite :

IV.2.1 Epaisseur :
La figure IV.4 montre l’épaisseur de la couche limite qui est définie comme la
distance à la paroi à partir de laquelle la vitesse devient supérieure à 99% de la vitesse
d’écoulement uniforme 𝑢(𝑦 = 𝛿) = 0.99𝑈, c'est-à-dire non perturbé par l’objet.
Quand on se déplace sur une perpendiculaire 𝑂𝑦 à la paroi, la vitesse 𝑢 se rapproche
rapidement d’une manière asymétrique de celle qui règne dans le fluide libre, on peut
considérer l’épaisseur de la couche limite comme étant la distance 𝛿 au-delà de
laquelle 𝑢 = 0.99𝑈.

𝑼∞

𝒖(𝒚 = 𝜹) = 𝟎. 𝟗𝟗𝑼

Fig (IV.4) Epaisseur de la couche limite.

44
 La couche limite

IV.2.2 Epaisseur de déplacement :

Pour définir l’épaisseur de déplacement, on évolue le flux manquant par rapport à
celui qu’on aurait dans l’hypothèse d’un écoulement uniforme jusqu’à la paroi.
Les lignes de courants en écoulement visqueux sont déplacées par rapport à leurs
positions en fluide non visqueux.
L’épaisseur de déplacement décrit le déficit de débit 𝑄𝑃 − 𝑄𝑉 .
Ce déplacement est exploité pour définir une épaisseur tels que les aires 𝐴 𝑒𝑡 𝐴′ soient
égales :

Fig (IV.5) Epaisseur de déplacement.

𝛿
𝑄𝑃 = ∫ 𝑈∞ 𝑑𝑦
0

𝛿
𝑄𝑉 = ∫ 𝑢𝑑𝑦
0

𝑄𝑃 : Débit volumique pour un fluide parfait.

𝑄𝑉 : Débit volumique pour un fluide visqueux.

𝛿
𝐴 = ∫ (𝑈∞ − 𝑢)𝑑𝑦 = 𝐴′ = 𝛿 ∗ 𝑈∞
0
Avec :

𝛿
𝑢
𝛿 ∗ = ∫ (1 − ) 𝑑𝑦
0 𝑈∞

L’épaisseur de déplacement 𝛿 ∗ représente le déficit en débit de fluide dû à la présence

de la couche Limite autour du corps.

45
 La couche limite

IV.2.3 Epaisseur de la quantité de mouvement :

On multiplie les débits par la vitesse 𝑢 et de la même manière on définit l’épaisseur de
la quantité de mouvement 𝜃.
𝛿 𝛿
∫ 𝑢. 𝑈∞ 𝑑𝑦 − ∫ 𝑢. 𝑢𝑑𝑦 = 𝜃. 𝑈∞ 2
0 0

𝛿
𝑢 𝑢
𝜃=∫ (1 − ) 𝑑𝑦
0 𝑈∞ 𝑈∞

L’épaisseur de Quantité de Mouvement 𝜃 représente le déficit en quantité de

mouvement dans la couche limite par rapport à ce qu’on aurait dans un écoulement
de fluide parfait.
Une autre définition est :
L’épaisseur de quantité de mouvement mesure la perte de quantité de mouvement
due à la friction des particules de fluide à la surface du corps.

Région laminaire Région

Région turbulente
transitoire

Fig (IV.6) Epaisseur de quantité de mouvement.

IV.2.4 Autres facteurs et coefficients :

 Facteur de forme de la couche limite :
𝛿∗
𝐻=
𝜃

 Contrainte tangentielle à la paroi :

𝜕𝑢
𝜏0 = 𝜇 ( )
𝜕𝑦

46
 La couche limite

 Coefficient de frottement local :

𝜏0
𝐶𝑓 =
1⁄2𝜌𝑈 2

 Coefficient de frottement global sur une surface S :

1
𝐶𝐹 = ∫ 𝐶𝑓 𝑑𝑆
𝑆

 Le flux d’énergie dissipée :

∞
𝜕𝑢 2
𝐹𝐷 = ∫ 𝜇 ( ) 𝑑𝑦
0 𝜕𝑦

 Equation de la couche limite (forme intégrale) :

𝑑𝜃 𝜃 𝑑𝑈∞ 𝐶𝑓
+ (𝐻 + 2) =
𝑑𝑥 𝑈∞ 𝑑𝑥 2

IV.3 Evolution de la couche limite :

IV.3.1 Plaque plane :
Quand un fluide s’écoule autour d’une plaque plane, l’action des contraintes de
cisaillement (forces de viscosité) retarde l’écoulement dans la direction latérale
(direction de l’écoulement), l’épaisseur de la couche limite augmente.
Dues aux contraintes de cisaillement, l’écoulement dans la couche limite est
rotationnel et il est irrotationnel en dehors de la couche limite (Prandtl).

Région
Région laminaire Région turbulente
transitoire

Fig (IV.7) Structure de l’écoulement dans la couche limite.

𝜕𝑢
𝜏=𝜇
𝜕𝑦

La valeur de 𝜏 en dehors de la couche est égale à 0 car :

47
 La couche limite

𝜕𝑢
=0
𝜕𝑦

𝜕𝑢
0≪𝑦≪𝛿 : 𝜕𝑦
est très élevé, l’effet de la viscosité est dominant.

𝑦>𝛿 : 𝜏 = 0 L’effet de viscosité n’est pas dominant.

La Formule de Blasius pour une couche limite laminaire :

𝛿 5
=
𝑥 √𝑅𝑒

IV.3.2 Une forme arrondie (profil, cylindre….) :

Considérant un écoulement autour d’un profil ou un obstacle avec une forme arrondi
(ex : un cylindre), il est clairement montré que le profil de vitesse change
progressivement due aux forces de viscosité et due aussi à la forme géométrique. La
figure IV.8 montre l’évolution de la couche limite jusqu’au point de séparation (𝑆) où
il y aura la formation de sillage.

𝝏𝑷 𝝏𝑷 𝝏𝑷 𝝏𝑷
𝝏𝒙
<𝟎 𝝏𝒙
=𝟎 𝝏𝒙
>𝟎 𝝏𝒙
>𝟎
𝝏𝑼 𝝏𝑼 𝝏𝑼 𝝏𝑼
𝝏𝒚
>𝟎 𝝏𝒚
>𝟎 𝝏𝒚
=𝟎 𝝏𝒚
<𝟎

(a) (b) (c) (d)

Fig (IV.8) Evolution du profil de vitesse.

48
 La couche limite

𝜕𝑃
< 0 ⇒ : Accélération dans le sens de l’écoulement ⇒ Gradient de pression
𝜕𝑥
favorable.

𝜕𝑃 𝜕𝑈
>0⇒ <0 : Décélération dans le sens de l’écoulement ⇒ Gradient de
𝜕𝑥 𝜕𝑦
pression défavorable.

Un gradient de pression défavorable peut donc provoquer un décollement au point S ;

Le décollement est accompagné d'une formation de tourbillons (figure IV.9) et peut
avoir de graves conséquences :

 Instabilité de l'écoulement ;
 Transition au régime turbulent ;
 Augmentation de trainée.

Fig (IV.9) Décollement de la couche limite (écoulement autour d’un cylindre).

La nature de l’écoulement pouvait être localement décrite au moyen du nombre de

Reynolds local :
𝑈𝑥
𝑅𝑒𝑥 =
𝜐

Quand 𝑅𝑒𝑥 > 5. 105 la couche limite devient turbulente, et la transition du laminaire
au turbulent se manifeste par un épaississement brutal de la couche limite.

49
 La couche limite

IV.4 Formules utilisées pour le calcul des paramètres de la couche

limite (cas d’une plaque plane) :
Le tableau IV.1 montre les formules utilisées dans le calcul des paramètres
caractérisant la couche limite pour différents régimes (cas d’une plaque plane).

Tab (IV.1) Les formules de Blasius et Prandtl pour le calcul des paramètres de la couche limite
(cas d’une plaque plane).
Paramètres Laminaire Turbulent
(solution de Blasius) (Approximation de Prandtl)
𝛿 5 𝛿 0.16
Epaisseur de la C.L = =
𝑥 √𝑅𝑒𝑥 𝑥 (𝑅𝑒𝑥 )1/7

𝛿∗ 1,72 𝛿∗ 0.02
Epaisseur de déplacement = =
𝑥 √𝑅𝑒𝑥 𝑥 (𝑅𝑒𝑥 )1/7

𝜃 0,664 𝜃 0.016
Epaisseur de quantité de = =
𝑥 √𝑅𝑒𝑥 𝑥 (𝑅𝑒𝑥 )1/7
mouvement

𝛿∗ 𝛿∗
Facteur de forme 𝐻= = 2,59 𝐻= = 1,25
𝜃 𝜃

0,664 0.027
Coefficient de frottement 𝐶𝑓 = 𝐶𝑓 =
√𝑅𝑒𝑥 (𝑅𝑒𝑥 )1/7
local

0,332. 𝜇 1⁄2 . 𝜌1⁄2 . 𝑈 3⁄2 0,0135. 𝜇 1⁄7 . 𝜌6⁄7 . 𝑈 13⁄7

Contrainte de cisaillement 𝜏0 = 𝜏0 =
sur la paroi 𝑥 1⁄2 𝑥 1⁄7

1,328 0,031
Coefficient de trainée 𝐶𝑑 = 𝐶𝑑 =
√𝑅𝑒𝐿 (𝑅𝑒𝐿 )1/7

50
 La couche limite

IV.5 Exercices avec réponses:

Exercice 01:
∗ ∗∗
Calculer les valeurs des épaisseurs adimensionnelles 𝛿 ⁄𝛿 , 𝛿 ⁄𝛿 ainsi que le
rapport de forme 𝐻 pour le profil de vitesse linéaire :

𝑢 𝑦
=
𝑈 𝛿
Sachant que :
𝑈: Vitesse du fluide libre.
𝑢: Vitesse à la distance y de la paroi.
𝛿: Epaisseur de la couche limite.
𝛿 ∗ : Epaisseur de déplacement.
𝛿 ∗∗ : Epaisseur en énergie.

Réponse :
∗
1. 𝛿 ⁄𝛿 = 1⁄2
∗∗
2. 𝛿 ⁄𝛿 = 1⁄4
3. 𝐻 = 3

Exercice 02 :
Soit un écoulement caractérisé par un profil de vitesse :

𝑢 (𝑥, 𝑦) 𝜋
= 𝑠𝑖𝑛 ( . 𝜂)
𝑈 2

Où :
𝑦
𝜂=
𝛿(𝑥)

∗
1. Exprimer les valeurs : 𝛿 ⁄𝛿 , 𝜃⁄𝛿 , et 𝜏0 . 𝛿
2. Vérifier qu’il s’agit bien toujours de constantes pour cet écoulement.
3. Calculer la loi d’épaisseur de la couche limite 𝛿(𝑥).

Réponse :
∗
1. 𝛿 ⁄𝛿 = 0,363 𝜃⁄ = 0,137 𝜋
𝛿 𝜏0 . 𝛿 = 𝜇. 2 . 𝑈
2. Il s’agit des constantes.
𝜈
3. 𝛿(𝑥) = 4,79√𝑈 . 𝑥

Exercice 03 :
La distribution de vitesse autour d’une plaque plane pour une couche limite
laminaire est donnée comme suit :

51
 La couche limite

𝑢
= 𝑎0 + 𝑎1 𝑦 + 𝑎2 𝑦 2 + 𝑎3 𝑦 3
𝑢1

1. Calculer l’épaisseur de quantité de mouvement 𝜃 ?

Réponse :
1. 𝜃 = 39𝛿⁄280

Exercice 04 :
Soit un écoulement laminaire autour d’une plaque mince et plate, on admet que le
profil de vitesse dans la couche limite répond à l’équation :

𝑢 3 𝑦 1 𝑦 3
= . − ( )
𝑈∞ 2 𝛿 2 𝛿

𝑈∞ : Vitesse d’écoulement libre.

𝑢: Vitesse à une distance 𝑦 de la paroi.
𝛿: Epaisseur de la couche limite.
𝜃: Epaisseur de quantité de mouvement : 𝜃 = 0.139𝛿

Déterminer :
1. Epaisseur de déplacement 𝛿 ∗ .
2. Facteur de forme 𝐻.
3. Contrainte de cisaillement 𝜏𝑚𝑎𝑥 .
4. Coefficient de frottement local 𝐶𝑓 en fonction de 𝛿.

Réponse :
1. 𝛿 ∗ = 0,375
2. 𝐻 = 2,7
3 𝑈∞
3. 𝜏𝑚𝑎𝑥 = 2 . 𝜇. 𝛿
3𝜇
4. 𝐶𝑓 =
𝛿.𝜌.𝑈∞

52
Chapitre V

ECOULEMENTS AUTOUR
DES OBSTACLES

V ECOULEMENTS AUTOUR DES OBSTACLES …………………………… 53

V.1 Introduction .………………………………………………….………………….……. 54
V.2 Evolution expérimentale de 𝐶𝑥 (𝑅𝑒 ) et description de l’écoulement
autour d’une sphère lisse ………………………………………………………….. 56
V.2.1 Evolution expérimentale de 𝐶𝑥 (𝑅𝑒 ) ……………………………………. 56
V.2.2 Description de l’écoulement ……………………………………………… 57
V.3 Evolution expérimentale de 𝐶𝑥 (𝑅𝑒 ) et description de l’écoulement
autour d’un cylindre ………………………………………..……………………….. 59
V.3.1 Evolution expérimentale de 𝐶𝑥 (𝑅𝑒 ) …………………………………… 59
V.3.2 Description de l’écoulement ……………………………………………… 60
V.4 Ecoulement autour d’un hydrofoil ………………..…………………………… 61
V.5 Ecoulement autour d’une hélice ………………………………………………... 62
V.6 Ecoulement autour d’une plaque plane (calcul du 𝐶𝑥 ) ………………….. 63
V.7 Quelques paramètres influents sur le coefficient de
Trainée 𝐶𝑥 ……………………………………………………………………………….. 65
V.7.1 Position de la plaque ………………………………………………………... 65
V.7.2 La surface …………………………………………………………………....... 65
V.8 Optimisation des formes géométriques ……………………………………… 65
V.9 Exercices avec réponses ……………………………………………………………. 67

53
 Ecoulements autour des obstacles

V.1 Introduction :
En aérodynamique et hydrodynamique, cas des hydrofoils les forces exercées par le
fluide sur le foil peuvent être décomposées en deux catégories qui se distinguent par
la nature du phénomène physique qui les provoque :
 La première catégorie comprend les forces de pression qui sont dues au
mouvement même du foil par rapport au fluide. Ces forces, perpendiculaires au
profil en tout point, peuvent être modélisées par une résultante unique,
perpendiculaire à la direction principale de l'écoulement et dirigée vers le haut : la
portance.
 La seconde catégorie comprend les forces de viscosité qui sont dues à la nature
visqueuse du fluide.

Considérons un profil d’objet non symétrique :

Les forces de pression sont normales et les forces de frottement sont tangentielles
sur la face supérieure il y a une dépression(𝑃 < 𝑃0 ), sur la face inférieure on
remarque qu’il y a une surpression (𝑃 > 𝑃0 ) comme montre la figure v.1.

Fig (V.1) Les forces de pression et les forces de frottement exercées sur un profil.

La résultante des forces de pression présente :

 Une composante verticale : la portance (Lift) ⃗⃗⃗⃗

𝐿𝑃 .
 Une composante horizontale : la trainée (drag) ⃗⃗⃗⃗𝐷𝑃 .
La résultante des forces de frottement présente :

 Une composante verticale : la portance (Lift) ⃗⃗⃗⃗

𝐿𝐹 .
 Une composante horizontale : la trainée (drag) ⃗⃗⃗⃗𝐷𝐹 .

La figure V.2 montre la trainée et la portance exercées sur le profil.

54
 Ecoulements autour des obstacles

⃗ = ⃗⃗⃗⃗
𝐿 𝐿𝑃 + ⃗⃗⃗⃗
𝐿𝐹

Et une force de trainée :

⃗ = ⃗⃗⃗⃗
𝐷 𝐷𝑃 + ⃗⃗⃗⃗
𝐷𝐹

Fig (V.2) La force de trainée et la force de portance exercées sur un profil.

Les coefficients de portance et trainée sont écrits respectivement comme suit :

𝐿
𝐶𝐿 =
0.5𝜌 . 𝑆 . 𝑉 2

𝐷
𝐶𝐷 =
0.5𝜌 . 𝑆 . 𝑉 2

Le coefficient de trainée 𝐶𝑥 dépend des caractéristiques suivantes :

𝑙
𝐶𝑥 = 𝑓 ( 𝐿𝑎 𝑓𝑜𝑟𝑚𝑒, 𝑅𝑒 , 𝑀𝑎 , 𝐹𝑟, )
𝜀

La forme de la géométrie a une grande influence sur la trainée.

Ex : Pour une sphère, et une géométrie qui ressemble à un casque d’un cycliste.
Pour un 𝑅𝑒 = 105 :

Géométrie (a) : 𝐶𝑥 = 0.5

Géométrie (b) : 𝐶𝑥 = 0.04 La trainée est réduite à 92%

Fig (V.3) Influence de la forme géométrique sur le coefficient de trainée.

(a) 𝐶𝑥 = 0.5 (b) 𝐶𝑥 = 0.04

55
 Ecoulements autour des obstacles

V.2 Evolution expérimentale de 𝑪𝒙 (𝑹𝒆 ) et description de

l’écoulement autour d’une sphère lisse:
V.2.1 Evolution expérimentale de 𝑪𝒙 (𝑹𝒆 ):
La figure ci-dessous reproduit l’évolution du coefficient de trainée 𝐶𝑥 en fonction du
𝜌2𝑈𝑅
𝑅𝑒 = pour une sphère de rayon 𝑅.
𝜇

Zone 01 Zone 02 Zone 03 Zone 04

𝑹𝒆𝒄 ≈ 𝟓. 𝟏𝟎𝟓
𝟐𝟒
𝑪𝒙 =
𝑹𝒆

Fig (V.4) L’évolution du coefficient de trainée en fonction du nombre de Reynolds 𝑅𝑒

(cas d’une sphère lisse)

Les valeurs numériques particulières sont indiquées sur le tableau V.1.

Tab (V.1) Valeurs numériques du coefficient de trainée en fonction du nombre de Reynolds 𝑅𝑒

(cas d’une sphère lisse)
𝑹𝒆 0.1 1 10 102 103 104 105 5.105 106
𝑪𝒙 245 28 4.4 1.1 0.47 0.43 0.47 discontinuité 0.14

Cette courbe fait apparaitre 4 zones caractéristiques :

Tab (V.2) Description des zones caractéristiques de la courbe du coefficient de trainée
(cas d’une sphère lisse)
Zone 1 Zone 2 Zone 3 Zone 4
𝑅𝑒 < 1 1 < 𝑅𝑒 < 103 103 < 𝑅𝑒 < 𝑅𝑒𝑐 𝑅𝑒 > 𝑅𝑒𝑐 ≈ 5. 105
𝑎
24 𝐶𝑥 = 𝐴 (1 + ) 𝐶𝑥 presque 𝐶𝑥 chute lorsque
𝐶𝑥 = √𝑅𝑒
𝑅𝑒 constant l’écoulement devient
Avec : 𝑎 = 9.06
turbulent. Il évolue
Et 𝐴𝑎2 = 24 assez peu ensuite

56
 Ecoulements autour des obstacles

Calculant la force de trainée 𝐹𝑥 :

1
𝐹𝑥 = 𝐶 𝜌𝑈 2 𝑆
2 𝑥
Avec :
𝜇
𝐶𝑥 = 24.
2𝜌𝑅𝑈

On a :

𝑆 = 𝜋𝑅 2

Donc :

𝜇 𝜌𝑈 2 2
𝐹𝑥 = 24. 𝜋𝑅
2𝜌𝑅𝑈 2

Et on aura :

𝐹𝑥 = 6 𝜋𝜇𝑅𝑈

L’approximation de cette force de trainée pour une sphère à bas de nombre de

Reynolds 𝐹𝑥 = 6 𝜋𝜇𝑅𝑈 constitue la formule de Stokes.

V.2.2 Description de l’écoulement :

Pour une sphère de rayon 𝑅 le coefficient de trainée 𝐹𝑥 est défini par :
1
𝐹𝑥 = 𝐶𝑥 𝜌𝑈 2 (𝜋𝑅 2 )
2

Rappelons la forme de l’évolution de coefficient de trainée 𝐶𝑥 (𝑅𝑒 ) (Figure V.4).

(a) Pour les valeurs de 𝑹𝒆 < 𝟏 :

L’écoulement est laminaire et approximativement linéaire. Les lignes de courant ont
l’allure représentée sur la figure. 𝐶𝑥 est inversement proportionnel à 𝑅𝑒 .

(b) Pour des valeurs de 𝑹𝒆 > 𝟏 (de l’ordre de 𝑹𝒆 ≈ 𝟐𝟎) :

Il apparait un tourbillon stable derrière la sphère comme montre la figure ci-dessous.
Les dimensions de ce tourbillon augmentent avec 𝑅𝑒 .

57
 Ecoulements autour des obstacles

(c) Pour des nombre de Reynolds de l’ordre de 300 à 450 :

Ce tourbillon finit par occuper toute la partie arrière de la sphère, comme il est
montré sur la figure.

(d) A partir de de 𝑹𝒆 voisinant 450 :

Le tourbillon se détache, en prenant une forme hélicoïdale. Ce tourbillon a pour
conséquence l’existence d’une force transversale tournante s’exerçant sur la sphère.

(e) Pour de 𝑹𝒆 ≈ 𝟏𝟎𝟎𝟎 :

L’écoulement n’est plus régulier, il se forme un sillage, zone turbulente et chaotique
derrière la sphère.

(f) Pour une valeur critique 𝑹𝒆𝒄 ≈ 𝟓. 𝟏𝟎𝟓 :

La couche limite effectue une transition laminaire-turbulente :
 Le point de décollement de la couche limite change de position (il recule).
 Le coefficient de trainée 𝐶𝑥 chute lors de cette transition.

Si 𝑅𝑒 augmente encore 𝑅𝑒 > 5. 105 ; le sillage se rétrécit, les tourbillons sont

chaotiques. Le coefficient de trainée 𝐶𝑥 évolue assez peu.

(a) (b)

(c) (d)

58
 Ecoulements autour des obstacles

(e) (f)
Fig (V.5) Description de l’écoulement en fonction du nombre de Reynolds 𝑅𝑒
(cas d’une sphère lisse)

V.3 Evolution expérimentale de 𝑪𝒙 (𝑹𝒆 ) et description de

l’écoulement autour d’un cylindre:
V.3.1 Evolution expérimentale de 𝑪𝒙 (𝑹𝒆 ):
L’allure de la variation de 𝐶𝑥 pour un cylindre en fonction de 𝑅𝑒 est défini par :

𝜌𝐷𝑈
𝑅𝑒 =
𝜇

L’allure de la variation de 𝐶𝑥 pour un cylindre en fonction du 𝑅𝑒 , est donnée sur la

figure ci-dessous.

𝑪𝒙

𝑹𝒆

Fig (V.6) L’évolution du coefficient de trainée en fonction du nombre de Reynolds 𝑅𝑒

(cas d’un cylindre)

59
 Ecoulements autour des obstacles

Les valeurs numériques particulières sont indiquées sur le tableau V.3.

Tab (V.3) Valeurs numériques du coefficient de trainée en fonction du nombre de Reynolds 𝑅𝑒

(cas d’un cylindre)
𝑹𝒆 10-2 10-1 1 10 102 103 104 105 5.105 106 107
𝑪𝒙 400 60 11 3 1.2 1.1 1.12 1.23 discontinuité 0.4 0.6

V.3.2 Description de l’écoulement :

(a) Pour de faibles valeurs de 𝑹𝒆 ; de l’ordre de l’unité :

L’écoulement est laminaire et approximativement linéaire. Les lignes de courant ont
l’allure représentée sur la figure. Nous sommes en présence d’un écoulement
laminaire visqueux non décollé.

(b) Pour des valeurs de 𝑹𝒆 de quelques dizaines :

Il apparait des tourbillons stables derrière le cylindre comme il est illustré sur la
figure. L’écoulement est toujours symétrique ; il existe un décollement laminaire
symétrique.

(c) Pour un Reynolds de l’ordre de 100 :

Des tourbillons se détachent périodiquement de cylindre, alors que d’autres se
forment. L’écoulement est périodique mais les lignes de courant sont encore
identifiables. Des tourbillons dits Van Karman apparaissent.

(d) Pour un 𝑹𝒆 compris entre 𝑹𝒆 = 𝟑𝟎𝟎 et 𝑹𝒆 = 𝟓. 𝟏𝟎𝟓 :

Le décollement de couche limite s’effectue toujours d’une couche limite laminaire.
Les tourbillons se détachent du cylindre et deviennent turbulents comme montre la
figure ci-dessous.

(e) Pour un Reynolds 𝟓. 𝟏𝟎𝟓 < 𝑹𝒆 < 𝟓. 𝟏𝟎𝟔 :

Le décollement s’effectue à partir d’une couche limite turbulente.

(f) valeur 𝑹𝒆 > 𝟏𝟎𝟔 :

Le 𝐶𝑥 reste quasiment constant.

La figure V.7 montre la structure de l’écoulement pour différents nombre de

Reynolds 𝑅𝑒 .

60
 Ecoulements autour des obstacles

(a) (b)

(c) (d)

(e) (f)
Fig (V.7) Description de l’écoulement en fonction du nombre de Reynolds 𝑅𝑒
(cas d’un cylindre)

V.4 Ecoulement autour d’un hydrofoil :

Le fluide s'écoule de part et d'autre du foil en suivant les lignes de courant qui se
séparent au bord d'attaque et se rejoignent au bord de fuite. Dans le cas d'une
incidence trop importante, la ligne de courant qui se trouve sur l'extrados ne peut pas
"coller" au profil sur toute sa longueur : On dit qu'il y a décollement. En conséquence,
il se crée une zone de recirculation due à la viscosité du fluide. Dans cette zone,

61
 Ecoulements autour des obstacles

l'écoulement devient turbulent, il y a une formation de tourbillons. C'est ce qu'on

appelle le décrochage. Le décollement de la couche limite laminaire provoque une
baisse de portance

Coefficient de
portance 𝑪𝒛

Angle d’incidence 𝒊

Fig (V.8) Coefficient de portance 𝐶𝑧 en fonction de l’angle d’incidence.

(cas d’un hydrofoil)

V.5 Ecoulement autour d’une hélice:

Pour calculer les performances hydrodynamiques d’une hélice marine, les valeurs de
la portance et la trainée sont estimées avant. Lorsqu’un écoulement est généré autour
d’une section d’une hélice marine, ce dernier crée un angle d’incidence avec la
section. La force de portance et la trainée sont générées. Les performances
hydrodynamiques sont calculées à la base de la portance et la trainée :

𝑑𝑇 = 𝑍[𝑑𝐿(𝑟). 𝑐𝑜𝑠𝛽(𝑟) − 𝑑𝐷(𝑟). 𝑠𝑖𝑛𝛽(𝑟)]

𝑑𝑄 = 𝑍[𝑑𝐿(𝑟). 𝑠𝑖𝑛𝛽(𝑟) + 𝑑𝐷(𝑟). 𝑐𝑜𝑠𝛽(𝑟)]𝑟𝑑𝑟

62
 Ecoulements autour des obstacles

Comme pour le cas des hydrofoils, le rendement de l’hélice chute à une certaine
vitesse d’écoulement (décrochage) et ça est dû à l’augmentation de l’angle
d’incidence.

Fig (V.9) Cas de l’écoulement autour d’une hélice marine.

(a) Forces agissantes sur une section d’hélice (b) Performances hydrodynamiques

V.6 Ecoulement autour d’une plaque plane (calcul du 𝑪𝒙 ):

On a la formule de la trainée :

𝐷 = 1⁄2 . 𝜌 . 𝑆 . 𝑈 2 . 𝐶𝑥

𝐷 : La trainée.
𝐶𝑥 : Coefficient de trainée.
𝑈 : Vitesse de l’écoulement relative à l’objet.

Pour un écoulement laminaire :

−1⁄
𝐶𝑥 = 1.328 𝑅𝑒𝑙 2 pour 𝑅𝑒 < 5. 105

Pour un écoulement transitoire :

−1⁄ 1740
𝐶𝑥 = 0.0743 𝑅𝑒𝑙 5 − pour 5. 105 < 𝑅𝑒 < 107
𝑅𝑒𝑙

63
 Ecoulements autour des obstacles

0.455 1740
𝐶𝑥 = − pour 5. 105 < 𝑅𝑒 < 109
(𝐿𝑜𝑔𝑅𝑒𝑙 ) 2.58 𝑅𝑒𝑙

Pour un écoulement turbulent :

(a) Turbulent lisse :
−1⁄
𝐶𝑥 = 0.0743 𝑅𝑒𝑙 5 pour 5. 105 < 𝑅𝑒 < 107

0.455
𝐶𝑥 = pour 5. 105 < 𝑅𝑒 < 109
(𝐿𝑜𝑔𝑅𝑒𝑙 )2.58

(b) Complètement turbulent :

𝐶𝑥 = [1.89 + 1.62𝑙𝑜𝑔 (𝜀 /𝑙)]−2.5 dépendant de la rugosité

Fig (V.10) Coefficient de trainée 𝐶𝑥 en fonction de nombre de Reynolds 𝑅𝑒

(cas d’une plaque plane à différents coefficient de rugosité)

64
 Ecoulements autour des obstacles

V.7 Quelques paramètres influents sur le coefficient de trainée 𝑪𝒙 :

V.7.1 Position de la plaque :
Quand la plaque est allignée parallèlement à l’écoulement,les lignes de courant sont
presque uniformes. La force de trainée est due aux forces de viscosités. La force de
pression est negligeable. Quand la plaque est perpendiculaire, La force de trainée est
due aux forces de pressions. . La force de viscosité est negligeable.

Fig (V.11) Ecoulement autour d’une plaque plane.

(a) une plaque parallèle à l’écoulement (b) une plaque perpendiculaire à l’écoulement

V.7.2 La surface :
Le coefficient de trainée dépend aussi de la forme géométrique, et la surface consiste
un paramètre important pour le calcul 𝐶𝑥 . La figure V.12 montre la surface utilisée
pour differents formes géométriques.

Fig (V.12) Choix de la surface pour le calcul du Coefficient de trainée 𝐶𝑥 .

(différentes formes géométriques)

V.8 Optimisation des formes géométriques :

Pour déplacer des obstacles sur une surface libre la traînée due aux vagues augmente.
Les formes des coques de navire sont conçues afin de réduire la traînée due aux
vagues. L’utilisation d’un bulbe réduit la résistance due aux vagues.

65
 Ecoulements autour des obstacles

Fig (V.13) Coefficnet de trainée 𝐶𝑥 en fonction du nombre de Froude 𝐹𝑟

(une coque avec et sans un bulbe)

66
 Ecoulements autour des obstacles

V.9 Exercices avec réponses:

Exercice 01:
Calculer la force de traînée due au mât de 20 mètres de haut d'un voilier dans un vent de
30 𝑛𝑜𝑒𝑢𝑑𝑠, si 𝑅 = 10 𝑐𝑚, sachant que 1 𝑛𝑜𝑒𝑢𝑑 nautique vaut 1.852 km. La courbe du
𝐶𝑥 (𝑅𝑒 ) pour le cas d’un cylindre et une sphère est illustrée (voir annexe V.1).

On donne :
𝜈𝑎𝑖𝑟 = 1,5 .10−5 𝑚2 /𝑠 , 𝜌𝑎𝑖𝑟 = 1.23 𝑘𝑔/𝑚3

Réponse :
1. 𝐹𝑥 = 744 𝑁

Exercice 02:
Une plaque plane mince de 14 𝑚 de longueur, et de 4 𝑚 de largeur est remorquée
parallèlement à elle-même dans le sens de sa longueur dans l’eau à la vitesse
de 120 𝑘𝑚/ℎ.
Déterminer la force de frottement s’exerçant sur la plaque ?
On remplace la plaque par un obstacle cylindrique mince de diamètre 𝐷 = 0.25 𝑚.
Calculer la force de frottement sur l’une des surfaces de cet obstacle ?

On donne :
1.328
𝑅𝑒 < 105 : 𝐶𝑥 =
√ 𝑅𝑒
0.074
105 < 𝑅𝑒 < 107 : 𝐶𝑥 = 1⁄
𝑅𝑒 5
𝑅𝑒 > 107 : 𝐶𝑥 = 0.455. (𝐿𝑜𝑔10 𝑅𝑒 )−2.58 , 𝜈𝑒.𝑚𝑒𝑟 = 1,2 .10−6 𝑚2 /𝑠 , 𝜌𝑒.𝑚𝑒𝑟 = 1025 𝑘𝑔/𝑚3

Réponse :
1. 𝐹𝑥 𝑝𝑙𝑎𝑞𝑢𝑒 = 56,43. 103 𝑁
2. 𝐹𝑥 𝑐𝑦𝑙𝑖𝑛𝑑𝑟𝑒 = 88,4 𝑁

Exercice 03:
On veut propulser à la vitesse de 72 𝑘𝑚/ℎ une vedette de masse égale à 10 tonnes
en la faisant supporter seulement par une aile sous-marine (hydrofoil) dont le
coefficient de portance vaut 0,7 et le coefficient de trainée 0,06.

On donne :
𝜌𝑒.𝑚𝑒𝑟 = 1025 𝑘𝑔/𝑚3

Calculer :
1. La surface de l’aile.
2. La puissance nécessaire à la propulsion.

67
 Ecoulements autour des obstacles

Réponse :
1. 𝑆 = 0,684 𝑚2
2. 𝑃𝑢𝑖 = 168,264 𝑘𝑤

Exercice 04:
Un sous-marin a une longueur 𝐿 = 84 𝑚 et une surface totale de coque de 𝑆 = 1800 𝑚2 .
Calculer la résistance due aux forces de frottement visqueux s’exerçant sur la coque,
sachant que le sous-marin a une vitesse 𝑉 = 15 𝑚/𝑠. Pour faire ce calcul, on admettra
qu’à la valeur 𝐶𝑥 du coefficient moyen de frottement établi pour une plaque plane est
applicable ici.

On donne :
𝑅𝑒 > 107 : 𝐶𝑥 = 0.455. (𝐿𝑜𝑔10 𝑅𝑒𝑙 )−2.58 𝜈 = 1,2. 10−6 𝑚2 /𝑠 𝜌 = 1025 𝑘𝑔/𝑚3

Réponse :
1. 𝐹𝑥 𝑠𝑜𝑢𝑠 𝑚𝑎𝑟𝑖𝑛 = 324,026. 103 𝑁

68
Chapitre VI

TURBULENCE

VI TURBULENCE ………………………………………………………………....……….. 69
VI.1 Introduction .………………………………………………………….………….……. 70
VI.2 Les equations moyennées ………………………………………………….……... 70
VI.2.1 Principe ……………………………………………….……………………….... 70
VI.3 Rappel des equations instantanées ……………………………………………. 71
VI.3.1 Equations principales ………………………..………………………….…. 71
VI.4 Critère de choix et classification ………………………………………………… 74

69
 Turbulence

VI.1 Introduction :
Quand le champ de vitesse subit d’importantes fluctuations spatiales et temporelles
on dit que l’écoulement est turbulent. Il est très compliqué de décrire ce champ, des
méthodes statistiques sont employées en remplaçant les propriétés instantanées par
des valeurs moyennes.

VI.2 Les equations moyennées :

VI.2.1 Principe :

 La résolution des équations de Navier Stokes peut se faire pour des nombres de
Reynolds faibles et pour des cas simples.
 Les problèmes de la mécanique des fluides sont pratiquement rencontrés dans les
écoulements turbulents c'est-à-dire pour de nombres de Reynolds importants qui
consistent à s’intéresser aux quantités moyennes.
 Pour cela l’opérateur moyen d’ensemble est appliqué aux équations de
mouvement en pratiquant une décomposition de Reynolds sur les inconnus du
problème.
 Les nouvelles equations sont appelées les equations moyennées par opposition
aux equations du mouvement dite instantanées.

Turbulent
Laminaire

t0 t0 +T

Fig (VI.1) Champs de vitesse en un point dans un Fig (VI.2) Structure de l’écoulement
régime turbulent. laminaire et turbulent
(cas d’une fumée de cigarette)

70
 Turbulence

1 𝑡0 +𝑇
𝑢̅ = ∫ 𝑢 𝑑𝑡
𝑇 𝑡0

On a :

𝑢 = 𝑢̅ + 𝑢́

𝑢: Vitesse instantanée.
𝑢:
̅ Vitesse moyenne.
𝑢́ : Fluctuation de la vitesse.
On obtient la vitesse moyenne en réalisant 𝑁 mesures.
𝑁
1
𝑢̅ = ∑ 𝑢𝑖
𝑁
𝑖=1

VI.3 Rappel des equations instantanées :

VI.3.1 Equations principales :
 Incompressibilité (continuité) :

𝜕𝑢𝑗
=0
𝜕𝑥𝑗

 Navier Stokes :

𝜕𝑢𝑖 𝜕𝑢𝑖 1 𝜕𝑃 𝜕 2 𝑢𝑖
+ 𝑢𝑗 =− +𝜈
𝜕𝑡 𝜕𝑥𝑗 𝜌 𝜕𝑥𝑖 𝜕𝑥𝑗 𝜕𝑥𝑗

Equation de continuité :

𝜕𝑢𝑗
=0
𝜕𝑥𝑗

𝜕(𝑢̅ + 𝑢′ ) 𝜕(𝑣̅ + 𝑣 ′ ) 𝜕(𝑤

̅ + 𝑤 ′)
+ + =0
𝜕𝑥 𝜕𝑦 𝜕𝑧

𝜕𝑢̅ 𝜕𝑢′ 𝜕𝑣̅ 𝜕𝑣 ′ 𝜕𝑤̅ 𝜕𝑤 ′

+ + + + + =0
𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑦 𝜕𝑧 𝜕𝑧

71
 Turbulence

𝜕𝑢̅ 𝜕𝑣̅ 𝜕𝑤̅ 𝜕𝑢̅𝑗

+ + =0 =0 (1)
𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑥𝑗
⇒
𝜕𝑢′ 𝜕𝑣 ′ 𝜕𝑤 ′ 𝜕𝑢𝑗′
+ + =0 =0 (2)
{ 𝜕𝑥 𝜕𝑦 𝜕𝑧 {𝜕𝑥𝑗

(1) et (2) vérifient l’équation de continuité.

Equation de quantité de mouvement (Navier Stokes):

𝜕𝑢𝑖 𝜕𝑢𝑖 1 𝜕𝑃 𝜕 2 𝑢𝑖
+ 𝑢𝑗 =− +𝜈
𝜕𝑡 𝜕𝑥𝑗 𝜌 𝜕𝑥𝑖 𝜕𝑥𝑗 𝜕𝑥𝑗

On a :

𝑢𝑖 = 𝑢̅𝑖 + 𝑢𝑖′

𝑃 = 𝑃̅ + 𝑃′

On remplace dans l’équation de quantité de mouvement :

𝜕(𝑢̅𝑖 + 𝑢𝑖′ ) ′
𝜕(𝑢̅𝑖 + 𝑢𝑖′ ) 1 𝜕(𝑃̅ + 𝑃′ ) 𝜕 2 (𝑢̅𝑖 + 𝑢𝑖′ )
+ (𝑢̅𝑗 + 𝑢𝑗 ) =− +𝜈
𝜕𝑡 𝜕𝑥𝑗 𝜌 𝜕𝑥𝑖 𝜕𝑥𝑗 𝜕𝑥𝑗

𝜕(𝑢̅𝑖 + 𝑢𝑖′ ) 𝜕(𝑢̅𝑖 + 𝑢𝑖′ ) 𝜕(𝑢̅𝑖 + 𝑢𝑖′ ) 1 𝜕(𝑃̅ + 𝑃′ ) 𝜕 2 (𝑢̅𝑖 + 𝑢𝑖′ )

+ 𝑢̅𝑗 + 𝑢𝑗′ =− +𝜈
𝜕𝑡 𝜕𝑥𝑗 𝜕𝑥𝑗 𝜌 𝜕𝑥𝑖 𝜕𝑥𝑗 𝜕𝑥𝑗

𝜕(𝑢̿𝑖 + 𝑢̅𝑖′ ) 𝜕𝑢̿𝑖 𝜕𝑢̅𝑖′ ′

𝜕𝑢̿𝑖 ′
𝜕𝑢̅𝑖′ 1 𝜕(𝑃̿ + 𝑃̅ ′ ) 𝜕 2 (𝑢̿𝑖 + 𝑢̅𝑖′ )
+ 𝑢̿𝑗 + 𝑢̿𝑗 + 𝑢̅𝑗 + 𝑢̅𝑗 =− +𝜈
𝜕𝑡 𝜕𝑥𝑗 𝜕𝑥𝑗 𝜕𝑥𝑗 𝜕𝑥𝑗 𝜌 𝜕𝑥𝑖 𝜕𝑥𝑗 𝜕𝑥𝑗

On a aussi :

𝑢̅𝑖′ = 0 et 𝑢̿𝑖 = 𝑢̅𝑖

On obtiendra :

𝜕𝑢̅𝑖 𝜕𝑢̅𝑖 ̅̅̅̅̅̅̅̅

′
𝜕𝑢𝑖′ 1 𝜕𝑃̅ 𝜕 2 𝑢̅𝑖
+ 𝑢̅𝑗 + 𝑢𝑗 =− +𝜈
𝜕𝑡 𝜕𝑥𝑗 𝜕𝑥𝑗 𝜌 𝜕𝑥𝑖 𝜕𝑥𝑗 𝜕𝑥𝑗

̅̅̅̅̅̅̅ 𝜕𝑢𝑖′ ̅̅̅̅̅̅̅̅

𝜕𝑢𝑖′ 𝑢𝑗′ ̅̅̅̅̅̅̅̅ 𝜕𝑢𝑗′
= 𝑢𝑗′ + 𝑢𝑖′
𝜕𝑥𝑗 𝜕𝑥𝑗 𝜕𝑥𝑗

̅̅̅̅̅̅̅
𝜕𝑢𝑖′ 𝑢𝑗′ ̅̅̅̅̅̅̅̅
𝜕𝑢𝑖′
′
= 𝑢𝑗
𝜕𝑥𝑗 𝜕𝑥𝑗

72
 Turbulence

𝜕𝑢̅𝑖 𝜕𝑢̅𝑖 ̅̅̅̅̅̅̅

𝜕𝑢𝑖′ 𝑢𝑗′ 1 𝜕𝑃̅ 𝜕 2 𝑢̅𝑖
+ 𝑢̅𝑗 + =− +𝜈
𝜕𝑡 𝜕𝑥𝑗 𝜕𝑥𝑗 𝜌 𝜕𝑥𝑖 𝜕𝑥𝑗 𝜕𝑥𝑗

𝜕𝑢̅𝑖 𝜕𝑢̅𝑖 1 𝜕𝑃̅ 𝜕 2 𝑢̅𝑖 ̅̅̅̅̅̅̅

𝜕𝑢𝑖′ 𝑢𝑗′
+ 𝑢̅𝑗 =− +𝜈 −
𝜕𝑡 𝜕𝑥𝑗 𝜌 𝜕𝑥𝑖 𝜕𝑥𝑗 𝜕𝑥𝑗 𝜕𝑥𝑗

On connait que :
𝜕𝑢𝑖
𝜏𝑖𝑗 = 𝜇
𝜕𝑥𝑗
𝜇
𝜈=
𝜌

On aura :

𝜕 2 𝑢̅𝑖 1 𝜕𝜏̅̅̅
𝑖𝑗
𝜈 =
𝜕𝑥𝑗 𝜕𝑥𝑗 𝜌 𝜕𝑥𝑗

̅̅̅̅̅̅̅
𝜕𝑢𝑖′ 𝑢𝑗′ 𝜌 ̅̅̅̅̅̅̅̅̅
1 𝜕𝜌𝑢 ′ ′
𝑖 𝑢𝑗 1 𝜕𝑅𝑖𝑗
− . =− . =−
𝜕𝑥𝑗 𝜌 𝜌 𝜕𝑥𝑗 𝜌 𝜕𝑥𝑗

𝑅𝑖𝑗 = −𝜌̅̅̅̅̅̅
𝑢𝑖′ 𝑢𝑗′ (Tenseur de Reynolds)

L’équation de Navier Stokes moyennée s’écrit :

𝜕𝑢̅𝑖 𝜕𝑢̅𝑖 1 𝜕𝑃̅ 1 𝜕

+ 𝑢̅𝑗 =− + (𝜏̅̅̅ + 𝑅𝑖𝑗 )
𝜕𝑡 𝜕𝑥𝑗 𝜌 𝜕𝑥𝑖 𝜌 𝜕𝑥𝑗 𝑖𝑗

Avec :

𝑢̅𝑖′ = 0
𝑢̅𝑗′ = 0
̅̅̅̅̅̅
𝑢𝑖′ 𝑢𝑗′ ≠ 0

 Les composantes 𝑅𝑖𝑗 sont appelées contraintes de Reynolds ou tensions de

Reynolds, elles sont analogues avec les contraintes visqueuses.
 L’équation de Navier Stokes moyennée différe de l’équation de Navier Stokes
instantanée par la présence d’un terme supplémentaire représentant l’effet du
champ fluctuant sur le champ moyen.
 Le tenseur de Reynolds est symétrique :

̅̅̅̅̅̅
𝑢′ 𝑢′ ̅̅̅̅̅̅
𝑢′ 𝑣 ′ ̅̅̅̅̅̅
𝑢′ 𝑤 ′
𝑅𝑖𝑗 = −𝜌 [ ̅̅̅̅̅̅
𝑣 ′ 𝑢′ 𝑣̅̅̅̅̅̅
′𝑣′ ̅̅̅̅̅̅
𝑣′𝑤′ ]
̅̅̅̅̅̅
𝑤𝑢 𝑤′ ′ ̅̅̅̅̅̅
′ ′
𝑣 𝑤̅̅̅̅̅̅̅
′𝑤′

73
 Turbulence

 En situation de turbulence développée (loin de paroi), le tenseur de Reynolds est

prépondérant devant le tenseur de contraintes visqueuses. 𝑅𝑖𝑗 ≫ ̅̅̅
𝜏𝑖𝑗

Tab (VI.1) Inconnus à calculer pour le cas des équations instantanées et moyennées.
Inconnues Equations Inconnus Nombre total
à calculer disponibles supplémentaires d’inconnus
(4) (4) (6) (10)
Champ de vitesse moyen 𝑢
̅𝑗 (3) Continuité Tensions de Champ de vitesse
Pression moyenne 𝑃̅ (1) Navier Stokes Reynolds ̅̅̅̅̅̅
𝑢𝑖′ 𝑢𝑗′ (6) moyen (4)
Pression moyenne
Tensions de
Reynolds (6)

 Il est donc nécessaire de modéliser les inconnues supplémentaires à l’aide de loi

« modèles de fermeture »
 La décomposition de Reynolds dans les équations de Navier Stokes a permis de
trouver un certain nombre d’équations pour les inconnus du problème.
 Le nombre d’inconnues est toujours supérieur aux nombre d’équations.
 On peut utiliser l’équation de Poisson à la place de l’équation de continuité.
 On obtient l’équation de Poisson en appliquant l’opérateur divergence à l’équation
de Navier Stokes.

VI.4 Critère de choix et classification :

Il existe deux grandes catégories de modélisation :
 Modèle du premier ordre : modèles à viscosité turbulente basé sur
l’hypothèse de Boussinesq qui consiste à modéliser directement les tensions de
Reynolds à l’aide de la viscosité turbulente 𝜇𝑡 .
o 0 équation (Longueur de mélange)
o 1 équation (𝑘 énergie cinétique turbulente)
o 2 équations (𝑘 − 𝜔, 𝑘 − 𝜀 )

 Modèles du 2 éme ordre : Les tensions de Reynolds sont calculées

directement, la modélisation se porte sur des moments d’ordre supérieur, la mise
en œuvre est plus délicate mais les résultats sont de meilleure qualité.
o 7 équations (𝑅𝑖𝑗 − 𝜀)

Le choix du modèle dépend du problème à traiter.

Les modèles à viscosité turbulente sont plus stables et couramment utilisés.
L’utilisation de 0, 1, 2 équations dépend de la complexité du problème à traiter.

74
Annexes

Anx (I.1) Anx (I.2)

Anx (I.3)

75
Anx (I.4) Anx (I.5)

Anx (II.1) Anx (II.2)

Anx (II.3) Anx (III.1)

76
Anx (V.1)

77
Références

[1] Mécanique des fluides 2 ème édition, Sakir Amiroudiene, Jean-Luc Battaglia.
[2] Fluid Mechanics Fundamentals and applications, Yunus A. Çengel, Jhon M.
Cimbala.
[3] Fluid Mechanics, Merl Potter, David C. Wiggert.
[4] Mécanique expérimentale des fluides, Tome 1 statique et dynamique des fluides
non visqueux, Raymond Comolet.
[5] Mécanique expérimentale des fluides, Tome 2 dynamique des fluides réels,
Turbomachines, Raymond Comolet.
[6] Introduction to fluid mechanics, Y. Nakayama.
[7] Turbulent flows, Stephen B. Pope.

78