Similitude de Froude en turbomachines
Similitude de Froude en turbomachines
(TERENCE,Phormio )
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:
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:
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
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 )
gH Q ρ ND 2 li
= f , , (design)
2
N D 2 ND
3
µ D
Q ρ ND 2 li
η = g , , (design)
ND
3
µ D
…
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
24
…
𝑄𝑄 𝛥𝛥𝛥𝛥⁄𝜌𝜌
Φ= Ψ= 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
∆(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 𝑁𝑁 𝐷𝐷
~1/(Mach)2
m 𝑚𝑚𝑅𝑅𝑇𝑇
̇ 01 𝑚𝑚̇ 𝑅𝑅𝑇𝑇01 𝑅𝑅𝑇𝑇01 𝑚𝑚̇ 𝑅𝑅𝑇𝑇01
= = 2
⋅ 2
→
ρ 01 ND 3
𝑝𝑝01 𝑁𝑁𝐷𝐷 3 𝑝𝑝01 𝐷𝐷 𝑁𝑁𝑁𝑁 𝑝𝑝01 𝐷𝐷 2
p01=ρ01RT01 ~(1/Mach)
η Ψ Φ Re M γ
P
Dans le contexte des turbomachines on identifie 7 paramètres
adimensionnels
Ψ = f( Φ, Re, M, γ, Géo)
η = g( Φ, Re, M, γ, Géo)
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
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
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
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à
𝑊𝑊̇ 𝑙𝑙𝑖𝑖
𝜂𝜂 = 𝑔𝑔 ,
3⁄2 𝐷𝐷
𝜌𝜌𝐷𝐷 2 𝑔𝑔𝑔𝑔
Les paramètres adimensionnels qui caractérisent l’opération
d’une turbine hydraulique sont:
Coeff. de charge gH Q
α= f 2 2 , 3
n D nD
α ψ Φ
η
Ψ
Coeff. de charge gH Q
η= g 2 2 , 3
n D nD
η ψ Φ
α d’ouverture rendement
η
Ψ
Coeff. de charge
𝑄𝑄 𝑄𝑄11 𝑄𝑄
Φ= = 𝑄𝑄11 = 2 1⁄2
𝑁𝑁 × 𝐷𝐷 3 𝑁𝑁11 × 𝟏𝟏𝟑𝟑 𝐷𝐷 𝐻𝐻
Puissance réduite P11
𝑊𝑊̇
𝑃𝑃11 = 2 3⁄2
𝐷𝐷 𝐻𝐻
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
Ψ Φ
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
H 0 , Q0 W0 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 W0 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
𝐻𝐻
ns = 35.5
ns = 130 ns = 1200
nq = NQ1/2 H −3/4
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
10000
15000
20000
1000
1500
2000
3000
4000
5000
6000
500
600
700
800
900
US
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 𝑾𝑾̇ = 𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐.
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.
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.
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
𝜋𝜋𝑫𝑫𝑁𝑁
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
φ 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
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
6400 mixte
2.33
D2 D1
D1
Ns basse Ns haute
Pompes
Q
Pompes
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
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 …
5
NQ 1/ 2
nq =
H 3/ 4 …
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
𝒄𝒄𝟐𝟐𝐮𝐮
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
g adimensionnelle
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%
10000
15000
20000
1000
1500
2000
3000
4000
5000
6000
500
600
700
800
900
US
c2=
ua U=2 472.39 m / s (90 0
)
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
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
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
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
AB, BC 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
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
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