République Algérienne Démocratique Populaire

Ministère de l’Enseignement Supérieur et Recherche

Scientifique

Université Saad Dahleb Blida-1

Institut d’Aéronautique et Etudes Spatiale
Option : Propulsion Avion
Mémoire fin d’étude en vue d’obtention du diplôme Master

Thème :

SIMULATION NUMERIQUE D’UNE FLAMME

TURBULENTE DE DIFFUSION DE COMBUSTION
D’HYDROGENE DANS L’AIR

Réaliser par :
BOUCHACHOUA Nouh
CHERFI Malik

Promoteur : Co-promoteur :
Mr. NECHE Ahmed Mr. LAZEB Sbaa

2020/2021
 Remerciment

Je tiens à remercier le bon Dieu (Allah) le tout

puissant de m’avoir attribué la faveur de réussir
mes études.

Je tiens à remercier mon promoteur Mr.

NECHE Ahmed pour son suivi et encore plus
Mr. LAZAB Sbaa pour son aide et son suivi
immense tout le parcours préparent ce travail.

Sans oublier, je tiens à remercier tous les

enseignants de l’IAES qui font partie de ma
formation et toute personne sert à faciliter les
choses au niveau de cet institut.

Enfin je tiens à remercier toute ma famille, mes

chers amis, mes collègues
 Résumé
L’objectif principal de cette étude est la simulation du comportement de
l’écoulement réactif de la flamme turbulente de diffusion dans le cas de la combustion
d’hydrogéné dans l’air, ainsi apporter une contribution a l’analyse de la structure de la flamme
, et connaitre pour une pression et une température de gaz frais données le comportement des
paramétrés thermodynamiques de combustion . Notre étude est basée sur des simulations dans
le cas bidimensionnel en utilisant le logiciel Ansys Fluent qui résout les équations de base de
combustion turbulente ou le modèle de turbulence K-e est choisi. Les résultats obtenus sont
ensuit comparés au résultats scientifique spécialisé.

Mots clés : Combustion turbulente, Flamme de diffusion, modèl K-e , simulation

numérique .

Abstract
The main objective of this study is simulation composting of the reactive flow of the
turbulent diffusion flame in the case of the combustion of hydrogen in air, this making a
contribution to the analysis of the structure of the flame, and to know, for a given fresh gas
pressure and temperature, the behaviour thermodynamic combustion parameters. Our study is
based on simulations in the two-dimensional case using Ansys Fluent software which solves
the basic turbulent combustion equations where the K-e turbulence model is chosen. Then the
results obtained are compared to the specialized scientific results.

Keywords: Turbulent combustion, Diffusion flame, K-e model, numerical simulation.

 ‫ملخص‬

‫الهدف الرئيسي من هذه الدراسة هو محاكاة سماد التدفق التفاعلي للهب االنتشار المضطرب في حالة احتراق‬

‫الهيدروجين في الهواء‪ ،‬مما يساهم في تحليل هيكل اللهب‪ ،‬ومعرفة بالنسبة لضغط ودرجة حرارة غاز طازج معينين‬

‫سلوك معامالت االحتراق الديناميكي الحراري‪ .‬تعتمد دراستنا على عمليات المحاكاة في الحالة ثنائية األبعاد باستخدام‬

‫برنامج انسيس فلونت الذي يحل معادالت االحتراق المضطرب األساسية أين تم اختيار نموذج اضطراب‬

‫‪ .‬ك ابسلون ثم يتم مقارنة النتائج المتحصل عليها بنتائج علمية متخصصة ‪.‬‬

‫الكلمات المفتاحية احتراق مضطرب لهب ذو متفاعالت منفصلة نموذج اضطراب ك ابسيلون محاكاة رقمية‬

‫ثم يتم مقارنة النتائج التي تم الحصول عليها بالنتائج العلمية المتخصصة‬
 Sommaire

Introduction .......................................................................................................................... 01

Chapitre 1 GENERALITE SUR LA THERMODYNAMIQUE

Introduction : ........................................................................................................................ 02
1.1 Équation d'état ................................................................................................................ 02
1.2 Mélanges de gaz parfaits ................................................................................................. 02
1.3 Stœchiométrie ................................................................................................................. 03
1.4.1 La loi zéro de la thermodynamique .............................................................................. 04
1.4.2 La première loi de la thermodynamique ....................................................................... 04
1.4.3 La deuxième loi de la thermodynamique ...................................................................... 05
[Link] Transformation réversible ......................................................................................... 06
[Link] Transformation irréversible ....................................................................................... 06
1.4.5 Troisième loi de la thermodynamique .......................................................................... 06
1.5 Enthalpie ........................................................................................................................ 07
1.6 Fonction Gibbs ............................................................................................................... 07
1.7 Constante d’équilibre ...................................................................................................... 07
1.8 Enthalpie standardisée et enthalpie de formation ............................................................. 08

Chapitre 2 Généralité sur la combustion

Introduction ......................................................................................................................... 09
2.1 Types de combustion ...................................................................................................... 10
2.2 Définition de la flamme ................................................................................................. 10
2.3 Le Mode d’introduction des réactifs ............................................................................... 11
Le régime d’écoulement des gaz dans le milieu réactionnel .................................................. 12
2.3.1 Flamme pré mélangée Laminaire ................................................................................. 13
[Link] Les Relations de Rankine-Hugoniot .......................................................................... 15
[Link] Lignes de Rayleigh ................................................................................................... 15
[Link] Ligne de Hugoniot .................................................................................................... 16
[Link] Ondes de détonation et de déflagration ...................................................................... 17
2.3.2 Flamme Non pré mélangée Laminaire .......................................................................... 19
2.3.3Les flammes turbulentes non-prémélangées…………………………………………….21
2.4 Fraction du mélange ....................................................................................................... 22
2.5 Équation de conservation pour la fraction de mélange ..................................................... 23

Chapitre 3 Modélisation de la combustion turbulente

3.1 Introduction .................................................................................................................... 24
3.2 Echelles de longueur et de temps en turbulence............................................................... 24
3.3 Les Modèles de turbulence.............................................................................................. 25
3.3.2 Equations Fondamentales............................................................................................. 25
3.3.3 La Viscosité D’Eddy .................................................................................................... 27
3.3.4 L’équation de transport de l’énergie cinétique .............................................................. 28
3.3.5 Le Modèle K-Epsilon ................................................................................................... 30
3.3.6 Le Modèle K-Omega ................................................................................................... 31
3.3.7 Le Modèle K-Omega SST ........................................................................................... 31
Les Avantages et les Inconvénients de chaque modèle .......................................................... 38

Chapitre 4 résolution numérique et résultats

4.1 Introduction .................................................................................................................... 40
4.2 PRESENTATION DE CODE DE CALCUL FLUENT ................................................... 40
4.2.1 Architecture Du Logiciel.............................................................................................. 40
4.3.1 Création de la géométrie .............................................................................................. 41
4.3.2 Génération de maillage ............................................................................................... 41
4.3.3 Critères de qualités de maillage ................................................................................... 42
4.3.4 Paramètres des données de calcul ................................................................................. 43
[Link] Choix de modelé de la turbulence.............................................................................. 45
4.3.4.2Choix de la méthode de solution ................................................................................ 46
4.4.1 Contour de température ................................................................................................ 47
4.4.2 Répartition de la température le long de l'axe ............................................................... 48
4.4.3 Contour énergie cinétique turbulente ............................................................................ 49
4.4.4 Fraction massique d’O2 ............................................................................................... 49
4.4.5 Fraction massique de H2 .............................................................................................. 50
4.4.6 Variation de fraction de mélange .................................................................................. 50
Conclusion .......................................................................................................................... 51
 Liste des figures
Chapitre 1 : généralité sur la thermodynamique
Figure 1.1 Le premier principe de la thermodynamique……………………………………...5

Chapitre 2 : généralité sur la combustion

Figure 2.1 Variation de température, combustible et de l’oxydant au cours d’une
Combustion…………………………………………………………………………………….9
Figure 2.2 Schéma d’une flamme pré-mélangée et non pré-mélangée…………………..........11
Figure 2.3 Structure de la flamme ………………………………………………………........12
Figure 2.4 Les différents régimes de l’écoulement……………………………………............13
Figure 2.5 Structure d’une flamme laminaire pré-mélangée…………………………....…….14
Figure 2.6 Courbe d’Hugoniot-Rankine pression/volume massique…………………...…….17
Figure 2.7 Comparaison de la dynamique des gaz de détonation et de déflagration……..........18
Figure 2.8 La forme de la flamme de diffusion…………………………………………..........19
Figure 2.9 Structure d’une flamme laminaire non pré-mélangée…………………….…….…19
Figure 2.10 représentation schématique des différents régimes de combustion turbulente
non pré mélangée……………………………………………………………………………...22

Chapitre 3 : modélisation de la combustion turbulente

Figure 3.1 La variation de la fonction de mélange par rapport a la distance…………………...33

Figure 3.2 La variation de la fonction de mélange pour différents modèles…………………..33
Figure 3.3 la variation de l’argument par rapport à la distance ……………………………….34
Figure 3.4 La variation de la vitesse de l’écoulement ………………………………………...35
Figure 3.5 comparaison entre le modèle linéaire et non linéaire ……………………………...36
Figure 3.6 La variation de la vitesse par rapport a y+…………………………………………36
Figure 3.7 la décomposition de la fonction de vitesse…………………………………..…….37
Figure 3.8 comparaison entre la fonction générée par DNS et Spalding…………………...…38

Chapitre 4 : résolution numérique et résultats

Figure 4.1 Schéma du domaine de calcul…………………………………………………….41
Figure 4.2 le maillage du domaine d’étude……………………………………………………42
Figure 4.3 Critères « skewness » et “orthogonal quality »……………………………...……..42
Figure 4.4 : vérification de Critères « orthogonal quality » dans notre maillage……….……..43
Figure 4.5 fenêtre fluent setup………………………………………………………………...44
Figure 4.6 fenêtre fluent hypothese………………………………………………………...…45
Figure 4.7 fenêtre fluent modèle de turbulence………………………………………….……45
Figure 4.8 fenêtre fluent species model…………………………………………………….…46
Figure 4.9 choix de la méthode de solution .……… ………………………………….........…46
Figure 4.10 contour de température…………………………………………………………..47
Figure 4.11 répartition de la température le long de l’axe………………………………….…48
Figure 4.12 Energie cinétique turbulente……………………………………………………..49
Figure 4.13 variation de la fraction massique d’O2 …………………………………………...49
Figure 4.14 variation de la fraction massique de H2…………………………………………50
Figure 4.15 variation de fraction de mélange…………………………………………………50
Introduction

La combustion est un phénomène extrêmement répandu dans la nature, et c’est grâce à elle, en
majeur partie, que l’activité humaine a pu se développer et se développe encore. Dans
l’antiquité et jusqu’au moyen âge, le feu était considéré comme l’un des quatre éléments de
base dont l’univers entier était constitué. C’est-à-dire que l’étude et la compréhension de ce qui
concerne la combustion et le feu a un intérêt tout à fait important.

Les systèmes industriels impliquant des phénomènes de combustion (fours, moteurs

automobiles ou aéronautiques, turbines à gaz...) sont soumis à des contraintes de plus en plus
importantes, tant sur le plan économique (réduction des coûts, amélioration des
performances,...) que sur le plan environnemental (réduction des émissions de polluants, des
émissions sonores,. . .). L’ensemble de ces considérations motive de nombreux travaux de
recherche lies à la combustion. En effet, la compréhension, la modélisation et éventuellement
le contrôle de la structure de flamme permettent non seulement l’amélioration des systèmes
actuels mais aussi le développement de nouvelles technologies performantes et propres.

L'utilisation de l'hydrogène dans les turbines à gaz stationnaires est en train de devenir l'une des
approches les plus efficaces pour réduire les émissions de carbone provenant de la production
d'électricité à grande échelle. Cela s'applique également, dans une synergie pratique, aux
systèmes de production d'électricité qui peuvent utiliser un flux constant d'hydrogène provenant
de la refonte à grande échelle du gaz naturel avec captage et stockage du carbone (CSC) [1] ou,
en exploitant l'énergie excédentaire provenant de ressources d'énergie renouvelable non
distribuables (éolien et solaire), un flux instable d'hydrogène produit à partir de l'électrolyse de
l'eau couplée à des solutions de stockage d'énergie à grande échelle (énergie-H-énergie).

Dans un premier temps, le premier chapitre est destiné à la présentation générale des relations
de la thermodynamique.

Dans le deuxième chapitre, on adresse les caractérisations de la combustion et de la flamme.

Dans le troisième chapitre, on établisse les modèles de turbulence nécessaires pour nos
simulations.

1
 CHAPITRE 1

GENERALITE SUR LA
THERMODYNAMIQUE
Introduction :

La thermodynamique est une branche de la physique qui traite de la chaleur, du travail et de la

température, et de leur relation avec l'énergie. La thermodynamique a quatre lois qui décrivent
toutes les interactions entre ces quantités, qui véhiculent une description quantitative à l'aide de
quantités physiques macroscopiques mesurables.

La thermodynamique s'applique à une grande variété de sujets en sciences, en particulier la

chimie physique, le génie chimique et le génie mécanique, mais aussi dans d'autres domaines
complexes tels que la météorologie. [1]

1.1 Équation d'état

Une équation d'état décrit en une seule équation la relation entre la pression, P, la température,
T et le volume V (ou volume spécifique v) d'une substance. Pour le comportement des gaz
parfaits, c'est-à-dire un gaz modélisé en négligeant les forces intermoléculaires et le volume des
molécules, les formes équivalentes suivantes de l'équation d'état s'appliquent [1] :

𝑃𝑉 = 𝑁𝑅𝑢 𝑇 (1.1)

𝑃𝑉 = 𝑚𝑅𝑇 (1.2)

𝑃𝑣 = 𝑅𝑇 (1.3)

𝑃 = 𝜌𝑅𝑇 (1.4)

1.2 Mélanges de gaz parfaits

Considérons un mélange de gaz à plusieurs composants composé de 𝑁1 moles d'espèce 1, 𝑁2

moles d'espèce 2, etc. La fraction molaire d'espèces 𝝌𝒊, 𝝌𝒊 est définie comme la fraction du
nombre total de moles dans le système qui sont des espèces 𝑖 [1] :

𝑁𝑖 𝑁𝑖
𝜒𝑖 = = (1.5)
𝑁1 +𝑁2 +⋯+𝑁𝑖 +⋯ 𝑁𝑡𝑜𝑡

2
1.3 Stœchiométrie

La quantité stœchiométrique est la quantité exacte d'oxydant nécessaire pour brûler une quantité
spécifique de carburant. Si plus d'une quantité stœchiométrique de comburant est fournie, le
mélange est dit pauvre en carburant, tandis que fournir moins que le comburant
stœchiométrique conduit à un mélange riche ou riche en carburant. Le rapport stœchiométrique
du combustible oxydant (généralement de l'air) est déterminé en équilibrant simplement les
atomes des deux côtés de l'équation de réaction chimique, en supposant une réaction parfaite.
Pour un combustible hydrocarboné donné par Cx Hy, la relation stœchiométrique peut être
exprimée par [2] :

𝐶𝑥 𝐻𝑦 + 𝑎(𝑂2 + 3.76𝑁2 ) → 𝑥𝐶𝑂2 + (𝑦⁄2)𝐻2 𝑂 + 3.76𝑎𝑁2

Ou :

𝑎 = 𝑥 + 𝑦 ⁄4

Le rapport air-carburant stœchiométrique est :

𝑚𝑎𝑖𝑟
(𝐴⁄𝐹 )𝑆𝑡𝑜𝑖𝑐 = ( ) (1.6)
𝑚𝑓𝑢𝑒𝑙

Le rapport d'équivalence

Le rapport d'équivalence, Φ, est utilisé pour indiquer quantitativement si un mélange carburant-

comburant est riche, pauvre ou stœchiométrique.

(𝐴⁄𝐹 )𝑆𝑡𝑜𝑖𝑐 (𝐹 ⁄𝐴)

𝜱= (𝐴⁄𝐹 )
= (𝐹⁄ (1.7)
𝐴)𝑆𝑡𝑜𝑖𝑐

Si Φ > 1, le mélange est riche en carburant.

Si Φ < 1, le mélange est pauvre en carburant.

Pour un mélange stœchiométrique, Φ est égal à l'unité.

3
Dans de nombreuses applications de combustion, le rapport d'équivalence est le facteur le plus
important pour déterminer les performances d'un système.

D'autres paramètres fréquemment utilisés pour définir la stœchiométrie relative sont le

pourcentage d'air stœchiométrique, qui se rapporte au rapport d'équivalence comme :

100%
% 𝑠𝑡𝑜𝑖𝑐ℎ𝑖𝑜𝑚𝑒𝑡𝑟𝑖𝑐 𝑎𝑖𝑟 =
𝛷

Pourcentage d'excès d'air :

(1−𝛷)
%𝑒𝑥𝑐𝑒𝑠𝑠 𝑎𝑖𝑟 = 100% (1.8)
𝛷

1.4 Les lois de la thermodynamique

1.4.1 La loi zéro de la thermodynamique

Si A est en équilibre thermique avec B et B est en équilibre thermique avec C, alors C sera en
équilibre thermique avec A.

La loi zéro implique que tout comme la pression est une propriété physique qui nous permet
d'anticiper quand les systèmes seront en équilibre mécanique lorsqu'ils seront réunis
indépendamment de leur composition et de leur taille, il existe une propriété qui nous permet
d'anticiper quand deux systèmes seront en équilibre thermique quelles que soient leur
composition et leur taille : nous appelons cette propriété universelle la température. [3]

1.4.2 La première loi de la thermodynamique

La première loi de la thermodynamique est une extension de la loi de conservation de l'énergie,

selon laquelle l'énergie ne peut être ni créée ni détruite. C'est-à-dire quelle que soit la quantité
d'énergie qu'il y avait au début de l'univers, il y aura donc cette quantité à la fin. [3]

Pour une masse fixe, c'est-à-dire un système, la conservation de l'énergie s'exprime pour un
changement fini entre deux états, 1 et 2, par :

4
𝑄12 − 𝑊12 = ∆𝑈1−2 (1.9)

𝑄12 : Chaleur ajoutée au système en passant de l'état 1 à l'état 2

𝑊12 : Travail effectué par le système sur l'environnement en passant de l'état 1 à l'état 2

∆𝑈1−2 : Variation de l'énergie totale du système en passant de l'état 1 à l'état 2

Fig 1.1 : Le premier principe de la thermodynamique

1.4.3 La deuxième loi de la thermodynamique

La deuxième loi de la thermodynamique stipule que l'entropie totale d'un système isolé ne peut
jamais diminuer avec le temps, et est constante si et seulement si tous les processus sont
réversibles.

Les systèmes de la vie réelle sont soumis à des frottements, ce qui en fait un processus
unidirectionnel, ce qui signifie que l'entropie augmente toujours. [3]

1.4.4 Réversibilité

En thermodynamique, la réversibilité est une caractéristique de certains processus (passages

d'un système d'un état initial à un état final spontanément ou en raison d'interactions avec

5
d'autres systèmes) qui est réversible, et restaure le système à son état initial, sans laisser d'effets
nets dans l'un des systèmes impliqués. [4]

[Link] Transformation réversible

𝑄𝑟𝑒𝑣
∆𝑆𝑠𝑦𝑠 = (1.10)
𝑇

[Link] Transformation irréversible

𝑄𝑟𝑒𝑣
∆𝑆𝑠𝑦𝑠 = + ∆𝑆𝑐𝑟é𝑒𝑟 (1.11)
𝑇

1.4.5 Troisième loi de la thermodynamique

La troisième loi de la thermodynamique s'énonce comme suit, concernant les propriétés des
systèmes fermés en équilibre thermodynamique :

L'entropie d'un système approche une valeur constante lorsque sa température approche le zéro
absolu. [4]

Cette valeur constante est indépendante de tout paramètre caractérisant le système fermé. Au
zéro absolu, le système sera dans un état avec le minimum d'énergie possible. L'entropie est liée
au nombre de micro-états accessibles, et il existe généralement un état unique (appelé état
fondamental) avec une énergie minimale.

𝑆 − 𝑆0 = 𝐾𝐵 ln 𝛺 (1.12)

𝐾𝐵 : La constante de Boltzmann.

𝛺 : Le nombre de micro états cohérents avec la configuration macroscopique. Le comptage des

états se fait à partir de l'état de référence du zéro absolu, qui correspond à l'entropie de S 0.

6
1.5 Enthalpie

Il est défini comme la somme de l'énergie interne du système et le produit de sa pression et de

son volume, dérivant sa formule peut être fait à partir de la première loi et en utilisant l'équation
de définition du travail effectué par le système à une pression constante comme suit [5] :

𝐻 (𝑈, 𝑃, 𝑉 ) = 𝑈 + 𝑃𝑉 (1.13)

1.6 Fonction Gibbs :

L'utilisation d'un système isolé de masse fixe et d'énergie et de volume fixes n'est pas
particulièrement utile pour de nombreux problèmes typiques impliquant l'équilibre chimique.
Il sera toujours nécessaire de calculer la composition d'un mélange à une température, une
pression et une stœchiométrie données. [1]

1.7 Constante d’équilibre :

Pour le système général, où :

𝑎𝐴 + 𝑏𝐵+. . → 𝑒𝐸 + 𝑓𝐹 + ⋯

La variation du nombre de moles de chaque espèce est directement proportionnelle à son

coefficient stœchiométrique :

𝑑𝑁𝐴 = −𝜅𝑎

𝑑𝑁𝐵 = −𝜅𝑏

𝑑𝑁𝐸 = +𝜅𝑒

𝑑𝑁𝐹 = +𝜅𝑓

L'argument du logarithme népérien est défini comme la constante d'équilibre :

(𝑃 ⁄𝑃 0)𝑒 (𝑃𝐹 ⁄𝑃 0 )𝑓
𝑲𝒑 = (𝑃 𝐸⁄ (1.14)
𝐴 𝑃 0)𝑎 (𝑃𝐵 ⁄𝑃 0)𝑏

Avec ces définitions, notre énoncé de l'équilibre chimique à température et pression constantes
est donné par :
7
 ∆𝑮𝟎𝑻 = −𝑅𝑢 𝑇 ln 𝐾𝑝 (1.15)

Avec:

𝑲𝒑 = exp(− ∆𝐺𝑇0 ⁄𝑅𝑢 𝑇) (1.16)

1.8 Enthalpie standardisée et enthalpie de formation

Pour toute espèce chimique, l'enthalpie standardisée est définie comme étant la somme d'une
enthalpie qui prend en compte l'énergie associée aux liaisons chimiques formant les espèces
elles-mêmes, appelée enthalpie de formation ℎ𝑓 . Et une enthalpie qui n'est associée qu'à la
température, terme décrivant le changement de température d'un état initial à un autre, appelé
changement d'enthalpie sensible ∆ℎ𝑠 . Ainsi, nous pouvons écrire l'enthalpie molaire normalisée
pour les espèces comme [5] :

ℎ̅𝑖 (𝑇) = ℎ̅𝑓,𝑖

0
(𝑇𝑟𝑒𝑓 ) + ∆ℎ̅𝑠,𝑖 (𝑇) (1.17)

ℎ̅𝑖 (𝑇) : Enthalpie normalisée à la température T

ℎ̅𝑓,𝑖
0
(𝑇𝑟𝑒𝑓 ): Enthalpie de formation à l'état de référence standard (𝑇𝑟𝑒𝑓 , 𝑃0 )

∆ℎ̅𝑠,𝑖 (𝑇): Changement d'enthalpie sensible en passant de 𝑇𝑟𝑒𝑓 à T

L'enthalpie standard de formation est une mesure de l'énergie libérée ou consommée lorsqu'une
mole d'une substance est créée dans des conditions standard à partir de ses éléments purs.

Ainsi, les enthalpies de formation ont une interprétation physique comme le changement net
d'enthalpie associé à la rupture des liaisons chimiques des éléments à l'état standard et à la
formation de nouvelles liaisons pour créer un autre composé d'intérêt.

L'enthalpie sensible est une mesure de l'énergie nécessaire pour passer d'une température de
référence à une température finale. [5]

8
 Chapitre 2

GENERALITE SUR LA
COMBUSTION
Introduction :

La combustion est une réaction chimique, globalement exothermique, au début lente mais
pouvant devenir ensuite rapide et même violente, avec émission des rayonnements et élévation
de température qui peut se traduire par un dégagement de chaleur. Elle se produit dans un
système « ouvert » ou dans un système « fermé », autrement dit avec ou sans échange de matière
et d’énergie entre le système et le milieu extérieur. Les réactions chimiques ne peuvent se
produisent que si des espèces chimiques, atomes, ions, radicaux, molécules, commencent par
se rencontrer ou du moins se trouver assez près pour se modifier mutuellement. La réaction de
combustion est globalement une réaction d’oxydo-réduction où l’oxydant est appelé comburant
(oxygène pur ou dilué ozone, chlore, nitrates, etc.) et le réducteur est appelé combustible (H2,
CO, hydrocarbure, etc.) se déroulant suivant un mécanisme complexe, comprenant un grand
nombre de réaction dites «élémentaires », les unes contre les autres libérant de la chaleur, le
dégagement de chaleur prenant le pas sur la consommation. Si cette réaction oxydait de façon
maximale le combustible, on pourrait écrire, globalement :

a carburant + b combustible c CO2+ d H2O +…

Cette écriture a pour objet d’indiquer simplement quels sont les produits initiaux et finaux du
processus, sans plus. On l’appelle la réaction de bilan stœchiométrique. Elle n’indique pas du
tout ce qui se passe au cours de la combustion. [5]

Fig 2.1 : Variation de température, combustible et de l’oxydant au cours d’une combustion

9
2.1 Types de combustion

Il existe 5 principaux types de combustion classés comme suit : [6]

 Combustion complète : lorsque l'air est suffisant pour épuiser tous les combustibles,
libérant 𝐻2 𝑂 et 𝐶𝑂2
 Combustion incomplète : lorsque l'air n'est pas suffisant pour brûler tous les
combustibles, dégageant 𝐻2 𝑂, 𝐶𝑂2 et CO
 Combustion lente : forme de combustion à basse température et sans flamme, entretenue
par la chaleur dégagée lorsque l'oxygène attaque directement la surface d'un
combustible en phase condensée.
 Combustion spontanée : se produit par auto-échauffement accélère rapidement jusqu'à
des températures élevées puis s'enflamme.
 Combustion rapide : comme la déflagration ou la combustion par détonation.

2.2 Définition de la flamme :

Une flamme est une région localisée de réaction chimique intense. Cette région de réaction
chimique intense est caractérisée par l'émission de chaleur et généralement aussi de lumière.
[7]

10
2.3 Le Mode d’introduction des réactifs :

Les flammes sont classées en deux grandes catégories : les flammes pré mélangées et non pré
mélangées.

Si le combustible et le comburant sont initialement séparés puis transportés l'un vers l'autre,
après allumage, la flamme se place toujours à une position où le combustible et le comburant
sont présents dans des proportions stœchiométriques. De telles flammes sont appelées flammes
non pré mélangées.

Si les réactifs sont mélangés a priori, après allumage, le front de réaction se propage dans
l'imbrûlé avec une vitesse définie. De telles flammes sont appelées flammes pré mélangées. [6]

Fig 2.2 : Schéma d’une flamme pré-mélangée et non pré-mélangée

11
 Fig 2.3 : Structure de la flamme

Type de Flamme Avantages Inconvénients

Pré-mélangé Peu polluant Instabilité de combustion
Stockage du mélange
Non Pré-mélangé Stockage du mélange en Production d’oxydes d’azote
dehors de la limite de
flamabilité

Le régime d’écoulement des gaz dans le milieu réactionnel

L’écoulement des gaz est défini différemment suivant son régime est que l’on va distinguer :

• les écoulements laminaires, qui se produisent à petits nombreux de Reynolds,

• les écoulements turbulents, qui se produisent à des Reynolds plus élevés. La turbulence est
due à l’instabilité de l’écoulement laminaire devant les petites perturbations expérimentales
possibles, et que l’état turbulent est la situation non stationnaire vers laquelle tend l’écoulement
instable.

La distinction se fait donc à l’aide du nombre de Reynolds 𝑅𝑒=𝜌𝑈𝑑μ avec 𝑈 et 𝑑μ une vitesse
et une dimension caractéristique du problème étudié. La transition d’un régime laminaire à un
régime turbulent se fait graduellement sur une plage critique de nombre de Reynolds, dont les
valeurs dépendent de chaque problème.

12
Pour un nombre de Reynolds inférieur à 2200, l’écoulement est dit laminaire. Supérieur à 3000,
le régime est turbulent et pour des valeurs comprises entre 2200 et 3000, l’écoulement est dit
transitoire. [9]

Fig 2.4 : Les différents régimes de l’écoulement

2.3.1 Flamme pré mélangée Laminaire

La flamme pré mélangée se propage à travers le mélange non brûlé à une vitesse définie. Ainsi,
I peut être considéré comme un front séparant le mélange brûlé et non brûlé. A ce niveau de
résolution, la flamme est une surface d'épaisseur infinitésimale à travers laquelle il y a un saut
de propriétés de mélange de la température des gaz non brûlés Tu à la température des gaz brûlés
Tb et de la concentration de combustible dans la région non brûlée yu à la concentration de
combustible dans la région brûlée yb =0. [6]

D’après la théorie de Mallard et Le Châtelier (1883), la zone intermédiaire entre les gaz frais et
les gaz brûlés peut être séparée en deux zones : la zone de préchauffage et la zone de réaction.

La zone de préchauffage : zone d’épaisseur 𝛿𝑝 dans laquelle la diffusion et la convection sont

les processus prédominants et où les réactions chimiques sont par conséquent négligeables.

La zone de réaction : zone d’épaisseur 𝛿𝑟 très mince dans laquelle se fait le dégagement de
chaleur dû aux réactions chimiques. Dans cette zone, les phénomènes de conduction sont alors
négligeables par rapport à ceux de diffusion.

13
 Fig 2.5 : Structure d’une flamme laminaire pré-mélangée

La définition de l’épaisseur de flamme laminaire (appelée épaisseur de Zeldovich) est obtenue

en supposant que la totalité de l’énergie produite par la combustion sert à élever la
température des gaz frais présents dans la zone de préchauffage. [8]

𝛬
𝛿𝐿 = (2.1)
𝜌𝑢 𝐶𝑝 𝑠𝐿0

𝛬, 𝜌𝑢 et 𝐶𝑝 sont respectivement la conductivité thermique, la masse volumiques et la capacité

calorifique massique à pression constante du mélange réactif. 𝑠𝐿0 Est la vitesse de combustion
laminaire d’une flamme plane, adiabatique et non-étirée.

Le front de flamme laminaire ne correspondant pas à une discontinuité brusque, l’épaisseur de

flamme laminaire peut être considérée comme la somme de l’épaisseur de la zone préchauffage
et de l’épaisseur de la zone de réaction

On peut alors définir un temps caractéristique de la combustion 𝜋𝑐 comme le rapport𝛿𝐿 /𝑠𝐿 .

Celui-ci intègre le temps de préchauffage et le temps de réaction. [6]

14
[Link] Les Relations de Rankine-Hugoniot

Une flamme pré mélangée peut être considérée comme une onde se propageant à travers un
mélange réactif et séparant la région non brûlée de la région brûlée.

En intégrant les équations de transport à travers la flamme sur une échelle de longueur
hydrodynamique, les relations de conservation suivantes sont obtenues. [6]

Masse :

𝜌𝑢 𝑉𝑢 = 𝜌𝑏 𝑉𝑏 = 𝑓 (2.2)

Moment :

𝜌𝑢 𝑉𝑢2 + 𝑝𝑢 = 𝜌𝑏 𝑉𝑏2 + 𝑝𝑏
(2.3)

Energie :

𝑉𝑢 ² 𝑉𝑏 ²
ℎ𝑢 + = ℎ𝑏 +
2 2

(2.4)

[Link] Lignes de Rayleigh

En réarrangeant les équations, on peut écrire :

1 1
𝑝𝑏 − 𝑝𝑢 = 𝜌𝑏 𝑉𝑏2 − 𝜌𝑢 𝑉𝑢2 = 𝑓² (𝜌 − 𝜌 ) = −𝑓²(𝑣𝑏 − 𝑣𝑢 ) (2.5)
𝑢 𝑏

Où :

𝑝̂ − 1 = −𝛾𝑀𝑢 ²(𝑣̂ − 1) (2.6)

Avec :

𝑝̂ = 𝑝𝑏 /𝑝𝑢 ; 𝑣̂ = 𝑣𝑏 /𝑣𝑢

Le nombre de Mach amont est défini comme :

15
 𝑉𝑢 ² 𝑉𝑢 ² 𝜌𝑢 𝑉𝑢 𝑣𝑢 𝑓²𝑣𝑢
𝑀𝑢 ² = = 𝛾𝑝𝑢 = = (2.7)
𝑐𝑢 ² 𝛾𝑝𝑢 𝛾𝑝𝑢
𝜌𝑢

𝑝̂−1
𝑀𝑢 ² = −
𝛾(𝑣̂−1)

(2.8)

(𝑝̂−1) 𝑣̂
𝑀𝑏 ² = − (2.9)
(𝑣̂−1) 𝑝̂

[Link] Ligne de Hugoniot

De l’équation de l’energie on peut écrire :

2𝑉 2
1 1 𝜌𝑢 𝑉𝑢2 𝜌𝑢 𝑏 1
ℎ𝑏 − ℎ𝑢 = 2 (𝑉𝑢2 − 𝑉𝑏2 ) = 2 ( − 2 ) = 2 𝑓²(𝑣𝑏 − 𝑣𝑢 )(𝑣𝑏 + 𝑣𝑢 ) (2.10)
𝜌𝑢 𝜌𝑢

En combinant l’équation de Rayleigh et l’équation précédente :

1
ℎ𝑏 − ℎ𝑢 = 2 (𝑝𝑏 − 𝑝𝑢 )(𝑣𝑏 + 𝑣𝑢 ) (2.11)

Pour le dégagement de chaleur qc par unité de masse, on peut écrire :

𝐶𝑝 𝑝 𝑝
ℎ𝑏 − ℎ𝑢 = −𝑞𝑐 + 𝐶𝑝 (𝑇𝑏 − 𝑇𝑢 ) = −𝑞𝑐 + (𝜌𝑏 − 𝜌𝑢 ) (2.12)
𝑅 𝑏 𝑢

En combinant les des équations précédentes, on peut écrire :

𝛾−1 𝛾−1 4𝛾 𝛾−1

(𝑝̂ + 𝛾+1) (𝑣̂ − 𝛾+1) = (𝛾+1)² + 2𝑞̂𝑐 (𝛾+1) (2.13)

𝑞𝑐
𝑞̂𝑐 =
𝑝𝑢 𝑣𝑢

16
 Fig 2.6 : Courbe d’Hugoniot-Rankine pression/volume massique

A cause de dégagement de chaleur au sein de milieu réagissant, la température est plus élevée
pour les produits que pour les réactifs, aussi la pression. La courbe 𝐻0 du milieu avant
combustion et la 𝐻 du milieu fluide après combustion. Du point 𝐼 de la première on peut par
combustion atteindre un point 𝐽, 𝐽′ou 𝐽" de la seconde. [7]

- Soit au-dessus du point B, donc avec (𝑝−𝑝0) positif et (1𝜌⁄−1𝜌0) négatif, c’est le cas de la
détonation.

- Soit au-dessus du point C, donc avec (𝑝−𝑝0) négatif et (1𝜌⁄−1𝜌0) positif, c’est le cas de la
déflagration.

[Link] Ondes de détonation et de déflagration

Les lignes de Rayleigh et les courbes d'Hugoniot sont obtenues à partir de la conservation
équations, toute solution réalisable doit satisfaire les deux. Ainsi, les solutions réalisables
réduire aux racines de l'ensemble des équations. Les résultats sont connus sous le nom de
relations Rankine-Hugoniot. [7]

17
Manipulations algébriques conduisant aux solutions des équations peuvent conduire à des
connaissances physiques importantes. Des manipulations algébriques considérables conduisent
aux racines suivantes des équations : [7]

1−𝑀 ²
𝑢 2(𝛾 2 −1) 𝑀𝑢 ²
𝑣̂ − 1 = (𝛾+1)𝑀 [1 ± √1 − ̂𝑐 ]
2 )² 𝑞 (2.14)
² 𝑢 𝛾 (1−𝑀𝑢

𝛾(1−𝑀𝑢 ²) 2(𝛾 2 −1) 𝑀𝑢 ²

𝑝̂ − 1 = [1 ± √1 − ̂𝑐 ]
2 )² 𝑞 (2.15)
𝛾+1 𝛾 (1−𝑀𝑢

A partir de ces équations, on voit que, pour chaque ensemble de paramètres, il peut y avoir deux
solutions.

2(𝛾 2 −1) 2
𝑀𝑢
Puisque 2) 2 𝑞̂𝑐 > 0 pour toutes les conditions, le terme sous la racine signe est
𝛾 (1−𝑀𝑢

toujours inférieur à 1. Par conséquent, le terme à l'intérieur de la parenthèse de boîte est toujours
positif. Ainsi, selon que 𝑀𝑢2 > 1 ou 𝑀𝑢2 < 1, on peut obtenir 𝑣̂ < 1ou 𝑣̂ > 1et 𝑝̂ < 1 ou 𝑝̂ >
1. La solution supersonique (𝑀𝑢 > 1) est connue sous le nom d'onde de détonation ; la solution
subsonique (𝑀𝑢 < 1) est appelée onde de déflagration.

Fig 2.7 : Comparaison de la dynamique des gaz de détonation et de déflagration

18
2.3.2 Flamme Non pré mélangée Laminaire :

Très généralement, dans les flammes de diffusion, les phénomènes de mélange sont beaucoup
plus lents que les phénomènes physico-chimiques qui accompagnent la [Link] en
résultent que le développement des flammes de diffusion est commandé essentiellement par
des facteurs aérodynamiques et l’aspect chimique lié à la nature des réactifs reste secondaire.
[9]

Ce type de flamme est connu sous le nom de flamme non pré mélangée. Comme le processus
de mélange est régi par la diffusion du carburant et de l'air, on l'appelle diffusion flamme.

Fig 2.8 : La forme de la flamme de diffusion

La hauteur de flamme est un paramètre important pour caractériser une flamme. La hauteur
peut être déterminée comme la distance le long de l'axe de la flamme

Fig 2.9 : Structure d’une flamme laminaire non pré-mélangée

19
On suppose une flamme de diffusion conique simple et que la diffusion n'est présente que dans
la direction perpendiculaire à l'écoulement, nous pouvons estimer la hauteur de la flamme grâce
à une analyse très simple. Le déplacement carré moyen d'une molécule (disons, l'oxygène) est
donné par [8] :

2
𝑦 = 2𝐷𝑡 (2.16)

Où, y est le déplacement moyen, D est le coefficient de diffusion et t est le temps requis. A
hauteur de flamme, les molécules d'oxygène sont diffusées jusqu'à la ligne médiane. Le
déplacement peut être pris comme le rayon de la buse de carburant, R. Le temps nécessaire à la
diffusion peut être écrit comme suit :

Le temps peut également être assimilé à :

𝐿𝐹
𝑡= (2.17)
𝑣𝑒

Où LF est la hauteur de flamme et ve est la vitesse moyenne du carburant à la sortie de la buse.

Nous pouvons écrire le débit volumique, Q, sous la forme :

𝑄 = 𝑉𝑒 𝜋𝑅² (2.18)

Et la hauteur de flamme peut être exprimée comme :

𝑣𝑒 𝑅²
𝐿𝐹 = = 𝑄/𝜋𝐷 (2.19)
2𝐷

Comme indiqué dans la section précédente, la définition d'un scalaire conservé est nécessaire.
Le scalaire conservé est une propriété scalaire, qui est conservée dans tout le champ
d'écoulement. Les scalaires conservés ne sont ni générés ni détruits dans le champ d'écoulement.
La masse peut être identifiée comme un tel scalaire qui ne peut pas être généré ou détruit dans
le champ d'écoulement. Il n'y a ni source ni puits pour la masse. Auparavant, nous avons dérivé
un certain nombre d'équations de conservation pour différentes variables, telles que la masse,
la quantité de mouvement, les espèces, l'énergie, etc. Rappelons que les équations de
conservation prennent une forme générique contenant principalement quatre termes : temporel,
advection, diffusion et source. Si nous pensons à une propriété scalaire φ, l'équation de
conservation peut s'écrire : [5]

𝜕𝜑
𝜌 𝜕𝑡 + 𝜌𝑉. ∇𝜑 = 𝛻. [𝜌𝛿 (𝛻𝜑)] + 𝑆 (2.20)

20
Où δ est le type de diffusivité associé au scalaire φ. Ici, les deux termes du côté gauche sont
respectivement temporels et advectifs. Le premier terme à droite est le terme de diffusion et le
dernier est le terme source. Pour les scalaires conservés, il ne devrait pas y avoir de terme source
dans l'équation de conservation ; à titre d'exemple, l'équation de conservation de la masse ne
contient pas de terme source non nul.

Pour les scalaires conservés, l'équation de conservation se présentera comme suit :

𝜕𝜑
𝜌 𝜕𝑡 + 𝜌𝑉. ∇𝜑 − 𝛻. [𝜌𝛿(𝛻𝜑)] = 0 (2.21)

Conceptuellement, nous pouvons identifier certains scalaires comme la masse, la masse de

carbone, la masse d'oxygène, etc. On peut aussi dériver des propriétés (par combinaison
appropriée de propriétés élémentaires) qui obéissent à la condition ci-dessus ou suivent le type
d'équation de conservation.

[Link] flammes turbulentes non-prémélangées

De manière similaire, il est possible de présenter les différents régimes de la combustion

turbulente de diffusion pouvant être rencontrés en pratique.

Cuenot et Poinsot ont tracé un diagramme en fonction du nombre de Damköhler de la flamme

et du nombre de Reynolds turbulent représentant les différents régimes de ce type de
combustion

On y distingue quatre domaines séparés par trois courbes correspondant des nombres de
Damköhler diffusifs constants ; la limite du régime de flammelettes par la droite (LFA pour
”Laminar Flamelet Assumption”) 𝐷𝑎 = 𝐷𝑎𝐿𝐹𝐴 et la limite d’extinction est représentée par la
droite 𝐷𝑎 = 𝐷𝑎𝑒𝑥𝑡 :

• Le premier domaine pour de faibles nombres de Reynolds (𝑅𝑒 < 1), la flamme est
laminaire.

Les trois domaines suivants sont pour des valeurs du nombre de Reynolds supérieur à l’unité.

• Le second régime correspond pour un Damköhler largement supérieur à 1 (𝐷𝑎 >

𝐷𝑎𝐿𝐹𝐴). Lorsque le temps chimique est largement inférieur au temps caractéristique de la
turbulence, la combustion peut être assimilée à un modèle de combustion infiniment rapide ou
régime des flamme lettes. On peut considérer la flamme turbulente comme étant constituée
d’une collection de flamme lettes.

21
• Le troisième domaine est compris entre les deux limites 𝐷𝑎𝐿𝐹𝐴 > 𝐷𝑎 > 𝐷𝑎𝑒𝑥𝑡. Les
temps caractéristiques de la chimie ne sont plus négligeables devant les temps de la
turbulence alors que les structures turbulentes agissent sur la structure de la flamme, il s’agit
du régime de chimie finie. Les effets instationnaires liés à la turbulence peuvent apparaître.

• Le quatrième domaine est au-delà de la limite d’extinction 𝐷𝑎 < 𝐷𝑎𝑒𝑥𝑡. Lorsque la

turbulence devient l’élément directeur, ils peuvent apparaître des extinctions locales de la
flamme.

Figure 2.10 : représentation schématique des différents régimes de combustion

turbulente non prémélangée .

2.4 Fraction du mélange :

La fraction de mélange est un scalaire conservé qui est très couramment utilisé pour la solution
de flammes non pré mélangées. Restreignons notre discussion à un système d'un seul flux
d'entrée de carburant pur. De même, nous supposons un seul flux d'oxydant pur. Nous prenons
également tous les produits de combustion comme une seule entité. La fraction de mélange, Z,

22
peut être définie comme le rapport de la masse de matière ayant son origine dans le flux de
carburant à la masse du mélange. [6]

1𝑘𝑔 𝑓𝑢𝑒𝑙 + 𝛾𝑘𝑔 𝑑 ′ 𝑎𝑖𝑟 → (1 + 𝛾)𝑘𝑔 𝑑𝑒 𝑝𝑟𝑜𝑑𝑢𝑖𝑡𝑠

Dans un petit volume, où le contenu est un mélange des trois espèces ci-dessus, si l'on veut
trouver la fraction de mélange, il aura des contributions des trois espèces :

1
𝑍 = (1)(𝑌𝐹 ) + (0)(𝑌0 ) + (𝛾+1)(𝑌𝑝𝑟 ) (2.22)

Où 𝑌𝐹 , 𝑌0 et 𝑌𝑝𝑟 sont les fractions massiques du combustible, du comburant et du produit,

respectivement. Dans notre système, la fraction de mélange dans le flux de carburant pur est de
1 et elle est de 0 dans le flux de comburant pur. Ainsi, la contribution des fractions massiques
de carburant et de comburant à la fraction de mélange doit être respectivement de 1 et 0. De
1
même, la fraction de masse de matière combustible par unité de masse de produits est de 1+𝛾.

Donc ça devrait être la contribution de la fraction massique des produits. La fraction de mélange
peut être déterminée comme suit :

1
𝑍 = 𝑌𝐹 + (𝛾+1) (𝑌𝑝𝑟 ) (2.23)

2.5 Équation de conservation pour la fraction de mélange :

Nous pouvons écrire les équations de conservation des espèces du combustible et du comburant
dans une dimension comme suit :

𝑑 𝑌
𝑝𝑟 𝑑 𝑑 𝑝𝑟𝑌
𝑚" 𝑑𝑥 (𝑌𝐹 + 𝛾+1 ) − 𝑑𝑥 [𝜌𝐷 𝑑𝑥 (𝑌𝐹 + 𝛾+1 )] = 𝑚"′𝐹 + 𝑚"′𝑃𝑟 /(𝛾 + 1) (2.24)

𝑚" : Le flux massique total

𝑚"′𝐹 ∶ Taux de production massique fuel

𝑚"′𝑃𝑟 : Taux de production massique produits

Donc :

𝑑𝑍 𝑑 𝑑𝑍
𝑚" 𝑑𝑥 − 𝑑𝑥 [𝜌𝐷 𝑑𝑥 ] = 0 (2.25)

23
 CHAPITRE 3

MODELISATION DE LA
COMBUSTION
TURBULENTE
3. Modélisation de la combustion turbulente

3.1 Introduction

L'étude des écoulements de fluides turbulents – comme dans toute autre étude scientifique –
nécessite une théorie quantitative traitable ou un modèle qui peut être utilisé pour évaluer
quantités qui présentent un intérêt et ont une pertinence pratique. Les caractéristiques inhérentes
des écoulements turbulents rendent cependant difficile l'élaboration d’une théorie ou modèle.

Le champ de vitesse est tridimensionnel, dépendant du temps et aléatoire. L'échelle des plus
grands mouvements turbulents est comparable à la caractéristique longueur de l'écoulement, et
par conséquent, ils sont considérablement affectés par la géométrie de l'écoulement.

En revanche, les plus petites échelles de longueur sont comparables aux échelles moléculaires.
Les échelles de temps associées sont également de nature aussi variée.

Des difficultés surviennent également en raison de la nature non linéaire du terme convectif et
le terme de gradient de pression dans l'équation de Navier-Stokes. Toute théorie ou modèle faire
face à la turbulence doit tenir compte de cette large gamme de longueur et échelles de temps et
devrait également être capable de gérer le caractère non linéaire des équations régissant.

Toutes les approches de modélisation de la turbulence qui ont été développées à ce jour sont
basées sur la résolution de l'équation de Navier-Stokes pour déterminer la vitesse d'écoulement
et les propriétés associées. [8]

3.2 Echelles de longueur et de temps en turbulence

L'une des caractéristiques de base de tout écoulement turbulent est l'existence d'une large
gamme d'échelles de longueur et de temps. Cette caractéristique est essentiellement ce qui
présente la difficulté de la modélisation numérique de la turbulence. Ainsi, il est
particulièrement important de comprendre les différentes échelles de durée et de temps
associées à l'écoulement turbulent.

Il y a, en général, quatre ensembles principaux d'échelles dans un écoulement turbulent. Ils sont
décrits dans les sous-sections suivantes. [8]

24
3.3 Les Modèles de turbulence

3.3.1Introduction:

La turbulence désigne l'état de l'écoulement d'un fluide dans lequel la vitesse présente en tout
point un caractère tourbillonnaire ,dont la taille, la localisation et l'orientation varient
constamment. Les équations de Navier-Stokes sont des équations aux dérivées partielles non
linéaires qui décrivent le mouvement des fluides. L'existence mathématique de solutions des
équations de Navier-Stokes n'est pas démontrée. Mais elles permettent souvent, par une
résolution approchée, de proposer une modélisation de nombreux phénomènes.

Reynolds a utilisé les équations de Navier-Stokes pour simplifier le problème de turbulence. Il

a utilisé les équations résultantes pour pouvoir modéliser l’écoulement turbulent afin de
déterminer les propriétés de ce dernier.

3.3.2 Equations Fondamentales

Equations de Navier et Stokes

Il est admis que le comportement de tout écoulement vérifiant l’hypothèse des milieux continus,
quelle que soit la nature du fluide (compressible ou non), et de l’écoulement (turbulent,
laminaire ou en transition), peut être représenté par les équations de Navier-Stokes qui
expriment la variation de la quantité de mouvement, auxquelles on y ajoute les équations de
conservation de la masse et de l’énergie. Ainsi, pour un écoulement de fluide visqueux,
compressible et conducteur de chaleur en négligeant les forces volumiques extérieures
(pesanteur, etc.), ces équations s’écrivent sous la forme instantanée suivante

Equation de continuité :

𝜕𝜌 𝜕𝑣
+ =0
𝜕𝑡 𝜕𝑥𝑗

Pour 𝜌𝑡 = 0

𝜕𝑣
=0
𝜕𝑥𝑗

En Remplace 𝑣 = 𝑣̅ + 𝑣′

25
 𝜕(𝑣̅ + 𝑣 ′ )
=0
𝜕𝑥𝑗

On fait la moyenne de l’équation précédente :

̅̅̅̅̅̅̅̅̅̅̅̅
𝜕(𝑣̅ + 𝑣 ′ )
=0
𝜕𝑥𝑗

L’équation devient :

𝜕𝑣̅
=0
𝜕𝑥𝑗

L’équation de conservation de quantité de mouvement :

𝜕𝑢 𝜕𝑢 𝜕𝑝 𝜕 𝜕𝑢
𝜌[ +𝑣 ]= − + (𝜇 )
𝜕𝑡 𝜕𝑥𝑗 𝜕𝑥𝑖 𝜕𝑥𝑗 𝜕𝑥𝑗

On remplace

𝑢 = 𝑢̅ + 𝑢′

𝑣 = 𝑣̅ + 𝑣′

𝑝 = 𝑝̅ + 𝑝′

𝜕(𝑢̅ + 𝑢′ ) 𝜕 𝜕(𝑝̅ + 𝑝′ ) 𝜕 𝜕(𝑢̅ + 𝑢′ )

𝜌[ + [(𝑣̅ + 𝑣 ′ )(𝑢̅ + 𝑢′ )]] = − + (𝜇 )
𝜕𝑡 𝜕𝑥𝑗 𝜕𝑥𝑖 𝜕𝑥𝑗 𝜕𝑥𝑗

On fait la moyenne de l’équation :

̅̅̅̅̅̅̅̅̅̅̅̅̅
𝜕(𝑢̅ + 𝑢′ ) ̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅
𝜕 ̅̅̅̅̅̅̅̅̅̅̅̅
𝜕(𝑝̅ + 𝑝′ ) ̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅
𝜕 𝜕(𝑢̅ + 𝑢′ )
𝜌[ + [(𝑣̅ + 𝑣 ′ )(𝑢̅ + 𝑢′ )]] = − + (𝜇 )
𝜕𝑡 𝜕𝑥𝑗 𝜕𝑥𝑖 𝜕𝑥𝑗 𝜕𝑥𝑗

̅̅̅̅̅̅̅̅̅̅
𝜕(𝑢̅ +𝑢 ′ ) ̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅
𝜕 ̅̅̅̅̅̅̅̅
𝜕(𝑝̅+𝑝′ ) ̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅
𝜕 𝜕(𝑢̅ +𝑢 ′ )
𝜌[ + 𝜕𝑥𝑗 [𝑢̅𝑣̅ + 𝑢̅𝑣′ + 𝑣̅ 𝑢′ + 𝑣′𝑢′]] = − 𝜕𝑥𝑖 + 𝜕𝑥𝑗 (𝜇 𝜕𝑥𝑗 )
𝜕𝑡

L’équation devient :

26
 𝜕𝑢̅ 𝜕𝑢𝑣
̅̅̅̅ 𝜕𝑝̅ 𝜕 𝜕𝑢̅ ̅̅̅̅̅
𝜌[ + ]= − + (𝜇 ) − 𝜌𝑢′𝑣′
𝜕𝑡 𝜕𝑥𝑗 𝜕𝑥𝑖 𝜕𝑥𝑖 𝜕𝑥𝑖

̅̅̅̅̅ est nommé le tenseur de Reynolds

Le terme 𝜌𝑢′𝑣′

Le Tenseur de Reynolds Peut être écrit sous la forme matricielle suivante :

̅̅̅̅̅
𝑢′ ²1 ̅̅̅̅̅̅̅̅
𝑢′1 𝑢′ 2 ̅̅̅̅̅̅̅̅
𝑢 ′1 𝑢 ′ 3
̅̅̅̅̅ = 𝜌 [̅̅̅̅̅̅̅̅
𝜌𝑢′𝑣′ 𝑢 ′ 2 𝑢 ′1 ̅̅̅̅̅
𝑢′ ² ̅̅̅̅̅̅̅̅
𝑢 ′ 𝑢′ ]
2 2 3
̅̅̅̅̅̅̅̅
𝑢 ′ 3 𝑢 ′1 ̅̅̅̅̅̅̅̅
′
𝑢 3𝑢 2′ ̅̅̅̅̅
𝑢′ ²3

Les équations 2 en (3D) sont : Les Equations RANS

Le nombre des inconnus est plus grand que le nombre des équations donc on a un problème de
fermeture. Pour fermer le système on a besoin des équations additionnelles

Les modèles de turbulence sont des équations additionnelles qu'on doit résoudre avec les
équations RANS pour trouver tous les inconnus du système

3.3.3 La Viscosité D’Eddy

̅̅̅̅̅̅
Dans le modèle de Eddy on doit changer le terme additionnel 𝜌𝑢 ′ 𝑣 ′ par l‘équation suivante :

̅̅̅̅̅
𝑢′ ²1 + ̅̅̅̅̅
𝑢′ ²2 + ̅̅̅̅̅
𝑢′ ²3 𝜕𝑢 𝜕𝑢
̅̅̅̅̅̅
𝜌𝑢 ′ 𝑣′ = 𝜌 + 𝜇𝑡 [ + ]
3 𝜕𝑥𝑖 𝜕𝑥𝑗

Cette égalité contient deux termes :

̅̅̅̅̅̅ ̅̅̅̅̅̅
𝑢 ′ ²1 +𝑢 ′² +𝑢̅̅̅̅̅̅
′²
2 3
𝜌 Isotrope
3

𝜕𝑢 𝜕𝑢
𝜇𝑡 [𝜕𝑥𝑖 + 𝜕𝑥𝑗 ] Anisotrope

On peut modifier le terme isotrope comme suit :

2 ̅̅̅̅̅
𝑢′ ²1 + ̅̅̅̅̅
𝑢′ ²2 + ̅̅̅̅̅
𝑢′ ²3 𝜕𝑢 𝜕𝑢
̅̅̅̅̅̅
𝜌𝑢 ′ 𝑣′ = 𝜌 + 𝜇𝑡 [ + ]
3 2 𝜕𝑥𝑖 𝜕𝑥𝑗

27
 2 𝜕𝑢 𝜕𝑢
̅̅̅̅̅̅
𝜌𝑢 ′𝑣′ = 𝜌𝑘𝛿𝑖𝑗 + 𝜇𝑡 [ + ]
3 𝜕𝑥𝑖 𝜕𝑥𝑗

Ou k est l‘énergie cinétique de l‘écoulement turbulent

Cette approche consiste à trouver la viscosité turbulente 𝜇𝑡

̅̅̅̅̅̅
En trouvant 𝜇𝑡 on peut remplacer le terme 𝜌𝑢 ′ 𝑣 ′ dans les équations RANS et le système va

devenir fermé. Les modèle k-Epsilon et k-Omega sont basés sur cette approche.

3.3.4 L’équation de transport de l’énergie cinétique

On peut trouver l‘équation de transport de l‘énergie cinétique de la manière suivante :

D'abord on fait la substruction de l‘équation (2) de l‘équation (a) :

𝜕(𝑢̅ + 𝑢′ ) 𝜕 𝜕(𝑝̅ + 𝑝′ ) 𝜕 𝜕(𝑢̅ + 𝑢′ )

𝜌[ + [𝑢̅𝑣̅ + 𝑢̅𝑣′ + 𝑣̅ 𝑢′ + 𝑣′𝑢′]] = − + (𝜇 )
𝜕𝑡 𝜕𝑥𝑗 𝜕𝑥𝑖 𝜕𝑥𝑗 𝜕𝑥𝑗

𝜕𝑢̅ 𝜕𝑢𝑣
̅̅̅̅ 𝜕𝑝̅ 𝜕 𝜕𝑢̅ ̅̅̅̅̅
𝜌[ + ]= − + (𝜇 ) − 𝜌𝑢′𝑣′
𝜕𝑡 𝜕𝑥𝑗 𝜕𝑥𝑖 𝜕𝑥𝑗 𝜕𝑥𝑗

On trouve l‘équation de conservation de quantité de mouvement en fonction des fluctuations :

𝜕𝑢′ 𝜕 ′ ′ ′ ′ ̅̅̅̅̅ 𝜕𝑝′ 𝜕 𝜕𝑢′

𝜌[ + [𝑢̅𝑣 + 𝑣̅ 𝑢 + 𝑣 𝑢 − 𝑢′𝑣′]] = − + (𝜇 )
𝜕𝑡 𝜕𝑥𝑗 𝜕𝑥𝑖 𝜕𝑥𝑗 𝜕𝑥𝑗

En multipliant par 𝑢′ on trouve :

𝜕𝑢′ 𝜕 𝜕 𝜕 𝜕
𝑢′ + 𝑢′ (𝑢̅𝑣 ′ ) + 𝑢′ (𝑣̅ 𝑢′ ) + 𝑢′ (𝑢′ 𝑣 ′ ) − 𝑢′ ̅̅̅̅̅̅̅
(𝑢 ′ 𝑣 ′)
𝜕𝑡 𝜕𝑥𝑗 𝜕𝑥𝑗 𝜕𝑥𝑗 𝜕𝑥𝑗

1 𝜕𝑝′ 𝜕 𝜕𝑢′
= −𝑢′ + 𝑢′ ( )
𝜌 𝜕𝑥𝑖 𝜕𝑥𝑗 𝜕𝑥𝑗

Chaque terme de cette équation va subir une modification mathématique comme suit :

𝜕𝑢 ′ 𝜕 1
Terme 1 : 𝑢′ = 𝜕𝑡 (2 𝑢′2 )
𝜕𝑡

28
 𝜕 𝜕
Terme 2 : 𝑢′ 𝜕𝑥 (𝑢̅𝑣 ′ ) = 𝑢′𝑣′ 𝜕𝑥 𝑢̅
𝑗 𝑗

𝜕 𝜕 1
Terme 3 : 𝑢′ 𝜕𝑥 (𝑣̅ 𝑢′ ) = 𝑣̅ 𝜕𝑥 (2 𝑢′2 )
𝑗 𝑗

𝜕 𝜕 1
Terme 4 : 𝑢′ 𝜕𝑥 (𝑢′ 𝑣 ′ ) = 𝜕𝑥 (𝑣′ 2 𝑢′²)
𝑗 𝑗

𝜕 𝜕
Terme 5 : 𝑢′ ̅̅̅̅̅̅̅
(𝑢 ′ 𝑣 ′ ) = 𝑢′ ̅̅̅̅̅̅̅
(𝑢 ′ 𝑣′ )
𝜕𝑥𝑗 𝜕𝑥𝑗

1 𝜕𝑝′ 1 𝜕
Terme 6 : 𝑢′ 𝜌 𝜕𝑥𝑖 = 𝜌 (𝜕𝑥 𝑢′𝑝′)
𝑖

𝜕 𝜕𝑢 ′ 𝜕 𝜕𝑢′ 𝜕𝑢′ 𝜕𝑢′

Terme 7 : 𝑢′ 𝜕𝑥𝑗 (𝜕𝑥𝑗 ) = 𝜕𝑥 (𝑢′ 𝜕𝑥 ) − 𝜕𝑥
𝑗 𝑗 𝑗 𝜕𝑥𝑗

En remplaçant dans l‘équation on trouve :

𝜕 𝜕
𝜕 1 ′2 ′ 𝑢
̅ 𝜕 1 ′2 𝜕 1 ′
( 𝑢 ) + 𝑢′ 𝑣 𝜕𝑥𝑗 + 𝑣̅ ( 𝑢 )+ (𝑣 ′2 𝑢′2 ) − 𝑢 𝜕𝑥𝑗 (̅̅̅̅̅̅
𝑢′ 𝑣 ′ )
𝜕𝑡 2 𝜕𝑥𝑗 2 𝜕𝑥𝑗

1 𝜕 𝜕 𝜕𝑢′ 𝜕𝑢′ 𝜕𝑢′

=− ( 𝑢′𝑝′) (𝑢′ )−
𝜌 𝜕𝑥𝑖 𝜕𝑥𝑗 𝜕𝑥𝑗 𝜕𝑥𝑗 𝜕𝑥𝑗

Puis on fait la moyenne de cette équation :

̅̅̅̅̅̅̅̅̅̅̅̅̅
𝜕 1 ′2 ̅̅̅̅̅̅̅̅̅̅̅
′
𝜕
𝑢
̅ ̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅
𝜕 1 ′2 ̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅
𝜕 1
( 𝑢 ) + 𝑢′ 𝑣 𝜕𝑥𝑗 + 𝑣̅ ( 𝑢 )+ (𝑣 ′ 2 𝑢′2 )
𝜕𝑡 2 𝜕𝑥𝑗 2 𝜕𝑥𝑗
̅̅̅̅̅̅̅̅̅̅̅̅̅̅
𝜕 ̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅ ̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅̅
𝜕𝑢 ′ ̅̅̅̅̅̅̅̅̅̅
′ 1 𝜕 ′ ′ 𝜕 ′ 𝜕𝑢′ 𝜕𝑢′
−𝑢 𝜕𝑥𝑗 (̅̅̅̅̅̅
′ ′
𝑢𝑣) = − ( 𝑢𝑝) + (𝑢 𝜕𝑥𝑗 ) − 
𝜌 𝜕𝑥𝑖 𝜕𝑥𝑗 𝜕𝑥𝑗 𝜕𝑥𝑗

𝜕 1 ̅2 𝜕 𝜕 1 ̅2 𝜕 1 ̅̅̅̅̅̅̅̅
( 𝑢′ ) + ̅̅̅̅̅
𝑢′𝑣′ 𝑢̅ + 𝑣̅ ( 𝑢′ ) + ( 𝑢′ 𝑢′𝑣 ′ )
𝜕𝑡 2 𝜕𝑥𝑗 𝜕𝑥𝑗 2 𝜕𝑥𝑗 2

1 𝜕 ′ ′ 𝜕 ̅̅̅̅̅̅̅̅̅
𝜕𝑢′ ̅̅̅̅̅̅̅̅̅
𝜕𝑢′ 𝜕𝑢′
=− ( ̅̅̅̅̅̅
𝑢 𝑝 )+ (𝑢′ )−
𝜌 𝜕𝑥𝑖 𝜕𝑥𝑗 𝜕𝑥𝑗 𝜕𝑥𝑗 𝜕𝑥𝑗

29
On réarrange l‘équation :

𝜕 1 ′2 𝜕 1 ′2
( 𝑢̅ ) + 𝑣̅ ( 𝑢̅ )
𝜕𝑡 2 𝜕𝑥𝑗 2

𝜕 1 ̅̅̅̅̅̅̅̅ 1 ̅̅̅̅̅ ̅̅̅̅̅̅̅̅̅

𝜕𝑢′ 𝜕 ̅̅̅̅̅̅̅̅̅
𝜕𝑢′ 𝜕𝑢′
=− ( 𝑢′𝑢′𝑣′ + 𝑣′𝑝′ − 𝑢′ ̅̅̅̅̅
) − 𝑢′𝑣′ 𝑢̅ − 
𝜕𝑥𝑗 2 𝜌 𝜕𝑥𝑗 𝜕𝑥𝑗 𝜕𝑥𝑗 𝜕𝑥𝑗

1 ̅2 = k
𝐴𝑣𝑒𝑐 : 𝑢′
2

Donc L’équation devient :

𝜕 𝜕 𝜕 1 ̅̅̅̅̅̅̅̅ 1 ̅̅̅̅̅ ̅̅̅̅̅̅̅̅̅

𝜕𝑢′ 𝜕 ̅̅̅̅̅̅̅̅̅
𝜕𝑢′ 𝜕𝑢′
𝑘 + 𝑣̅ 𝑘=− ( 𝑢′𝑢′𝑣′ + 𝑣′𝑝′ − 𝑢′ ) − ̅̅̅̅̅
𝑢′𝑣′ 𝑢̅ −
𝜕𝑡 𝜕𝑥𝑗 𝜕𝑥𝑗 2 𝜌 𝜕𝑥𝑗 𝜕𝑥𝑗 𝜕𝑥𝑗 𝜕𝑥𝑗

𝜕 𝜕 𝜕 1 ̅̅̅̅̅̅̅̅ 1 ̅̅̅̅̅ ̅̅̅̅̅̅̅̅̅

𝜕𝑢′ ̅̅̅̅̅̅̅̅̅
̅̅̅̅̅ 𝜕 𝑢̅ − 𝜕𝑢′ 𝜕𝑢′
𝑘 + 𝑣̅ 𝑘=− ( 𝑢′𝑢′𝑣′ + 𝑣′𝑝′ − 𝑢′ ) − 𝑢′𝑣′
𝜕𝑡 𝜕𝑥𝑗 𝜕𝑥𝑗 2 𝜌 𝜕𝑥𝑗 𝜕𝑥𝑗 𝜕𝑥𝑗 𝜕𝑥𝑗

Cette dernière représente l‘équation de transport de l‘énergie cinétique k

En utilisant un Benchmark on peut faire la modélisation suivante de l‘équation de k :

𝐷𝑘 𝜕𝑘 𝜕𝑘 𝜕 𝑡 𝜕𝑘
= + 𝑣̅ = [ ] + 𝑃 −
𝐷𝑡 𝜕𝑡 𝜕𝑥𝑗 𝜕𝑥𝑗 𝜎𝑘 𝜕𝑥𝑗

Avec : 𝜎𝑘 = 1

Et

𝑘 3⁄2
𝑉~𝑘 1⁄2 𝑙~ 𝜀
𝑡 ~𝑉𝑙

3.3.5 Le Modèle K-Epsilon

Le modèle k epsilon introduit aux équations RANS deux inconnus au lieu de 6

C'est deux inconnus sont :

k : l‘énergie cinétique

30
Epsilon : la dissipation de cette énergie cinétique

Ces deux inconnus vont être introduits comme conditions aux limites avec U la vitesse d'entré

Par analogie on peut trouver l‘équation de transport de epsilon comme suit :

𝜕(𝜌𝜀) 𝜕(𝜌𝜀) 𝜕 𝑡 𝜕𝜀 𝜀 𝜀²
+ 𝑣̅ = [( ) ] + 𝐶1𝜀 𝑃𝑘 − 𝐶2𝜀 𝜌
𝜕𝑡 𝜕𝑥𝑗 𝜕𝑥𝑗 𝜎𝑘 𝜕𝑥𝑗 𝑘 𝑘

Avec : 𝐶1𝜀 et 𝐶2𝜀 sont des constantes

Les k est les Epsilons vont être remplacé dans l’équation suivante :

𝑘²
𝜇𝑡 = 𝜌𝐶𝜇
𝜀

Donc 𝜇𝑡 est déterminée à partir de cette équation, 𝜇𝑡 est remplacer dans le tenseur de Reynolds
Modifié.

3.3.6 Le Modèle K-Omega

De même façon on peut trouver l‘équation de transport de Omega :

𝜕(𝜌𝜔) 𝜕(𝜌𝜔) 𝜕 𝑡 𝜕𝜔 𝛾 𝜔2
+ 𝑣̅ = [( ) ] + 𝑃𝑘 − 𝛽𝜌
𝜕𝑡 𝜕𝑥𝑗 𝜕𝑥𝑗 𝜎𝑘 𝜕𝑥𝑗 𝑡 𝑘

Avec :𝛾 , 𝜎𝑘 , 𝛽 sont des constantes

k et Omega vont être remplacer dans l’équation suivante :

𝑘
𝜇𝑡 = 𝜌
𝜔

-Donc 𝜇𝑡 est déterminée à partir de cette équation 𝜇𝑡 est remplacer dans le tenseur de Reynolds
Modifié

3.3.7 Le Modèle K-Omega SST

Pour le modèle mixte k omega SST on utilise la transition suivante

De l‘équation de transport de epsilon on remplace epsilon par :

31
 = 𝐶𝜇 𝑘𝜔

On trouve la même équation de omega mais avec un terme additionnel :

𝜕(𝜌𝜔) 𝜕(𝜌𝜔) 𝜕 𝑡 𝜕𝜔 𝛾 𝜔2 𝜌𝜎𝜔2

+ 𝑣̅ = [( ) ] + 𝑃𝑘 − 𝛽𝜌 +2 ∇k∇𝜔
𝜕𝑡 𝜕𝑥𝑗 𝜕𝑥𝑗 𝜎𝑘 𝜕𝑥𝑗 𝑡 𝑘 𝜔

Donc on peut changer entre les équations d’Epsilon et Omega en rajoutant or annulant ce terme

𝜌𝜎𝜔2
2 ∇k∇𝜔 : est le terme additionnel
𝜔

Pour cela on va introduire au terme additionnel la fonction F1 suivante :

𝜌𝜎𝜔2
2(1 − 𝐹1 ) ∇k∇𝜔
𝜔

La fonction F1 est nommée la fonction de mélange

La transition entre les deux équations ce fait de la manière suivante :

La fonction est choisie de telle façon à apparaitre près de l'obstacle qui va annuler le terme
additionnel

𝐹1 = 1 Donc

𝜌𝜎𝜔2
2(1 − 𝐹1 ) ∇k∇𝜔 = 0
𝜔

Donc le Solveur va utiliser l‘équation de transport de Omega

La fonction va disparaitre loin de l'obstacle qui va laisser le terme additionnel

𝐹1 = 0 Donc

𝜌𝜎𝜔2 𝜌𝜎𝜔2
2(1 − 𝐹1 ) ∇k∇𝜔 = 2 ∇k∇𝜔
𝜔 𝜔

Donc le Solveur va utiliser l‘équation de transport de Epsilon

32
 Fig 3.1 : La variation de la fonction de mélange par rapport a la distance

𝐹1 Prend la forme suivante :

𝑇𝑎𝑛ℎ(𝑎𝑟𝑔14 )

Fig 3.2 : La variation de la fonction de mélange pour différents modèles

33
 Fig 3.3 : la variation de l’argument par rapport à la distance

Le Code CFD va calculer d la distance à l’obstacle le plus proche avant le calcul des autres
paramètres. Le code va ensuite évaluer la valeur de l’argument

Donc F1 est calculé pour toute grille de maillage. Le Solveur va assumer quelle grille est
calculer en utilisant k-epsilon ou k-omega.

Dans le modèle k-omega SST un limiteur de viscosité est introduit.

Dans le modèle k-omega 𝜇𝑡 est calculée directement par la loi suivante :

𝑘
𝜇𝑡 = 𝜌
𝜔

Dans le cas du modèle k-omega SST 𝜇𝑡 est calculée avec cette loi :

𝑎1 𝜌𝑘
𝜇𝑡 =
max(𝑎1 𝜔, 𝑆𝐹2 )

Ou 𝐹2 est aussi une fonction de mélange

Si 𝑆𝐹2 est plus large donc la viscosité va être limitée

Le limiteur donne des résultats plus exacts et plus proches des mesures expérimentales

34
La fonction 𝐹2 est :

𝐹2 = tanh(𝑎𝑟𝑔22 )

2√𝑘 500𝜈
𝑎𝑟𝑔2 = max ( , )
𝛽 ∗ 𝑤𝑑 𝑤𝑑²

Le traitement près de l'obstacle

Ce traitement est nécessaire pour bien capter le profil de vitesse près de l'obstacle. On a besoin
donc d'un grand nombre de grille proche du notre géométrie.

Fig 3.4 : La variation de la vitesse de l’écoulement

L'approche utilisée dans ce cas est l'approche des fragments linéaires. Dans chaque grille on
pose use fonction qui décrit une petite portion du profil. Cette fonction est linéaire (elle donne
une droite). Plusieurs droites vont être montées pour construire le profil.

Une autre approche sera d'utiliser une seul grille qui va avoir une fonction non linéaire décrivant
le profil dans toute la région. Cela va rendre le calcul plus facile et aide la stabilité du schéma
numérique

35
 Fig 3.5 : comparaison entre le modèle linéaire et non linéaire

Cette fonction non-linéaire est déterminée à partir des résultats expérimentaux ou des résultats
du DNS

Fig 3.6 : La variation de la vitesse par rapport a y+

Le graph montre la vitesse tangentielle en fonction de la distance verticale à l'obstacle

Le graph est divisé de la manière suivante :

36
 Fig 3.7 : la décomposition de la fonction de vitesse

Le graph en noir est la fonction qu'on doit utiliser. Le graph bleu et vert correspond à des
fonctions empiriques qui sont mesurées dans les régions noté

Dans la région gauche ; le graph bleu est plus proche au graph noir.

Donc pour y+<5 on peut utiliser la fonction correspondante :

𝑈+ = 𝑦+

Dans la région droite ; le graph vert est plus proche au graph noir

Donc pour y+>30 on peut utiliser la fonction correspondante :

1
𝑈+ = log(𝐸𝑦 + )


Où Kappa et E sont des coefficients empiriques (0.4187 et 9.793 respectivement)

y+ et U+ sont la distance verticale à l'obstacle et la vitesse sans dimensions

𝑦𝑢𝜏
𝑦+ =

𝑈
𝑈+ =
𝑢𝜏

Avec : 𝑢𝜏 = √𝜏⁄𝜌 est une vitesse référentielle

Pour ne pas utiliser plusieurs fonctions pour décrire le développement du profil de vitesse, on
peut utiliser une autre approche

Cette approche est la fonction de Spalding :

37
 + 1 1
𝑦 + = 𝑈 + + 0.1108 [𝑒 𝑈 − 1 − 0.4𝑈 + − (𝑈 + )² − (𝑈 + )3 ]
2 6

Cette fonction est continue valable pour y+<300

Fig 3.8 : comparaison entre la fonction générée par DNS et Spalding

Les Avantages et les Inconvénients de chaque modèle

Le Modèle k-epsilon :

-Il prédit bien loin des obstacles, il est utilisé pour les écoulements turbulent seulement donc
epsilon décrit bien l’évolution de la dissipation des Eddy

-Il manque de sensibilité aux gradients de pression défavorables, à cause de la modélisation du

terme diffusif

𝜕 𝑡 𝜕𝜀
[( ) ]
𝜕𝑥𝑗 𝜎𝑘 𝜕𝑥𝑗

38
Le Modèle k-Omega :

-Prédit bien près du mur, le Omega contrôle mieux la dissipation de l’écoulement que epsilon,
donc le décollement (séparation) de la couche limite va être détecté plus précisément

-La séparation risque être excessive, car le modèle est très sensible à la valeur d’oméga.

Le Modèle k-omega SST :

-Il offre des avantages similaires à ceux du k-omega et k-epsilon

-Il donne des prédictions très précises du début et de la quantité de séparation de l'écoulement

-Il est recommandé pour les simulations de couche limite de haute précision.

-La dépendance à la distance du mur le rend moins adapté aux flux de cisaillement libres, et il
nécessite une résolution de maillage près du mur

39
 CHAPITRE 4

RESOLUTION NUMERIQUE
ET RESULTATS
4.1 Introduction

Dans ce chapitre on va présenter la discrétisation du domaine physique et la géométrie,

ainsi le maillage et sa sensibilité jusqu’à l`arrivé au maillage optimum. Nous allons utiliser ce
dernier pour simuler le comportement de la structure de la flamme turbulente de diffusion, ainsi
les différents paramètres thermodynamiques, à savoir la température, la pression, la vitesse et
les fractions massiques des réactifs et celles des espèces formées par la combustion. Ensuite on
va interprété les résultats de variation des paramètres précédents qui nous avons obtenu dans
notre calcul numérique.

4.2 PRESENTATION DE CODE DE CALCUL FLUENT

Il existe un certain nombre de codes industriels, permettant la prédiction d’écoulements

de fluides (FLUENT, CFX,…). Pour réaliser nos simulations, le choix s’est porté sur le code
de calcul Fluent.

Fluent est un logiciel de la CFD capable de modéliser les écoulements des fluides
impliquant des phénomènes physiques complexes tels que la turbulence, le transfert
thermique, les réactions chimiques, les écoulements multiphasiques pour des domaines
comportant des géométries complexes. Il est largement utilisé dans l’industrie aéronautique, il
offre une interface simple permet l'accès de l'utilisateur aux fonctions exigées pour le calcul
de la solution et l'affichage des résultats .

4.2.1 Architecture Du Logiciel

Il utilise une architecture qui lui permet de s’exécuter en tant que plusieurs processus
simultanés sur un seul ordinateur de travail ou sur des machines en parallèle, pour une
exécution plus efficace. Il est composé de trois éléments : le préprocesseur, le solveur et le
post-processeur.

 Préprocesseur : la génération du maillage est effectuée en Gambit qui est le

logiciel offert avec Fluent. Il permet de représenter la géométrie du système, de
définir les types des conditions aux limites des frontières du domaine .
 Solveur : le solveur permet de définir numériquement les conditions
opératoires (gravité, pression) dans lesquelles est effectuée la simulation, ainsi
que la spécification des conditions aux limites, la définition des propriétés du

40
 fluide et de spécifier le type du matériau. Enfin, il permet de choisir le
processus itératif pour le calcul.
 Post-processeur : le post-processeur est l’élément qui permet de visualiser la
géométrie et le maillage du domaine, mais surtout d’afficher les résultats
obtenus. Il est ainsi possible de visualiser les champs de vecteur vitesse, les
champs de pression, de turbulence ainsi que toutes les autres grandeurs
calculées sur un segment, une section du domaine ou sur tout le domaine. Il
offre aussi la possibilité de tracer des courbes et de visualiser les lignes de
courant ou la trajectoire de particules.

4.3 Les processus de la simulation

4.3.1 Création de la géométrie

La figure 1 présente le schéma qui est réalisé par Fluent du domaine de calcul montré
dans les vues latérale et arrière. Un mélange de l'hydrogène et l’azote sont utilisés comme
combustible. Le combustible entre dans le domaine avec une vitesse de 29,2 m/s. Le jet de
carburant accélère dans la conduite de carburant et s'éjecte à une vitesse proche de 34,8 m/s,
ce qui est le vitesse de sortie du brûleur mesurée dans l'expérience. La vitesse de l'air à co-
écoulement est fixée à 0,2 m/s. L'entrée d'air reçoit également une condition aux limites
d'entrée d'air (avec une vitesse d'entrée de 0,2 m/s). La pression et température sont
maintenues respectivement à 1bar et 300 K.

Fig 4.1 Schéma du domaine de calcul

4.3.2 Génération de maillage

La génération d’un maillage adapté à l’écoulement est une étape très importante dans la

41
Réalisation d’une simulation numérique. Approximativement, on constate que deux tiers du
temps est passé à la préparation de la géométrie et l’obtention d’un maillage acceptable.
Pour notre domaine de calcul le maillage est composé de 968230 cellules triangulaires de type
non structuré. Plusieurs tailles de maillage étendant entre. Le maillage ne doit être ni trop
lâche pour ne pas perdre d’informations, ni trop dense pour ne pas nécessiter trop de
ressources de calcul, il faut choisir le meilleur compromis entre la précision des résultats et le
temps de calcul.

La figure suivante présente le maillage de la chambre de combustion sous le logiciel

ANSYS

Fig 4.2 : le maillage du domaine d’étude

4.3.3 Critères de qualités de maillage

Il existe plusieurs critères pour vérifier la qualité de maillage tel que « orthogonal
quality » et « skewness », la figure ci-dessous présente ces deux derniers :

Fig 4.3 : Critères « skewness » et “orthogonal quality ».

42
 L’application de ces deux critères dans notre maillage nous donne les résultats suivants :

Fig 4.4 : vérification de Critères « orthogonal quality » dans notre maillage.

La majorité des éléments sont compris entre 0,88 et 1 d’après la figure 4.3 et 4.4 notre
maillage est de très bonne qualité.

4.3.4 Paramètres des données de calcul

la première des choses a faire c’est d’aller au menu setup, pour définir les étapes de calculs
chaque choix est justifié selon notre problème , fluent n’est plus une boite noir , mais ce sont
des étapes de calculs programmé et englobé dans le même logiciel.

43
 Fig 4.5 fenêtre fluent setup
Le menu général est destinée pour choisir le solver pressure based ou dansity based et le type
d’écoulement permanent ou transitoire, prendre en compte des effets de gravité.

 Régime permanent

la réaction de combustion étant une réaction rapide instantanée , on vas considérer un seule
instant stationnaire afin de capturer le phénomène complet.

 Solveur pressure based implcite

Pour notre domaine d’étude on doit par la suite faire choix a un modèle de combustion, on
modélise la flamme de diffusion en équilibre chimique.

Da << 1 : chimie lente le temps de combustion est inférieure aux temps caractéristiques de
turbulence.

Ka < 1 : même les plus petites structures ne peuvent pas altérer la structure de la flamme :
régime de « flamme lette ».

Fluent nous ne donne plus la main pour choisir dansity based solver donc on choisit pressure
based.

44
  Modèle de combustion ( species)


Hypothèses
Ecoulement compressible

Fluide Visqueux

Réactif

Fig (4.15) fenetre fluent hypothéses

[Link] Choix de modele du turbulence

Les écoulements turbulents sont affectés de façon significative par la présence de parois. Au
contact avec celle-ci, la vitesse du fluide est nulle. Au voisinage, la turbulence est fortement
amortie et les phénomènes dus à la viscosité moléculaire y sont prépondérants. La turbulence
augmente très rapidement lorsque l’on s’éloigne des parois. Certaines équations comme celle
de modèle K- e ne sont plus valables près des parois car l’écoulement turbulent n’est pas
pleinement établi. Une modélisation des différents phénomènes est alors nécessaire dans ce cas.
Fluent nous offre plusieurs modèles de turbulences.

Fig 4.7 fenêtre fluent modèle de

turbulence

45
 Fermeture de

L’équation de l’énergie

Fig 4.8 fenetre fluent species model.

4.3.4.2Choix de la méthode de solution

Algorithme
couplage
pression
vitesse

Fig 4.9 choix de la méthode de solution

46
4.4 résultats et commentaires

4.4.1 contour de température

Fig 4.10 contour de température

La figure 4.9 montre la propagation de température sur le plan de symétrie. A partir de la figure,
on peut voir la présence d'une région de couleur rouge à haute température qui atteint un
maximum de 2100 k autour de la zone de réaction puis une diminution progressive est observé
dans la direction d’écoulement jusqu’à l’extrémité a cause du mélange qui devient
progressivement pauvre. On voit aussi que la flamme a sa base positionnée au niveau de la lèvre
de la paroi du tuyau de carburant. Dans les expériences [10], le Laser-Doppler-Vélocimètre a
été utilisé pour obtenir les deux composantes de la vitesse simultanément. Spectroscopie
Raman-Rayleigh spontanée (SRRS/Raman) et Anti-Stokes cohérent Des techniques de
spectroscopie Raman (CARS) ont été utilisées séparément pour mesurer la température. SRRS
a été utilisé pour mesurer la fraction de mélange.

47
4.4.2 Répartition de la température le long de l'axe

Fig 4.11 répartition de la température le long de l’axe

La figure 4.11 montre le tracé de la température le long de l'axe obtenu dans les expériences et
le calcul. La sommet température dans le calcul est d'environ 16% plus élevée par rapport à
l'expérience, mais l'emplacement du pic de la température est proche de l’expérience et aussi a
l’extrémité du domaine d’étude.

48
4.4.3 contour énergie cinétique turbulente

Fig 4.12 Energie cinétique turbulente

d’après la figure 4.12 on remarque qu’au centre de la réaction l’énergie cinétique est a sa
valeur maximal et celle-ci est dû à l’augmentation de la température *valeur maximale* dans
ce point qui est d’environ 2100 k .

4.4.4 Fraction massique de O2

fig 4.13 variation de la fraction massique d’O2

49
4.4.5 Fraction massique de H2

Fig 4.14 variation de la fraction massique de H2

4.4.6 variation de fraction de mélange

Fig 4.15 variation de fraction de mélange

D’après les figures 4.13 , 4.14 et 4.15 on peut facilement remarquer que la flamme est de
diffusion, dont l’air et le carburant réactionnels introduits par des entrées séparées.

50
 On peut aussi subdiviser une flamme de combustion selon la variation de sa couleur :

 Une partie bleue allant vers le verre, comme on peut le voir dans la figure 4.13 cette
couleur est la conséquence de la concentration dominante du dioxygène.
 Une partie opaque brillante, plutôt jaune-blanc mais tirant vers le jaune-orangé à son
sommet : les produits de combustion de la zone bleue chassant l’air, la zone est donc
plus pauvre en dioxygène.
 Une zone rouge dans la figure 4.13 : cette région est également pauvre en dioxygène
par conséquent riche en dioxyde de Carbone.

Conclusion :

Ce dernier chapitre est consacré principalement à la présentation des différents résultats

obtenue à partir de logiciel ANSYS, qu’il nous a permis non seulement de mieux comprendre
la phénoménologie qui a lieu à l’intérieur de notre domaine d’étude mais aussi de vérifier les
différents paramètres telle que la distribution de la température. Ces résultats ont été comparés
avec la littérature scientifique, ou ils présentent relativement de bonne concordance et par
conséquent s’assurer de leurs validités.

51
 Liste des références
[1] Stephen R. Turns “An Introduction to Combustion Concepts and Applications” 3rd
edition Book. McGraw-Hill. New York, New York, United States. 1996
[2] Peter Atkins “The Laws of Thermodynamics_ A Very Short Introduction.” Book.
OXFORD. UK. 2007
[3] Jürg Waser “Basic Chemical Thermodynamics” Imperial College press. London,
United Kingdom. 1966.
[4] J. Wilks “The Third Law of Thermodynamics” Oxford University Press. London,
United Kingdom. 1961.
[5] “Enthalpy and internal energy liquids, solutions and vapours” Royal Society of
Chemistry The website.
[6] Achintya Mukhopadhyay "Fundamentals of Combustion Engineering" CRC Press
New York USA 2019
[7] Chung [Link] "Combustion Physics" Cambrigde University Press. United Kingdom.
2006
[8] Irvin Glassman "Combustion Fifth Edition" Elswevier. Oxford United Kingdom. 2015
[9] Maximilian Lackner "Combustion From Basics to Applications" WILEY- VCH.
Markono Print Media Pte Lt. Singapore. 2013
[10] Cheng T.-C., Fruechtel G., Neuber A., Lipp F., Hassel E. P., and Janicka J., 1995,
"Experimental data
base for numerical simulations of turbulent diffusion flames”, Forschung im Ingenieurwesen,
61 (6), pp.
165-171.