Polycope Gazodynamique Alliche2021 PART1

REPUBLIQUE ALGERIENNE DEMOCRATIQUE ET POPULAIRE
MINISTERE DE L’ENSEIGNEMENT SUPERIEUR ET DE LA RECHERCHE SCIENTIFIQUE

UNIVERSITE Dr YAHYA Farés DE MEDEA
COURS DE GENIE MECANIQUE
Titre :
INTRODUCTION A LA
DYNAMIQUE DES GAZ
(Cours et Exercices)
Auteur : Dr ALLICHE Mounir

Maitre de Conférences « A »
Université de Médéa
N° de Série :
Année : 2021
Approuvé par le Conseil Scientifique de la Faculté en date du : 29/09/2021
P a g e | 25
Semestre : 2
Unité d’enseignement : UEF 1.2.1
Matière : Dynamique des gaz
VHS : 45h00 (cours : 1h30, TD : 1h30)
Crédits : 4
Coefficient : 2
Objectifs de l’enseignement :
La dynamique des gaz est un domaine très vaste qui a pour objectif théorique, l’étude des
écoulements compressibles à grandes vitesses. Ces types d’écoulements sont le plus souvent
rencontrés dans le domaine pratique de l’industrie aéronautique et spatiale. Le présent
module traite seulement l’approche monodimensionnelle des écoulements compressibles des
gaz parfait
Connaissances préalables recommandées :
Thermodynamique et mécanique des fluides
Contenu de la matière :
Chapitre 1 : Introduction à la Gazodynamique (1 semaine)

1. Concepts et relations thermodynamiques
2. Relations isentropiques d’un gaz parfait
3. Compressibilité et propagation des ondes sonores
4. Expression générale de la vitesse du son
5. Nombre de Mach et ondes de Mach
6. Ecoulements subsonique, transsonique, supersonique et hypersonique
Chapitre 2 : Ecoulement Isentropique 1D en Conduit à Section Variable (5 semaines)

1. Équations de base (continuité, quantité du mouvement, énergie)
2. Lois générales de l’écoulement isentropique : état générateur et état critique
3. Ecoulement 1D dans une conduite de section variable et théorème d’Hugoniot
4. Étude d’un écoulement dans une tuyère : convergente et convergente-divergente
5. Aperçu sur les diffuseurs subsoniques et supersoniques
Chapitre 3 : Ondes de Choc (4 semaines)

I- Ondes de Choc Normales
1. Équations de base (continuité, quantité du mouvement, énergie) et relation de Prandtl
2. Relations de l’onde de choc normale en fonction du nombre de mach
3. Cas limites : ondes de choc faibles, ondes de choc fortes
4. L’onde de choc normale mobile
5. Tube de Pitot en supersonique
II. Ondes de Choc Obliques
1. Notion sur les ondes de choc obliques
2. Equations de base et relation de Prandtl
3. Réflexion des ondes obliques
Chapitre 4 : Ecoulement Non Isentropique 1D en Conduit à Section Constante (5 semaines)
I. Ecoulementadiabatique avec frottement : Ecoulement de Fanno
1. Analyse de l’écoulement de Fanno et équations de base
2. Variation des caractéristiques d’écoulement en fonction du nombre de Mach
P a g e | 26
3. Coefficient du frottement et variation d’entropie

4. Onde de choc dans l’écoulement de Fanno
II. Ecoulement sans frottement et avec échange de chaleur : Ecoulement de Rayleigh
1. Analyse de l’écoulement de Rayleigh et équations de base
2. Variation des caractéristiques d’écoulement en fonction du nombre de Mach
3. Variation d’entropie
III. Ecoulement avec frottement et avec échange de chaleur
Ecoulement isotherme avec frottement
Mode d’évaluation:
Contrôle Continu : 40%, Examen : 60%.
Références bibliographiques :
1- Patrick Chassaing. Mécanique des Fluides, 3 ème édition, Cépaduès, Toulouse, 2010.
2- André Lallemand. Ecoulement monodimensionnel des fluides Compressibles,
Techniques de l'ingénieur Génie énergétique, B- 8- 165
3- F. M. White. Fluid Mechanics, 5th edition, McGraw-Hill, New York, 2003.
4- R. W. Fox and A. T. McDonald. Introduction to Fluid Mechanics, 5th edition, New York: Wiley,
1999.
5- J. D. Anderson. Modern Compressible Flow with Historical Perspective, 3rd edition, New York:
McGraw-Hill, 2003.
6- H. Liepmann and A. Roshko. Elements of Gas Dynamics,Dover Publications, Mineola, NY, 2001.
7- Genick Bar–Meir, Fundamentals of Compressible Fluid Mechanics, Minneapolis, MN 55414-2411,
2009
8- Robert d. Zucker, Oscar Biblarz, Fundamentals Of Gas Dynamics,JOHN WILEY & SONS, 2002
9- Patrick Oosthuizen, William Carrascallen, Compressible Fluid Flow, McGraw-Hill, 1997
10- Klaus Hoffmann, Computational Fluid Dynamics, Volume II, EES, 4th edition, 2000

3UpIDFH
Préface :
/D dynamique des gaz FRQQXH VRXV ©gazodynamiqueª HVW XQ GRPDLQH WUqV
YDVWH TXL D SRXU REMHFWLI WKpRULTXH O¶pWXGH GHV pFRXOHPHQWV FRPSUHVVLEOHV j JUDQGHV
YLWHVVHV &HV W\SHV G¶pFRXOHPHQWV VRQW OH SOXV VRXYHQW UHQFRQWUpV GDQV OH GRPDLQH
SUDWLTXH GH O¶LQGXVWULH DpURQDXWLTXH HW VSDWLDOH 2Q SHXW WURXYHU VRQ DSSOLFDWLRQ HQ
7XUERPDFKLQHVGHV0DFKLQHVj)OXLGHV&RPSUHVVLEOHV
/HSUpVHQWPRGXOHWUDLWHVHXOHPHQWO¶DSSURFKHPRQRGLPHQVLRQQHOOH©'ªGHV
pFRXOHPHQWVFRPSUHVVLEOHVGHVJD]SDUIDLWV
'DQVOHFDGUHGHFRPSOpWHUOHXUEDJDJHWKpRULTXHHWPHWWUHHQWUHOHXUVPDLQVXQ
JXLGH IDFLOH j XWLOLVHU RQ SURSRVH GDQV OH SUpVHQW SRO\FRSLp SRXU QRV pWXGLDQWV GH
0DVWHU GH *pQLH 0pFDQLTXH XQ SUpFLV GH OD WKpRULH GH '\QDPLTXH GHV *D] 2Q D
HVVD\p GH WRXFKHU j WRXV OHV SRLQWVPHQWLRQQpVSDUOHSURJUDPPHRIILFLHOGX0DVWHU
*pQLH0pFDQLTXH2SWLRQ(QHUJpWLTXH &LMRLQW
8QUHFXHLOG¶H[HUFLFHVDpWpSURSRVpDILQG¶DLGHUO¶XWLOLVDWHXUGHFHSRO\FRSLpGH
PDLWULVHUOHVQRWLRQVSUpVHQWpHV'HVUpSRQVHVDEUpJpHVVHURQWGRQQpHVHQDQQH[H%
1pDQPRLQV FH SRO\FRSLp SHXW rWUH DXVVL XWLOLVp SDU OHV pWXGLDQWV GHV DXWUHV
ILOLqUHVV¶LQWpUHVVDQWjO¶pQHUJpWLTXH
3DU DLOOHXUV RQ QRWH j QRV FKHUV OHFWHXUV TXH FH FRXUV HVW GpMj SXEOLp FRPPH
©FRXUV HQ OLJQHª VXU OD SODWHIRUPH 0RRGOH GH O¶8QLYHUVLWp GH 0pGpD GHSXLV HW
UpDFWXDOLVpHQ&RXUV'\QDPLTXHGHVJD] XQLYPHGHDG]
(QILQMHVHUDLWUqVUHFRQQDLVVDQWjWRXVFHX[TXLYHXOHQWDPpOLRUHUODTXDOLWpGH
FHSRO\FRSLpSDUOHXUVFULWLTXHVHWUHPDUTXHV
'U$//,&+(0RXQLU
0DLWUHGH&RQIHUHQFHV$
8QLYHUVLWpGH0HGHD)DFXOWpGH7HFKQRORJLH
'HSDUWHPHQWGH*pQLH0pFDQLTXH
DOOLFKHPRXQLU#XQLYPHGHDG]
DOOLFKHBP#\DKRRIU

3DJH
,QWURGXFWLRQjODG\QDPLTXHGHVJD]'U$//,&+(0RXQLU8QLY0HGHD

3UpIDFH
3DJH

3UpIDFH

3DJH
Chapitre I Introduction
Chapitre I: Introduction Générale
Objectifs du chapitre:
Rappel :
x Des propriétés des gaz parfaits

x Des lois thermodynamiques
x Les équations régissantes les écoulements des fluides.
Page 3
Introduction à la dynamique des gaz Dr ALLICHE Mounir, Univ Medea
Toute matière qui ne peut pas rester en équilibre statique sous la contrainte de
cisaillement est appelée fluide. Ceux qui vérifient la loi de Newton de viscosité, on les
appelle fluides Newtoniens. L’autre type ne vérifie pas la loi de Newton, on les appelle
fluides non Newtoniens.
Les fluides peuvent être aussi divisés en deux classes : les gaz et les liquides. Un
liquide se caractérise d’un gaz par le fait qu’une masse liquide mise dans un récipient
(dont le volume est supérieur à celui de cette masse) aura une surface libre plane s’il y a
équilibre statique. Tandis qu’une masse de gaz occupera tout le volume du récipient
puisqu’un gaz ne possède pas de surface libre. En outre, la densité des liquides ne varie
pas d’une manière très sensible et ils peuvent être considérés comme incompressibles.
Par contre, la densité d’un gaz varie très sensiblement, on dit que les gaz sont
compressibles. D’autre part, un fluide peut être idéal ou réel. Un liquide idéal est
incompressible et ne possède pas de viscosité (cas imaginaire). Un gaz idéal a une
viscosité nulle, n’est pas conducteur thermique et il est compressible (cas imaginaire).
Quand un fluide se déplace à des vitesses comparables à sa vitesse du son. Le
changement de densité devient important et l’écoulement est nommé compressible.
L’effet de la compressibilité est négligeable pour les petits nombres de Mach.
Ce cours traitera les écoulements compressibles, qui ont un nombre de Mach
supérieur à 0.3 et alors les changements de densité ne peuvent pas être négligés. Le sujet
est aussi appelé dynamique des gaz, parce qu’un rapport de pression de seulement de 2/1
peut causer un écoulement sonique. Par contre, les liquides exigent des très hautes
pressions pour générer des vitesses soniques.
1.1 Notions de base :
a) Gaz parfait :
Un gaz parfait est un gaz qui vérifie, sous d’autres conditions, l’équation
suivante :
p
pv R T T (1.1)
U m
avec p : la pression du gaz, U : masse volumique, T : température absolue du gaz, v :
volume spécifique, R : constante caractéristique du gaz, : constante universelle des
gaz parfaits et m : la masse molaire. La relation (1.1) est appelée loi des gaz parfait.
Page 4
Quoique le gaz parfait soit un état imaginaire, on remarque que les gaz
rencontrés dans la réalité peuvent être considérés comme tels.
b/Gaz réels:
Un gaz réel est le gaz rencontré dans la nature et ne vérifie pas la loi (1) mais
s'approchent de l'état parfait dans les hautes températures et les basses pressions
relativement au point critique. Un gaz réel à une viscosité. Il est compressible. Sa
température spécifique est variable selon les conditions de température et de pression.
1-2-Propriétés des gaz :
On peut classer les propriétés des gaz en trois catégories:

- Propriétés de diffusion
- Propriétés physiques
- Propriétés thermodynamiques
1-2-1- Propriétés de diffusion:
Ce sont des propriétés qui indiquent la diffusion de quelques états déterminés

dans le gaz à cause de la progression d'autres propriétés dans ce gaz. Parmi ces
propriétés, on cite principalement;
A:Viscosité dynamique(μ) ou viscosité absolue:
C'est l'indice de mesure de la diffusion de la quantité de mouvement par la

contrainte de cisaillement τ entre les couches du gaz à cause de l'évolution de la vitesse
du/dy.
La relation qui lie ces quantités est connue sous l'appellation de«loi de
Newton»de viscosité;
W
P , μ[kg/m.s] ou [N.s/m²] (1.2)
§ du ·
¨¨ ¸¸
© dy ¹
μ: indique la résistance du gaz vis-à-vis de la contrainte de cisaillement.
B:Viscosité cinématique (ν):
C'est l'indice de diffus ion de l'énergie cinétique entre les couches du gaz à cause de
l'évolution de leur vitesse:
Page 5
W
U
Q (1.3)
§ du ·
¨¨ ¸¸
© dy ¹
Ou encore:
ν = μ/ρ (1.3bis)
ν[m²/s].
C:le coefficient de la conduction thermique spécifique(k):
C'est l'indice de la diffusion du flux thermique entre les couches du gaz à cause de
l'évolution de leur températures (dT/dy):
K= Q/A/(dT /dy) (1.4)
K[KJ/m.°K]
On peut connaître les propriétés de diffusion du gaz par des procédés expérimentaux.
1-2-2- Propriétés physiques:
Ce sont des propriétés qui montrent l'état physique du gaz. On les mesure
directement à l'aide d'instruments.
La pression (P):
C'est la résultante des forces de cohésion des molécules du gaz sur l'unité de
surface de la paroi solide injectée dans le gaz. Elle est proportionnelle avec la somme
des énergies cinétiques de ces particules dans une unité de volume du gaz(homogène).
On la mesure à l'aide d'un manomètre.
P[N/m²]
La température (T):
C'est la mesure du niveau de l'énergie interne du gaz.On la mesure par le
thermomètre. Son unité est le [°C].La température absolue est définie par:
T= t + 273, son unité [k].
La densité (ρ):
C'est la mesure de la concentration de la masse dans une unité de volume d'un
gaz. Sur unité est [kg/m³].Le poids spécifique est défini comme:
ω = ρg [N/m³]. (1.5)
Page 6
1-2-3- Propriétés thermodynamiques:

C'est l'ensemble des propriétés qui indiquent l'état thermodynamique d'un
gaz.On peut les calculer par la connaissance des propriétés physiques. On cite
principalement:
A:l'énergie interne spécifique (e) ou (u):
C’est l'énergie contenue dans une unité de masse d'un gaz en état statique. Ou
encore,c'est la somme des énergie cinétique aléatoire des particules d'une masse ou
encore c'est la quantité de chaleur transportée à une unité de masse d'un gaz sous un
volume constant.
Puisque la température du gaz est l'indice du niveau de l'énergie cinétique des
particules,on peut considérer que la température est aussi l'indice du niveau de l'énergie
interne càde↑si t↑ et e↓ si t↓. e[KJ/kg].
B:le travail d'écoulent spécifique:
C'est l'énergie qu'on doit fournir pour faire entrer une unité de masse à travers
des frontières du volume étudié dans le domaine de l'écoulement ou de la faire sortir de
ce volume sans autant changer la pression ou la température du gaz.
Il est exprimé par:
W=Pυ=( P/ρ ); [KJ/kg]. (1.6)
C:l'enthalpie (h):(l'énergie emmagasinée spécifique)
C'est la somme du travail d'écoulement et de l'énergie interne par unité de masse

d'un gaz en mouvement. Elle se mesure par la quantité de chaleur transmise sous une
pression constante à une unité de masse du gaz:
h = e +Pυ [kJ/kg]. (1.7)
D:la chaleur spécifique:
C'est la moyenne de la variation de l'énergie du fluide spécifique par rapport à la

variation de la température:
§ de · dq
*à volume constant : cv ¨ ¸ (1.8)
© dT ¹ dT v
§ dh · dq
*à pression constante c p ¨ ¸ (1.9)
© dT ¹ dT P
Page 7
Où : q=quantité de chaleur transportée du ou à l'unité de masse C [KJ/Kg.°C]

Pour un gaz parfait,on a:
JR R
h c pT .T (1.10) e cvT .T (1.11)
J 1 J 1
Avec : γ=Cp/Cv et R=Cp - Cv
E:L'entropie (S):
Selon le second principe de la thermodynamique, il existe une variable

d'états,appelée entropie, qui pour un processus réversible s'exprime par:
dS= dq/T (1.12)
ou encore;
TdS= de+Pd(1/ρ)=dh-(1/ρ)dP (1.13)
On peut déduire des relations générales entre S et h ou S et e. Ces relations sont
appelées relations de réciprocité:
§ wS · 1 § wh · cp
¨ ¸ ¨ ¸
© wT ¹ P T © wT ¹ P T
§ wS · 1 ª§ wh · 1º
¨ ¸ «¨ ¸ »
© wP ¹ T T ¬© wP ¹ T U ¼
§ wS · 1 § we · cv
¨ ¸ ¨ ¸ (1.14)
© wT ¹ U T © wT ¹ U T
§ wS · 1 ª§ we · Pº
¨¨ ¸¸ «¨¨ ¸¸ »
© wU ¹ T T ¬© wU ¹ T U ² ¼
Avec:
§ wh · 1 T § wU ·
®¨ ¸ ¨ ¸
¯© wP ¹ T U U ² © wT ¹ P
(1.15)
§ we · 1 ª § wP · º
®¨
¨ ¸
¸ «P T ¨ ¸ »
¯© wU ¹T U ² ¬« © wT ¹ U ¼
»
On a aussi:
Page 8
§ wh · 1 § wU ·
®¨ ¸ ¨ ¸
¯© wP ¹T U ² © wT ¹ P
(1.16)
§ wS · 1 § wP ·
®¨
¨ ¸
¸ ¨ ¸
¯© wU ¹T U ² © wT ¹ U
1-3-Application du premier principe de la thermodynamique:
D'après le 1er Principe relatif à la conservation d'énergie, on peut considérer

qu'unfluide recevant du travail et de la chaleur, soit ses énergies sous forme d'énergie
interne et cinétique:
1
6Q 6W u 2 u1 V2 ² V1 ² (1.17)
2
Chaleur +travail=variation d'énergie interne+ énergie cinétique
Pour cela considérons une conduite telle qu'elle indique la figure (1) suivante:
QC
P1 A1 '
V1 B1 '
U1
B2 '
A2 '
P2
A1 A1
V2
B1
S1 U2
WC
Z1 B2
A2
Position1 Z2 Position2
S2
Figure 1
1
P:pression, ρ:masse volumique, Q : volume massique, V:vitesse
U
S:section, Z:altitude
Pendant le temps "dt", le fluide contenu dans A1 A'1 A2 A' 2 vient en B1 B1 ' B2 B2 ' .
On sait que:
Page 9
.
dm .
m U1 S1v1 U 2 S 2 v2 (1.18)
dt
tout se passe comme si la masse "dm" toujours constante se trouve en position 2 mais
avec de nouvelles caractéristiques thermodynamiques.
Pour cette masse "dm", on pourra toujours calculer les diverses énergies;
- Celle causée par la pression (énergie de transvasement): dm P1v1 P2 v2
- Energie due aux forces de pesanteur: ce type d'énergie peut être négligé pour les
gaz dans la conception des machines ou le ∆Z est très petit.
Cependant pour des cas particuliers de variation des mines, forages et hautes
cheminées, cette variation ne peut en aucun cas être négligée: dm Z1 Z 2 g
- Travail mécanique entre l'entée et la sortie; dWm
- Quantité de chaleur reçue on dégagée par unité de masse; Q12 dm
-Variation de l'énergie interne; dm u 2 u1
§1 1 ·
-Variation de l'énergie cinétique; dm¨ V2 ² V1 ² ¸
©2 2 ¹
Ainsi, le 1erPrincipe appliqué à ce système s'écrira;

1
Q12 W12 P1v1 P2 v 2 g Z 2 Z 1 u 2 u1 V2 ² V1 ² (1.19)
2
Ou
1
Q12 W12 u 2 u1 P2 v 2 P1v1 g Z 2 Z1 V2 ² V1 ² (1.20)
2
Remarques :
a) Q12 ,W12 sont des valeurs algébriques {positif s'il y a appel de chaleur ou de travail
(compression ), négatif situation inverse(turbine)}.
Cette équation peut être rendue plus générale en y faisant figurer l'énergie chimique
C12 C1 C2 dans le membre gauche:
1 1
C1 u1 P1v1 gZ1 V1 ² We Qe C 2 u 2 P2 v 2 gZ 2 V2 ² WS QS (1.21)
2 2
Chacun des termes est exprimé en unité d'énergie [J/kg].
En ingénierie on préfère suivent regrouper u Pv h (Enthalpie)
Page 10
b) On notera que le 1er Principe de la thermodynamique se bornant à constater la

conservation de l'énergie sans se préoccuper de sa qualité (mécanique,
thermique…) ne fait aucune destination de qualité entre les deux forme des
énergies.
L'énergie thermique étant l'énergie d'agitation désordonnée des molécules quant à
l'énergie mécanique, elle correspond à une énergie de mouvement ou toutes les
particules sont orientées dans le même sens. La transformation du travail en chaleur est
facile car il suffit de laisser agir les chocs moléculaires pour que l'énergie parfaitement
ordonnée soit complètement désordonnée. La modification inversée étant au contraire
très délicate.On verra en turbomachines qu'elle ne peut s'effectuer qu'avec un rendement
inférieur à 1.
1-4- Application aux écoulements sans échanges d'énergie;

mécanique,chimique,..et de la chaleur avec l'extérieur:
L'équation générale sous ces hypothèses se simplifie pour devenir :
1 1
u1 P1v1 gZ1 V1 ² u 2 P2 v 2 gZ 2 V2 ² (1.22)
2 2
Ou sous la forme:
1
u Pv V ² gZ cte (1.23)
2
Cette équation est valable sur toute la longueur de la canalisation.
Cas d'un fluide incompressible:
Cette hypothèse peut être valable pour le gaz mais dans des conditions de
vitesses impliquant des variations de masse volumique négligeable permettant
d'admettre qu'on a une situation d'incompressibilité.
Dans ce cas, on a;
u cte
U cte (Cas particulier des ventilateurs)
'T 0
Ceci entraîne l'équation de BERNOULLI qui s'écrit sous forme de hauteur de fluide
considéré:
P V²
Z cte (1.24)
Ug 2g
Avec:
Page 11
P
: Hauteur équivalente à la pression statique
Ug
V²
: Hauteur équivalente à la pression dynamique
2g
P
Z : Hauteur piézométrique
Ug
P V²
Z : L'équivalent à la pression totale
Ug 2g
Remarques :
- Les pressions statistiques et totales sont mesurables.
-la pression dynamique n'est pas mesurable,elle est calculable.
-Quand un gaz est soumis à des faibles variations de pressions(cas du ventilateur) on le
considère comme incompressible.
De plus, pour les gaz les variations de cotes sont généralement négligeables
devant les variations dues aux vitesses. On néglige ainsi les termes en ρgz et l'équation
de Bernoulli prend la forme:
1
P UV ² cte Pt (1.25)
2
P:pression statique, Pt : Pression totale ou d'arrêt, ( 1/2)ρV²:pression dynamique
1-5-Application du second principe de la thermodynamique:
Selon le second principe de la thermodynamique, il existe une variable d'état S,

appelée entropie, qui pour un processus réversibles s'exprime par:
dS= dq/T (1.26)
Ou encore:
TdS= di+Pd(1/ρ)= dh-(1/ρ).dP(1.27)
On peut déduire des relations générales entre S et h ou S et i.
Ces relations, dites relations de réciprocité, sont données par (1.14), (1.15) et(1.16).
1-6-Rappels de physique du son:
1-6-1-Vitesse du son:
Ce qui différencie un écoulement compressible d'un écoulement incompressible,
principalement le fait que des petites variations de masses volumique, de pression et de
Page 12
vitesse se propagent à l'intérieure de fluide à une certaine vitesse qui est la célérité du
son.
Pour caractériser cette célérité, considérons une petite perturbation adiabatique
réversible(donc isentropique) qui se déplace dans le fluide initialement au repos contenu
dans un tube de section A constante(exemple système =cylindre –piston)
u =0
P+dP P
du ρ +dρ ρ
Onde acoustique
Figure 2
L'onde acoustique, en traversant le fluide à l'arrêt(après déplacement du piston), va

l'ébranler en lui induisant une vitesse V<<<a et en lui causant une variation de pression
et de masse volumique:
P o P dP
U o U dU
Ce système peut être décrit par les équations de continuité, de quantité de mouvement,
l'état et d'énergie.
- Equation de continuité: appliquée au système établi d'onde acoustique s'écrit;
div U a 0
Ou : (1.28)
§ wa · § wU ·
U¨ ¸ a¨ ¸ 0
© wu ¹ © wx ¹
- Equation de quantité de mouvements'écrit;

§ wa · 1 § wP ·
a¨ ¸ ¨ ¸ (1.29)
© wx ¹ U © wx ¹
L'équation (1.29) donne:
Page 13
§ wa · 1 § wP ·
¨ ¸ ¨ ¸ (1.30)
© wx ¹ Ua © wx ¹
D’où:
§ 1 wP · § wU · § wP ·
U ¨¨ ¸¸ a¨ ¸ 0 a² ¨¨ ¸¸ (1.31)
© Ua wx ¹ © wx ¹ © wU ¹
- Equation d'énergie:
Si on suppose le fluide comme gaz parfait et que l'écoulement est isentropique.
On aura:
§P ·
¨¨ rT ¸¸
©U ¹
(1.32)
§ P ·
¨¨ J cte ¸¸
©U ¹
La différentiation logarithmique de la seconde équation donne:
dP dU
J 0
P U
Donc : (1.33)
dP P
J
dU S
U
Ainsi la vitesse du son dans le gaz parfait équivalent à:
JP
a JrT (1.34)
U
1-6-2-Nombre de Mach:
Le nombre de Mach est par définition le rapport entre la vitesse locale de la
particule du fluide(V) et la vitesse du son dans le fluide au même point considéré(a)
V
M (1.35)
a
Le nombre de Mach"M" est un nombre additionnel permettant de caractériser en
plus du type de l'écoulement, la compressibilité du fluide.Les Aérodynamiciens font
Page 14
spécialement une distinction entre les différents intervalles du nombre de Mach.Ainsi, la

classification suivante est utilisée:
M<0.3 L'écoulement est incompressible, où les effets sont négligeable.

0.3<M<0.8 L'écoulement est subsonique, où les effets de la densité sont
importants mais les ondes de choc n'apparaissent pas.
0.8< M <1.2 L'écoulement esttranssonique,où les ondes de choc
apparaissent premièrement, divisant les régions subsonique et
supersonique de l'écoulement.
1.2< M <3 L'écoulement est supersonique, où les ondes de choc sont
présentées mais il n'y a pas de régions supersoniques.
3< M L'écoulement est hypersonique,où des ondes de choc,et d'autres
changements dans l'écoulement sont spécialement forts.
Les valeurs numériques données ci-dessus sont seulement des guides

approximatifs.Ces quatre catégories d'écoulement sont valables pour les écoulements
externes de hautes vitesses. Pour les écoulements internes dans les conduites, la plus
importante question est simplement quand l'écoulement est subsonique (M<1) ou
supersonique (M>1).
1-7-Les équations fondamentales appliquées aux fluides compressibles:

1-7-1-Equation de continuité:
L'équation de continuité pour un écoulement unidimensionnel d'un fluide
compressible est;(γ:poids spécifique)
G JAV cte (1.36)
ou
m U A V cte
Avec G et ݉ሶ sont respectivement débit poids et débit massique.
1-7-2-Equation d'énergie:
Pour un écoulement unidimensionnel d'un fluide compressible et il n'existe
aucune machine entre 1 et 2, l'expression de l’énergie est :
Page 15
V1 ² V2 ²
h1 QH h2 (1.37)
2g 2g
்
Où :݄ത ൌ ݅ ൅
ఊ
1-7-3-Equation de quantité de mouvement:

F U 2 Q2V2 U1Q1V1 (1.38)
1-7-4-Equation d'Euler:
Pour un écoulement d'un fluide compressible dans une conduite circulaire; on a :
dU 2Wds
VdV (1.39)
U Ur
Avec :
τ: constant de cisaillement
r: rayon de section
ds: longueur de la section de la conduite
ds
r τ
V Figure 3
Page 16
DEPARTEMENT GENIE MECANIQUEModule : Gazodynamique

Série N° 1
Ex.01 :
L’air circule à l’intérieur d’une conduite avec un débit de 0.15kg/s. une section
de cette conduite est refroidie à l’aide du nitrogène liquide. La chaleur cédée par l’air est
de 15kJ/s.
On donne : T1= 440K, P1=188kPa, V1=210m/s
T2= 351K, P2= 213Pa
Calculer : a) la section A
b) la variation ∆h, ∆u, ∆s
Ex.02 :
On a 5kg d’air qui se trouve dans un système fermé subit une detente
isentropique de P1=300MPa à T1=60°C jusqu’à une pression P2=450kPa.
Calculer T2 ?
Ex.03 :
Dans un système fermé, un gaz subit les transformations suivantes :
1-2 : compression isothermiquereversible
2-3 : apport de chaleur isochore (à v=cte)
3-4 : transformation isobarque (à P=cte)
4-1 : transformation isentropique adiabatique reversible.
Tracer les diagrammes P-v et T-s.
Ex.04 :
Considérons une compression infinitésimale d’un gaz parfait dans un
diagramme T-s. P+dP
T
1-2s : compression isentropique 2
2s
1-2 : compression polytropique
P
On donne : γ=k=Cp/Cv , ηpc= ∆hs/∆h = (h2s-h1)/(h2-h1)
ௗ் ௞ିଵ ௗ௉
a) Montrer que: ൌ
் ௞஗౦ౙ ௉ 1
ುమ
ೖషభ s
ሺ ሻ ೖ ିଵ
b) Montrer que: ߟ௖ ൌ ುభ
ೖషభ భ
ುమ Ǥ
ሺ ሻ ೖ ಏ౦ౙ ିଵ
ುభ
Page 17
Ex.05 :
Considérons une détente infinitésimale d’un gaz parfait dans un diagramme T-s.
3-4s : detente isentropique P3
T
3-4 : detente polytropique
3
On donne :γ=k=Cp/Cv , ηpd= ∆h/∆hs = (h3-h4)/(h3-h4s)
P4
ௗ் ௞ିଵ ௗ௉
a) Montrer que: ൌ Ʉ୮ୢ
் ௞ ௉
4
ೖషభ
ುర ಏ
ሺ ሻ ೖ ౦ౚ ିଵ 4s
b) Montrer que: ߟௗ ൌ ುయ
ುర
ೖషభ
ሺ ሻ ೖ ିଵ
Ǥ s
ುయ
Page 18
Chapitre III Ecoulement dans une conduite à Section variable
Chapitre 3. Ecoulement dans une conduite à Section variable –

Théorème d'Hugoniot
Objectifs du chapitre :
x Comprendre L’Effet de variation de la section de la conduite

sur la vitesse de l’écoulement
x Comprendre L’Effet de variation de la vitesse sur la pression
de l’écoulement
x Connaitre les conditions d’écoulement dans une conduite
convergente
x Connaitre les conditions d’écoulement dans une conduite
divergente
x Connaitre les conditions d’écoulement dans une tuyère
convergente-divergente
Page 31 Introduction à la dynamique des gaz Dr ALLICHE Mounir, Univ Medea

Ecoulement dans une conduite à Section variable –

Théorème d'Hugoniot
Introduction :
On présentera dans ce cours :
- Equations de base.
- Conditions à l'état générateur.
- Table de calcul.
- Tuyère convergente-divergente.
3.1- Equations de base:

On suppose qu'on a un écoulement isentropique dans des conduites à section
variable.
Soient les hypothèses suivantes:
- Ecoulement monodimensionnel et permanent.
- Pas de travail mécanique ou technique Wu 0
- Adiabatique: Q 0.
- On néglige l'énergie potentielle gz=0.
- Irréversible GW 0.
Les équations de base sont;

1- Equation de continuité : d U.VA 0. (3-1)
2- Equation d'énergie : dht = dh + d(V2/2) = 0 (3.2)
dU
3- Deuxième loi de la thermodynamique : TdS dh (3-3)
U
L'équation de continuité après différentiation logarithmique donne:
log( U V A) log U log V log A
D’où:
dU dV dA
0 (3-4)
U V A
L'équation de Barré-Saint Venant donne:
1 1
Ht H V² u Pv V ² (3.5)
2 2

En dérivant cette équation, on obtient :

dH t du Pdv vdP VdV (3.5a)
0 du Pdv vdP VdV (3.5b)

dq 0
Ainsi:
dP
VdV 0 (3-6)
U
D'autre part, on a: (voir chapitre 1)
dP JP
a ² dP a ² dU (3-7)
U S
U
D’où après substitution dans (3.6) :
dU dU dV V ² dV dV
a² VdV 0 V M ² (3.7b)
U U a² a² V V
Finalement, on obtient après injection dans (3.4):
dA dV
1 M ² 0 (3-8)
A V
La relation (3.8) représente le premier théorème d'Hugoniot qui lie la variation de la
section avec celle de la vitesse.
2eme Théorème d'Hugoniot:

Il lie la variation de la vitesse à celle de la pression:
dP dV 1 dP
3.6 VdV D’où:
U V V² U
JP
Or : V²=M²a² , Avec a²
U
D’où:
dV 1 dP

V M ²a² U
dV 1 1 dP 1 dP

V M ² JP U JM ² P
U

Finalement:
dV 1 dP
(3-9)
V JM ² P
Les relations d'Hugoniot permettent de déduire certaines règles
fondamentales pour les écoulements de fluide compressible et isentropique dans
n'importe quel type de canalisation (subsonique, supersonique):
1°/- En écoulement subsonique (M<1), les variations de vitesse et de section sont de
signes contraires: c.à.d. si A↑ alors V↓ et vis versa.
2°/- Dans le cas supersonique (M>1), dS et dv sont de même signe. Si V↑ alors A↑ et
vis versa.
3°/- Le second théorème d'Hygoniot montre clairement que les vitesses et les pressions
varient toujours en sens inverses indépendamment du type de l'écoulement : V↑ alors P↓
et vis versa.
4°/- La vitesse du gaz ne peut être égale à la vitesse du son (M=1) que pour une section
de canalisation où l'aire présente un max ou un minimum:
En effet : V=a pour dA/A=0.
Physiquement, seul le minimum convient.
Section A Densité ρ Vitesse V Pression P
Subsonique (M<1) ↓ ↓ ↑ ↓
Supersonique (M>1) ↑ ↓ ↑ ↓
3.2- Débit masse dans une canalisation:

On sait que le débit masse s'exprime par la relation suivante:
Q U V A (3.10a)
Qui peut s'écrire:
J 1
ª º
2J Pi « § P · J »
Q U A 1 ¨ ¸ (3.10b)
J 1 U i « ¨© Pi ¸¹ »
«¬ »¼
En faisant intervenir la loi de Poisson, on obtient:

1
§P· J
U U i ¨¨ ¸¸ ; Pi=connue, P=mesurable
© Pi ¹

1 J 1
ª º
§ P ·J 2J § Pi ·« § P · J »
Q U i A¨¨ ¸¸ ¨ ¸¸ 1 ¨¨ ¸¸
© Pi ¹ J 1 ¨© U i ¹ « © Pi ¹ »
¬ ¼
1 J 1
2J Pi § P ·J §P· J
Ui A ¨¨ ¸¸ 1 ¨¨ ¸¸ (3-10c)
J 1 Ui © Pi ¹ © Pi ¹
2J Pi
Avec VMAX (3.11a)
J 1 Ui
1 J 1
§P· J
§P· J
Et 8 ¨¨ ¸¸ 1 ¨¨ ¸¸ (3.11b)
© Pi ¹ © Pi ¹
Finalement:
Q U i A VMAX 8 (3-11c)
La section et la fonction Y varient en sens inverse : Y↑ , A↓ pour un même débit masse

et vis versa.
Y q/A
Fluide
1
J 1 J 1
2
col incompressible
§
¨ 2 J 1·
¸ §
¨ ·
© ¹ © J 1¸
¹
Fluide
supersonique compressible
subsonique
P P
Pi
0 Pc 1 Pi Pc 1
Pi Pi
Figure(1)-variation de la Figure(2)-débit d'une
Fonction Y Tuyère convergente
Si l'on maintient constante Pi , U i de l'état générateur (réservoir amant), le débit reste

constant tout le long de la conduite (voir figure 2).
Y atteint sa valeur max dans la section d'aire maximal ou col de la tuyère. Y=f (P/Pi)
s'annule pour P=0, P=Pi (voir figure 1).

Y passe par un maximum pour la valeur P Pc (pression critique obtenue par

w8
dérivation et annulation de la dérivée) on obtient :
wP
J
Pc § 2 · J 1
¨¨ ¸¸ (3-12)
Pi © J 1¹
Le théorème d'Hygoniot permet de déduire que la vitesse d'écoulement au col

est égale a la vitesse du son dans cette section:
Tc Tc
Vcol acol JrTc JrTi JrTi
Ti Ti
(3.13a)
1 1 1
§ J 1 2 · § J 1 · § J 1·
Vcol acol ai ¨1 Mc ¸ ai ¨ 1¸ ai ¨ ¸
© 2 ¹ © 2 ¹ © 2 ¹
D’où:
2 (3-13b)
Vcol acol ai
J 1
Le débit masse en régime critique s'écrit sous la forme:
qc U c Vc Ac (3.14a)
1
Ui § J 1 · J 1
Or: car ¨1 n² ¸
U © 2 ¹
D’où:
1
§ 2 · J 1 2
qc U i ¨¨ ¸¸ ai Ac
© J 1¹ J 1
1 1

§ 2 · J 1 2
qc Ac U i ai ¨¨ ¸¸
© J 1¹
J 1
§ 2 · 2 J 1
qc Ac U c ai ¨¨ ¸¸ (3-14b)
© J 1¹
Cette valeur représente un max pour le débit de la tuyère.

Caractéristiques de l'état critique: ( pour γ=1.405 càd pour l'air)

J
½
Pc § 2 · 1 J
¨ ¸ 0.527 °
Pi ¨© J 1 ¸¹ °
1
°
°
U c § 2 · J 1 °
¨¨ ¸¸ 0.634¾ (3-15)
Ui © J 1¹ °
1 °
ac § 2 · 2 °
¨¨ ¸¸ 0.912 °
ai © J 1 ¹ °¿
D’où : q c 0.685 Pi U i Ac (3-16)
Ou encore:
P
q c 0.685 i Ac (3-17)
rTi
Si vitesse du son est atteinte au col, le débit masse est donc proportionnel à la pression
du réservoir et inversement proportionnel à la racine carrée de la température du
générateur.
3-3- Ecoulement dans une tuyère convergente :
Soit une tuyère convergente alimentée par un réservoir ( Pi , U i , Ti ) la vitesse à

l'entrée de la tuyère est négligeable. La tuyère débite dans une enceinte ou régime une
pression PA que nous pouvons faire varier. (Figure 3).
x Si PA p le débit masse n et atteint un max q c au col pour une détente de rapport

J
PA Pc § 2 · J 1
¨¨ ¸¸ 0.527
Pi Pi © J 1¹ (3.18)
A ce moment, la vitesse de l'écoulement dans la section de sortie est la vitesse du son.
x Si PA continue à décroître ( plus bas que Pc ) , l'application de la relation du débit
( qi U i S VMAX 8 ) donnerait une diminution du débit mathématiquement ; ce

qui est physiquement incorrecte. Le débit reste constant et la pression P dans la
section de sortie ne suit pas l'abaissement de PA ; elle reste égale à Pc .
L'écoulement en aval du col n'est plus isentropique. Les relations précédentes ne
sont plus applicables. La détente de Pc à PA se fait de façon irréversible.

On dit que la tuyère amorcée quand les conditions au col sont les conditions critiques.
Pour les courbes (1, 2, 3), le régime est isentropique pour la courbe 4, le régime est
irréversible.
réservoir Tuyère enceinte valve

convergente
Pi Pt
Ti PA
PA Tt
Ui
V 0 V 0 Pc
P / Pt col
Pc PA
Valve
1fermé
Pc PA
2
3 4
m
4 3 2
qc
1
P 1 PA / Pi
0.527
Figure (3) Pt

Exemple 1 :
Soit une tuyère convergente dont la section au col Ac =0,001m², est alimentée
par un réservoir. La pression dans l'enceinte est PA 591KPa .

Calculer M c etm
c au col sachant que : Pt 1MPa, Tt 60qC, Vt 0
Solution :
On a Pa/Pt =0.591>0.527
D'après les tables d'un écoulement isentropique, pour PA Pc ; M c 0.9 .
షം
ఊିଵ ംషభ
On a Pc/Pt = 0591 = ቀͳ ൅ ‫ܯ‬௖ଶ ቁ Î Mc ≈0.9
ଶ
D’autre part, on a :
ఘ ఘ ඥఊ௥்೎
݉ሶ௖ ൌ ߩ௖ ܸ௖ ‫ܣ‬௖ ൌ ఘ೎ ߩ௧ ‫ܯ‬௖ ܽ௖ ‫ܣ‬௖ ൌ ఘ೎ ߩ௧ ‫ܯ‬௖ ඥܶ௧ ‫ܣ‬௖
೟ ೟ ඥ்೟
D’où :
ఘ ௉௧ ்
݉ሶ௖ ൌ ቀఘ೎ቁ Ǥ ቀ ቁ Ǥ ‫ܯ‬௖ Ǥ ටߛ‫் ݎ‬೎ Ǥ ඥܶ௧ Ǥ ‫ܣ‬௖
೟ ௥்௧ ௧
AN :
r=287 J/kg.K, Tt=333°K
ఘ ்
Mc=0.9 Î ఘ೎ ൌ ͲǤ͸ͺ͹Ͳ , ்೎ ൌ ͲǤͺ͸Ͳ͸, (table isentropique de l’air)
೟ ೟
On trouve :
ଵǤ଴‫כ‬ଵ଴ల
݉ሶ௖ ൌ ͲǤ͸ͺ͹Ͳ ‫ כ‬ቀଶ଼଻‫כ‬ଷଷଷቁ ‫Ͳ כ‬Ǥͻ ‫ כ‬ξͳǤͶ ‫ʹ כ‬ͺ͹ ‫Ͳ כ‬Ǥͺ͸Ͳ͸ ‫ כ‬ξ͵͵͵ ‫Ͳ כ‬ǤͲͲͳ
࢑ࢍ
࢓ሶࢉ ൌ ૛Ǥ ૚ૢ૞૜૜ ࢙
؆ ૛Ǥ ૛૙࢑ࢍȀ࢙

3.4- Etude d’un Divergent :

Soit une tuyère convergente-divergente. Les points 1 et 2 sont les sections
d’entrée et de sortie du divergent. Soit PA la pression de l’enceinte dans laquelle débite
la tuyère.
1 2
Figure 4
Premier cas : V1 < a1

L’écoulement est analogue à celui des fluides incompressibles. Le long du
divergent, la pression augmente et la vitesse diminue. La pression PA est égale à P2.
Toute discontinuité de pression à la sortie provoquerait un branlement qui
pourrait remonter le courant et modifier les conditions d’écoulement.
Deuxième cas : V1 = a1
D’après le théorème d’Hugoniot, V1 ne peut être supérieure à a1 , la section en
(1) étant minimale. Dans ce cas, d’apres le theoreme d’Hugoniot, la vitesse continue à
croitre et la pression à diminuer le long du divergent. Les conditions en (1) sont celles
de l’etat crtique. Selon la valeur de la pression PA, trois cas seront à envisager :
a) Si P2=PA : le nombre de Mach à la sortie peut etre calculé à partir du theoreme de
Saint Venant :
ംషభ
்೔ ఊିଵ ௉
‫ܯ‬஺ ଶ ൌ ቀ ೔ ቁ
ം
ൌͳ൅ (3.19a)
்ಲ ଶ ௉ಲ
A la sortie, on a les relations suivantes :

ംషభ
ଶ ௉೔ ം
‫ܯ‬஺ ൌ ඨఊିଵ ቈቀ௉ ቁ െ ͳ቉ (3.19b)
ಲ
ംశభ
ംషభ
஺ಲ ଵ ଵା ெಲ మ మሺംషభሻ
஺೎
ൌ ቈ మ
ംశభ ቉ (3.19c)
ெಲ
మ

b) Si P2 > PA : le jet se dilate à la sortie et crée une onde de dépression dont la

vitesse est inferieure à celle du son en (2). Cette onde ne peut donc remonter le
courant et modifier l’écoulement.
c) Si P2 < PA : il se produit une onde de recompression à la sortie, provoquant une
élévation de température.
Si le rapport PA /P2 est suffisamment grand, l’élévation de température peut être
assez élevée pour que la vitesse locale du son devienne plus grande que V2.
L’onde de compression remonte alors le divergent. Mais quand elle se rapproche
du col, le taux de compression de l’onde diminue ainsi que sa température et sa
vitesse. Elle s’arrête alors en une section où sa vitesse est égale à celle du
courant. On obtient alors une onde de recompression stationnaire appelée « onde
de choc ». A travers cette onde, il y’a une brusque variation irréversible des
propriétés du fluide suivie d’une variation isentropique dans la suite du
divergent.
3.5- Etude de l’écoulement à travers une tuyère convergente-divergente :

Considérons une tuyère convergente-divergente (dite de Laval) où on va faire
décroitre la pression Ps de l’enceinte avale où débite la tuyère. (voir Figure 5)
PD PS
1) 1 : l’écoulement est subsonique dans toute la tuyère (courbe 1)
Pi Pi
Dans le convergent : V et P jusqu’à une valeur minimale au col sans pour autant
être la pression critique
Dans le divergent : V et P jusqu’à la pression avale PA
PS PD
2) : cas limite où la vitesse du son est atteinte au col pour valeurs
Pi Pi
inferieures (courbe 2). L’écoulement est partout subsonique sauf au col (il est
sonique).
PB PS PD
3) : En faisant décroitre encore PS , le débit ne varie plus et reste égale
Pi Pi Pi
à qc . Dans cette zone, il se produit une onde de choc stationnaire dans le divergent.
L’écoulement est supersonique entre le col et l’onde de choc, subsonique ensuite
(courbes 3 et 4). Il y’a donc deux types d’écoulements séparés par une onde de choc

qui crée un processus irréversible. Il se produit une discontinuité brusque de

pression, de vitesse et de la température statique. L’expérience montre que cette
onde de choc stationnaire et normale à l’axe de la tuyère. Au fur et à mesure que PS
diminue, cette onde s’éloigne du col.
PS PB
4) : cas limite où l’onde de choc se produit juste à la sortie de la tuyère.
Pi Pi
PA PS PB P
5) : le rapport S est insuffisant pour que l’onde de choc puisse
Pi Pi Pi Pi
remonter la tuyère. L’expérience montre que le plus souvent ce régime est instable
et qu’il se forme des décollements dans le divergent. Le régime devient destructeur.
PS PA
6) = : l’écoulement est supersonique dans toute la tuyère. Cas limite où la
Pi Pi
pression à la sortie est juste égale à la pression de l’enceinte. C’est la valeur pour
laquelle une tuyère isentropique fonctionne correctement. L’écoulement est
entièrement isentropique dans la tuyère. Les possibilités maximales de la détente
sont utilisées. C’est le cas idéal de la propulsion. Malheureusement, il est rare
qu’une tuyère soit parfaitement adaptée.
PS PA
7) 0< : l'écoulement est supersonique dans toute la tuyère .il se forme une
Pi Pi
onde de détente irréversible a la sortie (onde de choc oblique).
En conclusion, une tuyère supersonique donnée convergente –divergente ne

ps pA
fonctionne correctement que pour une seule valeur .
pi Pi
Remarque :
Une tuyère fonctionne dans de très mauvaises conditions si la pression avale est
incorrecte .une fraction appréciable de l'énergie cinétique est alors perdue dans les
ondes irréversibles

Figure 5- Classification des écoulements dans une tuyère convergente-divergente

Exemple2:
Pt 1MPa
Tt 350q K Pb 954 KPa Pb 954 KPa

3
As 10 m²
Mc 0.68
Mc
Pt Ms
Ac
Tt As
Calculer Tc , Pc , Vc , M c et Ac ?
Solution:
T P U A
Pour M c 0.68 Table
o 0.9154, 0.7338, 0.8016, 1.110
Tt Pt Ut A
D’où;
Tc 0.9154 Tt 320q K
Pc 0.7338Pt 0.7338MPa 734KPa
Uc 0.8016U t 7.98Kg / m 3
Pt 1 10 6
Ut 9.955kg / m 3
rTt 287 350
Pb Ps 954 As
0.954 Table
o M s 0.26, 2.317
Pt Pt 10 3 A
D' ou;
AS 10 3
A 4.32 10 4 m²
2.317 2.317
Ac
1.110 o Ac 1.11A # 4.7952 10 4 m² | 4.8 10 4 m²
A
Vc M c JRTc 0.68 1.4 287 350 255m / s

Il y a quatre combinaisons de changement de la surface de la section et le

nombre de Mach.
Tableau 1-
Géométrie de la conduite Subsonique M<1 Supersonique
V↑ V↓
P↓ P↑
Jet subsonique Diffuseur supersonique
dA/A <0
V↓ V↑
dA/A >0 P↑ P↓
Diffuseur subsonique Jet supersonique
Pour M=1, puisque la vitesse infinie est impossible, l’équation (3-8) indique que
dV n’est finie que lorsque dA=0, ce qui correspond à une section minimale (col) ou
maximale.
Amax
Amin M<1
Subsonique subsonique
M=1 Supersonique supersonique

M>1
Figure 6-
Le col ou le convergent–divergent peut facilement accélérer l'écoulement

subsonique ou supersonique. Ceci est la seule façon qui peut être crée par la détente du
gaz de l'état stagné du réservoir.

L'autre section divergente- convergente n'est pas possible, par ce que le nombre
de Mach dans ce cas s'éloigne de la condition sonique au-lieu de s'approcher d'elle.
Malgré que l'écoulement supersonique en aval de la tuyère exige un col
sonique l'inverse n'est pas vrai. Un gaz compressible peut passer à travers un col sans
être sonique.
Col
M<1 M<1 M<1
Amin
Figure 7.
3.2. Relations des gaz parfaits :

On peut utiliser les équations de gaz parfait et l'écoulement isentropique pour
changer l'équation de continuité en une équation algébrique en fonction de la section et
le nombre de Mach.
On a :
U V A U V A
A U V
(3,20)
A U V
1
U U Ut 2 ª 1 º½ J 1
Alors: ® «1 J 1 M ²»¾ (3-21)
U Ut U ¯J 1 ¬ 2 ¼¿
D’où
1 1 1 1
V § JRT · 2
JRT 2 §T · 2 § Tt · 2
¨¨ ¸¸ ¨¨ ¸¸ ¨ ¸
V © 1 ¹ V © Tt ¹ ©T ¹
(3-22)
1
V 1 2 ª J 1 º½ 2
® «1 M ²»¾
V M ¯J 1 ¬ 2 ¼¿
Finalement on obtient en combinant (3-10) (3-11).
1 J 1
ª 1 º 2 J 1
A
1
1 « 2
J 1 M ² »
« 1 » (3-23)
A M « J 1 »
¬ 2 ¼

Pour l'air γ=1.4
A 1 1 0.2M ²
3 (3-24)
A M 1.728
A
A
1.0
Nombre de
1.0 1.5 Mach
0 0.5 2.0
Figure 8.
Cette figure montre que la section minimale dans un écoulement isentropique est la
section du col. Toutes les autres sections A ont une surface plus grande que A .

PS
P Ps
Pi Pi
1
Ecoulement subsonique
1 dans toute la tuyère
1
PD
0.8 2 Pi
Onde de choc stationnaire
3
0.6
4
0.527
PB
0.4
Pi Décollement de veine dans le
divergent
0.2 PA
Ecoulement supersonique
Pi dans toute la tuyère
0 0
x
C
Cc
4
2 3
Figure(5)-Classification des écoulements dans

Une tuyère convergente-divergente

DEPARTEMENT GENIE MECANIQUE Module : Gazodynamique

Série N°3
Ex. 01 :
De l’air circule dans une tuyère convergente avec les propriétés suivante :
A l’entrée : M1=0.3, A1=0.001m2, P1=650KPa, T1=62°C.
A la sortie : M2=0.8
Calculer Tt1, T2, Pt2, P2, U2, V2 et A2.
Ex. 02 :
De l’air entre dans une tuyère convergente à une vitesse négligeable. Les
propriétés dans le réservoir sont : Tt=95°C et Pt =406kPa. La pression dans l’enceinte
est de 240KPa. La section au col est de 0.01m2.
Calculer le débit de l’écoulement.
Ex. 03 :
De l’hélium qui se détend isentropiquement à travers une tuyère convergente.
Calculer la densité à la sortie de la tuyère sachant que : Pb=Patm, Tt=250°C, Pt=800Kpa,
Ac=0.002m2, J=1.66, r =2.077 KJ/kg K.
Ex. 04 :
De l’hélium circule dans une tuyère convergente, calculer le nombre de Mach au
x
col et la section au col sachant que : T1=20°C, P1=350Kpa, V1=201m/s, m C =0.15kg/s.
Ex. 05 :
On considère l’écoulement isentropique de l’air à travers une tuyère convergente
divergente. On donne :
A l’entrée : T1=27°C, P1=600Kpa, A1=1250mm2, M1=0.5.
A la sortie : M2=2.
Calculer la section A2.
Ex. 06 :
L’air à la température de 97°C se détend dans une tuyère de la pression 70N/cm2
à la pression 10 N/cm2. La détente dans la tuyère est supposée isentropique dans le
convergent et avec une perte de 5% dans le divergent. Connaissant le rapport de détente
critique (rapport dépression) Dc=1.89 et la section au col Sc=2.30cm2. Calculer :
a) Le débit de l’air
b) La vitesse et la température de l’air dans la section de sortie de la tuyère.
La vitesse à l’entrée de la tuyère est négligeable.
Ex. 07 :
On donne :
Pb=954 kPa
Pt=1.0MPa AS=10-3 m2
Tt=350°K
Mc=0.68
Calculer : Tc, Pc, Vc , Mc et Ac.

Ex.08 :
Un gaz avec une masse molaire de 39.9 et un rapport de
chaleur spécifique 1.67 est déchargé d’un grand réservoir dans
lequel la pression est 500Kpa et la température est 300C à travers
une tuyère. Supposant que l’écoulement est isentropique
unidimensionnel, trouvez :
1) Si la pression à une section de la tuyère est 80Kpa, le
nombre de Mach, la température et la vitesse à cette
section.
2) Si la tuyère à une section circulaire de diamètre 12mm au point précédent, le débit
massique à travers la tuyère.
Ex.09 :
Un gaz avec une masse molaire de 4 et un rapport de
chaleur spécifique de 1.3 s’écoule à travers une conduite de
section variable. A un point de la conduite, la vitesse est
150m/s, la pression est 100KPa et la température est 15°C.
Trouvez le nombre de Mach à ce point. A un autre point de la
conduite la température est -10°C. Trouvez le nombre de Mach,
la pression et la vitesse à ce second point en supposant 2
1
l’écoulement unidimensionnel et isentropique.
Ex.10 :
L’air circule à travers une tuyère qui a une
section d’entrée d’aire 10cm2. Si l’air a une vitesse de
80m/s, une température de 28°C, et une pression de
700KPa à l’entrée et une pression de 250KPa à la
sortie, Trouver le débit massique à travers la tuyère.
En supposant l’écoulement isentropique
unidimensionnel, trouver la vitesse à la section de
1
sortie. 2
Ex.11 :
Un tube de Pitot est placé dans un écoulement d’air subsonique. La température et la
pression dans l’écoulement sont 27°C et 96KPa respectivement. La différence entre la pression
statique et la pression de Pitot est 32KPa. Trouver la vitesse d’air (a) supposant l’écoulement
incompressible (b) supposant l’écoulement compressible.
Ex.12 :
L’air circule dans une tuyère convergente divergente avec une section d’entrée d’aire
5cm et une section de sortie d’aire 3.8cm2. A la section d’entrée la vitesse de l’air est 100m/s,
2
la pression est 680KPa, et la température est 60°C. Trouver le débit massique à travers la tuyère.
En supposant l’écoulement isentropique, trouver la vitesse et la pression à la section de sortie.

Ex.13 :
Les gaz d’échappement d’une fusée ont une masse molaire de 14. Ils sont supposés se
comporter comme gaz parfait avec un rapport de chaleur spécifique de1.25.
Ces gaz sont accélérés à travers une tuyère. A un point de la tuyère où l’aire de la
section est 0.7m2, la pression est 1000KPa, la température est 500°C et la vitesse est 100m/s.
Trouver le débit massique à travers la tuyère, et la pression et la température d’arrêts. Aussi,
trouver la vitesse maximale qui sera générée en détendant l’écoulement. Si la pression à un
autre point est 100KPa, trouver la température et la vitesse en ce point de l’écoulement en
supposant l’écoulement unidimensionnel isentropique.
Ex.14 :
Un gaz a une masse molaire 44 et un rapport de chaleur spécifique 1.3. A un certain
point d’écoulement, la température et la pression statique sont 15°C et 80KPa, respectivement et
la vitesse est 100m/s. Le gaz est isentropiquement détendu jusqu’à que la vitesse atteint la
valeur 300m/s. Trouver la pression, température et le nombre de Mach à la section de sortie.
Ex.15 :
L’hélium à une pression de 120KPa et une température de 20°C s’écoule avec une
vitesse de 800m/s. Trouver le nombre de Mach, la température d’arrêt et la pression d’arrêt.
Ex.16 :
Considérons l’écoulement unidimensionnel à travers une conduite. A un certain point de
la conduite, la vitesse est 360m/s, la température est 45°C, et la pression est 120KPa. Trouver le
nombre de Mach, la température d’arrêt et la pression d’arrêt à ce point d’écoulement. Si la
température à un autre point de l’écoulement est 90°C, trouver le nombre de Mach et la pression
dans ce point.
Ex.17 :
Une tuyère convergente de forme conique a une section d’entrée de diamètre 15cm. La
pression, la température et la vitesse à l’entrée de la tuyère sont 70KPa, 60°C et 180m/s
respectivement. Si la pression à la sortie de la tuyère est 78KPa, trouver le diamètre de la
section de sortie.
Ex.18 :
Un avion vole au niveau de la mer avec une vitesse de 220m/s. Quelle est la pression
maximale qui peut agir sur la surface de l’avion.
Ex.19 :
Considérons l’écoulement d’air avec une vitesse de 650m/s, une pression de 100KPa, et
une température de 20°C. Trouver la pression et la température d’arrêt.

Chapitre IV Ecoulement adiabatique avec frottement
Chapitre 4.
Ecoulement monodimensionnel adiabatique

Avec frottement dans une conduite à section constante
(Théorème de FANNO)
x Comprendre L’Effet des frottements de la paroi de la conduite

sur les propriétés de l’écoulement
x Comprendre Le théorème de FANNO
x Connaitre les critères de FANNO ;
x Savoir calculer la longueur maximale
x Savoir tracer le diagramme T-S

CHAPITRE 4: Ecoulement Monodimensionnel Adiabatique

Avec Frottement dans une Conduite à Section Constante
(Théorème de FANNO)
4.1- Equations de base:
a) Hypothèses:
1-Ecoulement monodimensionnel et permanent.
2-Adiabatique ( GQ 0 ).
3-Pas de travail technique (W tech =0)
4-Energie potentielle négligeable (gdz=0)
5- section constante (A 1 =A 2 =A)
b) Equations de base :
1- Equation de continuité:
UVA cte
m
d’où
U1V1 U 2V2 G (G: vitesse massique ) (4-1)
2- Equation d'énergie:
h t1 W Q ht 2 o ht1 ht 2 o Tt1 Tt 2 (pour le gaz parfait)
d’où
1 1
h 1 V1 ² h2 V2 ²
2 2
c'est-à-dire
1 1 G²
h+ V ² cte o h cte (4-2)
2 2 U²
3- Théorème d'Euler:
¦ Fext m V2 V1
On a:
F f P1 A1 P2 A2 m V2 V1
-F f A(P 1 -P 2 ) ( U V A)1, 2 (V2 V 1 )
D’où
Ff
( P1 1 P2 ) ( U 2V2 ² U1V1 ²)
A
Ft
P 1 U1 V 1 ² P2 U 2V2 ² (4-3)
A
4- Deuxième loi de la thermodynamique:
'S S 2 S1 ! 0 (fonction irréversible) (4-4)
D’où :
'S R log( Pt 2 / Pt1 ) ! 0 Pt1 ! Pt 2

Càd que le frottement réduit la pression.
5- Fluide compressible et gaz parfait:

P
RT (4-5)
U
h c pT (4-6)
Finalement, on a obtenu un système de six équations à sept inconnues (P, V, T,

h, S, Fx, ρ).
Pour le résoudre, on doit chercher une 7 ème équation en fixant un paramètre. Pour
cela, on utilisera la courbe de FANNO.
Exemple:
De l'air circule dans une conduite isolée wQ 0 . On donne :
P 1 =98,75 k Pa, T 2 =287°K, D=7,16 mm, P t1 101 KPa , T t1 296 q K
x
Trouver : m , Pt 2 , F t ?
Solution:
On a d'après la loi de conservation de masse :
. . P1
m U1 V1 A1 U 2 V2 A2 o m V1 A 1
RT1
3D ² 3,14 (7,16 10 3 )²
ou A= 40,24mm 4,024 10 5 m²
4 4
P1 98,75 T1
On a 0,9777 table

o M 1 0,18 0,9936
Pt1 101 Tt1
T1 0,9936 296 294,10 K
P1
D’où : U1 1,17 Kg / m 3
RT1
V1 M 1 JRT1 0,18 1,4 287 294,10 65,3m / s
Finalement, on obtient :
.
m U1 V 1 A=1,17 65,3 4,024 10 5 =3,075 10 3 Kg / s
D'autre part, on a: T t1 T t 2 (conduite isolée)
T2 287 P2
D’où 0,9696 o M 2 # 0,4 0,8956
Tt 2 296 Pt 2
Or : V 2 M 2 JRT2 0,4 1,4 287 287 135 ,83m / s
U1 V1 1,17 65,3

U2 0,5625Kg / m3
V2 135,83

D’où P2 U 2 RT2
0.5625 287 287 46.33KPa
46.33
Finalement, on aura : Pt 2 51.73KPa
0.8956
L'équation d'équilibre des forces donne:
P1 U1V1 ² t P2 U 2V2 ² F ft A> P1 P2 U1V1 ² U 2V2 ² @
F
A
>
F ft 4.024 10 5 98.75 46.33 10 3 1.17 65.3² 0.5625 135.83² @
F ft 1.89 N
4-2- Ecoulement dans une conduite à section constante:
On considère ici l'effet des frottements des parois dans un écoulement

adiabatique à travers une conduite à section constante. Ce type d'écoulement est connu
sous le nom "d'écoulement de FANNO".
Considérons un élément de volume de l'écoulement et examinons l'équilibre de
forces agissantes: Volume de contrôle
x Force de pression : PA P dP A
V V+dV x
x Force de frottement : W P AP P P+dP
dx
Avec : Figure 1- volume de control utilisé

dans l'analyse d'un écoulement avec
AP périmètre.dx = P dx frottement dans une conduite
WP contrainte de cisaillement de la paroi due aux frottements.
D’où;
PA P dP A W P AP V dV V
m (4-5)
D’autre part, le diamètre hydraulique est définie en général par:

4 Air 4A
DH (4-6)
périmètre P
Dans le cas d'une conduite à section circulaire:
3D ²
4
DH 4 D
3D

D’où;
dP A W P P dx UVAdV
En divisant par UAV ² , on obtient:
dP W P dV
P dx (4-7)
UV ² UV ² A V
Or le principe de conservation de masse implique:
UVA cte UV cte
D’où;
UV U dU V dV UV VdU UdV dUdV
En négligeant les quantités dρdV, on obtient:
UdV VdU 0 (4-8)
En divisant par UV , on aura:
dV dU
0 (4-9)
V U
L'équation de l'énergie donne:
1 V² V dV ²
ht h V² cte c p T c p T dT cte
2 2 2
En négligeant dV ² puisque dV est très petite, on obtient:
c p dT VdV 0 (4-10)
En outre, les équations d'état d'un gaz parfait donnent:

P
RT P URT
U
P dP
Et R T dT P dP U dU R T dT
U dU
D’où :
URT dP URT RdUT RUdT RdUdT
En négligeant dρdT, on obtient:
dP RTdU URdT (4-11)
V² V²
D'autre part, M²
a² JRT
V dV ² V dV ²
Et M dM ²
a da ² T dT JR

2
§ dV ·
2 V ² ¨1 ¸
§ dM · © V ¹
Ou encore; M ² ¨1 ¸
© M ¹ § dT ·
JRT ¨1 ¸
© T ¹
En négligeant les dérivées d'ordre supérieur à 1, on aura:
§ dM · V ² § dV ·§ dT ·
M ² ¨1 2 ¸ ¨1 2 ¸¨1 ¸
© M ¹ JRT © V ¹© T ¹
dM § dV ·§ dT ·
1 2 ¨1 2 ¸¨1 ¸
M © V ¹© T ¹
dM dV dT dV dT
1 2 1 2 2
M V T V T
dVdT | 0
dM dV dT
D' où ; 2 2
M V T
Ou
dM dV dT
(4-12)
M V 2T
L'équation (4-6) implique en divisant par c p T :
dT VdV
0 (4-13)
T c pT
2
V² 2 V dV
Or puisque M ² et M dM ,Il vient:
JRT JR T dT
2VdV V²
2MdM dT
JRT JRT ²
Ou encore:
V ² dT
VdV JRTMdM
2 T
D’où :
dT JR V ² dT
MdM 0
T cp 2c p T T
En exprimant c p et R en fonction de γ, on aura:

dT J 1 dT
J 1 MdM M² 0
T 2 T
On obtient finalement:
dT J 1 M ² dM
(4-14)
T >1 J 1 M ² / 2@ M
Considérons maintenant l'équation (4-11) en divisant par P URT :
dP dU dT

P U T

dU
En éliminant entre cette équation et l'équation (4-9), on obtient:
U
dU dV dT
(4-15)
U V T
D'autre part, on a d'après (4-8):
dV dM dT

V M 2T
D’où;
dP dM dT

P M T
En substituant (4-14) dans cette expression, on aura:
dp ª J 1 M ² / 2 º dM
«1 » (4-16)
p ¬ 1 J 1 M ² / 2¼ M
En fin, résonnons à l'équation (4-7);
dV 1 dP W P P
dx 0 (4-17)
V JM ² P UV ² A
En substituant (4-12) dans (4-16), on obtient;
dM 1 dT 1 dP W P P
dx 0
M 2 T JM ² P UV ² A
En substituant (4-15) et (4-16) dans cette équation:
dM J 1 M ² / 2 dM 1 J 1 M ² dM W P
P dx 0
M >1 J 1 M ² / 2@ M >1 J 1 M ² / 2@ M UV ² A
Ou encore:
dM JM ²>1 J 1 / 2M ² @ ª W P P º
« UV ² A dx» (4-18)
M 1 M ² ¬ ¼
D’où en substituant dans (4-14) et (4-16):
dT J J 1 M 4 ª W P P º
(4-19)
1 M ² «¬ UV ² A »¼
dx
T
dP JM ²>1 J 1 M ² @ ª W P P º
« UV ² A dx» (4-20)
P 1 M ² ¬ ¼
Remarque:
Puisque la contrainte de cisaillement W p , la vitesse et le nombre de Mach sont
toujours positifs, l'équation (4-18) indique que le signe de "dM" dépend du signe de la
quantité (1-M²):

Si M<1, les frottements implique que le nombre de Mach augmente. Cependant

si M>1, le nombre de Mach décroît en présence des frottements. La viscosité
pousse toujours le nombre de Mach à tendre vers 1.
La conduite joue le rôle d'une tuyère pour M<1 et le rôle d'un diffuseur pour
M>1.
Puisque une fois le nombre de Mach atteint la valeur 1, les changements dans les
conditions en aval de l'écoulement n'affectent pas l'écoulement en amant.
Il pourra y avoir une onde de choc comme conséquence des frottements.
Les équations (4-19) et (4-20) montrent que si M<1, dT et dP sont négatives
tandis que si M>1, elles sont positives.
Examinons, enfin, l'effet des frottements sur l'entropie. On a :

dT dP
ds cp R
T P
Ou encore:
ds dT J 1 dP

cp T J P
En utilisant les équations (4-15) et (4-16), on obtient:
ds ª Wp P º
J 1 M ²« dx» (4-21)
cp ¬ UV ² A ¼
Cette relation montre que l'entropie augmente toujours:
dM dV dP dT ds
M<1 ↑ ↑ ↓ ↓ ↑
M>1 ↓ ↓ ↑ ↑ ↑
On peut exprimer toute les équations précédentes en fonction de facteur de

frottement de FANNO « f » exprimé par la relation :
Wp
f (4-22)
1
UV ²
2

On peut aussi utiliser un autre facteur dit facteur de frottement de Darcy, f D ,

donné par la relation:
fD
f (4-12bis)
4
En général:
§ E ·
f f ¨¨ Re, , M ¸¸
© DH ¹
16
x Pour un écoulement laminaire: f , Re nombre de Reynolds
Re
2
§ H 5.74 ·
x Pour un écoulement turbulent: f 0.0625 / log¨¨ 0.9 ¸¸
© 3.7 DH Re ¹
La transmission du laminaire vers le turbulent se fait à Re=2300.
La rugosité ε dépend du type de matériau de la conduite. On prend pour des tubes en
acier de commerce ε ≈ 0.045 mm et pour des tubes étirés à froid ε ≈ 0.0015mm.
Le facteur de frottement de Darcy peut être donné par le diagramme de Mody.
En fonction de f, l'équation (4-14) exprimé sous la forme:
4 fdx 2 1 M ² dM
(4-23)
DH ª 1 º M
JM ² «1 J 1 M ² »
¬ 2 ¼
4A
Avec : DH
P
L'intégration de (4-19) donne:
M2
2 1 M ²
l
4 fdx dM
³0 DH ³ § 1 · M
M1
JM ²¨1 J 1 M ² ¸
© 2 ¹
D’où:
M2
4f 2 1 M ² dM
DH
l ³ ª 1 º M
(4-24)
M1
JM ² «1 J 1 M ² »
¬ 2 ¼
Où f "facteur de frottement moyen" sur toute la longueur "l" de la conduite.
En général, le changement de f est insignifiant; ce qui pousse à utiliser f au lieu de f .

ª 1 º
M 21 «1 J 1 M 2 2 »
4f 1§ 1 1 · J 1 ¬ 2 ¼
l ¨ 2 2 ¸ log (4-25)
DH J © M 1 M 2 ¹ 2J ª 1 º
M 2 2 «1 J 1 M 21 »
¬ 2 ¼
Puisqu'en frottements importants, le nombre de Mach tend vers 1, alors on
prend M 2 1 et la longueur l l max l :
4f
l max
§1 M ² · J 1
¨¨ ¸ log
J 1 M 2 > @ (4-26)
2 ¸
DH © JM ¹ 2J ª 1 º
2«1 J 1 M 2 »
¬ 2 ¼
De même, on obtient à partir de (4-15), (4-16):
ª º
P
dP
1 « J 1 M ² / 2 »» dM
³ P ³M ««1 ª 1 º» M
«¬1 2 J 1 M ² »¼ »
P
«
¬ ¼
Ou encore;
1
ª º 2
P 1 « J 1 / 2 »
« » (4-27)
P M «ª 1 º»
« «¬1 2 J 1 M ² / 2»¼ »
¬ ¼
T
dT
1
J 1 M
³T T ³
M
1 J 1 M ² / 2
dM
Ou encore:
T J 1 / 2
(4-28)
T 1 J 1 M ² / 2
Où P et T des valeurs de P et T à M=1.

Pour des pressions et températures totales, on obtient:
J 1
P0 1 ª1 J 1 M ² / 2 º 2 J 1
(4-29)
P0 M «¬ J 1 / 2 »¼
Puisque :
J
P0 ª J 1 º J 1
«¬1 M ² »¼
P 2
En conclusion, les équations précédentes (4-2bis),(4-27),(4-28) sont de la forme:

4f
l max f ( M )
D
P
f (M )
P
T
f (M )
T
Ces fonctions sont disponibles sous forme de graphes ou de tableau (tableau 2).
Pour dériver des formules pratiques, on reprend la formule (4-24) :
1
1 M ²
L
fdx
³0 DH ³ § 1 ·
dM ²
M JM
¨1 J 1 M ² ¸
2 4
© 2 ¹ l1
fL 1 M ² J 1 J 1 M ²
ln l2
DH JM ² 2J 2 J 1 M ² Figure 2-
f est le facteur de frottement moyen entre 0 et L . (L* correspond à M=1)

Plusieurs problèmes ont une petite longueur pour laquelle l’écoulement ne
devient pas sonique
'L § fL · § fL ·
f ¨¨ ¸¸ ¨¨ ¸¸ x Lmax
D © D ¹1 © D ¹2
Log f x
laminaire turbulent
1 1
ª J 1 º 2 ª J 1
1
º
M ²»
2
H /D
P 1 « 2 » , U V 1 « 2
« » «
M « J 1 U V M« J 1 »»
M ²»
P
1
¬« 2 ¼» ¬« 2 ¼» 64
J 1 U
ª J 1 º 2 J 1 Re
1 M ²»
Pt 1 « 2
« ,
Pt M« J 1 »» Log Re
¬« 2 ¼» Figure 3-
fLmax
DH
Toutes ces quantités sont obtenues par le tableau de FANNO.

Où f=f(Re , ε/D)
Diagramme de Mody
UVD VD
Re Force d'inertie/force de frottement
P X

UV ² UVD
Re
P
V P
D
μ= viscosité dynamique
υ= viscosité cinématique
dx x x max
Pour M=1, 0®
dM ¯x x min
Or: x=0 implique qu'il n'existe de conduite, ce qui veut dire que la solution sera pour
x xmax
Conclusion:
1. En subsonique, la conduite travaille comme tuyère et accélère l'écoulement
jusqu'à M=1 à x= Lmax .
2. En supersonique, la conduite travaille comme diffuseur et décélère l'écoulement
jusqu'à M=1 à x= Lmax .
Critère de FANNO:
Que ce passe t-il si la longueur de la conduite est plus grande ou plus petite que Lmax ?
FANNO d'écrit ce phénomène en utilisant le diagramme T-S.
On aura besoin des équations T=f(M) et s=f(M), on obtient :
J 1
ª 1 J 1 º
ss § 2 · 2J
§ · J § J 1 T · 2J
ln «¨¨ ¸ »
T
¸ ¨ ¸ ¨ (4.30)
cp «© J 1 ¸¹ © T ¹ © 2 T ¹ »
¬« ¼»

T M<1 (subsonique)
m
T0 G UV
A
M=1 M=1 M=1 x LMAX

T
(s , T , P , A ,V , U , M 1)
G1 > G2 > G3
M>1 (supersonique)
s
s1 s2 s3
Figure 4-
L'entropie S augmente en écoulement adiabatique.

L'entropie S augmente jusqu'à une limite critique à M=1, x=Lmax.
Les propriétés P, T, V, ρ croissent en écoulement subsonique et décroissent en
écoulement supersonique.
Les frottements font augmenter la vitesse et diminuer la pression P0 .

Figure 5- Digramme de Moody

Exemples :
Exemple 01 :
De l’air s’écoule dans une conduite de 5 cm de diamètre. L’air entre avec M=2.5
et sort à M=1.5. Quelle est la longueur de la conduite ? Quelle longueur prendra
la conduite si M=1 à la sortie ? on suppose que ݂ ҧ ൌ ͲǤͲͲʹ et que l’écoulement
est adiabatique.
Exemple 02 :
De l’air s’écoule à travers un tuyau de diamètre 0.3m avec un débit ‫ݍ‬ሶ ൌ
ͳͲͲͲm3/min à la température et la pression de 293K et 150KPa respectivement.
Si la longueur du tuyau est 50m et f=0.005, trouver le nombre de Mach à la
sortie du tuyau, la pression et la température d’entrée.
Exemple 03 :
De l’air s’écoule à travers un pipe en acier inoxydable de diamètre 5cm. L’air
entre à M=0.3 avec la température et la pression de 40°C et 150KPa
respectivement. En utilisant un facteur de frottement évalué à partir des
conditions d’entrée, trouver le nombre de Mach, la pression et la température
aux distances 4, 8, 12, 16, 18 et 19m de l’entrée.
En utilisant ces résultats, trouver les facteurs de frottement actuels en ces
positions dans le pipe et comparer ces valeurs à la valeur supposée constante.
Rappel : variation du facteur de frottement
࣋ࢂࡰ
f=fonction (Re, ε/DH) avec Re= ࣆ ࡴ
or ࢓ሶ ൌ ࣋ࢂ࡭
࢓ሶࡰ ૝࢓ሶ
d’où : ࡾࢋ ൌ ࡭ࣆࡴ ൌ ࡼࣆ (P= perimètre)
puisque ࢓ሶ ൌ ࢉ࢙࢚ࢋ et P=cste alors Re est proportionnel à ૚Ȁࣆ

le coefficient de viscosité ࣆ varie en fonction de la température :
ࣆ૚ ࢀ ࢔
ൌ ቀ ૚ቁ ( pour les gaz)
ࣆ૛ ࢀ૛
Où T est la température absolue

n index compris entre 0.5 et 0.8 pour les gaz.
Cependant le changement de Re le long d’une conduite est généralement petit.

En outre, puisque l’écoulement est souvent turbulent, « f » n’est pratiquement
dépendant que de Re. Donc, il est souvent adéquat de traiter « f » comme
constante et de l’évaluer en utilisant les conditions d’entrée.

Exemple 04 :
Un travailleur dans une mine souterraine est approvisionné en oxygène à un
débit de 0.06kg/s à travers une corde ombilicale qui a une longueur de 8m. La
pression et la température de l’oxygène à l’entrée de la corde sont maintenues à
250kPa et 10°C respectivement et la pression à la sortie est de 50kPa.
Trouver le diamètre de la corde ombilicale si elle est fabriquée d’une matière
donnant une rugosité effective ε de 0.005mm.
(on donne pour l’oxygène : R=260J/kg.°K, γ=1.4)


Série N°4
EX.01 :
De l’air s’écoule à travers une conduite à section constante. La pression, la température et le
nombre de Mach à l’entrée de la conduite sont respectivement 180KPa, 30°C et 0.25. Si le nombre de
Mach à la sortie atteint 0.75 à cause des frottements, déterminer la pression, la température et la vitesse de
sortie. Supposer que l’écoulement est adiabatique.
EX.02 :
L’air s’écoule à travers un tuyau bien isolé de 4in de diamètre avec un débit de 500lbm/min. la
pression baisse de 50 Psia à l’entrée jusqu’à la valeur de 40 Psia à la sortie. Si la température à l’entrée est
200°F, trouver le nombre de Mach à la sortie du tuyau.
Inch=2.54cm
PSI=6894.757 Pa
Lbm= 0.4535924 kg
EX.03 :
De l’air s’écoule à travers un tuyau de 5cm de diamètre. Les mesures indiquent que la vitesse à
l’entrée était 70m/s, la température est 80°C et que la pression est de 1Mpa. Trouver la température, la
pression et le nombre de Mach à la sortie du tuyau s’il a une longueur de 25m.
Supposer un écoulement adiabatique et que le facteur de frottement est 0.005.
EX.04 :
L’air traverse une conduite de 0.15m de diamètre. A son entrée, le nombre de Mach est égale à
0.1, la pression est 70KPa et la température est 35°C. si l’écoulement peut être considéré comme
adiabatique et que le facteur de frottement est 0.005, déterminer la longueur de la conduite pour que le
nombre de Mach sera égale à 0.6 à la sortie. Trouver la pression et la température à la sortie de la
conduite.
EX.05 :
De l’air s’écoule, à partir d’un large réservoir, à travers une conduite bien isolé de 12mm de
diamètre. Si l’air entre dans la conduite à M=0.2 et la quitte = M=0.6, trouver la longueur de la conduite.
Considérer que le facteur de frottement est de 0.005. Quelle doit être la longueur de la conduite si le
nombre de Mach à la sortie est 1.
EX.06 :
De l’air est stocké dans un réservoir à une pression de 1.6Mpa et une température de 20°C. Quel
est le débit massique maximum possible de l’écoulement (à partir du réservoir) à travers d’une conduite
de 1.2cm de diamètre et 30cm de longueur. La conduite décharge dans l’atmosphère où la pression

atmosphérique est de 101KPa. Le facteur de frottement moyen est de 0.006 et l’écoulement est supposé
adiabatique subsonique.
EX.07 :
L’air s’écoule à travers une conduite de 12m de long et 25mm de diamètre. A l’entrée de la
conduite, l’air a une vitesse de 80m/s, la pression est 350KPa et la température est 50°C. Si le facteur de
frottement est de 0.005, trouver la vitesse, la pression et la température à la sortie de la conduite.
Considérer l’écoulement adiabatique.
EX.08 :
De l’air se détend à partir d’un large réservoir dans lequel la pression et la température sont
respectivement 200KPa et 30°C, à travers une tuyère convergente qui donne à sa sortie un nombre de
Mach de 0.2. L’air s’écoule alors à travers une conduite de 25mm de diamètre, le nombre de Mach à la
sortie de la conduite est 0.8. En considérant que l’écoulement dans la tuyère est isentropique et qu’il est
adiabatique dans la conduite, trouver la longueur de la conduite et la pression à sa sortie. Le facteur de
frottement dans la conduite est supposé égale à 0.005.

Chapitre V ECOULEMENT DANS UNE CONDUITE AVEC TRANSFERT DE CHALEUR
Chapitre 5.
ECOULEMENT DANS UNE CONDUITE à SECTION CONSTANTE AVEC

ECHANGE DE CHALEUR (SANS FROTTEMENTS)
(Théorème de RAYLEIGH)
x Comprendre L’Effet des échanges de chaleur à travers la

paroi de la conduite sur les propriétés de l’écoulement
x Comprendre Le théorème de Rayleigh
x Savoir calculer les courbes de Rayleigh
x Savoir tracer le diagramme T-S (courbe de Rayleigh)

Chapitre 5:
ECOULEMENT DANS UNE CONDUITE à SECTION CONSTANTE AVEC

ECHANGE DE CHALEUR (SANS FROTTEMENTS)
5.1- INTRODUCTION :
On va considérer dans ce chapitre que l'écoulement n'est pas adiabatique,
mais il y’a un apport ou retrait de chaleur à travers les parois de la conduite à section
constante.
Exemple: chauffage on refroidissement des parois ou à partir des réactions chimiques
qui auront lieu dans l'écoulement comme est le cas de la combustion (chambre de
combustion).
On adoptera les hypothèses suivantes :
1 dQ 2 P dP
1- régime permanent
P
2- section constante (dA=0) V dV
V
. M dM
3- échange de chaleur ( q z 0) M
T dT
TP
4- absence de frottement (f=0) U dU
5- écoulement monodimensionnel Figure 5.1- volume de control
6- énergie potentielle négligeable ( ' Z=0)
1- Equation de continuité :
.
m
U1 V1 A1 U 2 V2 A2 U1 V1 U 2 V2 UV (5.1)
A
2- Théorème de quantité de mouvement :
.
U 2 A 2 U1 A1 m(V2 V1 ) P1 U1 V1 ² P2 U 2V2 ² (5.2)
3- Equation de l'énergie:
1 1
h 1 V1 ² gZ1 q Wu h2 V2 ² gZ 2 (5.3)
2 2
d’où :
V2 ² 1
q=h 2 h1 V1 ² q ht 2 ht1 c p (Tt 2 Tt1 ) (5.3b)
2 2

4- 2 ème loi de la thermodynamique :

dP dT dP
T dS=dh - dS cp R (gaz parfait) (5.4)
U T p
௉
Avec : ൌ ܴܶ
ఘ
5.2- Calcul et équations de base :

Considérons un écoulement d'un gaz à travers le volume de control de la figure
(5.1).
L'équation de conservation de masse donne:
U A V cte (5.5)
Puis que A=constante alors U V cste
d’où : U V ( U dU )(V dV )
U V U V U dV V dU VdU dU dV tel que d U dV 0
d’où U dV V dU 0 (5.5b)
En divisant par U V , on obtient:
dV dU
0
V U (5.5c)
Le théorème des mouvements appliqué au volume de control donne:

P A ( P dP) A U V A>(V dV ) V @ (5.6)
D’où:
dP P dP A UVA V dV V (5.6b)
D’où:
dP UVdV 0 (5.6c)
En fin, le principe de conservation d'énergie donne (étant donné qu'il n'y a pas de
travail produit) :
ª V dV ² V ²º
dQ m «c p T dT c pT » (5.7)
¬ 2 2¼
On note ici que l'échange de chaleur est, thermodynamiquement, strictement
le transfert d'énergie dû à la différence de températures.

Dans les systèmes de combustion, il est plus rigoureux d'écrire l'équation d'énergie telle
que;
ª V dV ² V ²º
m 'h f m «c p T dT c pT » (5.7b)
¬ 2 2¼
Où 'h f est la différence entre la chaleur de formation des réactants et la chaleur de
h f comme étant
formation des produits. Cependant, il est plus adéquat de considérer m'
la chaleur produite comme résultat de la combustion.

En outre, il est convenable de définir « dq » comme étant le rapport de la chaleur
échangée par unité de débit masse:
dQ
dq (5.8)
m
L'équation de conservation d'énergie (5-7) peut s'écrire donc sous la forme:
dq =c p dT V dV (5.9)
D'autre part, considérons la définition de la température d'arrêt T 0 : (de

stagnation)
V²
T0 T
2c p
(V dV )²
Ou T 0 dT0 (T dT )
2c p
V dV
D’où : dT 0 dT
cp
Alors, l'équation (5.9) peut prendre la forme :

dq=c p dT0 (5.9b)
L'équation (5.9) donne :

dq dT
cp V (5.10)
dV dV
On a pour un gaz parfait, on a:
P
RT P URT
U
P dP
Et R(T dT ) P dP ( U dU ) R(T dT )
U dU
D’où :

dP dU dT
(5.11)
P U T
Or d'après l'équation (5.5c), on a :
dU dV
(5.12)
U V
Et d'après (5-2), on a:
dP VdV
U (5.12b)
P P
Donc l'équation (5.11) peut être mise sous la forme :

VdV dV dT
-U (5.13)
P V T
D’où:
VdV dV dT

RT V T
dT T V
On encore (5.13b)
dV V R
En substituant dans (5.10), on obtient :
dq T cp
cp V ( 1) (5.14)
dV V R
R J 1
Or, puis que R=c p c v et , il vient que :
cp J
dq T V
cp (5.14b)
dV V J 1
T V
On remarque que si V n , le terme c p p par contre si V n , le terme n.
V J 1
A des très petites vitesses, le premier terme sera dominant et par conséquent puis
dq
que T et V sont positives, est donc positif c à d que (dV) aura le même signe
dV
que dq. Ceci signifie qu’à des petites vitesses, l'apport de chaleur diminue la vitesse.

Aux grandes vitesses, le second terme de l'équation (5.14b) sera dominant et

dq
sera négative. Donc, aux grandes vitesses, «dV» a un signe opposé à celui de
dV
« dq ».
Ce qui signifie que l'apport de chaleur fait décroître la vitesse par contre le retrait de
chaleur augmente la vitesse.
dq T V
La transition d'une forme à une autre se produira à 0 , d’où : cp
dV V J 1
c à d quand V= (J 1)c p T JRT a (5.15)
On conclut que :
Apport de chaleur Retrait de
chaleur
M<1 V↑ V↓
M>1 V↓ V↑
Tableau (5-1)
c à d quand M=1, on conclut que : (voir tableau (5.1))

dq
si M>1, !0
dV
dq
si M<1, O (5.16)
dV
dq
si M=1 0
dV
Introduisons maintenant la définition du nombre de Mach :
V=M a
V dV M dM a da Ma adM dM da Mda
Division par V et remplaçant par Ma=V
V dV Ma adM Mda dMda

V V V V V V
En négligeant dM.da, on obtient:
dV da dM
(5.17)
V a M

On a vue que:
VdV dV dT
U
P V T
Ou encore:
VdV dV dT § P·
J ¨¨ a² JRT J ¸
a² V T © U ¸¹
D’où:
dV dT
1 JM ² (5.18)
V T
Finalement, il vient que:
1 JM ² dV dM
(5.18b)
2 V M
On remarque que « dM » a le même signe que « dV ». Ce qui signifie que l'apport de
chaleur tendra M vers 1 par contre le retrait de chaleur le retrait de chaleur tend M à
s'éloigner de 1.
dT
L'équation (5.13b) montre aussi que est positive pour des faibles vitesses et
dV
dT T V
négative pour des grandes vitesses. sera nulle quand c à d quand:
dV V R
1
V²=RT M (5.19)
J
1
Ceci montre que la température maximale dans un écoulement existera quand M=
J
D'autre part, la relation d'entropie est donnée par :
dS dT J 1 dP
(5.20)
cp T J P
En utilisant (5-7) et (5-1), on obtient:

dS 1 dT (J 1) dV
(5.20b)
cp J T J V
En substituant (5.13b) dans cette équation, on obtient :

dS dV
(1 M ²) (5.20c)
cp V

Cette équation montre aussi que l'apport de chaleur augmente toujours l'entropie. Par
contre, le retrait de chaleur diminue toujours l'entropie. Le maximum d'entropie se
produit à M=1.
Apport de chaleur Retrait de chaleur

M<1 S↑ S↓
M>1 S↑ S↓
Tableau (5.2)- changement d'entropie convenable
Il serait aussi convenable d'exprimer les changements qui se produisent dans

l'écoulement en termes de nombre de Mach. Considérons deux points d'un écoulement.
1 2
dQ
P1 P2
V1 V2
M1 M2
Figure (5.2)- Volume de contrôle
L'équation d'équilibre des mouvements appliquée au volume de contrôle donne :

P1 U1V1 P2 U 2V2
2 2
(5.21)
JP
Puisque a ² , on peut écrire :
U
P1 1 JM 1 P2 1 JM 2
2 2
(5.21b)
D’où:
1 JM 1
2
P2
(5.22)
1 JM 2
2
P1
Par définition, la pression de stagnation peut être donc exprimée par :

J
P0 § J 1 · J 1
¨1 M ²¸ ; ( P0 o M 1) (5.23)
P © 2 ¹
Lorsqu'il y a retrait de chaleur, la pression de d'arrêt changera.

A partir des équations (5-22) et (5-23) vient:
> @
J
P2 1 J 1 M 2 / 2 J 1
2
P02
(5.23b)
> @
J
P01 P1
1 J 1 M 1 / 2 J 1
2
D'autre part, l'équation de continuité donne:

U2 V1
(5-24)
U1 V2
Et dans le cas d'un gaz parfait, on aura:
U2 P2 T1
(5-24b)
U1 P1 T2
D’où;
T2 P2 V2
(5-25)
T1 P1 V1
Or:
V2 M 2 a2 M2 T2
(5.25b)
V1 M 1 a1 M1 T1
En substituant (5-25b) et (5-22) dans (5-25), on obtient:
T2 >
M 2 1 JM 1
2
@ 2 2
(5.25c)
T1 M1
2
>1 JM @2
2 2
Et en substituant encore dans (5-23), on aura:

V2
2
>
M 2 1 JM 1 @ 2
V1 M1
2
>1 JM @2
2
(5-26)
Aussi, la substituant des équations (5-26) et (5-22) dans (5-24b) donne:

U2 >
M 1 1 JM 2
2
@ 2
U1 M2
2
>1 JM @1
2
(5-27)
De même façon, la définition de la température d'arrêt donne:

T0 J 1
1 M²
T 2
D’où:
T02 >
T2 1 J 1 M 2 / 2
2
@
T01 >
T1 1 J 1 M 1 2 / 2 @ (5.28)

En substituant (5-25c) da ns cette équation, on obtient:
T02 M 2 1 JM 1
2 2 2
>1 J 1 M 2
/2 @
>1 J 1 M / 2@
2
2
(5.28b)
2 2 2
T01 M 1 1 JM 2 1
Remarque:
Il est important de noter que l'équation (5.10) montrer que « T 0 » change
dans le même sens que « q » c à d que l'apport de chaleur augmente T 0 et le retrait
de chaleur de diminue T 0 .
D'autre part, on a vue que l'apport de chaleur fait varier M vers 1. Par
contre le retrait de chaleur éloigne M de 1. Alors, T0 doit avoir une valeur maximale
quand M=1 voir figure (5.3).
Pour une valeur quelconque de M 1 et pour M 2 =1, on a :
TO max 1 (1 JM 1 ²)²
(5-28c)
TO1 2(1 J ) M 1 ² >1 (J 1) M 1 ² / 2@
T0
retrait de chaleur
T0max
retrait de chaleur
Apport de chaleur Apport de chaleur
M=1 M
Figure (5.2)- variation de la température d’arrêt

en fonction du nombre de Mach

De même façon, la relation entre l'entropie et le nombre de Mach est donnée par :
§T · §P ·
S 2 S1 c p Ln¨¨ 2 ¸¸ R.Ln¨¨ 2 ¸¸ (5.29)
© T1 ¹ © P1 ¹
Ou encore :
S 2 S1 T2 (J 1) P
Ln( ) Ln( 2 )
cp T1 J P1
(5.29b)
S 2 S1 ª T P J 1 º
Ln « 2 ( 1 ) J »
cp ¬ T1 P2 ¼ (5.29c)
En utilisant les équations (5-25b) et (5-20c) , il vient que :

S 2 S1 ª M ² (1 JM 1 ²)² 1 JM 1 ² (J 1) J º
Ln « 2 ( ) »
cp ¬ M 1 ² (1 JM 1 ²)² 1 JM 1 ² ¼ (5.29d)
J 1
ª º
S 2 S1 M ² § 1 JM ² · J
D’où : Ln « 2 ¨¨ 1
¸¸ » (5.30)
cp « M 1 ² © 1 JM 2 ² ¹ »
¬ ¼
Rappelons l'équation (5.9b) : Q=c p (T02 T01 )
q TO 2
D’où 1 (5-31)
c p T01 T01
En utilisant le résultat (5.6), on obtient alors :

q M 2 ² >1 JM 1 ² @ >1 (J 1) M 2 ² / 2@
1 (5.31b)
c p T01 M 1 ² >1 JM 2 ² @ >1 (J 1) M 1 ² / 2@
Remarque:
q T02
On voit que puisque est fonction uniquement du rapport de
c p T01 T01
l'équation (5.31). Il est souvent plus pratique d'utiliser le rapport de température

d'arrêt pour mesurer l'échange de chaleur.
Finalement, il est convenable d'écrire les équations ci-dessus en terme de
conditions qui existent quand M 2 1 ; dans ce cas : P 2 P , T 2 T etc… On aura
alors les équations suivantes:

T0 2(J 1) M ²>1 (J 1) M ² / 2@
(5-32)
T0 (1 JM ²) 2

T (1 J )² M ²
(5.32)
T (1 JM ²)²
P 1 J
(5.33)
P (1 JM ²)²
J
P0 (1 J ) 2 ª (J 1) º½ J 1
®( ) «1 M ²»¾ (5.34)
P0 (1 JM ²) ¯ J 1 ¬ 2 ¼¿
V (1 J ) M ²
(5.35)
V (1 JM ²)
U (1 JM ²)
(5.36)
U (1 J ) M ²
SS ª 1 J J 1J º
lin« M ²( ) » (5.37)
cp ¬ 1 JM ² ¼
T0 T P
Les valeurs des équations des quantités , , … etc sont calculées pour
T0 T P
n'importe quelle valeur de M. Elles sont généralement données sous forme de tables
(table A12) "table de Rayleigh".
Remarque :
On note ici qu'aucune hypothèse n'est faite pour donner la manière dont
l'échange de chaleur se fait entre les sections considérées de l'écoulement.
L'échange se fait à travers une petite partie de l'écoulement uniquement ou à
travers toute la région entre les deux sections considérées de l'écoulement mais d'une
façon uniforme.
5.3- Relations entropie- température / courbe de Rayleigh :

On a déjà vu la relation :
dS dT J 1 dP

cp T J P
D’où:
S S § T · J 1 § P ·
ln ¨ ¸ ln ¨ ¸ (5.38)
cp ©T ¹ J ©P ¹

Or d'après la relation (5.33), on a:

1 J ª P º 1
M²
J «¬ P »¼ J
(5.39)
Et à partir des relations (5-32) et (5-33), on obtient :
0.5
P 1 §T ·
¨ ¸ (5.40)
P M ©T ¹
D’où:
2
§ P · ª P º T
¨ ¸ « 1 J 1» J (5.40b)
©P ¹ ¬ P ¼ T
Ou encore:
2
§ P · P T
¨ ¸ 1 J J 0 (5.40c)
©P ¹ P T
§ P ·
En résolvant cette équation quadratique de ¨ ¸ , on aura :
©P ¹
1 J 4J T
2
1 J
2
P T
r (5.40d)
P 2 2
Finalement, en substituant (5.40c) dans (5.40b), on obtient:
1 J
2
4J T ½
§ T · ª J 1º ° 1 J
2
ss T °
ln ¨ ¸ « » ln ® r ¾ (5.41)
cp ©T ¹ ¬ J ¼ ° 2 2 °
¯ ¿
La courbe représentative de l'équation (5.41) est connue sous le nom de "ligne
T
de Rayleigh" c'est le variation de en fin de s s / c p pour une valeur quelque de
T
J (Figure (5.3)).

T
T
Tmax B
T
1 A
Apport
M<1
Retrait M>1
Apport de chaleur
Retrait de chaleur s s / cp
Figure (5.3)- ligne de Rayleigh
L'écoulement en question, c à d unidimensionnel avec échange de chaleur sans

frottement est comme sous le nom "d'écoulement de Rayleigh".
>
Il est clair de la Figure (5.3) qu'à chaque valeur de ( S S ) / c p @ correspond deux
T
valeurs possibles de ; l'une correspond à l'écoulement subsonique et l'autre à
T
l'écoulement supersonique.
Le point A correspond à un nombre de Mach M=1 et T=T* .
Le point B est le point pour lequel T est atteint sur maximum.
dT
Pour obtenir ceci (maximum de T ) faisons la dérivée et égalons à Zéro, on
dM
obtient :
2M (1 J )² 4JM 3 (1 J )²
0 (5.42)
(1 JM ²)² (1 JM ²) 3
1
D’où M (5.43)
J
C'est le nombre de Mach auquel T=T max
En substituant dans (5-32), on obtient :
Tmax (1 J )²
(5.44)
T 4J

L'influence de l'apport de chaleur sur les différentes variables dans un

écoulement de Rayleigh est représentée par le tableau (5.3):
q T0 M P V S
M<1 n n n p n n
M>1 n n p n p n
M<1 p p p n p p
M>1 p p n p n p
Tableau (5-3)


Série N°5
Exercice 1-
De l’air s’écoulant à travers une conduite à section constante s’échauffe. La température
de l’air à une certaine section de la conduite est de 200°C.
Si au-delà d’une section voisine la température de stagnation augmente de 1%, estimer le
pourcentage d’augmentation de V et de M et la valeur de (dS/Cp) pour des nombres de Mach,
aux sections considérées, de 0.4, 0.8, 1.2 et 1.6.
Exercice 2-
De l’air s’écoule sans frottement entre deux sections (1) et (2). On donne : T1=333 K,
P1=135KPa, V1=732 m/s et M2=1.2.
Calculer : a) P2, T2, ρ2 , b) q12 , c) S2-S1
Exercice 3-
De l’air s’écoule à travers une conduite à section constante. La pression et la
température de l’air à l’entrée de la conduite sont respectivement 100KPa et 10°C et le nombre
de Mach à l’entrée est 2.8. De la chaleur est transférée à l’air pendant qu’il d’écoule à travers la
conduite et en sort avec un nombre de mach de 1.3. Trouver la pression et la température à la
sortie de la conduite.
S’il n’y a pas d’ondes de choc, trouver la quantité maximale de chaleur qui puisse être
transférée à l’air par unité de masse d’air. Trouver aussi la pression et la température qui peut
exister avec ce taux max de transfert de chaleur.
Supposer un écoulement permanent et que l’effet des frottements est négligeable et que l’air est
un gaz parfait.
Exercice 4-
De l’air s’écoule à travers une conduite de section constante dont les parois sont
maintenues à basses des températures. L’air entre avec un Mach de 0.52, une pression de
200KPa et une température de 350°C. Le taux de transfert de chaleur de l’air aux parois de la
conduite est estimé à 400KJ/Kg d’air.
Trouver le nombre Mach, la température et la pression à la sortie de la conduite. Supposer que
l’écoulement est permanent, les frottements sont négligeables et que l’air est un gaz parfait.
Exercice 5-
L’air entre dans une chambre de combustion à une vitesse de 80m/s sous une pression
de 180KPa et une température de 120°C. Trouver la quantité maximale de chaleur qui puisse
être générée dans la chambre de combustion par unité de masse d’air.

Si le carburant a une chaleur spécifique de 45 MJ/Kg, trouver le taux air-carburant. Si ce rapport

air/carburant est ajusté jusqu’à 90% de cette valeur, trouver la réduction dans le débit massique
de l’écoulement à travers la chambre de combustion qui doit avoir si la température et pression
d’arrêt restent les mêmes. Supposer que l’écoulement est permanent, les frottements sont
négligeables, que l’air est un gaz parfait et que l’effet de la masse du carburant est négligeable.

Chapitre VI Ondes de Choc
Chapitre 6.
LES ONDES DE CHOC
x Comprendre L’Effet des irréversibilités sur les propriétés de

l’écoulement
x Comprendre Le théorème de Rayleigh-Prandtl
x Distinguer entre les types d’ondes de choc (droite, oblique et
courbe).
x Distinguer entre Onde de Choc Statique ou en Mouvement.

Polycope Gazodynamique Alliche2021 PART1

Transféré par

Informations du documentcliquez pour développer les informations du document

Informations du documentcliquez pour développer les informations du document

Droits d'auteur :

Formats disponibles

Polycope Gazodynamique Alliche2021 PART1

Transféré par

Informations du document

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

Polycope Gazodynamique Alliche2021 PART1

Transféré par

Droits d'auteur :

Formats disponibles

REPUBLIQUE ALGERIENNE DEMOCRATIQUE ET POPULAIRE

MINISTERE DE L’ENSEIGNEMENT SUPERIEUR ET DE LA RECHERCHE SCIENTIFIQUE

COURS DE GENIE MECANIQUE

Auteur : Dr ALLICHE Mounir

Chapitre 1 : Introduction à la Gazodynamique (1 semaine)

Chapitre 2 : Ecoulement Isentropique 1D en Conduit à Section Variable (5 semaines)

Chapitre 3 : Ondes de Choc (4 semaines)

3. Coefficient du frottement et variation d’entropie

Chapitre I: Introduction Générale

x Des propriétés des gaz parfaits

1.1 Notions de base :

1-2-Propriétés des gaz :

On peut classer les propriétés des gaz en trois catégories:

1-2-1- Propriétés de diffusion:

Ce sont des propriétés qui indiquent la diffusion de quelques états déterminés

A:Viscosité dynamique(μ) ou viscosité absolue:

C'est l'indice de mesure de la diffusion de la quantité de mouvement par la

B:Viscosité cinématique (ν):

C:le coefficient de la conduction thermique spécifique(k):

1-2-2- Propriétés physiques:

1-2-3- Propriétés thermodynamiques:

A:l'énergie interne spécifique (e) ou (u):

B:le travail d'écoulent spécifique:

C:l'enthalpie (h):(l'énergie emmagasinée spécifique)

C'est la somme du travail d'écoulement et de l'énergie interne par unité de masse

D:la chaleur spécifique:

C'est la moyenne de la variation de l'énergie du fluide spécifique par rapport à la

Où : q=quantité de chaleur transportée du ou à l'unité de masse C [KJ/Kg.°C]

Selon le second principe de la thermodynamique, il existe une variable

1-3-Application du premier principe de la thermodynamique:

D'après le 1er Principe relatif à la conservation d'énergie, on peut considérer

- Travail mécanique entre l'entée et la sortie; dWm

- Quantité de chaleur reçue on dégagée par unité de masse; Q12 dm

-Variation de l'énergie interne; dm u 2  u1

Ainsi, le 1erPrincipe appliqué à ce système s'écrira;

b) On notera que le 1er Principe de la thermodynamique se bornant à constater la

1-4- Application aux écoulements sans échanges d'énergie;

Cas d'un fluide incompressible:

1-5-Application du second principe de la thermodynamique:

Selon le second principe de la thermodynamique, il existe une variable d'état S,

1-6-Rappels de physique du son:

L'onde acoustique, en traversant le fluide à l'arrêt(après déplacement du piston), va

- Equation de quantité de mouvements'écrit;

spécialement une distinction entre les différents intervalles du nombre de Mach.Ainsi, la

M<0.3 L'écoulement est incompressible, où les effets sont négligeable.

Les valeurs numériques données ci-dessus sont seulement des guides

1-7-Les équations fondamentales appliquées aux fluides compressibles:

1-7-3-Equation de quantité de mouvement:

DEPARTEMENT GENIE MECANIQUEModule : Gazodynamique

Chapitre 3. Ecoulement dans une conduite à Section variable –

x Comprendre L’Effet de variation de la section de la conduite

Page 31 Introduction à la dynamique des gaz Dr ALLICHE Mounir, Univ Medea

Ecoulement dans une conduite à Section variable –

3.1- Equations de base:

Les équations de base sont;

Page 32 Introduction à la dynamique des gaz Dr ALLICHE Mounir, Univ Medea

En dérivant cette équation, on obtient :

0 du  Pdv  vdP  VdV (3.5b)

-Variation de l'énergie interne; dm u 2 u1

0 du Pdv vdP VdV (3.5b)

( qi U i S VMAX 8 ) donnerait une diminution du débit mathématiquement ; ce

U1 V1 1,17 65,3