Cours de Turbomachine A Fluide Compressible

See discussions, stats, and author profiles for this publication at: https://www.researchgate.
net/publication/268401367
Cours de turbomachine à ﬂuide compressible
Article · May 2007
CITATIONS READS
2 5,276
1 author:
Xavier Ottavy
CNRS / Laboraroire de Mécanique des Fluides et d'Acoustique
66 PUBLICATIONS 557 CITATIONS
SEE PROFILE
Some of the authors of this publication are also working on these related projects:
AXIOOM: ANR+NSFC View project
Dynamics of the corner stall separation in a compressor cascade View project
All content following this page was uploaded by Xavier Ottavy on 21 January 2015.
The user has requested enhancement of the downloaded file.
Cours de turbomachine – Xavier Ottavy (CNRS – Laboratoire de Mécanique des Fluides et d’Acoustique UMR 5509 - ECL)
Cours de turbomachine
à fluide compressible
Xavier OTTAVY
CNRS UMR 5509
Laboratoire de Mécanique des Fluides et d’Acoustique
à l’École Centrale de Lyon
Plan du cours
Introduction
Analyse thermodynamique monodimensionnelle
Analyse de l’écoulement dans le plan circonférentiel
Équation de l’équilibre radial simplifié
Cas d’application : dessin d’un étage de compresseur axial

1ière partie - Introduction
Introduction
• Définition
• Fonctions et domaines d’utilisation des turbomachines
• Notion d’étage – échanges d’énergies
• Courbes caractéristiques
Introduction
• Approches 1D, 2D, 2.5D et 3D

Définition
Définition
Une turbomachine est une machine tournante qui réalise un transfert

d’énergie entre son arbre propre, et un fluide en mouvement. Ce transfert
peut s’effectuer dans les deux sens :
Introduction
• une récupération de l’énergie du fluide sur l’arbre de la machine

(fonction réalisée par les machines de type turbine)
• une augmentation de l’énergie du fluide par fourniture d’énergie

mécanique sur l’arbre de la machine (fonction réalisée par les
machines de type compresseur, ventilateur, pompe …)
Fonctions et domaines d’utilisation des turbomachines
Récupération de l’énergie d’un fluide (turbines)

liquide : récupération d’énergie potentielle hydraulique (barrages,…)
gaz : turbines de dentiste, turbocompresseurs, turbopompes, …
turbines associées à d’autres éléments (compresseurs, chambres de
combustion,…) pour la production d’énergie mécanique, ou pour la propulsion en
aéronautique.
Introduction
Compression de gaz (compresseurs)

fonction qui se présente dans des domaines très diversifiés : industrie chimique
(pression de réaction), industrie pétrolière (extraction du pétrole), ou simplement
création d’air comprimé.
compresseurs associés à d’autres éléments (turbines, chambres de combustion,…)
pour la production d’énergie mécanique, ou pour la propulsion en aéronautique.
Transport de fluide
élévation : fournir une énergie pour vaincre le champ gravitationnel (pompes)
transport horizontal : apport périodique d’énergie au fluide pour vaincre les pertes
de charges (boosters)
Ventilation
Production d’énergie mécanique à partir d’une source de chaleur

Production réalisée par des turbines à gaz ou des turbines à vapeur. Ces machines
associent dans un cycle thermodynamique turbines, compresseurs, sources de
chaleur, refroidisseurs,… Puissance variant de quelques kW à plusieurs dizaines de
MW.
Production d’énergie électrique (aérospatiale, avions, chars, réseau nationale,…)
Introduction
Production d’énergie mécanique : entraînement d’hélice de bateau, d’avion

(turbopropulseur), de rotor d’hélicoptère …
Turbines à vapeur essentiellement destinées à la production de forte puissance
d’énergie électrique dans les centrales thermiques.
Propulsion par réaction

Ces machines associent dans un cycle thermodynamique turbines, compresseurs,
chambres de combustions, tuyères…
Turboréacteurs
Turbofans (multiflux)
Introduction PW4156 - Pratt & Whitney (epower-propulsion.com)

Introduction CFM56 – Snecma Moteurs (epower-propulsion.com)

Introduction BR715 – BMW/Rolls-Royce (epower-propulsion.com)

Introduction GP7000 – EA (GE / Pratt & Whitney) (epower-propulsion.com)

Introduction PEGASUS – Rolls-Royce (epower-propulsion.com)

Introduction F404 – General Electric (epower-propulsion.com)

Notion d’étage – Échanges d’énergies
Notion d’étage – Échanges d’énergies
Géométries des turbomachines
Les géométries sont très diverses (de l’éolienne à la Pelton), mais une majorité
des turbomachines peut être répertoriée en 3 catégories :
Introduction
• Les machines axiales : le fluide entre et sort avec une vitesse débitante
approximativement axiale. Machines caractérisées par des débits importants,
mais des taux de pression limités (de l’ordre de 1,4 pour un compresseur
transsonique et de 2 pour un compresseur supersonique).
• Les machines centrifuges : le fluide sort approximativement dans un plan

radial, l’entrée pouvant ne pas être radiale. Machines caractérisées par des
débits limités et des taux de pression important (pouvant atteindre 10 grâce
au travail de la force de Coriolis et à l’augmentation de la pression statique
liée à l’action de la force centrifuge.
• Les machines mixtes

Notion d’étage – échanges d’énergies
Notion d’étage - échanges d’énergie
Un étage de turbomachine se compose d’une partie mobile appelée rotor (ou rouet) et
d’une partie fixe appelée stator (ou selon le cas : redresseur, distributeur,
diffuseur,…)
Introduction
• Le rotor :
Rôle : assurer le transfert d’énergie entre l’arbre de la machine et le fluide en

mouvement.
L’écoulement étant défléchi au passage de la roue, il existe donc une force exercée
par le fluide sur les aubages.
Le point d’application de la force se déplace du fait de la rotation des aubages, il y
a donc travail => échange d’énergie
• Le rotor (suite):
Énergie de pression : une turbomachine échange nécessairement de l’énergie de

pression avec le fluide (même si cela ne doit pas être sa fonction principale).
Cas compresseur : augmentation de la pression pour compenser les pertes de charge du
circuit.
Cas turbine : une partie de l’énergie récupérée l’est toujours sous forme de pression.
Introduction
Énergie cinétique : une turbomachine échange nécessairement de l’énergie

cinétique avec le fluide du fait de la giration de l’écoulement au passage de la roue
mobile.
Énergie calorifique : il n’y a pas d’énergie calorifique directement échangée entre

le fluide et la roue.
Cependant le fluide peut recevoir de la chaleur naissant de la dégradation d’une partie de
l’énergie cinétique <= travail des forces de frottement liées à la nature visqueuse du
fluide.
Phénomène de dissipation principalement localisé près des parois = transformation de la
forme d’énergie et non transfert de l’énergie (« pertes » => rendement).
Faible surface des parois en rapport avec les grands débits rendent les échanges de
chaleur avec l’extérieur négligeable => parois considérées comme adiabatiques
• Le stator :
Rôle : modifier la forme d’énergie (énergie cinétique en pression, ou

inversement).
Il existe comme pour la roue mobile une force exercée par le fluide sur les aubages,
liée à la déflection de l’écoulement.
Par contre l’aubage étant fixe, il n’y a pas de déplacement du point d’application de
la force. Donc pas de travail => pas d’échange d’énergie
Introduction
Redresseur de compresseur axial :

Situé en aval de la roue mobile
Rôle : redresser l’écoulement vers la direction axial, transformant ainsi l’énergie
cinétique de la composante giratoire de vitesse en pression statique.
« Orienter » le fluide dans une direction compatible avec le prochain étage.
Distributeur de turbine axiale :

Situé en amont de la roue mobile
Rôle : provoquer une giration de l’écoulement, transformant ainsi une partie de l’énergie
de pression statique disponible sous forme d’énergie cinétique. Cette énergie est ensuite
récupérée au niveau de la roue mobile.
Diffuseur de compresseur centrifuge :

Récupération de pression statique avec l’augmentation de la section de passage (rayon).
Courbes caractéristiques
• Compresseurs et turbines sont en général calculés pour un point de

fonctionnement (débit massique m et taux de pression Π) où le rendement est
maximal : c’est le point de fonctionnement nominal.
Introduction
• Il est cependant intéressant de connaître le comportement de la machine à

d’autres débits, d’où la notion de plage de fonctionnement. C’est la
fourchette de débit où la machine conserve un taux de pression acceptable
avec un rendement acceptable. Ce fonctionnement hors adaptation est illustré
sur les courbes caractéristiques.
• Actuellement, les recherches sont largement orientées sur l’extension des

plages de fonctionnement :
Interaction rotor/stator – rôle des effets potentiel dans l’amorce du décollement
tournant
Traitement du carter pour profiter de son interaction avec les écoulement de jeux.
But : repousser la zone de pompage.
Courbes caractéristiques des compresseurs

Débit Taux de pression Rendement isentropique Vitesse de rotation
γ −1
m& T pt 2 γ
Π −1 N
m& = t1
Π= ηis =
Tt 2
red
p pt1 −1 T
t1 t1
Tt1
Introduction
• Exemple de déclenchement d’un décollement tournant

Thèse : Nicolas Gourdin – ONERA/ECL - calcul elsA 2,5D - compresseur subsonique CME2 de Snecma Moteurs
Introduction
Courbes caractéristiques des turbines
Débit réduit
Rendement
Introduction
• Asymptote commune, indépendante de la vitesse de rotation => blocage sonique dans

distributeur (partie fixe)
•Plage de rendement très étalées => caractère accéléré de l’écoulement au passage des aubes
(contrairement aux compresseurs où l’écoulement décélère)
Approche 1D, 2D, 2.5D et 3D
Approche 1D, 2D, 2.5D et 3D

Écoulement réel dans une turbomachine complexe :
tridimensionnel
visqueux
instationnaire
Introduction
Définition de surfaces «méridienne» et «aubes-à-aubes»

Surface méridienne
Surface aubes-à-aubes
2ième partie - Analyse thermodynamique monodimensionnelle
Introduction
Analyse thermodynamique 1D

• Équations de conservation de base
• Bilan des différentes contributions
• Équation de l’entropie (équations de Gibbs)
• Travail mécanique et travail utile
• Conditions d’arrêt
• Intérêt du diagramme entropique et enthalpique
• Rendements isentropiques
• Rendements polytropiques

Équations de conservation de base

Équations de conservation de l’énergie totale
de d (V 2 ) dq dwe
+ = + forme locale générale

dt dt dt dt
e est l'énergie interne par unité de masse (qui n'inclue pas l'énergie cinétique du
mouvement d'ensemble macroscopique des molécules, mais uniquement l'énergie
cinétique liée à l'agitation de nature aléatoire de celles-ci)
q représente l'énergie calorique massique échangée avec l'extérieur
we est le travail des forces extérieures par unité de masse. Après utilisation de
l’équation de continuité, on a :
r
d (1 ρ )
( )
dwe r r 1 r V r
= f ⋅ V + div τ ⋅ V − grad p − p
r
dt ρ ρ dt
f représente les forces extérieures par unité de masse, appliquées sous forme volumique
τ est tenseur des contraintes visqueuses
Le travail des forces de pression comporte 2 termes :
r V 2
V r
1 − grad p qui est un terme de transport
ρ 1
d (1 ρ )
2 −p qui est un terme de compressibilité
dt
Expression adaptée à un système ouvert, en introduisant l’enthalpie (h):

p
h=e+
ρ
L’équation d’énergie devient :
dh d (V 2 / 2) 1 ∂ p 1
( )
r r r dq
+ = + div τ ⋅ V + f ⋅ V +
dt dt ρ ∂t ρ dt
En exprimant cette équation sous la forme :
dh d (V 2 / 2) dq dwT
+ = +
dt dt dt dt
On a donc par identification :
dwT r r 1 1∂p
( )
r
= f ⋅ V + div τ ⋅ V +
dt ρ ρ ∂t
où dwT/dt représente la puissance dite utile par unité de masse

Équations de conservation de l’énergie cinétique
Cette équation peut se substituer à l'équation de la dynamique ou à celle de la

quantité de mouvement. Elle s'exprime par :
d(V 2 2) dwe dwi puissance des forces extérieures

= + puissance des forces intérieures
dt dt dt
dw i σ: D τ :D p r τ:D d (1 ρ )
=− =− + divV = − +p
dt ρ ρ ρ ρ dt
rr t rr
où D est le tenseur des déformations : D =
1
2
grad V+ gra d V ( )
Soit, finalement :
r
( )
d (V 2 2) r r 1 r V r 1
= f{⋅ V + div τ ⋅ V − grad p − τ :D
dt ρ ρ 243 ρ
travail des 14243 14 123
forces travail des travail des dissipation d'
de volume forces forces de pression énergie cinétique
visqueuses (transport) par viscosité
( )
r
div τ V ρ = terme de production dû au travail des forces extérieures de viscosité
τ : D ρ = terme de dissipation interne (irréversibilité mécanique) dû au travail
intérieur de la viscosité
Équations de conservation de l’énergie interne
Cette équation s'obtient directement en retranchant membre à membre les deux

équations précédentes :
de dq dwi
= −
Soit : dt dt dt
de dq 1 d (1 ρ )
= + τ :D −p
dt dt
{ ρ
123 dt 3
1424
apport extérieur dissipation travail de
de chaleur mécanique compression
(dit volumique)
Remarque : dans l'équation de l'énergie totale les termes de pression interviennent

dans le travail surfacique extérieur : la pression est donc une pression extérieure.
Dans l'équation de l'énergie interne, ces termes proviennent du travail des forces
internes : la pression est donc une pression intérieure.
d(p ρ )
Expression adaptée à un système ouvert (en ajoutant aux 2 membres) :
dt
dh 1 dp dq 1
= + + τ :D
dt ρ dt dt ρ
Bilan des différentes contributions
Bilan des différentes contributions

Énergie totale
r
d (1 ρ )
( )
2 r r r r
de d (V 2) dq 1 V
+ = + { f ⋅ V + div τ ⋅ V − grad p − p
dt dt dt
{ ρ
travail des 14243
ρ
1424 3 1424 dt 3
apport extérieur forces travail des travail des travail de
de chaleur de volume forces forces de pression compression
visqueuses (transport) (dit volumique)
Énergie cinétique
r
( )
2
d (V 2) r r 1 r V r 1
= { f ⋅ V + div τ ⋅ V − grad p − τ :D
dt ρ ρ
travail des 14243 14243
ρ
123
forces travail des travail des dissipation d'
de volume forces forces de pression énergie cinétique
visqueuses (transport) par viscosité
Énergie interne >0 <0
de dq 1 d (1 ρ )
= + τ :D −p V
dt dt
{ ρ
123 dt 3
1424
apport extérieur dissipation travail de
de chaleur mécanique compression
(dit volumique)
Équation de l’entropie
Équation de l’entropie (équations de Gibbs)

En divisant l'équation de l'énergie interne par T en tenant compte de :
r r r1
div q  q  r
= div − q ⋅ grad
T  T T
d (1 ρ ) 
r r
1  de 1 dq 1τ :D r 1 −q q r1 1τ :D
 + p  = + = + div + ⋅ grad +
T  dt dt  T dt
{ T ρ
123 T 4
14 ρ244 T
3 ρ 243
14 T T ρ
123
apport de irréversibilité apport de chaleur irréversibilité irréversibilité
chaleur mécanique réversible thermique mécanique
Le deuxième membre de l'équation précédente s'identifie comme la variation

d'entropie du système par unité de masse ds
dt
D'où la première équation de Gibbs :
de d (1 ρ ) ds
+p =T
dt dt dt
La deuxième équation de Gibbs s'obtient en introduisant l'enthalpie :
dh 1 dp ds
− =T
dt ρ dt dt
Travail mécanique et travail utile

Interprétation physique du travail utile
(∂Df )
f
V1 V2
(∂D1 ) (∂D'1 ) Dm (∂D2 ) dm (∂D'2 )
dm
m = ρ ϑ(t)
m = ρ ϑ(t+dt)
masse de gaz m=ρϑ(t), comprise entre les sections d'entrée ∂D1 et de sortie ∂D2 et occupant
le domaine matériel Dm.
à l'instant t+dt cette masse de fluide m=ρϑ(t+dt) se sera déplacée vers l'aval.
Puissance utile échangée avec cette masse m :
dWT dw
P= = ∫ ρ T dϑ
dt Dm dt
Expression du travail utile, vue précédemment :

dwT r r 1 1∂p
( )
r
= f ⋅ V + div τ ⋅ V + (cf. énergie totale)
ρ ρ ∂t
dt
Donc :
∂p
( )
dwT r r r r
P=∫ ρ dϑ = ∫ ρ f ⋅ V dϑ + ∫ τ V ⋅ n dS + ∫ dϑ
Dm dt Dm ∂ Dm Dm ∂ t
Avec un écoulement permanent, une condition d'adhérence imposant une

vitesse nulle sur les parois fixes (∂Df) et des tensions visqueuses
négligeables sur ∂D1 et ∂D2 (au sein du fluide), on a :
dW T r r
dt
≈ ∫D m
ρ f ⋅ V dϑ
r
ou f représente les forces de volume autres que les forces de pesanteur,
c.à.d. dans notre cas les forces exercées par les aubages de la machine
sur le fluide.
Le travail utile représente donc le travail échangé avec la machine.

Autres expressions du travail mécanique et du travail utile
• Travail mécanique :
En ajoutant les variations d'énergie cinétique dans l'équation de l'énergie interne :
de d (V 2 ) dq 1 d (1 ρ ) d (V 2 )
2 2
+ = + τ :D− p +
dt dt dt ρ dt dt
et, en identifiant avec l'équation de l'énergie totale :
d (1 ρ ) d (V 2 )
2
dwe dwd
= −p +
dt dt dt dt
• Travail utile :
En effectuant la même opération avec l'équation d’enthalpie :
dh d (V 2 ) dq dwd 1 dp d (V 2 )
2 2
+ = + + +
dt dt dt dt ρ dt dt
et, en identifiant avec l'équation de l'énergie totale (enthalpie) :
dwT dwd 1 dp d (V 2 2)
= + +
dt dt ρ dt dt
Différence entre puissance mécanique et puissance utile
dw T dw e d( p ρ )
= +
{dt {dt 1
42 dt4
3
puissance puissance puissance de
utile mécanique transvasement
soit, en intégrant sur le domaine Dm défini précédemment :

dW T dWe d( p ρ )
− = ∫D ρ dϑ
dt dt m dt
d p ∂p r r
ρ dϑ = ∫ dϑ + ∫
dt ∫D m ρ
= p V ⋅ n dS
Dm ∂ t ∂ Dm
Si on fait l'hypothèse d'un écoulement permanent, et si nous tenons compte

du fait que la condition d'adhérence impose une vitesse nulle sur les parois
fixes (∂Df), il vient :
dW T dWe r r
dt
−
dt
= ∫∂D ∪∂ D
1 2
p V ⋅ n dS
Le travail utile diffère donc du travail mécanique par le travail des forces de pression
sur les surfaces libres d'entrée et de sortie que l'on appelle travail de transvasement.
Analyse thermodynamique 1D • Illustration : Compresseur et son entrée d’air
(∂D0) (∂D1) (∂D2)

V
D1
D0
domaine D1: travail utile échangé (par l'intermédiaire des pales mobiles) et
travail des forces de pression sur les surfaces de contrôle (∂D1) et (∂D2) (car
les pressions p1 et p2 sont différentes).
∆wT ≠ 0 et ∆we ≠ 0
domaine D0 : aucun travail échangé avec les parois matérielles
Mais : puissance des forces de pression nulle sur (∂D0) (si on la situe
suffisamment loin pour considérer que la vitesse du fluide y est nulle), elle ne
l'est pas sur (∂D1). ∆wT = 0 et ∆we ≠ 0
Conditions d’arrêt
Notion de variable d’arrêt absolue
Écoulement uniforme (p,T,ρ) et permanent de fluide

Pas d’échange de travail ni de chaleur avec l’extérieur
Plaçons y une sonde de température et de pression
Équation de conservation de l’énergie M1 M2

x
entre les points 1 et 2, sur ces 2 lignes
p, T, ρ, V
de courant
x
M'1 M'2
∆h + ∆(V 2) = ∆w T + ∆q = 0
2
L'enthalpie au point 2 englobe donc à la fois l'enthalpie statique et l'énergie

cinétique en 1, d'où son nom d'enthalpie totale. Mais, c'est aussi l'enthalpie
obtenue après avoir arrêté le fluide, d'où son nom d'enthalpie d'arrêt.
h0 = h + V 2 2
Définition des variables d’arrêt
Les conditions d'arrêt sont les conditions que l'on obtiendrait par une
transformation fictive ramenant isentropiquement (réversiblement et sans échange
de chaleur) et sans échange de travail utile, le fluide à vitesse nulle.
• Température d’arrêt
Dans le cas d’un gaz parfait : h = C p (T ) ⋅ T

h0 V2
=T +
C p (T ) 2 C p (T )
En considérant que : h0 = C p (T0 ) ⋅ T0 ≈ C p (T ) ⋅ T0 ≈ C p ⋅ T0
V2
ainsi : T0 = T +
2Cp
T est la température statique qui représente l'énergie cinétique moyenne d'agitation
moléculaire (de nature aléatoire).
V2/2Cp est la température dynamique qui représente l'énergie cinétique du
mouvement d'ensemble des molécules.
• Pression et masse volumique d’arrêt
En utilisant la définition précédente, qui exprime que l'on passe de l'état dynamique
(p, r, T, V) à l'état d'arrêt (p0, r0, T0, 0) par une transformation isentropique:
γ γ
p0  ρ0   T  γ −1
=  = 0
p ρ T 
γ 1
 V 2  γ −1  V2  γ −1
p0 = p  1 +  ρ0 = ρ  1 + 
 2 C T
p   2C p T 
• Autre formulation
V V
En introduisant la notion de nombre de Mach : M = =
a γ rT
V2 γ r V 2 γ −1 2
Il vient : = = M
2 Cp T 2 Cp a 2 2
1 γ
 γ −1 2   γ −1 2  γ −1  γ −1 2  γ −1
T0 = T  1 + M  ρ0 = ρ 1 + M  p0 = p  1 + M 
 2   2   2 
Approximation pour un écoulement incompressible
Les effets de compressibilité sont négligeables si le fluide se déplace à basse

vitesse. On considère généralement que la limite acceptable de l'approximation est
M < 0,3. Dans ce cas, pour l’air :

γ −1 0.4 2
M2 ≈ 0.3 ≈ 0.02 << 1
2 2
Soit en effectuant un développement limité au premier ordre :
 γ γ −1 2  γ p 2
p0 ≈ p  1 + ⋅ M  = p+ ⋅ V
 γ −1 2  2 γ rT
1
p0 ≈ p + ρ V
2
2
Cas de l’entropie
L'entropie statique s'exprime, en fonction des variables p et ρ et de la définition des

variables d'arrêt :
( )
s = C v Log p ρ γ = C v Log p0 ρ0 γ = s 0 ( )
s ≡ s0
Ce résultat pouvait être prévu sans aucun calcul, car l'entropie étant associée à une
quantité de chaleur (échangée ou produite par irréversibilités), elle ne dépend pas
du repère considéré.
Variables statiques et d’arrêt absolues ou relatives
p0 et T0 sont les pression et température « ressenties » par des sondes fixes

Pour mesurer les variables statiques, il ne faut pas modifier la dynamique du fluide.
Nécessité de sondes solidaires d'un référentiel lié au fluide (difficile).
En pratique: sondes de paroi parallèle à l'écoulement de fluide.
Si on se rappelle que le phénomène de pression résulte de l'intégrale des forces

d'impact des molécules sur la paroi, alors :
capteur sur paroi parallèle au mouvement du fluide, la pression mesurée n'inclut que les
chocs liés au mouvement aléatoire des molécules : c'est la pression statique
capteur sur paroi perpendiculaire au mouvement du fluide, la pression mesurée inclut à la
fois les chocs liés au mouvement aléatoire des molécules et ceux liés à leur mouvement
d'ensemble de vitesse : c'est la pression totale.
Lorsqu'une sonde est liée à un référentiel qui n'est ni fixe, ni solidaire des particules
de fluide, elle mesure une quantité que l'on appelle variable d'arrêt relative.
Exemple : pales d'un compresseur animées d'un mouvement de rotation uniforme :
Le fluide possède, dans le référentiel lié aux pales, une certaine vitesse W (dite vitesse
relative, différente de la vitesse absolue V). La température "ressentie" par la pale au point
d'arrêt est la température d'arrêt relative :
T 0 R = T + W 2 2C p
Équation thermique d’état
• Pour un gaz parfait, la définition des variables d'arrêt entraîne :

p = ρ rT
1
 γ −1 2  1− γ  γ −1 2  −1
= ρ0  1 + M  ⋅ r ⋅ T0 1 + M 
 2   2 
γ
γ −1 2  1− γ
= ρ0 r T 0  1 + M 
 2 
p0 = ρ0 r T 0
• Ce résultat est évident si on considère que l'état d'arrêt (même s'il est fictif) est
un état au sens thermodynamique : les variables y caractérisant le fluide sont
donc régies par l'équation thermique d'état.
Expression du travail utile pour une transformation adiabatique

• Expression du travail utile et de l’entropie établies précédemment :
1
dwT = dp + d (V 2 2) + dwd
ρ
Pour un écoulement adiabatique, les

dq δq δq irréversibilités de nature mécanique
ds = ext + ext + méca dwd sont les seules causes
T rév T irrév T irrév
d'augmentation d'entropie
144 42444 3
A
• Transformation adiabatique => A=0, donc

1
dwT = dp + d (V 2 2) + Tds
p ρ p + dp Conditions
ρ ρ + dρ statiques
1 dwT ≠ 0 2
T T + dT
r ds ≠ 0 r r
V V + dV
dwT ≡ 0 ds ≡ 0 dwT ≡ 0 ds ≡ 0
p0 p0 + dp0
dwT ≠ 0
ρ0 ρ0 + d ρ0
1’ ds ≡ ds0 2’
T0 T0 + dT0 Conditions
r 1 r d’arrêt
V =0 dwT = dp0 + 0 + T0 ds0 V =0
ρ0
• Rappel : Équation d’énergie exprimée avec l’enthalpie :

dh + d (V 2 / 2) = dq + dwT
dh0 = dq + dwT
• Cas adiabatique : dq=0, donc :
dh0 = dwT
1
dh0 = C p dT0 = dwT = dp0 + T0 ds0 Pertes par dissipation
ρ0
• Cas transfo isentropique (adiabatique + réversible) entre 1 et 2, avec un gaz parfait :

γ −1
 
 T02    γ
− 1 = C pT01   02  − 1
p
∆WT = C p (T02 − T01 ) = C pT01 
 p 
 T01    01  
 γ γ −1
  
∆WT = C pT01  Π − 1 Isentrop.
 
 
Compression : Π > 1 => T02>T01 => ∆WT12 > 0
Détente : Π < 1 => T02<T01 => ∆WT12 < 0
Importance de la température d’entrée
Diagrammes entropique et enthalpique
Diagrammes entropique (T,s) et enthalpique (h,s)

Diagrammes couramment utilisés pour représenter les transformations.
Si gaz parfaits évoluants dans une plage de températures limitée,
température et enthalpie sont identiques à une constante près :

h = h(T) = C p (T) ⋅ T ≈ C p T
Intérêt d’un tel diagramme

• les ordonnées représentent l'énergie du système. On peut donc directement y
visualiser :
pour une transformation adiabatique : ∆h = ∆w T + ∆q = ∆w T
Les échanges de travail utile
pour une transformation à p constante : ∆h = ∆w T + ∆q = ∆q
Les échanges de chaleur
pour un système isolé : ∆h = 0 :
Les transferts internes (énergie cinétique énergie de pression)
• dans le cas d'une transformation adiabatique, les abscisses représentent le degré

d'irréversibilité de la transformation.
compressions ou détentes usuelles : visualisation directe des transferts d'énergie,
quantitativement en ordonnées et qualitativement en abscisses.
Diagrammes entropique et enthalpique
Quelques iso-valeurs intéressantes

• Isentropiques :
les compressions et détentes usuelles réversibles sont représentées par
des droites verticales.

• Adiabatiques irréversibles :
selon le second principe ∆wd > 0, que ce soient des compressions ou des détentes,
l'entropie augmente.
• Isothermes (ou isenthalpiques) :
ce sont des droites horizontales.
• Isobares (cas des apports de chaleur usuels) :
Après intégration de l’équation de Gibbs, on montre que les isobares sont des
exponentielles croissantes se déduisant l'une de l'autre par translation horizontale.
Remarques :
l'écart vertical entre deux isobares augmente avec la température : ce résultat est
déterminant pour comprendre le fonctionnement d'une turbine à gaz.
dans le cas d'un gaz non parfait, il est indispensable d'utiliser le diagramme
enthalpique, où toutes ces courbes sont distordues.
ex : vapeur d'eau utilisée dans les turbines à vapeur (le diagramme enthalpique
correspondant est le diagramme de Mollier).
Rendements isentropiques
Ces rendements s'appellent isentropiques parce qu'ils comparent la
transformation réelle à une transformation isentropique fictive.
T T
T2 p2 p1
T2 is T1
p1 p2
T2
T1 T2 is
s s
compression détente
∆wT is h2is − h1 ∆wT h2 − h1
ηc = = η = =
∆wT h2 − h1 ∆wT is h2is − h1
T
Pour un gaz parfait à Cp constant, ces rendements peuvent aussi s'écrire :
T2is − T1 T2 − T1
ηc = η =
T2 − T1 T2is − T1
T
Rendements isentropiques total-à-total
Les rendements isentropiques total-à-total d'une roue mobile sont définis de

façon analogue à la formulation initiale, mais en utilisant les variables d'arrêt
pour conserver leur sens physique de rendements énergétiques globaux

(entrée-sortie).
h0 2 is − h0 1
ηc = pour une compression
h0 2 − h0 1
h0 2 − h0 1
ηT = pour une détente
h0 2 is − h0 1
Dans le cas où Cp constant => rendements définis avec les températures d'arrêt.
Étage de turbomachine :
la partie mobile : transfert d'énergie entre la machine et le fluide. En effet, les pales étant
mobiles, l'ensemble des forces de pression et visqueuses exercées sur le fluide travaillent.
la partie fixe ne réalise qu'une transformation interne de la forme d'énergie du fluide (pas
d'échange d'énergie avec la machine).
Existence de forces (fixes) entre les pales et le fluide => pas de travail.
le premier principe exprime que : ∆h0 = ∆wT + ∆q ≡ 0
il n'est pas possible de caractériser le degré d'irréversibilité de la transformation par ce type
de rendement. On utilise les rendements isentropiques statique à statique
Étage de compresseur
T
p02 is p02 p03

2 3
T0 3=T0 2 V3
2
2 Cp
2
V2
p3
T0 2 is 2 Cp
T02 is − T01
p2 ηc =
T02 − T01
1 p0 1
T01
V12
2 Cp p1
s
Étage de compresseur (suite)
Energie p -p >0
Energie totale
2 1
Roue fournie : mécanique : 2 2
mobile Cp (T0 2-T01) > 0 p0 2 - p01 > 0 V2 - V1 > 0
2 2
Dissipation
visqueuse :
p 02 is- p 0 2 > 0
Energie totale Energie p -p >0

3 2
Roue échangée : mécanique : 2 2
Cp (T0 3 -T0 2) ≡ 0 p0 3 - p02 < 0 V3 - V2 < 0
fixe 2 2
Dissipation
visqueuse :
p 0 2- p 03 > 0
Étage de turbine
T p0 1 p0
1 2
2
T0 2=T0 1
2
V1
2 Cp p1
2
V2
2 Cp
T02 − T01
p2 ηt =
T02is − T01
p03 is 3
T0 3 p0
2 3
V3
T0 3 is 2 Cp p3
s
Étage de turbine (suite)
p -p <0
Roue 1 Energie Energie totale

2
mécanique : échangée :
fixe 2
V2 V1
2
Cp (T0 2 -T0 1) ≡ 0
- > 0 p 0 2 - p 01< 0
2 2
Dissipation
visqueuse :
p 0 1- p 02 > 0
p -p <0 Energie Energie totale

2
Roue 3
mécanique : récupérée :
2 2
mobile V3 V2
- < 0 p 0 3 - p 02 < 0 Cp (T0 2-T0 3) > 0
2 2
Dissipation
visqueuse :
p 0 is- p 03 > 0
3
Rendements polytropiques
Processus polytropiques
• Une définition :
une transformation polytropique est une transformation au cours de laquelle le
rapport entre la chaleur totale échangée, et la variation d'enthalpie est égale à une
constante β.
dq + dwd
= β = constante
dh0
Pour un gaz parfait : dh0 = dwT + dq = dp ρ + d (V ) + dwd + dq

2
2
dh − dp ρ C p dT − dp ρ
(sans variation d’Ec) ⇒β = =
dh C p dT
1 1  1 
avec : p = ρ rT ⇒ dT =  dp + p d   
rρ  ρ 
d’où : Cp  1  1 1
 dp + pd    (β −1) = − dp
r ρ ρ  ρ
Cp 1  Cp 1
⇔ (β − 1) ⋅ p d   +  (β − 1) +1 dp = 0
r ρ   r ρ
Cp γ
soit, en utilisant la relation de Mayer : C p − Cv = r ⇒ =
r γ −1
et en multipliant toute l'équation par ρ ⋅ γ −1 ,

p γ β −1
dp γ (β − 1) d(1 ρ )
il vient : + ⋅ =0
p γ β −1 1 ρ
γ (β − 1)
Soit, en posant : n=
γ β −1
p
Et en intégrant => = constante
ρn
On appelle n exposant polytropique, qui est donc constant tout le long de la

transformation.
• Remarques sur la valeur de l’exposant polytropique :
Transformation isotherme : nc = 1
Transformation adiabatique réversible : nc = γ
Compression : les phénomènes dissipatifs créent de la chaleur et éloignent encore

plus la transformation adiabatique de l'isotherme (on a un sur-échauffement par
rapport à l'adiabatique réversible). Donc : nc ≥ γ
Détente: les phénomènes dissipatifs diminuent le refroidissement naturel de

l'adiabatique réversible (ils rapprochent la transformation réelle d'une
transformation isotherme). Donc : 1 ≤ nT ≤ γ
• Assimilation du processus adiabatique irréversible réel à un processus

polytropique:
dq + dwd
= β = constante = coéfficient de pertes
dh0
Avec des effets dissipatifs (pertes) répartis de manière "relativement" uniforme tout
le long de la transformation, on peut assimiler le processus réel à un processus
polytropique (cf. développement des couches limites et des pertes qu’elles génèrent,
de façon relativement constante tout au long de la surface si le profil turbulent est
établi).
Considérons une compression adiabatique irréversible de p1 à p2, et soit un

élément infinitésimal de cette transformation compris entre les pressions pi
et pi+1 = pi+dp.
T p2 Le rendement isentropique de cette
transformation élémentaire est, par
définition :
(i+1) p i+1
(i+1)' h(i +1)' − h(i ) dhis
=
h(i+1) − h(i ) dh
pi
(i) avec : dhis = dwT + dq = dp ρ
is
p1
s dp ρ
η Pc =
dh
Ce rendement local, constant tout au long de la transformation, est appelé

rendement polytropique.
Pour une détente, on définit de façon analogue : dh
η PT =
dp ρ
Relations entre rendements et exposants polytropiques

dq + dw d dh − dp ρ
Rappel : β = =
dh dh
d'où, en utilisant l'expression du rendement polytropique :
pour une compression : ηPc = 1 − β

1
pour une détente : η PT =
1−β
γ (β − 1)
Comme, par définition : n =
γ β −1
n −γ
β (n γ − γ ) = n − γ ⇒ β=
γ (n − 1)
d'où, pour un compresseur :
n ( γ − 1)
η Pc = 1 − β =
γ (n − 1)
Ou bien : n − 1 = γ − 1 ⋅ 1
n γ η Pc
Pour une détente, on montre de manière analogue, que :
n −1 γ −1
= ⋅ η PT
n γ
Relations entre rendements isentropiques et polytropiques
Travail utile échanger au cours d’une détente : ∆wT = C p (T2 − T1 )

γ −1
 
En utilisant le rendement isentropique :

  γ
∆wT = C pηT (T2is − T1 ) = C pηT T1   2  − 1
p
 p 
 1  
 
En utilisant le rendement polytropique : n −1
 
 T2   
∆wT = C pT1  − 1 = C pT1   2  − 1
p n
 p 
 T1   1  
 
n−1 γ −1
Avec : = ⋅ η PT
n γ
on montre que : ηT > η PT ηT = η PT (1
14 +2f4
3)
∞ , avec f∞ > 0
facteur de
récupération
Dans le cas d'une compression, on montre, d'une manière analogue, que :

η Pc
ηc < η Pc ηc =
(1+ f ∞ )
1424 3
facteur de
pertes
Interprétation
Comparaison entre compression isotherme et adiabatique
p
2 1
∆wT = ∫1 ρ dp
p2
adiabatique
rT  dp  p
p = ρrT = ⇒   =−
isotherme 1ρ  d(1 ρ ) Tct 1ρ
K  dp  p
p= ⇒   = −γ
(1 ρ )γ  d(1 ρ ) Isent . 1ρ
p
1
1/ρ 1/ρ
(∆wT )adiab. > (∆w T )isoth.
1
phénomènes dissipatifs (créateurs de chaleur) éloignent l'adiabatique de l'isotherme

dans le cas d'une compression et à la rapprocher dans le cas d'une détente. Outre
leur effet néfaste qui est la dégradation d'énergie mécanique, ils ont donc un
"effet thermodynamique" encore négatif dans le cas d'une compression, mais
favorable dans le cas d'une détente.
Le rendement polytropique, du fait de son caractère local, ne caractérise que l'effet
purement dissipatif : c'est un rendement aérodynamique.
Le rendement isentropique, du fait de son caractère global, englobe à la fois l'effet
dissipatif et l'effet thermodynamique qui en résulte.
Analyse de l’écoulement dans le plan circonférentiel 3ieme partie - Analyse de l’écoulement dans le plan circonférentiel
Introduction

• Notion de triangle de vitesse
• Équation d’Euler
• Écoulement en grille d’aubes de compresseur

Notion de triangle de vitesse
Analyse de l’écoulement dans le plan circonférentiel Notion de triangle de vitesse
• A partir de la relation :
r r r r r r
V =W +ω ∧ r =W +U
• Représentation dans le plan aube-à-aube qui permet de visualiser rapidement le

module et la direction du vecteur vitesse en amont et en aval :
d’une roue fixe (repère absolu – vitesse et angle absolus V et α)
d’une roue mobile (repère relatif – vitesse et angle relatifs W et β)
• Outil pratique pour prévoir le fonctionnement d’une roue et estimer :

la charge aérodynamique sur l’aubage
angle d’incidence au bord d’attaque
déflection imposée à l’écoulement
donc force d’aubage (id travail échangé si la roue est mobile)
fonctionnement hors adaptation (possibilité de décollement)
influence d’une variation de rayon, des couches limites pariétales
niveau d’accélération ou de décélération dans la roue,…
Analyse de l’écoulement dans le plan circonférentiel Cas compresseur

cas simple :
machine axiale (pas de variation de rayon et plan aube-à-aube = cylindre),
pas de variation de vitesse axiale axiale
Analyse de l’écoulement dans le plan circonférentiel Cas turbine

cas simple :
machine axiale (pas de variation de rayon et plan aube-à-aube = cylindre),
pas de variation de vitesse axiale
Analyse de l’écoulement dans le plan circonférentiel Vrillage des pales le long de l’envergure (compresseur)
En amont d’un rotor, une augmentation de R (du moyeu au carter) modifie la valeur
de U (=ω.R), et donc les triangles de vitesse.
βt > β p
Pour que l’incidence sur l’aubage

soit bien adaptée sur toute son
envergure, il faut modifier l’angle
de calage des aubages en fonction
du rayon, d’où leur forme vrillée.
Analyse de l’écoulement dans le plan circonférentiel Influence des zones visqueuse pariétales (compresseur)
Après plusieurs étages, le profil de vitesse de l’écoulement dans le plan méridien

met en évidence la diminution de la vitesse débitante Va près des parois.
Progression dans les étages
Dans les CL, l’angle d’incidence est donc

plus élevé => augmentation du travail (cf
Euler).
Mais : surcharge + phénomènes visqueux +
jeu en bout d’aubage (éventuellement)
=> pertes également plus importantes.
Problème critique : lorsque la roue est déjà

très chargée, cette surcharge peut amener à
des décollements des CL d’aubages
=> chute brutale des performances dans ces
zones (pouvant affecter l’ensemble de la
β1CL > β1 roue)
Analyse de l’écoulement dans le plan circonférentiel Allure des courbes caractéristiques d’un compresseur
• Diminution du débit à vitesse de rotation constante :
diminution de la vitesse axiale Va

donc augmentation du travail fourni
Augmentation du taux de pression (dans un premier temps)
Si diminution trop importante du débit :

décollement des couches limites d’aubages (augmentation des pertes)
donc chute du rendement
déflection ne se faisant plus correctement => diminution du travail
Chute importante du taux de pression
Analyse de l’écoulement dans le plan circonférentiel Allure des courbes caractéristiques d’un compresseur (suite)
• Diminution de la vitesse de rotation à débit constant : (Va=cte)
diminution de la vitesse d’entraînement U

Ainsi, diminution de l’angle d’incidence β1
donc diminution du travail fourni
diminution du taux de pression

Analyse de l’écoulement dans le plan circonférentiel Adaptation du débit à la vitesse de rotation et à la caractéristique du
circuit de charge (compresseur)
Si on place un compresseur donné dans un circuit donné et que la vitesse de rotation

N est fixée, le débit Q s’ajuste de lui-même
(de même que pour une turbine donnée dans un circuit donné et disposant d’un
certain taux de détente, la vitesse de rotation s’ajuste d’elle même).
Pourquoi ?
Point de fonctionnement 1 :
intersection de la courbe
caractéristique du compresseur à
N1=cte et de la courbe de pertes
de charge du circuit aéraulique.
(la pression délivrée par la

machine ne sert qu’à vaincre les
pertes de charges comprises
entre les deux infinis amont et
aval)
Analyse de l’écoulement dans le plan circonférentiel Si augmentation brutale de la vitesse de rotation de N1 à N1’ => augmentation de
l’angle β1 => augmentation du travail => augmentation du taux de pression (la
machine peut donc s’opposer à des pertes de charges dans le circuit plus importantes
que précédemment)
De ce fait, le débit va augmenter

jusqu’à ce que le niveau de pertes
revienne équilibrer le taux de
pression délivré. Ceci s’accompagne
d’une augmentation de Va qui
restitue l’allure du triangle de
vitesse comparable à la précédente
(triangles homothétiques).
Pour passer du point de fonctionnement 1 (Π1,Q1) à un autre point de

fonctionnement 1’ (Π1’,Q1’), il y a adaptation du triangle de vitesse.
Ce type d’adaptation se réalise lors de la montée en vitesse de la machine (qui

amène le débit de 0 à Q).
Équation d’Euler
Analyse de l’écoulement dans le plan circonférentiel Équation d’Euler

• L'équation d'Euler relie la quantité d'énergie échangée entre le fluide et les aubages de la
machine, aux caractéristiques aérodynamiques de l'écoulement en amont et en aval de la
roue.
• Cette équation est établie à partir de la projection sur l'axe de la machine de l'équation
intégrale du moment de quantité de mouvement, qui permet d'introduire et d'expliciter le
couple exercé sur l'arbre par le fluide, ou inversement.
Détermination de l’équation d’Euler

• Forme intégrale de l'équation du moment de quantité de mouvement :
δ r r r r r r r r r
ρ ϑ ρ τ
δ t ∫D ∫∂ D ∫∂ D
( r ∧ V ) d = ( r ∧ V )(U − V ) ⋅ n dS + r ∧ ⋅ n dS
(1)
r r r r
− ∫ pr ∧ n dS + ∫ (r ∧ ρ f ) dϑ
∂D D
• Couple exercé sur l’arbre par le fluide :

En projetant (1) dans la direction de l’axe machine, et en considérant le terme de
pesanteur comme globalement nul, on obtient le moment axial des forces (de
pression et de viscosité) exercées par les parois (fixes ou mobiles) sur le fluide.
δ
Ma =
r
δt
∫ D
RVθ dm + ∫
∂D1 ∪∂D2
RVθ dms − ∫
∂D1 ∪∂D2
R(τ ⋅ n )θ dS
Équation d’Euler
Analyse de l’écoulement dans le plan circonférentiel • Application à un tube de courant.
Hypothèses :
sauf le cas (très particulier) d'un écoulement fortement cisaillé, le terme provenant
des tensions visqueuses au sein du fluide est négligeable sur ∂D1 et ∂D2 .
le domaine d'intégration se limite à un tube de courant compris entre R et R+dR :
sur ∂D1 et ∂D2 les valeurs de R et Vθ (considérée comme moyennée dans la direction
θ) pourront donc être considérées comme constantes dans le plan méridien.
le moment cinétique contenu dans D est supposé constant avec le temps
=> les dérivées temporelles dans le repère lié à D sont nulles.
la conservation de la masse impose que le débit sortant de ∂D2 (dms2) soit égal au
débit entrant dans ∂D1 (- dms1).
Le moment axial élémentaire est alors : dM a = (R2Vθ 2 − R1Vθ1 )dms

r r
d’où puissance échangée avec la machine : dP = ω ⋅ dM a = (U 2Vθ 2 − U1Vθ 1 ) dms
dP
L'énergie apportée au fluide par unité de masse est : ∆WT = = U 2Vθ 2 − U1Vθ 1
dms
Dans le cas d'un écoulement adiabatique : ∆h0 = ∆WT + ∆q ≡ ∆WT
∆h0 = U 2Vθ 2 − U1Vθ 1

Équation d’Euler
Analyse de l’écoulement dans le plan circonférentiel Interprétation de l’équation d’Euler
• Rappel : Vθ = Wθ + ω.R
L’équation d’Euler peut donc aussi s’écrire :
∆h0 = (U 2Wθ 2 − U1Wθ 1 ) + (U 22 − U12 )
Travail Travail des forces

aérodynamique de Coriolis
• Travail aérodynamique dans une machine axiale :

Lorsque les variations de rayon sont négligeables le long d’une ligne de courant,
c’est le seul travail existant. On a alors :
∆h0 = U (Wθ 2 − Wθ 1 ) + O
∆h0 = U .Va ( tan β 2 − tan β1 ) (Si Va est constant)
Ce travail est lié à la déflection de l’écoulement au passage de la roue mobile

Équation d’Euler
Analyse de l’écoulement dans le plan circonférentiel • Cas d’un compresseur axial :
Convention de signe :
Vθ2 > 0, Vθ1 > 0 et Vθ2 > Vθ1 ⇒ Vθ 2 − Vθ 1 > 0 (augmentation de l’EC dans le rotor)
ou Wθ2 < 0, Wθ1 < 0 ⇒ Wθ 2 − Wθ 1 > 0 mais |Wθ2| < |Wθ1| (décélération de l’écoulement
relatif dans le rotor)
donc : ∆h0 > 0 cohérence avec la convention de signe thermodynamique.

Équation d’Euler
Analyse de l’écoulement dans le plan circonférentiel • Cas d’une turbine axiale :
Convention de signe :
Vθ2 > 0, Vθ3 < 0 et Vθ3 < Vθ2 ⇒ Vθ 3 − Vθ 2 < 0 (V3<V2 => récupération de
l’EC dans le rotor)
ou Wθ2 < 0, Wθ3 < 0 ⇒ Wθ 3 − Wθ 2 < 0 mais |Wθ2| < |Wθ3| (accélération de l’écoulement
relatif dans le rotor)
donc : ∆h0 < 0 cohérence avec la convention de signe thermodynamique.
Équation d’Euler
Analyse de l’écoulement dans le plan circonférentiel • Remarque sur le gradient de pression longitudinal :
Compresseur : écoulement décéléré (augmentation de la section de passage) :
Gradient de pression ∂p/∂s >0 défavorable pour les couches limites qui se développent
sur les aubages => déflection de l’ordre de 40°
Turbine : écoulement accéléré (diminution de la section de passage) :

Gradient de pression ∂p/∂s <0 favorable pour les couches limites qui se développent
sur les aubages => déflection de l’ordre de 100°
• Machine radiale ou mixte :

Travail de la force de Coriolisr: on montre que le terme U 22 − U12 représente le travail
r
de la force de Coriolis 2.ω ∧ W .
Ce travail n’est pas explicitement lié à la force d’aubage (déflection de l’écoulement)
et ne dépend donc pas du comportement des zones visqueuses se développant sur les
aubages : Il ne dépend que de la variation de R et n’est donc pas limité.
Attention au rôle défavorable de la force de Coriolis sur la nature de la turbulence.
r r r
Travail de la force centrifuge : les forces centrifuges −ω ∧ (ω ∧ r ) travaillent lors
d’un changement de rayon dans le repère relatif. Leur travail se traduit par une
variation de la pression statique (transfert d’énergie ne s’effectuant pas par
l’intermédiaire des aubages et se faisant sans pertes).
Ce travail n’est pas comptabilisé dans l’équation d’Euler, car dans le repère absolu les
forces passant par l’axe ne travaillent pas.
Écoulement en grille d’aubes de compresseur
Analyse de l’écoulement dans le plan circonférentiel Écoulement en grille d’aubes de compresseur axial
• Une fois déterminée la déflection (β2-β1) correspondant au travail désiré, il

reste à :
dessiner les aubages qui permettent cette déflexion avec le minimum de pertes,
déterminer le nombre et la position de ces aubages dans la roue.
• Problème complexe (3D et visqueux) => utilité de données empiriques

provenant de mesures sur des structures plus simples : les grilles d’aubes fixes.
• Les informations obtenues sont intéressantes sur 4 points :

prévision du travail maximum admissible par les aubages (cf critères de charges)
estimation de la déviation de l’écoulement par rapport à la direction du bord de
fuite des aubages
estimation des pertes de pression d’arrêt au passage de la roue
Détermination de l’angle d’incidence optimum (correspondant aux pertes minima)
• Informations utiles dans la phase de dimensionnement (ceci n’exclue pas des

calculs ultérieurs plus sophistiqués), mais transposables qu’à des machines
axiales.
Analyse de l’écoulement dans le plan circonférentiel Définition des paramètres géométriques
c : corde
i : angle d’incidence
g : pas
θ : déflexion
γ : angle de calage
δ: angle de déviation
σ: solidité = c/g
ϕ : cambrure
Analyse de l’écoulement dans le plan circonférentiel Principaux paramètres influant sur les performances de la grille
La nature du profile utilisé (cambrure, répartition d’épaisseur, état de surface)
L’angle d’incidence lié à β1 et au calage γ
Le calage γ : une augmentation de γ :

décroît le guidage car le recouvrement est moindre
augmente β1 à incidence équivalente donc donne plus de travail pour la même déflection θ
(Pour θ =10°, quand on passe de β1=30° à β1=75°, la portance est multiplié par 8)
La solidité σ (=c/g) : une augmentation de σ induit :

un meilleur guidage de l’écoulement
Un niveau de pertes par aubages inférieur du fait de la plus faible circulation par aubage (pic
de vitesse inférieure)
mais : un plus grand nombre d’aubages : d’où un risque de blocage et un niveau de pertes
total qui peut être supérieur.
Le nombre de Reynolds basé sur la corde (qui doit être > 2.5 105)
Le niveau de turbulence dans l’écoulement extérieur.
Le nombre de Mach (attention à la formation de poche supersonique).

Analyse de l’écoulement dans le plan circonférentiel back
Nombe de Reynolds basée sur

la corde et la vitesse relative
en amont de l’aubage doit
être supérieur à 2,5.105 :
W ⋅c
Rec = ≥ 2,5 ⋅105
υ
Pour des valeurs élevées de

ce nombre de Reynolds, les
couches limites sont
principalement en régime
turbulent, et les pertes sont
nettement plus faibles.
Analyse de l’écoulement dans le plan circonférentiel Pertes de pression d’arrêt

• Ces pertes sont dues :
aux couches limites se développant sur les aubages
aux mélange du sillage en aval
éventuellement aux ondes de choc (dans les cas supersoniques)
∆P0 R
• Pertes caractérisées par un coefficient de pertes défini par : ω =
1/ 2.ρ .W12
2H 2
2
θ*  σ  cos β1  3H 2 − 1
avec ω = 2  ⋅ ⋅  ⋅
 c 2 cos β 2 β 3
 cos 2    θ *  σ H2  terme souvent
1 −   ⋅  voisin de l’unité
  c  2 cos β 2 
θ 2* est l’épaisseur de quantité de mouvement du sillage

H2 est le facteur de forme en 2 : H=δ*/θ* (très voisin de 1.1)
.
• Pour calculer les pertes, il faut (si l’on ne dispose pas de calcul de couche
limite) estimer θ 2* / c dans le sillages. Pour cela il existe plusieurs corrélations
avec le niveau de décélération sur l’aubage.
Analyse de l’écoulement dans le plan circonférentiel Critères de charges pour les compresseurs
Nécessitée au moment du dessin d’un critère qui permette de savoir si la CL d’aubage
va ou non supporter le gradient de pression adverse (si le niveau de travail désiré va
pourvoir se faire avec un niveau de pertes acceptable).
Critère facilement utilisable s’il ne prend en compte que des paramètres

aérodynamiques d’entrée et de sortie de la roue.
• Critère de De Haller :
W2
≥ 0.72
W1
Valeur faible par rapport à celle obtenue avec un diffuseur 2D ou axisymétrique, car :
- effets 3D dus aux CL de paroi qui introduisent un blocage

- le gradient réel W2/Wmax est local et est inférieur à W2/W1
Analyse de l’écoulement dans le plan circonférentiel • Facteur de diffusion :
Il vient d’une forme du critère de Buri qui est un critère local de décollement de
couches limites turbulentes, incompressibles et bidimensionnelles, représentant la
tendance à l’accroissement rapide de la couche limite au voisinage du décollement.
Ce facteur a été explicité sous une forme entrée-sortie par Lieblein

Hypothèses : grille d’aubages NACA 65 d’épaisseur maximum 10%,
mais efficacité satisfaisante de ce facteur sur d’autres aubages.
W2 Wθ 2 − Wθ 1
D =1− + ≤ 0.6
W1 2 ⋅ σ ⋅ W1
Relation entre le facteur de diffusion et l’épaisseur de quantité de mouvement des CL
θ* 
  ≈ 0.0804.D − 0.0272.D + 0.0071
2
 c 
Analyse de l’écoulement dans le plan circonférentiel • Facteur de diffusion équivalent :

Le facteur de diffusion cité précédemment ne s’applique qu’au point de
fonctionnement nominal, d’où la nécessité d’un autre critère pour le
fonctionnement hors adaptation. C’est le facteur de diffusion équivalent :
cos β 2  * 1.43 cos 2 β1 

Deq = ⋅ 1.12 + a ⋅ i − i + 0.61 ⋅ tan β 2 − tan β1  ≤ 2.0
cos β1  σ 
a = 0.0117 pour les aubages de la série NACA 65 A10

a = 0.007 pour les aubages de la série C4
i* est l’angle d’incidence optimum
i est l’angle d’incidence réel
Équation de l’équilibre radial simplifié 4ieme partie - Équation de l’équilibre radial simplifié
Introduction

• Équation fondamentale de la dynamique et hypothèses
• Exemple de quelques distributions radiales

Équation fondamentale de la dynamique et hypothèses
Equation fondamentale
r
dV r 1
dt
= f + div σ
ρ () avec σ = τ − pI ( ) effets de pression
effets visqueux
tenseur des contraintes
Hypothèses
1. écoulement permanent
2. équation valable hors zone aubée
3. surfaces de courant cylindriques (machine axiale)
4. écoulement axisymétrique → Aθ = Aθ
5. tension de contrainte visqueuses τij négligeables
6. Forces volumiques négligeables
r
()
1 6 r 5
∂V r r 1 1
+ V .gradV = f + div τ − .grad ( p )
∂t ρ ρ
Repère cylindrique
 ∂VR 1  ∂VR  ∂VR   ∂p 

 ∂R  − Vθ 
 R  ∂θ  ∂z   ∂R 
VR   
 ∂Vθ 1  ∂Vθ  ∂Vθ  V  = 0 + 0 − 1  1 ∂p 
0 +  − V  . .
R  ∂θ
R  θ
 ∂R  ∂z  ρ  R ∂θ 
Vz 
  
 ∂Vz 1 ∂Vz ∂Vz   ∂p 
   ∂z 
 ∂R R ∂θ ∂z 
Projection sur la direction radiale

3 4 3
∂VR 1  ∂VR  ∂VR 1 ∂p
⋅VR +  − Vθ  ⋅ Vθ + ⋅ Vz = − ⋅
∂R R  ∂θ  ∂z ρ ∂R
3 : surface de courant cylindrique VR=0

4 : écoulement axisymétrique : VR(θ)=cte Vθ2 1 ∂p
− = − ⋅
R ρ ∂R
Équation de l’équilibre radiale simplifiée

1
dh − dp = Tds (Équation de Gibbs)
ρ 2
∂s ∂h0 ∂V
1 ∂p ∂s ∂h 2 = T ∂s − ∂h0 + 1 ∂ V 2 + V 2
−
ρ ∂R
=T −
∂R ∂R
=T −
∂R ∂R
+
∂R ∂R ∂R 2 ∂R
θ z ( )
∂s ∂h0 ∂V ∂V
=T − + Vθ ⋅ θ + Vz ⋅ z
∂R ∂R ∂R ∂R
1 ∂p Vθ2 ∂s ∂h0 ∂V ∂V
− ⋅ =− =T − + Vθ ⋅ θ + Vz ⋅ z
ρ ∂R R ∂R ∂R ∂R ∂R
Vθ ∂R Vθ ∂V ∂Vz ∂h0 ∂s
⇔ ⋅ Vθ ⋅ + ⋅ R ⋅ θ + Vz ⋅ = −T
R ∂R R ∂R ∂R ∂R ∂R
Vθ ∂RVθ ∂V ∂h ∂s
⋅ + Vz ⋅ z = 0 −T
R ∂R ∂R ∂R ∂R
( )
2
1 ∂ RVθ
Négligeable si pas de
∂Vz2 ∂h
E.R.S. ⋅ + = 2⋅ 0 gradient d’entropie
R 2
∂R ∂R ∂R selon R.
Bilan équations / inconnues

• Considérons l'écoulement au passage d'une roue. Si nous faisons les bilans du

nombre d'inconnues du problème et du nombre d'équations disponibles, nous
avons :
6 inconnues :
- les répartitions radiales de h0, Vz et Vθ dans une section à l'amont de la roue
- les répartitions radiales de h0, Vz et Vθ dans une section à l'aval de la roue
3 équations :
- l'équation d'équilibre radial simplifiée dans la section à l'amont de la roue
- l'équation d'équilibre radial simplifiée dans la section à l'aval de la roue
- l'équation d'Euler pour chaque rayon au passage de la roue
• Il y a donc la possibilité de choisir librement 3 répartitions.

Quelques lois d’évolutions classiques
1 2
°
m
i
Le tourbillon libre
RVθ = Cte
• Conditions amont h01 = Cte et R1Vθ1 = Cte
• Tourbillon libre à l’aval R2Vθ 2 = Cte = k
=> h02 , Vz1 et Vz2 constants.

Le tourbillon forcé
Vθ = k R
=> (
Vz2 = Vzi2 + 2⋅ k ⋅ (ω − k ) R 2 − Ri2 )
Vzi est calculé pour satisfaire la conservation du débit global.
L’angle absolu constant
α = Cte (tan α = Vθ Vz )
sin 2 α
R 
=> Vz = Vzi ⋅  i 
R
k2
Vθ = k1 R n ± Cours de turbomachine – Xavier Ottavy (CNRS – Laboratoire de Mécanique des Fluides et d’Acoustique UMR 5509 - ECL)
R
Quelques autres lois classiques

• Le vortex général
k2
Vθ = k1 R n ±
R
• L’étage exponentiel
k2
Vθ = k1 ±
R
• Le degré de réaction constant
Λ = Cte
View publication stats
Téléchargement du cours sur :

Cours de Turbomachine A Fluide Compressible

Transféré par

Informations du documentcliquez pour développer les informations du document

Droits d'auteur :

Formats disponibles

Cours de Turbomachine A Fluide Compressible

Transféré par

Informations du document

Description originale:

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Cours de Turbomachine A Fluide Compressible

Transféré par

Droits d'auteur :

Formats disponibles

See discussions, stats, and author profiles for this publication at: https://www.researchgate.

Analyse thermodynamique monodimensionnelle

Analyse de l’écoulement dans le plan circonférentiel

Équation de l’équilibre radial simplifié

Cas d’application : dessin d’un étage de compresseur axial

1ière partie - Introduction

• Approches 1D, 2D, 2.5D et 3D

Analyse thermodynamique monodimensionnelle

Analyse de l’écoulement dans le plan circonférentiel

Équation de l’équilibre radial simplifié

Cas d’application : dessin d’un étage de compresseur axial

Une turbomachine est une machine tournante qui réalise un transfert

• une récupération de l’énergie du fluide sur l’arbre de la machine

• une augmentation de l’énergie du fluide par fourniture d’énergie

Fonctions et domaines d’utilisation des turbomachines

Récupération de l’énergie d’un fluide (turbines)

Compression de gaz (compresseurs)

Production d’énergie mécanique à partir d’une source de chaleur

Production d’énergie mécanique : entraînement d’hélice de bateau, d’avion

Propulsion par réaction

Introduction PW4156 - Pratt & Whitney (epower-propulsion.com)

Introduction CFM56 – Snecma Moteurs (epower-propulsion.com)

Introduction BR715 – BMW/Rolls-Royce (epower-propulsion.com)

Introduction GP7000 – EA (GE / Pratt & Whitney) (epower-propulsion.com)

Introduction PEGASUS – Rolls-Royce (epower-propulsion.com)

Introduction F404 – General Electric (epower-propulsion.com)

Notion d’étage – Échanges d’énergies

Géométries des turbomachines

• Les machines centrifuges : le fluide sort approximativement dans un plan

• Les machines mixtes

Notion d’étage - échanges d’énergie

Rôle : assurer le transfert d’énergie entre l’arbre de la machine et le fluide en

Énergie de pression : une turbomachine échange nécessairement de l’énergie de

Énergie cinétique : une turbomachine échange nécessairement de l’énergie

Énergie calorifique : il n’y a pas d’énergie calorifique directement échangée entre

Rôle : modifier la forme d’énergie (énergie cinétique en pression, ou

Redresseur de compresseur axial :

Distributeur de turbine axiale :

Diffuseur de compresseur centrifuge :

• Compresseurs et turbines sont en général calculés pour un point de

• Il est cependant intéressant de connaître le comportement de la machine à

• Actuellement, les recherches sont largement orientées sur l’extension des

Courbes caractéristiques des compresseurs

• Exemple de déclenchement d’un décollement tournant

Courbes caractéristiques des turbines

• Asymptote commune, indépendante de la vitesse de rotation => blocage sonique dans

Approche 1D, 2D, 2.5D et 3D

Définition de surfaces «méridienne» et «aubes-à-aubes»

2ième partie - Analyse thermodynamique monodimensionnelle

Analyse thermodynamique monodimensionnelle

• Bilan des différentes contributions

• Équation de l’entropie (équations de Gibbs)

• Travail mécanique et travail utile

• Intérêt du diagramme entropique et enthalpique

Analyse de l’écoulement dans le plan circonférentiel

Équation de l’équilibre radial simplifié

Équations de conservation de base