Tuyère de Laval

La tuyère de Laval est un tube en forme de sablier utilisé pour accélérer des gaz chauds et sous pression qui le traversent jusqu'à ce qu'ils atteignent une vitesse supersonique. La tuyère convertit de manière optimale la chaleur des gaz en énergie cinétique. Elle permet de produire de grandes quantités d'énergie à partir de gaz de combustion. Des tuyères de Laval sont utilisées dans les moteurs-fusées, les turbines à vapeur et les turbines à gaz. Dans le cas d'un moteur-fusée ce type de tuyère joue un rôle fondamental dans l'optimisation de la poussée en maximisant la vitesse d'éjection des gaz. La tuyère de Laval doit son nom à l'ingénieur suédois Gustaf de Laval qui en a découvert le principe en 1887^[1].

Historique

Gustaf de Laval construisit en 1887 une petite turbine à vapeur pour prouver que de tels appareils peuvent être fabriqués dans des dimensions réduites, et en 1890 il développe une tuyère permettant d'augmenter la vitesse de la vapeur entrant dans la turbine.

Caractéristiques d'une tuyère de Laval

Une tuyère de Laval est un tube dans lequel circule un gaz ayant la forme d'un sablier : son diamètre commence par se réduire (dans le sens de circulation du gaz) puis augmente à nouveau. Il comprend trois parties :

le convergent : c'est la partie de la tuyère qui va en se rétrécissant,
le col est la section de la tuyère où le diamètre est minimum,
le divergent dont le diamètre s'accroit à nouveau.

Pour qu'une tuyère de Laval parvienne à accélérer des gaz de manière optimale, il est nécessaire que le convergent et le divergent (qui ne sont pas symétriques) aient des formes bien précises et que le diamètre du col prenne une valeur donnée. Tous ces paramètres sont déterminés à partir des caractéristiques du gaz entrant (pression, température, débit, masse moléculaire) et de la pression externe.

Principe de fonctionnement d'une tuyère de Laval

Le principe de fonctionnement d'une tuyère de Laval repose sur les propriétés des gaz lorsqu'ils circulent aux vitesses subsonique et supersonique. Lorsqu'un gaz circule à une vitesse subsonique dans un tuyau dont le diamètre se rétrécit, sa vitesse augmente. La vitesse du gaz ne peut toutefois pas dépasser celle du son (Mach 1). En effet en régime d'écoulement supersonique (vitesse supérieure à la vitesse du son) le comportement du gaz s'inverse : pour que sa vitesse augmente il faut que le diamètre du tuyau augmente (démonstration plus bas : Relation d'Hugoniot). Pour accélérer un gaz à des vitesses supersoniques, il faut donc qu'il circule d'abord dans une section de tuyau convergente jusqu'à ce qu'il atteigne la vitesse Mach 1 et à partir de cette section du tuyau, qu'on appelle le col, le gaz doit progresser dans un tuyau de diamètre croissant (le divergent) pour que la vitesse continue à augmenter.

La tuyère de Laval ne fonctionne selon ce principe que si la vitesse du gaz atteint la vitesse Mach 1 au niveau du col. Pour y parvenir il faut que la tuyère soit conçue de manière que la pression en sortie soit au minimum deux fois plus faible que celle en entrée (P/P_o=0,53 exactement pour de l'air). Si cette condition est remplie, la vitesse au col atteint Mach 1 et la tuyère est dite amorcée. Si la pression en sortie est plus forte que cette valeur, la tuyère ne s'amorce pas. Au contraire si le rapport est plus important le rendement augmente. Celui-ci est optimal lorsque la pression en sortie est égale à la pression ambiante (au niveau du sol (1 bar)) : on dit alors que la tuyère est adaptée. Pour un moteur-fusée le rapport de section du divergent doit donc être d'autant plus important que le moteur fonctionne à des altitudes élevées c'est-à-dire à des pressions ambiantes faibles.

Exemple de paramètres de fonctionnement d'une tuyère *adaptée*
Section de la tuyère	Pression	Température	Vitesse d'écoulement du gaz
Entrée du convergent	80 bars	2 700 °C	~ 0
Col de la tuyère	44 bars	2 400 °C	Mach 1
Sortie du divergent	1 bar	1 075 °C	3700 m/s

Calcul de la vitesse des gaz expulsés

La modélisation du comportement des gaz dans une tuyère de Laval repose sur les concepts et hypothèses simplificatrices suivants :

le gaz en entrée de la tuyère de Laval se comporte comme un gaz parfait ;
l'écoulement est isentropique, c'est-à-dire que l'entropie est constante, conséquence de l'hypothèse qu'il s'agit d'un fluide non visqueux et que le processus est adiabatique (c'est-à-dire qu'il n'y a pas d'échange de chaleur entre le fluide et la tuyère) ;
le débit de gaz est constant durant toute la période de fonctionnement ;
l'écoulement du gaz est non-turbulent et symétrique tout au long de l'axe de la tuyère ;
le gaz est compressible.

La vitesse des gaz expulsés peut être calculée en utilisant la conservation de l'enthalpie totale du fluide lors de son passage dans la tuyère qui provient du fait qu'il n'y a pas d'échange ni de travail, ni de chaleur entre le gaz accéléré et son environnement :

h+e_{c}=C_{p}T+{\frac {1}{2}}v^{2}=C_{p}T_{0}=cste

où :

$h$ est l'enthalpie du fluide (en tous points de sa trajectoire) (J),
$e_{c}$ est la vitesse cinétique du fluide (en tous points de sa trajectoire) (m².s^-2),
$C_{p}$ est la capacité thermique du gaz à pression constante (J.K^-1), supposée indépendante de la température,
T est la température statique en tous points à l'intérieur de la tuyère (K),
v la vitesse du gaz en tous points à l'intérieur de la tuyère (m.s^-1),
T₀ est la température à l'entrée de la tuyère (K),

En utilisant le fait que $C_{p}={\frac {\gamma R}{M(\gamma -1)}}$ , et que la transformation du gaz est isentropique qui implique que ${\frac {T}{T_{0}}}=\left({\frac {P}{P_{0}}}\right)^{(\gamma -1)/\gamma }$ , on trouve l'équation suivante fournissant la vitesse d'éjection du fluide ^[2]^,^[3]^,^[4] :

v_{e}={\sqrt {2C_{p}T_{0}\left[1-\left({\frac {P_{e}}{P_{0}}}\right)^{\frac {\gamma -1}{\gamma }}\right]}}={\sqrt {{\frac {T_{0}R}{M}}\cdot {\frac {2\gamma }{\gamma -1}}\cdot \left[1-\left({\frac {P_{e}}{P_{0}}}\right)^{\frac {\gamma -1}{\gamma }}\right]}}

avec :
$v_{e}$	= Vitesse des gaz à la sortie de la tuyère en m/s
$T_{0}$	= Température totale à l'entrée de la tuyère en kelvin
$R$	= Constante universelle des gaz parfaits
$M$	= Masse molaire du gaz
$\gamma$	= ${\frac {c_{p}}{c_{v}}}$ = Indice adiabatique
$C_{p}$	= Capacité thermique du gaz à pression constante
$C_{v}$	= Capacité thermique du gaz à volume constant
$P_{e}$	= Pression statique du gaz en sortie de tuyère
$P_{0}$	= Pression totale du gaz à l'entrée de la tuyère

Voici quelques valeurs typiques de vitesse de gaz en sortie de tuyères v_e de moteurs-fusées brulant différentes combinaisons d'ergols :

1 700 à 2 900 m/s pour des moteurs-fusées à monergol liquide
2 900 à 4 500 m/s pour des moteurs-fusées à ergols liquides
2 100 à 3 200 m/s pour des moteurs à propergol solide.

Lorsque la section de sortie tend vers l'infini la pression statique P_e tend vers zéro et la vitesse tend vers une valeur finie appelée vitesse de Crocco

v_{e,max}={\sqrt {2C_{p}T_{0}}}

On introduit parfois le paramètre de similitude ou nombre de Crocco

Cr={\frac {v_{e}}{v_{e,max}}}={\sqrt {1-\left({\frac {P_{e}}{P_{0}}}\right)^{\frac {\gamma -1}{\gamma }}}}=\left[1+{\frac {2}{(\gamma -1)M_{e}^{2}}}\right]^{-{\frac {1}{2}}}

Relation d'Hugoniot

Article détaillé : Équation d'Hugoniot.

Le comportement des gaz dans une tuyère de Laval repose sur le principe d'accélération des gaz décrit par l'équation d'Hugoniot (1885).

{\frac {dS}{S}}=({\mathcal {M}}^{2}-1){\frac {dv}{v}}

où S est la section du conduit, v la vitesse et ${\mathcal {M}}={\frac {v}{c}}$ le nombre de Mach.

À l'entrée de la tuyère, en régime subsonique, on reconnaît l'effet Venturi. Après le col, le gaz est en régime supersonique, et l'évasement du conduit augmente encore la vitesse d'éjection.

Démonstration

La conservation du flux de masse $q=\rho vS$ s'écrit :

${\frac {d\rho }{\rho }}+{\frac {dv}{v}}+{\frac {dS}{S}}=0$

La conservation de l'énergie s'écrit :

$d({\frac {v^{2}}{2}})+{\frac {dp}{\rho }}=0$

d'où vient :

$v\,dv+{\frac {dp}{d\rho }}\,{\frac {d\rho }{\rho }}=0$

On reconnaît l'expression de la vitesse des ondes sonores $c^{2}={\frac {dp}{d\rho }}$

$v\,dv+c^{2}\,{\frac {d\rho }{\rho }}=0$

Puis on substitue l'expression de $\rho$ de la conservation de masse, il vient :

$v\,dv=c^{2}\,({\frac {dv}{v}}+{\frac {dS}{S}})$

Finalement

$({\frac {v^{2}}{c^{2}}}-1){\frac {dv}{v}}={\frac {dS}{S}}$

Débit massique

Le débit de la tuyère ${\dot {m}}=\rho Sv$ peut être calculé au col d'aire $S_{c}$ où le nombre de Mach est égal à l'unité. Pour cela on utilise les relations d'écoulement isentropique depuis le réservoir caractérisé par une pression génératrice (totale) $P_{0}$ et une température génératrice (totale) $T_{0}$ . Pour un gaz parfait :

{\dot {m}}=P_{0}\,S_{c}\,{\sqrt {\frac {M\gamma }{RT_{0}}}}\left({\frac {\gamma +1}{2}}\right)^{\frac {\gamma +1}{-2(\gamma -1)}}

Toutes choses égales par ailleurs le col fixe le débit massique de gaz.

Notes et références

↑ A Brevet US 522066 A Steam Turbine, 1894
↑ Richard Nakka's Equation 12.
↑ Robert Braeuning's Equation 1.22.
↑ (en) George P. Sutton, Rocket propulsion elements : an introduction to the engineering of rockets, New York/Brisbane etc., Wiley-Interscience, 1992, 6^e éd., 636 p. (ISBN 0-471-52938-9)

Voir aussi

(en) John D. Anderson, Jr., Fundamentals of Aerodynamics, New York/St. Louis/Paris etc., McGraw-Hill Education, 1991, 772 p. (ISBN 0-07-001679-8)
Étude numérique de phénomène de décollement d'ondes de choc dans une tuyère supersonique de type CD convergente-divergente, Latreche Takieddine, Rachedi Amine [1]

Articles connexes

[1] A Brevet US 522066 A Steam Turbine, 1894

[2] Richard Nakka's Equation 12.

[3] Robert Braeuning's Equation 1.22.

[4] (en) George P. Sutton, Rocket propulsion elements : an introduction to the engineering of rockets, New York/Brisbane etc., Wiley-Interscience, 1992, 6^e éd., 636 p. (ISBN 0-471-52938-9)

[1]

[2]

[3]

[4]