Dans ce module on ne considère que des écoulements

compressibles internes dans des conduites rectilignes à section

droite à forme constante.
L’écoulement est pris à l’état permanent et on se limite à un fluide
parfait: non visqueux et non conducteur de la chaleur.
Les écoulements compressibles sont caractérisés par une
variation de la masse volumique ρ. Par conséquent, en plus des
équations de conservation de la masse et de la quantité de
mouvement, on doit considérer l’équation de conservation
d’énergie.
À ces équations, on ajoute l’équation d’état des gaz parfaits qui
permet de faire le couplage entre la pression et l’énergie (ou la
température).
Puisque la masse volumique est variable, l’ensemble des
propriétés thermodynamiques est affecté. On fait donc un court
rappel sur ce sujet.
 9.1 9.2 9.3 9.4 9.5 9.6

9.1 Thermodynamique

9.2 La vitesse du son

9.3 Écoulement adiabatique et isentropique

9.4 Écoulement isentropique avec changement de section

9.5 Choc normal

9.6 Tuyères convergente et convergente-divergente

Ecoulements incompressibles Ecoulements compressibles
Conservation de la masse Conservation de la masse
Conservation de la quantité de Conservation de la quantité de
mouvement
mouvement
Conservation de l’énergie
5 équations: 3 pour la quantité de
4 équations: 3 pour la quantité de mouvement + une pour la
mouvement + une pour la conservation de la masse + une pour
conservation de la masse la conservation de l’énergie
L’équation d’état qui lie la pression à
la masse volumique
 vitesse de l'objet (de l'écoulement)
Nombre de Mach =
vitesse du son

3.0< Ma
Écoulement hypersonique

1.2< Ma <3.0
Écoulement supersonique
0.8< Ma <1.2
Écoulement transsonique

0.3< Ma <0.8
Écoulement subsonique

Ma <0.3
Écoulement incompressible
 Rappel

p = ρ RT
R est le rapport entre la costante universelle Ru=8314.3 J/(kmol K)
et la masse moléculaire du gaz M : R=R u/M

 ∂u  du
=cv = T2

 ∂T v dT ∫ cv dT
u2 − u1 =
T1
 Rappel
p
h= u +
ρ

 ∂h  dh T2
=c p = 
 ∂T  p dT ∫ c p dT
h2 − h1 =
T1

kR R cp
c=
p − cv R; =cp ;=cv ;=k
k −1 k −1 cv
 Rappel

f
δ Qrev δ Qrev
∆s =∫ ds =
T T
i

δ q = δ Qrev

δq
ds = δ q = Td s
T

Le premier principe 𝑑𝑑𝑑𝑑 = 𝛿𝛿𝑞𝑞 − 𝑝𝑝𝑝𝑝𝑝𝑝 s’écrit donc

 Rappel

= Tds − p d v
du d=
h Tds + vd p

d’où,

Tds =du + pdv =dh − vdp

Formulation du premier principe en fonction de variables d'état

uniquement!
 Rappel

Sous la forme Tds, l’équation du premier principe s’écrit

= dh − vdP
Tds
dh = c p dT v = RT / P Pour un gaz parfait

dT dP
=
ds c p −R
T P
 Rappel

(Cp ≅ constante et Cv ≅ constante)

T2 P2
s2 −
= s1 c p ln − R ln
T1 P1

T2 V2
s2 −=
s1 cv ln − R ln
T1 V1
 Rappel

R
T2 P2 T2  p2 
s2 −
= s1 c p ln − R ln
Cp
= 
T1 P1 T1  p1 

k −1

cp =
kR  T2   p2  k Relation isentropique
k −1   =  lorsque cp=cste.
 T1  s =const .  p1 
Définition Relation Gaz Parfait
Relation pression-densite-temperature
= p ( ρ , T ) p ρ RT
f=
Énergie interne par unité de masse uˆ uˆ (T )
Enthalpie par unité de masse uˆ + p / ρ
h= h= uˆ + RT
∂uˆ duˆ
Chaleur spécifique à volume constant cv = cv =
∂T v dT
dh duˆ
∂h =
cp = +R
Chaleur spécifique à pression constante cp = dT dT
∂T p c p − cv =R
1
Volume spécifique
ρ
 = c p − cv
R Air
k = 1.4
ft 2 m2
= =
R 1716 287 2
c  1
s2  R s K
=R cv  p − 1  cv = R ft 2 m2
 cv  ⇒ k −1 = =
cv 4293
s2  R
718 2
 s K
cp  k ft 2 m2
k= c = R = =
c p 6009 1005
 p
k −1 s2  R s2 K
cv 

La variation d’énergie interne et d’enthalpie d’un

gaz parfait à chaleur spécifique constante
 =s1 s=
2, c p cste.

1
cv R cv / R
 T2   ρ2   ρ 2   T2  T  k −1
  =   ⇒  =   =  2
 T1   ρ1   ρ1   T1   T1 
k
cp R cp / R
 T2   p2   p2   T2  T  k −1
  =   ⇒  =   =  2
 T1   p1   p1   T1   T1 
k
k
 p2   T2   ρ2 
k −1
=
    =  
 p1   T1   ρ1 
k
 ρ 2  p1 p
  = 1 ⇒ = cst
ρ
 1 2 p ρ k
Exemple 9.1

De l’argon s’écoule dans un tube. La condition initiale est:

𝑝𝑝1 = 1.7 𝑀𝑀𝑀𝑀𝑀𝑀, 𝜌𝜌1 = 18 𝑘𝑘𝑘𝑘⁄𝑚𝑚3 , tandis que l’état final est 𝑝𝑝2 = 248 𝑘𝑘𝑘𝑘𝑘𝑘𝑎𝑎, 𝑇𝑇2 = 400𝐾𝐾

Estimer (a) la température initiale, (b) la masse volumique finale, (c) le

changement d’enthalpie, et (d) le changement d’entropie

Solution

Les propriétés =
de l’argon =
R 208 (
m 2 / s 2 K et k 1.67)
p1
a) Calcul de T1: =
T1 = 454 K
ρ1 R
p2
b) Calcul de ρ2: ρ2
= = 2.98 kg / m3
T2 R

Version adaptée du livre: Fluid Mechanics, par Frank M. White, 7e édition

Exemple 9.1 (suite) T1 = 454 K
T2 = 400 K
Solution

Les propriétés de
= l’argon =
R 208 ( )
m 2 / s 2 K et k 1.67

c p (T2 − T1 )
kR
c) Calcul de Δh: h2 − h= cp =
k −1
1

∆h =−28000 J / kg

T  p  p1 = 1.7 MPa
d) Calcul de Δs: =
s2 − s1 c p ln  2  − R ln  2 
 T1   p1  p2 = 288kPa
∆s =334 m 2 / ( s 2 K )
 9.1 9.2 9.3 9.4 9.5 9.6

9.1 Thermodynamique

9.2 La vitesse du son

9.3 Écoulement adiabatique et isentropique

9.4 Écoulement isentropique avec changement de section

9.5 Choc normal

9.6 Tuyères convergente et convergente-divergente

 Idée

Une forme permettant de caractériser la compressibilité d’un

fluide, est d’évaluer la vitesse avec laquelle se propagent des
petites perturbations au sein de celui-ci. Cette vitesse correspond
à la vitesse du son.
Pour évaluer la vitesse du son, on imagine une perturbation de
pression infinitésimale dp,créée par le déplacement d’un piston
avec une vitesse infinitésimale dc à l’extrémité d’un tube.
La perturbation se propage dans le tube sous la forme d’une onde
plane avec une vitesse c.
 ρ
D T

Vitesse c
 ρ
D T

Vitesse c
 ρ
D T

Vitesse c
 ρ ρ+Δρ ρ
D P p + Δp T
T T +ΔT
v=0 Δv

Vitesse c
 Vitesse du son

ρ ρ + Δρ
ρ
P p + Δp
D T
T T+ΔT
v=c v=c-Δv

Vitesse c

∆ρ
∆V =
C
ρ + ∆ρ
ρ AC= ( ρ + ∆ρ ) A ( C − ∆V )
 Vitesse du son

ρ ρ + Δρ
P p + Δp ρ
D T+ΔT
T T
v=c v=c-Δv

ρ AC ( −∆V ) = −∆pA ∆p= ρ C ∆V

m (Vout − Vin=
) pA − ( p + ∆p ) A
 Vitesse du son

ρ ρ +Δρ
P p + Δp ρ
D T+ΔT
T T
v=c v=c-Δv

∆ρ
∆p= ρ C ∆V ∆p =ρ C 2
ρ + ∆ρ
∆ρ
∆V =
C
ρ + ∆ρ ∂p ∆p ρ
= lim= C 2 lim = C 2
∂ρ ∆ρ →0 ∆ρ ∆ρ → 0 ρ + ∆ρ
La vitesse du son correspond à la vitesse de ∂p
propagation d’une perturbation dans un C= a=
fluide. ∂ρ s

Si les gradients de vitesse et de température

produits par l’onde deviennent si petits que
les contraintes de cisaillement et le transfert
de chaleur sont négligeables, il s’agit d’une
onde acoustique.

Par conséquent, le mouvement du fluide lié à

une onde acoustique est isentropique.
 Gaz parfait

Avec l'équation d'état et la relation isentropique, on a

p = ρ RT
p cnst × ρ
= k

 ∂p  ρk p  ∂p  p
  = const × k × ρ k −1
= k × const × = k  =
 =
RT
 ∂ρ s ρ ρ  ∂ρ T ρ

 ∂p   ∂p 
de sorte que   = k  
 ∂ρ s  ∂ρ T
1/2
 ∂p 
a =   a = kRT
 ∂ρ s
 Interprétation alternative

• Soit un milieu fluide non visqueux avec les propriétés p0 ,ρ0,s0 au

repos
• On considère une perturbation de pression infinitésimal et on
effectue un développement en série de Taylor
 ∂p  1  ∂2 p 
∆p = p − p0 =   ( ρ − ρ0 ) +  2  ( ρ − ρ0 ) + ...
2

 ∂ρ  s 2  ∂ρ  s

• On appelle vitesse du son la grandeur a telle que

 ∂p 
a=  
 ∂ρ  s
Exemple 9.2

On a du monoxide de carbone à: 𝑝𝑝 = 200 𝑘𝑘𝑘𝑘𝑘𝑘, 𝑇𝑇 = 3000𝐶𝐶 . On doit estimer la

vitesse du son.

Solution

Les propriétés du
= CO sont R 297
= ( )
m 2 / s 2 K et k 1.4

T = 300 + 273 = 573 K

Vitesse du son =a =
kRT 488 m / s

Version adaptée du livre: Fluid Mechanics, par Frank M. White, 7e édition

 9.1 9.2 9.3
9.3 9.4 9.5 9.6

9.1 Thermodynamique

9.2 La vitesse du son

9.3 Écoulement adiabatique et isentropique

9.4 Écoulement isentropique avec changement de section

9.5 Choc normal

9.6 Tuyères convergente et convergente-divergente

Lorsqu’on tient compte de l’énergie, on découvre une couche
limite thermique semblable à la couche limite dynamique.
L’épaisseur de la couche limite thermique, est en général
différente de celui de la couche limite dynamique. Tout comme
l'épaisseur de la couche limite hydrodynamique, elle est liée à la
vitesse de l'écoulement.
Le rapport entre les épaisseurs des deux couches limites est
donné par le paramètre adimensionnel connu sous le nom de
nombre de Prandtl.
Le fluide courant étant parfait, l’intérêt est porté iniquement sur
les phénomènes à l’extérieur ou sans une couche limite.
On suppose un écoulement compressible permanent au
voisinage d’une paroi adiabatique.
Entre deux points, l’équation unidimensionnelle d’énergie à
l’état permanent peut s’écrire:
V12 V22
h1 + + gz1 =h2 + + gz2 − q + wv
2 2

Si l’on néglige la variation d’énergie potentielle et en absence

d’un travail à l’arbre, on trouve:
 V12 V22 V12 V22
h1 + + gz1 =h2 + + gz2 − q + wv h1 + =h2 + =cst
2 2 2 2

Grandeur d’arrêt: on appelle état d’arrêt, totale, ou de stagnation,

l’état que prends toute variable de l’écoulement si on l’amenait au
repos de manière adiabatique et réversible, donc isentropique
On note par h0 l’enthalpie d’arrêt, totale ou de stagnation. Elle
représente la combinaison de l’énergie cinétique avec l’enthalpie
V2
h0= h +
2
 h0
V2
h0= h +
2
La valeur atteinte par h lorsque le fluide est emmené au repos selon un
processus adiabatique réversible.

h0 T0

∫ dh = ∫ c
h T
p dT h0 −=
h c p (T0 − T )

V2
T0= T +
2c p
 Remarque

V2 V2
h0= h + T0= T +
2 2c p

Il est important de noter que ces relations sont valables dans tout
écoulement adiabatique. Donc, aussi valides dans un écoulement
avec génération d’entropie, du moment qu’il n’y a pas transfert de
chaleur.
Ceci est particulièrement important dans un écoulement
présentant une onde de choc, à travers laquelle l’entropie
augmente, mais la température totale T0 demeure constante
 h0

v2
h0 = h +
2

h0
  V22 
gz1.  
 2 .

gz 2
  V22 
 
T01  2C 
 p .

T2
 V
Ma =
a
V2 T0 k −1
a = kRT T0= T + = 1+ Ma 2
2c p T 2
kR
cp =
k −1
 k / k −1 1/ k −1
p0  k − 1  ρ0  k − 1 2 
= 1 + Ma 2  = 1 + Ma 
p  2  ρ  2 

T0 k −1
= 1+ Ma 2
T 2
Pour un processus isentropique, l’entropie de l’état de stagnation
est la même que celle du fluide au repos
La température de stagnation Τ0 est supérieure à la température
statique Τ
La différence entre ces deux températures correspond à l’énergie
cinétique du fluide
V2
T0 − T =
2c p

La vitesse est maximale si T=0

Puisque
a =
kRT et a0 kRT0

T0 k −1 a0 k −1 T0
= 1+ Ma 2 =
1+ Ma 2 =
T 2 a 2 T
  k −1   k −1 
Pour Ma=1 =
T0 T 1 + Ma 2  =
T0 T * 1 + 
 2   2 
−1
T* température critique T *  k +1 
= 
T0  2 
V* vitesse critique 2kRT0
=V* = a*
k +1
a* vitesse du son crítique

Par définition, les quantités critiques correspondent à une section de

l’écoulement où Ma=1
 k / k −1 1/ k −1
∗
p  2  T∗  2  ρ∗  2 
= =  =

p0  k + 1  T0  k + 1  ρ0  k + 1 

p = ρRT
 Nous regarderons la
La thermodynamique
projection enthalpie-entropie
est 3D

Transformation irreversible

Transformation à s=cte
 p0
p0 − r
Si le processus est isentropique

V2
h+ =
h0 h0
2
Pour un processus réel:
2
V Même h0, différent p0
2

s
 = dh − vdp
Tds ds = 0 dh = dp / ρ

V2
h0= h + dh + VdV =
0
2

dp V 2 dp / ρ + VdV =
∫ρ+ 2 = cst 0
 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 

p0 k k
=
1 + Ma 2 + Ma 4 + O( Ma 6 )
p 2 8
 1 2 M2 
=
p0 − p pkM 1 + + O( M 4 ) 
2  4 

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

Pour M < 0.3

2 4 
2
1 M 1
p=
0 − p ρV 1 + + O( M )  p0 − p ≈ ρV 2
2  4  2
 14
c p cv = 1.4
12

10

8
ρ0 ρ
6

0
0 0.5 1 1.5 2 2.5 3
Ma
 Air (k=1.4) Valeurs critiques:
Ma=1, air (k=1.4)
T0 k −1 T0 T* 2
= 1+ Ma 2 = 1 + 0.2 Ma 2 = = 0.8333
T 2 T T0 k +1
a0 k −1 T0 a0 a* 2
=
1+ Ma 2 = = 1 + 0.2 Ma 2
= = 0.9129
a 2 T a a0 k +1
k k
k −1
p0   p0
(1 + 0.2Ma )  2 
k −1 *
3.5 p k −1
= 1 + Ma 2  = 2
= =  0.5283
p  2  p p0  k + 1 
ρ0
1 1
ρ0  k −1 (1 + 0.2Ma ) ρ *  2  k −1
2.5
 k −1
= 2
= =
= 1 +
ρ 0  k + 1 
Ma 2  ρ 0.6339
ρ  2 
V2  T0    ρ  2/5    p  2/7 
c pT + =
c pT0 Ma = 5  − 1 = 5  
2

0
 − 1 = 5   0  − 1
   V=
*
a=
*
kRT *
2 T   ρ    p  
 Exemple 9.3

De l’air s’écoule dans un tube sans transfert de chaleur (adiabatique). Au point 1

𝑝𝑝1 = 170 𝑘𝑘𝑘𝑘𝑘𝑘, 𝑇𝑇1 = 320𝐾𝐾, 𝑉𝑉1 = 240 𝑚𝑚/𝑠𝑠 . Au point 2 𝑝𝑝2 = 145 𝑘𝑘𝑘𝑘𝑎𝑎, 𝑉𝑉2 = 290 𝑚𝑚/𝑠𝑠

Estimer (a) T01, (b) p01 , (c) ρ01, (d) Ma1, (e) Vmax, (f) V* , et (g) p02.

Solution

Les propriétés de l’air =R 287 m 2 / ( s 2 K=

) , k 1.4,
= c p 1005 m 2 / ( s 2 K )

V12
a) Calcul de T01: T01 =T1 + =349 K T02 = T01
2c p
=
Vmax =
2c pT01 837 m / s

a* 2 2k
= V=
*
a=
*
RT01= 342 m / s
a0 k +1 k +1
Version adaptée du livre: Fluid Mechanics, par Frank M. White, 7e édition
Exemple 9.3 (suite) k
T0 k −1 p0  k −1  k −1
= 1+ Ma 2 = 1 + Ma 2 
T 2 p  2 
Solution

T 
p1 (1 + 0.2 Ma12 )
3.5
=
Ma1 5  01 −=
1 0.67 p01 = =
230kPa
 T1 
p01
ρ 01
= = 2.29 kg / m3
RT01
k
V22  T02 
Calcul de T2 ,p2
k −1
T2 =T02 − =307 K =p02 p=
2  211 kPa
2c p  2 
T

Forme alternative T 
=
Ma2 5  02 −=
1 0.83
 T2 
p2 (1 + 0.2 Ma )
3.5
p02 = 2
2 =
212kPa
p02 est 8% de moins que p01 car l’écoulement n’est pas isentropique entre 1 et 2
 9.1 9.2 9.3 9.4 9.5 9.6

9.1 Thermodynamique

9.2 La vitesse du son

9.3 Écoulement isentropique

9.4 Écoulement isentropique avec changement de section

9.5 Choc normal

9.6 Tuyères convergente et convergente-divergente

On examine l’effet de la variation de section d’une conduite
sur un écoulement compressible. Le problème est
unidimensionnel et les variables ne dépendent que de la
coordonnée x normale à la section, c’est-à-dire, le long de
l’axe de la conduite. À l’état stationnaire, 𝑚𝑚̇ = 𝜌𝜌𝑉𝑉𝑉𝑉 = 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐.,
de sorte que 𝑑𝑑 𝜌𝜌𝑉𝑉𝑉𝑉 = 0, donc

d ρ dA dV
+ + =
0
ρ A V
L’écoulement étant isentropique
dp
+ VdV =
0
ρ
et à partir de
 ∂p  dρ dp
=dp =  d ρ a 2
dρ =
 ∂ρ  s ρ ρ a2
on trouve dρ V dρ Ma 2
=
− 2 dV =
0 =
− dV =
0
ρ a ρ V
 dρ Ma 2 d ρ dA dV
=
− dV =
0 + + =
0
ρ V ρ A V
dA dV dp
=
( Ma 2 − 1) =
−( Ma 2 − 1)
A V ρV 2
Ce résultat révèle, selon que l’écoulent soit subsonique ou
supersonique, la forme à donner à la conduite pour effectuer une
compression ou une détente. Pour le cas d’une tuyère, conçue
pour accélérer un fluide, elle doit être convergente (dA < 0)⤇dp < 0,
si l’écoulement est subsonique (Ma < 1) et divergente (dA > 0) ⤇
dp < 0, si l’écoulement est supersonique( Ma > 1).
On peut illustrer des cas décrits par l’équation “aire-vitesse-Mach”
dA dV
= ( Ma 2 − 1)
A V

A diminue A augmente
M <1 M >1
V augmente V augmente
Subsonique Supersonique p diminue
p diminue

A augmente A diminue
M <1 M >1
V diminue V diminue
Subsonique Supersonique
p augmente p augmente
 dA dV
= ( Ma − 1)
2

A V
dA -- -- + +
M <1 >1 <1 >1
dV/V + -- -- +
dp/p -- + + --

En régime subsonique la section et la vitesse varient en sens

opposés, tandis qu’en régime supersonique, la section et la
vitesse varient dans le même sens.
On voit que l’équation précédente donne une accélération infinie
au point sonique (Ma=1). Mathématiquement, il s’agit d’une
singularité!
dV dA 1
=
V A Ma 2 − 1

La seule solution est dA=0 à Ma=1. C’est-à dire pour un minimum

ou un maximum A.
Pour la combinaison d’un divergent suivi d’un convergent, donc a
section maximale, le divergent entraine la diminution du nombre
de Mach lors d’un écoulement subsonique, ou fait augmenter le
nombre de Mach, s’il s’agit d’un écoulement supersonique. Il ne
peut donc pas y avoir un écoulement sonique pour cette
combinaison à section maximale.
Dans le cas d’une combinaison a section minimale, soit d’un
convergent suivi d’un divergent, l’effet contraire se produit sur
l’écoulement et un écoulement sonique peut avoir lieu au col.
La seule configuration physique qui permet d’accélérer un
écoulement est le col
 Conditions de Conditions soniques Conditions de
stagnation stagnation
V = V0 = 0 V = V0 = 0
V = V *= a*
h0, T0, ρ0, p0, s0

Réservoir h*,T* ρ*, p*, s* Réservoir

Grandeurs de stagnation ou d’arrêt Grandeurs soniques

Une propriété ou grandeur de stagnation est Une propriété ou grandeur sonique est définie en
définie en un point de l’écoulement où la vitesse un point où la vitesse est égale à la vitesse locale
est nulle du son
Notation h0, T0, ρ0, p0, s0 Notation h*, T*, ρ*, p*, s*
Caractéristique V = V0 = 0 Caractéristique V = V* = a*
Il est pratique de relier le nombre de Mach à une section donnée,
à l’aire de cette même section. Ceci peut s’obtenir en considérant
la conservation du débit massique 𝒎𝒎̇ le long de la conduite

Au col, M = 1 (conditions soniques) Reste de la tuyère, M ≠ 1

ρ*, V*, A* m ρ=
= VA ρ *V * A* ρ, V, A

A ρ * V * ρ * a* ρ * ρ 0 a*
= = =
A *
ρ V ρ V ρ0 ρ V
 1
ρ *
 2  k −1
=
A ρ * V * ρ * a* ρ * ρ 0 a* ρ0  k + 1 
= = =
ρ V ρ V ρ0 ρ V
1
*
A ρ0  k − 1 2  k −1
= 1+
ρ 
Ma 
2 

k +1
 k −1 2 
2( k −1)
1 + Ma
A 1  2 
=  
A*
Ma  1 + k − 1  Aire adimensionnelle
 2  référée à l’aire critique
Branche subsonique Branche supersonique

A
A
A*
 Au col, M = 1 (conditions soniques) Reste de la tuyère, M ≠ 1

ρ*, V*, A* ρ, V, A
1 γ 1
V A ρ ∗ ∗
ρ  p  γ  2  γ −1  p  γ
= =  ∗  =    
V∗ A ρ ρ  p   γ + 1   p0 
∗

k +1
−
 k −1 2 
2( k −1)
k
−
1
 1 + Ma 
V  2  k −1  k −1 2
k −1
2
+ Ma 
    k −1 
1 Ma
V∗  k +1  2   1+  Vitesse adimensionnelle
 2  référée à la vitesse critique
Toutes les variables sans dimension T0 T
A A*
peuvent être exprimées en fonction
d’une seule: le nombre de Mach Ma
V V*
Le rapport A/A* est minimum si Ma=1
Un même rapport A/A* peut
correspondre à un écoulement
subsonique et à un écoulement
supersonique.
1

p p0

1 Ma
 Subsonique Supersonique
M A/A* p/p 0 ρ/ρ 0 Τ/Τ 0 V/V* M A/A* p/p 0 ρ/ρ 0 Τ/Τ 0 V/V*
0 Inf 1.0000 1.0000 1.0000 0 1.0000 1.0000 0.5283 0.6339 0.8333 1.0000
0.0500 11.5914 0.9983 0.9988 0.9995 0.0548 1.0500 1.0020 0.4979 0.6077 0.8193 1.0411
0.1000 5.8218 0.9930 0.9950 0.9980 0.1094 1.1000 1.0079 0.4684 0.5817 0.8052 1.0812
0.1500 3.9103 0.9844 0.9888 0.9955 0.1639 1.1500 1.0175 0.4398 0.5562 0.7908 1.1203
0.2000 2.9635 0.9725 0.9803 0.9921 0.2182 1.2000 1.0304 0.4124 0.5311 0.7764 1.1583
0.2500 2.4027 0.9575 0.9694 0.9877 0.2722 1.2500 1.0468 0.3861 0.5067 0.7619 1.1952
0.3000 2.0351 0.9395 0.9564 0.9823 0.3257 1.3000 1.0663 0.3609 0.4829 0.7474 1.2311
0.3500 1.7780 0.9188 0.9413 0.9761 0.3788 1.3500 1.0890 0.3370 0.4598 0.7329 1.2660
0.4000 1.5901 0.8956 0.9243 0.9690 0.4313 1.4000 1.1149 0.3142 0.4374 0.7184 1.2999
0.4500 1.4487 0.8703 0.9055 0.9611 0.4833 1.4500 1.1440 0.2927 0.4158 0.7040 1.3327
0.5000 1.3398 0.8430 0.8852 0.9524 0.5345 1.5000 1.1762 0.2724 0.3950 0.6897 1.3646
0.5500 1.2549 0.8142 0.8634 0.9430 0.5851 1.5500 1.2116 0.2533 0.3750 0.6754 1.3955
0.6000 1.1882 0.7840 0.8405 0.9328 0.6348 1.6000 1.2502 0.2353 0.3557 0.6614 1.4254
0.6500 1.1356 0.7528 0.8164 0.9221 0.6837 1.6500 1.2922 0.2184 0.3373 0.6475 1.4544
0.7000 1.0944 0.7209 0.7916 0.9107 0.7318 1.7000 1.3376 0.2026 0.3197 0.6337 1.4825
0.7500 1.0624 0.6886 0.7660 0.8989 0.7789 1.7500 1.3865 0.1878 0.3029 0.6202 1.5097
0.8000 1.0382 0.6560 0.7400 0.8865 0.8251 1.8000 1.4390 0.1740 0.2868 0.6068 1.5360
0.8500 1.0207 0.6235 0.7136 0.8737 0.8704 1.8500 1.4952 0.1612 0.2715 0.5936 1.5614
0.9000 1.0089 0.5913 0.6870 0.8606 0.9146 1.9000 1.5553 0.1492 0.2570 0.5807 1.5861
0.9500 1.0021 0.5595 0.6604 0.8471 0.9578 1.9500 1.6193 0.1381 0.2432 0.5680 1.6099
2.0000 1.6875 0.1278 0.2300 0.5556 1.6330
 Formule très utilisée

P Pk P
m ρ=
= VA VA =m = VA kaMaA
RT kRT kRT

P P P
m = k kRT MaA kMaA m = kMaA
kRT kRT kRT
 Formule très utilisée

P T0  k − 1  p0  k − 1
k 1

m = kMaA = 1 + Ma 2  = 1 +
 k −1
Ma 2 
ρ0  k − 1 2 
= 1+
k −1

ρ 
T  Ma 
kRT 2  p  2  2 

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

Si dans une conduite l’écoulement est sonique au col, alors le débit

massique est maximal: on dit que l’écoulement est étranglé (choked).
 Écoulement étranglé

La formule générale du calcul du débit, se simplifie lorsque

l’écoulement est étranglé (Ma=1 dans la formule précédente):
k +1 k +1
 2  2( k −1) p0  2  2( k −1)
k= ρ0 ( RT0 ) A k 
1/2
m max 1/2
 
*
 A*
 k + 1  a0  k + 1 
Si le fluide est de l’air (k=1.4), cette dernière devient:

Dans ces relations: A*=aire du col= aire critique

Exemple 9.4

De l’air s’écoule isentropiquement dans une conduite. Au point 1 on a

𝑝𝑝1 = 500 𝑘𝑘𝑘𝑘𝑘𝑘, 𝑇𝑇1 = 470𝐾𝐾, 𝑉𝑉1 = 180 𝑚𝑚/𝑠𝑠 𝐴𝐴1 = 0.05𝑚𝑚2 .
a)Estimez 𝑝𝑝01 , 𝑇𝑇01 , 𝑀𝑀𝑀𝑀1 , 𝐴𝐴∗ , 𝑚𝑚̇ , b) Au point 2E(subsonique) 𝐴𝐴 = 0.036𝑚𝑚2 , on doit
trouver 𝑝𝑝2 , 𝑀𝑀𝑀𝑀2 , c) Au point 2F(supersonique) 𝐴𝐴 = 0.036𝑚𝑚2 , on doit trouver
𝑝𝑝2 , 𝑀𝑀𝑀𝑀2

2𝐸𝐸 2𝐹𝐹

T1 = 470 K
p1 = 500kPa
V1 = 180m / s

𝐴𝐴 = 0.05𝑚𝑚2 𝐴𝐴 = 0.036𝑚𝑚2 𝐴𝐴 = 0.036𝑚𝑚2

Version adaptée du livre: Fluid Mechanics, par Frank M. White, 7e édition

Exemple 9.4 (suite I) 𝑝𝑝1 = 500 𝑘𝑘𝑘𝑘𝑘𝑘, 𝑇𝑇1 = 470𝐾𝐾,
𝑉𝑉1 = 180 𝑚𝑚/𝑠𝑠 𝐴𝐴1 = 0.05𝑚𝑚2 .

Solution

Les propriétés de l’air =R 287 m 2 / ( s 2 K=

) , k 1.4,
= c p 1005 m 2 / ( s 2 K )

a) Calcul de T01: T01 =T1 +

V12
=486 K p1 (1 + 0.2 Ma
p01 = 1)
2 3.5
=
563kPa
2c p
=a1 =
kRT1 435 m / s ρ1
=
RT1
= 3.71kg / m3
p1
V1
Ma=
1 = 0.414
a1
Exemple 9.4 (suite II) 𝑝𝑝1 = 500 𝑘𝑘𝑘𝑘𝑘𝑘, 𝑇𝑇1 = 470𝐾𝐾,
𝑚𝑚
𝑉𝑉1 = 180 , 𝐴𝐴1 = 0.05𝑚𝑚2 , 𝑀𝑀𝑀𝑀1 = 0.414
𝑆𝑆
Solution

1 (1 + 0.2Ma )
2 3
A1 A1
=
1
1.547 =
A* = 0.0323m 2
A* Ma1 1.728 1.547

Calcul du débit massique au point 1: méthode 1 =m ρ=

1 AV
1 1 33.4kg / s
*
Calcul du débit massique au point 1: méthode 2
p A
= = 33.4kg / s
m 0.6847 01

RT01
Exemple 9.4 (suite III)

Solution

Pour déterminer le nombre de Mach au point 2, on doit résoudre cette

équation non linéaire. Pour ce faire, on peut appliquer une méthode numérique
ou bien utiliser des tables

A2 0.036
= = =
1.115
1 (1 + 0.2Ma )2
2 3

A* 0.0323 Ma2 1.728

La réponse exacte en p01

régime subsonique est: Ma2 = 0.6758 =p2 = 415kPa
(1 + 0.2Ma2 )
2 3.5
Exemple 9.4 (suite IV)

Solution au point 2E(subsonique): avec la table

Subsonique
M A/A* p/p 0 ρ/ρ 0 Τ/Τ 0 V/V*
0 Inf 1.0000 1.0000 1.0000 0
0.0500 11.5914 0.9983 0.9988 0.9995 0.0548
0.1000 5.8218 0.9930 0.9950 0.9980 0.1094
0.1500 3.9103 0.9844 0.9888 0.9955 0.1639
0.2000 2.9635 0.9725 0.9803 0.9921 0.2182
0.7 − 0.65
Ma2 =0.65 + (1.0944 − 1.115) =0.675 0.2500 2.4027
0.3000 2.0351
0.9575
0.9395
0.9694
0.9564
0.9877
0.9823
0.2722
0.3257
1.0944 − 1.1356 0.3500 1.7780 0.9188 0.9413 0.9761 0.3788
0.4000 1.5901 0.8956 0.9243 0.9690 0.4313
0.4500 1.4487 0.8703 0.9055 0.9611 0.4833
0.5000 1.3398 0.8430 0.8852 0.9524 0.5345
0.5500 1.2549 0.8142 0.8634 0.9430 0.5851
0.6000 1.1882 0.7840 0.8405 0.9328 0.6348
A/A*=1.115 0.6500 1.1356 0.7528 0.8164 0.9221 0.6837
0.7000 1.0944 0.7209 0.7916 0.9107 0.7318
0.7500 1.0624 0.6886 0.7660 0.8989 0.7789
0.8000 1.0382 0.6560 0.7400 0.8865 0.8251
0.8500 1.0207 0.6235 0.7136 0.8737 0.8704
0.9000 1.0089 0.5913 0.6870 0.8606 0.9146
0.9500 1.0021 0.5595 0.6604 0.8471 0.9578
Exemple 9.4 (suite V)

Solution au point 2F(supersonique): avec la table

Supersonique

1 (1 + 0.2Ma2 )
3 2 M A/A* p/p 0 ρ/ρ 0 Τ/Τ 0 V/V*
A2
= 1.115 1.0000 1.0000 0.5283 0.6339 0.8333 1.0000
A* Ma2 1.728 1.0500 1.0020 0.4979 0.6077 0.8193 1.0411
1.1000 1.0079 0.4684 0.5817 0.8052 1.0812
1.1500 1.0175 0.4398 0.5562 0.7908 1.1203
1.2000 1.0304 0.4124 0.5311 0.7764 1.1583
1.2500 1.0468 0.3861 0.5067 0.7619 1.1952
1.3000 1.0663 0.3609 0.4829 0.7474 1.2311
1.3500 1.0890 0.3370 0.4598 0.7329 1.2660
1.4000 1.1149 0.3142 0.4374 0.7184 1.2999
A/A*=1.115
Ma2 = 1.4 1.4500 1.1440 0.2927 0.4158 0.7040 1.3327
1.5000 1.1762 0.2724 0.3950 0.6897 1.3646
1.5500 1.2116 0.2533 0.3750 0.6754 1.3955
1.6000 1.2502 0.2353 0.3557 0.6614 1.4254
p01 1.6500 1.2922 0.2184 0.3373 0.6475 1.4544
p2 = 177 kPa 1.7000 1.3376 0.2026 0.3197 0.6337 1.4825
(1 + 0.2Ma2 )
2 3.5 1.7500 1.3865 0.1878 0.3029 0.6202 1.5097
1.8000 1.4390 0.1740 0.2868 0.6068 1.5360
1.8500 1.4952 0.1612 0.2715 0.5936 1.5614
1.9000 1.5553 0.1492 0.2570 0.5807 1.5861
1.9500 1.6193 0.1381 0.2432 0.5680 1.6099
2.0000 1.6875 0.1278 0.2300 0.5556 1.6330
 =R 287 m2 / ( s 2 K=
) , k 1.4,
= c p 1005 m2 / ( s 2 K )
Exemple 9.5

On désire accélérer de l’air (k=1.4) du repos à 𝒑𝒑𝟎𝟎 = 𝟐𝟐𝟐𝟐𝟐𝟐 𝒌𝒌𝒌𝒌𝒌𝒌, 𝑻𝑻𝟎𝟎 = 𝟓𝟓𝟓𝟓𝟓𝟓𝟓𝟓, au travers
un col, jusqu’à Ma =2.5 à la sortie(e), Le débit massique est 𝒎𝒎̇ = 𝟑𝟑𝟑𝟑𝟑𝟑/𝒔𝒔. Considérez
un écoulent isentropique et estimez 𝑨𝑨𝒄𝒄𝒄𝒄𝒄𝒄 , 𝒑𝒑𝒆𝒆 , 𝑻𝑻𝒆𝒆 , 𝑽𝑽𝒆𝒆 , 𝑨𝑨𝒆𝒆

Solution

Sortie supersonique => col sonique => Acol=A*

p0 A*
= 3kg / s
m 0.6847 = A=
Acol *
0.00830m 2
RT0

( )
3
p0 A 1 1 + 0.2 Ma 2
pe = =
11.7 kPa e
= 2.64
( )
3.5
1 + 0.2 Ma 2
A*
Ma 1.728
= =
Ae 2.64 A* 0.0219m 2
Version adaptée du livre: Fluid Mechanics, par Frank M. White, 7e édition
 =R 287 m2 / ( s 2 K=
) , k 1.4,
= c p 1005 m2 / ( s 2 K )
Exemple 9.5 (suite)

Solution

Sortie supersonique => col sonique => Acol=A*

T0 = = =
Ve a=
Te = 222 K ae kRTe 299m / s e Ma 747 m / s
1 + 0.2 Ma 2
Dans la seconde partie
Nous allons regarder le choc
normal et les écoulements
dans les tuyères
convergentes-divergentes