0% ont trouvé ce document utile (0 vote)
357 vues183 pages

Similitude de Froude en turbomachines

Ce document présente les paramètres de similitude utilisés dans l'analyse des turbomachines. Il décrit des nombres adimensionnels comme le coefficient de charge et de débit qui permettent de caractériser les machines et prédire leur performance à partir de modèles de similitude.

Transféré par

MarwaEnnouri
Copyright
© © All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats PDF, TXT ou lisez en ligne sur Scribd
0% ont trouvé ce document utile (0 vote)
357 vues183 pages

Similitude de Froude en turbomachines

Ce document présente les paramètres de similitude utilisés dans l'analyse des turbomachines. Il décrit des nombres adimensionnels comme le coefficient de charge et de débit qui permettent de caractériser les machines et prédire leur performance à partir de modèles de similitude.

Transféré par

MarwaEnnouri
Copyright
© © All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats PDF, TXT ou lisez en ligne sur Scribd

Si vous avez connu un, vous les avez connu tous.

(TERENCE,Phormio )

Citation du livre Fluid Mechanics and


Thermodynamics of Turbomachinery de [Link]
• Présenter les paramètres de similitude employés dans les
turbomachines

• Décrire des formes particulières de ces coefficients pour divers


types machines

• Regarder les cartes utilisées pour étudier les compresseurs et


les turbines

• Présenter le concept de vitesse spécifique

3
La spécification fondamentale d’une turbomachine concerne:

• le débit
• la nature du fluide(gaz, liquide)
• l’état à l’entrée
• la quantité de travail à échanger
À partir de cette spécification on peut:

• sélectionner le type de machine (axiale, radiale, mixte) ;


• estimer sa performance
• effectuer un dimensionnement préliminaire

La mise à l’échelle d’une machine existante, en réponse à une


nouvelle spécification, est une démarche utilisée en industrie
visant ce dimensionnement préliminaire. Cette étape permettra la
réduction de la durée des cycles de développement
L’analyse dimensionnelle fournit une support théorique à cette
première étape de conception. Elle incorpore des nombres
adimensionnels dont leur valeur peut être interprétée comme le
ADN de chaque appareil
La richesse du formalisme des nombres sans dimension permet
la génération de familles de machines dites similaires, la
prédiction de leur performance et un guide pour la construction de
prototypes
Dans cette démarche on utilise des quantités telles que les
coefficients de charge Ψ et de débit Φ, parmi d’autres
Le théorème π de Vaschy-Buckingham affirme que toute loi
physique peut se réduire à un nombre minimal d’inconnues en
introduisant des variables sans dimension (notées πi )

Si un phénomène est décrit par N variables physiques, exprimées


à l’aide de J dimensions fondamentales, il est alors possible de
définir un problème adimensionnel équivalent comme une
fonction de N-J grandeurs sans dimension

7
Modèle virtuel Maquette de laboratoire Prototype

8
Pour effectuer une conception agile, on vise le transport des
caractéristiques connues d’une turbomachine vers une autre.
Pour ce faire on a besoin que:

les dimensions respectives soient à l’échelle


les rapports de vitesses soient égaux (les lignes du
courant soient similaires)
les rapports des forces inertie/autres soient égaux

On doit alors satisfaire trois conditions de similitude différentes:


géométriques, cinématiques et dynamiques
 l   l 
  = 
Rapport de longueurs  D m  D  p
modèle prototype

 Cx   Cx 
φ =  = 
Rapport de vitesses
 U m  U  p

 Fi   Fi 
  =  Fx = ?
Rapport de forces
 Fx  m  Fx  p
La similitude géométrique, en principe la plus simple, se voit
confrontée à des difficultés pour la satisfaire pleinement dans
certains endroits d’une turbomachine, comme par exemple le jeu
de tête (séparation entre la pale en sortie et le carter) ou la
rugosité relative de la surface.
Pour cette raison, des corrections empiriques (souvent
confidentielles) sont appliquées par les manufacturiers, pour
transposer les données obtenus lors des essais sur le modèle
vers le prototype
En mécanique des fluides, la similitude dynamique regarde
principalement le rapport des forces d’inertie à d’autres forces
En génie on s’intéresse fréquemment aux forces visqueuses,
gravitationnelles et élastiques
Avec l’inertie, ces trois forces donnent lieu respectivement, aux
nombres adimensionnels de Reynolds, Froude et Mach
On peut déjà noter que l’écoulement dans turbomachine n’est pas
à surface libre, de sorte la force gravitationnelle n’occupe pas un
rôle prépondérant dans l’analyse
Reynolds
𝑉𝑉𝑉𝑉𝑉𝑉
rapport des forces d’inertie aux forces visqueuses
𝜇𝜇
Froude
𝑉𝑉
rapport des forces d’inertie aux forces gravitationnelles
𝑔𝑔𝑔𝑔
Mach
𝑉𝑉
rapport des forces d’inertie aux forces élastiques
𝑐𝑐
cx , m 𝑚𝑚̇ ⁄𝜌𝜌 𝐴𝐴
φ= =
U 𝑁𝑁 𝐷𝐷⁄2

𝑁𝑁: vitesse de rotation


m
φ= 𝐷𝐷: diamètre du rotor
ρND 3
We
ψ= 2
U

 We 
ψ= 2 2 
N D 
Pompes


m m 𝑄𝑄
Q= φ= Φ=
ρ ρND 3
𝑁𝑁𝐷𝐷3

We 𝑔𝑔𝑔𝑔
We = gH ψ= 2 2 Ψ= 2 2
N D 𝑁𝑁 𝐷𝐷
Lors de l’étude des machines centrifuges on a trouvé la relation:

Ψ = (1 − Φ𝑡𝑡𝑡𝑡𝑡𝑡𝛽𝛽2 )

Ce cas particulier indique que Ψ = 𝑓𝑓 (Φ, 𝛽𝛽2 ) , Φ étant le


coefficient de charge, Ψ le coefficient de débit Φ et 𝛽𝛽2 l’angle en
sortie (condition géométrique)
Dans le développement de ce modèle le fluide était considéré
sans viscosité, mais dans un contexte plus large il faudra en tenir
compte. On peut écrire alors:
Ψ =f (Φ Re Géo)

gH  Q ρ ND 2 li 
= f , , (design) 
2
N D 2  ND
3
µ D 

Une formulation similaire peut être imaginée pour le rendement

 Q ρ ND 2 li 
η = g , , (design) 
 ND
3
µ D 

À partir de quelle vitesse on ne peut plus considérer l’écoulement


d’un fluide compressible comme étant incompressible?
Le consensus établie indique qu’un écoulement dont le nombre
de Mach satisfait 𝑴𝑴 < 𝟎𝟎. 𝟑𝟑 , peur être traité comme étant
incompressible
Ce critère, peut être expliqué si l’on regarde la relation entre la
pression statique et la pression totale (d’arrêt ou de stagnation):
𝑘𝑘 ⁄𝑘𝑘−1
𝑝𝑝0 𝑘𝑘 − 1 𝑘𝑘 = 𝛾𝛾
= 1+ 𝑀𝑀𝑎𝑎2
𝑝𝑝 2
22

Série
n(n − 1) 2
(1 + x) =+
n
1 nx + x + O( x3 )
2! k k −1
n= x= Ma 2
k / k −1 k −1 2
p0  k − 1 2
= 1 + Ma 
p  2  k
nx = Ma 2
2
p0 k k
=
1 + Ma 2 + Ma 4 + O( Ma 6 )
p 2 8

23

Pression
1 2 Ma 2 
=
p0 − p pkMa 1 + + O( Ma 4 ) 
2  4 

1 1 V 2 1 pk 2 1
=
pkMa 2
=pk 2 = V ρV 2
2 2 a 2 kRT 2

Pour Ma < 0.3


2 4 
2
1 Ma 1
p=
0 − p ρV 1 + + O( Ma )  p0 − p ≈ ρV 2
2  4  2
Pression dynamique

24

L’écoulement dans un ventilateur est à faible nombre de Mach


𝜌𝜌 ≈ 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐. → 𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖. . Le coefficient de débit Φ sera ainsi basé
sur le débit volumique 𝑄𝑄, et le coefficient de charge Ψ sur la
variation de pression ∆𝑃𝑃 entre l’amont et l’aval du rotor

𝑄𝑄 𝛥𝛥𝛥𝛥⁄𝜌𝜌
Φ= Ψ= 2 2
𝑁𝑁𝐷𝐷 3 𝑁𝑁 𝐷𝐷

25
Ψ =f ( Φ Re Géo)

∆P  Q ρ ND 2 li 
= f , , 
ρN D
2 2
 ND
3
µ D

et pour le rendement
 Q ρ ND 2 li 
η = g , , 
 ND
3
µ D
Le débit d’air circulant par un ventilateur qui opère à 800 rpm est de 425
m3 /min. L’augmentation de la pression statique est de 7.6 cm d’H2O,
tandis que et celle de la pression totale est de 10 cm d’H2O. Le
rendement total-à-total est de 75%. Les conditions de stagnation (d’arrêt)
à l’entré sont: T01= 200C et P01 =1bar

On dispose d’un deuxième ventilateur, géométriquement similaire, dont


sa grandeur est ½ fois celle du premier. La vitesse de rotation de ce
ventilateur est de 1000 rpm et il opère sur un point de similitude dit
homologue. Sachant que les conditions thermodynamiques à l’entrée
son les mêmes pour les deux cas, on doit trouver:

Le débit d’air, l’augmentation de la pression statique et la variation de


pression totale dans le ventilateur à échelle réduite
 Q   Q  ➁
=φ = 3   3 
 N D 1  N D 2
3
 D2   N 2 
Q2 = Q1    
➀  D1   N1 
 ∆p0   ∆p0 
=ψ = 2   2 
 ρN 2
D 1  ρN 2
D 2
2 2
 ρ 2   D2   N 2 
∆p02 =
∆p01      
ρ
 1  1   1 
D N
 ∆ ( p + ρV 2 / 2)   ∆ ( p + ρV 2 / 2) 
  =  (faible nombre de Mach)
 ρ N 2
D 2
2  ρ N 2
D 2
1
 ∆p   ∆(V 2 / 2)   ∆p   ∆(V 2 / 2) 
 2 
+  = 2 
+ 
 ρ N D  2  N D  2  ρ N D 1  N D 1
2 2 2 2 2 2

 ∆(V 2 / 2)   ∆(V 2 / 2) 
 2 2  = 2 2  Similitude cinématique
 N D  2  N D 1

2 2
 ρ  D   N  ∆p Q
∆p2 =∆p1  2   2   2  W = 0
 ρ1   D1   N1  η
Dans un compresseur ou dans une turbine on est en présence de
l’écoulement d’un gaz à haute vitesse. On doit alors tenir compte
de la compressibilité qui se manifeste dans le nombre de Mach
Ce paramètre est ainsi partie intégrante des relations de
similitude
La nature du gaz devra également être incluse. Ceci sera pris en
compte par le rapport des chaleurs spécifiques 𝜸𝜸
Le visage des paramètres adimensionnels peut être varié, mais
certaines formes sont plus pratiques que d’autres
le débit

η, ∆h0 s , W= f ( m , ρ01 , T01 , γ, N , µ, D ) le fluide
l’état à l’entrée
la travail à échanger
la grandeur
la vitesse de rotation

We ∆h0 m W
η, =
ψ = φ= P = φψ =
( ND ) ρ01 N D 3 ρ01 N 3 D 5
2
U2
charge débit puissance

∆h0 W  m ρ01 ND 2 ND 
η, , = f , , ,γ
( )
2
ρ N 3 5
D  ρ01 ND 3 µ RT 
ND 01  01 
η Ψ P Φ Re M γ
W
P=
ρ01 N 3 D 5
Paramètre 𝛾𝛾𝑅𝑅
adimensionnel 1/Φ
𝛾𝛾 − 1

𝑾𝑾̇ 𝜌𝜌01 𝑁𝑁𝐷𝐷3 𝑊𝑊̇ 𝜌𝜌01 𝑁𝑁𝐷𝐷3 𝑚𝑚𝑐𝑐̇ 𝑝𝑝 𝛥𝛥𝑇𝑇0 𝑐𝑐𝑝𝑝 𝛥𝛥𝑇𝑇0
⇒ ⋅ = ⋅ = 2 2
𝟑𝟑
𝝆𝝆𝟎𝟎𝟎𝟎 𝑵𝑵 𝑫𝑫𝟓𝟓 𝑚𝑚̇ 3
𝜌𝜌01 𝑁𝑁 𝐷𝐷 5 𝑚𝑚̇ 𝜌𝜌01 𝑁𝑁 3 𝐷𝐷 5 𝑁𝑁 𝐷𝐷

𝛾𝛾 𝑅𝑅𝑇𝑇01 𝛥𝛥𝑇𝑇0 𝜟𝜟𝑻𝑻𝟎𝟎


= ⋅ 2
⋅ →
𝛾𝛾 − 1 𝑁𝑁𝑁𝑁 𝑇𝑇01 𝑻𝑻𝟎𝟎𝟎𝟎

~1/(Mach)2

À la place de W/(ρ01N3D5) on peut considérer ∆T0/T01 !


m
φ=
ρ01 N D 3

m 𝑚𝑚𝑅𝑅𝑇𝑇
̇ 01 𝑚𝑚̇ 𝑅𝑅𝑇𝑇01 𝑅𝑅𝑇𝑇01 𝑚𝑚̇ 𝑅𝑅𝑇𝑇01
= = 2
⋅ 2

ρ 01 ND 3
𝑝𝑝01 𝑁𝑁𝐷𝐷 3 𝑝𝑝01 𝐷𝐷 𝑁𝑁𝑁𝑁 𝑝𝑝01 𝐷𝐷 2

p01=ρ01RT01 ~(1/Mach)

𝒎𝒎̇ 𝒎𝒎̇ 𝑹𝑹𝑹𝑹𝟎𝟎𝟎𝟎


Au lieu de 𝟑𝟑 on peut considérer
𝝆𝝆𝟎𝟎𝟎𝟎𝜨𝜨𝑫𝑫 𝒑𝒑𝟎𝟎𝟎𝟎𝑫𝑫𝟐𝟐

Aussi, ∆𝒉𝒉 = 𝒇𝒇 ∆𝑻𝑻 = 𝒇𝒇(𝒑𝒑𝟎𝟎𝟎𝟎 ⁄𝒑𝒑𝟎𝟎𝟎𝟎 ), alors 𝑾𝑾̇ 𝒑𝒑𝟎𝟎𝟐𝟐


𝟑𝟑 𝟓𝟓 →
𝝆𝝆𝟎𝟎𝟎𝟎 𝑵𝑵 𝑫𝑫 𝒑𝒑𝟎𝟎𝟎𝟎
p02  m RT01 ρ01 ND 2 ND 
η, f , , ,γ
p01  P01D 2
µ RT 
 01 

η Ψ Φ Re M γ

P
Dans le contexte des turbomachines on identifie 7 paramètres
adimensionnels

Ψ = f( Φ, Re, M, γ, Géo)

η = g( Φ, Re, M, γ, Géo)

En fonction du type de machine, certains seront appropriés et


d’autres pourront être négligés
η, Ψ Φ Re M γ
p02  m RT01 ρ01 ND 2 ND 
η, f , , ,γ
p01  P01D 2
µ RT 
 01 
– Reynolds élévé (turbulent)

Simplifications – Même fluide (γ=cnste.)

– Même machine D=cnste.


η, Ψ Φ Re M γ
p02  m RT01 ρ01 ND 2 ND 
η, f , , ,γ
p01  p01D 2
µ RT 
 01 
– Reynolds élévé (turbulent): au-delà d’un nombre de Reynolds
limite, la performance ne varie plus (similitude restreinte)

– Même fluide (γ=cnste)

– Même machine D=cnste


Même machine, même fluide, écoulement turbulent

p02  m T01 N 
η, =
f , 
p01 
 p01 T01 

η Ψ Φ M
 M= f ( Ψ , Φ )
F (η , Φ, Ψ , M ) =0 
 η = g ( Ψ, Φ )
L’évaluation de la performance d’un compresseur est caractérisée
au travers de l’ensemble de paramètres 𝜼𝜼, 𝚿𝚿, 𝚽𝚽, 𝐌𝐌.
Le champ d’opération du compresseur est représenté par des
courbes de nombre de Mach 𝐌𝐌 constant dans le plan 𝚿𝚿, 𝚽𝚽. Il est
d’usage de superposer dans carte du compresseur, des contours
d’isorendement 𝜼𝜼
Selon les simplifications présentées, en pratique on utilise le taux
de compression 𝒑𝒑𝟎𝟎𝟎𝟎 ⁄𝒑𝒑𝟎𝟎𝟎𝟎 au lieu d’une autre forme pour le
coefficient 𝚿𝚿 et le débit corrigé 𝒎𝒎̇ 𝑻𝑻𝟎𝟎𝟎𝟎 �𝒑𝒑𝟎𝟎𝟎𝟎 à la place de 𝚽𝚽
Aussi, le nombre de Mach 𝑴𝑴 est substitué par la vitesse de
rotation corrigée 𝑵𝑵⁄ 𝑻𝑻𝟎𝟎𝟎𝟎
Dans cette carte on distingue deux limites qui contraignent la
plage d’opération d’un compresseur: le pompage et le blocage
sonique
Le pompage: il s’agit d’un phénomène qui a lieu des faibles
débits. Il est associé à un chute soudaine du taux de compression
suivie par des fortes oscillations de la pression, parfois
accompagnée d’une inversion de l’écoulement. Sous cette
condition le compresseur est inopérable

Le blocage sonique: ce phénomène à lieu des forts débits. Il


s’agit de la l’apparition d’une section sonique dans l’étage qui
limite le débit massique maximal. Dans cette limite, la
performance du compresseur chute rapidement
p02
p01

Courbe de

Rapport de Pression
vitesse
corrigée
constante

𝑵𝑵
𝑻𝑻𝟎𝟎

Débit m T0
p0
η
Courbes de
rendement

p02
p01

Rapport de Pression

Débit m T0
p0
Courbes de
rendement

p02
p01
Courbes de

Rapport de Pression
vitesse
angulaire
constante

𝑵𝑵
𝑻𝑻𝟎𝟎

Debit
m T0
p0
Sur la carte d’un compresseur on trace également la ligne
d’opération ou ligne de fonctionnement du compresseur
Elle correspond à la ligne reliant chacun des points de rendement
maximal sur chaque isovitesse
Chaque compresseur a un point d’opération optimal sur la ligne
d’opération
Cette position dans la carte c’est le point de design où le
compresseur fonctionnera la plupart du temps
Ligne d’opération

p02
p01 Point de
design idéal

Rapport de Pression

𝛈𝛈

𝑵𝑵
𝑻𝑻𝟎𝟎

Debit
m T0
p0
Un compresseur centrifuge opère au point nominal. Le rotor a un diamètre
de 40 cm et le rapport de pression totale est p02/p01=2. Les conditions à
l’entrée du compresseur sont T01 =200C P01=1 bar. Utilisez la carte du
compresseur et déterminez:

Le débit massique

La puissance requise

La vitesse angulaire
Un compresseur centrifuge opère au point nominal. Le rotor a un diamètre de 40 cm et le rapport
de pression totale est p02/p01=2. Les conditions à l’entrée du compresseur sont T01 =200C p01=1
bar. Utilisez la carte du compresseur et déterminez:
3.4

Pp2002
ηc = 0.72
Pp1001 0.76
3.0
0. 70
0. 78

0.60
2.6
0.80
0. 50

2.2
𝜼𝜼 = 𝟎𝟎. 𝟖𝟖𝟖𝟖 1.2
0.81
𝟐𝟐 1.1

1.8 1.0

0.9

0.8
𝟎𝟎. 𝟖𝟖𝟖𝟖
1.4
uu
0.7
γRT10
γRT01
1.0
0. 02
𝟎𝟎.0.04
𝟎𝟎𝟎𝟎 0.06 0. 08
m T
0.10
m RT1010
p01DD22
P
m RT01 u
= =
0.04, η 0.81, = 0.85 Le débit massique
γRT01
2
p01 D La puissance requise
La vitesse angulaire

D= 40 cm , p02/p01=2, T01 =200C, p01=1 bar

p01 D 2
m = 0.04
m RT01

105 × (0.4) 2
m =
0.04 × =
2.21 kg / s
293 × (8314 / 28.97)
Opère au point nominal. D= 40 cm, p02/p01=2,T01 =200C, P01=1 bar.
γ−1 γ
 p02  Le débit massique
= T01  = = ( 2)
0.4 1.4
T02 s  T02 s 293 357 K La puissance requise

 p01  La vitesse angulaire

T02 s − T01 m γR
η= =
T02 − T W (T02 − T 01 )
T02 − T01
01
( γ − 1)

DN / 2
= 0.85 DN N=DN/D
γRT01
Remarque:Pour un rapport de pression autre que celui du point nominal, le débit cherché sera
celui correspondant à un point sur la ligne d’opération
Ouverture variable
On constate que dans les turbines hydrauliques on retrouve des
parties mobiles pour ajuster le point d’opération.
On devra alors prêter une attention particulière à la similitude
géométrique
D’autre part, on note que le nombre de Mach 𝐌𝐌 et le rapport des
chaleurs spécifiques 𝜸𝜸 ne sont pas pertinents pour un
écoulement d’eau
Le champ d’opération d’une turbine est représentée par des
courbes d’isorendement dans le plan 𝚿𝚿, 𝚽𝚽 dont l’ensemble est
appelé ‘‘la colline de rendement’’
À cette représentation on superpose plusieurs courbes
correspondant à des positions d’organes mobiles
En pratique industrielle on utilise des formes particulières
associées au système d’unités pour les coefficient 𝚿𝚿 et 𝚽𝚽
Ψ = f( Φ, Re, M,
X Xγ, Géo)

η = g( Φ, Re, M, γ,
X X Géo)

Ψ = f( Φ, Re, Géo)

η = g( Φ, Re, Géo)
Ψ Φ Re Géo
gH  Q ρ ND li 
2
= f , , 
2
N D 2
 ND
3
µ D

 Q ρND 2 li 
η= g , , 
 ND
3
µ D
Dans l’univers des turbines hydrauliques l’intérêt est porté sur la
puissance produite 𝑾𝑾̇ = 𝜼𝜼𝜼𝜼𝒈𝒈𝒈𝒈𝒈𝒈
Il est alors possible d’exprimer le débit 𝑸𝑸 en fonction de la
puissance 𝑾𝑾̇ et de trouver un coefficient de débit 𝚽𝚽 avec un
visage différent
À ce stade, ce n’est pas indispensable de suivre de près les
étapes de ce développement qui est présenté par la suite
Ψ Φ Re Géo
gH  Q ρ ND li 
2
= f , , 
2
N D 2
 ND
3
µ D
 W 
Q= 
 ηρ gH 
 Q ρND 2 li 
η= g , , 
 ND
3
µ D
 W  gH  Q ρND 2 li 
Q=  =ψ = f , , 
 ρ gH  N 2 D2 µ
3
 ND D

Φ Re Géo
gH  W 1 ρND 2 li 
=ψ = f , , 
N 2 D2  ρ µ
3
gH ND D
Q
gH  W 1 ρND 2 li 
=ψ = f , , 
N 2 D2  ρgH ND µ
3
D
écoulement turbulent

gH  W ( gH )1/ 2 ρND 2 li 
=ψ = f , , 
N 2 D2  ρD ( gH ) µ
2 3/ 2
ND D
ND gH
est l ' inverse de
( gH )1/2 N 2 D2
gH  W ( gH )1/ 2 li 
= f , 
 ρ D ( gH )
2 2 2 3/ 2
N D ND D
Il existe déjà

𝑔𝑔𝑔𝑔 𝑊𝑊̇ 𝑙𝑙𝑖𝑖


2 2
= 𝑓𝑓 ,
3⁄2 𝐷𝐷
𝑁𝑁 𝐷𝐷 𝜌𝜌𝐷𝐷 2 𝑔𝑔𝑔𝑔

𝑊𝑊̇ 𝑙𝑙𝑖𝑖
𝜂𝜂 = 𝑔𝑔 ,
3⁄2 𝐷𝐷
𝜌𝜌𝐷𝐷 2 𝑔𝑔𝑔𝑔
Les paramètres adimensionnels qui caractérisent l’opération
d’une turbine hydraulique sont:

𝑔𝑔𝑔𝑔 𝑊𝑊̇ 𝑙𝑙𝑖𝑖


2 2
= 𝑓𝑓 ,
3⁄2 𝐷𝐷
𝑁𝑁 𝐷𝐷 𝜌𝜌𝐷𝐷 2 𝑔𝑔𝑔𝑔
𝑊𝑊̇ 𝑙𝑙𝑖𝑖
𝜂𝜂 = 𝑔𝑔 ,
3⁄2 𝐷𝐷
𝜌𝜌𝐷𝐷2 𝑔𝑔𝑔𝑔
η, Φ, Ψ, Géo
.
W ND li
η, , ,
ρD ( gH )
2 3/2
( gH )1/2
D
Formulation explicite
Formulation implicite
Géo = f( Φ, Ψ)
F (η , Φ, Ψ , Géo ) =0
η = g( Φ, Ψ )
Courbes
α d’ouverture

Coeff. de charge  gH Q 
α= f  2 2 , 3 
 n D nD 

α ψ Φ

Coeff .de débit


Φ
Courbes de
rendement

η
Ψ

Coeff. de charge  gH Q 
η= g  2 2 , 3 
 n D nD 

η ψ Φ

Coeff .de débit


Φ
Courbes Courbes de

α d’ouverture rendement

η
Ψ

Coeff. de charge

Coeff .de débit


Φ
En pratique industrielle, des variables réduites sont définies pour
la construction de cartes destinées au passage des données du
modèle vers le prototype

Ces grandeurs, notées avec un double indice 1, correspondent à


un fonctionnement en similitude sous une chute H de 1 m (1pi) et
un rotor de diamètre D de 1 m (1pi)
Vitesse angulaire réduite 𝑵𝑵𝟏𝟏𝟏𝟏
2
𝑁𝑁 2 × 𝐷𝐷 2 𝑁𝑁11 × 𝟏𝟏𝟐𝟐 𝑁𝑁𝑁𝑁
Ψ= = 𝑁𝑁11 = 1⁄2
𝑔𝑔 × 𝐻𝐻 𝑔𝑔 × 𝟏𝟏 𝐻𝐻

Débit réduit 𝑸𝑸𝟏𝟏𝟏𝟏

𝑄𝑄 𝑄𝑄11 𝑄𝑄
Φ= = 𝑄𝑄11 = 2 1⁄2
𝑁𝑁 × 𝐷𝐷 3 𝑁𝑁11 × 𝟏𝟏𝟑𝟑 𝐷𝐷 𝐻𝐻
Puissance réduite P11

𝑊𝑊̇
𝑃𝑃11 = 2 3⁄2
𝐷𝐷 𝐻𝐻

On note que les quantités N11, Q11 et P11 ne sont pas


adimensionnelles. Alors, leurs grandeurs dépendent du système
d'unités utilisé!
La vitesse spécifique et le diamètre spécifique, sont deux
concepts issues de l'étude des lois de similitude des
turbomachines. Ils sont très utiles pour la conception de
machines appartenant à des familles similaires (homothétiques)
et constituent une base pour le classement de chaque appareil
.
Avec ces chiffres adimensionnels, on peut déterminer la
turbomachine la plus adéquate pour une application donnée
Remarques:
La vitesse spécifique fournit par le milieu industriel se réfère au
point d'efficacité maximum.

Roue dite lente : Vitesse spécifique faible


Roue dite rapide Vitesse spécifique élevée

Cependant, les roues lentes tournent généralement vite tandis les


roues rapides tournent généralement lentement
Les deux nouvelles quantités adimensionnelles, vitesse
spécifique et diamètre spécifique, associées à la vitesse et au
diamètre de la machine, respectivement, peuvent être obtenues
par la combinaison des coefficients de débit et de charge

Bien qu’il s’agit de notions générales, elle seront présentées


dans le cadre de machines opérant avec des fluides
incompressibles
 Q   Q 
φ = 3 
N = 3 
 N D  φD  Q /φ Q1/ 2 Ψ1/ 4
D =
2
D=
(We / Ψ ) (We ) Φ1/ 2
1/ 2 1/ 4

1/ 2
 W   W 
ψ = 2 e 2  N = e2 
N D   ΨD 
𝐷𝐷 𝑊𝑊𝑒𝑒 1⁄4 Ψ 1⁄4
𝐷𝐷𝑠𝑠 = 1⁄2
= 1⁄2
𝑄𝑄 Φ
Diamètre spécifique
1/3
 Q   Q 
φ = 3  D= 
 N D   Nφ D  (We )
3/ 4
Φ1/ 2
N=
1/2 Q1/ 2 Ψ 3/ 4
 W   W 
ψ = 2 e 2  D= e2 
N D   ΨN 
𝑁𝑁𝑄𝑄1⁄2 Φ1⁄2
𝑁𝑁𝑠𝑠 = 3⁄4
= 3⁄4
𝑊𝑊𝑒𝑒 Ψ
Vitesse spécifique
Φ1⁄2 Ψ 1⁄4
𝑁𝑁𝑠𝑠 = 3⁄4 𝐷𝐷𝑠𝑠 = 1⁄2
Ψ Φ

Les quantités adimensionnelles, vitesse et diamètre spécifique 𝑁𝑁𝑠𝑠


et 𝐷𝐷𝑠𝑠 , respectivement, permettent de caractériser les
turbomachines par des variables dans lesquelles n’apparaissent
que le travail (spécifique), le débit et vitesse de rotation
On retrouve des formes particulières pour les divers types de
turbomachines. Certaines, même ayant des unités
 φ 1/ 2  1/ 2
N s =  3 / 4  NQ 𝑁𝑁 𝑄𝑄
ψ 
 Q 
Ns = 3/ 4
𝑛𝑛𝑞𝑞 = 3⁄4 3⁄4
𝑔𝑔 𝐻𝐻
φ = 
 N D

3 

(gH )
gH
ψ= Vitesse spécifique adimensionnelle Vitesse dimensionnelle
N 2D2
 ψ 1/ 4 
Ds =  1/ 2  D( gH )1/ 4
φ 
Ds = 1/ 2
Q
Diamètre spécifique adimensionnel
NQ1/ 2 •

Ns = Qeq =
W
(gH ) 3 / 4 ρ gH

ρ, g
  •  1/ 2

N W  
  ρ gH   𝑁𝑁𝑊𝑊̇ 1/2
N s =     𝑛𝑛𝑠𝑠 =
( gH )3/ 4  𝐻𝐻5⁄4
 
 
  Vitesse dimensionnelle
 
Vitesse spécifique adimensionnelle
Les deux formes de vitesse spécifique, 𝒏𝒏𝒔𝒔 et 𝒏𝒏𝒒𝒒 , peuvent
également être obtenues de manière pratique en comparant la
machine à caractériser avec un machine d’étalonnage
Cette idée sera présentée dans le cadre d’une pompe
 Q  gH W
φ = 3 
ψ= 2 2 P = φψ =
 ND  N D ρN 3 D5

Pompe de référence Pompe d’étalonnage

H 0 , Q0 W0 H = 1m
W = 1CV
N0

2 2
H  N   D
ψ =      W   H 0  𝑵𝑵𝟎𝟎 𝑾𝑾̇ 𝟎𝟎
1/2 5/4
H0  0  0
N D N0 = N     𝒏𝒏𝒔𝒔 = ⁄𝟒𝟒
 0
W  H  𝑯𝑯𝟓𝟓𝟎𝟎
W
3 5
 N   D
𝑃𝑃� =    
W 0  0  0
N D
𝑯𝑯 = 𝟏𝟏𝒎𝒎, 𝑾𝑾̇ = 𝟏𝟏𝑪𝑪𝑪𝑪, 𝑵𝑵 = 𝒏𝒏𝒔𝒔
 Q  gH W
φ = 3 
ψ= 2 2 P = φψ =
 ND  N D ρN 3 D5
Pompe de référence Pompe d’étalonnage

H 0 , Q0 W0 H = 1m
W = 1CV
N0

2 2
𝐻𝐻 𝑁𝑁 𝐷𝐷
ψ =
𝐻𝐻0 𝑁𝑁0 𝐷𝐷0 1/2
 Q   H0 
3/4 𝑵𝑵𝟎𝟎 𝑸𝑸𝟎𝟎
N0 = N 𝒏𝒏𝒒𝒒 =
   𝑯𝑯𝟑𝟑𝟎𝟎
⁄𝟒𝟒
 Q0   H 
3
ϕ  Q   D  N 
 =   
 0   D0   N 0 
Q 𝑯𝑯 = 𝟏𝟏𝒎𝒎, 𝑸𝑸 = 𝟏𝟏 𝒎𝒎𝟑𝟑 ⁄𝒔𝒔 , 𝑵𝑵 = 𝒏𝒏𝒒𝒒
L’interprétation de la vitesse spécifique sur la base d’une
machine d’étalonnage, dont la chute, le débit ou la puissance sont
unitaires, permet de simplifier virtuellement des unités et
d’imaginer ainsi une vitesse spécifique en rpm.
Cette approche est pratique en génie, mais pour des systèmes
d’unités différents, la valeur numérique de la vitesse
spécifique sera différente. En d’autres mots, il y aura des rpm
associées au système anglais et des rpm jumelées au système
international
𝑁𝑁 𝑄𝑄
𝑛𝑛𝑞𝑞 = 3⁄4
𝐻𝐻

nq correspond à la vitesse spécifique d’une machine, en rpm, qui


est géométriquement similaire à une machine d’étalonnage
opérant avec une tête de H = 1 m (1pi) et par laquelle circule un
débit Q = 1 m3/s (1 pi3/s)
𝑁𝑁𝑊𝑊̇ 1/2
𝑛𝑛𝑠𝑠 =
𝐻𝐻5⁄4

ns correspond à la vitesse spécifique d’une machine, en rpm, qui


est géométriquement similaire à une machine d’étalonnage
opérant avec une tête de H = 1 m (1pi) et qui consomme (produit)
une puissance de 𝑾𝑾̇ = 𝟏𝟏𝟏𝟏𝟏𝟏(𝟏𝟏, 𝒌𝒌𝒌𝒌, 𝟏𝟏 𝑯𝑯𝑯𝑯)
ns = 100

ns = 35.5
ns = 130 ns = 1200

On remarque que la vitesse spécifique caractérise le type d’aubage:


axial radial ou mixte, qu’une roue devra posséder en fonction de la
charge, le débit et la vitesse de rotation
nq = n11 Q11

nq = NQ1/2 H −3/4

ns = NW 1/2 H −5/4

ns = 3.65nq Le coefficient 3.65 n’est valable que pour des unités


métriques avec la puissance mesurée en chevaux-vapeur.
Bâche spirale
Aubes directrices

Aspirateur Rotor
Francis
0.5

22.5

93%

92 17.5
90
W
P11 = 2 3 / 2
P11

88
86 D H
15

84 12.5 ND
82 N 11 =
78 80 10 H 1/ 2

5 7.5

0
100 N11 275
za
• dd : diamètre de la buse
• ds : diamètre du jet
• D : diamètre de référence
• za : nombre d’augets
• zi : nombre d’injecteurs

zi
Pelton
W
P11 =
D 2 H 3/ 2
Kaplan
P11
89% 87 85 83 81

W
P11 = 2 3 / 2
D H

α=200 ND
N 11 =
H 1/ 2
α=100

92% 91 89 87 85 83 81

N11
Remarque: on associe une colline à chaque position des aubes mobiles
Les turbines hydrauliques, peuvent être classifiées selon la valeur
de la vitesse spécifique:

𝑁𝑁𝑊𝑊̇ 1/2
𝑛𝑛𝑠𝑠 =
𝐻𝐻5⁄4
Turbine adéquate pour des grandes chutes H et des faibles débits
Q. Les valeurs de ns sont faibles
Turbine adéquate pour des grandes chutes et des débits moyens
Les valeurs de ns sont moyens
Turbine adéquate pour des faibles chutes et des débits élévés. Elle
opère à des ns élévées
ns

Pelton Francis Kaplan

• La géométrie des turbines varie en fonction de ns

• Au fur et à mesure que ns augmente la forme de ces machines


change de radiale vers axiale. Ceci se produit lorsque le débit Q
augmente et la chute H diminue.
L’évolution de la géométrie découverte pour les turbines
hydrauliques en fonction de la vitesse spécifique, est aussi
valable pour l’ensemble des turbomachines. Toutes les
turbomachines, peuvent alors être classifiées selon cette quantité
adimensionnelle
En particulier, on a trouvé que les machines se transforment de
radiales à axiales au fur et mesure que la vitesse spécifique
augmente
La nature radiale-axiale d’un rotor dépende de la vitesse spécifique
Q
nq = N 3
4
H
SI 20 40 60 80 100 150 200 300

10000

15000

20000
1000

1500

2000

3000

4000

5000
6000
500
600
700
800
900

US

Rotor radial Rotor Francis Rotor mixte Rotors axial


Vitesse spécifique nq
Q
nq = N 3
4
H

nq (US gpm, pi) = 0.861 nq (m3/h, m)


nq (Imp. gpm, pi) = 0.787 nq (m3/h, m)
avec
1 US gpm = 0.2271 m3/h
1 Imp. gpm = 0.2728 m3/h
1 pi = 0.3048 m
Rotor Pelton Rotor Pelton Turbine Francis Turbine à vapeur Pompe centrifuge Compresseur Centrifuge

ns
Turbine à vapeur Turbine à vapeur Pompe centrifuge Turbine à gaz Compresseur radial Compresseur axial

ns
Turbine à vapeur Turbine Francis Turbine Francis Pompe Mixte Turbine Kaplan Compresseur axial

ns
Compresseur axial Pompe hélice Soufflante axiale Pompe hélice Turbine Kaplan

ns
On propose la construction d’une turbine de type Pelton ayant les mêmes
caractéristiques que celles d’un design existant. Les paramètres de vitesse
et de puissance sont données par une carte de rendement. Sur l’axe des
abscisses on trouve les paramètres 𝑵𝑵𝟏𝟏𝟏𝟏 = 𝑵𝑵𝑵𝑵⁄𝑯𝑯𝟏𝟏/𝟐𝟐 , tandis que sur l’axe
des ordonnées on trouve le coefficient 𝑷𝑷𝟏𝟏𝟏𝟏 = 𝑾𝑾̇ ⁄(𝑯𝑯𝟑𝟑/𝟐𝟐 𝑫𝑫𝟐𝟐 ) Le rendement η
et la vitesse spécifique dimensionnelle dans le système métrique sont
représentés par des isocontours. La charge ou chute nette pour
l’aménagement hydroélectrique est de 𝑯𝑯 = 𝟑𝟑𝟑𝟑𝟑𝟑𝟑𝟑 et la puissance produite
est 𝑾𝑾̇ = 𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐.

Considérez un seul injecteur et sur la base du point de design (le point


nominal), déterminez: la vitesse de rotation, le diamètre du jet de
l’injecteur et le diamètre de la roue.
patm
réservoir

H
V j ≈ 2 gH

z=0
U
W2 = W1

C1 U W=
1 C1 − U

C1

U
C1 • C1 : vitesse absolue du jet à l’entrée
• U : vitesse tangentielle de la roue
• (C1-U): vitesse relative à l’entrée
• (C1-U): vitesse relative à la sortie (pas des pertes)
𝒏𝒏𝒔𝒔 =30
Remarque: ns a été calculée avec
NW 1/2
la puissance en CV, la vitesse en ns =
rpm et la hauteur en mètres H 5/4

H=300m
̇ 20 000 kW
𝑾𝑾=
η =0.9

ND
=η 0=
.9, ns 30, 1/ 2
= 40
H
…sur la base du point de design (le point
nominal), déterminez…

40
Un seul injecteur et sur la base du point de design (le point nominal), déterminez:
la vitesse de rotation, le diamètre du jet de l’injecteur et le diamètre de la roue.

Remarque: ns a été calculée avec


ND NW 1/2
la puissance en CV, la vitesse en
=η 0=
.9, ns 30, = 40 ns =
rpm et la hauteur en mètres
H 1/ 2 H 5/4
H=300m
̇ 20 000 kW
𝑾𝑾=
Facteur de conversion

N × W 1/2 (en CV ) 𝑵𝑵 × 20000 × 1.359 1⁄2



s n= =
30 = N=227 rpm
s
H 5/4
300 5⁄4

V ≈ 2 gH = 2 × 9.81× 300 W = ηρQgH Q


d: diamètre du jet ? Q = VA Aire du jet ?
Un seul injecteur et sur la base du point de design (le point nominal), déterminez:
la vitesse de rotation, le diamètre du jet de l’injecteur et le diamètre de la roue.

� 𝒔𝒔 =ns a été calculée


Remarque: 𝑵𝑵
ND
avec la puissance en CV, la
vitesse en rpm et la hauteur en
=η 0=
.9, ns 30, 1/ 2
= 40
mètres H
H=300m
̇ 20 000 kW
𝑾𝑾=

Q Q = VA A = πd2 / 4 C1

d d: diamètre du jet

ND
N 1/ 2
= 40 D
H
diamètre de la roue
Une turbine Francis (simple), opère avec une charge H de 60 m et un
débit Q de 30 m3/s. Le rendement est η=88%.Utilisez la carte et estimez
les rpm et le diamètre D de la roue.

Remarque: Pour le calcul de ns la puissance est en HP


Supposez Hs=+3m
1000
ns 900
800 Hs =+ 1m KAPLAN
700
Hs = -0.5m
600
500 Hs = -2 m FRANCIS

H = 60m,
Hs = -4 m DOUBLES
400

300
Q = 30m / s,
3
260
Hs =3 m
Hs =5 m

η = 0.88 200
FRANCIS
SIMPLES

Hs =+5 m

100 Hs =+3 m
90

NW 1/2
80
70
ns = 5/4 60 6 jets

H 50
4 jets
40

3 jets
30
PELTON 2 jets
20
1 jet

10
2 4 6 8 10 20 40 60 100 200 400 800 1600
600 1000
60 H
=
H = 60m, =
Q 30 m3 / s, η 0.88, ns = 260

NW 1/2  Estimez les rpm et


ns = 5/4  le diamètre D de la roue.
 ns H 5/4
H
 N = 1/2
ηρgQH  W ( Hp)

W=
745.7 

1 HP=745.7 Watts

𝜋𝜋𝑫𝑫𝑁𝑁
2𝑔𝑔𝑔𝑔 ≈ D
60
Approximation grossière, mais..
Le diagramme de Cordier est une représentation dans le
plan 𝑵𝑵𝒔𝒔 − 𝑫𝑫𝒔𝒔 (diamètre spécifique, vitesse spécifique) qui aide à
choisir le type de machine et fournit une première estimation du
rendement
Ce graphique est très avantageux puisqu’il a été construit à
partir de données expérimentales (machines existantes). Cette
représentation donne de plus une première estimation du
diamètre du rotor
50
40
30

Φ1/2 20

N s = 3/4
Ψ 10
8

Ψ1/4 6

Ds = 1/2 4

Φ 3

2
.
1.0
0.8
0.6

0.4
0.3

0.2

0.1
0.3 0.5 1.0 2 3 4 5 10 20 30
=N s 9.00 Ds−2.103 Ds ≤ 2.8
=N s 3.25 Ds−1.126 Ds ≥ 2.8

=Ds 2.84 N s−0.476 N s ≥ 1.0


=Ds 2.84 N s−0.888 N s ≤ 1.0

 φ 1/ 2   Q 
N s =  3 / 4  φ =  
3
ψ   ND 
 ψ 1/ 4  ψ=
gH
Ds =  1/ 2  N 2D2
φ 
Pour sélectionner une machine pour une tâche donnée, on
considère le débit, la charge et la vitesse de rotation

Ces trois variables permettent de calculer la vitesse spécifique 𝑵𝑵𝒔𝒔

Avec 𝑵𝑵𝒔𝒔 et l’utilisation du diagramme de Cordier, on détermine le


diamètre spécifique 𝑫𝑫𝒔𝒔

Finalement, à partir du 𝐷𝐷𝑠𝑠 , on trouve une valeur du diamètre 𝑫𝑫 de


la machine
NQ1/2
NS = pompe
( gH ) 4
3
N, Q, ΔP, H Ns
1/2
NQ
NS = ventilateur
( ∆P0 / ρ ) 4
3

D ( ∆P0 / ρ )
1
4
DS =
Ns Q1/2
Ds D
D ( gH )
1
4
DS =
Ds Q1/2
Le facteur de forme n’est que le concept nord-américain pour la
vitesse spécifique adimensionnelle, Il est défini pour le point
nominal ηmax , et l’on note par le symbole S
C’est la vitesse spécifique
adimensionnelle

𝜔𝜔 𝑄𝑄
𝑆𝑆 =
𝑔𝑔𝑔𝑔 3⁄4
D
nq S
500 0.18 radial: pression élevée faible débit

radial
1000 0.37

3400 1.25 mixte

6400 mixte
2.33

axial: débit élevé faible pression


10000 3.67

*N en rpm, Q en gpm, H en pi nq = 2732S


D2
𝑁𝑁 𝑄𝑄 𝑟𝑟𝑟𝑟𝑟𝑟
𝑆𝑆 = 𝑁𝑁𝑠𝑠 = 𝑁𝑁
𝑔𝑔𝑔𝑔 3⁄4 𝑠𝑠

D2 D1

D1

Ns basse Ns haute
Pompes

• La plage de fonctionnement est confortable nq faible

• La courbe de H-Q est toujours (presque)


décroissante. H

• La puissance P augmente avec le débit.


P

Q
Pompes

• Lerendement maximal peut être plus élevé nq élevée


que celui des pompes centrifuges, mais la
plage de fonctionnement est plus étroite
H
• La courbe de H-Q est toujours décroissante.
P
• La puissance P décroît avec le débit.

Q
Pompes
= 0.95
Centrifuges: nq = 10 -100
0.90 (ηmax pour nq ≈ 50)
Mixtes: nq = 75 - 200
0.85 (ηmax pour nq ≈ 130)
0.80 Axiales nq = 200 - 320
(ηmax pour nq ≈ 250)
0.75
pompe
pompe centrífuge pompe mixte
axiale
0.70

0.65
0.60
10 15 20 25 30 40 50 60 70 100 150 200 250 300 nq
Q
nq = N 3
Pompes
4
d2
W2 C2 H , P ,η H
C2u

d1 U2 Q
d 2 / d1= 2 − 3.5

d2
W2 C2 H , P ,η

C2u
d1 U2
d 2 / d= 1.5 − 2 Q
1

d2

H , P ,η
W2 C2

C2u
Q
d 2 / d1 = 1 d1 U2
PROBLÈMES
Exemple: H = 70 pi Q = 5.35 pi3/s g=32.2 pi2/s
n = 870 rpm β2 = 22.5 o

Calculer la valeur de la vitesse spécifique (en fonction de Q)


adimensionnelle Ns et dimensionnelle nq (débit Q en gpm)

2𝜋𝜋 𝑄𝑄 𝜋𝜋 5.35
𝑵𝑵𝒔𝒔 = 𝑛𝑛 = 870 = 0.6442
60 𝑔𝑔𝑔𝑔 3⁄4 30 32.2 × 70 3⁄4

2401 3⁄ 1/2
𝑄𝑄 3 𝒓𝒓𝒓𝒓𝒓𝒓 − 𝑝𝑝𝑝𝑝 𝑠𝑠
𝒏𝒏𝒒𝒒 = 𝑛𝑛 = 870 × 3 ⁄ 4
= 1.761210
𝐻𝐻3⁄4 70 𝑝𝑝𝑖𝑖 3⁄4
1 pi3/s=448.8 gpm
Dans cette expression
3⁄ 1/2
𝑄𝑄 3 𝒓𝒓𝒓𝒓𝒓𝒓 − 𝑝𝑝𝑝𝑝 𝑠𝑠
𝒏𝒏𝒒𝒒 = 𝑛𝑛 = 1.761210
𝐻𝐻3⁄4 𝑝𝑝𝑖𝑖 3⁄4
Les caractères en gris pâle ne sont qu’un aide-mémoire pour rappeler
que la valeur numérique de la vitesse spécifique 𝒏𝒏𝒒𝒒 , en 𝒓𝒓𝒓𝒓𝒓𝒓, a été
obtenue avec un débit en 𝒑𝒑𝒑𝒑𝟑𝟑 ⁄𝒔𝒔 (initialement en 𝒈𝒈𝒈𝒈𝒈𝒈) et une
hauteur en 𝒑𝒑𝒑𝒑
NQ 1/ 2
nq =
H 3/ 4 …

Pour une pompe centrifuge, on a les données suivantes:


𝒏𝒏 = 𝟖𝟖𝟖𝟖𝟖𝟖 𝒓𝒓𝒓𝒓𝒓𝒓, 𝐙𝐙 = 𝟔𝟔 pales, 𝜼𝜼𝑯𝑯 = 𝟎𝟎. 𝟗𝟗𝟗𝟗(rendement hydraulique)
𝐐𝐐 = 𝟏𝟏𝟏𝟏 𝟎𝟎𝟎𝟎𝟎𝟎 𝐠𝐠𝐠𝐠𝐠𝐠 (gallons/ minute), 𝑫𝑫𝟐𝟐 = 𝟑𝟑𝟑𝟑 𝒑𝒑𝒑𝒑, (le diamètre en
sortie) 𝛃𝛃𝟐𝟐𝒂𝒂 = 𝟔𝟔𝟔𝟔. 𝟒𝟒𝟎𝟎 , (l’angle à la sortie)
Pour cette pompe, on a trouvé une relation industrielle empirique
𝚽𝚽𝟐𝟐 = 𝒏𝒏𝒒𝒒 /𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏 et on peut supposer que 𝒄𝒄𝟏𝟏𝟏𝟏 = 𝟎𝟎
On vous demande de calculer: 𝐔𝐔𝟐𝟐 , 𝒘𝒘𝟐𝟐𝒎𝒎 , 𝒄𝒄𝟐𝟐𝒖𝒖 , 𝐇𝐇𝒕𝒕𝒕𝒕𝒕𝒕𝒕 , 𝐇𝐇. Tenez
compte du glissement !

5
NQ 1/ 2
nq =
H 3/ 4 …

On vous suggère d’estimer, dans un premier temps, la vitesse


spécifique, 𝒏𝒏𝒒𝒒 = 𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏
Cette hypothèse doit être vérifiée à l’aide de vos résultats

147
U2 ,W2m ,C2u ,Hth, H ?
1 pi =12 po
π D2 n
U2 = 146.7 pi / s
60
Relation industrielle empirique pour Φ2 𝑼𝑼𝟐𝟐
 nq  1000
=
Φ2  =  = 0.082
 15900  15900
c2 m w2 m
Φ=
2 = 𝒄𝒄𝟐𝟐𝒎𝒎 = 𝒘𝒘𝟐𝟐𝒎𝒎 = 146.7 × 0.082 = 𝟏𝟏𝟏𝟏 𝒑𝒑 𝒊𝒊⁄𝒔𝒔
U2 U2
Z =6 D2 = 38 po
c2u = ? Q = 10 000 gpm
β2 = 68.4 o
n = 885 rpm
ηH=0.92 nq = 1000
β2 Z =6
Q = 10 000 gpm
D2 = 38 po
n = 885 rpm
C2 W2 β2 = 68.4 o ηH=0.92
C2m
W2m
C2u W2u
U2

𝑼𝑼𝟐𝟐 = 𝟏𝟏𝟏𝟏𝟏𝟏. 𝟕𝟕 𝒑𝒑 𝒊𝒊⁄𝒔𝒔


W2
𝒘𝒘𝟐𝟐𝒎𝒎 = 𝟏𝟏𝟏𝟏 𝒑𝒑 𝒊𝒊⁄𝒔𝒔
C2
U2
𝒄𝒄𝟐𝟐𝒖𝒖𝒂𝒂 = 𝑈𝑈2 − 𝑊𝑊2𝑢𝑢 = 𝑈𝑈2 − 𝑊𝑊2𝑚𝑚 × tan𝛽𝛽2a = 𝟏𝟏𝟏𝟏𝟏𝟏. 𝟒𝟒 𝒑𝒑 𝒊𝒊⁄𝒔𝒔

𝒄𝒄𝟐𝟐𝐮𝐮
U2 ,W2m ,C2u ,Hth, H ?
𝒄𝒄𝟏𝟏𝒖𝒖 = 𝟎𝟎 𝒄𝒄𝟐𝟐𝒖𝒖 𝑼𝑼𝟐𝟐
𝑐𝑐2𝑢𝑢 𝑈𝑈2 116.4 × 146.7
𝐻𝐻𝑡𝑡𝑡 = = = 531 𝑝𝑝𝑝𝑝 Tenez compte du glissement
𝑔𝑔 32.2
𝒈𝒈(𝒑𝒑𝒑𝒑⁄𝒔𝒔𝟐𝟐 )

cos β 2 a c2uf
σs =
1− =
0.8269 𝛽𝛽2 = 68.40 , σ s =
Z 0.7 c2ua

c2uf = σ s c2ua c2uf = 0.8269 ×116.4 = 96.25 pi / s


𝒄𝒄𝟐𝟐𝒖𝒖𝒖𝒖
U2, W2m, C2u, Hth, H?
𝒄𝒄𝟐𝟐𝒖𝒖𝒖𝒖 𝑼𝑼𝟐𝟐
c2uf U 2 96.25 ×146.7
H ηH
= = 0.92 × = 403.43 𝑝𝑝𝑝𝑝
g 32.2
g
Vérification de la valeur estimée de 𝑛𝑛𝑞𝑞

NQ1/2 885 × (10000)0.5


nq = 3/4
=0.75
983.14 rpm
H (403.43)
Pas trop loin de la valeur proposée nq= 1000
Estimez le diamètre, le type et le rendement d’une pompe
satisfaisant les conditions suivantes
H = 18 pi, Q = 250 gpm, n = 1500 rpm pi 3 s  0.00223 × Q [ gpm ]
Q =

 2π  Q 𝜋𝜋 250 × 𝟎𝟎. 𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎


=Ns n  × = 1500 × = 𝟎𝟎. 𝟗𝟗𝟗𝟗𝟗𝟗
 60  ( gH )
34
30 32.2 × 18 3 ⁄ 4

g adimensionnelle

Un graphique 𝜼𝜼 vs. 𝑫𝑫𝒔𝒔 est disponible


50
40
30

20

10
8
6

4
3

N s = 0.992 1.0
0.8
0.6

0.4
0.3

0.2

0.1
0.3 0.5 1.0 2 3 4 5 10 20 30
DS = 3.1
Diamètre et rendement

D ( gH )
14 3.1 × 0.5575 1⁄2
Ds = 3.1 𝐷𝐷 = ⁄
= 𝟎𝟎. 𝟒𝟒𝟒𝟒𝟒𝟒 𝒑𝒑𝒑𝒑
12 32.2 × 18 1 4
Q
Diagramme disponible
η ≈ 90%

𝑁𝑁 𝑄𝑄 1500 250 𝒓𝒓𝒓𝒓𝒓𝒓 − 𝑔𝑔𝑔𝑔𝑚𝑚1⁄2


𝑛𝑛𝑞𝑞 = 3 = 3� = 𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐
𝐻𝐻 �4 18 4 𝑝𝑝𝑖𝑖 3⁄4

H = 18 pi, Q = 250 gpm, n = 1500 rpm


Q = 0.5575 p3/s
Type de pompe
Q
nq = N 3
4
H

SI 20 40 60 80 100 150 200 300

10000

15000

20000
1000

1500

2000

3000

4000

5000
6000
500
600
700
800
900

US

Rotors radiales Rotors Francis Rotors mixtes Rotors axiales


Vitesse spécifique nq
D2
Pour un compresseur centrifuge on a les données suivantes:
D1S
n = 60 000 rpm D2 = 5.92 po, D1r =1.35 po
D1s =3.84 po T01 = 60 F p01 = 14.7 psia
p02 = 61,74 psia Z =33 pales radiales (β2 = 90 o ) , D1r
m = 2.2 lb/s Tenez compte du glissement

-Calculez le travail spécifique We en utilisant des conditions cinématiques et


géométriques (ne passez pas par la thermodynamique).

-Calculez le travail spécifique Ws idéal. Par la suite, trouvez le rendement η=


Ws/ We.
D2
Pour un compresseur centrifuge on a les données suivantes:
D1S
n = 60 000 rpm D2 = 5.92 po, D1r =1.35 po
D1s =3.84 po T01 = 60 F p01 = 14.7 psia
p02 = 61,74 psia Z =33 pales radiales (β2 = 90 o ) , D1r
m = 2.2 lb/s Tenez compte du glissement

-[Link] le nombre de Mach M1, la température statique T1 ainsi


que la masse volumique ρ1 à l’entrée (1). L’aire correspond à l’anneau limité
par D1r et D1s

-Calculez le débit volumique Q1 , la ‘hauteur’ idéale H=W s/g, la vitesse et le


diamètre spécifique nqet Ds (dimensionnels) et le rendement η. Vérifiez les
calculs à l’aide du tableau suivant. Supposez Q2= Q1
Rendement η
0.4 0.5 0.6 0.7 0.8
50 2.42 2.65 2.91
60 1.94 2.15 2.26

Vitesse spécifique dimensionnelle


65 1.77 1.92 2.02 2.14
70 1.66 1.82 1.89 1.96
80 1.44 1.55 1.63 1.68
85 1.36 1.48 1.53 1.57 1.70
90 1.30 1.39 1.43 1.46 1.59
100 1.16 1.25 1.29 1.32 1.41
110 1.07 1.14 Ds 1.17 1.21 1.29
120 1.00 1.06 1.10 1.15 1.22
130 0.91 1.00 1.03 1.08 1.18
140 0.87 0.96 1.00 1.06
150 0.83 0.94 1.00 1.07
160 0.80 0.91 1.00 1.04
170 0.80 0.91 1.00 1.11
180 0.80 0.91 1.00
Ns 190 0.79. 0.91 1.01
200 0.79 0.91
Calculez le travail spécifique We 1po = 0.0254m
Facteur de conversion

π D2 N π × 5.92 × 0.0254 × 60000


=
U2 = = 472.39 m / s
60 60
2 2
σ = 1− (900 ) σ =−
1 =0.939 (c1u = 0)
z 33

c2=
ua U=2 472.39 m / s (90 0
)

c2uf =σ × c2ua =443.76 m / s

We c=
2 uf U 2 2.0983 × 10 5
m 2
/ s 2
Calculez le travail spécifique Ws idéal. Par la suite, trouvez le rendement η= Ws / We

p02
=
rp = 4.2
p01 Conversion degrés R K

5
T01 = 460 + 60 = 520 R =
T01 ( K ) = T01 ( R) 288.89 K
9

γR
Ws T01 ( rp (γ −1 γ ) − 1) = 1.471×105 m 2 / s 2 plus tard
γ −1
Ws
η = 0.7016
= plus tard
We
[Link] le nombre de Mach M1, à l’entrée, la température statique T1 ainsi que la
masse volumique ρ1 à l’entrée (1). L’aire correspond à l’anneau limité par D1r et D1s

γ +1

RT0  γ −1 2  2( γ −1)
m= M γ 1 + M 
p0 A  2 
π ( D 12s − D 12r ) 1 po = 0.0254m
A=
4
π (3.84 2 −1.35 2 ) × 0.02542
A = 0.0065 m 2
4
[Link] le nombre de Mach M1, à l’entrée, la température statique T1 ainsi que la
masse volumique ρ1 à l’entrée (1). L’aire correspond à l’anneau limité par D1r et D1s

T01 = 288.89 K
γ +1

287 × 288.89  γ −1 2  2( γ −1)
1× − M γ 1 + M  =
0
101320 × 0.0065  2 
p01 = 14.7 psia
m =
2.2 lb / s 1 kg / s
Le nombre de Mach M1, à l’entrée
γ +1

RT0  γ −1 2  2( γ −1)
m − M γ 1 + M  =
0
P0 A  2 

M = 0.4039

f =@(M)m*sqrt(R*T01)/(P01*A) - M*sqrt(gama)*(1+0.5*gama1*M*M)^expop

gama1 = gama-1;
expop = -0.5*(gama+1)/(gama1);
Calculez le débit volumique Q1, la ‘hauteur’ idéale Hs=W s/g, 𝑇𝑇01 (𝐾𝐾) = 288.89𝐾𝐾
𝑝𝑝01 = 101.3 𝑘𝑘𝑘𝑘𝑘𝑘
T01 γ −1 2 T1 = 279.76 K
= 1+ M
T1 2

c1 = M 1 γ RT1 c1 = 135.41 m / s

( γ−1) / γ
 p1   T1  =p1 9.055 × 104 Pa
 = 
 p01   T01 

 p1 
ρ1 =  ρ1 =1.1278 kg / m3
 1
RT
Calculez le débit volumique Q1, la ‘hauteur’ idéale Hs=W s/g,

 m (kg / s )  1 pi = 0.3048 m
Q1 ( pi / s ) =  ( 0.3048)
3 3
3 
ρ
 1 ( kg / m )

 1(kg / s ) 
3  ( )
Q1 ( pi / s ) =
3 3 3
0.3048 31.31 pi /s
 1.1278(kg / m ) 
ρair 𝑸𝑸𝟏𝟏 = 𝟑𝟑𝟑𝟑. 𝟑𝟑𝟑𝟑𝒑𝒑 𝒊𝒊𝟑𝟑 ⁄𝒔𝒔
 Ws (m 2 / s 2 ) 
H s ( pi ) =  2  ( 0.3048)
 9.8( m / s )
plus tard
=
Ws 1.471×105 m 2 / s 2

 1.471 × 105 (m 2 / s 2 ) 
H s (=
pi )  2  ( =
0.3048 ) 4.92 × 10 4
pi
 9.8( m / s ) 
𝑯𝑯𝒔𝒔 = 𝟒𝟒. 𝟗𝟗𝟗𝟗 × 𝟏𝟏𝟏𝟏𝟒𝟒 𝒑𝒑𝒑𝒑
Calculez le débit volumique Q1 , la ‘hauteur’ idéale H=W s/g, la vitesse et
le diamètre spécifique nqet Ds (dimensionnels) et le rendement η. Vérifiez
les calculs à l’aide du tableau suivant. Supposez Q2= Q1
Q1 = 31.31( pi 3 / s )
NQ 60000 × (31.31)
1/ 2 0.5
=
Ns = = 101.57 =
H 4.92 × 104 pi
H 3/ 4
(4.924 ×10 )
4 0.75

D2 H 1/4 5.92 /12 × (4.92 ×104 )0.25


Ds = 1/2 0.5
= 1.31
Q (31.31)

Calculez le débit volumique Q1 , la ‘hauteur’ idéale H=W s/g, la vitesse et


le diamètre spécifique nqet Ds (dimensionnels) et le rendement η. Vérifiez
les calculs à l’aide du tableau suivant. Supposez Q2= Q1
Rendement η η = 0.7016
0.4 0.5 0.6 0.7 0.8
plus tard
50 2.42 2.65 2.91
60 1.94 2.15 2.26

Vitesse spécifique dimensionnelle


65 1.77 1.92 2.02 2.14
70 1.66 1.82 1.89 1.96
80 1.44 1.55 1.63 1.68
85 1.36 1.48 1.53 1.57 1.70
90 1.30 1.39 1.43 1.46 1.59
100 1.16 1.25 1.29 1.32
1.32 1.41
N s = 101.57
110 1.07 1.14 Ds 1.17 1.21 1.29
120 1.00 1.06 1.10 1.15 1.22
130 0.91 1.00 1.03 1.08 1.18
=
Ds 1.32 ≈ 1.31
140 0.87 0.96 1.00 1.06
150 0.83 0.94 1.00 1.07
160 0.80 0.91 1.00 1.04
170 0.80 0.91 1.00 1.11
180 0.80 0.91 1.00
Ns 190 0.79. 0.91 1.01
200 0.79 0.91
H gH
Ψ=
n2 D 2

n1 n2 n3

Q Q
Φ=
nD 3
Si on compare une pompe avec elle même pour différentes
vitesses de rotation N0 et N, on a
2 2 3
H  N   D  Q   D  N 
=      =   
H0  0  0
N D  0  0  0
Q D N
1 1

H
2
 N  Q 
2
 H0  2
= =
    H =  2 Q
H0  0  0
N Q  Q0 
𝐻𝐻 = 𝑘𝑘0 𝑄𝑄 2
2 𝐻𝐻0
𝐻𝐻 = 𝑘𝑘0 𝑄𝑄 𝑘𝑘0 =
𝑄𝑄02

Pour une pompe qui tourne à une vitesse N0 et qui opère au point
(Q0,H0), il existe un ensemble de points ayant le même rendement
sur la parabole H=k0Q2. Ces points correspondent à vitesses de
rotation différentes et reçoivent le nom de points homologues
H
H = k 0Q 2

Point d’opération
H0 ( H 0 , Q0 ) Points homologues

H1 Les points sur la courbe H=k0Q2


ont le même rendement
H2 N0

Q2 Q1 Q0 Q
Si on compare différentes pompes similaires qui tournent à la
même vitesse de rotation N= N0
2 2 3
H  N   D  Q   D  N 
=      =   
H0  N 0   D0   0  0  0
Q D N
1 1

H
2
 D Q 
2/3
 H0  2/3
= =
    H =  2/3 Q
H0  0  0
D Q  Q0 
𝐻𝐻 = 𝐾𝐾0 𝑄𝑄2⁄3
2⁄3 𝐻𝐻0
𝐻𝐻 = 𝐾𝐾0 𝑄𝑄 𝐾𝐾0 = 2⁄3
𝑄𝑄0
Pour une pompe opérant sur un point (Q0,H0) à une vitesse de
rotation N0 , il existe un ensemble de points sur la courbe H=K0Q2/3
qui ont le même rendement. Ces points homologues correspondent
à des pompes similaires avec des diamètres différents, mais qui
tournent à la même vitesse de rotation N0
H

( H 0 , Q0 ) H = K 0Q 2/3
Point d’opération
H0 Points homologues

H1 Les points sur la courbe H=K0Q2/3


ont le même rendement
N0

Q1 Q
H = k 0Q 2
Lorsque le diamètre H Même D,
différentes vitesses
ainsi que la vitesse N
de rotation changent,
on passe du point
d’opération A au B
nouveau point C
suivant les parcours: H = K 0Q 2 / 3
H0 Même N, différents
AB, BC A diamètres D
C

N0

Q0 Q
Les points sur la courbe H=K0Q2/3 ont le même rendement
Les points sur la courbe H=k0Q2 ont le même rendement
La hauteur de charge H, le rendement η et la caractéristique du
système hydraulique Hs associé à une pompe sont donnés par les
équations:
H =20 + 0.8333Q − 0.1667Q 2
η 29.643Q − 3.2143Q 2
Hs 10 + 2.1164Q 2

La vitesse de rotation de la pompe est de N=1800 rpm, la charge


H est donnée en mètres, le débit Q est exprimé en lt/s et le
rendement η en %.
Si la vitesse du rotor du même système est augmentée à 3600
rpm, quel sera le débit et quelle sera la puissance au point
d’opération?

Quelle est la vitesse de rotation (en rpm) nécessaire pour


augmenter le débit 1.7 fois la valeur obtenue à 1800 rpm?
H Si la vitesse du rotor du même
 Q   Q 
= Φ =   3 
système est augmentée à 3600 rpm,
 D N  a  D N b
3 quel sera le débit et quelle sera la
puissance au point d’opération?

 Nb   3600 
=Qb =  Qa  =  2.283 4.566( lt / s ) = + 0.8333Q −
 a
N  1800  H 20 0.1667Q 2
=η 29.643Q − 3.2143Q 2
Hs = 10 + 2.1164Q 2
H = 21.034 m,
Q = 2.283 l / s N a = 1800 rpm
N a = 1800 rpm

H=
20 + 0.8333Q − 0.1667Q 2
H=
s 10 + 2.116Q 2

Q
H Si la vitesse du rotor du même
 gH   gH 
=Ψ = 2 2   2 2
système est augmentée à 3600 rpm,
 N D  a  N D b quel sera le débit et quelle sera la
puissance au point d’opération?
2
 Nb 
2
 3600 
=H b =  Ha  = 21.034 84.136 m H = 20 + 0.8333Q − 0.1667Q 2
 a
N  180 0  =η 29.643Q − 3.2143Q 2
Hs = 10 + 2.1164Q 2
H = 21.034 m,
Q = 2.283 l / s N a = 1800 rpm
N a = 1800 rpm

H=
20 + 0.8333Q − 0.1667Q 2
H=
s 10 + 2.116Q 2

Q
H Si la vitesse du rotor du même
= =
H b 84.136 =
m , Qb 4.566( l / s ), N b 3600 rpm système est augmentée à 3600 rpm,
quel sera le débit et quelle sera la
W = ρ gHQ / η = 1000 × 9.8 × 84.136 × 4.566 × 10−3 / η puissance au point d’opération?

𝜼𝜼 = 𝟓𝟓𝟓𝟓. 𝟏𝟏𝟏𝟏%
H =20 + 0.8333Q − 0.1667Q 2
=η 29.643Q − 3.2143Q 2
=Hs 10 + 2.1164Q 2
H = 21.034 m,
Q = 2.283 l / s N a = 1800 rpm
N a = 1800 rpm

H=
20 + 0.8333Q − 0.1667Q 2

H= 10 + 2.116Q 2 𝜼𝜼 = 𝟐𝟐𝟐𝟐. 𝟔𝟔𝟔𝟔𝟔𝟔𝑸𝑸 − 𝟑𝟑. 𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝑸𝑸𝟐𝟐

Q
Le point d’opération
Qb =
1.7 × 2.283 =
3.882 H=41.88m,Q=3.882 l/s
n’est pas un point similaire
Hb =
10 + 2.116Qb2 =
41.88

H = 21.034 m,
Pour une pompe opérant à
Q = 2.283 l / s
différentes vitesses de
N a = 1800 rpm rotation, il y a un ensemble
de points avec le même
H=
20 + 0.8333Q − 0.1667Q 2 rendement sur la parabole
H=k0Q2.
H=
b 10 + 2.116Qb2

Quelle est la vitesse de rotation (en rpm) nécessaire pour augmenter le débit 1.7 fois la valeur
obtenue à 1800 rpm?
H = k 0Q 2

Qb =
1.7 × 2.283 =  H0 
3.882
k0 =  2 
Hb =
10 + 2.116Qb2 =
41.88  Q0 
 41.88 
H = 21.034 m, = = 2 
2.779
 (3.882) b
Q = 2.283 l / s
N a = 1800 rpm
H = 2.279Q 2
Courbe de points similaires
H=
20 + 0.8333Q − 0.1666Q 2

H=
b 10 + 2.116Qb2

Quelle est la vitesse de rotation (en rpm) nécessaire pour augmenter le débit 1.7 fois la valeur
obtenue à 1800 rpm?
=
H k=
0Q
2
2.279Q 2
Courbe de points similaires
Hb =
10 + 2.116Qb2 =
41.88 m
Qb =
1.7 × 2.283 =
3.882 l / s

H = 21.034 m ,
Q = 2.283 l / s Q Q
  = 
H=
20 + 0.8333Q − 0.1666Q 2 N a = 1800 rpm  N a  N b

H=
b 10 + 2.116Qb2 H a = 21.031 m
Qa = 2.7507 l / s
Nb

Le point a décrit par Ha =21.031m, Qa =2.7507 l/s, pour une vitesse de rotation Na =1800 rpm, est
similaire par rapport au point b avec Hb =41.88m ,Qb =3.882 l/s et Nb=2540.29 rpm

Vous aimerez peut-être aussi