Chapitre II Ecoulement Isentropique Unidimensionnel En Conduit A Section Variable

R.Lanani-Benchabi

CHAPITRE 2
Ecoulement isentropique unidimensionnel dans un conduit à section variable

II-1 Introduction

Les écoulements isentropiques unidimensionnels compressibles de fluides non visqueux se

rencontrent quand on envisage des écoulements dans des conduits ou des écoulements
externes. On peut citer notamment le cas de la tuyère convergente- divergente sans onde de
choc, ni transfert de chaleur ou frottements ou l’écoulement autour d’un profil d’aile d’avion.

II-2 Etat générateur

Cette notion d’état générateur est particulière aux fluides compressibles.

L’état générateur est l’état du fluide lorsque sa vitesse est nulle,
par exemple les conditions régnant dans le réservoir alimentant une tuyère.
Aux grandeurs caractérisant cet état nous affectons l’indice 0 : P0, 0, T0….

II-3 Propriétés isentropiques et de stagnation d’un écoulement

II-3-1 Température totale ou de stagnation

Cette température se produit en un point où la vitesse d’écoulement est nulle et si en ce point
l’écoulement est adiabatique on peut écrire :

𝑼𝟐
CP .T0 = CP .T + (2-1)
𝟐

La température d’arrêt est :

𝑼𝟐
T0 =T + 𝟐𝑪 (2-2)
𝒑

𝑈
D’après la définition du nombre de Mach, on peut écrire : Ma = 𝑎

𝑻𝟎 𝜸−𝟏
=1+ 𝐌𝐚𝟐 (2-3)
𝑻 𝟐

II-3-2 Pression totale ou de stagnation

La pression d’arrêt est la pression qui est atteinte de manière isentropique. En utilisant les
relations d’isentropiques on trouve :

𝑷𝟎 𝜸−𝟏
= (1 + 𝐌𝐚𝟐 ) /-1 (2-4)
𝑷 𝟐

II-3-3 densité totale ou de stagnation

1
Chapitre II Ecoulement Isentropique Unidimensionnel En Conduit A Section Variable
R.Lanani-Benchabi

𝟎 𝜸−𝟏
= (1 + 𝐌𝐚𝟐 ) 1/-1 (2-5)
 𝟐
II-4 Ecoulement isentropique dans des canalisations

II-4-1 Equation de Barré Saint Venant

Pour un fluide compressible, sans transfert de chaleur, et en introduisant le concept de
l’enthalpie on a :
𝑼𝟐
h+ = cte (2-6)
𝟐

En appliquant cette équation entre l’état générateur et un état quelconque du fluide on a :

𝑼𝟐
Cp T0 = Cp T + (2-7)
𝟐

De cette équation on tire la vitesse d’écoulement et on obtient :

𝑻
𝐔𝟐 = 2 Cp (T0 – T) = 2 Cp T0 (1- 𝑻 ) (2-8)
𝟎

𝑇 𝑃0 -1/
En introduisant les relations isentropiques ( 𝑇0 = ( ) ) et la relation du gaz parfait (P0 =
𝑃
0.r.T0) On obtient :

𝜸 𝑷 𝑷
𝐔𝟐 = 2 𝜸−𝟏  𝟎 (1- (𝑷 )-1/ ) (2-9)
𝟎 𝟎

Cette équation permet le calcul de la vitesse d’écoulement en un point connaissant la pression

statique en ce point et la pression d’arrêt.

II-4-2 Vitesse maximale (vitesse limite)

Il est à noter qu’il existe une vitesse maximale pour l’écoulement. Si on effectue une détente
du gaz depuis un réservoir jusqu’au vide (P= 0), ce qui entraine T = 0, On obtient d’après
l’équation de Barrée Saint Venant :

𝟐𝜸 𝑷𝟎
Umax = = 𝟐𝐂𝐩𝐓𝟎 (2-10)
𝜸−𝟏 𝝆𝟎

Remarque
Cette vitesse maximale ne dépend que des conditions régnant dans l’état générateur.Elle ne
peut être atteinte que pour une détente isentropique.
Cette vitesse est intéressante , elle est la vitesse maximale que peuvent atteindre les gaz
d’échappement d’une fusée dans le vide, elle conditionne donc la poussée maximale des
moteurs de fusée.

2
Chapitre II Ecoulement Isentropique Unidimensionnel En Conduit A Section Variable
R.Lanani-Benchabi

II-5 Conditions soniques

Le point de l’écoulement où Ma = 1 est appelé point critique. Les caractéristiques du fluide

en ce point sont appelées grandeurs critiques et sont notées a*, *, P*, T*, A*.On peut les
calculer à partir des conditions d’arrêt.

𝑻𝟎 𝜸+𝟏
= (2-11)
𝑻∗ 𝟐

𝑷𝟎 𝜸+𝟏 /-1
=( ) (2-12)
𝑷∗ 𝟐

𝟎 𝜸+𝟏 1/-1
=( ) (2-13)
∗ 𝟐
Pour l’air  =1,4 :
𝑻∗
= 0,8333 (2-14)
𝑻𝟎

𝑷∗
= 0 ,5283 (2-15)
𝑷𝟎

∗
= 0,6339 (2-16)
𝟎

II-6 Ecoulement isentropique avec changement de section Théorème d’Hugoniot

Soit l’écoulement unidimensionnel stationnaire d’un fluide parfait compressible dans une
conduite à section variable.

II-6-1 Equations de base de l’écoulement isentropique

En combinant les relations isentropiques d’un écoulement compressible avec l’équation de
continuité, on peut étudier des problèmes pratiques des écoulements compressibles.

1- Equation de continuité
Pour l’écoulement unidimensionnel stationnaire, l’équation de continuité exprime la
conservation du débit :
.
m = .S.U = cste (2-17)

ou par différentiation on obtient :

𝐝𝛒 𝐝𝐒 𝐝𝐔
+ + =0 (2-18)
𝛒 𝐒 𝐔

𝐝𝛒 𝐝𝐒 𝐝𝐔
=- - (2-19)
𝛒 𝐒 𝐔

3
Chapitre II Ecoulement Isentropique Unidimensionnel En Conduit A Section Variable
R.Lanani-Benchabi

2- Equation de quantité de mouvement exprimé par l’équation d’Euler

L’équation d’Euler décrivant le mouvement stationnaire est :

𝐝𝐏
UdU + =0 (2-20)
𝛒

𝟏 𝐝𝐔
= −𝑼 𝐝𝐏 (2-21)
𝛒

En combinant les équations (2-19) et (2-21) on obtient :

𝐝𝐔 𝐝𝐒 𝐝𝐔
−𝑼 .d = - - (2-22)
𝐝𝐏 𝐒 𝐔

𝑈 𝜕𝑃
Puisque Ma = 𝑎 et la vitesse du son est a2 = 𝜕𝜌 , ce qui donne :

𝐝𝐔 𝟏 𝐝𝐒
=( ) (2-23)
𝐔 𝑴𝒂𝟐 −𝟏 𝐒

C’est la première équation de la théorie d’Hugoniot donnant une relation entre la variation de la
section d’une conduite et de la vitesse.

On a aussi d’après l’équation d’Euler (2-20) :

𝐝𝐏 𝟏 𝐝𝐒
=( ) (2-24)
𝐔 𝟐 𝟏−𝑴𝒂𝟐 𝐒

En combinant Les équations (2-19) et (2-23) on trouvera :

𝐝𝛒 𝐌𝐚𝟐 𝐝𝐒
= (𝟏−𝐌𝐚𝟐) (2-25)
𝛒 𝐒

II-6-2 Effet de la variation de la section

Les équations d’Hugoniot sont très importantes car elles permettent de voir l’effet de la
variation de la section sur les variations des grandeurs physiques en fonction du régime
d’écoulement.
- Ecoulement subsonique (Ma < 1):
Lorsque l’écoulement est subsonique (𝐌𝐚𝟐 − 𝟏) < 0, une augmentation de la section entraine
une diminution de la vitesse et inversement. Ainsi pour un écoulement subsonique, pour
augmenter la vitesse il faut un convergent et pour la diminuer, le canal doit être un divergent.
- Ecoulement supersonique (Ma > 1):
Lorsque l’écoulement est supersonique (𝐌𝐚𝟐 − 𝟏) > 0, une augmentation de la vitesse est
associée à une augmentation de la section et inversement. Ainsi pour un écoulement
supersonique, pour augmenter la vitesse il faut un divergent et pour la diminuer, le canal doit
être un convergent.

4
Chapitre II Ecoulement Isentropique Unidimensionnel En Conduit A Section Variable
R.Lanani-Benchabi

Ma  1 U augmente Ma 1 U augmente

P diminue Supersonique P diminue
Subsonique S diminue S augmente

Tuyére convergente Tuyére divergente

U diminue U diminue
P augmente P augmente
S augmente S diminue
Ma  1 Ma 1
subsonique supersonique

Diffuseur

- Ecoulement sonique Ma = 1 :
Puisque l’accélération infinie est physiquement impossible, donc dS = 0, ce qui signifie
mathématiquement une section minimale ou maximale, mais la section maximale n’a pas de
sens physique et en plus si on tient compte de la discussion précédente, l’écoulement sonique
ne peut avoir lieu en section minimale de la conduite appelée col.

Figure 2-1 Ecoulement compressible dans une conduite à section variable

Remarque :
Si l’écoulement est sonique au col, alors l’évolution en aval peut être soit subsonique, soit
supersonique, c’est la condition aux limites en aval du divergent qui détermine lequel des
deux comportements se réalise.

II-7 Débit massique d’une conduite de section variable

Le point essentiel dans l’étude des écoulements compressibles dans les conduites à section
variable est la détermination du débit en fonction du nombre de Mach, de la section et des
conditions d’arrêt.

5
Chapitre II Ecoulement Isentropique Unidimensionnel En Conduit A Section Variable
R.Lanani-Benchabi

On a :
.
m = .S.U = cste (2-26)

P
En substituant  = rT et U =a. Ma, l’équation du débit devient :
.
𝐏 𝑷
m = 𝐫𝐓 .S. a. Ma = Ma. 𝒓𝑻 . 𝑺. 𝛄. 𝐫. 𝐓 (2-27)
.
𝛄
m = Ma . S.P. (𝒓𝑻) 1/2 (2-28)

En introduisant les relations isentropiques, on aura :

.
𝛄 𝜸−𝟏
m = Ma . S.P0. (𝒓𝑻 ) ½.(1 + 𝐌𝐚𝟐 ) +1/2(1-) (2-29)
𝟎 𝟐

II-8 Equation du rapport de section

En combinant l’équation du débit massique avec les relations de l’écoulement isentropiques,

on peut tirer une expression reliant seulement la section S et le nombre de Mach.
Puisque le débit est conservé sous les conditions soniques :
.
m = .S.U = *.S*.U* (2-30)

où :
𝐒 𝛒∗ 𝐔∗
= . (2-31)
𝐒∗ 𝛒 𝐔

D’après les relations trouvées précédemment pour l’écoulement isentropique on a :

𝛒∗ 𝛒∗ 𝛒𝟎 𝟐 𝜸−𝟏
== . = [(𝜸+𝟏)(1 + 𝐌𝐚𝟐 )] 1/-1 (2-32)
𝛒 𝛒𝟎 𝛒∗ 𝟐

Aussi en considérant les relations isentropiques pour les températures on aura :

𝐔∗ 𝐚∗ 𝐌𝐚∗ 𝟏 𝟐 𝜸−𝟏
= 𝐚. = 𝐌𝐚 . [( 𝜸+𝟏)(1 + 𝐌𝐚𝟐 )] 1/2 (2-33)
𝐔 𝐌𝐚 𝟐
En combinant les équations (2-32) (2-33) on aboutit à

𝐒 𝟏 𝟐 𝜸−𝟏
= 𝐌𝐚 . [( 𝜸+𝟏)(1 + 𝐌𝐚𝟐 )] +1/2(-1) (2-34)
𝐒∗ 𝟐

6
Chapitre II Ecoulement Isentropique Unidimensionnel En Conduit A Section Variable
R.Lanani-Benchabi

T.D N°1 Ecoulement isentropique

Exercice1 :
La figure ci contre représente une chaudière qui produit de la vapeur d’eau à un débit massique de
13,4Kg/s.
Par une canalisation cylindrique, la vapeur arrive dans une
section S de diamètre 10cm à une pression de 15bar et une
température de 541K.
On donne les caractéristiques de la vapeur d’eau :
𝛾 = 1,3 ; r = 462J/Kg.°K.
1- On suppose que la vapeur est un gaz parfait, calculer la
masse volumique de la vapeur d’eau en sortie de la chaudière.
2- Déterminer la vitesse d’écoulement.
3- Calculer la vitesse du son.
4- En déduire le nombre de Mach et préciser la nature de l’écoulement.
5- Calculer la pression d’arrêt à l’intérieur de la chaudière.

Solution

1- La masse volumique de la vapeur

La vapeur est un gaz parfait :

𝑃 15.10 5
P = 𝜌. 𝑟. 𝑇 ⤇ 𝜌 = = = 6Kg/m3
𝑟.𝑇 462.541

2- La vitesse d’écoulement
. . . .
m m 4.m 4.13,4
m = 𝜌. 𝑠. u ⤇ u = = 𝜋 .𝑑 2
= = = 284,5m/s
.𝑠 . .𝜋 .𝑑 2 6.3,14.0,12
4

3- La vitesse du son

c = γ. r. T = 1,3.462.541 = 570,021m/s

4- Le nombre de Mach

u 284,5
Ma = = = 0,449≈ 0,5
c 570,021

Ma < 1 ⤇ l’écoulement est subsonique

5- La pression d’arrêt

1,3
𝛾 𝛾
𝑃0 𝛾−1 𝛾 −1 1,3−1 2 1,3−1
=[1 + 𝑀𝑎2 ]𝛾 −1 ⤇ P0 = P. [1 + 𝑀𝑎2 ]𝛾 −1 5
= 15.10 . 1 + 0,5
𝑃 2 2 2

7
Chapitre II Ecoulement Isentropique Unidimensionnel En Conduit A Section Variable
R.Lanani-Benchabi

P0 = 16,041.105 Pa = 16,041 bar

Exercice2 :
Un réservoir contient de l’air comprimé à une pression de 4 bar, supposé pression d’arrêt à l’état
initial. L’ouverture d’une vanne dans ce réservoir provoque la détente de l’air vers l’extérieur sous
forme d’un jet ayant un diamètre de 5mm. Les paramètres extérieurs du jet d’air à l’état final sont :
P = 1 bar ; T= 25C.
1- Calculer la vitesse du son à l’extérieur du réservoir.
2- Déterminer la masse volumique de l’air à l’extérieur du réservoir.
3- Ecrire l’équation de Saint – Venant, en terme de rapport de pression, entre un point d’arrêt et un
point sur le jet d’air.
4- En déduire le nombre de Mach au niveau du jet d’air. Quel est la nature de l’écoulement.
5- Calculer la vitesse d’écoulement du jet d’air.

Solution

1- La vitesse du son à l’extérieur du réservoir

c = γ. r. T = 1,4.287. (25 + 273) = 346m/s

2- La masse volumique

𝑃 10 5
P = 𝜌. 𝑟. 𝑇 ⤇ 𝜌 = = = 1,169Kg/m3
𝑟.𝑇 287.298

3- Equation de Barré saint-Venant

𝛾
𝑃0 𝛾−1
=[1 + 𝑀𝑎2 ]𝛾 −1
𝑃 2

4- Le nombre de Mach

𝛾 −1
2 𝑃 𝛾
Ma = .[ 0 − 1] = 1,558≈ 1,56
𝛾−1 𝑃

Ma  1 ⤇ L’écoulement est supersonique

5- La vitesse d’écoulement

u
Ma = ⤇ u = Ma .c = 1,56.346 = 539,76m/s
c

Exercice3 :
De l’hélium entre dans une tuyère convergente divergente à des conditions spécifiques.
Déterminer la température et la pression au col sachant que la température, la pression et la vitesse à
l’entrée sont respectivement T = 800K ; P = 0, 7MPa et U = 100m/s.
On donne : 𝛾 = 1,667 ; Cp = 5,1926KJ/Kg.°K.

8
Chapitre II Ecoulement Isentropique Unidimensionnel En Conduit A Section Variable
R.Lanani-Benchabi

Solution

- La température d’arrêt :

u 21 100 2
T0 = T1 + = 800 + = 800,86 𝐾 ≈ 801K
2.C p 2.5192,6

- La pression d’arrêt :

𝛾 𝛾 1,667
𝑃0 T T 801
=[T 0 ]𝛾 −1 ⤇ P0 = P1. [T 0 ]𝛾 −1 = 0,7.106 [800 ]1,667 −1 = 0,702MPa
𝑃1 1 1

- La température au col :

𝑇0 𝛾+1 2 2
On a :
𝑇∗
=
2
⤇ T* = T0.( 𝛾 +1) = 801.( 1,667+1) = 600,67K≈ 601K

𝛾 𝛾 𝛾 1,667
𝑃0 T 𝛾+1 2 2
=[ T0∗ ]𝛾 −1 = [ ]𝛾 −1 ⤇ P* = P0. [𝛾+1 ]𝛾 −1 = 0,702.106. [1,667+1 ]1,667 −1
𝑃∗ 2

P* = 0,342MPa

Exercice4 :
Du gaz carbonique est détendu isentropiquement dans un canal de la pression P 1= 125KPa et la
température T1= 100°C à la pression P2= 80KPa et une vitesse U2 = 325m/s.
Calculer :
1- La température T2 et le nombre de Mach à l’état final.
2- La température de stagnation.
3- La pression de stagnation.
4- La vitesse U1 et le nombre de Mach à l’entrée du canal.
On donne : 𝛾 = 1,3 ; r = 189J/Kg.°K.

Solution

1- La température de l’état 2

On a :

γ −1 γ −1 1,3−1
T2 P2 γ P2 γ 80 1,3
= ⤇ T2 = T1. = (100+273). = 336K
T1 P1 P1 125

u2 u2 325
Ma2 = = = = 1,13
c2 γ.r.T2 1,3.189.336

L’écoulement est transonique

2- La température de stagnation

9
Chapitre II Ecoulement Isentropique Unidimensionnel En Conduit A Section Variable
R.Lanani-Benchabi

𝑇0 𝛾 −1 1,3−1
=1+ 𝑀𝑎2 ⤇ T0 = 336. [ 1 + . (1,13)2] = 401K
𝑇2 2 2

3- La pression de stagnation

𝛾 γ 𝛾 1,3
𝑃0 T γ−1 T 3 401
= [ 0 ]𝛾 −1 = [1 + Ma2 ]γ −1 ⤇ P0 = P2. [ 0 ]𝛾 −1 = 80.10 . [ ]1,3−1 = 171KPa
𝑃2 T2 2 T2 336

4- La vitesse et le nombre de Mach à l’entrée du canal

u 21
T0= T1 + ⤇ u1 = 2. 𝐶𝑝 (𝑇0 − 𝑇1 ) = 2.819(401 − 373) = 214m/s
2.C p

𝛾.𝑟 1,3.189
Cp = = = 819J/Kg.°K
𝛾−1 1,3−1

u1 u1 214
Ma1 = = = = 0,71
c1 γ.r.T1 1,3.189.373

Ecoulement subsonique à l’entrée du canal.

Exercice 5 :
De l’oxygène s’écoule dans un canal de 5cm de diamètre. A une section donnée la pression P est égale
à 120KPa et le débit d’écoulement vaut 0,4kg/s.
Calculer à cette section :
1- La vitesse et le nombre de Mach.
2- La masse volumique de stagnation sachant que T0 = 300°C.
On donne : 𝛾 = 1,4 ; r = 260J/Kg°K.

Solution

𝛾.𝑟 1,4.260
Cp = = = 910J/Kg.°K
𝛾−1 1,4−1

1- La vitesse et le nombre de Mach

On a :

u2
T0= T + ⤇ u1 = 2. 𝐶𝑝 (𝑇0 − 𝑇)
2.C p

Et on a aussi

. . . . .
P m.r.T m.r.T 4 .m.r.T 4 .m .260.T
m = 𝜌. 𝑠. u = . s. u ⤇ u= = 𝜋 .𝑑 2
= ==
r.T 𝑃.𝑠 𝑃. 4 𝑃.𝜋.𝑑 2 120.10 3 .𝜋.(0,05)2

2 équations à 2 inconnus, on obtient :

10
Chapitre II Ecoulement Isentropique Unidimensionnel En Conduit A Section Variable
R.Lanani-Benchabi

T = 542K et u = 239m/s

u 239 239
Ma = = = = 0,538
c γ.r.T 1,4.260.542

2- La masse volumique de stagnation

1 1 1
0 γ−1 γ−1 1,4−1
⤇ 0 = 0,852. [1 +
1,4−1
= [1 + Ma2 ]γ −1 Ma2 ]γ −1 = 0,852. 1 + (0,538)2
 2 2 2

0 = 0,98Kg/m3

𝑃 120.10 3
P = 𝜌. 𝑟. 𝑇 ⤇ 𝜌 = = = 0,852Kg/m3
𝑟.𝑇 260.542

Exercice 6 :
Dans la veine d’essai d’une soufflerie supersonique, les conditions de l’air sont :
P= 105 Pa, T= -23°C, et Ma= 3
Calculer la température et la pression du fluide dans le réservoir.
Si l’aire de la section d’essai de test est 10 -2 m2.
Trouver le débit massique dans la soufflerie, la vitesse et la section au col.

Solution

1- La température et la pression dans le réservoir

𝑇0 𝛾 −1 1,4−1 1,4−1
=1+ 𝑀𝑎2 = 1 + . (3)2 ⤇ T0 = T.[ 1 + . (3)2 ] = 700K
𝑇 2 2 2

𝛾
𝛾 𝛾 1,4
𝑃0 T 𝛾−1 T 700
=[ T0 ]𝛾 −1 = [1 + 2
𝑀2𝑎 ]𝛾−1 ⤇ P0 = P. [ 0 ]𝛾 −1 = 105. [ ]1,4−1 = 280KPa
𝑃 T 250

2- Le débit massique dans la soufflerie

.
P P 𝛾 1,4
m = 𝜌. 𝑠. u = . s.Ma.c = . s.Ma. γ. r. T = P.s.Ma = 105.10-2.3. = 13,25 Kg/s
r.T r.T 𝑟.𝑇 287.250

La vitesse au col

u∗ u∗
Ma = 1 = = ⤇ u* = γ. r. T ∗
c∗ ∗
γ.r.T

𝑇0 𝛾+1 2 2
𝑇∗
=
2
⤇ T* = T0.( 𝛾 +1) = 700. ( 1,4+1) = 583,3K

u* = 1,4.287.583,3 = 484m/s

11
Chapitre II Ecoulement Isentropique Unidimensionnel En Conduit A Section Variable
R.Lanani-Benchabi

La section au col
. . .
𝑚 𝑚 .𝑟 .𝑇 ∗
m = ρ∗ .s ∗ . u∗ ⤇ s ∗ = =
ρ ∗ .u ∗ P ∗ .u ∗

On a :

𝛾 𝛾 1,4
𝑃0 𝛾+1 2 2
=[ ]𝛾 −1 ⤇ P* = P0. [𝛾+1 ]𝛾 −1 = 280.103. [1,4+1 ]1,4−1 = 147,92KPa
𝑃∗ 2

D’où :
. .
𝑚 𝑚 .𝑟.𝑇 ∗ 13,25.287.583,3
s∗ = = = = 2360mm2
ρ ∗ .u ∗ P ∗ .u ∗ 147,92.484

Exercice7 :
De l’air à 403°K et 1 atm pénètre dans une tuyère convergente à une vitesse de 150 m/s et subit une
détente isentropique jusqu’à une pression de sortie de 76 KPa.
Si la section d’entrée de la tuyère est 5.10-3m2
Trouver :
1- La température, la pression et l’enthalpie de l’état générateur.
2- Le nombre de Mach à l’entrée.
3- La température, le nombre de Mach et la section de sortie.

Solution

1- La température de l’état générateur

u 21 150 2
T0= T1 + = 403 + = 414,25K
2.C p 2.1000

La pression de l’état générateur

𝛾 𝛾 1,4
𝑃0 T T 5 414,25
= [ 0 ]𝛾 −1 ⤇ P0 = P1. [ 0 ]𝛾 −1 = 1,013.10 . [ ]1,4−1 = 111560Pa
𝑃1 T1 T1 403

L’enthalpie de l’état générateur

h0 = Cp.T0 = 103. 414,25 = 414250J/Kg

2- Le nombre de Mach

u1 u1 150
Ma1 = = = = 0,372
c1 γ.r.T1 1,4.287.403

3- La température de sortie

12
Chapitre II Ecoulement Isentropique Unidimensionnel En Conduit A Section Variable
R.Lanani-Benchabi

γ −1 γ −1 1,4−1
T0 P0 γ P2 γ 76000 γ 1,4
= ⤇ T2 = T0.. = 414,25 = 371,18K
T2 P2 P0 111560

u 22
T0= T2 + ⤇ u2 = 2. 𝐶𝑝 (𝑇0 − 𝑇2 ) = 2.1000(414,25 − 371,18) = 293m/s
2.C p
Le nombre de Mach

u2 u2 293
Ma2 = = = = 0,758
c2 γ.r.T2 1,4.287.371,18

La section de sortie
. . .
m m 0,687
m = ρ1 . s1 . u1 = ρ2 . s2 . u2 ⤇ s2 = = P2 = 76000 = 31,38.10-4m2
ρ 2 .u 2 .u .293
r.T 2 2 287 .371 ,18

.
P1 1,013.10 5
m = ρ1 . s1 . u1 = .s .u = . 5. 10−3 . 150= 0,687Kg/s
r.T 1 1 1 287.403

Exercice8:
L’air provenant d’un réservoir s’écoule dans un conduit adiabatique, la température et la pression
sont : T= 400°K et P = 500KPa.
Déterminer à une section du conduit où P= 430 KPa :
1- Le nombre de Mach.
2- La Température.
3- La masse volumique.
4- La vitesse d’écoulement.

Solution

1-Le nombre de Mach

γ−1
𝛾
𝑃0 𝛾−1 𝟐 P γ
=[1 + 𝑀𝑎12 ]𝛾 −1 ⤇ Ma1 = .[ 0 − 1] = 0,447
𝑃1 2 𝜸−𝟏 P2

2- La température

𝑇0 𝛾−1 𝑇0 400
=1 + 𝑀𝑎12 ⤇ T1 = 𝛾 −1 = 1,4−1 = 384,62K
𝑇1 2 1+ 2 𝑀𝑎 12 1+ 2 .(0,447)2

3- La masse volumique

𝑃1 430.10 3
P1 = ρ1 . r. T1 ⤇ ρ1 = = = 3,89Kg/m3
𝑟.𝑇1 287.384,62

4- La vitesse d’écoulement
u
Ma = ⤇ u = Ma.c = Ma. γ. r. T = 0,447. 1,4.287.384,62 = 175,72m/s
c

13
Chapitre II Ecoulement Isentropique Unidimensionnel En Conduit A Section Variable
R.Lanani-Benchabi

T.D N°2 Ecoulement isentropique

Exercice1 :
On veut décomprimer de l’air à P 0 = 200KPa, T0 = 500K en utilisant une conduite (convergente-
divergente) qui produit un nombre de Mach 2,5 à la sortie et donne un débit de 3kg/s. Déterminer :
1- L’expression du débit maximum au col.
2- La section au col.
3- La section de sortie.
4- La pression de sortie.
5-La vitesse de sortie.

Solution

1- Expression du débit maximum

.
m = ρ∗ . s ∗ . u∗

on a

1 1
0 γ+1 γ−1 2 γ−1
= [ ] ⤇ 𝜌 ∗ = 0 . [ ]
𝜌∗ 2 γ+1

T∗
et u* = c* = γ. r. T ∗ = γ. r. .T
T0 0

𝑇0 𝛾+1 2
𝑇∗
=
2
⤇ T* = T0.( 𝛾 +1)

Donc le débit d’étranglement est égal à :

.
1
2 γ−1 2
m = 0 . [ ] . 𝑠∗. γ. r. .T
γ+1 𝛾 +1 0

P0
pour l’air, on a P0 = 0 .r.T0 ⤇ 0 =
𝑟.T0
.
P0 .𝑠 ∗ .0,6847
m= 1
(𝑟.𝑇0 )2

2- La section au col
. .
1 1
P0 .𝑠 ∗ .0,6847 m .(𝑟.𝑇0 )2 3 .(287.500 )2
m= 1 ⤇ 𝑠∗ = = 0,0083m2
)2
0,6847 .P0 0,6847 .200.10 3
(𝑟.𝑇0

3- La section de sortie

14
Chapitre II Ecoulement Isentropique Unidimensionnel En Conduit A Section Variable
R.Lanani-Benchabi

𝐒𝟐 𝟏 𝟐 𝜸−𝟏
= 𝐌𝐚 . [( 𝜸+𝟏)(1 + 𝟐
𝐌𝐚𝟐 )] +1/2(-1)
𝐒∗ 𝟐 𝟐
⤇ S2 = 218,8cm2

4- La pression de sortie

𝛾
𝑃0 𝛾−1 𝑃0 200.10 3
=[1 + 𝑀𝑎22 ]𝛾 −1 = 17,08 ⤇ 𝑃2 = = = 11,7KPa
𝑃2 2 17,08 17,08

𝑇0 𝛾−1
=1+ 𝑀𝑎22 ⤇ 𝑇2 = 222,2K
𝑇2 2

5- La vitesse de sortie

u2
Ma2 = ⤇ u2 = Ma2.c2 = Ma2. γ. r. T2 = 2,5. 1,4.287.222,2 = 747m/s
c2

Exercice2 :
De l’air à la pression de 1MPa et la température de 800K entre dans un convergent-divergent.
L’écoulement est stable, unidimensionnel et isentropique avec 𝛾 = 1,4.
Pour un nombre de Mach à la sortie Ma2 = 2 et une section au col A* =20cm2.
Déterminer en utilisant la table thermodynamique de l’écoulement isentropique :
1- Les grandeurs physiques critiques.
2- Les grandeurs physiques à la sortie.
3- Le débit massique.

Solution

1- Les grandeurs physiques critiques

𝑃0 10 6
P0 = ρ0 . r. T0 ⤇ ρ0 = = = 4,355Kg/m3
𝑟.𝑇0 287.800

D’après la table de l’écoulement isentropique pour Ma =1, on a:

La pression critique

𝑃∗
= 0,5283 ⤇ 𝑃∗ = 106. 0,5283 = 0,5283MPa
𝑃0

La température critique

𝑇∗
= 0,8333 ⤇ 𝑇 ∗ = 800. 0,8333= 666,6K
𝑇0

La masse volumique critique

𝜌∗
= 0,6339⤇ 𝜌 ∗ = 4,355. 0,6339= 2,761Kg/m3
𝜌0

15
Chapitre II Ecoulement Isentropique Unidimensionnel En Conduit A Section Variable
R.Lanani-Benchabi

u∗ u∗
Ma = 1 = = ⤇ u* = γ. r. T ∗ = 1,4.287.666,6 = 517,5m/s
c∗
γ.r.T∗

2- Les grandeurs physiques à la sortie

d’après la table de l’écoulement isentropique, on a :

𝑃2
= 0,1278 ⤇ 𝑃2 = 0,1278MPa
𝑃0

𝑇2
= 0,5556⤇ 𝑇2 = 444,5K
𝑇0

𝜌2
= 0,2300⤇ 𝜌2 = 1,002Kg/m3
𝜌0

u2 u2
Ma2 = = ⤇ u2 = Ma2 . γ. r. T2 =2. 1,4.287.444,5 = 845,2m/s
c2 γ.r.T2

On a :

𝐴2
= 1,6875 ⤇ 𝐴2 = 𝐴∗ . 1,6875 = 20. 1,6875 = 33,75cm2
𝐴∗

3- Le débit massique
.
m = 𝜌 ∗ . 𝑠 ∗ . 𝑢∗ = 2,761.20.10-4.517,5 = 2,86Kg/s

Exercice3 :
Du Nitrogène entre dans une tuyère convergente avec une pression de 100KPa et une température de
400K. En supposant un écoulement isentropique et que le gaz est accéléré depuis M a = 0,3.
Calculer les grandeurs physiques du gaz à la sortie du convergent sachant que le rapport entre la
section de sortie et la section critique vaut 1,6281.

Solution

On a :

𝐴2
= 1,6281
𝐴∗

De la table thermodynamique, ce rapport correspond à :

Le nombre de Mach à la sortie

0,38 1,6587

Ma2 1,6281 Ma2 = 0,391

0,4 1,5901

16
Chapitre II Ecoulement Isentropique Unidimensionnel En Conduit A Section Variable
R.Lanani-Benchabi

Le rapport des températures pour Ma2 = 0,391

0,38 0,9719

T T
0.391 = 0,9701
T0 T0

0,4 0,9690

Le rapport des pressions pour Ma2 = 0,391

0,38 0,9052

P P
0.391 = 0,8993
P0 P0

0,4 0,8956

Les grandeurs physiques à la sortie sont :

𝑇2
T2 𝑇0 0,9701
= 𝑇1 = ⤇ T2 = 395K
T1 0,9823
𝑇0

𝑝2
p2 𝑝0 0,8993
= 𝑝1 = ⤇ p2 = 95,7KPa
p1 0,9395
𝑝0

17
Chapitre II Ecoulement Isentropique Unidimensionnel En Conduit A Section Variable
R.Lanani-Benchabi

18
Chapitre II Ecoulement Isentropique Unidimensionnel En Conduit A Section Variable
R.Lanani-Benchabi

19
Chapitre II Ecoulement Isentropique Unidimensionnel En Conduit A Section Variable
R.Lanani-Benchabi

20
Chapitre II Ecoulement Isentropique Unidimensionnel En Conduit A Section Variable
R.Lanani-Benchabi

21
Chapitre II Ecoulement Isentropique Unidimensionnel En Conduit A Section Variable
R.Lanani-Benchabi

22
Chapitre II Ecoulement Isentropique Unidimensionnel En Conduit A Section Variable
R.Lanani-Benchabi

23
Chapitre II Ecoulement Isentropique Unidimensionnel En Conduit A Section Variable
R.Lanani-Benchabi

24
Chapitre II Ecoulement Isentropique Unidimensionnel En Conduit A Section Variable
R.Lanani-Benchabi

25
Chapitre II Ecoulement Isentropique Unidimensionnel En Conduit A Section Variable
R.Lanani-Benchabi

26
Chapitre II Ecoulement Isentropique Unidimensionnel En Conduit A Section Variable
R.Lanani-Benchabi

27
Chapitre II Ecoulement Isentropique Unidimensionnel En Conduit A Section Variable
R.Lanani-Benchabi

28
Chapitre II Ecoulement Isentropique Unidimensionnel En Conduit A Section Variable
R.Lanani-Benchabi

29
Chapitre II Ecoulement Isentropique Unidimensionnel En Conduit A Section Variable
R.Lanani-Benchabi

30
Chapitre II Ecoulement Isentropique Unidimensionnel En Conduit A Section Variable
R.Lanani-Benchabi

31
Chapitre II Ecoulement Isentropique Unidimensionnel En Conduit A Section Variable
R.Lanani-Benchabi

32
Chapitre II Ecoulement Isentropique Unidimensionnel En Conduit A Section Variable
R.Lanani-Benchabi

33