DYNAMIQUE DES FLUIDES COMPRESSIBLES

I - INTRODUCTION AUX FLUIDES COMPRESSIBLES

Jusqu’à maintenant, les fluides étudiés sont supposés incompressibles. Cette

hypothèse se justifie dans le cas où les fluides sont des liquides. Dans le cas des gaz, la
masse volumique varie en fonction de la pression et de la température selon la loi des gaz
parfaits. Pour étudier les fluides compressibles, une revue thermodynamique est
nécessaire.

II - RAPPEL THERMODYNAMIQUE

La pression, la température et la masse volumique des gaz sont liées par une
équation d’état. Dans le cas où le gaz est parfait cette équation est :

P = rrT (1)

avec r = R/ M où R est la constante des gaz parfaits R= 8314 N.m.mol -1.K-1 et M est la
masse molaire du gaz. Dans le cas de l’air r = 287 N.m.kg -1.K-1.

II-1 Energie interne

les gaz parfaits ont d’autres caractéristiques simples et utilisables. En général,

l’énergie Interne d’une substance peut être exprimée comme U = u (v, T) alors :

(2)

avec v : 1/ : volume spécifique.

La chaleur spécifique à volume constant est définie comme alors :

(3)

Vous avez déjà vue que pour certaine substance satisfaisant à P = rrT (gaz parfait)
que , Par conséquent :

Dynamique des fluides compressibles 1

II-2 Enthalpie
L’enthalpie est définie comme : h = U +P/r. Pour un gaz parfait : h = U + r.T .
Or U = f(T), donc h est une fonction uniquement de la température : h = g(T)
Pour obtenir une relation entre h et T, on exprime h sous la somme générale
comme
h = h(P, T); donc :

et
(4)

avec

Nous avons vu que pour un gaz parfait h = f(T) uniquement. Par conséquent
, alors : .

Or ; donc

Soit : CP - Cv = r

Le rapport des chaleurs spécifiques est défini comme : ; on obtient alors :

Pour un gaz parfait, la chaleur spécifique est fonction de la température. Dans une
marge de température raisonnable, la chaleur spécifique d’un gaz parfait peut être
considérée constante. Sous cette considération on a :

(5)

et

(6)

II - 3 Entropie

L’entropie est très utilisable dans l’étude des fluides compressibles. Elle est définie
par :

Dynamique des fluides compressibles 2

 Or (7)

L’inégalité de clausius, déduite du second principe de la thermodynamique, nous

donne que : . Comme conséquence du second principe, ce résultat peut être

étendue a :
ou

Si la Transformation est Si la T Transformation est.

réversible on a irréversible

Si la transformation est en plus adiabatique , alors :

dS = 0 (réversible + adiab)
et
dS > 0 (irréversible + adiab )

Par combinaison du premier et du second principe de la thermodynamique, on peut

obtenir :
TdS = dU + P.dv (8)

d’autre part, dU = d(h - Pv) = dh - v.dP - P.dv donc : T.dS = dh - v.dP (9)

Pour un gaz parfait, la valeur d’énergie peut être évaluée ainsi :

(10)

ou

(11)

Dans le cas où la chaleur spécifique est constante, l’intégration des équations (10) et (11)
nous donne :

(12)

ou

(13)

Dynamique des fluides compressibles 3

Si la transformation est isentropique (adiabatique + réversible) dS = 0, donc :

dU + P.dv = 0 Cv.dT + P.dv = 0

et pour un gaz parfait on a :
dh - v.dP = 0 CP.dT - v.dP = 0

donc (14)

Or d’après les équations donnant dS on a :

Sachant que CP - Cv = r et , donc : .

En remplaçant dT par son expression (équation (14)), on a :

(15)

après intégration de l'équation (15), on obtient :

Ln(P) + k.Ln(v) = Ln(c) = Cte où Ln(P) + Ln(vk) = Ln(c) = Cte

Cette équation nous permet d’écrire que : P.vk = Cte ou P/k = Cte.

II-4 Travail de compression ( Transformation irréversible )

Par définition, le travail de compression fourni au fluide est donné par l’intégrale :

(16)
Ce travail est fourni par l’extérieur.

Dynamique des fluides compressibles 4

 Dans le cas où le gaz est parfait :
- Si la transformation est isotherme ( T = cte ), on a :

P.v = Cte Þ (17)

- Si la transformation est polytropique on a :

P.vk = Cte Þ (18)

II-5 Travail de transvasement (Transformation Réversible )

Soit la transformation suivante :

Figure 1

dans cette transformation, le travail total de compression reçu par le fluide est donné par :

(19)
Ce travail se compose comme suit : - Travail fourni par le piston P1 (P1v1)
- Travail fourni par la pression P2 (- P2v2)
- Travail wT fourni par le passage de machine T.

donc : Þ (20)

III- PROPAGATION D'UNE ONDE SONORE

III -1 Vitesse de Son

La célérité du son est un paramètre fondamental de la dynamique des fluides

compressibles.

Dynamique des fluides compressibles 5

 Soit une onde sonore qui se déplace dans une conduite suivant x où le fluide est
stationnaire (Vx = 0) :

Figure 2

On s’intéresse à la propagation de l’onde à la vitesse C et la variation des propriétés

du fluide de part et d’autre de l’onde. Avant l’onde, la pression et la densité du fluide sont
P et r. Le passage de l’onde provoque une variation de P et r , de dP et dr. La vitesse
après l’onde devient Vx + dVx soit dVx car Vx = 0.

Pour une observateur fixe, il voit que l’onde évolue dans la conduite et le système
parait instationnaire. Cependant, on peut placer un volume de contrôle autour de l’onde et
que l’observateur est fixé à l’onde. Pour celui-ci, l’onde apparaît stationnaire et que le
fluide s’approche du volume de contrôle à la vitesse (C) et s’éloigne à la vitesse C - dV x.

Figure 3

Les équations de base appliquées à ce volume de contrôle sont :

* Equation de continuité :
Hypothèses (1) écoulement permanent
(2) vitesse uniforme dans une section.

Equation de base :

(21)
d'après (1) cette
intégrale = 0

après intégration de l'équation (21), on obtient :

Dynamique des fluides compressibles 6

ou :
On peut négliger ddVx par rapport à d ou dVx, on obtient alors :
;

ou : (22)

* Equations de conservation de quantité de mouvement :

Equation de base : (23)

avec Fv : force de volume Ff : force de frottement FP : force de Pression

Hypothèses (3) force de volume négligeable Þ Fv = 0,

(4) force de frottement négligeable Þ Ff = 0.

D'après la conservation du débit masse, on a : .

La force de pression s'écrit : FP = P.A - (P + dP).A = - A.dP

L’équation (23) de la quantité de mouvement devient :

donc : (24)

En combinant les équationbs (22) et 24), on obtient :

Soit (25)

La vitesse de propagation de l’onde sonore est alors donnée par :

La pression et la masse volumique dans un écoulement polytropique sont liées par :

donc : Ln(P) - kLn() = Cte

Dynamique des fluides compressibles 7

Soit :

Alors :

Pour un gaz parfait , donc :

(26)

Si on compare la vitesse du son à la vitesse du fluide, on introduit un nouveau

paramètre sans dimension nommé nombre de Mach : .

- Si M >1 : l’écoulement est dit supersonique

- Si M = 1 : l’écoulement est sonique
- Si M < 1 : l’écoulement est subsonique .

IV- ETUDE D'ECOULEMENT MONODIMENSIONNEL D'UN FLUIDE

COMPRESSIBLE

Les propriétés du fluide dans un écoulement compressible sont affectées par :

- La surface d’échange.
- Les forces de frottement visqueuses
- et les transferts de chaleur avec le milieu extérieur.

Dans cette partie, nous étudions en premier temps les écoulements isentropiques
(viscosité et échange de chaleur négligeables) dans une conduite à section variable. Puis,
nous étudions les écoulements adiabatique des fluides visqueux dans une conduite à
section constante. Enfin on tient compte des échanges de chaleur avec le milieu extérieur.

IV-1 Ecoulement isentropique ( s = cte)

IV-1-1 Equations de base

On considère un écoulement isentropique de fluide incompressible à travers une

conduite.

Dynamique des fluides compressibles 8

 Figure 4

* Equation de continuité :
Hypothèses (1) écoulement permanent.
(2) écoulement unidirectionnel .

(27)

* Equation de quantité de mouvement :

(28)

Hypothèses (3) effet de la gravité négligeable Fv = 0

(4) force de Frottement Ff = 0

En considérant les hypothèses (3) et (4), l'équation (28) devient :

(29)

er
*1 Principe de la Thermodynamique :

Q + wfourni = E
Chaleur échangée Travail Fourni au Variation de l’énergie
avec l’extérieur Système du système

Hypothèses (5) : transformation isentropique et adiab Þ Q = 0

Wfourni = P1 v1 - P2 v2 + WT (30)

Dynamique des fluides compressibles 9

wT est nul car il n’y a pas de machine, E est l'énergie mécanique totale : .

En considérant l'hypothèse (3), l'énergie mécanique totale devient : .

L’équation énergie s'écrit alors :

(31)

En introduisant la notion d’enthalpie h = U + P.v, on obtient : .

Soit : : Formule de Zéner (32)

* Second Principe de la Thermodynamique :

L’écoulement étant isentropique donc : S1 = S2 = S = Cte.

IV-1-2 Cas où le fluide est un gaz parfait :

Les équations de base sont :

Continuité :
Quantité de Mouvement :

Premier principe :

Second principe : S1 = S 2 = S

Equation d’état :

Processus isentropique :

Si on connaît les propriétés à l’état (1), on aurait sept inconnues (2, A2, V2, P2,
h2, S2, T2) dans six équations. Cependant, on connaît une relation entre h et T. Pour un
gaz parfait :
dh = CP.dT
Pour CP = cte, on a : h = h2 - h1 = CP.T = CP.(T2 - T1)

L’équation d’énergie devient, pour un gaz parfait :

Dynamique des fluides compressibles 10

 (33)

IV-1-2-1 Conditions génératrices :

Si dans une région de l’espace de l’écoulement on a V = 0, l’enthalpie h 0 = CP.T0.

Cette région est appelée état générateurs. Dans cette état, les propriétés seront menues
d’indice 0. En appliquant l’équation (33) d’énergie entre l’état générateur et un point
quelconque on a :

Sachant que , on a : (34)

ou

D'après la définition de la vitesse du son , on a:

Or , on a :

On obtient finalement : (35)

(36)

(37)

Dynamique des fluides compressibles 11

IV-1-2-2 Relation de Barré et Saint -Venant :

L’équation d’énergie (32) est donnée par la formule de Zéner . Pour

un gaz parfait, on a : . En appliquant cette équation entre l’état

générateur et un état quelconque on a :

Donc

En introduisant la loi d’isentropique , on obtient :

équation de Saint Venant

(38)

IV-1-2-3 Etat critique :

Les conditions critiques correspondent à l’état où la vitesse de l’écoulement est

égale à celle du son, Soit M = 1. En appliquant les équations (35), (36) et (37) entre l’état
générateur et l’état critique, on obtient pour k = 1,4 :

La vitesse critique est donnée pour :

Dynamique des fluides compressibles 12

 Il est possible d’exprimer la section le long de la conduit en fonction de la section
où le nombre de Mach M = 1 ( Section Critique ). En effet ,

donc :

Þ (39)

Le traçage du Rapport en fonction de M montre que pour une valeur de

correspond à deux nombres de Mach. Le minimum est relatif au M = 1.

Figure 5 : Variation de en fonction du nombre de Mach

(cas d'une transformation isentropique k = 1.4)

IV-1-3 Effet de la variation de la section sur les propriétés de l’écoulement

(Transformation isentropique)

Dynamique des fluides compressibles 13

 Nous essayons de déterminer l’influence de la variation dA de la section sur la
vitesse V et la pression P.
Est ce que, pour une variation de dA de la section, dV et dP sont positif ou négatif?
Pour répondre à cette question, il est commode de travailler avec la forme
différentielle des équations régissant l’écoulement. Pour cela, on considère un tube de
section variable selon la direction x .

Figure 6

Hypothèses considérées : (1) Régine permanent

(2) Ecoulement uniforme dans une section.
(3) Force de gravité négligeable
(4) Force de frottement négligeable.

équation de continuité : (40)

conservation de la quantité de mouvement :

(41)

Calcul de FP?

On suppose que sur l’élément de A est appliquée la pression moyenne P + dP/2.

Le développement de cette équation s'écrit :

Dynamique des fluides compressibles 14

or , donc . (42)

En remplaçant FP par son expression dans l'équation (41) de conservation de la quantité de

mouvement, on obtient :

(43)

soit

(44)

En devisant l'équation (44) Par V2, on obtient : (45)

D’autre part, d’après l’équation de continuité on a : . En appliquant le

Logarithme sur cette équation on a :

Ln() + Ln(Vx) + Ln(A) = Cte

La dérivée de cette équations donne :

ou
. (46)

D'après l’équation (45) de conservation de qualité du mouvement, on a :

Or , donc : (47)

D’après cette équation, on voit que :

Dynamique des fluides compressibles 15

 Si M < 1 les variations de la pression dP et de la section dA ont le même signe
Si M > 1 dP et dA ont des signes opposés.

L’équation (47) s’écrit autrement :

: relation de HUGONIOT (48)

D’après l'équation (48), on remarque que :

- Si M > 1 dV et dA de même signe,

- Si M < 1 dV et dA de signe opposé.
- Si M = 1, dA /dV = 0 : ceci veut dire que la section passe pour un
extremum, si M = 1 d'après la figure 5, on voit que M = 1, peut être
obtenu uniquement au col ou à une section minimale.

Ces Résultats sont résumes sur la figure 7.

Figure 7

IV-1-4 Ecoulement isentropique dans une tuyère convergente

Dynamique des fluides compressibles 16

 Dans cette partie, on s’intéresse au comportement du rapport P/P 0 dans une tuyère
convergente débouchante dans un milieu à pression fixée et contrôlée par une vanne. Les
résultats sont montrés dans la figure 8 :

Figure 8

On remarque que :
i la vanne est fermée, il n'y a pas d’écoulement dans la tuyère et P = P 0 tout le long
de la tuyère.

ii Pb < P0, la vanne est ouverte; il y a écoulement subsonique à la sortie M < 1 et

Pe = P b .

iii Si Pb continue à diminuer par rapport à P 0; la vitesse de l’écoulement croit et la

pression à la sortie diminue. L'écoulement est subsonique et Pe = Pb.

iv Si Pb diminue de sorte à avoir un écoulement sonique au col (Me=1 et Pe=Pb=P*).

v Si Pb < P*, on a Me = 1 et Pe = P*. Ceci n’affecte pas l’écoulement dans la tuyère

(P et Pe) . On dit que la tuyère est choquée et l’écoulement n’est plus considéré
comme isentropique à la sortie de la tuyère.

L’écoulement à travers une tuyère convergente peut être diviser en deux régimes :
Régime (1) : 1 ³ Pb/P0 ³ P*/P0 : écoulement isentropique tout le long de la tuyère.
Dynamique des fluides compressibles 17
 Régime (2) : Pb/P0<P*/P0 : écoulement isentropique jusqu’au col; mais non
isentropique après la sortie.

L’écoulement correspondant au régime (2) peut être montré sur le diagramme (T,
S) (figure 9) :

Figure 9

IV-1-5 Ecoulement isentropique dans une tuyère convergent-divergente

Soit une tuyère convergente-divergente débouchant dans un grand réservoir où la

pression Pb est contrôlée par une vanne. L’état générateur est considéré constant. La
pression à la sortie de la tuyère est notée P e. On s’intéresse à étudier l’effet de la variation
de Pb sur la distribution de la pression à l’intérieur de la tuyère. Les résultats sont
représentées sur la figure-10 :

Figure 10
Initialement, la vanne est fermée; il n’y a pas écoulement dans la tuyère et la
pression est P0 dans la tuyère .

Dynamique des fluides compressibles 18

i: En ouvrant légèrement la vanne, on obtient la courbe (i). En tout point de cette
courbe l’écoulement est subsonique et le fluide peut être considéré comme
incompressible.

ii : En ouvrant d’avantage la vanne, le débit croit et on obtient la courbe (ii). Dans ce

cas l’effet de la compressible devient important, l’écoulement reste encore
subsonique en tout point et la décélération du fluide prend lieu dans le divergent.

iii : Si la vanne s’ouvre encore, on obtient la courbe (iii). Au col, le nombre de Mach
atteins l’unité M = 1 .

Dans les trois cas i, ii, iii l’écoulement est isentropique.

iv : en diminuant encore la Pression Pb, on obtient toujours M = 1 au col, mais le fluide

devient accéléré au divergent et on obtient un écoulement supersonique.
Si Pb = Pe l’écoulement est isentropique (courbe iv).
Ce résultat a été déjà vue dans la courbe donnant A/A* en fonction de M, à savoir
on a deux valeurs de M pour un même rapport A/A* (un en subsonique et l'autre
supersonique)

v: Si Pb diminue encore de sorte à avoir Pb ≠ Pe. Dans cette condition, il y a une

discontinuité entre la sortie de la tuyère et le réservoir; l’écoulement n’est plus
isentropique à la sortie.

IV-2 Ecoulement visqueux dans une conduite adiabatique à section constante

Considérons une conduite à section constante adiabatique parcourue par un fluide

compressible visqueux.

Figure 11
Dans cette paragraphe, on considère les hypothèses suivantes :

Dynamique des fluides compressibles 19

 - Régime permanent.
- Vitesse uniforme dans une section
- Effet de gravité négligeable
- Ecoulement adiabatique

En considérant ces hypothèses les équations de bases s'écrivent :

(1) Continuité : (49)

(2) Quantité de Mouvement : (50)

(3) Premier principe : Þ h01 = h02 (51)

car la conduite est adiabatique

(4) Second principe : (gaz parfait) (52)

(5) Equation d’état : (53)

(6) Dans le cas des gaz parfaits : h = CP.T = CP.(T2 - T1) (54)

Si on connaît les propriétés de l’état (1), on a six équations avec 7 inconnues ( T 2,

P2, 2, V2, h2, S2, Fx, ).

IV-2-1 Ligne de Fanno :

La résolution de ce système montre qu’il y a un nombre infini de solutions possible

de l’état (2). En traçant toutes les solutions de l’état (2), Calculées à partir de la
connaissance de l’état (1), sur le diagramme TS on obtient la ligne de Fanno représentée
sur la figure 12.

La plus simple méthode pour déterminer cette courbe est la suivante : on fixe
différentes valeurs de T2 et on calcule les autres propriétés à partir du système d’équations
précédent (49) à (54).

Dynamique des fluides compressibles 20

 Figure 12 :Schématisation de la ligne de Fanno représentée sur le diagramme TS

On remarque :
* Pour S maximale M = 1
* La branche supérieure M < 1
* La branche inférieure M > 1
* Si l’écoulement initial est subsonique, état (1), l’effet de la viscosité fait
augmenter le nombre de Mach (M 1); l'entropie S augmente aussi sous
l'effet de la viscosité.
* Si l’écoulement initiale est supersonique, état (1’), l’effet de la viscosité fait
diminuer le nombre de Mach (M 1).
* En tout point de la conduite, on peut imaginer un état générateur entre deux
points différents (1) et (2) P01 > P02 pour contre T01 = T02 car Q = 0.

Finalement l’effet de la viscosité est résumé dans le tableau suivant :

Subsonique Supersonique
Propriétés M<1 M>1 Obtenu à partir de
Température génératrice Constante Constante Equation d'énergie
T0
Entropie S Augmente Augmente Second principe
Pression Diminue Diminue T0 = constante;
génératrice P0 S augmente
Température T Diminue Augmente ligne de Fanno
Vitesse V Augmente Diminue Equation d'énergie et la
tendance de T
Nombre de Mach M Augmente Diminue Tendance de V, T et la
définition de M
Masse volumique  Diminue Augmente Equation de
Continuité
Pression P Diminue Augmente Equation d'état et
effet sur  et T

Dynamique des fluides compressibles 21

 IV-2-2 Détermination des tables de Fanno
(pour un fluide obéissant à la loi des gaz parfait )

La variable primaire indépendante dans la ligne de Fanno est la force de frottement

Ff. La connaissance de la force de frottement totale entre deux points sur la ligne de Fanno
permet de prédire les propriétés de l’état (2) si on connaît celles de l’état (1). La force de
fortement totale est obtenue en intégrant tout le long la surface latérale de la conduite.
Puisque cette force de frottement varie le long de la conduite, on doit développer une
équation différentielle et intégrer pour trouver l’évolution Ff. Pour établir cette équation
différentielle, on utilise un volume de contrôle élémentaire (figure 13) :

Figure 13

En considèrent les hypothèses citées précédemment, on obtient :

* Equation de continuité :

on suppose que :

on obtient : (55)

* Equation de qualité de mouvement :

En considérant la conservation du débit masse, on obtient :

(56)

* Premier principe de la Thermodynamique :

Dynamique des fluides compressibles 22

 Soit en substituant h = U + Pv, on obtient : (57)

Pour compléter notre formulation, on doit relier la force de frottement au variable

d’écoulement. Or :

(58)

avec Pw est le périmètre mouillé de la conduite.

Calcul de w : Soit autour du point P, on construit un volume de contrôle.

Figure 14

En appliquant le théorème de quantité de mouvement sur le domaine de contrôle, on a :

(59)

On suppose que dx est très faible de sorte à avoir V2 = V1.

Fv : force de volume négligeable donc .

donc : (60)

Soit : (61)

Dans le cas où la conduite n’est pas cylindrique, on utilise le diamètre hydraulique

Dh=4A/P. A est la section de passage et P est le périmètre mouille. En introduisant la

notion de coefficient de perte de change  on a : ; donc : .

Dynamique des fluides compressibles 23

En remplaçant par son expression dans dFf, on a :

ou (62)

en remplaçant l’expression de dFf dans l’équation (56) de conservation de la quantité de

mouvement, on a :

(63)

En divisant par P l'équation (63), on a :

Or pour un gaz parfait et , on obtient :

(64)

Dans le but d’obtenir une équation liant M et x on doit éliminer et de

l’équation (64). D’après la définition du nombre de Mach :

donc :

(65)

A partir de l’équation de continuité : et (66)

Pour un gaz parfait , on peut écrire :

(67)
Dynamique des fluides compressibles 24
En combinant les trois équations (65), (66) et (67), on a : (68)

En substituant l’expression de dans l’équation (64) de quantité de mouvement finale,

on a :

Soit (69)

Cependant afin d’obtenir une expression de M en fonction de x on doit avoir une

expression de en fonction de M.

Or à partir de l’équation (35) d’énergie et pour un gaz parfait, on a :

Puisque sur la ligne de Fanno, la température génératrice est constante donc :

alors

(70)

en substituant l’expression de dans l’équation (69), on obtient :

Soit :

(71)

Dynamique des fluides compressibles 25

 L’intégration de l'équation (71) entre l’état (1) et l’état (2) permet d’obtenir une
expression compliquée liant M1 et M2. Cette expression peut être obtenue numériquement
pour chaque nouvelle combinaison M1 et M2. Le calcul de cette intégration peut être
simplifié considérablement en utilisant les conditions critique (où M = 1); c’est à dire en
considérant les bornes d’intégration M et 1.

Un nombre de Mach égale à 1 est obtenu quand la longueur de la conduite est très
grande (figure 16).

Figure 16

(72)

Peut être calculé par On définit un coefficient moyen l

intégration par partie avec

Remarque l est fonction de Reynolds

comme v = cte le long de la conduite
donc l est fonction uniquement de µ
L’intégration de l'équation (72) donne :

(73)

Cette équation donne la valeur maximale du terme correspondant à la

valeur du nombre du Mach M.

Dynamique des fluides compressibles 26

 Puisque est fonction uniquement de M, et sachant que la conduite est de

longueur L, le nombre de Mach initiale M1 et le nombre de Mach final M2, la valeur de

correspondante est obtenue comme suit :

(74)

Puisque T0 est contante, on peut écrire :

(75)

de même : (76)

A partir de la conservation de la masse (.V = *.V*), on a :

Pour un gaz idéal , on a: (77)

Le rapport entre la pression génératrice d’un état et la pression génératrice de l’état

critique est donné par :

donc : (78)

Dynamique des fluides compressibles 27

 Les rapports sont tablés en fonction du nombre de

Mach selon les formules établies précédemment.

IV-3- Ecoulement de fluide non visqueux (parfait) dans une conduite à section
constante échangeant de la chaleur

Afin d’étudier l’influence d’échange de chaleur, avec le milieu extérieur, sur

l’écoulement des fluides compressibles, nous établissons des équations de bases relatives à
l'écoulement de gaz parfait, monodimensionnel. Pour cela, considérons dans une conduite
à section constante un volume élémentaire :

Figure 17

Hypothèses : (1) régime permanent

(2) écoulement uniforme dans une section
(3) force de gravité négligeable
(4) force de viscosité négligeable.
En considérant ces hypothèses, les équations de base s’écrivent :

(1) Continuité : (79)

(2) Quantité de Mouvement : (80)

(3) Premier principe : Þ (81)

Dynamique des fluides compressibles 28

(4) Second principe : (gaz parfait) (82)

(5) Equation d’état : (83)

(6) Dans le cas des gaz parfaits : h = CP.T = CP.(T2 - T1) (84)

Si on connaît les propriétés de l’état (1), on a six équations avec 7 inconnues ( T 2,

P2, 2, V2, h2, S2, ).

IV-3-1 Ligne de Rayleigh

En raisonnant de la même manière que pour la détermination de la ligne de Fanno,

on fixe différentes valeurs de T 2 et on calcule les autres propriétés à partir du système
précédent. Les résultats de ce calcul sont tracés sur le diagramme Ts, on obtient une
courbe appelée "ligne de Raleigh" (figures 18 et 19).

On remarques :

1- au point où la température est maximale le nombre de Mach ,

2- au point ou l’entropie est maximale, M = 1;

3- sur la branche supérieure de la courbe, M est toujours inférieur à 1 et augmente en
se déplaçant à droite.
4- en tout point du branche inférieure, M est toujours supérieure à 1 et décroît en se
déplaçant à droite.
5- Par rapport à un point initiale de la courbe, en additionnant de la chaleur on se
déplace vers la droite le long de la ligne de Réglage et inversement si on
refroidisse.

Dynamique des fluides compressibles 29

 Figure 18 : Représentation sur le diagramme TS
de l'effet de l'apport ou perte de chaleur

Figure 19 : Réduction de la pression génératrice à cause du chauffage

pour deux cas d'écoulements

Dynamique des fluides compressibles 30

 Propriétés Chauffage Refroidissement Obtenu à partir
M<1 M>1 M<1 M>1
Température Augmente Augmente Diminue Diminue 1erprincipe et
génératrice T0 h0=CP T0

Entropie S Augmente Augmente Diminue Diminue Second principe

Pression Diminue Diminue Augmente Augmente
génératrice P0

Augmente Diminue
Température T A partir de la
Augmente Diminue ligne de
Diminue Augmente Rayleigh

Vitesse V Augmente Diminue Diminue Augmente Equation de la

quantité de
mouvement et
effet sur P
Nombre de Mach Augmente Diminue Diminue Augmente Tendance de la
M ligne de
Rayleigh
Masse volumique Diminue Augmente Augmente Diminue Equation de
 continuité et
effet sur V
Pression P Diminue Augmente Augmente Diminue Tendance de la
ligne de
Rayleigh

IV-3-2 Détermination des tables de Rayleigh pour un gaz parfait :

Comme dans le cas de la détermination des tables de Fanno, on considère un état de

référence (état critique où M = 1) et on détermine les différents rapports des propriétés
entre l’état critique et un état quelconque.

Le rapport est obtenu en utilisant l’équation de conservation de la quantité de

mouvement :

ou (85)

Pour un gaz parfait : , on a :

(86)

Or : donc :

Dynamique des fluides compressibles 31

Finalement : (87)

Pour un gaz parfait :

d'après l'équation de continuité :

alors : , donc : (88)

de même : (89)

Pour la température génératrice, on peut écrire :

(90)

de même pour la pression génératrice :

(91)

Dynamique des fluides compressibles 32