Ecole Sup. Polytech.

de Kaya 2020-2021 GC1 et GHR1

MECANIQUE DES FLUIDES

Notes de cours et Exercices

Licence 1 Génie Civil et Génie Hydraulique

École Supérieure Polytechnique Kaya (ESPK)

Ouagadougou - Burkina Faso

24 novembre 2021

Dr. OUEDRAOGO Soumaïla

Maitre de Conférences en Physique Appliquée
Université Joseph Ki-ZERBO/ Institut des Sciences -
Ouagadougou

©Dr OUEDRAOGO 2 [email protected]

 TABLE DES MATIÈRES

1 Introduction à la mécanique des fluides : propriétés générales des fluides 4

1.1 Définition . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
1.1.1 Qu’est ce qu’un fluide ? . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
1.1.2 Particule fluide . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.2 Caractéristiques physiques des fluides . . . . . . . . . . . . . . . . . . . . . . . . 5
1.2.1 Masse volumique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.2.2 Compressibilité . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.2.3 Viscosité . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
1.3 Forces subies par un fluide . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
1.3.1 Force de volume : force de pesanteur . . . . . . . . . . . . . . . . . . . . . 9
1.3.2 Forces de surface : force de pression et force de frottement . . . . . . . . 10

2 Statique des fluides 13

2.1 Équation fondamentale de la statique des fluides . . . . . . . . . . . . . . . . . . 13
2.1.1 Autre méthode plus générale . . . . . . . . . . . . . . . . . . . . . . . . . . 15
2.2 Application aux fluides incompressibles : l’hydrostatique . . . . . . . . . . . . . 15
2.2.1 Équation fondamentale de l’hydrostatique . . . . . . . . . . . . . . . . . . 15
2.2.2 Conséquences . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
2.3 Forces hydrostatiques . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
2.3.1 Résultante des forces hydrostatiques . . . . . . . . . . . . . . . . . . . . . 19
2.3.2 Pression sur le fond horizontal des vases . . . . . . . . . . . . . . . . . . . 22
2.3.3 Détermination du point d’application d’une force hydrostatique . . . . . 22
2.4 Poussée d’Archimède . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
2.5 Stabilité d’équilibre d’un corps totalement immergé . . . . . . . . . . . . . . . . 25
2.6 Corps flottants . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
2.6.1 Définitions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
2.6.2 Différentes positions d’équilibre (cas d’un bateau) . . . . . . . . . . . . . 26
2.7 Application aux fluides compressibles . . . . . . . . . . . . . . . . . . . . . . . . . 27

1
Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

3 Dynamique des fluides parfaits incompressibles 28

3.1 Équation de continuité . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
3.1.1 Débits massique et volumique . . . . . . . . . . . . . . . . . . . . . . . . . 28
3.1.2 Bilans de masse : équation de continuité . . . . . . . . . . . . . . . . . . . 29
3.1.3 Équation de continuité pour divers écoulements . . . . . . . . . . . . . . 30
3.2 Équation d’Euler . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
3.2.1 Principe fondamental de la dynamique appliqué à un fluide parfait . . . 31
3.2.2 Application à un écoulement unidirectionnel d’un fluide incompressible 33
3.3 Équation de Bernouilli . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
3.3.1 Hypothèses - Énoncé . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
3.3.2 Démonstration . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
3.3.3 Équation de Bernouilli d’un écoulement avec échange de travail . . . . . 37
3.4 Application du Théorème de Bernouilli . . . . . . . . . . . . . . . . . . . . . . . . 38
3.4.1 Écoulement par orifice : formule de Torricelli . . . . . . . . . . . . . . . . 38
3.4.2 Effet Venturi . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
3.4.3 Tube de Pitot . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
3.4.4 Phénomène de cavitation . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
3.5 Théorème d’Euler . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42
3.5.1 Hypothèses - Énoncé . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42
3.5.2 Application . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44

4 Dynamique des fluides réels incompressibles 45

4.1 Régime d’écoulement - Nombre de Reynolds . . . . . . . . . . . . . . . . . . . . 45
4.2 Pertes de charges . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
4.2.1 Définition . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
4.2.2 Pertes de charge singulières . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
4.2.3 Pertes de charges linéaires . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50

5 Annexe 1 : Analyse Vectorielle 53

5.1 Notion de champ . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
5.2 Gradient . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
5.2.1 Définition . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
5.2.2 Propriétés . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
5.2.3 Exemples . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
5.3 Divergence . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
5.3.1 Définition . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
5.3.2 Interprétations physiques . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
5.3.3 Propriétés de l’opérateur divergence . . . . . . . . . . . . . . . . . . . . . 56
5.3.4 Théorème de la divergence . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
5.4 Rotationnel . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
5.4.1 Définition . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57

©Dr OUEDRAOGO 2 [email protected]

Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

5.4.2 Interprétation physique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58

5.4.3 Propriétés . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58
5.5 Laplacien . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58
5.5.1 Définition . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58
5.6 Vecteur formel « nabla » . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59
5.6.1 Gradient, divergence, rotationnel . . . . . . . . . . . . . . . . . . . . . . . 59
5.6.2 Laplacien . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59

©Dr OUEDRAOGO 3 [email protected]

CHAPITRE 1
INTRODUCTION À LA MÉCANIQUE DES FLUIDES :
PROPRIÉTÉS GÉNÉRALES DES FLUIDES

La mécanique des fluides est la science des lois de l’é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 histori-
quement 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’hydrau-
lique, l’hydrodynamique, l’aérodynamique,etc. 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 très puissants : les supercalculateurs.

1.1 Définition
1.1.1 Qu’est ce qu’un fluide ?
C’est un milieu matériel :
– continu : ses propriétés varient d’une façon continue, propriétés considérées comme
caractéristiques non d’un point sans volume mais d’une particule, volume de fluide
extrêmement petit autour d’un point géométrique ; par exemple, on affecte à chaque
point P, pour chaque instant t, une masse volumique ρ représentative de la population
des molécules intérieures au volume dV de la particule ;

4
Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

– déformable (il n’a pas de forme propre) : les molécules peuvent facilement glisser les
unes sur les autres ; cette mobilité fait que le fluide prendra la forme du récipient qui
le contient ;
– qui peut s’écouler, mais tout fluide peut s’écouler plus ou moins facilement d’un ré-
cipient à un autre ou dans une conduite. Des forces de frottements qui s’opposent au
glissement des particules de fluide les unes contre les autres peuvent apparaître : on
parle de viscosité.
L’état fluide englobe deux des trois états de la matière : le liquide et le gaz. Les liquides et
gaz habituellement étudiés sont isotropes, c’est-à-dire que leurs propriétés sont identiques
dans toutes les directions de l’espace.

1.1.2 Particule fluide

La particule fluide est une portion de fluide à laquelle correspondent, à un instant t, une
vitesse, une pression, une température, une masse volumique, etc. Le volume envisagé est
très petit à notre échelle, mais doit contenir encore un très grand nombre de molécules pour
que les chocs moléculaires puissent être remplacés par la pression moyenne.
Les particules fluides ne sont pas des particules microscopiques sur lesquelles le mou-
vement brownien dû à l’agitation moléculaire est très perceptible ; la notion de continuité
repose sur celle de la compacité du réseau moléculaire intrinsèquement lacunaire.
Chaque particule d’un fluide est soumise à des forces de volume qui sont des forces à
longue distance induites par des champs de forces - le plus banal étant le champ de pesan-
teur - et à des forces de surface, forces de contact transmises à la surface de la particule
par les éléments environnants. On peut dire qu’un fluide est un corps homogène et continu
dont les diverses particules peuvent se déplacer ou se déformer sous l’action d’une force très
faible.

1.2 Caractéristiques physiques des fluides

1.2.1 Masse volumique
Définition
Considérons un milieu continu fluide à l’intérieur d’un volume V, et soit dV un volume
élémentaire défini autour d’un point M du volume V. Désignons par d m la masse de fluide
contenue dans le volume dV .
Le rapport ρ = d m/dV représente la masse volumique moyenne du fluide contenu dans
le volume dV .
On définit la masse volumique au point M par :

dm
ρ = lim
dV →0 dV

La masse m du fluide contenue dans le volume V est alors :

©Dr OUEDRAOGO 5 [email protected]

Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

Ñ
m= ρdV
V

La masse volumique s’exprime en kg .m −3 .

Ordres de grandeur des masses volumiques (à 20o C)

Corps ρ ( en kg .m −3
eau 1000
huile 914
mercure 13400
air 1,2

Les liquides sont caractérisés par une masse volumique relativement importante : ρ g az ¿
ρ l i qui d e
Pour les gaz, la masse volumique dépend de la température et de la pression.
Pour un gaz parfait, l’équation d’état donne
p
pV = nRT ⇒ ρ =
rT
R
où r = M est la constante massique des gaz parfaits (R = 8, 314 J .mol −1 .K −1 et M la masse
molaire du gaz.

Densité

ρ f l ui d e
d=
ρ f l ui d e/r é f ér ence
Pour les liquides, l’eau est le fluide de référence tandis que pour les gaz c’est l’air.

Poids volumique

P mg
γ=$= = = ρg
V V

1.2.2 Compressibilité
La propriété physique qui permet de faire la différence entre un liquide et un gaz est
la compressibilité. Un liquide est un fluide occupant un volume déterminé, ou du moins
ce volume ne peut varier que très peu, et seulement sous l’action de fortes variations de
pression ou de température. Un gaz, au contraire, occupe toujours le volume maximal qui
lui est offert : c’est un fluide essentiellement compressible (ou expansible).

©Dr OUEDRAOGO 6 [email protected]

Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

Définition
La compressibilité traduit la diminution de volume en réponse à un accroissement de
pression. Pour quantifier cet effet on introduit le coefficient de compressibilité isotherme
défini par :

1 ∂v
µ ¶
χ=−
v ∂p T
où v = ρ1 (m 3 kg −1 ) est le volume massique.
χ s’exprime en P a −1 .
Un accroissement de pression entraîne une diminution de volume, et inversement ; d’où
la nécessité de mettre un signe moins devant le coefficient de compressibilité.
Pour les gaz parfaits : χ = p1 .

Relation entre masse volumique et compressibilité

Le volume (et donc la masse volumique) peut varier sous l’effet de la pression ou de la
température. En plus du coefficient de compressibilité isotherme, on définit donc un coeffi-
cient de dilatation thermique à pression constante :

1 ∂v
µ ¶
α=
v ∂T p

Dans un fluide en mouvement les trois grandeurs p, v = 1/ρ et T ne sont pas uniformes et
l’équilibre thermodynamique n’est réalisé que localement, à l’échelle de la particule. L’équa-
tion différentielle d’état liée à ces grandeurs

∂v ∂v
µ ¶ µ ¶
dv = dp + dT
∂p T ∂T p

peut être transformée en faisant apparaître les deux coefficients χ et α :

d v = −χvd p + αvd T

Nous n’étudierons que des écoulements de liquides ou de gaz dans lesquels la tempé-
rature peut être considérée comme constante (dT = 0). L’approximation suivante sera donc
faite :
Liquide = fluide incompressible (χ = 0) ⇒ ρ = cte : fluide isovolume (dv = 0)
En pratique, il est raisonnable de dire qu’un fluide est isovolume si ∆ρ/ρ ≤ 4% au cours
de son mouvement. La variation pour l’eau est ∆ρ/ρ = 5.10−4 pour une variation de tem-
pérature ∆T = 1 K et ∆ρ/ρ = 5.10−4 pour une variation de pression ∆p = 1 bar . On peut
donc souvent traiter l’eau comme un fluide incompressible et utiliser dans les équations du
mouvement une masse volumique ρ = c t e.

©Dr OUEDRAOGO 7 [email protected]

Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

1.2.3 Viscosité
C’est une grandeur qui caractérise les frottements internes du fluide, autrement dit sa
capacité à s’écouler. Elle caractérise la résistance d’un fluide à son écoulement lorsqu’il est
soumis à l’application d’une force. C’est à dire, les fluides de grande viscosité résistent à
l’écoulement et les fluides de faible viscosité s’écoulent facilement. Elle peut être mesurée
par un viscosimètre à chute de bille, dans lequel en mesure le temps écoulé pour la chute
d’une bille dans le fluide. Elle peut également être mesurée par un récipient dont le fond
comporte un orifice de taille standardisée. La vitesse à laquelle le fluide s’écoule par cet ori-
fice permet de déterminer la viscosité du fluide.
On distingue la viscosité dynamique et la viscosité cinématique.

Viscosité dynamique
La viscosité dynamique exprime la proportionnalité entre la force qu’il faut exercer sur
une plaque lorsqu’elle est plongée dans un courant et la variation de vitesse des veines de
fluide entre les 2 faces de la plaque.
Elle est exprimée par un coefficient représentant la contrainte de cisaillement nécessaire
pour produire un gradient de vitesse d’écoulement d’une unité dans la matière.
Considérons deux couches de fluide adjacentes distantes de ∆z. La force de frottement
F qui s’exerce à la surface de séparation de ces deux couches s’oppose au glissement d’une
couche sur l’autre. Elle est proportionnelle à la différence de vitesse des couches soit ∆v, à
leur surface S et inversement proportionnelle à ∆z :

∆v
F = µS
∆z
Le facteur de proportionnalité µ est le coefficient de viscosité dynamique du fluide.

Exercice : Trouver l’unité du coefficient de viscosité dynamique

Ordre de grandeurs de la viscosité dynamique

Fluide µ ( en P a.s)
Eau (0o C) 1, 787 × 10−3
Eau (20o C) 1, 002 × 10−3
Eau (100o C) 0, 2818 × 10−3
Huile d’olive (20o C) 100 × 10−3
Glycérol (20o C) ≈ 1, 00
Oxygène 1, 95 × 10−5

©Dr OUEDRAOGO 8 [email protected]

Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

Viscosité cinématique

µ
ν=
ρ
L’unité légale de la viscosité cinématique est le m 2 .s −1 . Cependant, on utilise le Stokes
(St) comme unité pratique de mesure de la viscosité cinématique.

1St = 10−4 m 2 .s −1

Remarque :
– La viscosité des liquides diminue beaucoup lorsque la température augmente. Par contre,
celle des gaz augmente avec la température.
– Pour tous les fluides, la viscosité augmente avec la pression.
– 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

Conclusion
Les fluides peuvent être classés en :
– fluides parfaits (sans frottement). On dit aussi fluide non-visqueux
– fluides réels : avec frottement ou fluides visqueux,
– fluides incompressibles : liquides où ρ = c t e et
– fluides compressibles (gaz).
L’écoulement d’un fluide est caractérisé par :
– un champ de vitesse ;
– la pression ;
– certaines propriétés de ce fluide telles que la masse volumique et la viscosité.

1.3 Forces subies par un fluide

L’un des buts de la mécanique étant de définir la position ou le mouvement des particules
matérielles sous l’action des forces qui les sollicitent, il faut donc définir le genre de forces
que nous aurons à considérer en mécanique des fluides.

1.3.1 Force de volume : force de pesanteur

Les champs de force (de pesanteur, magnétique, électrique, etc.) exercent sur les parti-
cules fluides des actions à distance qui sont proportionnelles aux volumes des particules. Ce
sont les forces de volume.

©Dr OUEDRAOGO 9 [email protected]

Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

En supposant les particules fluides électriquement neutres, on ne considère que la den-

sité des forces exercées par la gravité sur ce milieu fluide continu.
Considérons un petit volume élémentaire dV et soit d F la force élémentaire qui s’exerce
sur dV (les forces de pesanteur) :

~ = d m~
dF g = ρ~
g dV

~:
Pour tout le volume de fluide, la résultante des forces de pesanteur définie le poids P
Z Z
~ = dF
P ~ = ρ~ g dV

1.3.2 Forces de surface : force de pression et force de frottement

Imaginons une surface S fictive qui, au sein du fluide, sépare le fluide en deux domaines
D 1 et D 2 . Les particules qui se trouvent du côté de D 2 , mais contiguës à S, agissent sur les
particules de D 1 qui le touchent. Ce sont des actions à courte distance proportionnelles à
l’aire de contact et on les appelle forces de surface.
Dans ce milieu, la force que le domaine 1 exerce sur le domaine 2 à travers un élément
−→
de surface réel ou fictif d S a une direction quelconque. Cette force d F peut être décomposée
en :
−−→
– une composante tangentielle d F t ;
−−→
– une composante normale d F n ;

– La quantité ddFSt représente la contrainte tangentielle. Ces contraintes n’apparaissent

qu’en dynamique des fluides et correspondent aux frottements visqueux des couches
fludies en mouvement par rapport aux autres et par rapport à la paroi de la conduite.
– La quantité ddFSn qui représente la contrainte normale est appelée la pression. C’est la
seule composante qui apparait en statique des fluides.

©Dr OUEDRAOGO 10 [email protected]

Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

La pression
La pression p est une grandeur scalaire (positive) définie en tout point du fluide. Dans le
SI, la pression s’exprime en pascal (P a = N .m −2 ).
Autres unités :
– le bar : 1 bar = 105 P a
– l’atmosphère : 1 at m = 1, 013.105 P a
– 1 at m = 760 mmH g = 760 t or r = 10, 33 mH2O
Mesure de la pression
– le baromètre pour la pression de l’air (pression atmosphérique) ;
– le manomètre avec pour référence la pression atmosphérique ;
– les capteurs numériques ou mécaniques

Représentation de la force de pression

La force de pression est une force normale. Pour définir le champ de pression qui s’ap-
plique sur une surface élémentaire d S, on choisit un vecteur unitaire ~ n orienté vers le milieu
qui agit.
Observons ces expériences élémentaires qui consistent à percer un petit trou dans un
récipient rempli de liquide. On constate que, quelle que soit la forme du récipient et la po-
sition de l’orifice, le liquide jaillit toujours perpendiculairement à la paroi. On admet que ce
qui a lieu sur la frontière du récipient se produit encore à l’intérieur. En d’autres termes, si S
est une surface non matérielle, qui sépare un domaine D de fluide en deux sous-domaines
D 1 et D 2 , alors le fluide dans D 2 exerce sur D 1 une force normale à S en tout point M de S,
et vice versa, le fluide dans D 1 exerce sur D 2 une force égale et opposée, donc normale elle
aussi à S (principe de l’action et de la réaction).
Pour exprimer la force exercée par D 2 sur D 1 on introduit le vecteur unitaire → −
n orienté
vers le milieu qui agit (ici D 2 ) et on écrit, pour un élément d S tracé sur S et entourant un
point M de S :
~ = −p~
dF nd S

dF~ la force exercée sur l’élément de sur-

face dS et p, la pression qui règne en un
point M
Par conséquent D 2 exerce sur D 1 , par l’intermédiaire de S, une densité surfacique de
force (force par unité de surface) −p →−
n . La force de pression agit toujours vers l’intérieur du
volume délimité par l’élément de surface.
La pression est indépendante de la surface et de l’orientation de cette surface. Si on fait
0
passer une autre surface S par le point M , on obtient la même pression p. La notion de pres-
sion, totalement indépendante de la frontière du fluide (récipient), décrit les efforts exercés
à l’intérieur du fluide par une partie D 2 sur son complémentaire D 1 . On dit que l’on a donné
une description des efforts intérieurs (que les particules exercent les unes sur les autres).

©Dr OUEDRAOGO 11 [email protected]

Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

Pression atmosphérique, pression absolue, pression relative, dépression et vide

– La pression atmosphérique p at m est due à la présence de l’atmosphère, elle dépend de
l’altitude, de la température et de la latitude. Elle est mesurée à l’aide d’un baromètre.
– La pression relative p r est prise par rapport à la pression atmosphérique. Elle est me-
surée à l’aide d’un manomètre.
– La pression absolue p a est définie par rapport au vide parfait. Elle est la somme des
pressions atmosphérique et relative (ou manométrique) : p a = p r + p at m
– Lorsque la pression est inférieure à la pression atmosphérique, on parle de dépression
ou de vide.

©Dr OUEDRAOGO 12 [email protected]

CHAPITRE 2
STATIQUE DES FLUIDES

La statique des fluides est la science qui étudie les conditions d’équilibre des fluides au
repos. Plus précisément, elle concerne toutes les situations dans lesquelles il n’y a pas de
mouvement relatif entre les particules fluides :
– fluides au repos
– fluides uniformément accélérés
Il n’y a pas de contraintes dues aux frottements entre particules. Les forces en jeu sont
uniquement des forces de volume dues au poids et de forces de surface dues à la pression.

2.1 Équation fondamentale de la statique des fluides

Considérons un élément de volume de forme parallélépipédique à l’intérieur d’un fluide
en équilibre, de volume dV = d xd yd z, dans un repère cartésien, et faisons le bilan des forces
qui s’appliquent sur cet élément de volume :

z La force de volume : le poids du fluide donné par

~ = d m~
dP g = ρdV ~
g

13
Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

z Les forces de surface dues à la pression ; on peut décomposer la résultante en trois

composantes :

−→
d F = d Fx −
e→x + d F y −
e→y + d F z →
−
ez

Puisque les forces de surface sont nécessairement normales, la composante suivant z

correspond aux forces de pression s’exerçant sur les surfaces perpendiculaires à l’axe z. Donc :
£ ¤
d F z = p (z) − p (z + d z) d xd y

Par un développement au premier ordre, on a :

∂p
p (z + d z) = p (z) + dz
∂z
D’où :
∂p ∂p
d Fz = − d xd yd z = − dV
∂z ∂z
Par analogie sur les deux autres axes :

∂p
d Fx = − dV
∂x

∂p
d Fy = − dV
∂y
La force de surface se résume alors à :

∂p − ∂p − ∂p →
µ ¶
−→ → → −
³−−−−→ ´
dF = − ex + ey + e z dV = − g r ad p dV
∂x ∂y ∂z
z La condition d’équilibre de la particule fluide permet d’écrire :
−→ −→ → −
dF + dP = 0
h −−−−→ →
−
→
−
i
−g r ad p + ρ g dV = 0

D’où la relation vectorielle fondamentale de la statique des fluides :

−−−−→ →
−
g r ad p = ρ g
→−
En supposant que : g = −g →
−
e z on a alors pour un fluide au repos :

∂p
 
 ∂x   
  0
 ∂p   
 = 0 
  
 ∂y 
 
 ∂p  −ρg
 

∂z
ou

©Dr OUEDRAOGO 14 [email protected]

Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

∂p ∂p
= =0
∂x ∂y
∂p
= −ρg
∂z
Cette dernière relation signifie que la pression ne dépend que de l’altitude z : p(x, y, z) =
p(z).
L’équation différentielle à résoudre pour connaitre la pression en tout point d’un fluide
au repos est alors :

dp
= −ρg
dz

2.1.1 Autre méthode plus générale

L’équation fondamentale de la statique peut s’établir d’une façon plus générale, sans faire
intervenir un repère particulier. Considérons une portion quelconque de fluide de volume V
et écrivons que la somme des forces qui s’exerce sur cette portion du fluide est nulle.
– Forces de volume :
−→
Ñ
→− →
− →
−
d P = ρdV g ⇒ P = ρdV g
V
– Forces de surface :

−−−−→
Ï Ñ
~ = −p~ →
−
dF nd S ⇒ F = − p~
nd S = − g r ad pdV
S V

en utilisant la formule du gradient (Voir fiche Calcul vectoriel)

−−−−→
Ñ Ñ
→
− →− →− →
− →
−
P +F = 0 ⇔ ρdV g − g r ad pdV = 0
V V
−
− − −→
Ñ ³
→
−
´ →
−
⇔ ρ g − g r ad p dV = 0
V
soit quelque soit le volume V :

− −−−−→
→ →
−
ρ g − g r ad p = 0

2.2 Application aux fluides incompressibles : l’hydrostatique

2.2.1 Équation fondamentale de l’hydrostatique
La masse volumique du fluide est en tout point la même : ρ = c t e (fluide incompressible).
Par ailleurs, on peut considérer que l’accélération de la pesanteur est une constante :
g = ct e.
Par conséquent :

©Dr OUEDRAOGO 15 [email protected]

Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

dp
= −ρg = c t e
dz
Et par intégration :

dp
Z Z
p (z) = dz =− ρg d z = − ρg z + c t e
dz
Soit :
p (z) + ρg z = c t e = p 0 + ρg z 0

où p 0 est la pression à l’altitude z 0

Donc :

p (z) = p 0 + ρg (z 0 − z) = p 0 + ρg h

où h est la hauteur de fluide sous le niveau de référence.

Un tel champ de pression, affine en z, est appelé champ de pression hydrostatique, et
l’équation ci-dessus encadrée est l’équation de l’hydrostatique (équation fondamentale de
la statique des fluides dans le cas particulier d’un fluide isovolume dans le champ de pesan-
teur).
La plupart du temps, on prendra z 0 = 0 le niveau de référence correspondant à la sur-
face libre du fluide où p 0 = p a . Pour les applications numériques, on prendra la pression
atmosphérique standard : p a = 1, 013 × 105 P a.

©Dr OUEDRAOGO 16 [email protected]

Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

Autre forme plus générale

p 1 + ρg z 1 = p 2 + ρg z 2 = c t e

p1 p2
+ z1 = + z 2 = c st e
ρg ρg
Ainsi, à tout point quelconque de cote z, on a :
p
+ z = H = cte
ρg

En hydraulique, cette constante (homogène à une longueur) est appelée la hauteur pié-
zométrique.

Retenons : Entre deux points A et B de cotes et pressions respectives (z A , p A ), (z B , p B ) ;

l’équation fondamentale de l’hydrostatique permet d’écrire pour un fluide de masse volu-
mique ρ :

p A − p B = ρg (z B − z A ) ou p A − p B + ρg (z A − z B ) = 0
¡ ¢

2.2.2 Conséquences
Les conséquences de l’équation fondamentale de l’hydrostatique sont les suivantes :
z Les surfaces d’égale pression dans un fluide homogène sont des plans horizontaux
(plans isobares). En effet, quand p = c t e nous avons z = c t e. Réciproquement la pres-
sion est constante dans un plan horizontal quelconque.

z Si nous avons deux fluides différents, de masses volumiques différentes, non miscibles,
la surface de séparation est horizontale. Le fluide le plus lourd est en dessous (équilibre
stable). En particulier, la surface libre d’un liquide surmonté d’un gaz au repos (comme
l’air atmosphérique) est un plan horizontal (plan où la pression est la pression atmo-
sphérique constante).
En effet :

0
Soit S l’interface entre deux fluides de masses volumiques respectives ρ et ρ . Soient A
0 0
(respectivement A ) et B (respectivement B ) deux points distincts sur S. On a :
0 0
p(A) = p(A ) et p(B ) = p(B ) en vertu du principe de l ?action et de la réaction,
0 0
z(A) = z(A ) et z(B ) = z(B ) puisque ce sont les mêmes points géométriques,

©Dr OUEDRAOGO 17 [email protected]

Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

0 0 0 0 0
p(A) − p(B ) = ρg [z(B ) − z(A)] et p(A ) − p(B ) = ρ g [z(B ) − z(A )] d’après l’équation
fondamentale de la statique.

0 0 0 0
D’où, ρg [z(B ) − z(A)] = ρ g [z(B ) − z(A )] soit (ρ − ρ )g [z(B ) − z(A)] = 0.
0
Si ρ 6= ρ , alors z(B ) = z(A), et ceci quels que soient les points A et B . Il en résulte que,
sous l’action de la gravité, la surface libre d’un liquide ou la surface de séparation de
deux liquides non miscibles en équilibre est un plan horizontal.

On en déduit le principe des vases communicants : dans plusieurs vases de forme quel-
conque, communiquant entre eux et contenant un seul liquide en équilibre, les surfaces
libres dans les différents vases sont dans le même plan horizontal.

Ce point est à la base de la mesure de la différence de pression entre deux gaz à l’aide
du manomètre en U. Il suffit de mesurer la différence de niveau du liquide dans les
deux branches et de connaître la masse volumique de ce liquide.
z La différence de pression p A − p B entre deux points quelconques A et B pris à l’inté-
rieur du fluide ne dépend que de la distance verticale entre les deux points. Elle est
égale au poids d’une colonne de fluide ayant comme base l’unité de surface et comme
hauteur la différence de niveau entre les deux points. En effet :

p A + ρg z A = p B + ρg z B
p A − p B = ρg (z B − z A )
z Dans un fluide incompressible en équilibre, les variations de pression se transmettent
intégralement en tout point du fluide. La différence (p B − p A ) calculée précédemment
reste constante quelles que soient les pressions. Si p A varie, p B varie simultanément de
la même quantité. Ceci constitue le théorème de Pascal : "Dans un fluide incompres-
sible en équilibre, toute variation de pression en un point entraine la même variation de
pression en tout autre point"

©Dr OUEDRAOGO 18 [email protected]

Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

2.3 Forces hydrostatiques

On intéresse à la détermination des forces s’exerçant sur les surfaces solides immergées.
Puisque l’équation fondamentale de la statique des fluides permet de déterminer la pression
en tout point d’un fluide, il est possible de calculer les forces s’exerçant sur les parois d’un
solide totalement ou partiellement immergé.
−→
Rappel : On sait que d F = −p → −
n d S est la force de pression élémentaire s’exerçant sur la
surface élémentaire d S. Pour connaître la force totale s’exerçant sur une surface S, il suffit
−→
d’intégrer d F sur cette surface S.
Ï
→
−
F = −p →
−
n dS
S

2.3.1 Résultante des forces hydrostatiques

Cas d’une plaque plane totalement immergée
Si l’on considère le cas particulier simple d’une plaque plane totalement immergée (voir
figure ci-dessous), la force hydrostatique s’exerçant sur la face supérieure est portée par la
normale → −n à cette surface :

©Dr OUEDRAOGO 19 [email protected]

Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

Si l’on considère le cas particulier simple

d’une plaque plane totalement immergée
(voir figure ci-dessous), la force hydrostatique
s’exerçant sur la face supérieure est portée par
la normale → −
n à cette surface :
dh
Z Z
→
−
F = −→−n p (h)d S = −→ −n p (h) ∆l
cos θ
Z h2
→
− →
− 1
p 0 + ρg h d h
¡ ¢
F =−n
cos θ h1

où ∆l est la largeur de la plaque (suivant la

perpendiculaire au plan du schéma, supposée
égale à l’unité : ∆l = 1), et θ est l’angle d’incli-
naison de la plaque par rapport à la verticale.
Si h est la profondeur à partir d’un niveau de
référence (surface libre) où la pression est la
pression atmosphérique p 0 .
h 2 −h 1
Ce qui donne après calcul et sachant que la surface totale S de la plaque vaut : S = cos θ
et hG = h1 +h
2
2
profondeur du barycentre (centre de gravité) de la surface S :

→
−
F = −→
−
n p 0 + ρg hG S
¡ ¢

Ce résultat peut être généralisé à n’importe quelle surface immergée, quelle que soit sa
→
−
forme : pour évaluer F il suffit donc de connaitre sa surface S et la profondeur de son bary-
centre hG .

Cas d’une plaque plane séparant un liquide de l’atmosphère

Considérons maintenant une paroi plane de surface S située entre les profondeurs h 1 et
h 2 au-dessous de la surface libre du liquide. Soit p 0 la pression atmosphérique. Cette surface
−
→ −
→ −
→
est soumise aux deux forces F 1 et F 2 normales à S et de directions opposées : F 1 la force
−
→
exercée par le fluide sur la paroi et la force F 2 exercée par le milieu extérieur sur la paroi.

©Dr OUEDRAOGO 20 [email protected]

Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

Les poussées élémentaires −p → −

n d S reçues par la paroi sont des vecteurs parallèles qui se
composent en une résultante normale à la surface ; son intensité est, puisque p = p 0 + ρg h :
Z µ Z ¶
−
→ −
→ −
→
p 0 + ρg h d S = −n 1 p 0 S + ρg hd S
¡ ¢
F 1 = −n 1
S S
R
Or : S hd S = hG S avec hG profondeur du barycentre (centre de gravité) de la surface S :
hG = h1 +h
2
2

on en déduit que :

−
→
F 1 = −−
→S ¡p + ρg h ¢
n 1 0 G

Par ailleurs :
−
→
F 2 = −−
→Sp = −
n 2 0
→Sp
n 1 0

D’où la résultante de ces deux actions :

→
− − → − →
F = F 1 + F 2 = −ρg hG S −
→ = ρg h S −
n 1
→
G n2

La pression atmosphérique s’applique de part et d’autre de la paroi et n’intervient donc

plus sur la résultante.
Ainsi, dans le cas d’une paroi séparant un liquide de l’atmosphère, la résultante des
forces de pression s’affranchit du terme dû à la pression atmosphérique.

©Dr OUEDRAOGO 21 [email protected]

Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

Cas d’une plaque à surface courbe

−→
d F = d Fx →−
x + d Fy→
−y
avec les projections suivantes sur les axes en
fonction de α
d F x = cosαpd S = cosαρg hd S
d F y = si nαpd S = si nαρg hd S
L’intégration de d F x et de d F y sur toute la
surface de l’élément courbe permet d’évaluer
la résultante F . On exprimera d’abord d S en
fonction de d α selon la configuration.

2.3.2 Pression sur le fond horizontal des vases

Quelle que soit la forme des vases, s’ils sont

remplis du même liquide à la même hauteur
h, le fond de même surface S est soumis à la
même force de pression, égale au poids d’une
colonne verticale de fluide de base S, de hau-
teur h :
F = ρg hS

Ce résultat a été appelé le paradoxe hydrostatique : la poussée verticale sur le fond d’un
récipient est indépendante de la forme de la paroi à l’intérieur de son contour qui est sup-
posé fixe.
La force de pression sur la surface inférieure est donc la même, bien que les vases contiennent
des poids de liquide différents, mais il faut remarquer que si dans le cas du vase cylindrique
la paroi inférieure supporte effectivement tout le poids de l’eau qui se trouve dans le vase, il
n’en est pas de même dans les deux autres cas ; par exemple, dans le vase tronconique évasé
vers le haut les parois latérales supportent des forces de pression qui ont une résultante ver-
ticale non nulle et dirigée vers le bas ; par contre, pour le vase resserré vers le haut, les parois
latérales inférieures supportent des forces de pression qui ont une résultante non nulle, mais
dirigée cette fois-ci vers le haut.

2.3.3 Détermination du point d’application d’une force hydrostatique

→
−
Si F est la force hydrostatique s’exerçant sur une surface S, alors on peut avoir besoin de
connaître le point d’application A de cette force. Ce point s’appelle centre de poussée.

©Dr OUEDRAOGO 22 [email protected]

Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

→
−
Le centre de poussée se détermine en écrivant que le moment de F par rapport à un
point O quelconque est égal à la résultante des moments élémentaires par rapport à ce
même point O.

−−→ −→
Z
−−→ → −
O A ∧ F = OM ∧ d F
S

Exemple : cas d’une paroi verticale

La paroi verticale possède un axe de

→
−
symétrie (G, Z ). G est son centre de
surface. D’un coté de la paroi il y a
un fluide de masse volumique ρ , de
l’autre coté, il y a de l’air à la pres-
sion atmosphérique p 0 . On désigne
par p G la pression au centre de sur-
face G du coté fluide.

z la force de pression en M est :

−→ →
− ¡ ¢¢ →− ¡ ¢ → −
d F = p M d S X = p G + ρg yG − y M d S X = p G − ρg y d S X
¡

→
−
La résultante F vaut :

−→
Z
→
− →
−
F = d F = pG S X
S
R →
−
puisque S yd S = yG S = 0 (Moment statique de la surface S par rapport à l’axe (G, Z ).

→
−
z Le moment par rapport à l’axe (G, Z ) des forces de pression :
−−→ −→
Z
−
→
MG = G M ∧ dF
Z
→
− ¡ ¢ → −
= y Y ∧ p G − ρg y d S X
·S Z Z ¸³
2 →
−´
= p G yd S − ρg y d S − Z
S S

Sachant que :
R
yd S = yG S = 0
RS 2 ³ − ´ (Moment quadratique de la surface S par rapport à l’axe (G, →
−
S y dS = I → G, Z
Z ) passant
par le centre de surface G. Donc :

©Dr OUEDRAOGO 23 [email protected]

Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

−
→ →
−
M G = ρg I ³ −´Z
→
G, Z

z Le centre de poussée G 0 est telle que :

−−→ −→ −−→ →
Z
−−→ − −−→ ~
~ G +G
MG = G M ∧ d F = GG 0 ∧ F ⇒ M 0G ∧ F = ~
0
S
Avec les résultats précédents, on obtient :
°−
ρg I ³ °→ °
°
→
−´
G, Z °M G °
y0 = − = − °→
°− °
°
pG S °F °

Remarque : Le centre de poussée est toujours au-dessous du centre de surface G.

2.4 Poussée d’Archimède

On cherche l’effort exercé sur un objet immergé, c’est-à-dire la force totale exercée par le
fluide sur l’obstacle qui occupe le volume V totalement entouré par le fluide.
→
− Î
On sait que cette force s’exprime par F = S −p → −
n dS
Î → − Ð −−−−→
La formule du gradient ( f n dS =
S g r ad f dV ) permet de passer d’une intégrale de
V
surface à une intégrale de volume soit :

−−−−→
Ñ
→
−
F = g r ad pdV
V

−−−−→ →
−
Or, l’équation fondamentale de la statique des fluides permet d’écrire g r ad p = ρ f l ui d e g
→
− →
−
D’où F = − V ρ f l ui d e g dV , et, en supposant g constant sur tout le volume,
Ð

Ñ
→
− →
− →−
F =−g ρ f l ui d e dV = −ρ f l ui d e V g
V

ρ f l ui d e V est la masse de fluide déplacé par le volume solide.

→
−
La poussée F n’a pas de composante horizon-
tale et sa composante verticale est égale et op-
posée au poids du fluide déplacé par le corps
(c’est la poussée d’Archimède).

Le point d’application de la poussée d’Archi-

mède, appelé centre de poussée, est le centre
de masse de la partie immergée du solide
(centre de gravité du fluide déplacé et non ce-
lui du solide immergé).

©Dr OUEDRAOGO 24 [email protected]

Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

2.5 Stabilité d’équilibre d’un corps totalement immergé

→
−
Pour que le solide immergé soit en équilibre il faut que son poids P soit égal et opposé à
→
− → − →
− →
− →
−
la poussée d’Archimède F : F = − P ou ρ sol i d e Vsol i d e g = −ρ f l ui d e Vsubmer g é g
→
− →
−
z Si le solide est homogène, le poids P et la poussée F sont appliqués au même point :
le centre de gravité G du solide coïncide avec le centre de poussée A.
– Si ρ sol i d e = ρ f l ui d e le solide est en équilibre indifférent dans le liquide.
– Si ρ sol i d e > ρ f l ui d e le solide descend jusqu’à toucher le fond du récipient.
– Si ρ sol i d e < ρ f l ui d e le solide monte et lorsqu’il émerge à la surface du liquide, la
poussée d’Archimède diminue et le solide flotte.
z Si le solide est hétérogène, le centre de gravité G et le centre de poussé A sont des points
distincts.
– Si ρ̄ sol i d e = ρ f l ui d e le solide va tourner et osciller jusqu’à ce que les points G et
A soient sur une même verticale, G étant en dessous de A (condition d’équilibre
stable).
– Si ρ̄ sol i d e > ρ f l ui d e le solide tourne et descend à la fois.
– Si ρ̄ sol i d e < ρ f l ui d e le solide tourne et monte à la fois et vient finalement flotter à la
surface.

2.6 Corps flottants

Les résultats établis précédemment s’appliquent intégralement si la surface S limite un
corps solide en équilibre, immergé complètement ou flottant. C’est le « principe d’Archi-
mède » : le corps solide subit une poussée égale et opposée au poids du volume de fluide
déplacé. Elle passe par le centre de gravité du volume déplacé, ou centre de poussée.

2.6.1 Définitions
Un flotteur est un solide de forme quelconque, généralement fermé, en équilibre dans
un liquide. L’équilibre n’est possible que si le poids du flotteur est inférieur ou égal au poids
du volume de liquide qu’il peut déplacer.

Le plan de flottaison est le plan de la sur-

face libre du liquide. Ce plan coupe le flotteur
suivant une surface appelée flottaison dont le
contour s’appelle ligne de flottaison.

La carène est le volume immergé du flotteur

(volume situé au dessous du plan de flottai-
son).
La poussée du liquide passe par le centre de gravité du volume de la carène. C’est le
centre de poussée ou centre de carène.

©Dr OUEDRAOGO 25 [email protected]

Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

2.6.2 Différentes positions d’équilibre (cas d’un bateau)

Imaginons un bateau de centre d’inertie G ( G est un point lié au bateau), et notons G x ,
G y , G z trois directions cartésiennes (figure ci-dessus pour un Bateau partiellement immergé
dans l’eau). L’axe G z est verticale et G y est dans la direction d’avancement du bateau.
Les mouvements de rotation autour des axes G x , G y , G z sont classiquement appelés
«Tangage», «Roulis», «Lacet».
Les mouvements de translation suivant les directions G x , G y sont associés aux termes de
«Cavalement», «Embardée», «Pilonnement».

Appliquons le théorème d’Archimède au bateau totalement immergé dans l’air et l’eau,

mais partiellement immergé dans l’eau.
la résultante de la poussée d’Archimède est
−
→ →
−
P a = − (m eau + m ai r ) g

où m eau est la masse d’eau déplacée et m ai r est la masse d’air déplacé.

La masse d’air déplacé est très petite devant la masse d’eau déplacée. En effet la masse
volumique de l’air qui est 1, 29 kg .m −3 est très petite devant celle, 1000 kg .m −3 , de l’eau.
Pour le bateau, la résultante de la poussée d’Archimède est donc :
−
→ →
− →
−
P a = − (m eau + m ai r ) g ≈ −m eau g

Considérons le bateau comme un corps allongé cylindrique présentant un plan de sy-

métrie. On étudie ce qui se passe dans un plan de section droite. Si le centre d’inertie G du
bateau est au-dessous du centre de poussée P 0 , il y a stabilité (comme pour les corps com-
plètement immergés). Mais il peut y avoir également stabilité avec le centre d’inertie G au-
dessus du centre de poussée P 0 . On ne donne ici qu’une présentation intuitive du résultat.

©Dr OUEDRAOGO 26 [email protected]

Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

Pour une démonstration mathématique, on peut introduire le centre de courbure I (appelé

métacentre) de la ligne de carène en P 0 . Ce sont les positions relatives de P 0 , I et G qui per-
mettent de conclure à la stabilité ou non. Si G est au-dessous de I , le bateau est stable ; si G
est au-dessus de I , le bateau est instable.

2.7 Application aux fluides compressibles

Les systèmes gazeux sont une illustration simple des cas dans lesquels la compressibilité
du fluide ne peut être négligée. La masse volumique dépendant directement de la pression,
celle-ci n’est évidemment plus indépendante de l’altitude z.

Exercice
Soit un gaz parfait d’équation d’État pV = nRT .
1. A partir de l’équation fondamentale de la statique des fluide sous sa forme différen-
tielle, exprimer la pression p à l’altitude z en fonction de la masse molaire M du gaz, la
constante des gaz parfait R, l’accélération de la pesanteur g , la pression p 0 à l’altitude
z = z0 .
2. Déterminer la hauteur critique à laquelle la variation de la pression p par rapport à la
pression p 0 = p(z 0 = 0) = p at m est de 0, 1% à la température de 25o C .

©Dr OUEDRAOGO 27 [email protected]

CHAPITRE 3
DYNAMIQUE DES FLUIDES PARFAITS INCOMPRESSIBLES

Dans ce chapitre, nous allons étudier la dynamique les fluides, c’est-à-dire de considérer
non seulement des mouvements, mais aussi des efforts, une loi de comportement et d’ap-
pliquer le principe fondamental de la dynamique.
Dans les problèmes de mécanique des fluides, il importe souvent de comprendre la struc-
ture de l’écoulement et de calculer les répartitions de variables comme la pression, la vitesse,
la température et la masse volumique dans le milieu fluide. Dans les applications technolo-
giques, le fluide s’écoule autour de corps solides ou dans des conduites ; la connaissance des
distributions de pression et de vitesses au voisinage des parois est particulièrement utile. On
cherche aussi dans beaucoup de cas à déterminer les contraintes ou des quantités intégrées
comme la force et le moment qui s’exercent sur un corps solide fixe ou en mouvement.
Nous nous limiterons dans ce chapitre aux mouvements des fluides parfaits, c’est-à-dire
sans frottement (fluides non visqueux) et sans échange de chaleur, ou encore un fluide dont
les transformations sont thermodynamiquement réversibles. Nous étudierons tout particu-
lièrement le cas de fluides incompressibles.

3.1 Équation de continuité

3.1.1 Débits massique et volumique
Le débit massique à travers une section S est la quantité de fluide (fluide de masse volu-
mique ρ animé d’une vitesse d’écoulement → −v qui traverse la section S par unité de temps. À
travers la surface S, le débit massique de fluide, noté Q m ou ṁ en kg .s −1 , est donné par :

28
Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

Ï
ṁ = Q m = ρ→
−
v .→
−
n dS
S
Les débits sont généralement comptés positi-
vement dans le sens de l’écoulement : →−
n est
donc orienté dans le sens de l’écoulement.

Le débit volumique à travers la surface S, noté QV ou V̇ en m 3 .s −1 est donné par :

Ï
V̇ = QV = →
−
v .→
−
n dS
S

Le théorème de Green-Ostrogradski (voir Calcul Vectoriel) permet de transformer l’inté-

grale de surface en intégrale de volume :
Ï Ñ
ṁ = Q m = →
− →
−
ρ v . n dS = d i v ρ→
¡ −¢
v dV
S V

3.1.2 Bilans de masse : équation de continuité

L’équation de continuité traduit le principe de conservation de la masse : la variation
de masse pendant un temps dt d’un élément de volume fluide doit être égale à la somme des
masses de fluide entrant diminuée de celle du fluide sortant.
L’équation de continuité peut s’établir à partir d’un bilan global. Lorsqu’on effectue un
bilan non plus sur un volume élémentaire mais sur un volume fini de fluide, les quantités
étudiées apparaissent sous forme d’intégrales de volume. On parle de bilan global.
Ainsi la masse d’un volume de fluide contenue dans le volume V est égale à la somme
les masses entrantes (soit le débit massique Q m ) pendant dt à travers la surface S limitant le
volume V : Ñ Ñ
d i v ρ→
¡ −¢
m= ρdV = − Q m d t = −d t v dV
Ñ V V
d
Ñ
d i v ρ→
¡ −¢
ρdV = −
¡ ¢
v dV
dt V V
puisque le volume de contrôle V est fixe :

∂ρ
Ñ · ¸
¡ →
−
+ d i v ρ v dV = 0
¢
V ∂t
Cette équation corresponde à la loi de conservation de la masse sous forme globale.
Cette relation étant valable quelque soit V, elle entraîne l’équation de continuité, qui est une
équation locale valable en tout point du fluide :

∂ρ
+ d i v ρ→
¡ −¢
v =0
∂t

©Dr OUEDRAOGO 29 [email protected]

Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

∂ρ
La divergence est négative si le fluide est en phase de compression ( ∂t > 0 ), et la diver-
∂ρ
gence est positive si le fluide est en phase de dilatation ( ∂t < 0 ). S’il n’y a ni compression, ni
dilatation (fluide isovolume) la divergence est nulle.

3.1.3 Équation de continuité pour divers écoulements

Écoulement stationnaire
Un écoulement est dit permanent (ou stationnaire) lorsque toutes les grandeurs caracté-
ristiques du mouvement sont invariables dans le temps (vitesse, masse volumique, pression,
température, etc.), ce qui se traduit symboliquement par :

∂
Mouvement stationnaire ⇔ =0
∂t
Sur le plan cinématique, le champ de vitesse ne varie pas dans le temps. Dans ce cas :
– les lignes de courant sont fixes dans l’espace
– les trajectoires coïncident avec les lignes de courant.
En effet, toutes les particules passées en M à un instant donné ont en ce point le même
vecteur vitesse et donc la même trajectoire. La description et la prévision de l’écoulement
sont beaucoup plus simples et de ce fait mieux connues.
Les écoulements qui ne sont pas stationnaires sont appelés tout naturellement instation-
naires. Nous ne les traiterons pas ici. Leur description est mathématiquement complexe et
dans la pratique on recherchera un éventuel régime pseudo-stationnaire en définissant, si
elle existe, une période T au bout de laquelle les paramètres de l’écoulement reprennent des
valeurs identiques. À défaut, ils seront définis par leur valeur moyenne sur la durée T. C’est
le cas des écoulements turbulents dont nous parlerons plus tard.

Conséquences
z En régime permanent l’équation de continuité devient : d i v ρ→
¡ −¢
v =0
z Si l’écoulement est permanent, alors le débit massique est conservé à travers toute
section droite du tube de courant :

ṁ (S 1 ) = ṁ (S 2 )

Toutes les particules intérieures à ce tube

à une certaine époque y resteront toujours
intérieures du fait que leurs trajectoires sta-
tionnaires ne peuvent à aucun moment
percer la surface latérale du tube de cou-
rant.

NB : Si le fluide est incompressible, alors le débit volumique est aussi conservé. Pour
la plupart des applications en hydraulique où l’eau est considérée comme un milieu

©Dr OUEDRAOGO 30 [email protected]

Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

incompressible, le débit est pris au sens volumique. En revanche, en aéraulique la den-

sité du gaz varie, seul le débit massique est significatif du flux de matière.

Écoulement incompressible
Un écoulement est incompressible si le volume de toute particule de fluide reste constant
au cours de son mouvement. Les particules de fluide ayant une masse constante par construc-
tion, leur masse volumique est donc constante au cours de leur écoulement : ρ = c t e.

Conséquences
z Dans le cas d’un écoulement incompressible, l’équation de conservation de la masse
devient :
di v →
¡− ¢
v =0

z Le débit volumique est le même à travers tout tube de courant (de même que le débit
massique).
Ainsi pour un fluide incompressible en écoulement permanent (unidirectionnel) :

ṁ (S 1 ) = ṁ (S 2 )
Ï Ï
→
− →
−
ρ v . n dS = ρ→ −v .→
−
n dS
S1 S2
ρv 1 S 1 = ρv 2 S 2
v1S1 = v2S2
Cette dernière relation donne le rapport entre les vitesses d’écoulements et les sections
d’une conduite :
v1 S2
=
v2 S1

3.2 Équation d’Euler

La quantité de mouvement d’une particule matérielle de masse m et de vitesse instanta-
née →
−
v est m →
−
v.
Son accélération est →
−
a = d→
−
v /d t .

3.2.1 Principe fondamental de la dynamique appliqué à un fluide parfait

Le principe fondamental de la dynamique exprime que la variation dans le temps de la
quantité de mouvement d’un système matériel à nombre constant de particules est égale à
la somme des forces extérieures qui lui sont appliquées.
Pour un fluide, on va appliquer la loi fondamentale de la dynamique à un domaine ma-
tériel que l’on suit dans son mouvement. Cela se traduira par : La dérivée particulaire de la
quantité de mouvement est égale à la résultante des forces extérieures appliquées au domaine.

©Dr OUEDRAOGO 31 [email protected]

Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

Considérons un volume matériel fini de fluide V (t ) regroupant un ensemble donné de

particules fluides d’un écoulement quelconque, de frontière S(t ) et de normale extérieure
unitaire →−
n.
Au cours du mouvement, ce domaine va changer de position, de forme, de dimension,
mais restera constitué des mêmes particules fluides (il n’échange pas de masse avec le reste
du fluide).
La variation dans le temps de la résultante des quantités de mouvement du volume V
s’écrit :

d d→
−
v
µÑ ¶ Ñ
→
−
ρ v dV = ρ dV (3.1)
dt V (t ) V (t ) dt
Inventaires des forces exercées sur V :
– les forces de volume s’exerçant à distance en tout point M du volume : la densité des
→
−
forces par unité de volume exercées par le champ de gravité est ρ g ; la résultante sur
V , c’est-à-dire le poids, vaut : Ñ
→
−
ρ g dV (3.2)
V
– les forces de surface qui sont transmises par le fluide extérieur au domaine en tout
point P de la surface de contact S limitant ce domaine : la densité des forces par unité
de surface est −p →
−
n ; sa résultante, la force de pression, vaut :
Ï
−p →
−
n dS (3.3)
S
La variation dans le temps de la quantité de mouvement s’exprime alors par :

d→
−
v
Ñ Ï Ñ
→
−
ρ dV = −p →
−
n dS + ρ g dV (3.4)
V (t ) dt S V

ou bien encore, en appliquant le théorème d’Ostrogradski

d→
−
v −−−−→
Ñ µ ¶
→
−
ρ + g r ad p − ρ g dV = 0 (3.5)
V (t ) dt
La relation est valable quel que soit le domaine fluide V . En outre, nous pouvons sans
problème supposer l’intégrand continu. Alors l’intégrand est identiquement nul, et il vient
l’expression locale du bilan de quantité de mouvement, nommée équation d’Euler :

d→−v −−−−→ →
−
ρ= −g r ad p + ρ g (3.6)
dt
C’est une équation locale dont les différents termes sont homogènes à une force par unité
de volume et s’expriment donc en N .m −3 .

©Dr OUEDRAOGO 32 [email protected]

Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

3.2.2 Application à un écoulement unidirectionnel d’un fluide incompressible

Lorsque l’écoulement est rectiligne, on parle d’écoulement unidirectionnel. Il est tel que
la vitesse est en tout point parallèle à une direction unique fixe.
Les trajectoires sont des droites parallèles à l’axe x de l’écoulement. C’est le cas, par
exemple, d’un écoulement dans une canalisation cylindrique, ou d’un écoulement libre dans
les zones non perturbées par un obstacle.

L’équation d’Euler s’écrit pour un fluide incompressible dans le champ de pesanteur sa-
→
− −−−−→
chant que g = −g g r ad z :

d→
−
v −−−−→ ¡
ρ = −g r ad p + ρg z
¢
(3.7)
dt
où z est l’altitude (comptée positivement suivant la verticale ascendante) et → −v (v, 0, 0)
→
−
pour un écoulement unidimensionnel (les composantes de v dans le plan perpendiculaire
à la direction de cet écoulement sont nulles).
L’hypothèse d’écoulement rectiligne implique la colinéarité de → −
v et d →
−v /d t . donc de → −v
−−−−→ ¡
et g r ad p + ρg z . Ainsi la composante normale cette dernière expression est nulle.
¢

Par conséquent, l’expression (p + ρg z) reste constante quand on se déplace normale-

ment aux trajectoires ; la répartition des pressions est hydrostatique.
Ce résultat peut sembler évident car lorsque la vitesse du fluide n’intervient pas, l’hydro-
statique remplace la dynamique.
Si, de plus, l’écoulement est permanent, les composantes de la quantité de mouvement
s’écrivent :
∂v ∂ ¡
ρ p + ρg z
¢
=−
∂x ∂x
∂ ¡
p + ρg z
¢
0=−
∂y
∂ ¡
p + ρg z
¢
0=−
∂z
Compte tenu de l’équation de conservation de la masse :

©Dr OUEDRAOGO 33 [email protected]

Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

∂v
d i v→
−
v = =0
∂x
on conclut que :
−−−−→ ¡
g r ad p + ρg z = 0
¢

soit p + ρg z = c t e dans tout l’écoulement unidirectionnel.

3.3 Équation de Bernouilli

L’équation de Bernouilli traduit le bilan énergétique du fluide.

3.3.1 Hypothèses - Énoncé

Sous les hypothèses suivantes :
– le fluide est incompressible (ρ = c st e),
– le fluide est parfait,
→−
– les seules forces extérieures sont les efforts de pesanteur avec l’accélération g = −g →
−
ez
constante,
– l’écoulement est stationnaire (la vitesse →−v et la pression p exprimées en variables eu-
lériennes ne dépendent pas explicitement du temps),
on a, v désignant le module du vecteur vitesse → −
v :

v2 p
+ + g z = cte
2 ρ

sur chaque ligne de courant (ou chaque trajectoire).

→
−
NB : Il est à remarquer, que si le fluide est au repos, alors →
−
v = 0 , et on retrouve la loi de
l’hydrostatique.

3.3.2 Démonstration
On peut partir de l’équation d’Euler ou d’un bilan global (conservation de la masse,
conservation de l’énergie).

Cas d’un écoulement sans travail

On considère une veine d’un fluide incompressible de masse volumique ρ animée d’un
écoulement permanent.
On désigne, à partir de la figure ci-dessous, par :
– S 1 et S 2 respectivement la section d’entrée et la section de sortie du fluide à l’instant t,
0 0 0
– S 1 et S 2 respectivement les sections d’entrée et de sortie du fluide à l’instant t = (t +
d t ),

©Dr OUEDRAOGO 34 [email protected]

Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

– −→ et −
v 1
→ les vecteurs vitesse d’écoulement respectivement à travers les sections S et S
v 2 1 2
de la veine.
– d x 1 et d x 2 respectivement les déplacements des sections S 1 et S 2 pendant l’intervalle
de temps d t ,
0
– d m 1 : masse élémentaire entrante comprise entre les sections S 1 et S 1 ,
0
– d m 2 : masse élémentaire sortante comprise entre les sections S 2 et S 2 ,
– M : masse comprise entre S 1 et S 2 ,
0
– dV1 : volume élémentaire entrant compris entre les sections S 1 et S 1 ,
0
– dV2 : volume élémentaire sortant compris entre les sections S 2 et S 2 ,
→
−
On considère un axe Z vertical dirigé vers le haut. On note Z1 , Z2 et Z respectivement
les altitudes des centres de gravité des masses d m 1 , d m 2 et M .
On désigne par F 1 et F 2 respectivement les normes des forces de pression du fluide agis-
sant au niveau des sections S 1 et S 2 .

A l’instant t le fluide de masse (d m 1 + M ) est compris entre S 1 et S 2 . Son énergie méca-

nique est :

1 S2 d mv 2 ¡
Z
E m = E c + E p = d m 1 v 12 +
¢
0
+ d m 1 g Z1 + M g Z
2 S1 2
0 0 0
A l’instant t = (t +d t ) le fluide de masse (M +d m 2 ) est compris entre S 1 et S 2 . Son énergie
mécanique est :

S2 d mv 2 1
Z
0 0 0
+ d m 2 v 22 + d m 2 g Z2 + M g Z
¡ ¢
Em = Ec + E p = 0
S1 2 2

©Dr OUEDRAOGO 35 [email protected]

Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

0
On applique le théorème de l’énergie mécanique au fluide entre t et t : la variation de
l’énergie mécanique est égale à la somme des travaux des forces extérieures.

0
³→
−´
Em − Em = W F

1 1
d m 2 g Z2 − d m 1 g Z1 + d m 2 v 22 − d m 1 v 12 = F 1 d x 1 − F 2 d x 2 = p 1 S 1 d x 1 − p 2 S 2 d x 2
2 2

1 1 p1 p2
d m 2 g Z2 − d m 1 g Z1 + d m 2 v 22 − d m 1 v 12 = p 1 dV1 − p 2 dV2 = d m1 − d m2
2 2 ρ1 ρ2

Par conservation de la masse : d m 1 = d m 2 = d m et puisque le fluide est incompressible :

ρ 1 = ρ 2 = ρ, On aboutie à l’équation de Bernoulli :

v 12 − v 22 p1 − p2
+ + g (Z1 − Z2 ) = 0
2 ρ
L’équation ci-dessus peut être mise sous la forme :

v 12 − v 22 p1 − p2
+ + (Z1 − Z2 ) = 0 ⇔
2g ρg

v 12 p1 v 2 p2
+ + Z1 = 2 + + Z2
2g ρg 2g ρg
donc sur chaque ligne de courant d’un fluide en écoulement :

v2 p
+ + Z = cte (3.8)
2g ρg
ou
v2 p
+ + g Z = cte (3.9)
2 ρ
ou

v2
ρ + p + ρg Z = c t e (3.10)
2

Interprétation graphique de l’équation de Bernouilli

L’expression 3.8 est homogène à une à une hauteur de liquide :
z Z est l’altitude
p
z ρg + Z est appelée la charge piézométrique (ou hauteur piézométrique),
v2 p
z 2g
+ ρg + Z = H est la charge totale.

©Dr OUEDRAOGO 36 [email protected]

Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

Le théorème de Bernoulli peut alors s’interpréter graphiquement à partir des évolutions des
différentes hauteurs le long du circuit. Comme l’illustre la figure, quand on suit l’unité de
poids de fluide dans son mouvement le long de la trajectoire, on peut tracer trois lignes dif-
férentes :
z la ligne d’altitude qui représente la trajectoire du fluide,
z la ligne piézométrique, distante de la trajectoire de la quantité p/ρg ,
z la ligne de charge, obtenue en ajoutant v 2 /2g à la ligne piézométrique. Le théorème de
Bernoulli conduit à une ligne de charge horizontale (charge totale = cte). Il n’y a pas de
perte de charge dans l’écoulement d’un fluide parfait.

Interprétation en termes de pression de l’équation de Bernouilli

L’expression 3.10 est homogène à une à unité de pression (énergie par unité de volume) :
z p est la pression statique (elle existe même s’il n’y a pas de mouvement),
z p + ρg Z est la pression motrice ; elle génère le mouvement (ρg z est la pression de pe-
santeur),
2
z p + ρg Z + ρ v2 = p t est la pression totale.
L’équation de Bernoulli montre alors que la pression totale reste constante le long d’une
même ligne de courant (il n’y a pas de perte de pression dans l’écoulement d’un fluide parfait).

3.3.3 Équation de Bernouilli d’un écoulement avec échange de travail

Reprenons le schéma de la veine fluide du paragraphe précédent avec les mêmes nota-
tions et les mêmes hypothèses. On suppose en plus qu’une machine hydraulique est pla-
cée entre les sections S 1 et S 2 . Cette machine est caractérisée par une puissance nette P net
échangée avec le fluide, une puissance sur l’arbre P a et un certain rendement η. Cette ma-
chine peut être soit une turbine soit une pompe.
Le rendement est donné par :

©Dr OUEDRAOGO 37 [email protected]

Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

– pour une pompe

P net
η=
Pa
– pour une turbine
Pa
η=
P net
0
Entre les instant t et t = (t + d t ), le fluide a échange un travail net Wnet = P net d t avec
la machine hydraulique. Wnet est supposé positif s’il s’agit d’une pompe et négatif s’il s’agit
d’une turbine.

En considérant cette fois ci le travail de la machine hydraulique, l’équation de Bernouilli

se met sous la forme :

v 12 − v 22 p1 − p2 P net
+ + g (Z1 − Z2 ) =
2 ρ qm

3.4 Application du Théorème de Bernouilli

3.4.1 Écoulement par orifice : formule de Torricelli
Considérons un réservoir de grande section, rempli d’un liquide qui s’écoule à travers un
orifice de section faible (devant celle du réservoir) et situé à une hauteur sous la surface libre

©Dr OUEDRAOGO 38 [email protected]

Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

(voir figure ci-dessous). La pression atmosphérique s’exerce à la fois sur la surface libre et
sur le jet à la sortie de l’orifice. En supposant que le liquide est incompressible et non vis-
queux, il est possible d’utiliser l’équation de Bernoulli pour déterminer la vitesse de vidange
du réservoir.

Considérons un point A de la surface libre et un point B situé dans le jet. Il existe une ligne
de courant passant par ces deux points et le long de laquelle on peut appliquer l’équation
de Bernoulli pour déterminer l’expression de la vitesse en B qui est la formule deTorricelli
compte tenu du rapport de section entre le réservoir et l’orifice, on peut y négliger v A devant
v B . Par ailleurs, le point A étant situé sur la surface libre, on a donc p A = p at m .
p
vB = 2g H

La vitesse d’écoulement du liquide est alors indépendante de sa masse volumique ; le

seul paramètre dont dépend la vitesse de vidange est donc la hauteur de liquide séparant
l’orifice de la surface libre.
L’évaluation du débit de vidange nécessite de prendre en compte la contraction du jet.
Comme le montre l’encart de la figure ci-dessus, la vitesse évaluée précédemment corres-
pond à une section plus faible que celle de l’orifice. En conséquence, le débit volumique de
vidange s’obtient en calculant :

Q B = v B σ = v B SC c

σ est la section du jet où les lignes de courant peuvent être considérées rectilignes et
parallèles. On peut ainsi définir un coefficient de contraction C c = σS , lequel dépend essen-
tiellement du type de paroi ainsi que du profil de l’orifice dans la paroi.
La figure ci-dessous donne de manière non exhaustive quelques valeurs typiques de ce
coefficient de contraction.

©Dr OUEDRAOGO 39 [email protected]

Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

3.4.2 Effet Venturi

Associé au principe de conservation du débit d’un écoulement en conduite, la conserva-
tion de la pression totale le long des lignes de courant (équation de Bernoulli) conduit natu-
rellement à observer l’effet Venturi : un élargissement (rétrécissement) local de la conduite
provoque localement une surpression (dépression). Cet effet est à la base d’un dispositif
simple permettant la mesure du débit d’un écoulement. Il suffit d’aménager un rétrécisse-
ment le long d’une conduite cylindrique dans laquelle s’écoule un fluide supposé parfait et
incompressible.

TAF : Utiliser l’équation de Bernouilli sur les données de la figure ci-dessus et déterminer
le débit Q v en fonction de g , ∆h et α = SS BA .

Remarque Ce dispositif est couramment utilisé comme élément intégré aux circuits hy-
drauliques et permet une mesure simple du débit d’écoulement. Il génère toutefois des pertes
de charge : c’est donc un dispositif de mesure qui dissipe de l’énergie. La perte de charge peut
facilement être évaluée en plaçant en amont du rétrécissement un troisième tube manomé-
trique. Compte tenu du principe de conservation du débit, à section égale, la vitesse y est la
même que v A . L’application de l’équation de Bernoulli entre l’amont et l’aval impose donc
d’avoir la même pression si le fluide est parfait (non visqueux). Or, une mesure réelle montre
que la pression en aval est inférieure à celle en amont : cela signifie simplement que lorsque

©Dr OUEDRAOGO 40 [email protected]

Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

les frottements ne peuvent être négligés (viscosité non nulle) l’équation de Bernoulli établie
pour un fluide parfait n’est pas valable ; le rétrécissement génère en pratique une chute de la
pression totale que l’on qualifiera de perte de charge singulière (voir chapitre à venir sur les
pertes de charge).

3.4.3 Tube de Pitot

Dans son principe, il s’agit d’un dispositif extrêmement simple qui permet une mesure de
la vitesse d’écoulement d’un fluide. L’objet présente une forme profilée, est creux afin d’être
rempli du fluide dans lequel il est immergé, et doit être muni de deux prises de pression
(tubes manométriques). Comme le montre le schéma de la figure ci-dessous, l’un des deux
tubes manométriques est relié au front d’attaque de l’objet (point d’arrêt caractérisé par une
vitesse d’écoulement nulle), alors que l’autre est en prise avec le fluide statique remplissant
l’objet.

En supposant que le fluide est non visqueux, incompressible et que l’écoulement est sta-
tionnaire et uniforme en amont de l’objet, on va pouvoir identifier un certain nombre de
lignes de courant et y appliquer l’équation de Bernoulli. On supposera par ailleurs que toutes
ces lignes de courant sont approximativement à la même altitude. Le calcul (à faire par les
étudiants) donne avec les données de la figure ci-dessus :
q
U= 2g ∆h
La fiabilité de cette mesure repose essentiellement sur la qualité de la prise de pression
statique (point B) : la paroi doit être lisse et sans imperfection au voisinage de l’orifice (évi-
tant ainsi une surévaluation ou sous-évaluation de la pression). Bien qu’extrêmement simple
et rudimentaire, ce dispositif équipe notamment les avions modernes pour évaluer leur vi-
tesse de vol.

3.4.4 Phénomène de cavitation

Compte tenu de l’équation de Bernoulli, lorsqu’au sein d’un écoulement la vitesse atteint
localement des valeurs élevées, la pression du liquide chute et peut tomber sous le seuil de
pression de vapeur saturante. Dans ce cas, il y a localement transition de phase et vaporisa-
tion du liquide. Les bulles de vapeur ainsi formées perturbent l’écoulement et nuisent géné-
ralement au bon fonctionnement des dispositifs dans lesquels elles apparaissent (pompes,

©Dr OUEDRAOGO 41 [email protected]

Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

hélices propulsives...). Notamment, la migration de ces bulles vers des régions de plus haute
pression provoque leur implosion : le choc mécanique ainsi produit constitue un facteur
d’endommagement prématuré des hélices.
A 20o C la pression de vapeur saturante de l’eau est de2, 34 kPa.
La cavitation peut également être générée par la propagation d’une onde acoustique (gé-
néralement des ultrasons de puissance) ; on parle alors plutôt de cavitation acoustique, par
distinction avec la cavitation hydrodynamique décrite précédemment.

Conclusion
Les lois et les équations établies ci-dessus ont un intérêt pratique considérable du mo-
ment ou elles permettent de comprendre le principe de fonctionnement de beaucoup d’ins-
truments de mesure de débits tels que le tube de Pitot, le tube de Venturi et le diaphragme,
etc. Réservées aux fluides incompressibles, ces lois et équations peuvent être employées
dans certains cas particulier pour les fluides compressibles à faible variation de pression.
Une telle variation existe dans plusieurs cas pratiques. Cependant, lorsqu’on veut prendre
en considération la compressibilité dans les calculs, il est nécessaire d’employer les formules
appropriées.

3.5 Théorème d’Euler

3.5.1 Hypothèses - Énoncé
Sous les hypothèses suivantes :
– le fluide est incompressible (ρ = c st e),
– le fluide est parfait,
→
−
– les seules forces extérieures sont les efforts de pesanteur avec l’accélération g = −g → −
ez
constante,
– l’écoulement est stationnaire (la vitesse →−v et la pression p exprimées en variables eu-
lériennes ne dépendent pas explicitement du temps), de plus, → −
v et p sont de classe C 1
sur le domaine fluide D
on a, la loi de la conservation de la masse étant satisfaite :
Ï Ï
→
− ¡→
− →
− →
−
ρ v v • n dS = −p →−
¢
n d S + m(D) g (3.11)
S S

pour tout domaine fluide D limité par la surface S, le vecteur → −

n étant unitaire normal à
S et dirigé vers l’extérieur de S, et m(D) étant la masse de fluide contenu dans D.
Le domaine D est appelé volume de contrôle et la surface S surface de contrôle Ce même
théorème est connu aussi sous le nom de Théorème des quantités de mouvement.
Autrement dit :

©Dr OUEDRAOGO 42 [email protected]

Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

F IGURE 3.1 – Volume de contrôle, surface de contrôle

P→
−
La résultante ( F ext ) des actions mécaniques extérieures exercées sur un fluide isolé (fluide
contenu dans une enveloppe limitée par S 1 et S 2 ) est égale à la variation de la quantité de
mouvement du fluide qui entre en S 1 à une vitesse −
→ et sort par S à une vitesse −
v 1 2
→.
v 2

X→
−
F ext = q m −
¡→ − →¢
v2 − v 1 (3.12)

Les forces extérieures sont ici les forces de surface dont la résultante est la force de pres-
sion et les forces de volume dont la résultante est le poids de la masse m(D) de fluide contenu
dans le domaine D.
Ï Ï Ñ
~ext = →
− →
− →− →
−
ρ g dV
X
F −p n d S + m g = −p n d S + (3.13)
S S V

Démonstration
On utilise l’équation d’Euler avec les hypothèses ci-dessus-citées.

Remarque
Il est à remarquer que les deux intégrales figurant dans la relation 3.11 sont des intégrales
de surface. Pour appliquer le théorème d’Euler, il suffit de connaître → −v et p sur la surface S
de contrôle, et non en les points intérieurs à D. Pour appliquer ce théorème, il est inutile de
chercher à déterminer → −
v et p dans l’intérieur de D. Cette remarque est très importante pour
résoudre des problèmes dans les applications qui suivent.

©Dr OUEDRAOGO 43 [email protected]

Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

3.5.2 Application
Une application directe du théorème d’Euler est l’évaluation des forces exercées par les
jets d’eau. Celles-ci sont exploitées dans divers domaines : production de l’énergie électrique
à partir de l’énergie hydraulique grâce aux turbines, coupe des matériaux, etc. Ce théorème
permet de déterminer les efforts exercés par le fluide en mouvement sur les objets qui les
environnent.

– Efforts sur une canalisation hydraulique

– Efforts sur un embout (Poussée d’un réacteur)
– Efforts de pression exercées par un jet sur une plaque plane
– etc

©Dr OUEDRAOGO 44 [email protected]

CHAPITRE 4
DYNAMIQUE DES FLUIDES RÉELS INCOMPRESSIBLES

Dans le chapitre précédent nous avons supposé que le fluide était parfait pour appliquer
l’équation de conservation de l’énergie. L’écoulement d’un fluide réel est plus complexe que
celui d’un fluide idéal. En effet, il existe des forces de frottement, dues à la viscosité du fluide,
qui s’exercent entre les particules de fluide et les parois, ainsi qu’entre les particules elles-
mêmes. Pour résoudre un problème d’écoulement d’un fluide réel, on fait appel à des ré-
sultats expérimentaux, en particulier ceux de l’ingénieur et physicien britannique Osborne
Reynolds. Une méthode simplifiée de calcul des pertes de charge basée sur ces résultats ex-
périmentaux est proposée. Elle est indispensable pour le dimensionnement des diverses ins-
tallations hydrauliques (de pompage, de turbines, de machines hydrauliques et thermiques
dans lesquelles est véhiculé un fluide réel, etc.)
Un fluide est dit réel si, pendant son mouvement, les forces de contact ne sont pas per-
pendiculaires aux éléments de surface sur lesquelles elles s’exercent (elles possèdent donc
des composantes tangentielles qui s’opposent au glissement des couches fluides les unes sur
les autres). Cette résistance est caractérisée par la viscosité.

4.1 Régime d’écoulement - Nombre de Reynolds

Les expériences réalisées par Reynolds en 1883 lors de l’écoulement d’un liquide dans
une conduite cylindrique rectiligne dans laquelle arrive également un filet de liquide coloré,
ont montré l’existence de deux régimes d’écoulement : régime laminaire et régime turbu-
lent :
z Régime laminaire : Les filets fluides sont des lignes régulières, sensiblement parallèles
entre elles.

45
Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

z Régime turbulent : Les filets fluides s’enchevêtrent, s’enroulent sur eux-mêmes.

Des études plus fines ont montré qu’il existe encore une subdivision entre :
– les écoulements turbulents lisses et
– les écoulements turbulents rugueux.
La limite entre ces différents types d’écoulements est évidemment difficile à appréhen-
der.
En utilisant divers fluides à viscosités différentes, en faisant varier le débit et le diamètre
de la canalisation, 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 donné par l’expression suivante :

vd
Re =
ν
où v est la vitesse moyenne d’écoulement à travers la section considérée en (m/s)
d : Diamètre de la conduite ou largeur de la veine fluide en (m).
ν : Viscosité cinématique du fluide (m 2 /s).

Résultats empirique à titre indicatif :

z Si R e < 2000 l’écoulement est laminaire
z Si R e > 2000 l’écoulement est turbulent :
– Lisse si 2000 < R e < 105
– Rugueux si Re > 105

4.2 Pertes de charges

4.2.1 Définition
Considérons un écoulement entre deux points (1) et (2) d’un fluide réel dans une conduite,
tel que entre les points (1) et (2) il n’y ait pas de machine hydraulique. Reprenons le schéma
de la veine fluide du chapitre précédent avec les mêmes notations et les hypothèses sui-
vantes :
– Le fluide est réel et incompressible : cela suppose l’existence de forces élémentaire de
frottement visqueux d τ qui contribue dans l’équation de bilan par un travail négatif et
donner naissance à des pertes de charges.
– L’écoulement est permanent.

©Dr OUEDRAOGO 46 [email protected]

Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

On reprend l’application du théorème de l’énergie mécanique en prenant en compte ici

le travail Wd τ des forces de frottement visqueux d τ.
P

L’équation de Bernouilli prend la forme :

v 12 − v 22
P
p1 − p2 Wd τ
+ + g (Z1 − Z2 ) =
2 ρ dm

On défini la perte de charge entre les points (1) et (2) par

P
Wd τ
J 12 =
dm
qui est la perte d’énergie par frottement visqueux par unité de masse qui passe.
L’expression
v 12 − v 22 p 1 − p 2
+ + g (Z1 − Z2 ) = J 12
2 ρ
est telle que l’unité de chaque terme est le joule par kilogramme (J /kg )
On peut réécrire cette expression en des termes homogènes à une hauteur soit :

v 22p2 v 12 p 1 J 12
+ + Z2 = + + Z1 +
2g ρg 2g ρg g
Elle peut être interprétée graphiquement de la manière suivante :

©Dr OUEDRAOGO 47 [email protected]

Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

En l’absence de pertes de charge, la ligne de charge est confondue avec le plan de charge.
Ce plan de charge donne une représentation graphique de la constance tirée de l’équation
de Bernoulli pour un fluide parfait. La perte de charge totale exprimée en hauteur de liquide
depuis le début de l’écoulement, est égale à la distance entre la ligne de charge et le plan de
charge, mesurée sur la verticale passant par le point G 1 . La perte de charge entre deux points
G 1 et G 2 de l’écoulement est donnée par la différence de cote de la ligne de charge sur les
verticales passant par les points précédents.

La perte de charge J 12 peut être due à une perte de charge linéaire et une perte de charge
singulière.

J 12 = J S + J L

Par exemple, dans le circuit représenté dans la figure ci-dessous :

– les tronçons BC , DE , F G, H I et J K sont des coudes de différents angles, donc elles
présentent des pertes de charge singulières.
– Les tronçons AB , C D, E F , G H , I J et K L sont des conduites rectilignes, donc elles pré-
sentent des pertes de charge linéaires.

©Dr OUEDRAOGO 48 [email protected]

Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

4.2.2 Pertes de charge singulières

Quand la conduite subit de brusque variation de section ou de direction, il se produit
des pertes de charges dites singulières, elles sont généralement mesurable et font partie des
caractéristiques de l’installation.
On les exprime par :

v2
J S = −K s
2
où s est l’indice de l’accident de forme de la conduite.
K s : Coefficient (sans unité) de pertes de charge. Il dépend de la nature et de la géométrie
de l’accident de forme.
Les valeurs de K s sont données par les constructeurs dans leurs catalogues.
La figure ci-dessous donne des formules de calcul du coefficient de pertes de charge se-
lon le type de singularité.

©Dr OUEDRAOGO 49 [email protected]

Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

4.2.3 Pertes de charges linéaires

Les pertes de charges linéaires, sont des pertes de charge réparties régulièrement le long
des conduites. En chaque point d’un écoulement permanent, les caractéristiques de l’écou-
lement sont bien définies et ne dépendent pas du temps. La représentation graphique de
l’écoulement prend l’allure ci-dessous.

La vitesse étant constante, la ligne piézométrique et la ligne de charge sont parallèles.

La variation de hauteur piézométrique, évaluée en hauteur de liquide est égale à la perte de
charge linéaire entre les deux points de mesure.
Les pertes de charge linéaires sont proportionnelles à la longueur L(m) de la conduite,
inversement proportionnelles à son diamètre d (m), proportionnelle au carré de la vitesse
débitante v(m/s) du fluide.

v2 L
µ ¶
J L = −λ
2 d
où λ est le coefficient de perte de charge linéaire. Il dépend du régime d’écoulement et
notamment du nombre de Reynolds R e .
z Dans un régime d’écoulement laminaire où R e < 2000, on utilise la formule de Poi-
seuille :
64
λ=
Re
z Dans un régime d’écoulement turbulent lisse (2000 < R e < 105 ) : Formule de Blasius)

λ = 0, 316R e−0,25

z Dans un régime d’écoulement turbulent rugueux (R e > 105 ) : Formule de Blench)

ε
r
λ = 0, 79
d
avec
ε : rugosité de la surface interne de la conduite (mm)
d : diamètre intérieur de la conduite (mm)

©Dr OUEDRAOGO 50 [email protected]

Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

Parfois, on lit la valeur de λ sur un abaque établi par Moody.

Conclusion
Considérons un écoulement entre deux points (1) et (2) d’un fluide réel dans une conduite.
On suppose éventuellement, qu’il existe entre (1) et (2) des machines hydrauliques.
On note :
J 12 : Somme de toutes les pertes de charge, singulière et linéaires entre les sections (1) et
(2).
P n : Puissance mécanique échangé entre le fluide et les machines éventuellement placées
entre (1) et (2).
Le Théorème de Bernoulli prend la forme générale suivante :

1¡ 2 ¢ 1¡ Pn
v 2 − v 12 +
¢
p 2 − p 1 + g (Z2 − Z1 ) = J 12 +
2 ρ qm

©Dr OUEDRAOGO 51 [email protected]

 BIBLIOGRAPHIE

1. Mécanique appliquée : Résistance des matériaux, Mécanique des fluides, Thermo-

dynamique - 2ème édition : Pierre Agati, Frédéric Lerouge, Nicolas Mattéra. Dunod
2002.
2. Mécanique expérimentale des fluides : Raymond Comolet, Jacques Bonnin Dunod.
(2002) .
3. CAPES de Sciences physiques : Tome 1 – Physique – cours et exercices. Nicolas BILLY,
Jean DESBOIS, Marie-Alix DUVAL, Mady ELIAS, Pascal MONCEAU Aude PLASZCZYNSKI,
Michel TOULMONDE. 3e édition BELIN 2004.
4. IIntroduction to Thermal Systems Engineering : Thermodynamics, Fluid Mecha-
nics, and Heat Transfer. : Michael J. Moran, Howard N. Shapiro, Bruce R. Munson,
David P. DeWitt. John Wiley & Sons, Inc. (2002)
5. Fundamentals of Fluid Mechanics, 6th edition : B.R. Munson, D.F. Young, and T.H.
Okiishi. John Wiley & Sons, Inc. (2009) .

52
CHAPITRE 5
ANNEXE 1 : ANALYSE VECTORIELLE

5.1 Notion de champ

En physique, on est souvent amené à traiter des propriétés différentes en tout point de
l’espace, qui apparaissent aussi bien en physique classique qu’en certains domaines de la
physique quantique. Les champs proposent un formalisme intuitif pour modéliser ces phé-
nomènes.
→
−
Pour bien démarrer considérons un vecteur E qui dépend d’une seule variable t, et no-
→
− → − →
−
tons E = E (t ) cette relation de dépendance. Alors E est une application vectorielle de la
variable t . On remarquera que si t varie, alors la norme ou la direction (ou les deux à la fois)
→
−
du vecteur E peuvent varier.
En généralisant ce cas particulier d’application dans le plan R2 , l’espace R3 puis à un es-
pace Rn à n dimensions, on arrive à la définition suivante :

Définition Un champ est une application M : Rn → Kp continue, où n et p sont deux en-

tiers naturels non-nuls et où K désigne le corps des réels ou des complexes.
Il faut remarquer qu’un tel objet est défini en tout point de l’espace. Dans certains cas,
on peut restreindre les champs à une zone seulement de l’espace, mais cela est bien souvent
inutile, voire faux physiquement. Nous traiterons donc uniquement le cas général.
On distingue généralement deux cas particuliers :
– On appelle champ scalaire un champ qui, à tout point de l’espace, associe un réel ou
un complexe.
– On appelle champ vectoriel un champ qui, à tout point de l’espace, associe un vecteur
réel ou complexe.

53
Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

5.2 Gradient
Le premier outil vectoriel, car le plus simple, est l’opérateur gradient. Pour l’introduire,
imaginons que le champ est un champ scalaire réel, et qu ?il représente la hauteur d’un ter-
rain en chaque point. On veut savoir dans quelle direction aller pour monter : l’opérateur
gradient donne cette information.
Dans notre exemple, il s’agit d’un champ scalaire de dimension deux. Considérons indé-
pendamment chaque direction. On peut dériver notre champ par rapport à x : si le champ
est localement décroissant, on aura un nombre négatif ? s’il est croissant, on aura un nombre
positif. Ainsi, un vecteur d’abscisse ce nombre ira dans la direction vers laquelle le champ
croît. Il en est de même selon y. Ainsi le vecteur :

∂M
 
 ∂x 
 
 
 ∂M 
 

∂y
indique la direction vers laquelle le champ M est croissant. Ce vecteur est appelé vecteur
gradient, et généralise la notion de dérivée.

5.2.1 Définition
L’application qui à un champ scalaire M associe un champ vectoriel dont la valeur en
chaque point est donnée, en coordonnées cartésiennes, par :

∂M
 
 ∂x 
 
 
 ∂M 
 
 
 ∂y 
 
 
 
 ∂M 
∂z
−−−→
est appelée opérateur gradient et noté : grad

5.2.2 Propriétés
On distingue les propriétés du vecteur gradient, valeur du gradient d’un champ en un
point :

Propriétés du vecteur gradient

La direction du vecteur gradient indique la « pente » de plus grande inclinaison – on peut
dire aussi qu’elle pointe vers les zones de valeurs plus grandes.

©Dr OUEDRAOGO 54 [email protected]

Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

La norme du vecteur gradient est proportionnelle à la variation : plus celle-ci est violente,
plus le vecteur est long.

Propriétés de l’opérateur gradient

L’opérateur gradient est un opérateur linéaire.
Pour tous champs scalaires A, B, la relation suivante est vraie :
−−−→ ¡ −−−→ −−−→
∀λ, µ ∈ K, grad λA + µB = λgradA + µgradB
¢

En particulier, si le champ scalaire M dépend du temps, la relation suivante est vraie :

∂ ³−−−→ ´ −−−→ ∂
µ ¶
gradM = grad M
∂t ∂t
L’opérateur gradient n’est invariant que par des transformations orthogonales.

Théorème du gradient

a et ~
Soit M un champ scalaire, ~ b deux points de l’espace, L une courbe reliant ces deux
points. Alors :

→
− −−−→
Z
→
−
M ( b ) − M ( a ) = gradM · d→
−
r
L

5.2.3 Exemples
On dit d’une force qu’elle est conservative s’il existe une énergie potentielle E p telle que :
→
− −−−→
F = −grad E p (parce que le travail est nul sur une courbe fermée)
→
− −−−→
Si on note V le potentiel électrique, alors, en électrostatique : E = −gradV

5.3 Divergence
L’opérateur divergence est un outil d’analyse vectorielle qui mesure, pour faire simple,
si un champ vectoriel « rentre » ou « sort » d’une zone de l’espace, comme ce que l’on peut
observer sur un diagramme de lignes de champ. Il donne donc une information très liée aux
sources qui créent le champ.
Comme nous le préciserons, l’opérateur divergence est l’équivalent local de la mesure
d’un flux.

©Dr OUEDRAOGO 55 [email protected]

Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

5.3.1 Définition
~ associe le champ scalaire dont la valeur en tout
L’application qui à un champ vectoriel M
point est :

∂M x ∂M y ∂M z
+ +
∂x ∂y ∂z
~ ) est appelé opérateur divergence et noté : div
(où M i représente la coordonnée i de M

5.3.2 Interprétations physiques

On peut interpréter le signe de la divergence en termes de « sources ». Cela est notam-
ment flagrant en électromagnétisme, où le lien est direct :
→
− ρ
div E = (équation de Maxwell-Gauss)
ε0
Ici, la divergence du champ électrique est liée à la densité de charges, ?. Si les charges
sont positives, la divergence est positive. De même si les charges sont négatives.
Une autre interprétation possible est en termes de « compression » : si le champ repré-
sente le mouvement de particules, une divergence non nulle implique que des particules
se concentrent (ou s’éloignent) d’une zone de l’espace. En mécanique des fluides, un fluide
incompressible est donc caractérisé par : div →−
v =0
Enfin, une dernière interprétation est en termes de « volume » : le volume formé par des
particules dans le flot du champ augmente si la divergence est positive et diminue si elle est
négative.

5.3.3 Propriétés de l’opérateur divergence

Linéarité de l’opérateur divergence
Soit ~ ~ deux champs vectoriels. Alors, pour tous scalaires λ, µ, on a :
A et B
³ →
− →
−´ →
− →
−
div λ A + µ B = λdiv A + µdiv B

~ un champ vectoriel, on a :
De même, si f est une fonction scalaire et M
³ −→´ −−−→ − → −
→
div f M = grad f · M + f div M

Enfin, il est utile de savoir que :

∂ ∂−
µ ¶
−
→ →
div M = div M
∂t ∂t

©Dr OUEDRAOGO 56 [email protected]

Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

5.3.4 Théorème de la divergence

Le théorème de la divergence, ou théorème de Green-Ostrogradski, qui n’est qu’une forme
particulière du théorème de Stokes, trace un lien entre le flux et la divergence.
Soit M~ un champ vectoriel. Soit V un volume de l’espace, délimité par sa surface (fermée)
S. Alors :
ZZZ ³ Z Z
−
→´ −
→ → −
div M dV = ⊂⊃ M · d S
V S
Le théorème de la divergence relie flux sur une surface et divergence sur un volume.

5.4 Rotationnel
Il existe encore un outil vectoriel élémentaire, qui permet de mesurer la « rotation » du
champ.

5.4.1 Définition
~ un champ vectoriel. L’opérateur rotationnel est l’application qui à F
SoitF ~ associe le
−→ →
−
champ vectoriel rot F , qui vérifie :

→
− ³→
− →
−´ −→ →
− → −
∀ X 0, div X 0 ∧ F = −rot F · X 0

En coordonnées cartésiennes on a :

∂F z ∂F y
 
−
 ∂y ∂z 
 
→
− −→ →
−   ∂F x ∂F z 

∀F , rot F =  ∂z − ∂x 
 
 
∂F y
 
∂F x
∂x
− ∂y

©Dr OUEDRAOGO 57 [email protected]

Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

5.4.2 Interprétation physique

Comme son nom l’indique, l’opérateur rotationnel donne une mesure de la « rotation »
du champ. La direction d’un vecteur de ce champ donne l ?axe de rotation, son intensité la
vitesse de rotation autour de cet axe. S’agissant de vecteurs, on ne connait cependant pas le
centre de la rotation. Un champ vectoriel de rotationnel nul est dit « irrotationnel ».
1 −→ −
Par exemple, en mécanique des fluides, on définit le vecteur vorticité par :→−
ω = rot → v.
2

5.4.3 Propriétés
Comme l’opérateur divergence et l’opérateur gradient, l’opérateur rotationnel est linéaire
et vérifie les mêmes propriétés concernant la dérivation et les combinaisons linéaires.
Le théorème de Stokes justifie la relation suivante appelé Théorème du rotationnel :
Pour toute surface S, délimitée par le contour fermé C , pour tout champ vectoriel M~ , on
a:

→ → −
I Ï
− −→ −→ → −
M ·d l = rot M · d S
C S
Autrement dit, le flux du rotationnel d’un champ vectoriel à travers une surface égale la
circulation de ce champ sur le contour de cette surface.

5.5 Laplacien
5.5.1 Définition
Soit M un champ scalaire. Alors l’opérateur laplacien est l’application qui à M associe la
divergence du gradient de M . Le laplacien est noté ∆ .
Formellement :
³−−−→ ´
∆M = div grad M

Expression explicite
L’expression complète du laplacien dépend du système de coordonnées choisies. Pre-
nons l’exemple utile des coordonnées cartésiennes dans l’espace de dimension 3 :

∂2 ∂2 ∂2
∆= + +
∂x 2 ∂y 2 ∂z 2

©Dr OUEDRAOGO 58 [email protected]

Ecole Sup. Polytech. de Kaya 2020-2021 GC1 et GHR1

5.6 Vecteur formel « nabla »

Il est possible de regrouper tous les outils vectoriels abordés jusqu’ici sous un même
cadre formel, qui facilite bien souvent leur mémorisation et leur utilisation. On introduit
ainsi le vecteur formel nabla.
Le vecteur formel nabla est l’opérateur :

∂
 
∂x
→
− ∂
∇ =  ∂y 
∂
∂z
On retrouvera les expressions des autres outils vectoriels à partir de lui, mais il est im-
portant de garder à l’esprit qu’en toute rigueur, d’une part, son utilisation est un abus de
notation et d’autre part elle n’a de sens qu’en coordonnées cartésiennes et en dimension 3.
Néanmoins, on l’utilise pour des raisons de simplicité bien au-delà de ce seul cadre, mais il
faudra se rappeler qu’il ne s’agit alors que d’un moyen mnémotechnique.

5.6.1 Gradient, divergence, rotationnel

~ un champ vectoriel, soit N un champ scalaire.
Soit M

−−−→ →
−
grad N = ∇ N

−
→ → − − →
div M = ∇ · M

−→ −→ → − − →
rot M = ∇ ∧ M

En dimension 3, où ces notations ont un sens, on retrouve quelques propriétés, concer-

nant par exemple le fait qu’un produit vectoriel est orthogonal à ses deux arguments.

5.6.2 Laplacien
Le laplacien est défini comme la divergence du gradient. Soit N un champ scalaire :
³−−−→ ´ → − ³→− ´ → −2
∆N = div grad N = ∇ · ∇ N = ∇ N

©Dr OUEDRAOGO 59 [email protected]