Mécanique des Fluides et Pompes Rotodynamiques
Mécanique des Fluides et Pompes Rotodynamiques
Yann MARCHESSE
ECAM LaSalle
40 Montée Saint-Barthélemy
69321 Lyon Cedex 05
[Link]
Yann Marchesse
Pôle Énergétique / LabECAM
ECAM LASALLE
40 Montée St.-Barthélemy
69321 Lyon Cedex 05
1 Définitions et rappels 3
1.1 Généralités . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
1.2 Conservation de la masse . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
1.3 Bilan de quantité de mouvement . . . . . . . . . . . . . . . . . . . . . . . 12
1.4 Théorème de Bernoulli . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
1.5 Pertes de charge . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
1.6 Présence d’une machine dans le réseau hydraulique . . . . . . . . . . . . 28
1.7 Mesure du débit dans les conduites . . . . . . . . . . . . . . . . . . . . . 29
2 Réseaux hydrauliques 35
2.1 Choix des conduites . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
2.2 Constitution d’un réseau hydraulique . . . . . . . . . . . . . . . . . . . . 37
2.3 Courbe caractéristique d’un circuit hydraulique . . . . . . . . . . . . . . 37
3 Les vannes 47
3.1 Exemples de vannes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
3.2 Caractéristiques hydrauliques des vannes . . . . . . . . . . . . . . . . . . 48
3.3 Course de la vanne . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
3.4 Modification de l’écoulement au passage d’une vanne . . . . . . . . . . 53
3.5 Mesure des coefficients des vannes . . . . . . . . . . . . . . . . . . . . . . 54
3.6 Application au cas de l’écoulement de la section 3.4 . . . . . . . . . . . . 56
3.7 Recommandations pour le dimensionnement des vannes . . . . . . . . . 58
3.8 Application . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59
Le but de ce chapitre est d’une part de rappeler les deux expressions mathémati-
ques obtenues par application de la conservation de la masse et du bilan de quantité
de mouvement, deux principes fondamentaux de la mécanique des fluides. Ces ex-
pressions seront écrites sous les formes intégrales et locales. On rappellera ensuite le
théorème de Bernoulli utile à toute étude d’écoulements dans les canalisations. En-
fin, le phénomène de pertes de charge sera abordé. Toutes ces notions seront utilisées
par la suite dans le document. On attachera une importance non négligeable aux
équations de conservations de masse et de quantité de mouvement, celles-ci après
quelques hypothèses avérées se déclinant en de multiples expressions non moins
utiles. Notons qu’écrire ces grands principes est une première étape, les comprendre
en est une autre. En effet, la théorie insuffisante pour construire des machines permet
au moins de s’en servir intelligemment et les transformer si besoin est.
1.1 Généralités
Cette partie présente brièvement quelques aspects liés à la physique des écoule-
ments. Ceci nous permettra essentiellement de définir le concept de particule fluide
basé sur l’hypothèse d’un milieu continu. Celle-ci est très utile car le mouvemement
d’un fluide peut alors être décrit sans faire l’analyse détaillée de la dynamique des
particules discrètes qui le composent.
4 Définitions et rappels
Hypothèse de continuité : les propriétés d’un milieu, qu’il soit solide ou fluide, ne
sont pas uniformément distribuées. Cette distribution apparaı̂t d’autant moins uni-
forme que l’échelle d’observation est petite. Dans le cadre de notre étude, on se place
à une échelle telle qu’un élément de volume appelé particule fluide, suffisamment pe-
tit pour que la mesure puisse être considérée comme locale, soit suffisamment grand
pour contenir un grand nombre de molécules. La matière apparaı̂t alors comme un
milieu continu. Les quantités associées à la matière, telles que la vitesse, la pression
et la température sont considérées comme réparties sur tout le domaine d’étude. On
les représente par des fonctions continues. L’hypothèse d’un milieu continu ne per-
met donc pas l’étude des écoulements à basse densité, des écoulements à travers les
milieux poreux et des phénomènes d’ondes de choc (ceux-ci se caractérisant par des
discontinuités des propriétés de l’écoulement telles que la température, la pression,...).
Homogénéité et isotropie : une spécificité du fluide est d’être un état condensé mais
désordonné de la matière. Il est alors naturel de supposer ses propriétés physiques
comme étant identiques en tout point du domaine d’étude (homogénéité) et demeu-
rant indépendantes d’un changement de direction (isotropie).
Les particules fluides d’un écoulement ne sont pas isolées, mais interagissent entre
elles en échangeant de la quantité de mouvement. Ce paragraphe se propose de ca-
ractériser les contraintes mises en jeu. Considérons pour cela une surface élémentaire
−
→
dS située dans un écoulement. Les efforts dF exercés sur cette surface par les parti-
cules fluides voisines permettent de définir un vecteur unitaire d’effort (ou densité
−→ −
→ − →
surfacique d’effort) T tel que dF = T dS appelé vecteur contrainte (Figure 1.2), et lié
au tenseur des contraintes σ, tenseur d’ordre deux, tel qu’en tout point M on ait :
−
→
T = σ−
→
n ou Ti = σij n j (1.1)
6 Définitions et rappels
Le tenseur des contraintes rassemble tous les efforts s’exerçant sur les surfaces des
particules fluides 1 . Ces contraintes sont caractérisées par des efforts normaux liés à
la pression et tangents à la surface appelées contraintes de cisaillement.
Figure 1.3 : Volume de fluide parallelépipédique initialement au repos (gauche) et comprimé sous l’effet de
la pression (droite).
1. À ces efforts, il faut rajouter ensuite les efforts s’appliquant sur le volume de la particule fluide
(les efforts de pesanteur par exemple).
Généralités 7
(Figure 1.4.b). Cette déformation est caractérisée par un taux de déformation exprimé
de façon généralisé en fonction du champ de vitesse :
!
1 ∂ui ∂u j
dij = + . (1.2)
2 ∂x j ∂xi
Figure 1.4 : Volume de fluide parallelépipédique initialement au repos (gauche) et cisaillé (droite).
Dans beaucoup de problèmes où les contraintes visqueuses doivent être prises en
compte, il est plus commode d’utiliser la viscosité cinématique ν = µ/ρ (m2 /s).
Pour les gaz, la loi de Sutherland permet d’estimer la viscosité à partir de sa
température :
3/2
µ( T ) T T0 + S
=
µ0 T0 T+S
où µ0 est la viscosité à la température T0 et S la température de Sutherland égale à
110,4 K pour l’air.
Nous traiterons ici que de fluide newtonien pour lequel le tenseur des contraintes
visqueuses s’écrit :
!
∂ui ∂u j 2 ∂u
τij = µ + − µ k δij (1.4)
∂x j ∂xi 3 ∂xk
où le deuxième terme devient nul dans le cas d’écoulement incompressible, caractérisé
−
→
par la relation div V = 0. Le terme −2/3 provient de l’hypothèse de Stokes 3 . Si on
considère ce type d’écoulement, la matrice du tenseur des contraintes visqueuses
s’écrit alors dans la base de vecteurs (−
→
x ,−
→
y ,−
→z):
∂u ∂v ∂u ∂w ∂u
2 + +
∂x ∂x ∂y ∂x ∂z
∂u ∂v ∂v ∂w ∂v
∂y + ∂x
τ = µ 2 +
∂y ∂y ∂z
∂u ∂w ∂v ∂w ∂w
+ + 2
∂z ∂x ∂z ∂y ∂z
Notons en premier lieu qu’un fluide au repos ne présente aucune contrainte visqueuse
(i.e. τ = 0), son champ de vitesse étant nul. D’autre part, les contraintes normales et
de cisaillement dépendent du gradient de vitesse de l’écoulement, et non pas de la
vitesse seule. Ainsi, des couches de particules fluides ayant une vitesse relative non
nulle présenteront des échanges de quantité de mouvement entre elles (i.e. les plus
rapides auront tendance à accélérer les plus lentes, et inversement les plus lentes
auront tendance à freiner les plus rapides). Nous allons l’illustrer sur un exemple
donné sur la figure 1.7. Considérons un écoulement caractérisé par le champ de vitesse
suivant :
−
→
V = u(y)− →x
On s’intéresse aux efforts visqueux exercés par les particules situées au-dessus d’un
élément de surface dS de vecteur normal sortant unitaire −→n (= − →
y d’après la figure).
3. Afin de satisfaire
certaines conditions mathématiques sur ce tenseur, celui-ci s’écrit τij =
µ ∂ui /∂x j + ∂u j /∂xi + λ∂uk /∂xk δij . L’hypothèse de Stokes entraı̂ne la relation 3λ + 2µ = 0. Insis-
tons sur le fait qu’il s’agit d’une pure hypothèse qui est conservée car confirmée par de nombreuses
expériences.
10 Définitions et rappels
Dans la plupart des écoulements visqueux, les contraintes normales (i.e. τ11 ,
τ22 , et τ33 ) sont faibles devant les contraintes de cisaillement et sont la plupart
du temps négligées. Les contraintes normales peuvent devenir très importantes
lorsque le gradient longitudinal de la vitesse est grand, ceci est le cas dans les
ondes de choc.
Tenseur des contraintes : comme nous l’avons dit plus haut, on note des efforts de
pression, ainsi que des efforts liés à la viscosité du fluide. On écrira alors le tenseur
des contraintes :
!
∂ui ∂u j 2 ∂u
σij = − pδij + µ + − µ k δij (1.6)
∂x j ∂xi 3 ∂xk
Rappelons que le fluide newtonien est caractérisé par un taux de déformation linéaire
avec la contrainte de cisaillement (cf. relation 1.3).
La suite du document présente l’application de deux principes fondamentaux sur
des volumes de contrôle ou sur des volumes élémentaires, conduisant à des expres-
sions intégrales ou locales. Un volume de contrôle est de taille assez importante pour
Conservation de la masse 11
contenir beaucoup de fluide, alors que le volume élémentaire est par définition plus
petit, et n’autorise donc pas de bilan intégral. Il ne s’agira néanmoins que de plusieurs
écritures de la “même” équation.
Figure 1.8 : Représentation d’un volume de contrôle. Dans le cas (a), le volume de contrôle est fixe et est
traversé par le fluide en mouvement ; dans le cas (b), le volume de contrôle est en mouvement,
et est alors constitué des mêmes particules fluides.
−
→→
Z Z
ρ V .−
∂
ρdv + n dS = 0 (1.7)
∂t
D S
D ∂ → −−→
−
= + V .grad (1.8)
Dt ∂t
La conservation de la masse s’écrit alors :
D
Z
ρdv = 0 (1.9)
Dt
Dt
4. Bien que Newton ait postulé, en 1687 dans son Principia, le fait que la masse spécifique d’un
objet soit constante, ce principe ne fut appliqué à la mécanique des fluides qu’en 1749 par d’Alembert
sous forme dérivée dans un article intitulé Essai d’une nouvelle théorie de la résistance des fluides. Huit
années plus tard, Euler proposa les équations de continuité et de bilan de quantité de mouvement.
Deux des trois équations de bilan utilisées de nos jours ont donc vu le jour avant la fin du 18ème
siècle. La troisième équation de bilan, énergétique elle, ne viendra qu’avec le développement de la
thermodynamique au 19ème siècle.
12 Définitions et rappels
avec l’indice t pour rappeler le fait que le volume de contrôle est en mouvement, et
peut donc voir sa forme évoluer au cours du temps.
Les relations mathématiques (1.7) et (1.9) sont des écritures intégrales du prin-
cipe de conservation de la masse. Des écritures locales peuvent être atteintes en
considérant un volume élémentaire au lieu d’un volume de contrôle.
Considérons maintenant un volume élémentaire dV de dimension dx × dy × dz
(Fig. 1.9) et fixe dans l’espace. Le bilan du débit massique sur toutes ses surfaces
permet d’écrire une première forme locale de la conservation de la masse :
∂ρ ∂
+ (ρu j ) = 0 (1.10)
∂t ∂x j
Dρ −
→ ∂ρ − → −−→ −
→ ∂ρ −
→
+ ρ × div V = + V .gradρ + ρ × div V = + div ρ V
∂t ∂t ∂t
correspondant ainsi au membre de gauche de la relation (1.10).
(a) (b)
Figure 1.11 : Application du bilan de quantité de mouvement pour le calcul de la charge délivrée par une
pompe centrifuge à un écoulement. a, vue de deux aubes successives d’une pompe centrifuge ;
b, domaine de contrôle associé à ces deux aubes.
Encore une fois, il est important de comprendre que les expressions (1.12), (1.14), et
(1.17) sont des représentations différentes d’une unique équation caractérisant le bilan
de quantité de mouvement.
Ces équations expriment le fait que les variations au cours du temps de la vitesse
peuvent être créées par l’application de forces s’exerçant sur le fluide (la gravité ou la
différence de pression par exemple) mais aussi que la réponse du fluide va dépendre
de deux phénomènes de nature très différentes : (1) la friction entre des éléments de
fluide voisins se déplaçant à des vitesses différentes et (2) l’effet d’entraı̂nement du
fluide par son propre mouvement (effet d’inertie). Cet effet d’inertie aura tendance
à dominer l’évolution du fluide lorsqu’il se déplace à grande vitesse. L’importance
relative de ces deux contributions détermine le régime de l’écoulement (laminaire ou
turbulent).
Les équations de Navier-Stokes sont très compliquées et difficiles à résoudre,
ceci pour plusieurs raisons : (i) il s’agit d’une équation vectorielle dans laquelle les
équations scalaires associées (une pour chaque direction) sont couplées de façon com-
plexe entre elles ; (ii) le terme de convection Vj ∂Vi /∂x j entraine la non linéarité des
équations ; (iii) le couplage pression-vitesse renforce cette difficulté (lorsque le terme
de pression devient négligeable, la solution est atteinte).
1
p + ρV 2 + ρgz = Cste (1.18)
2
5. Daniel Bernoulli, physicien suisse (1700-1782). Son importance en mécanique des fluides est
centrée sur son livre Hydrodynamica (1738) dans lequel il introduit le terme hydrodynamique. C’est
dans ce livre qu’il essaie de définir les relations entre la pression et la vitesse de l’écoulement. Mal-
heureusement, son aboutissement reste obscure, et sa formulation est très éloignée de celle que nous
connaissons aujourd’hui. Il faut attendre la contribution d’Euler pour atteindre la relation différentielle
(1.18) utilisée de nos jours : dp = −ρUdU.
16 Définitions et rappels
En hydraulique, la pratique veut que cette équation soit exprimée en énergie par
unité de poids, obtenue en divisant l’expression (1.18) par ρg :
p V2
+ + z = Cste (1.19)
ρg 2g
Ces termes
sont homogènes à des longueurs et définissent des hauteurs, dynamique
2 ∗
V /2g et piézométriques (z + p/ρg = p /ρg). L’ensemble est la charge hydraulique
H du fluide, elle représente en hauteur de fluide, l’énergie mécanique totale par unité
de poids du fluide. On peut écrire en chaque point i d’un réseau hydraulique :
Vi2 p
Hi = zi + + i
2g ρg
Hi représente l’altitude qu’atteindrait le fluide s’il était ralenti sans frottement jusqu’à
une vitesse nulle. De même s’il n’existe pas de perte dans l’écoulement alors la conser-
vation de l’énergie s’écrit H = Cste. Si ce n’est pas vérifé expérimentalement, il existe
alors des pertes de charge dans le réseau (cf. paragraphe suivant). Le tracé de Hi en
fonction de positions particulières dans le circuit représente la ligne de charge. Dans
un écoulement, cette dernière a une pente négative consécutivement à des pertes dans
le circuit sauf si on interpose un organe capable de fournir de l’énergie au fluide : une
pompe.
1 1
p1 + ρgz1 + ρV12 = p2 + ρgz2 + ρV22 + ρg∆H ou H1 = H2 + ∆H
2 2
La quantité ∆H est toujours positive et définit en hauteur d’eau la perte d’énergie
mécanique, appelée perte de charge, subie par l’écoulement entre les deux sections.
L V2 L 8Q2
∆Hl = λ =λ 5 2 (1.20)
D 2g D π g
Lorsque la conduite a une section de forme non-circulaire, la section peut alors être
ramenée à un diamètre équivalent, appelé diamètre hydraulique. Chézy a établi une
relation théorique pour des gaines de forme quelconque Dh = 4 × Surface/Périmètre.
La valeur du coefficient de pertes de charge dépend du régime de l’écoulement, ca-
ractérisé par le nombre de Reynolds qui est le rapport des efforts d’inertie et de
viscosité :
Force d’inertie ρV 2 /D ρVD
Re = = = (1.21)
Force visqueuse µV/D2 µ
Lorsque ce nombre croı̂t, les différents régimes suivants s’établissent successivement :
(1) un régime laminaire où λ est inversement proportionnel au nombre de Reynolds,
(2) un régime turbulent où λ est d’abord uniquement fonction du nombre de Rey-
nolds, on parlera de conduite lisse, (3) une deuxième phase du régime turbulent dite
phase de transition, où λ est à la fois fonction de Re et de la rugosité relative ϵ/D,
18 Définitions et rappels
avec ϵ caractérisant la “taille” de la rugosité des parois, et (4) une phase où λ est
indépendant de Re et uniquement fonction de ϵ/D. Dans cette région, les pertes de
charge sont strictement proportionnelles au carré de la vitesse (on parle de régime de
turbulence pleinement développée). La figure 1.12 représente de façon schématique
l’évolution du coefficient λ avec le nombre de Reynolds. Les trois régimes sont iden-
tifiés par des nombres de Reynolds distincts, Rec , Re1 et Re2 . Alors que la valeur
numérique du nombre de Reynolds critique (Rec ) est toujours la même, celles de Re1
et Re2 est fonction de la rugosité de la paroi. La succession de ces différentes phases
peut être expliquée par la présence d’une couche de fluide proche de la paroi dans
laquelle l’écoulement voit sa vitesse passer brutalement d’une vitesse nulle en paroi
(adhérence du fluide) à la vitesse de l’écoulement. Il existe donc dans cette couche
limite, une sous-couche dans laquelle le mouvement est laminaire (sous-couche vis-
queuse). Quand le nombre de Reynolds est faible, la région laminaire occupe toute la
section de la conduite. Quand le nombre de Reynolds augmente, l’écoulement central
turbulent apparaı̂t, mais la couche périphérique reste laminaire et est encore assez
épaisse pour noyer toutes les irrégularités de la surface ; l’écoulement s’effectue donc
en tuyau lisse. Puis le nombre de Reynolds croissant, l’épaisseur de la sous-couche
laminaire diminue, et les irrégularités de surface émergent. Quand elles sont en ma-
jorité émergées, elles créent une turbulence telle qu’elles masquent complètement les
effets de la viscosité et le phénomène devient indépendant du nombre de Reynolds.
Figure 1.12 : Évolution du coefficient de pertes de charge régulières, λ, en fonction du nombre de Reynolds.
de grains de sable tamisés. Les résultats de ses expériences mettent en évidence les
différents régimes d’écoulements cités précédemment via des évolutions différentes
du coefficient de pertes de charge (Fig. 1.13). Cependant cet état de rugosité de sur-
face régulier n’est pas observé dans le cas des conduites industrielles pour lesquelles
des hétérogénéités apparaissent. Les évolutions du coefficient λ avec le nombre de
Reynolds, ceci pour des rugosités différentes, sont alors quelque peu différentes de
celles notées par Nikuradse. L’évolution du coefficient de pertes de charge pour les
conduites industrielles est donnée sur le diagramme de Moody, figure 1.14. L’estima-
tion du coefficient de pertes de charge régulières est précise à ± 15% sur la gamme
du nombre de Reynolds tracée sur cette figure.
Il est possible d’estimer la valeur numérique du coefficient de perte de charge λ à
partir de modèles. La formulation à utiliser dépend du régime d’écoulement traduit
par la valeur numérique du nombre de Reynolds :
⋆ Écoulement laminaire : lorsque Re < 2000, l’écoulement est calme et régulier,
et les petites perturbations qui pourraient être créées sont amorties et ne peuvent
pas se développer plus. Le régime d’écoulement est dı̂t laminaire et le coeffi-
cient de perte de charge régulière est évalué selon la formule de Hagen 6 [49] -
Poiseuille 7 [69] :
64
λ= (1.22)
Re
Le régime laminaire est rare dans l’industrie (écoulement de pétrôle dans les
pipe-lines, canalisation d’huile de graissage−coussinet, etc.). Cependant, cette
configuration peut être utile dans l’évaluation expérimentale de la viscosité dy-
namique d’un fluide puisqu’il est possible d’écrire après transformation, dans
ce cas laminaire, ∆p = 32µVL/D2 . Ainsi, la mesure de la vitesse et des pertes
de charge permet de remonter à la valeur de la viscosité dynamique, µ.
⋆ Écoulement transitoire : lorsque 2000 < Re < 2500, le régime est transitoire au
sens où il est situé entre le régime laminaire et celui turbulent présenté dans la
suite. La valeur numérique du nombre de Reynolds critique, en-dessous duquel
le régime laminaire est stable, dépend généralement des conditions à l’entrée
[66]. Il n’y a pas de valeur numérique suffisamment précise dans cette région,
et il n’y a donc pas de formulation proposée.
⋆ Écoulement turbulent : lorsque Re > 2500, l’écoulement est agité et ne s’effec-
tue plus par des filets rectilignes. Le régime d’écoulement est dı̂t turbulent et
le comportement est différent suivant que la conduite est lisse ou rugueuse (du
point de vue du fluide).
▶ Écoulement turbulent lisse : les aspérités seront, dans ce cas, complètement
noyées dans la sous-couche laminaire et la rugosité ne joue aucun rôle dans
les pertes de charge. Cependant le profil de vitesse longitudinale évolue
avec le nombre de Reynolds et il devient alors compliqué de proposer
un critère stipulant d’impliquer ou non la rugosité. Parmi les nombreuses
6. Gotthilf Hagen (né le 3 mars 1797 à Königsberg et mort le 3 février 1884 à Berlin) est un ingénieur
hydraulicien prussien spécialiste des voies navigables
7. Jean-Léonard-Marie Poiseuille (né le 22 avril 1797 à Paris et mort le 26 décembre 1869 à Paris) est
un physicien et médecin français
20 Définitions et rappels
Figure 1.13 : Évaluation du coefficient de pertes de charge linéaires en fonction du nombre de Reynolds
et de la rugosité relative de la conduite d’après l’expérience de Nikuradse [64] (Dans cette
expérience la rugosité des parois est maitrisée par l’utilisation de grains calibrés et elle ne
correspond donc pas strictement aux écoulements dans les conduites industrielles présentant
quant à elles des tailles de rugosités diverses.).
Figure 1.14 : Diagramme de Moody pour des conduites industrielles (évaluation du coefficient de pertes
de charge linéaires en fonction du nombre de Reynolds et de la conduite) pour lesquelles les
rugosités sont hétérogènes. On remarque ici l’absence du décrochement en V noté dans le
cas des conduites ayant des rugosités homogènes pour des nombres de Reynolds proches de
2500 (i.e. log(Re) = 3,4) − d’après White [92].
Pertes de charge 21
études, nous prendrons ici le critère proposé par la référence [66] stipulant
que la rugosité n’a aucun rôle dans les pertes de charge régulières tant que
Re7/8 × ϵ/D < 19, 25. La valeur numérique du nombre de Reynolds Re1
peut alors être évaluée à partir de cette équation (i.e. Re7/8
1 × ϵ/D = 19, 25).
Plusieurs formulations sont dès lors proposées en fonction de la valeur
numérique du nombre de Reynolds. Lorsque Re< 105 , une première for-
mulation possible est proposée par Blasius 8 [20] :
0, 316
λ= (1.23)
Re0,25
Cette formulation est implicite, et donc le coefficient λ y est estimé par voie
itérative.
▶ Écoulement turbulent rugueux − dépendant du nombre de Reynolds : la plus
grande partie des aspérités est située cette fois-ci au-delà de la sous-couche
laminaire et la rugosité joue donc un rôle dominant dans les pertes de
charge. Ceci est observé lorsque Re1 < Re < Re2 = 24 (Re1 − 1, 528 × D/ϵ)
[66]. La formule de Colebrooke 11 , faisant intervenir la rugosité relative ϵ/D,
peut être utilisée :
1 ϵ/D 2, 51
√ = −2 × log +√ (1.25)
λ 3, 7 λ × Re
Swamee et Jain [81] ont proposé une version simplifiée explicite de cette
équation :
2
ϵ/D 5, 74
λ = 0, 25/ log + 0,9 (1.26)
3, 7 Re
▶ Écoulement turbulent rugueux − indépendant du nombre de Reynolds : ce type
d’écoulement est rencontré dès lors que Re > Re2 . Le coefficient de perte
de charge régulière ne dépend plus du nombre de Reynolds et peut alors
être évalué par la formulation de Nikuradse :
1 ϵ/D
√ = −2 × log (1.27)
λ 3, 7
8. Heinrich Blasius (né le 9 août 1883 à Berlin et mort le 24 avril 1970 à Hambourg) est un ingénieur
hydraulicien et professeur allemand.
9. Theodore von Karman (né le 11 mai 1881 et décédé le 6 mai 1963) est un ingénieur et physicien
hongrois et américain.
10. Ivan (naturalisé en Johann) Nikuradze (né le 20 novembre 1894 à Samtredia et mort le 18 juillet
1979) est un ingénieur hydraulicien allemand d’origine géorgienne.
11. Cyril Frank Colebrooke (né le 26 juillet 1910 à Swansea (Pays de Galles) et mort le 12 janvier 1997
à Worthing, en Angleterre.) est un physicien britannique.
22 Définitions et rappels
Dénomination ϵ (mm)
Conduite en acier neuve 0,15
Conduite en acier usagée 0,30
Conduite en fonte neuve 0,25
Conduite en fonte usagée 0,50
Conduite en béton lisse 0,15 à 0,20
Application
Q 4Q 4 × 1, 55 × 10−3
V= = = 1,52 m/s
S πD2 π × 0, 0362
L’écoulement est dans le régime turbulent, et reste à savoir si la rugosité est à considérer
dans notre calcul. Il faut pour cela calculer la valeur numérique de la quantité Re7/8 ×
ϵ/D :
(7/8)
Re7/8 × ϵ/D = 5, 4 × 104 × 0, 15/36 = 57, 6
Ainsi Re7/8 × ϵ/D > 19, 25 et la rugosité influence les pertes de charge régulières. La
formulation de Colebrooke peut être utilisée et on obtient pas voix itérative λ = 0,036.
Le diagramme de Moody propose, quant à lui, la valeur √ approchée λ = 0,030 et la
formulation pour les conduites industrielles λ = 0, 79 ϵ/D = 0,051. Nous garderons
λ = 0,036 dans nos calculs. Les pertes de charge linéaires, ∆Hl , valent donc :
8L 8 × 150
−3
2
∆Hl = λ 2 5 Q = 0, 036 2
2
1, 55 × 10 = 17, 73 m.
π gD π × 9, 81 × 0, 0365
La perte de charge vaut, si elle est exprimée cette fois-ci en terme de pression :
de cette phase. D’autre part, le point de transition est difficilement déterminé par
l’observation. On préfèrera alors effectuer des mesures et étudier l’évolution d’un pa-
ramètre qui dépend du nombre de Reynolds. Dans le cas des conduites, on s’intéresse
à l’évolution du coefficient de pertes de charge en fonction du nombre de Reynolds
(Fig. 1.17). Celui-ci voit ses valeurs décroitrent tant que le nombre de Reynolds reste
inférieur à une première valeur critique (notée Rec1 ). Le coefficient de perte de charge
augmente brutalement au passage de la transition. Pour des nombres de Reynolds
supérieurs à une deuxième valeur critique (notée Rec2 ), le comportement turbulent se
traduit par une nouvelle loi de décroissance différente de celle observée en laminaire.
On pourra alors trouver un critère de transition, Rec tel que
1
Rec = [Rec1 + Rec2 ] (1.29)
2
V2
∆Hs = ξ (1.30)
2g
ξ dépend de la forme géométrie de la singularité, et est déterminée la plupart du
temps par voie expérimentale. V est la vitesse de l’écoulement et dans le cas d’un
changement de section (i.e. élargissement ou rétrecissement brusque) la vitesse la plus
importante des deux sera considérée dans la formulation. L’ouvrage “Memento des
pertes de charge − Idel’Cik” [82] propose des formulations d’évaluation des coefficients
de pertes de charge dans diverses configurations simples. Voici quelques singularités
régulièrement observées dans une installation industrielles.
Élargissement brusque
Dans le cas d’un élargissement brusque (Figure 1.18.a), le coefficient de pertes de
charge singulières peut être évalué à l’aide de la formule de Borda-Carnot :
2
S
ξ= 1− 1 (1.31)
S2
Rétrécissement brusque
L’écoulement subit cette fois-ci une décroissance brutale de sa section de pas-
sage (Figure 1.18.b). S1 représente la section de passage la plus faible comme c’était
26 Définitions et rappels
(a) (b)
Figure 1.19 : coefficient de perte de charge d’un rétrécissement brusque lorsque Re < 104 et pour différents
rapports S1 /S2 (d’après Idel’Cik [82]).
(a) (b)
Figure 1.20 : Débouché d’une conduite dans un réservoir (a) et départ d’une conduite d’un réservoir (b).
Coude
Dans le cas de conduites circulaires coudées, ayant des diamètres identiques en
entrée et en sortie de coude, pour lesquelles Rc /D > 1, Rc étant le rayon de cour-
bure, le coefficient de pertes de charge peut être évalué à l’aide de la formulation de
Weisbach : " #
D 7/2 θ
ξ = 0, 13 + 1, 85 (1.33)
2Rc 90
avec l’angle θ l’angle du coude exprimé en degré. Lorsque le rayon de courbure est
nul (i.e. coude brusque) et θ = 90 degré, ξ = 1,2.
Les coefficients de pertes de charge sont donnés par des modèles semi-empiriques.
Une fois les pertes de charges totales du réseau estimées, on prend l’habitude de
majorer de 10% la valeur estimée.
V2 Le V 2 D
ξ =λ conduit à la relation Le = ξ (1.34)
2g D 2g λ
Les valeurs de la conduite virtuelle, D et λ, sont prises égales aux valeurs de la
conduite composant le circuit hydraulique. Tout se passe donc comme si la présence
de singularités augmentait la longueur de conduite du réseau. Cette notion est sou-
vent utilisée dans les catalogues de fournisseurs de matériels hydrauliques.
p1 V2 p V2 Pv Pu
+ z1 + 1 = 2 + z2 + 2 + −
ρg 2g ρg 2g ρgQ ρgQ
Mesure du débit dans les conduites 29
avec Pv la puissance dissipée par effet visqueux et donc associée aux pertes de charge.
Dans le cas des turbines, des éoliennes ou des moulins (Fig. 1.21.a), les aubes
prélèvent de l’énergie mécanique à l’écoulement, i.e., Pu < 0. L’énergie est alors ini-
tialement sous forme d’énergie potentielle (cas des retenues d’eau des stations hy-
droélectriques) ou cinétique (cas des éoliennes ou des moulins). L’énergie mécanique
est ensuite convertie en puissance électrique ou alors en une autre forme d’énergie
mécanique, comme cela était le cas pour les moulins broyant les céréales. Dans le cas
des pompes (Fig. 1.21.b), la roue, clavetée sur un arbre mis en rotation par un mo-
teur électrique, dispose d’aubes. Ces dernières fournissent l’énergie à l’écoulement, la
charge en sortie de la machine se trouvant dès lors plus importante qu’à son entrée.
Le chapitre 4 est consacré à l’utilisation des pompes centrifuges dans les réseaux hy-
drauliques.
(a) (b)
Figure 1.21 : a, Roue à eau conçue par Jacob Leupold. L’originalité de cette roue repose sur la présence
de tuyaux tuyauxcourbés alimentant en eau la région centrale de la roue ; b, Installation
hydraulique avec une pompe de type centrifuge et montée en charge (d’après [2]).
Le Venturi
Le tube de Venturi est composé de deux parties consécutives, un convergent suivi
d’un divergent (Fig. 1.23.a). La fabrication de ce débitmètre nécessite une grande
Mesure du débit dans les conduites 31
précision et de faibles tolérances expliquant ainsi son coût élevé. Il est principale-
ment placé dans des canalisations sous pression. Le Venturi classique détaillé dans la
norme ISO 5167-1 possède un angle au sommet de la partie convergente valant 21 ˚ ,
et celui de la partie divergente valant 7 ˚ . Les pertes de charge par un tel débitmètre
restent faibles.
(a) (b)
Figure 1.23 : Représentation d’un Venturi (a) et d’un diaphragme (b) (d’après A. Dupont [34]).
Les mesures de pression sont réalisées sur deux sections placées légèrement en
amont et en aval du convergent. Les deux prises de pression sont reliées à un ma-
nomètre différentiel à mercure et la différence de hauteur ∆h enregistre la quantité
( p1 − p2 )/(ρg). Le débit massique est formulé de la façon suivante
q
Q = KV Sc 2ρ ( p1 − p2 )
Le diaphragme
est alors augmentée. Un chanfrein est alors réalisé avec un angle égal à 45° [6]. Étant
donnée sa simplicité son coût est très faible et l’installation ou le remplacement reste
aisé. Contrairement au Venturi pour lequel la variation de section de passage est
contrôlée, ici cette section de passage passe brutalement d’une section de diamètre D à
une section plus faible de diamètre d. De ce fait, les pertes de charge sont supérieures.
D’autre part, la place occupée par les diaphragmes devient plus faible.
Le principe de mesure étant identique à celui du Venturi, la relation débit massique-
hauteur est de la même forme
q
Q = K D Sc 2ρ ( p1 − p2 )
Le rotamètre
Le Rotamètre (de la société Rota) est un débitmètre constitué d’un tube vertical
(Fig. 1.24.a) dont la section de passage augmente du bas vers le haut et dans lequel un
organe, appelé de façon impropre flotteur, est animé d’un mouvement de translation
(et parfois de rotation). L’altitude du flotteur est fonction de son poids propre, de
l’effort d’Archimède et de la traı̂née générée par le passage de l’écoulement autour de
l’obstacle. Seul ce dernier effort dépend de la vitesse et donc du débit de l’écoulement.
Une relation entre ce débit et l’altitude observée est alors possible.
L’application du théorème de Bernoulli et de la conservation de la masse permet
d’écrire la relation entre le débit et la section de passage locale, S2 :
Q = KS2
La turbine
La turbine est composée d’une hélice libre de tourner autour d’un axe aligné avec
la conduite de l’écoulement (Fig. 1.24.b). La vitesse de rotation peut être mesurée par
un capteur magnétique.
Mesure du débit dans les conduites 33
(a) (b)
Figure 1.24 : Représentation d’un rotamètre (d’après A. Dupont [34]) (a) et d’une turbine (b) (d’après F.
Kreith [52]).
Un débitmètre à effet vortex est basé sur la génération d’instabilités par un cylindre
placé dans un écoulement. La fréquence de cette génération est caractérisée par le
nombre de Strouhal St = f D/V où f est la fréquence des lâchers de tourbillons, D
le diamètre du cylindre circulaire, et V la vitesse de l’écoulement à caractériser. Le
nombre de Strouhal dépend du nombre de Reynolds de l’écoulement mais demeure
constant sur la plage Re = [300 - 100 000]. L’estimation de la vitesse d’écoulement
repose donc sur la mesure de cette fréquence réalisée la plupart du temps à partir de
jauges de déformations (Fig. 1.25.a).
(a) (b)
Figure 1.25 : Représentation d’un débitmètre à effet vortex (a) et d’un débitmètre à effet Coriolis. Respec-
tivement d’après H. Yamazaki et al. [95] et d’après D. Liu [18].
avec c le coefficient du tube de Pitot égal à 1 pour des tubes standards. La valeur de
ce coefficient est différente de l’unité lorsque les effets de compressibilités deviennent
importants.
13. Henri de Pitot, physicien français (1695-1771), mis en œuvre une “machine pour mesurer la
vitesse des eaux courantes et le sillage des vaisseaux”.
Chapitre 2
Réseaux hydrauliques
Les réseaux hydrauliques ont pour rôle d’acheminer, dans des conditions mai-
trisées, un fluide d’un endroit à un autre pour y être utilisé, lors de son parcours ou
à l’extrémité du réseau. Les réseaux hydrauliques dans les bâtiments permettent par
exemple de distribuer la chaleur et le froid via des échangeurs thermiques. Le fluide
doit alors avoir suffisamment d’énergie pour circuler et lutter contre les dissipations
générées par les composants du réseau. Le but de ce chapitre est de présenter la
démarche permettant d’obtenir la courbe du système hydraulique à alimenter. Cette
courbe sera tracée sur le même graphique que celle caractérisant le système de pom-
page (présenté au chapitre suivant).
Tableau 2.1 : Vitesse en m/s pour diverses tuyauteries (d’après Oniga [66]).
Constitution d’un réseau hydraulique 37
débits identiques, les pertes de charges s’accumulent alors. La courbe résultante est
obtenue en considérant différents débits, et en additionnant les pertes de charge des
deux conduites (Fig. 2.3.b).
L’assemblage des courbes caractéristiques de chaque conduite nécessite une bonne
analyse du réseau − qu’il ne faut pas négliger − à partir du schéma.
Figure 2.3 : Courbes caractéristiques résultantes de conduites placées en parallèle (a) ou en série (b).
Lorsque les conduites sont disposées en parallèle, pour différentes valeurs de pertes de charge,
la courbe résultante est atteinte par addition des débits traversant chaque conduite ; lorsque
les conduites sont disposées en série, pour des différentes valeurs de débits, on additionne les
pertes de charge.
Courbe caractéristique d’un circuit hydraulique 41
Figure 2.4 : Prise en compte de la différence d’altitude entre le bassin d’alimentation et le réservoir. À
gauche, le réseau à alimenter ; à droite, la courbe caractéristique de la conduite débouchant
sur le réservoir.
Figure 2.5 : Prise en compte de la différence d’altitude entre plusieurs bassins d’alimentation dont les
surfaces libres sont toutes localisées à la même altitude et plusieurs réservoirs. À gauche, le
réseau à alimenter ; à droite, les courbes caractéristiques des conduites débouchant sur les
réservoirs.
42 Réseaux hydrauliques
Ceci est illustré sur la figure 2.6.a avec deux réservoirs à alimenter étant situés à des
altitudes différentes, la courbe résultante notée C1 ||C2 étant tracée en gras.
Figure 2.7 : Courbes caractéristiques d’une conduite AC composée d’une pompe et de sa conduite d’aspi-
ration (Les pertes de charge dans la conduite BC sont négligées). Les autres parties du réseau
ne sont pas tracées.
est tracée à partir de l’origine et ne fait intervenir que ses pertes de charge. Dans ce
cas précis, après avoir tracé la courbe résultante des deux conduites desservant les
réservoirs, notée C1 ||C2 , on ajoute les pertes de charge de la conduite C3 .
Figure 2.8 : Réseau présentant un débit de fuite imposé dans une conduite (La courbe “Réseau” correspond
aux conduites CAsp. et C1 ).
(a) (b)
Figure 2.9 : Collecteur en appui d’une pompe (a) et collecteur à remplir (b).
Les vannes 1 sont des éléments importants d’un réseau hydraulique et sont utilisées
pour de nombreuses actions. Nous nous intéresserons ici qu’aux vannes dont le rôle
est de réguler le débit qui diffèrent des vannes de contrôle par le fait qu’elles ne sont pas
actionnées de façon pneumatique ou hydraulique permettant alors un contrôle avec
précision de l’écoulement. Ainsi le terme vanne de contrôle est plus associé à la fonction
de la vanne qu’à son type. Si les vannes ne sont pas correctement dimensionnées ou
manipulées, elles peuvent être à l’origine d’incidents graves. Fermer ou ouvrir, par
exemple, une vanne trop rapidement entraine l’apparition de phénomènes transitoires
destructeurs, ce qui sera évoqué dans un prochain chapitre. Si elles sont soumises à
de la cavitation, elles peuvent s’endommager grandement et être à l’origine de fuites,
ce qui nécessitera alors leur remplacement.
3.1 Exemples de vannes
Les vannes de régulation sont généralement utilisées dans des boucles de contrôle
dont le but est de maintenir une variable du procédé (débit, pression, température,
etc.) dans une plage de valeurs souhaitées. L’organe de la vanne assurant la modifica-
tion du débit est le clapet obturant plus ou moins le passage de l’écoulement par une
modification de sa position dans la partie interne de la vanne. Le corps de vanne est
le système mécanique assurant le réglage de la position du clapet (et donc du débit).
Le clapet est déplacé en translation (corps droit) ou en rotation (corps d’angle). Parmi
les nombreuses vannes de régulation, on note :
1. Le terme vanne est traduit en anglais par le terme valve mais a une signification beaucoup plus
large puisqu’elles concernent les vannes de régulation (throttling valves), les vannes d’arrêt (on-off valves)
et les clapets (nonreturn valves).
48 Les Vannes
Q
Kv = p
∆p/d
2. Le système métrique ne correspond pas forcément aux unites du système international.
Caractéristiques hydrauliques des vannes 49
(a) (b)
(c) (d)
Figure 3.1 : (a) Vue en coupe d’une vanne à tige montante dans laquelle la dissipation énergétique est
assurée par la perturbation de l’écoulement à travers la section de passage laissée par le clapet
(d’après Fisher [40]) ; (b) vanne papillon dans laquelle le clapet est un disque qui peut tourner
à 90° autour d’un axe. En fonction de l’angle que fait le clapet, la surface de passage laissée
à l’écoulement diffère (voir section 3.4) et modifie alors les dissipations énergétiques (d’après
[Link]) ; (c) Vanne à cône (le cône est localisée dans la partie droite et peut être
déplacé vers la gauche − cette vanne est utilisée pour libérer un écoulement vers l’atmosphère
− d’après [Link]) ; (d), vanne d’arrêt (le clapet peut être déplacé verticalement
obstruant ainsi de façon plus ou moins importante le passage de l’écoulement − d’après
[Link]).
50 Les Vannes
Ainsi une vanne n’aura pas les mêmes valeurs numériques des coefficients Cv et Kv
pour une même ouverture. Le coefficient ξ v est le coefficient de perte de charge sin-
gulière tel que nous l’avons déjà vu dans le chapitre de rappel.
Lorsque l’on compare des coefficients de vannes identiques mais de tailles différen-
tes, les coefficients ξ v et Cd sont identiques (puisqu’ils sont sans dimension) alors que
les valeurs de Cv dépendront de la taille. Néanmoins les rapports d’échelles de deux
vannes de diamètres D1 et D2 entraineront Cv1 = Cv2 ( D1 /D2 )2 .
Certaines vannes sont caractérisées par des pertes de charge nulles (ou quasiment
nulles) lorsqu’elles sont totalement ouvertes. Il est alors préférable d’utiliser le coeffi-
cient Cd qui varie entre 0 et 1 alors que les deux autres coefficients tendent vers l’infini
lorsque l’ouverture est totale.
Il est possible d’évaluer la valeur numérique d’un coefficient connaissant la valeur
numérique d’un autre à partir de ces formulations :
1 1
ξv = −1 et Cd = √ (3.4)
Cd2 1 + ξv
et
S2
Cv = 8, 146 × 104 √
1 + ξv
avec S l’aire de la section de passage exprimée en m2 . Les valeurs numériques de
ces coefficients sont tracées généralement en fonction du pourcentage d’ouverture
qu’offre le clapet à l’écoulement. Néanmoins, pour certaines vannes le fait d’être tota-
lement ouverte ne signifie pas que l’aire de passage correspond à celle de la conduite.
Une vanne papillon par exemple obture légèrement le passage même quand elle est
pleinement ouverte. La valeur 100% d’ouverture correspondra alors au pourcentage
vis-à-vis de l’aire maximale de passage. Des tracés représentatifs de coefficient de
pertes de charge singulières sont proposées sur la figure 3.2. On note des rapports de
coefficients très élevés entre une ouverture totale (100%) et une ouverture très faible
(10%).
(a)
(b)
Figure 3.2 : Évolution des coefficients Cd et de perte de charge singulière des vannes à soupape, papillon
et à cône (Données extraites de la référence [84]).
52 Les Vannes
Figure 3.3 : Forme de la cage d’un clapet de vanne en fonction du type d’évolution du coefficient Cv relatif
(d’après Emerson−Fisher [41]).
utilisé dans le cas où la vanne produit la majorité des pertes de charge du
réseau total.
— la course à ouverture rapide : le débit augmente plus rapidement à l’ouverture
de la vanne comparativement à la région de la fermeture. Ce comportement
est utile lorsque l’on souhaite un fonctionnement proche du tout-ou-rien (on-off
valve) puisque sur une large gamme d’ouverture la vitesse de l’écoulement est
proche de la vitesse maximale.
— la course à égal pourcentage : la variation de Cv par unité d’ouverture est direc-
tement proportionnelle à la valeur de Cv avant changement, i.e. ∆Cv /∆O = kCv .
On peut réécrire la formulation de la sorte : ∆Cv /Cv = k∆O. Ce type de vanne
permet un réglage fin sur la première moitié d’ouverture de la vanne puisqu’il
faut déplacer de façon significative le clapet pour faire varier le débit. D’autre
part l’augmentation du débit est élevée dans la seconde moitié permettant alors
au débit d’atteindre la valeur souhaitée. Parmi les trois types d’évolution, la
course à égal pourcentage est la plus utilisée dans les vannes de contrôle.
Ainsi le type de vanne doit être choisi scrupuleusement afin de satisfaire les usages
de celle-ci.
Modification de l’écoulement au passage d’une vanne 53
presque 10 fois le diamètre de la conduite, soit 500 mm avant de retrouver des lignes
de courant alignées avec la conduite. Cette distance est nettement plus pour l’incli-
naison égale à 5 ˚ puisque les lignes de courant semblent retrouver leur forme dès
le clapet passé (Figure 3.6). Nous le verrons dans le paragraphe suivant, cette obser-
vation est très importante à comprendre pour placer convenablement un capteur de
pression dans l’évaluation des pertes de charge générées par une vanne.
La figure 3.7 présente l’évolution de la variation de la pression statique de l’écou-
lement par rapport à celle observée un diamètre en amont du clapet et le long de la
conduite. On note en premier lieu une augmentation de la pression suite à un ralen-
tissement de l’écoulement à l’approche du clapet. La diminution de l’aire de passage
conduite ensuite à une augmentation de la vitesse de l’écoulement par conservation
du débit volumique à l’origine d’une dépression élevée (−500 Pa sur la figure 3.7).
On observe ensuite une décroissance linéaire de la pression statique après une dis-
tance de deux diamètres en aval du clapet. Celle-ci est associée aux pertes de charge
régulières que subit l’écoulement une fois celui-ci établi.
Figure 3.6 : Modification de l’écoulement se déplaçant de la gauche vers la droite par la vanne lorsqu’elle
est très ouverte (haut, 5˚) ou peu ouverte (bas, 65˚).
Figure 3.8 : Schéma d’un banc de mesure des performances d’une vanne (d’après Bosserman [22]).
nomètre doit être placé au moins à cinq ou huit diamètres en aval de la vanne (Figure
3.8). Cette longueur peut atteindre parfois trente diamètres. L’accumulation des pertes
de charge régulières devient donc non négligeable et doit alors être ôtées pour estimer
uniquement les pertes de charge singulières de la vanne, ∆Hv . Celles-ci peuvent être
évaluées par la formulation suivante :
L V2
∆Hv = ∆Hmes. − λ
D 2g
avec ∆Hmes. les pertes de charge mesurées à partir des deux manomètres, λ le coeffi-
cient de pertes de charge régulières de la conduite de diamètre D et L l’espacement
entre les deux manomètres.
Figure 3.9 : Schéma de la conduite de l’étude numérique illustrant la section 3.4 et définition des plans à
partir desquels sont estimées des variations de pression statique et totale.
aval du clapet (elles devraient être très proches). La différence des pressions totales,
ptot = p + 1/2ρV 2 , évaluées à l’entrée et à la sortie du domaine permet d’estimer
les pertes de charge totales dissipées sur la conduite, ∆pc,tot. 3 . Les valeurs obtenues
numériquement lorsque l’inclinaison du clapet vaut 65 ˚ sont rassemblées dans le ta-
bleau 3.1.
Le coefficient de perte de charge singulière vaut donc pour cette inclinaison du clapet
(i.e. 65 ˚ ) :
∆psing 28347, 19
ξv = = = 56, 8
1 2 499
ρV
2
En appliquant cette démarche pour toutes les inclinaisons, on trouve les valeurs
numériques du coefficient de pertes de charge singulières du tableau 3.2.
Inclinaison (°) 5 15 25 35 45 55 65
∆psing (Pa) 53 252 771 1948 4713 11 469 28 347
ξv 0,11 0,50 1,54 3,90 9,44 22,98 56,80
Tableau 3.2 : Valeurs numériques des pertes de charge singulières et du coefficient de pertes de charge
singulières de la vanne en fonction de l’inclinaison de son clapet.
3.8 Application
On souhaite dimensionner une vanne à soupape capable de générer une perte de
charge égale à 3,0 bar lorsque le débit vaut 52 l/s et au maximum 0,55 bar à 380 l/s.
Nous ferons l’hypothèse ici que les diamètres nominaux des vannes sont distancés
tous les 5 cm. Seuls les critères d’ouverture seront utilisés dans le dimensionnement
de la vanne, ceux étant liés à la cavitation et aux phénomènes transitoires ne seront
pas abordés (mais devraient l’être). On utilisera les caractéristiques de vanne de la
figure 3.2 pour la résolution de ce problème.
Critère #1 Pour le débit maximum (lorsque la vanne est totalement ouverte) la perte
de charge vaut 0,55 bar (ou ∆H = ∆p/(ρg) = 5,6 m). On va déterminer ici le diamètre
de la vanne permettant de satisfaire ce critère. D’après les figures 3.2.a et 3.2.b, les
coefficients Cd et ξ v valent respectivement :
Cd = 0, 41 et ξ v = 4, 95
∆H 0,5 5, 6 0,5
V = 2g = 2 × 9, 81 × = 4, 7 m/s
ξv 4, 95
L’aire de la section de passage, satisfaisant à la fois un débit QMax = 380 l/s et V = 4,7
m/s, s’écrit
QMax 380 × 10−3
S= = = 8, 06 × 10−2 m2
V 4, 7
√
d’où un diamètre D = 4S/π = 32,0 cm. Étant donnée la normalisation des diamètres
des vannes stipulée par l’énoncé, nous avons le choix entre les diamètres de vanne
60 Les Vannes
Critère #4 Celui-ci stipule une vanne ne devant pas être utilisée avec une ouverture
inférieure à 10%. On doit donc l’évaluer lorsque Qmin = 52 l/s et ∆p = 3,0 bar. En
nous basant sur le diamètre évalué d’après le critère précédent, la vitesse vaut V =
4Qmin /(πD2 ) = 0,54 m/s. En utilisant la relation (3.1), on obtient ξ v = 2054 et d’après
le graphique de la figure 3.2.b, l’ouverture de la vanne devrait être proche de 15%,
valeur acceptable. La vanne à soupape de diamètre égal à 35 cm satisfait donc tous
les critères utilisés ici.
Chapitre 4
Principe de fonctionnement des
pompes rotodynamiques
Historiquement, le principe de pompage de l’eau est apparu dès lors que l’homme
a su construire des habitations artificielles pour se protéger des éléments naturels.
Le besoin en eau l’oblige alors à trouver un système de transport de cette eau, du
puits ou de la rivière à son habitat. Les Grecs et les Romains fûrent les premiers à
utiliser des systèmes (rotatifs) de pompage dont la fonction était de relever l’eau dans
des aqueducs alimentant ensuite les villes et les bains. Il faudra attendre la fin du
xviiième siècle pour que les premières pompes soient réellement construites et utilisées
de façon industrielle. En effet, avant cette date, le moulin à eau restait le principal
moyen de production d’énergie et les besoins économiques étaient satisfaits par la
croissance démographique qui fournissait une main-d’œuvre bon marché. L’explosion
démographique, qui apparût dans la seconde moitié de ce siècle, obligea l’outil de
production à se mécaniser. Cette mécanisation mise en place grâce à l’apparition de
la machine à vapeur entraina le besoin en minerai situé à des profondeurs de plus
en plus importantes. Dès lors l’utilisation de pompes (à pistons) fût obligatoire afin
d’évacuer l’eau s’infiltrant dans les galeries. Il s’agit probablement là de la première
utilisation industrielle des pompes. L’un des acteurs de cette évolution est Thomas
Newcomen qui en 1712 mis au point un système utilisant la puissance de la vapeur
en puissance mécanique permettant la mise en mouvement des organes des pompes.
Ce système de pompage fût utilisé à travers l’Angleterre principalement pour évacuer
l’eau des mines. Les pompes centrifuges, dont le principe avait déjà été énoncé en 1689
par Denis Papin dans son “Acta Eruditorium”, ne furent industrialisées qu’à la fin du
XVIIIème siècle (Fig. 4.1) même si la première pompe semble être celle de Le Demours
62 Principe de fonctionnement des pompes rotodynamique
par cette vitesse de rotation. Un débit variable est présent dans le cas d’une pompe à
pistons mais l’est aussi pour des pompes à engrenages. Le fluide étant enfermé dans
le creux entre deux dents, lorsqu’il y a engrènement, la dent agit comme un piston
et chasse le fluide. A chaque position d’engrènement correspond un point de contact
différent, d’où la fluctuation de débit. Les caractéristiques principales des pompes
volumétriques sont données dans le tableau 4.1.
Performance Avantages
Quelques 100aines de m3 /h Q fixe ∀ la pression de refoulement
Pas d’émulsion des produits pompés (bière,. . .)
Adapté pour les produits visqueux
Performance Avantages
jusqu’à 105 m3 /h Plage de fonctionnement très large
Bon rendement (70 à 80 %)
Vitesse de rotation de 750 à 3000 tr/min
Inconvénients
Inadaptées aux fluides très visqueux
Production d’une pression différentielle faible
Généralités sur les pompes rotodynamiques 65
-(i) Les pompes centrifuges pour lesquelles l’écoulement entre dans la partie interne
de façon alignée avec l’axe de rotation, l’écoulement ressortant de la partie interne
selon un plan dont la normale est alignée avec l’axe de rotation (Figure 4.4).
-(ii) Les pompes axiales ont la particularité quant à elle, de part leur géométrie, de
laisser l’écoulement principalement de façon alignée avec l’axe de rotation tout au
long de son parcours dans le rotor (Figure 4.5).
-(iii) Les pompes hélico-centrifuges sont caractérisées par une entrée axiale de l’écoule-
ment mais un écoulement en sortie ayant à la fois des composantes de vitesses radiale
et axiale (Figure 4.6).
Nous le verrons plus loin, ces trois types de pompe n’ont pas les mêmes usages en
terme de débit et de quantité d’énergie délivrée à l’écoulement. Une pompe rotodyna-
mique est composée d’une roue clavetée sur un arbre, tournant à l’intérieur d’un corps
de pompe. Dans le cas des pompes centrifuges, celui-ci forme une canalisation pro-
gressivement divergente appelée volute (Fig. 4.4). Dans le cas des pompes axiales, le
corps de pompe reste généralement cylindrique. Cette roue est pourvue d’aubes dont
les formes diffèrent en fonction du type de pompe rotodynamique. Pour les pompes
centrifuges, la concavité de la courbure des aubes est dirigée à l’opposé du sens de
rotation, imprimant au liquide un mouvement de rotation. Pour les pompes axiales, la
forme en hélice oblige l’écoulement à aller de l’avant. Une roue est caractérisée par les
angles de ses aubes, son diamètre, sa largeur, et sa vitesse de rotation. Généralement
les aubes des pompes centrifuges suivent la forme d’une spirale logarithmique, les
autres formes plus simples comme par exemple des portions d’arcs circulaires condui-
sant à des rendements plus faibles [90]. Un bouchon de remplissage est placé sur la
partie supérieure de la volute permettant à la pompe de ne pas démarrer vide d’eau :
une pompe centrifuge n’est pas autoamorçante. De même, un bouchon de vidange est
situé dans la partie inférieure de la volute des pompes centrifuges.
Le grand avantage des pompes rotodynamiques réside dans leur simplicité ; leur
construction ne pose pas de problème majeur, car nulle part une grande précision
n’est demandée. En effet, on peut accepter jusqu’à un jeu de quelques millimètres
sans inconvénient.
Une pompe est généralement placée entre deux canalisations placées en amont
et en aval de la pompe appelées respectivement conduites d’aspiration et de re-
foulement (Fig. 4.7). Celles-ci n’ont généralement pas les mêmes diamètres puisque
la vitesse souhaitée à l’aspiration est plus faible que celle présente au refoulement
[66]. Cette pratique provient du fait qu’une réduction de la vitesse tend à éloigner le
phénomène de cavitation par une diminution des pertes de charge à l’aspiration [?].
66 Principe de fonctionnement des pompes rotodynamique
(a) (b)
Figure 4.4 : Pompe centrifuge − a, vue ouverte ; b, vue de coupe (D’après Bergeron [19]).
(a) (b)
Figure 4.5 : Pompe axiale − a, dessin [4] ; b, vue de coupe [7] − l’écoulement a une trajectoire rectiligne
au passage dans le rotor puis est évacué dans la conduite coudée derrière laquelle est placé le
moteur.
(a) (b)
Figure 4.8 : a, gravure d’Anthony Turner d’une vis d’Archimède [85] ; b, vis d’Archimède utilisée pour pom-
per l’eau des polders à Kinderdijk (photo d’après M.A. Wijngaarden - page internet wikipedia
concernant Kinderdijk).
met pas d’améliorer ces performances. Pour cela une analyse fine de l’écoulement au
passage du rotor et du stator doit être menée. Les méthodes numériques, de type CFD
(pour Computational Fluid Dynamics), sont utilisées depuis peu et permettent de mieux
comprendre la physionomie de l’écoulement. Cette méthode retrouve généralement
avec satisfation la hauteur manométrique apportée par les pompes (Fig. 4.10). Une
identification des améliorations de la géométrie est alors réalisée et est souvent à
l’origine d’une augmentation des performances de la pompe.
sieurs pompes sont associées entre elles pour augmenter les performances du système
de pompe, il convient alors de déterminer graphiquement la courbe équivalente “hau-
teur manométrique / débit” de ces pompes utilisées ensemble.
On appelle hauteur manométrique, la charge fournie au fluide par la pompe.
Cette charge, notée Hm , permet d’obtenir le débit souhaité à un point souvent situé à
une altitude supérieure au plan de pompage.
Ho = H1 + ∆Ha et H2 = H3 + ∆Hr
Hm = H2 − H1 = H3 − Ho + ∆Hr + ∆Ha
p3 − p o V 2 − Vo2
Hm = + z3 − z o + 3 + ∆Ha + ∆Hr
ρg 2g
∆p ∆V 2
= + ∆z + + ∆Ha + ∆Hr
ρg 2g
∆p ∆V 2
= + ∆z + + ∆H
ρg 2g
avec ∆H les pertes de charge générées par le réseau hydraulique complet. Cette re-
lation met en évidence le rôle de la charge manométrique fournie par la pompe à
travers ces quatre termes :
— la hauteur de pression, ∆p : quantité d’énergie nécessaire pour vaincre une
pression aval ou une vide amont à la pompe. Cette hauteur est généralement
mesurée sur les surfaces libres du bassin d’alimentation et du réservoir ali-
menté. Si les deux surfaces libres sont caractérisées par la même pression (e.g.
la pression atmosphérique), la pompe n’a pas à fournir d’énergie pour lut-
ter contre cette contribution. Ceci est aussi le cas pour un réseau hydraulique
fermé.
— la hauteur géométrique, ∆z : il s’agit de la hauteur géométrique que devra
réaliser le fluide durant son trajet. Cette quantité est généralement estimée via
Courbes de la charge manométrique en fonction du débit 71
(a) (b)
Figure 4.10 : Détermination des performances d’une pompe centrifuge à partir d’une méthode numérique
(a, maillage des domaines rotor et stator ; b, courbe de la charge manométrique en fonction
de la vitesse angulaire des pales), d’après Stel [79].
(a) (b)
Figure 4.12 : Courbes caractéristiques d’une pompe centrifuge (a, forme classique ; b, autre forme possible).
Figure 4.13 : Couplage de deux pompes en série (a, couplage en série d’une pompe immergée et d’une
pompe de surface ; b, couplage en série de deux pompes de surface), d’après Duclos [33] ; c,
Couplage de deux pompes en parallèle (chaufferie ECAM).
Dans la pratique, on choisit une pompe pour qu’elle fournisse un débit à une pres-
sion donnée (que l’on convertit en hauteur manométrique). La pompe et le circuit
étant caractérisés par leurs courbes, on choisira la pompe la plus proche du point
de fonctionnement souhaité. Si celui-ci est trop éloigné des caractéristiques de notre
pompe, on peut alors soit utiliser une pompe “surdimensionnée” (i.e., délivrant un
débit trop important) et introduire une vanne dans le circuit permettant d’introduire
une nouvelle perte de charge et baisser ainsi le débit. Une autre solution est de rogner
les aubes de la pompe (Cf. paragraphe 4.13.1). Cette opération est possible sur la plu-
part des pompes, c’est pourquoi les constructeurs de pompes associent non pas une
courbe mais une zone de hauteur, définissant alors une plage d’utilisation (Fig. 4.25).
(a) (b)
on ne pourra pas négliger les pertes de charge. L’écoulement créé par la pompe ali-
mente un réservoir situé à une altitude plus élevé que le plan de pompage. Toutes les
courbes des conduites du réseau sont illustrées sur la figure 4.16.b.
Les deux approches permettant d’estimer le point de fonctionnement de la pompe
sont illustrées sur la figure 4.17. Dans le cas (a), le système alimentant est composé de
la pompe et de sa conduite d’aspiration. Les pertes de charges correspondant à une
dissipation irréversible de l’énergie de l’écoulement, celles-ci sont donc ôtées de la
charge délivrée par la pompe. Le sytème à alimenter est composé alors de la conduite
de refoulement seule. Le point d’intersection entre les courbes des deux systèmes,
tracées en gras, correspond au point de fonctionnement du système. On trouve le
point de fonctionnement de la pompe, en ajoutant à partir de ce point les pertes de
charges de l’aspiration (i.e. on se déplace sur une droite de débit constant jusqu’à at-
teindre la courbe caractéristique de la pompe). La deuxième approche, illustrée sur la
figure 4.17.b, propose de composer le système à alimenter des conduites d’aspiration
et de refoulement. Les deux conduites sont traversées par le même débit, on addi-
tionne alors les pertes de charge. Le système alimentant est dans ce cas caractérisé
seulement par la courbe de la pompe. Cette fois, le point de fonctionnement est direc-
tement estimé par l’intersection des deux courbes.
(a)(b)
Figure 4.16 : a, réseau hydraulique de l’estimation du point de fonctionnement d’une pompe à partir
de deux approches différentes ; b, courbes caractéristiques des conduites d’aspiration et de
refoulement, et de la pompe.
(a) (b)
Figure 4.17 : Détermination du point de fonctionnement de la pompe utilisée dans le réseau hydraulique de
la figure 4.16 à partir de deux méthodes. a, la courbe du système alimentant est construite en
enlevant les pertes de charge de la conduite d’aspiration de la charge délivrée par la pompe,
et la courbe du système à alimenter correspond à la courbe caractéristique de la conduite de
refoulement ; b, la courbe du système alimentant correspond à la courbe seule de la pompe,
et le système à alimenter correspond aux pertes de charge des conduites d’aspiration et de
refoulement. Dans les deux approches, on retrouve les mêmes points de fonctionnement de
la pompe (•).
Puissances et rendements 77
4.5.1 Puissances
Les puissances impliquées dans l’étude d’une pompe sont (1) la puissance commu-
niquée au fluide ou puissance utile (Pu = ρgHm Q), (2) la puissance de la pompe
ou puissance hydraulique (P p = ρgH p Q), (3) la puissance absorbée par la roue
(Pr = ρgH p ( Q + q)), et (4) la puissance absorbée sur l’arbre (P a = Cm ω) avec Cm
le couple délivré par le moteur entrainant la pompe, q le débit de fuite.
4.5.2 Rendements
Le rendement de la pompe ou rendement global est le rapport de la puissance
communiquée au fluide et de la puissance absorbée sur l’arbre moteur :
Pu ρgHm Q
η= =
Pa Cm ω
Ce rendement est donné dans les catalogues des fabricants, et peut être affecté par
différentes pertes dans la pompe : (1) les pertes hydrauliques, fonctions du débit, (2)
les pertes par fuite s’il existe un débit de fuite q (le débit dans la roue doit être de
Q + q). Le débit de fuite augmente avec la pression et diminue avec la viscosité. (3)
les pertes mécaniques fonctions de la vitesse de rotation. On peut alors décomposer
le rendement global :
— Rendement manométrique : rapport de la puissance communiquée au fluide
et de la puissance de la pompe
Pu ρgHm Q Hm H p − ∆H p
ηm = = = =
Pp ρgH p Q Hp Hp
∆H p
= 1−
Hp
Pp ρgH p Q 1 Qréel
ηv = = = ≡
Pr ρgH p ( Q + q) 1 + q/Q Qthéorique
— Rendement mécanique : rapport de la puissance absorbée par la roue et de la
puissance du moteur. Ce rendement tient compte des pertes dans les paliers
78 Principe de fonctionnement des pompes rotodynamique
de l’arbre de rotation qui supportent la pompe et des frottements qui ont lieu
entre les aubes et les parois fixes.
Pr ρgH p ( Q + q)
ηméca = =
Pa Cm ω
Le rendement global peut dès lors s’écrire η = ηm × ηv × ηméca avec un rendement
caractéristique des performances internes à la pompe (ηm × ηv ), et un rendement lié
à la partie externe à la pompe (ηméca ). Pour des débits nuls, aucune puissance n’est
communiquée au fluide (Pu = 0) et pour un débit très important la puissance utile
est nulle car les pertes de charge sont trop importantes. Finalement, la courbe de
rendement a une allure en cloche. La pompe devra travailler évidemment dans des
conditions pour lesquelles le rendement est maximum, afin entre autre de baisser
les coûts d’utilisation. Ce dernier aspect est important et à ne pas négliger car les
pompes consomment près de 20 % de l’énergie électrique utilisée par l’industrie [37].
Un abaissement de la puissance consommée peut être atteint par amélioration de la
pompe elle-même mais aussi par le moteur associé.
Ainsi d’après cette relation, si on considère une famille de pompes, telle que cela a
été défini précédemment, les valeurs numériques de la charge manométrique avec le
débit dépendront des dimensions des pompes et de leur vitesse de rotation, alors que
2 3
les valeurs numériques de gH/( Dω ) avec Q/ D ω seront identiques quelles que
soient les dimensions de la pompe et sa vitesse de rotation. La courbe µ = f (δ) ca-
ractérise donc une famille de pompe, quelles que soient la vitesse de rotation et les di-
mensions des pompes de cette famille. La figure 4.20 présente en ce sens des résultats
numériques de l’évolution du coefficient de hauteur pour une pompe centrifuge [79].
Les hauteurs manométriques ont été obtenues initialement, de façon numérique, à
partir de différentes vitesses de rotation de la pompe, 612 tr/min., 806 tr/min., 1000
tr/min. et 1150 tr/min. (Figure 4.10.b). Lorsque les coefficients de Rateau sont uti-
lisés cette fois, on note des points faisant parti de la même courbe caractétistique,
2 3
gH/( Dω ) = f (Q/ D ω ). Ceci confirme ce qui a été dit plus haut. Si une autre
pompe de la même famille était étudiée, sa courbe µ= f (δ) serait confondue avec les
quatre présentes sur la figure 4.20. Ainsi l’intérêt à utiliser les coefficients de débit et
de hauteur repose sur la manipulation d’une unique courbe, représentant alors une
famille de pompes similaires quelles que soit leur vitesse de rotation et leur dimen-
sion. Il est toujours possible évidemment de retrouver la courbe caractéristique d’une
pompe particulière de la famille de pompe, à partir de la courbe µ = f (δ), connaissant
ses dimensions et la vitesse de rotation d’utilisation.
µ3 ( gH )3 δ1/2 ( gH )3/4
= ou ω= (4.2)
δ2 ω 4 Q2 Q1/2 µ3/4
Par définition, la vitesse spécifique de cette famille de pompes est la vitesse de rota-
tion d’une pompe délivrant un débit de 1 m3 /s à une hauteur de 1 m. Ainsi, d’après
la relation précédente, ωs = δ1/2 g3/4 /µ3/4 qui peut être réécrite à partir de la relation
(4.2) :
Q1/2
ωs = ω 3/4 (4.3)
H
avec ω exprimée en rad/s, H en m, et Q en m3 /s. Si on prête attention à l’expression
(4.4), on note que la vitesse spécifique n’est pas sans dimension (i.e. m−3/4 × s−3/2 ). Sa
valeur numérique dépend donc du système d’unités utilisé pour exprimer la vitesse
angulaire, le débit et la charge manométrique. Cela rend d’autant plus confus sa ma-
nipulation que les unités varient en fonction des pays. Ainsi aux Etats-Unis, la vitesse
angulaire, le débit et la charge manométrique sont exprimés en tr/min, en gal US (=
3,785 l) et en pied (= 30,48 cm) respectivement alors que ces unités deviennent tr/min,
gal Imp (= 4,546 l) et pied au Royaume Unis. Il est donc parfois préférable d’ajouter
l’accélération de pesanteur dans l’expression précédente et exprimée dans les unités
adaptées afin de rendre l’expression de la vitesse spécifique sans dimension. On par-
lera généralement de coefficient de vitesse spécifique (Figure 4.21), il s’écrit dans ce
cas :
Q1/2
ωs = ω (4.4)
( gH )3/4
Q1/2
NS = N (4.5)
H 3/4
avec N exprimée en tr/min, H en m, et Q en m3 /s. Les valeurs numériques de la
vitesse spécifique correspondent donc aux gammes proposées sur la figure 4.22. À
chaque valeur numérique de la vitesse spécifique Ns , il lui correspond une forme
de roue particulière, indépendante de la taille de cette roue. On note par exemple,
d’après la relation (4.5), qu’à valeurs égales de débit et de vitesse de rotation, les
Similitudes 83
Figure 4.22 : Formes de la roue en fonction de la vitesse spécifique, N × Q1/2 /H 3/4 , estimée à partir
de la relation (4.5) (les valeurs de NS sont estimées à partir de N exprimé en tr/min., Q en
m3 /s et H en m). On reconnait la pompe centrifuge pour les valeurs numériques faibles de
la vitesse spécifique et la pompe axiale pour laquelle la valeur numérique est élevée.
pompes de grande vitesse spécifique fonctionnent avec une hauteur d’élévation faible
(i.e. les pompes axiales). Il existe une progression continue entre la pompe centrifuge
pure radiale jusqu’à l’hélice de type axial. Notons que la forme n’est ni un avantage,
ni un inconvénient en soi pour l’utilisateur ou le constructeur.
On peut se faire une idée de la vitesse spécifique d’une pompe en regardant sa
roue. En effet, les roues fortement radiales ont des aubes à faible passage de liquide.
La longueur des aubes est importante (le rapport des diamètres extérieur et d’aspi-
ration) excède deux et peut atteindre jusqu’à cinq. Les roues héliocentrifuges ont un
diamètre extérieur semi-conique et des aubes larges. Le rapport des diamètres est
voisin de l’unité. Les hélices quant à elles ont une largeur de passage très importante
pour un rapport de diamètre égal à l’unité. Ainsi pour des rapports de diamètres
diminuant et des largeurs de passage augmentant, la vitesse spécifique croı̂t.
Lorsque la valeur numérique de la vitesse spécifique augmente, i.e. lorsque l’on
passe de la famille des pompes centrifuges pour se diriger vers la famille des pompes
axiales, on observe une modification de la courbe caractéristique et une variation dans
la valeur maximale possible du rendement (Figure 4.23).
4.7 Similitudes
Les lois de similitude, pour les pompes rotodynamiques, sont généralement uti-
lisées pour estimer la courbe caractéristique d’une pompe à une vitesse ω ′ connaissant
celle mesurée à la vitesse ω. Ceci est rendu possible dès lors que les écoulements ob-
servés pour ces deux vitesses, i.e. ω et ω ′ , sont similaires, c’est à dire ayant des formes
identiques. Des lois de similitude existent aussi pour des modifications du diamètre
d’entrée mais ne seront pas étudiées ici.
Le débit éjecté par les aubes d’une pompe étant linéaire avec la vitesse
de rotation,
3
si on considère une pompe rotodynamique, la fonction f Q/ D ω de la relation
84 Principe de fonctionnement des pompes rotodynamique
Figure 4.23 : Évolution de la courbe de charge manométrique et de la valeur maximale que peut atteindre
une pompe de la famille considérée avec la vitesse spécifique (Document SIHI).
(4.1) devient alors constante. En conséquence la quantité gH/ω 2 reste aussi constante.
Ainsi pour une variation de la vitesse de rotation de ω à ω ′ , un point localisé par ses
coordonnées (Q,H) sera déplacé aux nouvelles coordonnées (Q′ ,H ′ ) tout en satisfai-
sant les relations de similitude suivantes :
′ 2
Q′ ω′ H′ ω
= =
Q ω H ω
Ces relations pouvant être évidemment remplacées par d’autres faisant intervenir la
vitesse de rotation exprimée en tr/min, notée N. Intéressons-nous à l’utilité de ces
formulations et considérons pour cela la courbe caractéristique d’une pompe tour-
nant à une vitesse de rotation N = 1450 tr/min (Fig. 4.24). On souhaite connaı̂tre sa
caractéristique pour une vitesse plus faible, N’ = 1200 tr/min. Le rapport des vitesse
et ce même rapport au carré valent respectivement :
2
N′ N′
= 0, 828 et = 0, 685
N N
Q′ = 0, 828 × Q, et H ′ = 0, 685 × H
Tableau 4.3 : Points de la caractéristique de la pompe de la figure 4.24 pour deux vitesses de rotation (le
premier nombre représente le débit, le deuxième la hauteur manométrique).
1450 tr/min (30 ; 38) (25 ; 72) (20 ; 100) (15 ; 122) (10 ; 138)
1200 tr/min (24,8 ; 25,7) (20,7 ; 49,2) (16,6 ; 68,5) (12,4 ; 83,5) (8,3 ; 94,2)
Le débit : celui-ci est spécifié par le responsable du projet et découle des besoins si-
multanés à assurer (besoins domestiques, appareils à alimenter, puissance calorifique
d’une installation de chauffage dont la chute de température est fixée,...). Malgré le
fait que les statistiques (Tableau 4.4) ne font apparaı̂tre que des données moyennées,
86 Principe de fonctionnement des pompes rotodynamique
il ne faut tenir compte que des débits instantanés. Pour une utilisation domestique,
on admet en général comme débit de la pompe 1/3 de la consommation totale jour-
nalière.
Les catalogues des constructeurs : on peut maintenant s’intéresser aux pompes pro-
posées par les constructeurs dans la famille la plus adaptée. Ces derniers proposent
aux clients des plages d’utilisation de leurs pompes pour lesquelles le rendement est
Accessoires pour une installation hydraulique 87
correct et ceci pour des vitesses de rotation fixées. La figure 4.25 représente les plages
d’utilisation d’une pompe centrifuge fonctionnant à une vitesse de rotation égale à
1450 tr/min. Il s’agit ici d’une plage et non pas d’une courbe car un rognage des
aubes, c’est à dire une réduction du diamètre, est toujours possible.
Il est très difficile d’estimer à l’avance les pertes de charge d’un réseau hy-
draulique. Les points de fonctionnement sur site et par le calcul seront pro-
bablement différents. La pompe risque alors de ne pas fonctionner à son rende-
ment le meilleur. Pour une optimisation, on peut adapter la pompe au réseau en
réduisant par exemple le diamètre de la roue à la valeur calculée après essais.
− Clapet anti-retour : il doit être installé sur la conduite de refoulement pour éviter
tout reflux de fluide en cas d’arrêt soudain du moteur. Il est aussi préférable de choi-
sir des types de clapet avec ressort intérieur et avec obturateur à ogive, aidant ainsi à
atténuer le phénomène de coup de bélier.
− Vanne : non seulement elle permet le démontage de la pompe sans vider l’installa-
tion mais elle sert aussi pour le démarrage du groupe et pour le réglage de débit.
88 Principe de fonctionnement des pompes rotodynamique
Figure 4.26 :
Clapet de pied avec
crépine.
4.10 La cavitation
La cavitation est le phénomène qui se produit quand, en un point de l’écoulement,
la pression absolue locale devient égale à la pression de vapeur saturante du liquide,
notée pv ( T ). La baisse locale de cette pression peut provenir :
— d’une augmentation de la vitesse, généralement consécutive à une modification
locale de l’aire de la section de passage de l’écoulement ;
— d’une augmentation d’altitude entrainant une diminution de la pression at-
mosphérique ;
— d’une perte de charge générée par une singularité (e.g. vanne, diaphragme,
etc.).
Ce phénomène n’est pas à prendre à la légère tant l’impact à la fois sur le rendement
de fonctionnement, mais surtout sur la survie du matériel, est important. Afin d’éviter
son apparition, il existe un critère basé sur la connaissance du réseau hydraulique
amont à la pompe et sur les caractéristiques internes de celle-ci.
La suite de ce paragraphe se propose dans une première partie de détailler les
changements de phase impliquées dans ce processus. On terminera enfin par la présen-
tation du critère cité plus haut.
Changements de phase
(a) (b)
Figure 4.27 : (a), Nomenclature des changements de phase ; (b) Courbe de saturation.
Notion de NPSH
Afin d’éviter les nuisances citées plus haut, il faut qu’en tout point du circuit, la
pression absolue totale reste supérieure à la tension de vapeur saturante, la première
s’écrivant :
1
pabs, tot = ps + ρV 2 (4.6)
2
où ps , ρ et V sont la pression statique, la masse volumique et la vitesse de l’écoulement
au point étudié. Généralement la masse volumique du bassin d’alimentation est considé-
rée dans cette formulation, quelle que soit la localisation à laquelle on souhaite évaluer
la pression absolue totale. Néanmoins dans le cas de pompe à eau chaude, il peut
exister des écarts de température entre la tuyauterie et le bassin de telle sorte que les
valeurs numériques des masses volumiques ne soient plus confondues.
La pression absolue étant quantifiable en tout point du circuit, on est donc capable
d’estimer la marge que le fluide dispose par rapport à la cavitation, i.e. pabs, tot − pv ( T ).
Cette quantité variable dans le circuit, est appelée NPSH, pour Net Positive Suction
Head (Hauteur de Charge Nette Absolue). Cette valeur traduit l’énergie encore dis-
ponible dans le fluide, au-dessus du seuil de tension de vapeur saturante : c’est une
réserve d’énergie. À l’entrée de la pompe, cette quantité va atteindre une valeur mi-
nimale. En ce point du circuit (et uniquement), on parle de NPSH disponible. Cette
La cavitation 91
Figure 4.29 : Évolution de la pression dans la pompe et trajectoire des particules fluide dans la pompe (les
évolutions à travers les différents organes de la pompe ne sont pas respectées ici).
pabs, tot − pv ( T ) p o − p1
− >0 ou NPSHdisponible > NPSHrequis
ρg ρg
cas a cas b
Figure 4.30 : Réseaux hydrauliques. cas a, pompe montée en charge ; cas b, pompe montée en aspiration.
94 Principe de fonctionnement des pompes rotodynamique
ps − pv ( T )
NPSHdispo = + h − ∆Ha
ρg
=⇒ Cas b : Pompe montée en aspiration : un liquide est prélevé d’un bassin, pour
lequel cette fois-ci la pression à sa surface est égale à la pression atmosphérique et
situé à une hauteur h en-dessous de l’entrée de la pompe. Les indices b et e seront
utilisés pour signifier le bassin et l’entrée de la pompe respectivement. L’application
du théorème de Bernoulli nous permet encore une fois d’exprimer la pression absolue
totale à l’entrée de la pompe, c’est-à-dire :
pe V2 p
+ e = b + (zb − ze ) − ∆Ha ,
ρg 2g ρg
patm − pv ( T )
NPSHdispo = − h − ∆Ha
ρg
=⇒ Cas général : d’après les deux cas précédents, on peut exprimer de façon générale
le NPSH disponible :
ps − p v ( T )
NPSHdisp = + hgéo − ∆Ha (4.9)
ρg
où ps est la pression à la surface du réservoir alimentant la pompe. Dans le cas d’un
bassin ouvert, cette pression est égale à la pression atmosphérique. S’il s’agit d’un
réservoir fermé pressurisé, il faut alors additionner la pression (relative) lue sur le
manomètre monté sur le réservoir à la pression atmosphérique. hgéo correspond à
la différence de hauteur entre le niveau d’alimentation et l’entrée de la pompe. La
pompe fonctionnant en charge, cette hauteur sera considérée comme positive alors
que dans le cas pour lequel un liquide est pompé d’un niveau plus bas, la hauteur
devient négative. ∆Ha représente les pertes de charge dans la conduite d’aspiration.
On note d’après la relation (4.9) que le niveau de pression disponible fait intervenir
trois composantes, qui ajoutent ou au contraire qui ôtent de l’énergie au système.
Ainsi pour augmenter le NPSH disponible, on peut agir sur différents paramètres :
(1) se placer si possible dans le cas d’une pompe montée en charge et augmenter au
plus la hauteur d’alimentation ; (2) pressuriser si possible l’alimentation ; (3) réduire
au maximum les coudes dans la conduite d’aspiration et optimiser le diamètre de
celle-ci afin de réduire les pertes de charge ; (4) diminuer la température du liquide
afin d’abaisser la pression de vapeur saturante (Fig. 4.27.b).
Pour l’eau, aux conditions de température et de pression normales, i.e. 20 ˚ C et
1 013 mbar, et lorsque le bassin a sa surface libre à pression atmosphérique, on a
La cavitation 95
Une autre solution pour éviter la cavitation dans la roue repose sur l’augmentation
de la pression de l’écoulement dans sa région d’entrée à l’aide d’un dispositif appelé
inducteur (ou parfois roue de gavage). Il s’agit-là d’une roue axiale placée à l’aspiration
d’une roue centrifuge et montée sur le même arbre (Figure 4.31). Cette roue axiale
augmente la charge de l’écoulement à son passage (et donc la pression) et rendant
donc plus favorable la roue centrifuge par rapport à la cavitation. Par translation
des effets, l’inducteur est souvent soumis à la cavitation et son dimensionnement
doit alors être réalisé attentivement. La courbe de la figure 4.31 met en évidence une
influence nette de l’inducteur sur les performances de la roue en régime cavitant.
L’inducteur diminue de façon non négligeable les valeurs du NPSH requis sur une
plage de débit allant de 0, 3 × Qn à 1, 4 × Qn où Qn est le débit nominal. À ce débit, le
gain de NPSH requis représente 58% du NPSH requis de la pompe centrifuge seule.
Figure 4.31 : Gauche : photo d’un inducteur placée en amont d’une turbo-pompe de moteur-fusée (d’après
Morel [59]) ; droite, influence de la présence d’un inducteur sur le NPSH requis d’une pompe
centrifuge [67].
Analyse interne d’une pompe centrifuge 97
Figure 4.33 : Caractéristiques de la cavitation dans une pompe centrifuge pour différentes valeurs
numériques du coefficient de débit (d’après Govinda et al. [47]).
(a) (b)
Figure 4.36 : Coupes frontale (a) et méridienne (b) d’une pompe centrifuge.
Analyse interne d’une pompe centrifuge 101
(a) (b)
Figure 4.37 : a, vue frontale de la roue d’une pompe centrifuge et présentation des angles β ; b, trajectoires
absolue et relative de l’écoulement dans les aubes de la roue.
La vitesse du fluide en un point M de l’écoulement peut être définie (1) par rapport
−
→
à un trièdre de référence fixe (vitesse absolue, V ), et (2) par rapport à un trièdre lié
−
→
aux organes mobiles (vitesse relative, W ) (Fig. 4.38.a). Le passage du domaine fixe au
domaine mobile est effectué à l’aide de la composition des vitesses :
−
→ − → −→
V = U +W (4.10)
−
→
où U est la vitesse d’entraı̂nement orientée dans le sens de rotation. Cette vitesse est
circonférencielle et est égale en valeur absolue à rω, où r est la distance de la particule
considérée par rapport à l’axe de rotation et ω la vitesse angulaire de rotation. La
composition des vitesses donne lieu au tracé d’un triangle des vitesses (Fig. 4.39.a).
On notera α et β respectivement les angles que font les vitesses absolue et relative
avec la vitesse d’entraı̂nement.
−
→− →→
On définit un repère orthonormé ( i , j ,− z ) lié aux aubes en mouvement (Fig.
−
→
4.39.b). Le vecteur i est orienté radialement sortant et est donc placé dans un plan
−
→ −
→ −→ − →
méridien 4 , le vecteur j est aligné selon U , et −→z tel que − →
z = i ∧ j . Ces vecteurs
sont utiles à des décompositions des vitesses introduites plus haut.
L’entrée et la sortie de la roue sont généralement indiquées respectivement par
les indices 1 et 2. Ainsi, les vitesses d’entrainement en ces points valent U1 = ωr1 et
U2 = ωr2 .
−
→ −
→
Les vitesses V et W sont décomposées en deux termes : (1) une composante
−→ −
→
projetée sur j et donc alignées avec U dite tangentielle ou orthoradiale, Vu ou Wu ,
mesurée positivement dans le sens de rotation ; et (2) une composante projetée sur
−
→
i dı̂te radiale ou méridienne, Vr et Wr , ou Vm et Wm . Il faut bien comprendre que
−→
toute composante projetée selon le vecteur U est liée au transfert d’énergie, alors que
4. Un plan méridien d’une géométrie axisymétrique contient l’axe de rotation.
102 Principe de fonctionnement des pompes rotodynamique
(a) (b)
Figure 4.38 : Cinématique dans la roue d’une pompe centrifuge. (a), de façon généralisée ; (b), en entrée
et en sortie des aubes.
(a) (b)
Figure 4.39 : Définition du triangle des vitesses (a) et des vecteurs unitaires (b).
les composantes méridiennes sont quant à elles liées au débit (massique) traversant
−
→ →
les roues mobiles, ρS2 W 2 .−n 2 avec − →
n 2 sortant radialement de la surface d’éjection S2
placée au rayon r2 . Les différentes relations entre ces projections sont les suivantes :
( −
→ −→
Vr = Wr Vr = Vsinα = Wsinβ
−
→ −→ − → ou
Vu = Wu + U Vu = Vcosα = Wcosβ + U
avec U1 un vecteur α1 très proche de 90 ˚ , la composante Vu1 étant donc très faible.
D’autre part, l’angle β 2 est appelé angle de sortie de la roue. Sa valeur est comprise
entre 130 ˚ et 165 ˚ .
La suite du document se propose de formuler analytiquement la charge délivrée
par la roue. On s’attend, d’après ce qui a été dit plus haut, à noter dans cette expres-
sion uniquement des termes de transfert d’énergie.
dW 2 dp
= − gdz − + rω 2 dr
2 ρ
dp dU 2
= − gdz − +
ρ 2
Il vient donc
d (W 2 − U 2 ) dp
+ dz + =0
2g ρg
L’équation de Bernoulli en mouvement relatif (fluide incompressible) est donc obte-
nue à partir d’un bilan sur un filet de fluide :
W 2 − U2 p
+z+ = constante
2g ρg
Cette quantité est conservée entre l’entrée (indice 1) et la sortie (indice 2) des aubes
mobiles, en considérant d’autre part la relation W 2 − U 2 = V 2 − 2UVu :
V22 − V12 p − p1 U2 Vu2 − U1 Vu1
+ 2 + z2 − z1 =
2g ρg g
Le terme de gauche caractérise des transferts d’énergies cinétique (∆V 2 /2g) et d’énergie
potentielle (∆p/ρg + ∆z), lors du passage du fluide dans la roue mobile. Cette quan-
tité correspond donc à la charge d’Euler telle qu’elle a été définie précédemment.
Par identification le terme du membre de droite est une autre forme d’écriture, plus
pratique, de la charge d’Euler :
U2 Vu2 − U1 Vu1
HE =
g
Analyse interne d’une pompe centrifuge 105
Le bilan de quantité de mouvement sera cette fois encore appliqué dans un repère
relatif, repère dans lequel l’écoulement est permanent. Dans le cas où ce bilan est
appliqué dans un repère absolu, les termes d’instationnarités difficiles à estimer dans
ce cas doivent être pris en compte (dérivée temporelle dans la relation (1.14) page
13). La charge d’Euler traduisant la puissance communiquée au fluide (ρgQHE ) par
la roue, et cette puissance étant directement liée à la puissance absorbée par l’arbre
(Cω, si on considère la puissance mécanique totalement transférée au fluide), il nous
faut donc exprimer le couple via l’effort exercé par les aubes sur le fluide pour obtenir
in fine la charge d’Euler.
Considérons le schéma de la figure 4.41, et intéressons-nous à l’effort exercé par
l’aube de gauche sur le fluide au point P. Dans le cas où les effets visqueux sont
négligés, cet effort élémentaire s’écrit
− p−
→
n p dS
où −
→
n p est la normale sortante à la paroi, et dS la surface élémentaire autour de ce
point. Cet effort local crée un couple élémentaire autour de l’axe de rotation
−
→ →
dC = −
r ∧ − p−
→
n p dS
avec ||−
→
r || la distance radiale au point P. Le couple total exercé par les aubes en mou-
vement sur le fluide s’écrit donc en intégrant l’expression précédente sur la surface
des deux aubes en contact avec le fluide :
−
→
Z
C = −
→
r ∧ − p−
→
n p dS
Sp
Nous nous plaçons dans cette démonstration dans l’hypothèse énoncée plus haut.
Toutes les particules fluides sont dès lors soumises aux mêmes efforts extérieurs quelle
que soit leur position. Les effets que pourraient avoir le voisinage des aubes sont donc
négligés. Ceux-ci ont une répercussion sur l’écoulement, et leurs conséquences sont
discutées dans le paragraphe 4.11.9. Le bilan de moment de quantité de mouvement
s’écrit dans notre cas :
Z
−−→ −
→ − →− →
Z
−−→ −→
Z
−−→ −
→ −−−−→ −−−−→
OM ∧ ρ W ( W. n )dS = OM ∧ σ n dS + OM ∧ ρ f v dv + M(Fe ) + M(Fc )
S S D
(4.11)
Ce bilan est appliqué sur le volume de contrôle proposé sur la figure 4.42 défini par la
région de passage des particules entre deux aubes consécutives. L’entrée et la sortie du
domaine sont caractérisées respectivement par des indices 1 et 2 ayant des normales
unitaires sortantes −→
n1 et −
→
n2 purement radiales, la paroi étant identifiée par l’indice p.
Le premier terme de l’expression (4.11) représente le moment du débit de quantité
de mouvement, alors que les termes du membre de droite sont respectivement les
moments des efforts surfaciques, des efforts volumiques, et des efforts massiques
d’inertie d’entraı̂nement et de Coriolis définis par les expressions suivantes :
−−−−→ Z
−−→ −
Z
M(Fe ) = − OM ∧ γe dm = − −
→ →
r ∧ [−
→
ω ∧ (−
→
ω ∧−
→
r )]dm (4.12)
D D
−−−−→ Z
−−→ −
Z
−
→
M(Fc ) = − OM ∧ γc dm = −2 −
→ →
r ∧ [−
→
ω ∧ W ]dm (4.13)
D D
avec les termes placés aux extrémités du membre de gauche nuls car −
→
r colinéaires à
−
→
n1 et −
→
n2 . Le moment du débit de quantité de mouvement de la relation (4.11) s’écrit,
en tenant compte des notations de la figure 4.42 :
−−→ −→ − →→ −→ −→ → −→ −→ →
Z Z Z
OM ∧ ρ W ( W.−
n )dS = −
→
r ∧ ρW1 (W1 .−
n 1 )dS + −
→
r ∧ ρW2 (W2 .−
n 2 )dS
S S1 S2
−→ −→ →
Z
+ −
→
r ∧ ρW p (W p .−
n p )dS
Sp
L’intégrale en paroi devient nulle car la vitesse est tangente à cette paroi. Les deux
−
→
autres termes font apparaı̂tre la composante − →r ∧ Wi = ri Wui − →z constante sur les
surfaces d’intégration S1 et S2 . On fait apparaı̂tre d’autre part le débit massique“dans
le sens de l’écoulement” et défini par la relation suivante :
−→−→
Z
Qm = ρ W. N dS
S
−
→ −
→
où N 1 = −−
→
n 1 et N 2 = −
→
n 2 . Le moment du débit de quantité de mouvement s’écrit
donc :
−−→ −
→ −→→
Z
OM ∧ ρ W ( W.− n )dS = Qm (r2 Wu2 − r1 Wu1 ) −
→
z
S
Les moments des efforts d’inertie s’écrivent :
−−−−→ − → −−−−→ Z
M(Fe ) = 0 et M(Fc ) = −2ω rWr −
→
z dm
D
Notons d’autre part que dm = ρdV = ρA(r )dr, avec A(r ) la surface de révolution
située à la distance radiale r de l’axe de rotation (Wr est normale à cette surface).
L’expression du moment de quantité de mouvement projeté selon z s’écrit finalement :
Z r2
C = Qm (r2 Wu2 − r1 Wu1 ) + 2ρω A(r )rWr dr
r1
U2 Vu2 − U1 Vu1
HE = (4.15)
g
108 Principe de fonctionnement des pompes rotodynamique
Pour une inclinaison d’aubes donnée et à vitesse de rotation fixée, il apparaı̂t que
la charge d’Euler est une fonction linéaire décroissante du débit. La comparaison avec
l’expérience fait apparaı̂tre des différences (Fig. 4.43) dont l’origine est connue. Dans
la zone 1, des phénomènes de recirculation dans la pompe dégradent son rendement,
et des pertes hydrauliques existent dans la zone 3.
vitesse de rotation donnée (ω), il peut être intéressant de connaı̂tre celle-ci pour une
autre vitesse de rotation (ω ′ ). La vitesse d’entraı̂nement s’écrit U2 = r2 ω, d’où :
2 2 BQ
gHE = Aω + BQω ou gHE /ω = A + (4.17)
ω
Ainsi nous retrouvons les lois de similitude abordées dans la section 4.7.
compris dans l’intervalle angulaire [0, π ]. La figure 4.44 représente les deux cas pour
lesquels l’angle β est obtu ou aigu. Pour les angles situés dans l’intervalle [0, π/2[,
la fonction tangente est positive, alors que dans l’intervalle ]π/2, π ] elle est négative.
De ce fait, dans le premier cas, la charge augmente avec le débit, alors qu’elle est
décroissante dans le deuxième cas (Fig. 4.45). On notera d’ailleurs que la figure 4.45
représente la charge d’Euler et ne tient pas compte de la volute qui est généralement
adaptée à un angle d’inclinaison obtu. Un sens de rotation désadapté entraı̂ne alors
une mauvaise circulation du fluide générant de fortes pertes de charge, diminuant de
façon catastrophique le rendement. D’autre part, pour des raisons non développées
ici, il est préférable d’utiliser des courbes de pompes décroissantes plus stables que
celles ayant des courbes de pompes croissantes. Pour toutes ces raisons, les pompes
ayant des aubes incurvées vers l’arrière sont généralement préférées.
(a) (b)
Figure 4.44 : Influence du sens de rotation sur le triangle des vitesses : a, β est obtu ; b, β est aigu.
Dans les deux cas, la vitesses d’entraı̂nement est orthoradiale et est orientée dans le sens de
rotation des aubes, le vecteur relatif est tangente à la trajectoire relative (en pointillé sur la
figure). La vitesse absolue est ensuite obtenue par composition des deux vecteurs précédents.
orthoradiale de la vitesse absolue d’aubes, Vu2 , comparativement à celle qui serait ob-
servée si l’écoulement était parfaitement guidé, Vu′ 2 (Fig. 4.47). Ainsi, le glissement
a pour conséquence une charge d’Euler, gHE = U2 Vu2 , diminuée par rapport à celle
observée, encore une fois, si l’écoulement était parfaitement guidé, gHE = U2 Vu′ 2 .
Ce phénomène est quantifié par le facteur de glissement, proposé par Stodola [80]
en 1927 :
Vu
σ = ′2 (4.18)
Vu2
évoluant entre 0 et 1, les valeurs numériques étant généralement proches de 0,8−0,9.
Le facteur de glissement représente donc la qualité du transfert d’énergie entre les
aubes et l’écoulement. Il représente un critère important à considérer lors de la concep-
tion interne d’une pompe. Plusieurs modèles ont été proposés permettant d’approcher
au mieux la valeur numérique du facteur de glissement. Stodola [80] propose en 1927
la formulation suivante :
π cosβ′2
σ = 1− (4.19)
Z 1 − ϕ2 tanβ′2
avec ϕ2 = Vr2 /U2 . Wiesner [94] a rassemblé l’ensemble des modèles connus en 1967
et a proposé la formulation du facteur de glissement suivante :
q
cosβ′2 /Z0,7
σ = 1− (4.20)
1 − ϕ2 tanβ′2
avec Z le nombre d’aubes. Cette formulation est utilisable lorsque la relation suivante
est observée :
−8, 16 × cosβ′2
D1
< exp (4.21)
D2 Z
Lorsque ceci n’est pas observé, l’auteur propose une correction à apporter à cette
formulation. D’après l’auteur, la formulation (4.20) prévoit des valeurs de facteur de
112 Principe de fonctionnement des pompes rotodynamique
glissement proches de celles prévues par les autres modèles et a l’avantage d’être
simple d’utilisation. Néanmoins le modèle de Wiesner n’évalue pas correctement le
facteur de glissement pour toutes les géométries de rotor [71]. Une autre formulation
[72] très utilisée est la suivante :
1
σ= (4.22)
Km 1 + cosβ′2
1+
Z 1 − ( D1 /D2 )2
avec Km = 0, 02 × Ns + 0, 94 où Ns est la vitesse spécifique de la pompe (Cf. section
4.6).
Cependant les modèles les plus anciens ne sont fonction que de paramètres géomé-
triques (i.e. β′2 , Z, D1 /D2 , etc.) et aucunement de paramètres liés à l’écoulement (par
exemple le coefficient ϕ2 comme cela est le cas dans le modèle de Wiesner). De ce fait
ces modèles peuvent être inadaptés sur des points de fonctionnement à rendement
plus faible (i.e. off-design). Qiu et al. [71] propose en 2011 alors une nouvelle formu-
lation du facteur de glissement pouvant être utilisée lorsque le rendement n’est pas
maximal et pour lequel les écoulements dans le rotor n’est pas bien conditionné. Ce
facteur s’écrit de façon général :
où σrad. représente les effets de la rotation , σturn. représente les effets de l’inclinaison
des aubes et σpassage représente la variation de la surface de passage. Le deuxième
terme semble être prépondérant dans les écoulements off-design. Dans son étude, Qiu
teste son modèle sur 90 rotors de formes différentes avec 1650 débits différents. Ce
modèle semble évaluer des valeurs numériques du facteur de glissement très proches
de celles mesurées, ceci pour des pompes de rotor de type très différent. Il est à noter
cependant que ce modèle nécessite des valeurs numériques de paramètres parfois
difficiles à évaluer. Ceci est le cas par exemple pour la variation radiale de l’angle
d’inclinaison, i.e. dβ/dm avec m la distance méridienne, alors que ce paramètre semble
être prépondérant dans une bonne évaluation du facteur de glissement.
Une fois le facteur de glissement estimé, on peut alors s’intéresser à la charge
réelle délivrée par le rotor, Hm , différente de celle prévue par la théorie d’Euler, HE . La
première est plus faible que cette dernière, on peut définir un rendement hydraulique,
ηh , caractérisant les dissipations énergétiques présentes dans le rotor uniquement :
Hm gHm
ηh = = (4.23)
HE U2 Vu2
En utilisant le triangle des vitesses, il vient :
ηh σU22
ϕ2
Hm = × 1+
g tan β′2
Ainsi, d’après cette formulation, dès lors que le rendement hydraulique est connu et
que le facteur de glissement est évalué à l’aide de la formulation (4.20), il est possible
de prévoir la charge réelle délivrée par le rotor à l’écoulement.
Analyse d’une pompe axiale 113
caractérisées par une vitesse spécifique élevée (i.e. NS > 200 ou ωS > 4). En effet,
elles sont principalement utilisées dans le cadre de forts débits (103 à 105 m3 /h) pour
des faibles hauteurs manométriques (Hm ∼ 10 m). Cette plage des débits est aussi
couverte par les pompes centrifuges, mais la hauteur manométrique associée est 10
voire 100 fois plus élevée. Nous supposerons par ailleurs que les vitesses à l’entrée
à la pompe sont exemptes de toute composante giratoire. Les pompes axiales ont
un certain nombre d’avantages [83] : (i) des dimensions transversales réduites par
rapport aux pompes centrifuges, (ii) une modification possible et lors de l’utilisation
de l’inclinaison des aubes. Par contre l’élévation en hauteur reste faible (comme cela
a été dit plus haut) et une courbe de rendement très ascendante lorsqu’on s’éloigne
du point de meileur rendement.
4.12.1 Constitution
Une pompe axiale est composée d’un convergent d’entrée, d’une roue à aubes
mobile, d’un redresseur et d’un diffuseur axial. Le convergent permet d’uniformiser
les vitesses moyennes à l’entrée tout en diminuant le taux de turbulence. Le redresseur
est principalement destiné à ramener les filets fluides suivant l’axe de la pompe en
permettant une récupération et une transformation partielle en pression de l’énergie
cinétique. Différentes coupes sont dessinées sur les figures 4.50 et 4.51.
Nous supposerons dans notre cas que les particules de fluide entrant dans la partie
mobile de la pompe axiale conserveront leur distance radiale à l’axe de rotation. Ainsi
la composante radiale de la vitesse de ces particules est nulle. Il s’agit là de l’hypothèse
d’écoulement en équilibre radial traduisant l’équilibre entre les efforts centrifuges et les
efforts de pression. La composante radiale pourrait néanmoins être non-nulle dans le
cas de faibles débits pour lesquels un phénomène de recirculation existe.
L’écoulement généré par les aubes mobiles fait apparaı̂tre une vitesse d’entraı̂nement
−
→ −
→
( U ), et une vitesse relative ( W ), toutes deux reliées par la composition des vitesses
(relation 4.10).
Comme nous l’avons vu précédemment, il est très utile de décomposer la vitesse
d’une particule fluide en une composante débitante ou axiale (Va ), et une composante
Analyse d’une pompe axiale 115
tangentielle (Vu ) :
Va = V sin α = W sin β : Composante débitante ;
Vu = V cos α = U + W cos β : Composante de transfert d’énergie.
116 Principe de fonctionnement des pompes rotodynamique
(a) (b)
Figure 4.51 : Vues frontale (a) et de dessus (b) des aubes mobiles d’une pompe axiale - définition des
vitesses.
La conservation du débit met en évidence le fait que la vitesse axiale est un invariant.
Si le fluide à l’entrée n’est pas animé de pré-rotation, l’évolution du triangle des
vitesses est alors donné par la figure 4.52.
On note bien le fait que la vitesse axiale reste constante, ainsi que la vitesse d’en-
traı̂nement (si on demeure sur le cylindre de rayon constant). On remarque d’autre
−
→
part que les variations des projections sur U des vitesses absolues et relatives entre
la sortie et l’entrée sont égales,
Comme plus haut les efforts de pesanteur seront négligés. L’expression précédente
devient après développement 5 :
−
→ −
→ − → → −→ − → → −
→ − → →
r ∧ ρ W 1 W 1 .− n 1 dS1 + − →r ∧ ρ W 2 W 2 .− n 2 dS2 + −
→
r ∧ ρ W p W p .−
n p dS p
= − →r ∧ − p1 −
→n 1 dS1 + − →
r ∧ − p2 −
→
n 2 dS2
+− →r ∧ − p−→n p dS1 − −→r ∧−→
γe dm
−
→ −
→
− r ∧ γ dm (4.25)
c
L’hypothèse d’un fluide parfait ayant été posée, le fluide glisse le long des parois sans
adhérence avec un vecteur−vitesse perpendiculaire à la normale locale. Le troisième
5. On remarque dans cette expression l’absence des contributions des surfaces de normales
extérieures purement radiale. Ces contributions étant nulles, on a choisi ici de les ignorer.
118 Principe de fonctionnement des pompes rotodynamique
terme du membre de gauche est donc nul. On introduit le débit élémentaire “dans le
sens de l’écoulement”
−
→ − →
dQ = W i . N i dSi
−
→ −
→
avec N = −−
1
→n et N = −
1
→
n . Les produits vectoriels impliquant les efforts d’iner-
2 2
tie sont nuls. Enfin le troisième terme du membre de droite correspond au couple
élémentaire exercé par les aubes sur le fluide dC et orienté selon −
→
z . La relation (4.25)
devient donc
−
→ −
→ −
→ −
→ −
→
ρdQ − →r ∧W −− →
2 r ∧ W = − p r j dS + p r j dS + dC
1 1 1 2 2
UVu2
HE = (4.26)
g
la charge théorique délivrée (i.e., la charge d’Euler ) à un filet de fluide situé à une dis-
tance r de l’axe de rotation. On retrouve une expression similaire à celle obtenue dans
le cas des pompes centrifuges. Ceci est logique car les deux types de pompes font parti
de la même famille des pompes rotodynamiques. Il n’est pas souhaitable cependant
que la charge délivrée par les aubes dépende de la distance radiale. En effet les gra-
dients de pression dans la direction radiale seraient à l’origine d’une déformation des
lignes de courant et donc de pertes hydrauliques. Pour éviter cela, la charge doit être
identique quelle que soit la distance radiale du filet de fluide. Il faut donc Vu2 = K/r,
étant donné le fait que U = rω. Les aubes sont dessinées la plupart du temps pour
satisfaire cette condition.
(a) (b)
Figure 4.53 : Modification du point de fonctionnement à partir (a) d’un vannage sur la conduite de refou-
lement, (b) d’un changement de la vitesse de rotation de la pompe.
120 Principe de fonctionnement des pompes rotodynamique
Réseau sous-estimé
Ce cas n’est pas fréquent car les pertes de charge sont généralement surestimées
dans le calcul. Ce problème provient généralement d’une hauteur géométrique ou
des pertes de charge réelles supérieures à celles estimées. Ce dernier point peut être
consécutif à un encrassement des conduites.
On observe alors sur place un débit plus faible que celui souhaité (Fig. 4.55.a).
Une possibilité consiste alors à réduire les pertes de charge (1) si cela est possible. Si
cette modification n’est pas possible, la vitesse de rotation peut être augmentée (2)
pour déplacer le point de fonctionnement vers des débits plus élevés. Si aucune de
ces modifications ne peut être apportée, le responsable de l’installation devra alors
changer la pompe et adapter ainsi le diamètre des aubes.
Réseau surestimé
Ce cas est plus fréquent que le précédent, et est lié à une hauteur géométrique ou
des pertes de charge réelles inférieures à celles estimées.
Le débit observé sur place est dans ce cas plus élevé que celui souhaité (Fig. 4.55.b).
Il y a un risque de cavitation si la pompe est placée en aspiration. Les pertes de charge
peuvent être augmentées (1) par l’utilisation d’une vanne ou d’un diaphragme. On
peut d’autre part rogner les aubes ou diminuer la vitesse de rotation (2), ces deux
opérations ayant le même effet à savoir une diminution de la charge fournie par la
pompe. Le point de fonctionnement voit alors sa valeur de débit chuter.
122 Principe de fonctionnement des pompes rotodynamique
(a) (b)
Figure 4.55 : Modification du point de fonctionnement lorsque le réseau est mal-dimensionné. a, réseau
sous-estimé ; b, réseau surestimé (illustrations d’après la référence [43]).
Chapitre 5
Analyse Dimensionnelle
L’analyse dimensionnelle est fondée sur le principe fondamental très simple se-
lon lequel une relation entre grandeurs physiques doit être dimensionnellement ho-
mogène, c’est à dire indépendante du système d’unités. Tout particulièrement en
mécanique des fluides, l’analyse dimensionnelle est une étape importante car elle per-
met d’exprimer des lois entre plusieurs grandeurs physiques. Son utilité se fait sentir
dès lors que l’on veut simplifier les termes d’une équation afin d’en alléger l’écriture.
Elle permet aussi d’aboutir à une réduction du nombre de variables limitant ainsi
l’étude paramétrique associée. Nous verrons d’autre part que l’analyse dimension-
nelle guide l’ingénieur dans l’élaboration de tests sur maquettes en lui indiquant les
produits adimensionnels à respecter.
∂ui
=0
∂xi
!
∂ui ∂u ∂p ∂2 u i
ρ + uj i = ρgi − +µ
∂t ∂x j ∂xi ∂x j ∂x j
On considère le problème comme étant caractérisé par différents paramètres (i) une
longueur L, (ii) une vitesse U, (iii) une accélération g, et (iv) une pression ρU 2 . La
vitesse mentionnée ici peut par exemple correspondre à celle observée dans une
124 Analyse Dimensionnelle
région très amont d’un obstacle à étudier. La longueur doit quant à elle caractériser
cet obstacle (diamètre d’une conduite, longueur d’une plaque, etc.). Les variables de
l’écoulement peuvent ainsi être exprimées à partir de ces paramètres de référence et
de variables sans dimension (notées ici avec le symbole + ) : xi+ = xi /L, ui+ = ui /U,
t+ = tU/L, g+ = gi /g, et p+ = p/ρU 2 . Les deux équations précédentes peuvent donc
être manipulées pour faire apparaı̂tre ces variables sans dimension :
U ∂ui+
× + =0
L ∂xi
U 2 ∂ui+ U2 + ∂ui
+
+ U 2 ∂p+ U ∂ui+
ρ × + +ρ × uj = ρg × gi − ρ × + +µ 2 × + +
L ∂t L ∂xi+ L ∂xi L ∂x j ∂x j
On définit les produits sans dimension suivant Fr2 = U 2 /gL et Re = ρUL/µ. Ces
groupes sans dimension sont développés plus bas. Les équations de mouvement adi-
mensionnées s’écrivent finalement :
∂ui+
=0 (5.1)
∂xi+
∂ui+ + ∂ui
+
1 + ∂p+ 1 ∂ui+
+ uj = 2 gi − + + (5.2)
∂t+ ∂xi+ Fr ∂xi Re ∂x + +
j ∂x j
ce qui leur donne alors un sens physique alors qu’ils sont dépourvu pourtant de
dimension. D’autre part les particularités qui affectent un écoulement donné, loca-
lement ou globalement, ne sont pas définies par la valeur de grandeurs mesurables
considérées isolément, vitesse, longueur. . .Ces grandeurs n’ont d’intérêt que dans la
mesure où on considère leurs valeurs dans leur ensemble par l’utilisation des produits
sans dimension. On se propose ici d’énumérer quelques produits sans dimension ob-
servés régulièrement.
ρU 2 /L ρUL
2
= = Re (5.3)
µU/L µ
Ainsi lorsque les efforts d’inertie dominent très fortement les efforts visqueux, le
nombre de Reynolds caractéristique de l’écoulement aura une valeur très importante
(Re >> 1), et l’écoulement sera dans ce cas caractérisé par des mouvements très
désordonnés. On parlera de régime turbulent. Au contraire les écoulements à nombre
de Reynolds très faibles auront des trajectoires déterministes et seront dit laminaires.
dans la région amont absente du phénomène transitoire, ( L/U )/T. On préfère généra-
lement manipuler la fréquence de présence du phénomène transitoire f = 1/T (par
exemple la fréquence des lâchers de tourbillons d’un cylindre totalement immergé
dans un écoulement permanent), et le nombre de Strouhal s’écrit finalement :
fL
St = (5.5)
U
Le temps caractéristique d’un écoulement permanent est une valeur très élevé et ten-
dant mathématiquement vers l’infini. La fréquence devient donc nulle, entrainant une
valeur nulle du nombre de Strouhal. Le concept de nombre de Strouhal est repris plus
loin, page 157, dans le chapitre concernant l’aérodynamique.
p/L p
2
= (5.6)
ρU /L ρU
Les efforts de pression sont généralement importants dans le cas d’écoulements d’air
à grande vitesse, la relation des gaz parfaits (i.e. p/ρ = rT = c2 ) permet alors de
transformer la relation précédente :
Efforts de pression
= γM2 (5.7)
Efforts d’inertie
avec M = U/c le nombre de Mach. Il est couramment admis que lorsque sa valeur
devient plus élevée que 0,2 − 0,3 les effets de compressibilité deviennent importants.
On note au passage que le fait d’introduire la célérité de l’air dans le problème en-
traine irrémédiablement la nécessité d’une équation de l’énergie pour l’estimation du
champ de température dont dépend la célérité.
5.3.1 Principe
La mise en forme adimensionnelle fait appel au théorème π mis en place par
Vaschy [87] et Buckingham [24]. Considérons un phénomène physique décrit par une
relation de la forme
f ( x1 , x2 , . . . , x n ) = 0 (5.8)
α
π1 = x1α1 x2α2 . . . xk k xk+1
β β β
π 2 = x 1 1 x 2 2 . . . x k k x k +2
..
.
η η η
π n − k = x1 1 x2 2 . . . x k k x n
— Le théorème π suppose que l’on est capable de lister toutes les grandeurs
physiques qui interviennent dans le phénomène physique ;
— Les paramètres xi ne sont pas obligatoirement une entité seule mais
peuvent être un groupement de paramètres (par exemple g∆ρ pour mettre
en évidence la poussée d’Archimède dans un fluide) ;
— Le théorème π nous indique le nombre de produits sans dimension que
l’on peut construire, il n’y a pas en général unicité de la solution. L’uti-
αj
lisateur fait le choix des paramètres xi en cherchant si possible à faire
apparaı̂tre des produits sans dimension habituels (Re, Fr, Nu, Pr, M,...) ;
— Il faut aussi comprendre que le théorème π ne donne pas la forme de la
relation qui lie les différents πi , la recherche de cette relation fait l’objet
du dépouillement des résultats expérimentaux.
Il vient donc (α1 , α2 , α3 ) = (−2, −1, −2). D’où le premier nombre sans dimension,
π1 = T/(ρd2 U 2 ). On fait de même pour le deuxième nombre sans dimension :
β2 + 1 = 0 : exposants de M ;
π2 β − 3β 2 + β 3 − 1 = 0 : exposants de L ;
1
− β3 − 1 = 0 : exposants de T.
Similitude 129
d’où ( β 1 , β 2 , β 3 ) = (−1, −1, −1) entrainant π2 = µ/(ρdU ). Il s’agit là (de l’inverse)
du nombre de Reynolds. Le théorème π indique donc que le produit sans dimension
impliquant la traı̂née est lié au nombre de Reynolds :
T
= f (Re)
ρd2 U 2
On fait généralement apparaı̂tre ici la pression dynamique (i.e. 1/2ρU 2 ) et une sur-
face de référence A dans ce produit correspondant alors l’écriture du coefficient de
traı̂née, CD = T/(1/2ρU 2 A). Dans le cas d’une sphère la surface de référence est le
maı̂tre-couple A = πd2 /4 c’est à dire la surface projetée de la sphère sur un plan
perpendiculaire à l’écoulement. A posteriori l’expérience montre que l’évolution du
coefficient de traı̂née avec le nombre de Reynolds est similaire pour des sphères de
diamètres différents (Fig. 5.1). Nous notons qu’au-delà d’un nombre de Reynolds
égal à 1000, celui-ci n’a quasiment plus aucune influence sur le coefficient de traı̂née.
On observe tout de même une chute de ce dernier pour des nombres de Reynolds
proches de 2 × 105 pour redevenir stable. Cette chute est consécutive à la présence
d’une couche limite turbulente autour de la sphère.
On a fait l’hypothèse dans l’exemple précédent d’une surface lisse, ce qui n’est
pas toujours le cas. En effet la rugosité peut avoir une influence de l’écoulement
dans la région proche de la paroi et de ce fait influencer l’évolution du coefficient
de traı̂née avec le nombre de Reynolds. Ceci est illustré sur la figure 6.19.
5.4 Similitude
Un champ d’application de l’analyse dimensionnelle est celui des problèmes de si-
militude entre un modèle (i.e., maquette à échelle réduite) et un prototype. Cette tech-
130 Analyse Dimensionnelle
nique est couramment utilisée lorsque les essais réels sont très complexes et lorsque
les simulations numériques deviennent trop coûteuses en temps de calcul.
L’analyse dimensionnelle est là essentiellement pour donner des informations sur
les conditions de similitude pour lesquelles les mesures effectuées à échelle réduite
sont représentatives du système réel. En effet, il serait trivial de penser que tous les
phénomènes mesurés sur maquette sont proportionnels à ceux observés en réalité.
De plus, il peut y avoir des différences notables du fait que les essais sont effectués à
échelle réduite, cet effet est appelé effet d’échelle.
Malgré cela, il existe des domaines pour lesquels les essais sur modèles sont
nécessaires : l’hydraulique (études de barrages, de canaux, d’aménagements côtiers),
l’hydrodynamique (effet de la houle sur les structures marines), l’aérospatial (concep-
tion de maquettes d’avions, de fusées).
Voyons maintenant ce qui caractérise les études de similitudes.
φ ( π1 , π2 , . . . , π n − k ) = 0
La similitude sera donc respectée si tous les nombres sans dimension sont identiques
pour la maquette et pour le prototype.
Tableau 5.1 : Prédiction numérique du couple de frottement d’un disque de diamètre 300 mm en rotation
autour de son axe de révolution dans l’air.
ω (rad/s) 300 500 700
C (N.m) 4,67 ×10−2 1,20 ×10−1 2,21 ×10−1
C ρR2 ω
Cm = et Re =
1/2ρR5 ω 2 µ
Le coefficient de couple de frottement est estimé pour les trois vitesses (Tab. 5.2), et
son évolution avec le nombre de Reynolds peut être approchée par une loi de type
puissance :
0, 0939
Cm =
Re0,161
L’estimation du couple de frottement exercé par l’huile sur un disque de diamètre
égal à 120 mm est obtenue à partir de l’estimation du nombre de Reynolds dans
cette nouvelle configuration, puis du couple de frottement adimensionnel, et enfin le
couple de frottement. Les résultats obtenus à chacune de ces étapes sont rassemblés
dans le tableau 5.3. On remarque que les nombres de Reynolds atteints dans le cas du
disque en rotation dans l’huile (Tab. 5.2) sont dans la gamme de ceux observés dans
la configuration numérique (Tab. 5.3) nous permettant ainsi d’utiliser l’approximation
en puissance.
Tableau 5.2 : Coefficient de couple de frottement d’un disque de rayon 150 mm en rotation autour de son
axe de révolution dans l’air.
ω (rad/s) 300 500 700
Re = ρR2 ω/µ 4,37 ×105 7,28 ×105 1,02 ×106
Cm = C/(1/2ρR5 ω 2 ) 1,15 ×10−2 1,01 ×10−1 0,95 ×10−2
Tableau 5.3 : Couple de frottement exercé par l’huile sur un disque de rayon 120 mm en rotation autour
de son axe de révolution.
ω (rad/s) 300 500 700
Re = ρR2 ω/µ 2,16 ×105 3,60 ×105 5,04 ×105
Cm = 0,0939/Re0,161 1,30 ×10−2 1,20 ×10−1 1,13 ×10−2
C = Cm × (1/2ρR5 ω 2 ) (N.m) 13,40 34,34 63,40
132 Analyse Dimensionnelle
5.5 Récapitulatif
Afin d’y voir plus clair, récapitulons un peu ce que nous venons de voir. Il apparaı̂t
donc plusieurs types d’écoulements pour lesquels certains processus physiques do-
minent. Et lorsque l’on désire effectuer une étude sur maquette, l’ensemble des pro-
duits sans dimension peuvent être ou tous respectés (similitude complète), ou alors
la théorie rencontre une incohérence. Dans cette dernière configuration, il importe
alors de respecter en priorité certains nombres sans dimension, caractéristiques du
phénomène physique dominant.
1. L’écoulement doit alors présenter une turbulence fortement dévelopée, entraı̂nant dès lors un
problème pour les écoulements surfaciques de faible profondeur. Ce problème peut néanmoins être
contourné en introduisant une distorsion d’échelle.
Chapitre 6
Introduction à l’aérodynamique
Figure 6.1 : Soufflerie de l’ONERA à Modane dans les Alpes françaises (d’après [89]).
µ R a
π1 = ≡ Re π2 = π3 = ≡M
ρLV ρL2 V 2 V
1. On pourrait impliquer la pression, en lieu et place de la célérité, ce qui conduirait à un résultat
identique en utilisant la loi des gaz parfaits et la définition de la célérité.
La résultante aérodynamique 135
R
D’après les propriétés de la variationnelle φ(π1 , π2 , π3 ) = 0, entraı̂nant = f(Re,
ρL2 V 2
M). Rappelons que cette analyse a été effectuée pour une incidence donnée. Si celle-ci
est amenée à varier, la résultante sera aussi changée. On aura finalement en introdui-
sant un coefficient C :
R
= C (Re, M, i) (6.1)
ρL2 V 2
avec C dépendant des nombres de Reynolds et de Mach, ainsi que de l’incidence.
Nous verrons plus loin que ce coefficient dépend aussi et surtout de la forme de
l’obstacle. On remarque qu’au dénominateur, la dimension de longueur est à la puis-
sance 2 : on choisira donc une surface de référence S au lieu d’une longueur L. Une
discussion quant au choix de cette surface sera ménée au paragraphe 6.5 (p. 155). On
écrira alors
1
R = CρSV 2 (6.2)
2
La présence du coefficient 1/2 permet de mettre en évidence la pression dynamique,
1/2ρV 2 .
On définit un système d’axes orthonormés (− →
x, −
→z ) de telle sorte que −
→x soit
−
→
aligné selon l’écoulement et z perpendiculaire à l’écoulement. La décomposition
de la résultante des efforts dans ce système d’axes donne respectivement la traı̂née T
et la portance P (Fig. 6.2). On peut dès lors exprimer ces deux efforts identiquement
à la relation (6.2) :
1
T = Cx ρSV 2 (6.3)
2
1
P = Cz ρSV 2 (6.4)
2
Le fait d’écrire les efforts de traı̂née et de portance sous cette forme permet d’autre
part d’isoler les termes pour lesquels on connait bien l’influence (ρ, S et V). Par contre
ceux qui ont une influence plus subtile (inclinaison, l’état de surface, . . .) sont inclus
dans les coefficients mis en place. Ces coefficients, Cx et Cz , sont appelés respective-
ment coefficients de traı̂née et de portance, avec
(a) (b)
Figure 6.5 : Efforts de pression sur l’intrados et l’extrados. a, principe ; b, mesures effectuées par le Labo-
ratoire d’Études Aérodynamiques (Poitiers) sur une aile inclinée.
il existe un point pour lequel la vitesse est perpendiculaire à la paroi : le point d’arrêt,
A. En ce point, la vitesse est minimale, et la pression maximale (p∞ + 1/2ρV 2 , avec
p∞ la pression statique à l’infini). À partir de ce point, la vitesse augmente et la pres-
sion diminue. Il détermine donc le point de séparation des écoulements intrados et
extrados. En ce qui concerne les deux surfaces : (i) sur la surface inférieure (intrados)
la vitesse croı̂t et la pression décroı̂t jusqu’à ce que chacune retrouve leur valeur à l’in-
fini amont. La surpression est donc située entre p∞ + 1/2ρV 2 et p∞ ; (2) sur la surface
extérieure (extrados) la vitesse croı̂t jusqu’à une valeur supérieure à celle des condi-
tions à l’amont infini, de ce fait la pression diminue d’autant plus. Il existe alors une
forte zone d’aspiration. On remarque donc que l’écoulement autour du profil tend
à aspirer la surface supérieure et à pousser la paroi inférieure. La portance est ainsi
créée. Notons tout de même, que les efforts intrados et extrados n’ont pas la même
intensité. L’effet aspirant est plus important que l’effet de poussée et représente 2/3
voire 4/5 de la portance totale.
138 Introduction à l’aérodynamique
Néanmoins, il ne faut pas poser de suite qu’il est préférable d’avoir un écoulement
laminaire dans la couche limite. La traı̂née de frottement n’est pas la seule traı̂née
présente, et la traı̂née de pression qui dépend fortement de la forme de l’obstacle
dans l’écoulement peut elle devenir très importante, même dans le cas laminaire.
2. Notons que le terme de couche limite est artificiellement introduit. La transition de la région
proche de la paroi vers l’écoulement externe n’est pas brutale mais prend place de façon continue. Le
terme “limite” n’est donc a priori pas adapté.
Origine des forces de résistance (traı̂née) 139
(a) (b)
Figure 6.8 : a, décollement de la couche limite (Photo ONERA) ; b, distribution de la pression autour d’un
obstacle selon le type d’écoulement (attaché ou décollé).
l’aval (Fig. 6.9). Il en résulte l’ajout d’une résistance à l’avancement et une diminution
de la portance consécutivement à une longueur efficace de l’aile plus faible. Il est à
noter d’autre part que ces tourbillons sont très persistants, et le fait qu’un avion situé
en aval rencontre ces perturbations peut être très dangereux. Face à des avions deve-
nant de plus en plus lourds, augmenter la portance devient nécessaire. La longueur
des ailes peut être augmentée mais ceci n’est souvent pas possible pour des raisons
d’encombrement sur les pistes d’aéroport 3 . Une solution proposée par Richard T.
Whitcomb et son équipe du centre de recherche de la NASA (1976) consiste alors à
amoindrir les tourbillons marginaux par l’utilisation d’ailettes marginales (Fig. 6.10).
Ceci a pour autre conséquence de générer une portance latérale légérement dirigée
vers l’avant, diminuant alors la traı̂née de l’aile. Plus récemment, un passionné d’avia-
tion, Christian Hugues [60], a conçu un dispositif appelé minix (Fig. 6.11). Ce der-
nier en forme de cylindre à fente hélicoı̈dale “fluidifie” l’écoulement ayant pour une
conséquence à la fois une diminution de la traı̂née proche de 8% pour les écoulements
incompressibles et une augmentation de la portance proche de 5,5%.
Cette traı̂née porte mal son nom, car même si elle est bien induite par la portance
(c’est le prix à payer pour avoir de la portance), elle ne réagit pas comme une traı̂née
au sens 1/2CρV 2 . Plus une aile porte, plus la traı̂née induite est importante. Si une
aile ne porte pas, comme cela peut être le cas lors d’une chute libre, celle-ci n’existe
pas.
(a) (b)
Figure 6.10 : a, ailette marginale en bout d’aile d’un A380 (Airbus S.A.S. 2008 - Photo exm company / H. Goussé) ;
b, ailette marginale testée en soufflerie (d’après Whitcomb [91]).
tement connue ; (2) la couche limite fait apparaı̂tre un décollement, le sillage devient
alors important. La résistance totale est très importante, et la résistance de frottement
devient négligeable devant celle de pression. Cette résistance de pression (en présence
de décollement) est appelée Traı̂née de forme. Dans le cas d’un cycliste par exemple
pour lequel celui-ci essaye d’adopter la position la plus aérodynamique possible Mar-
tin et al. montrent que la traı̂née de pression contribue à 90% de la pression totale [57].
Il devient alors judicieux pour le cyliste d’adapter au mieux sa position afin que la
pression le pénalise le moins possible.
4. Ludwig Prandtl, mécanicien allemand (1875 - 1953). Prandtl est le père fondateur de
l’aérodynamique moderne en lui donnant de nombreuses théories. Il fût d’autre part le professeur
de Theodore Von Karman (1881 - 1963) célèbre chercheur hongrois en aéronautique et de Paul Richard
Heinrich Blasius (1873 - 1970) mécanicien allemand qui donna une première application à la présence
de la couche limite sur une plaque plane (cf. paragraphe 6.4.5 page 147).
Étude simple de la couche limite 143
fluctuations de vitesses. L’effort de frottement devient alors plus grand. Nous verrons
néanmoins qu’il est préférable d’avoir une couche limite turbulente.
Une définition qualitative vient d’être donnée, or on définit souvent l’épaisseur de
couche limite δ comme étant la distance à la paroi à laquelle u est très voisin de la
vitesse de l’écoulement extérieur Ue , plus précisemment :
u(y = δ) = 0, 99 × Ue
La couche limite est donc la région de l’écoulement proche d’une paroi caractérisée
par des vitesses au moins inférieures à 1% de celle de l’écoulement général.
Cette frontière demeure néanmoins très difficile à estimer expérimentalement car
le gradient de vitesse y est très faible. De plus, elle est fluctuante dans le temps
consécutivement à des effets d’intermittence et non pas statique.
∂2 u ∂2 u
∂u ∂u 1 ∂p
u +v =− +ν + 2 (6.6)
∂x ∂y ρ ∂x ∂x2 ∂y
∂2 v ∂2 v
∂v ∂v 1 ∂p
u +v =− +ν + (6.7)
∂x ∂y ρ ∂y ∂x2 ∂y2
En raison de la faible épaisseur de la couche limite, le mouvement du fluide dans
la couche limite est principalement parallèle à l’élément plan. L’évolution de la vi-
tesse selon y est très grande (ces vitesses évoluent d’une valeur nulle à une vitesse
égale à celle de l’écoulement externe à la couche limite sur une distance très courte,
l’épaisseur de couche limite). Suivant la direction x, la vitesse évolue très faiblement,
une variation notable est observée sur une longueur caractéristique du problème
(par exemple, les dimensions du corps). Pour toutes ces raisons les dérivées dans
la direction transversales sont prépondérantes devant celles selon l’écoulement, ainsi
∂2 /∂x2 ≪ ∂2 /∂y2 .
D’autre part les ordres de grandeur des termes du membre de gauche de l’équation
projetée selon x valent
2
∂u ∂u U
u = v ∼
∂x ∂y x
et ceux de l’équation projetée selon y :
2
∂v ∂v U δ
u = v ∼
∂x ∂y x2
De ce fait, pour des points localisés dans la couche limite et placés dans une région
éloignée du point d’arrêt l’épaisseur de la couche limite reste très négligeable devant
la valeur de la localisation, x. Les termes de convection de l’équation projetée selon x
sont très importants devant ceux de l’équation projetée selon y.
Enfin, l’équation de continuité permet d’atteindre l’ordre de grandeur de la vitesse
longitudinale, [v] ∼ [Uδ/x ]. Les termes de viscosité de l’équation projetée selon x sont
dès lors prépondérants devant ceux de l’équation projetée selon y. Cette dernière ne
fait plus apparaı̂tre qu’un seul terme, le gradient de pression.
Les équation de la couche limite s’écrivent finalement
∂u ∂u 1 dp ∂2 u
u +v =− +ν 2 (6.8)
∂x ∂y ρ dx ∂y
∂p
=0 (6.9)
∂y
La relation (6.9) tend à dire que la pression est constante le long de l’épaisseur
de la couche limite. De ce fait, la pression qui est appliquée sur cette couche par
l’écoulement est intégralement restituée à l’obstacle, et peut être estimée à partir de
l’équation de Bernoulli, p + 1/2ρU 2 = Cste. La pression de l’écoulement externe
à la couche limite joue alors un rôle primordiale dans la formation de celle-ci. Par
exemple, le point de transition laminaire/turbulent dépend fortement de cela. Si la
Étude simple de la couche limite 145
pression décroı̂t fortement dans le sens de l’écoulement, ce qui se passe à l’arrière des
obstacles, il est alors possible que la couche limite se sépare (cf. § 6.4.3).
Les équations (6.8) et (6.9) représentent les équations de Prandtl auxquelles on
rajoute l’équation de conservation de la masse. Il faut d’autre part indiquer les condi-
tions aux limites : (1) en paroi, y = 0 : u = 0 ; v = 0, (2) hors de la couche limite,
y −→ ∞ : u = Ue ; v = 0.
∂2 u ∂2 u
1 dp dp
0=− +ν ou =µ (6.10)
ρ dx ∂y2 y =0 dx ∂y2 y =0
La figure 6.13 propose des formes de profil de vitesse pour différents signes du
gradient de pression, dp/dx. On note que dans la région caractérisée par dp/dx > 0,
il existe une zone de retour de l’écoulement (appelée zone de recirculation). On dit
alors qu’en ce point, la couche limite se sépare ou “décolle”. La tangente au profil de
vitesse y est normale à la surface de la paroi, (∂u/∂y)y=0 = 0. La contrainte pariétale,
τp , devient donc nulle en ce point :
∂u
τp = µ =0 (Décollement).
∂y y =0
∂p ∂U
= −U (6.12)
∂x ∂x
Ce gradient de pression adverse étant valable aussi dans la couche limite (d’après
l’expression 6.9), s’il est suffisamment important, il peut renverser l’écoulement et
provoquer un décollement de la couche limite.
5. Cet effet se rajoute à l’épaississement de la couche limite provoqué par la diffusion de quantité
de mouvement consécutive à la viscosité.
Étude simple de la couche limite 147
∂u ∂v
+ =0 (6.13)
∂x ∂y
∂u ∂u ∂2 u
u +v =ν 2 (6.14)
∂x ∂y ∂y
∂p
=0 (6.15)
∂y
L’expérience montre une évolution du profil des vitesses u(y) le long de l’obstacle
(i.e., en fonction de x). On peut cependant écrire de façon identique ces profils en
faisant√intervenir une nouvelle variable η qui fait intervenir les deux autres x et y,
η = y U/νx.
Dans ce cas, le profil de vitesse u(η ) reste inchangé (Fig. 6.14). On définit alors une
fonction de courant ψ uniquement fonction de la nouvelle variable :
√ Ux
ψ( x, y) = νxU f (η ) = √ f (η )
Rex
6. Paul Richard Heinrich Blasius (1873-1970), mécanicien allemand. La relation (6.17) a été intro-
duite par Blasius lors de son travail de thèse en 1908. Il fut un des élèves de Prandtl, et cette rélation
représente la première application de la théorie de la couche limite énoncé par Prandtl en 1904.
148 Introduction à l’aérodynamique
2 f ′′′ + f f ′′ = 0 (6.17)
Cette équation est appelée équation de Blasius. Remarquons comme il est surprenant
à partir de cette méthode transformant les variables de départ (x et y) d’obtenir une
équation différentielle très simple dont les solutions sont données dans le tableau 6.1
et sur la figure 6.15.
η f f ′ = u/U f ′′ η f f ′ = u/U f ′′
0 0 0 0,33206 2,6 1,07252 0,77246 0,20646
0,2 0,00664 0,06641 0,33199 3,0 1,39682 0,84605 0,16136
0,4 0,02656 0,13277 0,33147 4,0 2,30576 0,95552 0,06424
0,6 0,05974 0,19894 0,33008 5,0 3,28329 0,99155 0,01591
0,8 0,10611 0,26471 0,32739 6,0 4,27964 0,99898 0,00240
1,0 0,16557 0,32979 0,32301 7,0 5,27926 0,99992 0,00022
1,4 0,32298 0,45627 0,30787 8,0 6,27923 1,00000 0,00001
1,8 0,52952 0,57477 0,28293 8,8 7,07923 1,00000 0,00000
2,2 0,78120 0,68132 0,24835
√ √
Soit η = δ U/νx = 5, ce qui nous permet d’écrire : δ/x = 5/ Rex . On connaı̂t ainsi
l’évolution de l’épaisseur de couche limite pour une plaque plane avec une incidence
nulle. La solution de l’équation de Blasius faisant intervenir la grandeur f ′ en paroi,
on peut donc formuler la contrainte pariétale, et ensuite le coefficient de frottement
pariétal, respectivement :
ρU 2 τ0 ( x ) 0, 664
τ0 ( x ) = 0, 332 √ et c f (x) = 1
= √ (6.18)
Rex 2 Rex
2 ρU
vers la turbulence est observée pour des valeurs du nombre de Reynolds Re = 3,8
×105 pour la sphère (Fig. 6.18) et Re = 3,0 ×105 pour le cylindre (Fig. 6.3). La taille
du sillage étant fortement diminuée, le coefficient de traı̂née est alors caractérisé par
des valeurs plus faibles. On parlera de crise de la traı̂née.
L’état de surface des obstacles influence la localisation dans la couche limite de
la transition vers la turbulence. Le nombre de Reynolds critique devient plus faible
lorsque la rugosité de la surface augmente [13, 14] (Fig. 6.19). La figure 6.20 propose
l’évolution du coefficient de traı̂née d’une balle de golf, d’un ballon de football, et
d’une sphère lisse. Parmi ces trois obstacles, la chute de la traı̂née se produit pour
une valeur plus faible du Reynolds pour la balle de golf comparativement à la sphère
lisse. La longueur parcourue par cette balle est susceptible d’être plus importante dans
ce cas. De ce fait, les équipementiers sportifs s’intéressent grandement à ces effets afin
de proposer des produits de plus en plus performants. Les fabricants de ballons de
football ont notamment travaillé l’état de surface, par le biais entre autre des coutures,
dans le but d’améliorer leur stabilité et permettant d’autre part de diminuer de près
de 30% la traı̂née (Fig. 6.21).
0, 37x
δ= (6.19)
Re1/5
x
0, 074
cf = (6.20)
Re1/5
x
expression à comparer avec la relation 6.18. Il apparaı̂t donc que les frottements vis-
queux sont plus importants lorsque la couche limite devient turbulente.
Étude simple de la couche limite 151
Modélisation de la turbulence
Les expressions des coefficients de frottement vus plus haut sont basées sur des
configurations simples d’écoulement (plaque plane). Dans des cas de formes d’obs-
tacles plus complexes, celles-ci sont inutilisables. Les codes de calculs numériques,
permettant alors d’obtenir des solutions quelle que soit la configuration, impliquent
les équations de conservation de la masse, le bilan de quantité de mouvement, et
l’équation d’énergie (dans le cas d’écoulements compressibles). Ces calculs sont ef-
fectués à l’aide de techniques numériques (différences finies, éléments finis, volumes
finis,. . .) dont il ne sera pas discuté ici.
Elles sont basées aussi pour certaines sur une modélisation de la turbulence. Dans
ce type d’approche, les effets de la turbulence sont simplement impliqués par l’addi-
tion d’un terme, la viscosité turbulente µ T , à la viscosité “moléculaire” µ. Le bilan de
quantité de mouvement selon x s’écrit alors :
∂2 u
∂u ∂u ∂p
ρ u +v = − + (µ + µ T ) 2 (6.21)
∂x ∂y ∂x ∂y
où les différents paramètres de vitesse et pression sont des valeurs moyennes (la
turbulence est mis de côté pour ne s’intéresser qu’à l’écoulement moyen). La difficulté
réside ici dans l’estimation de la valeur locale de la viscosité turbulente à partir des
propriétés de l’écoulement. Il faut en effet que ce paramètre représente au mieux les
particularités de la turbulence en chacun des points de calculs. Il existe pour celà
plusieurs types de modèles : algébriques pour les plus simples, et avec une ou deux
équations supplémentaires représentant le transport de l’énergie cinétique turbulente
et du taux de dissipation pour les plus complexes. Le plus connu à ce jour est le
modèle kϵ qui donne des résultats satisfaisants pour des écoulements en absence de
recirculation.
profil “laminaire” NACA 63-215 ayant un coefficient de traı̂née plus faible mais ca-
ractérisé par un domaine d’utilisation de portance plus réduit (Fig. 6.23). Plus tard le
profil NASA NLF(I)-0115 [78] fût proposé et possèdant un coefficient de traı̂née aussi
faible que le profil NACA 63-215 tout en ayant une gamme d’utilisation de la portance
aussi grande que le profil NACA 23015.
Par contre, si l’obstacle n’est plus profilé, le sillage est tel que la traı̂née de pres-
sion devient importante et celle de frottement très négligeable. La transition vers la
turbulence est alors une solution pour à la fois réduire la zone de sillage et la traı̂née.
La surface de référence doit rester constante lorsque l’obstacle voit son inclinaison
varier. Elle représente de façon intrinsèque l’obstacle, quelle que soit la configura-
tion. Les clichés d’Étienne-Jules Marey 8 illustrent très bien ce propos (Fig. 6.26). Les
maı̂tres-couples du prisme dans les deux cas sont identiques, mais le sillage est plus
étendu dans la deuxième configuration. La traı̂née sera donc plus importante. On uti-
lisera dans cet exemple la surface mouillée pour les raisons évoquées plus haut dans
ce paragraphe.
Notons à cette occasion, le rôle crucial des arrêtes vives dans l’aérodynamisme des
corps.
Dans le domaine de l’automobile, la quantité SCx est plus utilisée que le coefficient
de traı̂née, celui-ci désignant la qualité d’un véhicule indépendamment de sa taille.
En divisant ce paramètre par la surface frontale du véhicule (comprise entre 1,75 m2
et 2,25 m2 ), on obtient le Cx (de l’ordre de 0,28 et 0,40 pour les véhicules actuels). Pa-
reillement, dans le domaine de construction des grues, cette même quantité (i.e. SCx )
appelée surface au vent est étudiée [88]. Il est très difficile dans cette application d’es-
timer sa valeur numérique puisque l’obstacle est en treillis rendant alors l’écoulement
dans son voisinage très complexe et des effets de sillage très prononcés.
se comporter de façon turbulente bien que la couche limite sur ce cylindre reste lami-
naire. Dans la région sous-critique, 300 ≤ Re ≤ 3,5 × 105 , la couche limite laminaire
se sépare à 80 degré du point d’arrêt, et la production de tourbillons est cette fois
bien définie et périodique. Dans le régime de transition, 1,5 × 105 ≤ Re ≤ 3,5 × 106 ,
la couche limite du cylindre est turbulente et décolle à une valeur plus en aval, 140
degré. Ce régime, appelé crise de la traı̂née, est caractérisé par une chute de la traı̂née.
Dans le régime supercritique, Re ≥ 3,5 × 106 , la génération périodique de tourbillons
reprend.
Les lâchers d’instabilités par l’obstacle sont caractérisés par le nombre de Strouhal,
fD
St = ,
U
où U est la vitesse amont au cylindre, D son diamètre, et f le nombre d’instabilités
générées par une même région du cylindre (i.e. supérieure ou inférieure) pendant une
seconde. Ce nombre correspond donc à la fréquence adimensionnelle des lâchers de
tourbillon par le solide. Les efforts aérodynamiques sont sensibles à ces instationna-
rités et sont caractérisés à leur tour par des valeurs périodiques. Ainsi la portance
évolue autour d’une valeur nulle avec une fréquence identique à celle des lâchers de
tourbillons d’une même rangée. La traı̂née est sensible quant à elle à chaque tour-
billon généré, et donc au double de la fréquence précédente. Le nombre de Strouhal
158 Introduction à l’aérodynamique
dépend du régime d’écoulement (Fig. 6.28), mais aussi de la forme de l’obstacle. Ainsi
le nombre de Strouhal pour des formes prismatiques est affecté par les dimensions
relatives de l’obstacle, les formes des coins et est compris dans la gamme 0,1 < St <
0,2 (Fig. 6.29). Pour des obstacles arrondis, cette gamme devient 0,2 < St < 0,3 (Fig.
6.30).
Figure 6.30 : Nombre de Strouhal pour des obstacles fermés (d’après Bachmann [36]).
160 Introduction à l’aérodynamique
pour la section circulaire. On observe d’autre part que lorsque le rayon du congé aug-
mente, une crise de la traı̂née existe et ceci pour un nombre de Reynolds critique Rec se
rapprochant de celui observé pour la section circulaire (Rec = 3, 0 × 105 ), au fur et à
mesure que le congé voit sa valeur augmenter. Ainsi toutes ces observations montrent
que les valeurs numériques des coefficients dépendent de la forme de l’obstacle plus
que de la taille de cet obstacle.
figure 6.32 que le fait d’incliner le profil augmente la dépression sur l’extrados. Ce-
pendant, cette évolution est rompue pour une inclinaison trop importante de l’aile.
Étude globale
Les variations des deux coefficients vis à vis de l’angle d’attaque sont données sur
la figure 6.33.a. On note que le coefficient de portance croı̂t de façon linéaire avec α
jusqu’à une valeur (appelée angle de décrochage) pour laquelle l’avion voit sa por-
tance terriblement chuter. Ceci est consécutif au décollement de la couche limite, étant
donnée la forte réduction de la vitesse proche du profil. Le coefficient de traı̂née fait
apparaı̂tre quant à lui un minimum. On voit donc que pour des angles différents la
portance est maximale et la traı̂née minimale. Il convient alors de trouver un compro-
mis sur la valeur de l’angle, qui donnera à l’aile sa meilleure efficacité.
Figure 6.32 : Coefficient de pression le long du profil d’aile. (a), répartition sur les deux surfaces ; (b),
influence de l’angle d’attaque sur le coefficient de pression le long de l’extrados. Résultats
numériques obtenus à partir du code de calculs Fluent - le profil est caractérisé par une corde
d’un mètre et l’écoulement par une vitesse de 119,7 m/s (M = 0,4, T = −50°C).
Portance 1/2ρU 2 Cz Cz
= = (6.22)
Traı̂née 1/2ρU 2 Cx Cx
162 Introduction à l’aérodynamique
(a) (b)
Figure 6.33 : Influence de l’angle d’attaque sur les coefficients de traı̂née et de portance (a), et sur la
finesse (b).
que l’on cherchera ainsi à faire tendre vers un maximum. Ce rapport est appelé finesse
et nous lui donnerons un autre sens physique plus loin. On peut dès lors tracer la
finesse en fonction de l’angle d’attaque (Fig. 6.33.b) et se rendre compte de la présence
d’un angle (unique) pour lequel l’aile est la plus efficace.
6.5.5 Polaire
Une représentation, souvent utilisée pour caractériser une aile, est celle du coeffi-
cient de portance en fonction du coefficient de traı̂née : polaire d’une aile (Fig. 6.34).
Celle-ci est graduée en incidence, la branche inférieure correspondant aux angles d’in-
cidence faibles.
Les principaux points caractéristiques d’une polaire sont (i) Mo le point de por-
tance nulle, (ii) M1 le point de finesse maximale (pente maximale pour i dans l’inter-
valle [io ; im ]), (iii) le point de portance maximale, (iv) Md le point de décrochage. On
peut remarquer qu’il existe un lien entre l’angle de finesse maximale et une consom-
mation minimale. En vol, à altitude et vitesse constantes, le poids est équilibré par la
portance et la propulsion équilibre la traı̂née :
1
Cz ρV 2 S = mg (6.23)
2
1
Cx ρV 2 S = Fp (6.24)
2
avec Fp la force de propulsion. Le rapport des deux équations précédentes donnent
Cz mg mg
= f = ou Fp = (6.25)
Cx Fp f
D’après la relation (6.25), la force de propulsion, pour un poids donné, est minimale
pour une finesse maximale. En ce point, la consommation est la plus faible.
On peut effectuer un raisonnement similaire pour connaı̂tre la configuration d’une
vitesse maximale, ceci pour une force de propulsion donnée. D’après la relation (6.24),
on peut écrire
s s
2Fp 2Fp
V= d’où VMax = (6.26)
ρSCx ρSCx,min.
Coefficients de traı̂née
On note de prime abord que le Cx est fortement lié à la taille du sillage, ce qui
confirme les propos précédents (cf. page 150).
La plaque plane placée perpendiculairement à l’écoulement génère le coefficient de
traı̂née le plus important (Cx = 2). Cette configuration est en effet la plus défavorable
à l’avancement des corps.
La figure (b) montre le cas d’un cylindre, avec pour coefficient Cx = 1,2 beaucoup
plus faible que le cas précédent. Cependant, celui-ci peut être très fortement diminué
en travaillant le profil en aval de l’obstacle. Cx = 0,12 pour le cas (c).
164 Introduction à l’aérodynamique
Figure 6.35 : Coefficients de traı̂née et répartition des traı̂nées de frottement et de pression pour différentes
formes de profil (d’après Anderson [15]).
Considérons les cas (b) et (d) présentant les mêmes valeur de Cx . Le dernier cas
considère un cylindre de diamètre d/10 entrainant donc un nombre Reynolds dix fois
moins important (on considère les autres paramètres inchangés). L’expérience montre
que Cx ne varie pas pour des nombres de Reynolds compris dans l’intervalle [104 ,
105 ] ; d’autre part la forme est conservée, le coefficient de traı̂née reste donc identique
dans les deux cas. Cependant la traı̂née dans le deuxième cas est dix fois moins élevée
que le cas (b) car la surface de référence est dix fois moins grande.
Une autre comparaison est illustrée sur les cas (c) et (d). Nous comparons ici un
large corps profilé avec un cylindre de diamètre faible. Pour le corps profilé, la traı̂née
vaut :
1 1
T = ρU 2 SCx = 0, 12( ρU 2 d)
2 2
alors que pour le cylindre, on a
1 2 1 1
T= ρU SCx = ρU 2 (0, 1d)(1, 2)Cx = 0, 12( ρU 2 d)
2 2 2
Les traı̂nées sont identiques ! Ceci illustre alors le fait que la traı̂née d’un cylindre est
équivalente à celle d’un obstacle profilé dix fois plus grand.
dans le cas de la plaque plane, alors que la traı̂née de frottement domine pour l’obs-
tacle profilé. On apportera alors grand soin à l’état de surface sur une aile d’avion
profilé, afin qu’un défaut ne vienne augmenter le faible frottement consécutivement à
un décollement prématuré.
Figure 6.37 : Influence d’un véhicule en manœuvre de dépassement sur le véhicule dépassé − Cy et Cn
représentant les coefficients d’effort latéral et de moment de lacet (d’après Noger et Széchényi
[65]).
sont toujours d’actualité et ont été complétées par des vitesses plus importantes (su-
personiques, hypersoniques,. . .).
La soufflerie utilise le principe de réciprocité, à savoir que si les conditions aux li-
mites sont respectées et l’écoulement stable, les efforts aérodynamiques de cet écoulement
sur un obstacle immobile sont identiques à ceux agissant sur un corps se déplaçant
dans un fluide immobile. On notera cependant une différence sur les différents efforts,
liée à des taux de turbulence souvent différents dans les deux cas et difficilement mai-
trisables en soufflerie.
Figure 6.39 : Étude aérodynamique d’une voiture de compétition. Photo de la soufflerie S4 de l’Institut
AéroTechnique du CNAM (Une balance à 6 composantes, est disposée sous le plancher de la chambre
d’expériences. Le véhicule à qualifier repose sur ses roues, sur quatre patins de forme rectangulaire affleu-
rant le plancher et liés au bâti mobile de la balance. Trois dynamomètres à jauges de contraintes, dont les
axes de mesure sont dirigés perpendiculairement au plancher, transmettent les poussées verticales à des
butées hydrostatiques. Trois autres dynamomètres répartis dans un même plan horizontal et orientés sui-
vant des directions orthogonales privilégiées, sont utilisés pour déterminer les composantes horizontales du
torseur des efforts.).
Description
Une soufflerie est un large venturi dans lequel l’écoulement est généré par un
ventilateur associé à un moteur. Les pales du ventilateur sont identiques à celles d’un
avion à propulsion à hélices. La soufflerie peut être ouverte (type Eiffel, Fig. 6.40.a),
ou fermée (type Götingen, Fig. 6.40.b).
En aérodynamique, pour que les mesures soient fiables, il est indispensable que
les dimensions du flux d’air que délivre la soufferie soient très supérieures à celles de
l’objet à tester. En cas contraire, se crée un effet de “bouchon” qui induit des erreurs
de mesures. À veines d’essais identiques, la soufflerie Eiffel est plus petite que la
soufflerie Götingen, et de ce fait plus simple et moins chère. Cependant, la deuxième
permet de faire varier les conditions expérimentales (pression, gaz,. . .), et reste moins
consommatrice.
Dans le cas de tests effectués sur des véhicules, les nouvelles soufferies présentent
un sol défilant et permettant ainsi de reproduire le mouvement de la route par rapport
au véhicule. Dans les souffleries classiques, plus anciennes, le flux d’air est freiné au
168 Introduction à l’aérodynamique
(a) (b)
Figure 6.40 : Souffleries ouverte type Eiffel (a) et fermée type Götingen (b).
niveau du sol, ce qui n’est pas le cas dans la réalité. Seul ce dispositif est capacle de
restituer de manière réaliste les écoulements d’air sous le plancher et sous le capot, et
ainsi de prendre en compte leur impact réel sur l’aérodynamique globale.
Mesures de la vitesse
La vitesse moyenne de l’écoulement dans la veine d’essais est estimée soit à partir
de tubes de Pitot, soit à partir de la variation de pression mesurée entre les sections
A1 et A2 : s
2( p1 − p2 )
V2 =
ρ[1 − ( A2 /A1 )2 ]
Si on désire d’autre part mesurer des vitesses locales (près des parois par exemple),
d’autres moyens de mesure, intrusifs ou non, sont disponibles.
1−→
u .( −
→
es − −
∆f = →
ei )
λ
où λ est la longueur d’onde du laser dans le milieu considéré, −
→
ei et −
→
es respective-
ment les directions du faisceau incident et la direction de l’observateur (i.e., photo-
multiplicateur).
(a) (b)
Figure 6.42 : Configurations classiques de la vélocimétrie laser à effet Doppler (a) et à franges (b).
▷ La vélocimétrie par image de particules (PIV en anglais pour Particule Image Ve-
locimetry) mesure la vitesse de particules entraı̂nées par l’écoulement à partir de deux
photos successives prises à un intervalle de temps très court. Les deux illuminations
successives de l’écoulement par deux flashs sont enregistrées par une caméra CCD,
et un post-traitement par technique de corrélation permet de déterminer les posi-
tions successives de la même particule. Un exemple de champ de vitesse à l’aval d’un
diffuseur est donné sur la figure 6.43. On y voit de façon très claire la recirculation
générée par la présence du gradient de pression adverse à l’origine du décollement
de la couche limite.
lame ou une poutre sous l’effet de la composante désirée. Le type d’agencement des
jauges permet alors la mesure de différents efforts.
Figure 6.44 : Balances à 6 composantes, la soufflerie se situant dans le plan supérieur (document Institut
AéroTechnique du CNAM).
Les frottements pariétaux peuvent être mesurés soit à l’aide d’un élément flottant
inséré sur la surface à étudier (Fig. 6.46.a), soit en utilisant des tubes (de Preston), ou
à partir de la théorie du coin d’huile (Fig. 6.46.b).
172 Introduction à l’aérodynamique
(a) (b)
Figure 6.46 : (a), mesure directe de la force de frottement à l’aide d’un élément flottant (balance de
frottement) ; (b) mesure de la forme d’un film d’huile déposé sur la paroi (d’après la théorie
de la lubrification, coin d’huile).
(a) (b)
Figure 6.47 : a, lignes de courant de l’écoulement dans la région des dents d’une roue cylindrique dentée ;
b, pertes aérodynamiques estimées expérimentalement et numériquement (d’après Pallas et
al. [68]).
10. On pourra lire le document “Modélisation de la turbulence”, proposé dans le cadre du module
éponyme en 3ème année.
174 Introduction à l’aérodynamique
réel du véhicule, contre 0,002 m2 de SCx pour une soufflerie. Ceci suffit à justifier
l’utilisation des souffleries. En fait, les deux approches sont complémentaires. Le cal-
cul, moins précis, permet ainsi d’avoir une visualisation complète de l’écoulement,
et de les analyser sous tous les angles : pressions, vitesses, etc. Il est certes possible
de visualiser les écoulements en soufflerie mais pas de façon aussi complète qu’en
numérique. D’autant que le calcul autorise des observations de l’écoulement dans des
zones difficilement accessibles en soufflerie, comme le sous-capot.
Concernant le temps, d’un côté lancer un calcul ne prend que le temps de récupérer
les données géométriques du véhicule à tester, et quelques manipulations logicielles.
Soit environ une journée pour une nouvelle configuration. Alors que la réalisation
d’une maquette nécessite au moins une semaine et coûte cher. Toutefois, une fois
passée ce cap, le bilan s’inverse. Les temps de calcul sont encore pour l’heure très im-
portants, tandis qu’une fois la maquette installée dans la soufflerie, il est possible de
tester un grand nombre de configurations et de réaliser des modifications en quelques
heures. Une facilité qui favorise l’échange entre les aérodynamiciens et les dessina-
teurs.
Ainsi, le calcul permet de préparer le travail en soufflerie, celle-ci étant réservée
aux séances de finalisation des formes.
Chapitre 7
Phénomènes transitoires dans les
conduites
réservoir par exemple, et donne alors naissance à une onde de dépression se dirigeant
vers la vanne. Il existe alors une alternance de surpressions et de dépressions (par
rapport à la pression de service) qui sont par la suite amorties par frottement. Ce
phénomène est appelé coup de bélier.
Il est important de noter que ces conditions d’alternance en pression peuvent en-
dommager le matériel ce qui n’est pas visible de suite. Après une période assez longue
d’utilisation la canalisation peut se rompre ou des brides se desserrer avant d’être ar-
rachée. L’origine du dommage est alors difficilement identifiable.
∂ ∂
(ρS) = − (ρSu) (7.1)
∂t ∂x
Cette équation généralisée laisse la possibilité aux trois variables (masse volumique,
section de passage et vitesse) d’évoluer au cours du temps et dans l’espace. Elle pourra
néanmoins être simplifiée selon les hypothèses acceptées lors des coups de bélier de
masse ou d’ondes.
Nous considérerons dans la suite l’écoulement aligné selon un axe x horizontal. Les
équations selon les deux autres directions seront donc ignorées dans les résolutions
suivantes. D’autre part les efforts volumiques de pesanteur sont alignés selon l’axe
z et orienté dans le sens négatif de ce dernier (i.e. f z = − gz). Enfin les efforts de
viscosité étant supposés négligeables ici, la composante selon x du bilan de quantité
de mouvement se résume à l’équation d’Euler :
∂u ∂u ∂p ∂
ρ +u = − − ρ ( gz)
∂t ∂x ∂x ∂x
∂p∗
= −
∂x
avec p∗ = p + ρgz. Si on ne considère que des écoulements à l’horizontal, le terme de
pression se simplifie et on obtient :
∂u ∂u ∂p
ρ +u =− (7.2)
∂t ∂x ∂x
Maintenant que les équations “générales” sont posées, intéressons-nous à différentes
configurations pratiques. Nous procéderons dans la suite par étape en traitant en
premier lieu le cas simple d’un écoulement incompressible lors de l’ouverture d’une
vanne dans une conduite de section strictement constante. Ensuite, nous examinerons
l’influence de la compressibilité du fluide et introduirons alors le concept d’onde.
Pour terminer, la conduite sera considérée comme étant élastique et capable de se
déformer sous l’influence des variations de pression.
tées par le module d’Young E. Les valeurs numériques de ces grandeurs sont ras-
semblées dans le tableau 7.1. La surface libre du bassin d’alimentation est localisée à
une altitude égale à h = 160 m au-dessus du centre de la conduite BC. Celle-ci dis-
pose d’une vanne en C initialement ouverte et générant des pertes de charge selon
l’évolution :
∆H = 640 × Q2 . (7.4)
La vanne étant partiellement ouverte et n’ayant pas encore été manipulée, il existe
un régime permanent. La vanne débouche à l’air libre en D, il règne donc en ce
point une pression atmosphérique comme cela est le cas sur la surface libre du bassin
d’alimentation (au point A). Le théorème de Bernoulli peut être appliqué entre les
points A et D, en vue d’une estimation du débit permanent :
avec ∆H les pertes de charge générées par la vanne partiellement ouverte, estimées à
partir de l’expression (7.24). La vitesse de la surface libre du réservoir est supposée
nulle, de plus h = z A − z D , il vient alors :
r
1 h
h= + 640 Q2 ≈ 640Q2 ou Q= = 0, 5 m3 /s
2gS2 640
énoncées plus haut (i.e. ρ = Cste et S = Cste). Elles s’écrivent respectivement pour la
conduite BC :
∂u ∂u ∂p
=0 et ρ =−
∂x ∂t ∂x
La première relation indique que la vitesse de l’écoulement est uniforme le long de la
conduite BC. On peut intégrer la seconde relation entre les deux extrémités, il vient
alors :
∂u
ρL = p B − pC
∂t
Le point B est placé au pied du réservoir imposant une hauteur importante, nous
supposerons alors que la pression en ce point ne varie donc pas dans le temps lorsque
la vanne est manipulée. La variation de la pression dans la conduite BC s’écrit donc
∆p = pC − p B , qui devient finalement
∂u
∆p = −ρL (7.5)
∂t
On remarquera que le terme temporel est négatif consécutivement à la décélération
de l’écoulement, générant alors une augmentation de la pression dans la conduite (i.e.
∆p > 0). On note aussi que l’évaluation de la variation de pression est fortement liée
à l’estimation de l’évolution temporelle de la vitesse. Cette dernière est généralement
obtenue à l’aide de la connaissance des pertes de charge à l’origine du changement
du point de fonctionnement.
(a) (b)
Figure 7.3 : Cheminée d’équilibre. a, installation dans un réseau hydraulique ; b, vue en coupe (d’après
Dupont [34]).
∂u ∂u g
ρLc = p B − pS = −ρghc ou = − hc (7.6)
∂t ∂t Lc
ayant pour solution un terme sinusoı̈dal. La relation (7.6) correspond presque à ce que
l’on recherche, mis à part le fait que la vitesse intervienne. Afin de faire apparaı̂tre le
terme h, utilisons la conservation du débit massique de la conduite vers la cheminée :
∂zc ∂hc
ρSu = ρSc = ρSc
∂t ∂t
∂u Sc ∂2 h c
=
∂t S ∂t2
∂2 h c Sg ∂2 h c
2
+ hc = 0 ou 2
+ Ω2 h c ( t ) = 0 (7.7)
∂t Sc L c ∂t
le niveau
p de l’eau dans la cheminée laisse donc apparaı̂tre des oscillations de pulsation
Ω = Sg/Sc Lc . La solution de cette équation est de la forme hc (t) = A × sin(Ωt + Φ).
Les paramètres A et Φ sont déterminés à l’aide de conditions à l’instant initial :
— les niveaux des deux réservoirs sont égaux (hc (0) = 0), d’où Φ = 0.
∂hc
— le débit est transmis à la cheminée, soit suo = S , entraı̂nant ainsi A =
p ∂t t=0
Suo /SΩ = uo SL/Sc g.
Ainsi, la solution devient hc (t) = uo SL/Sc g × sin (Ωt) et le maximum de surpres-
p
s
SLc
∆pMax ≡ ρghc,Max = ρg uo .
Sc g
Application numérique : si l’on reste dans le cas de la conduite précédente, avec en plus
une cheminée de diamètre trois fois plusp grand que celui de la conduite (i.e. Sc = 9S),
il vient ∆pMax = 1000 × 9,81 × 2,5 × 767/(9 × 9, 81) = 0,7 bar, au lieu des 1,5 bar
dans le cas de la conduite sans cheminée. Ce système est largement utilisé dans le cas
d’alimentation de turbine afin d’amortir les brusques variations dans la manœuvre
du distributeur de la turbine.
Coup de bélier d’ondes 185
Ces résultats ne sont valables que pour des temps de fermeture ou d’ouverture
lents. L’application de ces lois dans le cas d’ouverture plus courte indiquerait
des surpressions beaucoup plus grandes que celles observées. En effet, si les va-
riations de pression sont suffisamment rapides et importantes, l’eau devient un
fluide compressible susceptible de propager des ondes à grande vitesse. Le reste
du chapitre s’intéresse aux effets de cette compressibilité.
1 dV
χ=− (7.8)
V dp
1 dρ
χ= (7.9)
ρ dp
2 ∂p 1
c = = (7.10)
∂ρ ρχ
Cette expression met en évidence des célérités différentes dans l’air, l’eau ou l’acier
(Tableau 7.2).
Le prochain paragraphe se propose d’inspecter la propagation d’une onde issue
d’une perturbation locale, mais toujours en considérant la conduite pourvue d’une
section constante (n’allons pas trop vite).
Une fois celle-ci arrêtée, la deuxième tranche subit la même décélération, et la même
compression. Cette compression est véhiculée de proche en proche et ne s’apparente
plus alors à un mouvement en bloc mais plutôt à la propagation d’une onde dont on
va estimer la célérité 2 .
Le paragraphe 7.3 nous a permis de mettre en place les équations générales qui
gouvernent le mouvement des particules. Dans le cas de fluide compressible, les
équations de continuité et d’Euler deviennent respectivement (en tenant compte du
fait que la section reste constante) :
∂ρ ∂u ∂ρ ∂u ∂u ∂p
+ρ +u =0 et ρ +u + =0
∂t ∂x ∂x ∂t ∂x ∂x
On supposera que les termes de convection sont négligeables devant les autres termes :
∂ρ ∂ρ ∂u ∂u
u ≪ u ≪
∂x ∂t ∂x ∂t
Ceci est légitime à partir de l’instant où la vanne est fermée, le déplacement des
particules est très inférieur à la variation des paramètres au passage de l’onde. Il reste
alors :
∂ρ ∂u
+ρ =0 (7.11)
∂t ∂x
∂u ∂p
ρ + =0 (7.12)
∂t ∂x
Les équations (7.11) et (7.12) font apparaı̂tre trois paramètres, ρ, u et p. On peut en
éliminer un à l’aide de la relation (7.10) :
∂ρ ∂ρ ∂p 1 ∂p
= × = 2
∂t ∂p ∂t c ∂t
1 ∂p ∂u
2
+ρ =0 (7.13)
c ∂t ∂x
Intéressons-nous maintenant à la différence ∂(7.12)/∂x − ∂(7.13)/∂t, on trouve la re-
lation :
2. Dans ce cas, on parlera de la célérité au lieu de la vitesse de l’onde, la dénomination de vitesse
étant associée à un déplacement de matière.
188 Phénomènes transitoires dans les conduites
∂2 p 1 ∂2 p
− =0 (7.14)
∂x2 c2 ∂t2
∂2 u 1 ∂2 u
− =0 (7.15)
∂x2 c2 ∂t2
Nous reconnaissons ici l’équation de propagation d’une onde. Cette relation intervient
aussi dans les études concernant les cordes vibrantes, les ondes acoustiques, le rayon-
nement électromagnétique, etc. Il s’agit dans tous les cas d’une équation différentielle
du second ordre du paramètre étudié ϕ :
∂2 ϕ 1 ∂2 ϕ
− =0
∂x2 c2 ∂t2
La solution contient deux termes, F (t − x/c) et f (t + x/c) . Le premier terme ca-
ractérise une onde qui se propage vers les x positifs, c’est à dire dans le sens de
l’écoulement permanent (onde progressive) ; le deuxième une onde se propageant
dans le sens inverse à l’écoulement (onde régressive).
La vitesse d’agitation des particules fluides u et la pression p ont donc des solu-
tions équivalentes :
x x
p − pini = F t − + f t+ (7.16)
c c
1 h x x i
u − uini = F t− − f t+ (7.17)
ρc c c
où l’indice ini indique des conditions initiales, c’est-à-dire les conditions observées
avant le départ de l’onde. Ces conditions peuvent donc être différentes d’une extrémité
à l’autre d’une conduite. Elles peuvent même évoluer dans le temps pour une même
extrémité. Les fonctions F et f doivent être homogènes à une pression. Ces expres-
sions, appelée dans la littérature Formules d’Allievi 3 , peuvent être vérifiées par la sa-
tisfaction des équations (7.12) et (7.13). C’est à partir de ces équations que la méthode
graphique de Bergeron a été introduite (§ 7.8).
Les fonctions F et f représentant des ondes de pression se propageant sans déforma-
tion respectivement dans les sens positif et négatif, à une abscisse x dans la conduite,
la pression p − pini est la superposition de ces deux ondes.
3. Lorenzo Allievi (1856-1941), ingénieur italien ayant travaillé sur les problèmes transitoires dans
les centrales hydroélectriques.
Coup de bélier d’ondes 189
∂2 p 1 ∂2 p
− =0
∂x2 a2 ∂t2
avec a la célérité des ondes telles que :
v
u 1
a=u (7.18)
u
D
ρ χ+
t
eE
Prenons le temps d’examiner les expressions (7.10) et (7.18), expressions des célérités
d’une onde dans une conduite de section strictement constante et dans une conduite
élastique. Tout se passe comme si le coefficient de compressibilité était augmenté
dans le deuxième cas par D/eE, de ce fait la célérité de l’onde dans la conduite est
inférieure à celle de cette même onde dans le même milieu libre. De plus, l’influence
de la déformation de la conduite n’est appréciable que pour les fluides pratiquement
incompressibles, pour lesquels χ et D/eE sont à peu près du même ordre de grandeur.
Elle est négligeable pour les gaz.
Application numérique : χ = 0, 5.10−9 m2 /N, ρ = 103 kg/m3 , e=2 mm, D=200 mm,
E = 2, 0 × 1011 N/m2 :
célérité (m/s)
Section constante Section variable
1414 1000
une pression po . Ainsi on peut remplacer dans Les solutions d’Allievi les grandeurs
initiales, i.e. pini = po et uini = uo . Ces solutions s’écrivent, si on tient compte en plus
du fait qu’il n’y a pas d’onde progressive F et que la vanne impose u = 0 :
x 1 h x i
p − po = f t + = ∆p et 0 − uo = −f t +
a ρa a
L’onde régressive est donc une onde de surpression f = ρauo remontant le long de la
conduite (Fig. 7.5.a), la pression dans la conduite derrière l’onde valant p = po + ρauo .
La surpression observée par la conduite par rapport à la pression nominale, i.e. po ,
vaut donc p − po = ρauo ou ∆h = auo /g lorsqu’elle est exprimée en mètre de colonne
d’eau.
L’onde régressive est donc une onde de dépression f = −ρauo remontant le long de la
conduite (Fig. 7.5.e), la pression dans la conduite derrière l’onde valant p = po − ρauo .
La dépression observée par la conduite par rapport à la pression nominale, i.e. po , vaut
donc p − po = −ρauo ou ∆h = − auo /g lorsqu’elle est exprimée en mètre de colonne
d’eau.
192 Phénomènes transitoires dans les conduites
↬ À t = 3τ, l’onde régressive atteint le réservoir (Fig. 7.5.f). Toute la conduite est en
dépression, et l’écoulement s’arrête. Une dépression régnant dans la conduite, l’eau
du bassin rentre dans la conduite. Il y a naissance d’une onde progressive au niveau
du bassin :
po − ( po − ρauo ) = F ⇒ F = ρauo
1 x
u−0 = F t− ⇒ u = uo
ρa a
↬ À t = 4τ, le cycle est fini, l’onde progressive atteint la vanne. Un nouveau cycle
prend place.
4L
T= (7.19)
a
L’ensemble de ce qui vient d’être dit est résumé sur la figure 7.6. Pour connaı̂tre les
conditions de pression et de vitesse dans la canalisation au point P d’abscisse x à
l’instant t, il faut se souvenir que celles-ci correspondent à la superposition d’ondes
progressive et régressive. La première est partie du réservoir à l’instant t − x/a avec
une vitesse + a, et la deuxième de la vanne à l’instant t + ( L − x )/a avec une vitesse
− a. La résolution de ce problème passe par l’écriture des deux équations d’onde et
de leur condition limite respective.
Application numérique : on note dans cette application des supressions et des dépres-
sions ayant des valeurs ±ρauo = ± 29,4 bar. Cette valeur numérique très élevée est
liée à la production d’une seule onde produite par fermeture instantanée, ce qui
n’est jamais observé dans les installations réelles. La présence de plusieurs ondes,
progressives et régressives, entraı̂ne en effet des interférences diminuant l’effet tran-
sitoire. D’autre part, on observe généralement une atténuation de ce phénomène
consécutivement aux pertes de charge et aux imperfections de l’élasticité de la conduite
(figure 7.7). Néanmoins cette décroisance est lente car la transformation d’énergie en
chaleur de frottement suite aux phénomènes précédents est peu importante, et les
variations de pression sont alors persistantes.
Coup de bélier d’ondes 193
Figure 7.5 : Coup de bélier lors d’une fermeture instantanée de la vanne située à l’extrémité de la conduite.
Figure 7.6 : Descriptif du comportement des ondes au cours du coup de bélier (notation : h = p/(ρg).)
Q0 − Q1
∆p = ρa(Vo − V1 ) = ρa (7.20)
S
Cette relation est simple d’utilisation pour des évolutions consécutives de la vitesse
tant que la première onde générée, réfléchie par l’autre extrémité, n’est pas de retour
sur la vanne.
Considérations pratiques − identification du type de coup de bélier 195
vrai pour des fermetures rapides de la vanne. En effet, le coup de bélier en masse ne
fait pas intervenir le matériau de la conduite.
Si les analyses montrent qu’un coup de bélier peut apparaı̂tre dans une installa-
tion, l’ingénieur peut faire appel, en premier lieu, à différentes solutions : (i) augmen-
ter le temps de modification de l’écoulement, (ii) augmenter la classe de la conduite
permettant une meilleure résistance à la pression, (iii) limiter la vitesse de l’écoulement,
(iv) réduite la vitesse de propagation de l’onde en changeant le type de conduite et
(v) placer des soupapes de sécurité.
Il existe d’autre part des moyens de protections dont le rôle est de limiter la
dépression et la surpression, le but n’étant pas de supprimer le coup de bélier car
ceci est impossible. Les moyens les plus régulièrement utilisés sont le volant d’iner-
tie, la soupape de décharge, les réservoirs d’air et les cheminées d’équilibre (étudiées
précédemment).
Dans le cas du volant illustré sur la figure 7.9, le moment d’inertie de ce volume de
révolution par rapport à l’axe de rotation O−
→z vaut
R22 + R21
I=m
2
De ce fait, au moment de l’arrêt brusque du moteur, le volant d’inertie a accumulé
une énergie cinétique égale à 1/2 × Iω 2 . D’après le théorème de l’énergie cinétique, la
variation de l’énergie cinétique est égale à la somme des travaux des efforts extérieurs,
soit ici Cωdt pendant le temps dt. Par conséquence, d’après ce théorème :
1 2
d Iω = Cωdt
2
En développant, et en passant aux différences finies, cette relation s’écrit finalement :
C
∆ω =
∆t (7.21)
I
Ainsi, à partir de la vitesse angulaire au moment de la disjonction, il est possible
d’estimer les vitesses successives (ω − ∆ω) à des temps consécutifs (t + ∆t), et les
nouvelles caractéristiques de la pompe en fonction de ces nouvelles vitesses de ro-
tation. On choisit généralement la valeur de ∆t égale au temps d’un aller-retour de
l’onde entre la pompe et le réservoir, soit ∆t = 2L/a.
(a) (b)
Figure 7.10 : a, soupape de sureté anti-bélier (produit Dresser Consolidated) ; b, Principe de la dispo-
sition d’un réservoir d’air.
battant percé.
Si notre imagination est assez vaste pour pouvoir considérer un observateur se déplaçant
dans le sens des x positifs à la vitesse a (Fig. 7.15), son déplacement vaut :
x
x = at + Cste ou t− = Cste
a
x
Dans ce cas, F t − = Cste, entrainant donc la relation importante
a
Méthode graphique de Bergeron 203
a
h+ Q = Cste (7.22)
gS
Ici la constante représente les conditions de départ du coureur et peut s’exprimer sous
la forme hinit. + aQinit. /gS, avec hinit. et Qinit. les conditions de hauteur et de débit au
moment du départ du relayeur considéré. Si notre imagination a résisté jusqu’ici, il
nous sera facile de considérer cette fois-ci un autre observateur se déplaçant en sens
inverse (Fig. 7.15), mais toujours à la vitesse a, son déplacement vaut :
x
x = − at + Cste ou t+ = Cste
a
x
Dans ce cas, f t + = Cste, entrainant donc une deuxième relation tout aussi im-
a
portante
a
h− Q = Cste (7.23)
gS
Nous voyons donc que pour un observateur se déplaçant à une vitesse a (quel que
soit le sens), la charge et le débit volumique sont reliés linéairement.
Dans cet exemple, nous traiterons de façon simultanée le problème de façon ana-
lytique en utilisant les relations portées par les relayeurs et de façon graphique en
traçant l’épure de Bergeron basée sur des segments.
nous allons calculer plus loin. À un instant infiniment petit avant la fermeture, un
premier relayeur quitte la vanne en direction du réservoir, il court à une vitesse a
sur une droite d’équation h − aQ/( gs) = Cste (il s’agit du relayeur à gauche portant
la relation 7.23). Lorsqu’il quitte la vanne, ses conditions initiales ( Qinit. , hinit. ) étaient
(0,5 m3 /s, 160 m), d’où l’équation qu’emprunte le relayeur 1 : h − aQ/( gs) = −140.
↬ À l’instant t = 4L/a, le relayeur 4 arrive sur la vanne pour laquelle le débit est
nul. La pression alors h = −140 m d’après le tracé ou l’équation du relayeur. Nous
obtenons ainsi le point (0 m3 /s, −140 m). À partir de cet instant, le quatrième va
s’empresser d’inverser la surpression, on retrouve alors les conditions du premier
relayeur. Un cycle complet a été effectué sur la vanne. Le dernier relayeur quant à
lui, partant de la vanne pour se diriger vers le réservoir lui permettra d’observer la
condition (0,5 m3 , 160). Lorsqu’il inversera la tendance, un cycle complet sera alors
effectué à l’embouchure du réservoir. Le relai est fini.
Notons que nous retrouvons une période de 4L/a que ce soit au niveau de la
vanne mais aussi à l’embouchure du réservoir. Nous retrouvons d’autre part, comme
il a été dit plus haut, les mêmes estimations de surpression que celle prédı̂te par les
formulations d’Allievi. Les résultats obtenus à partir des relayeurs de Bergeron sont
rassemblés sur la figure 7.17.
Notons tout de même que tracer les droites caractéristiques des relayeurs est une
première étape, et non une fin en soi. En effet, il est plus important dans une deuxième
étape de fournir les renseignemments de débits et de hauteurs aux extrémités pour
des temps caractéristiques du problème. C’est là le but de la méthode de Bergeron.
Pour cela, on simplifie en plus généralement l’écriture temporelle, c’est à dire qu’un
temps noté t = 1 équivaut au temps t = L/a. De même, on renseigne sur l’épure de
Bergeron les extrémités considérées. Ainsi la notation ] − ∞, 1] R signifie que pour des
instants antérieurs à t = L/a, le débit vaut 0,5 m3 /s et la surpression à l’embouchure
du réservoir est nulle d’après la figure 7.17.
On note donc une surpression égale à aQo /gS sur la vanne entre les temps 0 et
L/a, comme il est indiqué sur la figure 7.5 obtenue à partir des solutions d’Allievi.
Insistons sur le fait que tout ce qui vient d’être dit dans ce chapitre est basé sur les
hypothèses de réversibilité complète du phénomène : absence de perte de charge,
élasticité parfaite de la conduite, et écoulement monophasique 4 . Dans le cas réel, les
dissipations atténuent cet effet dans la durée. D’autre part, les canalisations présentent
la plupart du temps des changements de direction à l’orgine de réflexions d’ondes.
En tout point de la canalisation la pression est ainsi la superpostion de deux ondes,
ce qui diminue généralement l’amplitude de la pression.
Figure 7.17 : Représentation des surpressions et dépressions au niveau de la vanne, et des débits au niveau
du réservoir dans le cas d’une fermeture instantanée d’une vanne placée à l’extrémité d’une
conduite d’adduction.
2
Q
∆H = 640 (7.24)
t
1−
T
Le réservoir est de grande taille, on peut donc considérer son niveau comme étant
constant, les pertes de charge seront négligées dans la conduite BC. Comme il a été
montré précédemment, le débit permanent vaut Qo = 0,5 m3 /s.
On ne s’intéresse qu’aux caractéristiques de la vanne aux instants 0, 2τ, 4τ, 6τ, etc.
correspondant à des multiples de l’aller-retour de la première onde générée, avec 2τ
= 2L/a = 1,3 s. Les caractéristiques de la vanne sont rassemblées dans le tableau 7.3 et
les courbes sont tracées sur la figure 7.18 sur laquelle la condition du réservoir a été
rajoutée (h = 160 m).
Tableau 7.3 : Courbes caractéristiques de la vanne pour des temps multiples d’un aller-retour de la première
onde générée.
Instant Caractéristique
0 ∆H = 640 × Q2
2τ ∆H = 1372 × Q2
4τ ∆H = 4782 × Q2
6τ ∆H = 268 960 × Q2
Au-delà Q=0
La méthode graphique de Bergeron consiste à faire partir des relayeurs des extrémités
de la conduite BD. L’orientation de l’axe des droites est de B vers C étant donnée
l’orientation de l’écoulement permanent. La constante des relations portées par les
relayeurs sera déterminée par les conditions de l’extrémité de départ. Les conditions
imposées par l’extrémité d’arrivée seront soit une hauteur (pour le réservoir), soit une
valeur de perte de charge et de débit (pour la vanne).
210 Phénomènes transitoires dans les conduites
La pente des relations des relayeurs vaut (Cf. paragraphe 7.8.4) 600 m/(m3 /s), les
relayeurs courent sur des droites de pente égale à ± 60 m par 0,1 m3 /s.
Les conditions imposées par le réservoir entraı̂ne Q = 0,12 m3 /s et h = 160 m, que l’on
retrouve graphiquement.
Figure 7.19 : État de la pression et du débit lors d’une fermeture progressive d’une vanne.
214 Phénomènes transitoires dans les conduites
La solution de cette application est résumée dans le tableau 7.6. Toutes les lignes de
ce tableau ont pour but d’estimer la décroissance de la vitesse de rotation de l’arbre-
moteur, ∆ω. On recommence ensuite une nouvelle ligne à partir d’une nouvelle valeur
ω − ∆ω. On effectue un calcul tous les deux fois ∆t (en prenant pour unité de temps
Méthode graphique de Bergeron 215
∆ω = 0, 194 × C
Figure 7.21 : Résolution graphique de l’application utilisant un volant d’inertie sur un système de pompage.
a pV ( T )
ho − Qo < − 10, 33 (7.25)
gS ρg
où la valeur 10,33 correspond à la pression atmosphérique exprimée en mètre de
colonne d’eau. Pour des conditions de température modérée (i.e. 20 ˚ C), la pression
de vapeur saturante est proche de 2340 Pa entrainant alors pV ( T )/ρg − 10, 33 ≈ −10
m. Il s’agit là d’un ordre de grandeur pour la quantité ho − aQo /gS à ne pas atteindre.
On note que le débit au point de cavitation est négatif durant la période [2,4].
Le signe négatif provient de l’orientation de l’axe des relayeurs orienté du réservoir
vers la vanne. Ce signe devient positif si on considère le cas d’une alimentation d’un
réservoir par une pompe pour lequel l’axe des relayeurs est alors orienté de la pompe
vers le réservoir. Dans les deux cas de figure, le volume de la poche gazeuse croı̂t du
point de changement de régime et remonte vers le réservoir.
Méthode graphique de Bergeron 217
(a) (b)
Figure 7.22 : Tracé de l’épure de Bergeron en présence de cavitation lors d’une fermeture rapide d’une
vanne. a, intersection de la droite d’un relayeur avec les conditions de cavitation ; b, épure
finale (Flèches en trait normal, premier cycle ; flèches en trait gras, deuxième cycle).
où t f représente le temps de fermeture tel que t f > 2L/a. Considérons le coureur de
Bergeron partant du réservoir au temps L/a, l’équation empruntée par ce coureur est
a
h = ho − ( Q − Qo ). Lorsqu’il arrive sur la vanne, à l’instant t = 2l/a, la pression
gS
vaut d’après la relation précédente :
2L a 2L
h = ho − Q − Qo
a gS a
218 Phénomènes transitoires dans les conduites
2LQo
Après simplification, on peut écrire la valeur de la surpression, ∆h = h − ho = .
gSt f
On obtient alors la formule de Michaud :
2ρLuo
∆p =
tf
Rappelons-nous cependant que cette formule est basée sur une fermeture linéaire de
la vanne, ce qui n’est pas le cas dans la pratique. La formule de Michaud conserve
uniquement un intérêt théorique.
Bibliographie sommaire
— Les Pompes
⋄ Know and Understand Centrifugal Pumps
L. Bachus & A. Custodio, Ed. Elsevier
⋄ Techniques de l’Ingénieur
⋄ Les pompes centrifuges - Conceptions, règles d’utilisation
Encyclopédie des sciences et des techniques industrielles, Section techniques
modernes. Édition Nathan.
⋄ Machine à fluide - Principe et fonctionnement
M. Pluviose, Ed. Ellipse
⋄ Pompes, ventilateurs, compresseurs
A. de Kováts et G. Desmur, Ed. Dunod
⋄ Les pompes centrifuges
M. Chicha, Les presses de l’Université de Montréal
⋄ Machines hydrauliques
L. Bergeron, Dunod
⋄ Mémento Technique
lmt, Division Électro-Hydraulique
⋄ Les pompes centrifuges - Conceptions, règles d’utilisation
Pompes Guinard/Nathan
⋄ Hydraulique urbaine - Tome II
220 Bibliographie sommaire
— Les Vannes
⋄ Hydraulics of pipelines − Pumps, valves, cavitation, transients
J. P. Tullis, 1989, John Wiley & Sons, Inc.
⋄ Valve Handbook (2nd Edition)
P. L. Skousen, 2004, McGraw-Hill Handbooks.
[17] S. Barber, S. Haake, and M. Carré. Using CFD to understand the effects of seam
geometry on soccer ball aerodynamics- The engineering of Sport 6 - Volume 2 : develop-
ments and disciplines. Springer, 2006.
[18] D. Bergeron. Environmental Engineers’ Handbook. CRC Press, 2nde edition, 1997.
[19] L. Bergeron. Machines hydrauliques. Dunod, 4ème edition, 1928.
[20] H. Blasius. Grenzschichten in Flüssigkeiten mit kleiner Reibung. Z. Math. Physik,
56 :1–37, 1908.
[21] W. Bohl. Stromungsmaschinen Berechnung und Konstruction. Vogel Buchverlag
Wurzburg.
[22] B. E. Bosserman, A. Ali, and I. M. Schuraytz. Butterfly valves : torque, head loss
and cavitation analysis − Manual of water supply practices. M49 Second Edition,
2001.
[23] C. E. Brennen. Hydrodynamics of pumps. Concepts ETI 1994 / Oxford Science
Publications, 1994.
[24] E. Buckingham. On physically similar systems. Illustrations of the use of dimensional
equations. Physical Review, 4 :345–376, 1914.
[25] R. Camerer. Z. ges. Turbinenwes,. 1915.
[26] A.D.S. Carter and H.P. Hughes. A Theoretical Investigation Into the Effect of Profile
Shape on the Performance of Airfoils in Cascade. ARC Report No. 2384, 1946.
[27] DS. Chamieh. Forces on a whirling centrifugal pump-impeller. PhD thesis, Calif.
Inst. of Tech. (Division of Engineering and Applied Science), 1983.
[28] R. Comolet. Mécanique expérimentale des fluides - Tomme II : Dynamique des fluides
réels. Masson, 4ème edition, 1994.
[29] P. Coussot and J.L. Grossiord. Comprendre la rhéologie, de la circulation du sang à
la prise du béton. EDP Sciences, 2001.
[30] T.N. Crouch, D. Burton, N.A.T. Brown, M.C. Thompson, and J. Sheridan.
Flow topology in the wake of a cyclist and its effect on aerodynamic drag. J. Fluid Mech.,
748 :5–35, 2014.
[31] Collectif d’auteurs. Memento Technique. Edition Le Matériel Téléphonique −
Division Électro-hydraulique, 1975.
[32] S.L. Dixon. Fluid Mechanics and Thermodynamics of Turbomachinery. Butterworth-
Heinemann, 4ème edition, 1998.
[33] L. Duclos. Pompes hydro et appareils d’élévation.
[34] A. Dupont. Hydraulique urbaine. Eyrolles, 2ème edition, 1969.
[35] L. P. Erm. Dynamic testing of aricraft models in water tunnel. 15th Australasian Fluid
Mechanics Conference, Sidney, Australia, 13-17 Dec., 2004.
[36] H. Bachmann et al. Vibration problems in structure : practical guidelines. Birkhauser,
1995.
[37] J. Fay. Industrie et Technologie, 910, 2009.
BIBLIOGRAPHIE 223
[56] E. marchand. Nouvelle théorie des pompes centrifuges − Étude théorique et pratique.
Librairie E. Bernard & Cie, 53ter Quai des Grands-Augustins, 2nde edition, 1896.
[57] J.C. Martin, C.J. Davidson, and E.R. Pardyjak. Understanding sprint-cycling per-
formance : the integration of muscle power, resistance, and modeling. Int. J. of Physiol.
Perform., 2(1) :5 pages, 2007.
[58] P.J. McNulty and I.S. Pearsall. Cavitation inception in pumps. In Proceedings
of ASME Int. Symp. on Cavitation Inception, 1979.
[59] P. Morel. Vers une meilleure maitrise des écoulements dans les inducteurs. La houille
blanche, 3/4, 2001.
[60] Y. Moulin. Christian hugues, minix - dispositif de réduction de vortex et de
traı̂née induite, http ://[Link], 2016.
[61] Y. Moulin. Yoann moulin, photographe, http ://[Link], 2016.
Accédé le 01/02/2017.
[62] Y. Nakayama. Visualised Flow. Pergamon Press, Oxford, 1988.
[63] Y. Nakayama and R. F. Boucher. Introduction to fluid mechanics. Butterworth-
Heinemann, 1999.
[64] J. Nikuradse. Strömungsgesetze in Rauhen Rohren. Forschungsheft, Série B (n°361),
1933.
[65] C. Noger and E. Széchényi. Influence of vehicule shape on transient aerodyna-
mic forces generated during passing manoeuvre. In Proceedings of SIA − Fluid
Dynamics Applications in Ground Transportation, Lyon, Oct. 2005.
[66] T. Oniga. Calcul des tuyaux. Matemine−France, 1949.
[67] P. Pagnier, R. Morel, J.Y. Champagne, and C. Henry. Prédiction de l’apparition de
la cavitation dans un inducteur placé en amont d’une roue centrifuge. La houille blanche,
5 :13–21, 1999.
[68] S. Pallas, Y. Marchesse, C. Changenet, and F. Ville. A windage power loss model
based on CFD study about the volumetric flow rate expelled by spur gears. Mechanics
and Industry, 13 :317–323, 2012.
[69] J.L.M. Poiseuille. Recherches expérimentales sur le mouvement des liquides dans les
tubes de très petits diamètres. Comptes Rendus, 11 :961–967, 1840.
[70] E. Portier. Étude de l’écoulement au travers de diffuseurs équipés de vannes de gui-
dage ; analyse des performances et application au contrôle des débits pulsés. Thèse de
l’Université de Poitiers, 2000.
[71] X. Qiui, D. Japikse, J. Zhao, and M.R. Anderson. Analysis and Validation of a
Unified Slip Factor Model for Impellers at Design and Off-Design Conditions. J. of
Turbomachinery, 133 :041018–1 041018–9, 2011.
[72] R. Rey, F. Bakir, and J. Poulain. pompes rotodynamiques − Fonctionnement. Tech-
niques de l’Ingénieur, BM4302, 2012.
[73] O. Reynolds. An experimental investigation of the circumstances which determine
whether the motion of water shall be direct or sinuous, and the law of resistance in
parallel channels. Royal Society of London − Philosophical Transactions, 174, 1883.
BIBLIOGRAPHIE 225
[74] S. Ristic. Flow Visualisation Techniques in Wind Tunnels, Part II Optical Methods.
Scientific Technical Review, 57(2) :38–49, 2007.
[75] S. Ristic. Flow Visualisation Techniques in Wind Tunnels, Part I Non Optical Methods.
Scientific Technical Review, 57(1) :39–49, 2007.
[76] R. H. Sabersky, A. J. Acosta, and E. G. Hauptmann. Fluid Flow. McMillan Publ.
Co., 3ème edition, 1989.
[77] C. Scruton. An Introduction to Wind Effects on Structures. Oxford University
Press, 1981.
[78] M. S. Selig, M.D. Maughmer, and D.M. Somers. Natural-Laminar-Flow airfoil for
General-Aviation Applications. Journal of Aircraft, 32(4) :pp. 710–715, 1995.
[79] H. Stel, G.D.L. Amaral, C.O.R. Negrao, S. Chiva, V. Estevam, and R.E.M. Mo-
rales. Numerical Analysis of the Fluid Flow in the First Stage of a Two-Stage Centri-
fugal Pump With a Vaned Diffuser. Journal of Fluids Engineering, 135 :pp. 071104–1
071104–9, 2013.
[80] A. Stodola. Steam and Gas Turbines. McGraw-Hill New York, 1927.
[81] P.K. Swamee and A.K. Jain. Explicit equations for pipe-flow problem. J. of Hydraulic
Division, 102(5) :657–654, 1976.
[82] traduit du russe par Mme M. Meury. Memento des pertes de charge − Coefficients
de pertes de charge singulières et de pertes de charge par frottement. Eyrolles, 1960.
[83] A. Troskolanski. Les turbopompes : théorie, tracé et construction. Eyrolles, 1977.
[84] J. P. Tullis. Hydraulics of pipelines − Pumps, valves, cavitation, transients. John
Wiley & Sons, Inc., 1989.
[85] A. Turner. Leçons de Physique expérimentale de l’abbé Nollet (livre III, leçon IX).
[86] M. Van Dyke. An album of fluid motion. Parabolic Press, Stanford, California.,
1982.
[87] A. Vaschy. Sur les lois de similitude en physique. Annales Télégraphiques, 19 :25–28,
1892.
[88] D. Voisin. Étude des effets du vent sur les grues à tour. Thèse de l’Université de
Nantes, 2003.
[89] Brüel & Kjaer (Magazine WAVES). Tunnel vision, 2016.
[90] W.L. Weiss. Suggested Design Change for a Centrifugal Pump Impeller handling Dred-
ged Mud. Report for C.E. 422 - Hydraulic Research, 1959.
[91] R. Whitcomb. A design approach and selected wind-tunnel results at high subsonic
speeds for winglets wing-tip mounted winglets. Nasa Technical Note (TN D-8260),
1976.
[92] F.M.. White. Fluid Mechanics − Mcgraw-Hill Series In Mechanical Engineering (4th
Edition). McGraw-Hill, 1998.
[93] C. Wieselberger. A sphere drag. Zeitschrift fur Mechanisch, 5 :140–144, 1914.
226 BIBLIOGRAPHIE
[94] F.J. Wiesner. A review of slip factor for centrifugal compressors. Journal of Engineering
Power, Transactions of the American Society of Mechanical Engineers, 89 :558–572,
1967.
[95] H. Yamazaki. . Journal of Instrumentation and Control, 10 :173, 1971.