1st National Conference on Applied Energetics.

NCAE 2020, February 11 and 12, 2020, NAAMA, Algeria.

L’EFFET DE LA VARIATION DE LA GÉOMÉTRIE SUR LE

COMPORTEMENT THERMIQUE D’UNE TUYÈRE CONVERGENTE
DIVERGENTE
Rai Naceur Bouharkat1, Medjahed Driss Meddah2, Mahammedi Abdelkader3
1
Laboratoire de Technologie Industrielle Univ Ibn Khaldoun Tiaret . rainaceurbouharkat@[Link]
2
Laboratoire LaRTFM, Département de Génie Mécanique, ENP d’Oran .medjaheddriss@[Link]
3
Universié de Ziane Achour Djelfa. [Link]@[Link]

RÉSUMÉ
L'étude des flux dans les tuyères est un sujet de recherche motivée par la nécessité de comprendre et maîtriser
les phénomènes induits, Nous nous intéressons principalement à l'identification des caractéristiques des
écoulements compressibles et visqueux dans une tuyère convergente-divergente.
Notre contribution est une étude numérique de l'écoulement de l'air dans une tuyère avec variation de la géométrie
(deux tuyères avec section de sortie de ex47.6cm² et ex33.6cm² pour la longueur de tuyère de 63cm et deux
tuyères avec les mêmes sections précédentes pour la longueur de 83cm) pour étudier l'influence des conditions
thermiques pariétales sur le comportement des fluides.

Mots-clés : tuyère convergent divergent, Aval, Nombre de Mach.

NOMENCLATURE

Symbols : C vitesse, m/s

𝑚̇𝑔 débi global, Kg/s F la pousséé, N
𝑚̇𝑔 débi de l’air, Kg/s t le temps, s
𝑣𝑠 Vitesse des sortie des gaz, m/s T temperature, K
𝑣𝑒 Vitesse du véhicule, m/s M Nombre de Mach
P pressure, Nm-2 e epaisseur, m
Greek Letters: Indices / Exponents :
ρ density, kg.m-3 o grandeur totale
A surface, m²

1. INTRODUCTION

L’étude des écoulements gazeux est un domaine de recherche qui s’est récemment beaucoup développé,
notamment grâce à des avancées majeures dans le domaine d’usinage.

1.1 Les tuyères

Une tuyère ou tuyère propulsive, dans le domaine de l'astronautique, est un conduit (appelé aussi divergent) de
section conique (fusées de feu d'artifice) ou oblongue (tuyères de missiles, et lanceurs), placé à la sortie d'un
propulseur et qui sert à transformer en énergie cinétique l'énergie des gaz de combustion.

1.2 Forme des tuyères

Il y a plusieurs formes de tuyères nous citons parmi les plus connues et les plus utilisées celles de formes
convergentes, divergentes, et convergente-divergente, comme l’indique la figure (1).

Fig 1: La tuyère supersonique [1].

1.3 Applications des tuyères dans l’industrie

1.3.1 Microcanaux
Le concept des microcanaux n'est pas nouveau puisqu'il a été introduit vers 1980 par les chercheurs Tuckerman
et Pease [2], L’utilisation des canaux pour le transport de fluides (liquides, gaz) est fondamentale lors de
l’association de plusieurs microsystèmes et joue le même rôle que les fils électriques dans les circuits intégrés
conventionnels.
Il est aujourd’hui possible de fabriquer des canaux dans la silice à partir de nombreuses techniques de gravure
mais on peut aussi disposer de microcanaux en verre ou avec d’autres matériaux.
Les critères de choix les plus importants sont les contraintes géométriques (sections compatibles avec les autres
composants du système, longueur, et surface de la conduite) et le choix des matériaux constituant la surface
interne du canal (qui doit être compatible avec le fluide utilisé[3].

1.3.2 Turboréacteurs
La mécanique de propulsion s’intéresse surtout par les turboréacteurs, leurs fonctionnements et leurs types
(figure2), qui nous intéresse par cette étude, la formule de la force de poussé dans les turboréacteurs est calculée:
𝐹 = 𝑚̇𝑔 𝑣𝑠 − 𝑚̇𝑎 𝑣𝑒 (1)
Si on a au sol la vitesse de véhicule est nulle donc la poussée devient :
𝐹 = 𝑚̇𝑔 𝑣𝑠 (2)
Cette formule démentre l’importance de l’augmentation de la vitesse des gaz à la sortie.

Fig2 : Schéma d’un turboréacteur simple flux, sec mono corps [4]
1.4 travails dans le domaine
A cause de la grande utilité des tuyères dans les différents domaines, De nombreux chercheurs ce sont intéressé,
à faire des études théoriques et expérimentales. Ce domaine de recherche date d’une vingtaine d’années puisque
le premier article remonte à 1983, et depuis, le nombre d’articles sur le sujet croît de manière exponentielle.
Omer San , Ilhan Bayraktar , Tuba Bayraktar [5] se sont intéressés dans cette étude, à l’effet de taille sur le
comportement de l’écoulement. Deux méthodes sont utilisées pour analyser les paramètres de l'écoulement dans
une tuyère convergente divergente, l’une est basée sur les équations de Navier stockes et l’autre c’est l’équation
du Brunette augmentée.
Ameur Djilali, et ses colaborateurs [6] ont modélisé l’écoulement d’un gaz (azote) dans des micro-filtres. La
région du filtrage est une membrane qui contient des orifices, en fait des microcanaux, dont le diamètre est de
quelques micromètres. Au voisinage de chaque orifice, l’écoulement du gaz supposé normal à la membrane, est
en régime de raréfaction dit de transition, intermédiaire entre le régime continu relevant des équations de Navier-
Stokes et le régime moléculaire libre où les molécules du gaz évoluent sans interactions entre elles.
Jinliang Xu et Chuangxin Zhao[7] L’objectif de leur étude est l’utilisation des équations de Navier stokes avec
les conditions de glissement sur les parois pour simuler le débit supersonique dans une tuyère convergente
divergente. La première observation montre que les ondes de choc peuvent avoir lieu à l’intérieur ou à l’extérieur
de la micro tuyère sous la pression atmosphérique. La structure de L'onde de choc de l’écoulement aux
caractéristiques suivantes: Il y a un point de séparation de couche limite en aval de la gorge et il existe une
interface séparant les ondes de couche limite visqueux avec un flux de retour et le débit de base non visqueux.
Le nombre de Mach est d’environ 0.8 à la gorge supersonique dans le divergent, et il n’atteint pas l'unité
exactement à l'emplacement de la gorge, mais un peu en aval de la gorge. La couche limite visqueuse est très
épaisse et peut couvrir la moitié de la largeur de la tuyère sur le plan de sortie.
L’effet visqueux est important pour affecter et former la structure des flux d'ondes de choc dans les micro-
échelles à cause de la grande valeur rapport du volume sur surface. Au plan de sortie de la tuyère la couche limite
visqueuse peut couvrir la moitié de la sortie de l’ensemble de largeur de la tuyère.

L’objectif principal de ce travail est de voir l’influence des conditions thermiques pariétales sur le comportement
dynamique et thermique, d’un écoulement bidimensionnel de fluide visqueux dans une tuyère. Nous nous
sommes intéressées par l’augmentation de la vitesse des gaz à la sortie de la tuyère comme l’indique l’équation
(2), et pour réaliser cet objectif dans cette étude on a varié la géométrie de la tuyère (deux tuyères avec section
de sortie de ex47.6cm et ex33.6cm pour la longueur de tuyère de 63cm et deux tuyères avec les mêmes sections
précédentes pour la longueur de 83cm) pour voir quelle est la géométrie optimale, les différentes géométries
étudiées sont présentées au (figure 3).

Fig 3: les différentes géométries de l’étude

 2. FORMULATION DU PROBLÈME

L’étude de l’écoulement d’un fluide compressible devient plus compliquée que celle d’un fluide incompressible.
En effet, les variations de température ou de pression qui peuvent apparaître dans l’écoulement d’un liquide ne
modifient en rien les volumes mis en jeu car la dilatation ou la compression sont généralement négligeables. En
revanche, ces phénomènes prennent une grande importance lorsqu’il s’agit de vapeur ou de gaz.
L’étude de l’écoulement des fluides compressible ne peut être abordée sans avoir fixé au préalable un certain
nombre d’hypothèses simplificatrices (nature du gaz : parfait, type d’évolution : isotherme ou
adiabatique,…etc.[8].

2.1 Les équations générales:

Pour un écoulement bidimensionnel, visqueux, et compressible, les équations s’écrivent sous la forme suivante:
𝜕𝜌
1. Continuité : + ∇. 𝜌𝑣 (3)
𝜕𝑡
𝜕(𝜌𝑣)
2. Quantité de mouvement 𝜕𝑡
+ ∇. 𝜌𝑣𝑣 = −∇p + ∇τ (4)
𝜕(𝜌𝐸)
3. Conservation d’énergie 𝜕𝑡
+ ∇. 𝜌𝐸𝑣 = −∇(𝜌𝑣) + ∇(τ. 𝑣) − ∇q (5)
La dynamique des gaz est une étude de l’écoulement de fluides avec une vitesse comparable à la vitesse locale
sonique. C’est le cas quand la vitesse de fluide est 30% plus que la vitesse locale sonique, i.e. M > 0.30. Dans
ces cas, la densité n’est plus une constante.

Grandeurs totales :
Un point d’arrêt de l’écoulement est un point grandeur définie en un point d’arrêt.

Grandeurs soniques :
Un point sonique d’un écoulement est un point ou la vitesse est égale à la célérité locale du son. Une
grandeur sonique est une grandeur définie en un point aux conditions soniques[9].

2.2. Écoulement unidimensionnel isentropique avec variation de section

2.2.1 Effet de la variation de section :
Examinons l’effet de la variation de section dans le cas général d’un fluide d’équation quelconque Les équations
de conservation de la masse et de l’énergie sont respectivement :
𝑑𝜌 𝑑𝑐 𝑑𝐴
+ + =0 (6)
𝜌 𝑐 𝐴

𝑐2
𝑑 (ℎ + 2
) = 𝑑ℎ + 𝑐𝑑𝑐 = 0 (7)
2
L’écoulement étant isentropique 𝑑ℎ = 𝑣𝑑𝑝 = 𝑣𝑎 𝑑𝜌 Par conséquent, l’équation de l’énergie peut se
réécrire :
𝑑𝑝
𝑎2 𝑝
+ 𝑐𝑑𝑐 = 0 (8)
𝑑𝑝 𝑐 2 𝑑𝑐 𝑑𝑐 𝑐2 𝑑𝐴
𝑝
+ 𝑎2𝑐
= 0 => 𝑐 (1 − 𝑎2 ) + 𝐴
=0 (9)
Résolvant pour dρ/ρ et insérant dans l’équation de conservation de la masse on obtient :
𝑑𝑐 𝑑𝐴 𝑑𝑐 𝑑𝐴 −𝑑𝑝 𝑑𝐴
𝑐
(1 − 𝑀2 ) + 𝐴
= 0 => 𝑐
(1 − 𝑀2 ) = − 𝐴
= 𝑐
(1 − 𝑀2 ) = − 𝑝𝑐 2 (1 − 𝑀2 ) (10)
L’équation (10) à une grande importance, car elle permet de déterminer la forme à donner à la conduite :
Si Ma <1(vitesses subsoniques) la vitesse varie en sens inverse de la section 𝒅𝒗 > 𝟎 𝒔𝒊 𝒅𝑨 < 𝟎
Si Ma>1(vitesses supersoniques), la vitesse varie dans le même sens que la section 𝒅𝒗 > 𝟎 𝒔𝒊 𝒅𝑨 > 𝟎
 P diminue M>1 P diminue
M<1 A augmente
A diminue Supersonique
Subsonique

M<1 P augmente P augmente

A augmente
M>1
Subsonique Supersonique A diminue

Fig 4 Forme de tuyère en fonction du nombre de Mach [10].

2.4 Evolution des paramètres de l’écoulement :

Pour un écoulement isentropique, l’équation d’énergie entre deux états dans le champ d’écoulement devient :
𝑢12 𝑢22
ℎ1 + 2
= ℎ2 + 2
= ℎ0 = 𝑐𝑜𝑛𝑠𝑡 (11)
L’enthalpie d’un gaz parfait est donnée par h = CpT. En substituant dans l’équation d’énergie de l’écoulement
permanent, on trouve ho = Cp To = constante, Où ho appelée l’enthalpie de stagnation, reste constante dans
le champ d’écoulement
𝑇0 𝑦−1 2
𝑇
=1+ 2
𝑀 (12)

Donc, la température de stagnation reste constante dans un écoulement isentropique et la relation de la

température locale à la température de stagnation est une fonction seulement du nombre de Mach local.
En incorporant l’équation de vitesse sonique, l’équation d’état de gaz parfait et l’équation d’énergie on trouve
les relations utiles pour l’écoulement isentropique permanent des gaz parfaits [11].

𝑎0 𝑇 1/2 𝑦−1 2 1/2

𝑎
= ( 𝑇0 ) = (1 + 2
𝑀 ) (13)
𝑝0 𝑇 𝑦/(𝑦−1) 𝑦−1 2 𝑦/(𝑦−1)
𝑝
= ( 𝑇0 ) = (1 + 2
𝑀 ) (14)
𝜌0 𝑇 1/(𝑦−1) 𝑦−1 1/(𝑦−1)
𝜌
= ( 𝑇0 ) = (1 + 2 𝑀2 ) (15)

[Link]élisation du problème :
La complexité des géométries rencontrées dans les applications pratiques en mécaniques des fluides sont telles
qu’il existe un besoin bien défini pour le traitement de frontières irrégulières. Cependant, les approches
généralisées pour traitement de géométries arbitraire en 3 dimensions remontent en de début des années 70 avec
les travaux du professeur Joe [Link] [12] sur la génération des maillages adaptés.

3.1 Présentation du domaine d'étude

Dans cette étude nous avons étudié quatre domaines d’étude d’une tuyère dont les caractéristiques géométriques
sont indiquées sur la figure (3).
- Demi angle du convergent : 4.76°.
-La longueur : 83cm et 63cm
- Section de sorties : ex47.6cm² et ex33.6cm²
- section du col : ex10.5 cm².
- Le convergent et le divergent sont raccordés par un arc de cercle.
3.2 Conditions aux limites
Notre écoulement est délimités par des parois imperméables et les fluide est de l’air, donc c’est un fluide
visqueux, ce qui nous conduit à une condition aux limites de non glissement (u=v=0) sur les parois. Soient les
paramètres d’arrêt de l’écoulement : P0=238828.19Pa, T0=308.39 k.
Les propriétés thermiques du gaz sont représentées au (tableau1), l’étude est menée par ANSYS.

Table1 : Les propriétés physiques de l’air [13]

Propriétés Valeurs numériques

Masse volumique 1.225 Kg/m3

Viscosité dynamique 1.7894.10-5 Kg/ms

3.RÉSULTAT ET DISCUSSIONS

Vue la difficulté de traiter ce genre d’écoulement (compressible, visqueux, présence de couche limites), nous
commençons notre étude par une validation des résultats de simulation par les résultats analytiques (figure5).

3.1 Validation des resultats

a b

c d
Fig 5 : Variation des paramètres (a : vitesse, b : température, c : pression, d : densité) le long de la tuyère

Ces graphes nous démontre la bonne concordance entre les résultats de simulation et les résultats calculées pour
les déférents paramètres (a : vitesse, b : température, c : pression, d : densité), donc l’outil de simulation est
fiable pour l’utilisées dons les prochaines simulations.
3.2 Les contours
3.2.1 Contours de température
Avant de présenter les résultats de notre simulation en forme de graphes, nous avons présenté les différents
contours de températures pour voir la distribution de la température le long de la tuyère ainsi que les valeurs
thermiques à la sortie de la tuyère (figure 6).

a : la longueur de 83cm et surface de ex33.6cm² b : la longueur de 83cm et surface de ex47.6cm²

c : la longueur de 63cm et surface de ex33.6cm² d : la longueur de 63cm et surface de ex47.6cm²

Fig 6 : Contours de température

3.2.2 Contours de pression

a : la longueur de 83cm et surface de ex33.6cm² b : la longueur de 83cm et surface de ex47.6cm²

c : la longueur de 63cm et surface de ex33.6cm² d : la longueur de 63cm et surface de ex47.6cm²

Fig 7 : Contours de pression

3.2.3 Contours de vitesse

a : la longueur de 83cm et surface de ex33.6cm² b : la longueur de 83cm et surface de ex47.6cm²

 c : la longueur de 63cm et surface de ex33.6cm² d : la longueur de 63cm et surface de ex47.6cm²
Fig 8 : Contours des vitesses

Ces graphes montre que la géométrie qui nous donne une rapide et grande chute de température (figure6) et de
pression (figure7) est pour la géométrie (d : la longueur de 63cm et surface de ex47.6cm²), c’est un résultat
logique parce que dans cette géométrie (b) on a augmenté la section de sortie après une petite distance (rapide).

Rappelant que notre objectif est de trouver une vitesse de sortie maximale pour trouver une force de poussée
maximale comme l’indique (l’equation2), la géométrie qui vérifier cette vitesse est aussi la géométrie (d : la
longueur de 63cm et surface de ex47.6cm²) mais cette vitesse est répartie sur les côtés, si on prend en
considération les vitesses horizontales on trouvera une autre géométrie optimale comme nous montrons par la
suite.

3.4 Comparaison graphique des résultats

a b

c d
Fig 9 : Contours des vitesses (a : vitesse, b : température, c : pression, d : densité) le long de la tuyère
Dans ces graphes les paramètres thermiques sont étudiés dans la ligne horizontale de la tuyère et c’est normale
que les valeurs aux extrémités sont plus grandes que les valeurs aux centre, ces graphes nous montre que la
vitesse maximale est enregistrée par la tuyère (la longueur de 83cm et surface de ex47.6cm²).
La tuyère courte (la longueur de 83cm et surface de ex47.6cm²) nous donne la deuxième grande vitesse comme
on indiqué dans la partie des contours précédente.

[Link]

Ce travail de simulation des écoulements dans les tuyères a été dicté par le besoin de comprendre les
mécanismes qui régissent ces écoulements utilisant comme outil de calcul, le code Fluent.

Les résultats obtenus montrent l’influence conditions thermiques et pariétales sur la structure de l’écoulement,
en particulier dans les zones où les effets visqueux sont significatifs. Nous nous sommes intéressées aux effets
de changement de section et des changements des conditions d’échange par convection sur le comportement
dynamique et thermique de l’écoulement, les résultats de la simulation montre que plus nous augmentons la
section de sortie et plus nous augmentons le transfert chaleur nous remarquons une diminution de la température
donc une diminution l’enthalpie à la sortie qui est présentée par la relation h=cpT (Pour les gaz parfaits).

Pour le nombre de Mach (ou la vitesse) nous avons trouvé que ce dernier diminue pour les petites géométries,
est augmenté lorsque la température diminue à la sortie de la tuyère. Nous remarquons aussi que la pression
diminue avec l’augmentation de la section est avec l’apport de chaleur.

Nous avons remarqué que la grande section à la sortie génère plus de pertes, une grande vitesse donc une
géométrie optimale.

REFERENCES

[1] Yohann Perrot, thèse, étude mise au point et validation de modèles de turbulence compressible, Thèse de
doctorat, INSA de Rouen 2006.
[2] D.B. Tuckerman, R.F.W. Pease, Hight-performance heatsinking for VLSI. Electron Device Letters, 1981.
[3] David dersiaux, [Link] microstructures108/3/.(2004)).
[4] Cours de thermopropulsion ii, dr henni mansour zoubir, universite saad dahleb blida page p28
[5] Omer san , ilhan bayraktar , tuba bayraktar : size and expansion ratio analysis of micro nozzle gas flow
international communications in heat and mass transfer 36 (2009) 402–411
[6] Ameur djilali., modélisation analytique et simulation numérique par la méthode de Monte carlo d’un
écoulement de gaz dans des microcanaux 18 ème congrès français
[7] Jinliang xu et chuangxin zhao., two-simulation numerical simulation of shock waves in micro convergent-
divergent nozzles ,international journal of heat and mass transfer [Link] 2434-2438, (2007)
[8] Riadh ben hamouda; notions de mécanique des fluides, cours et exercices corrigés, centre de publication
universitaire, tunis (2008)
[9] Alain drotz marc a. Habisreutinger ; mécanique des fluides compressibles, ecole polytechnique fédérale de
lausanne laboratoire d’ingénierie numérique lausanne, (2007)
[10] Wolf gregory,2005 the n2k consortium. Ii. A transiting hot saturn around hd 149026 with a large dense core
1,2 arxiv: astro-ph/0507009v1 30 juin 2005
[11] Sébastien Candel, Cours mécanique des fluides, Livre papier, 2001 (2eme édition).
[12] Thompson, joe f . Proceedings of a symposium on the numerical generation of curvilinear coordinate
systems and their use in the numerical solution of partial differential equations, nashville, tennessee, defense
technical information center, held april 1982.
[13] Aoudjera farida. Calcul d’un ecoulement dans une tuyere miniaturisee, these de majister, universite el-hadj
lakhdar-batna.2011.