‫الجمھوریة الجزائریة الدیمقراطیة الشعبیة‬

République Algérienne Démocratique et Populaire

‫وزارة التعلیم العالي و البحث العلمي‬
Ministère de l’Enseignement Supérieur et de la Recherche Scientifique

Université Larbi Ben M’hidi Oum-El-Bouaghi

Faculté des Sciences et Sciences Appliquées
Département de Génie Mécanique

Mémoire de Fin d'Etudes en vue de l’obtention du diplôme de Master

En Génie Mécanique
Option: Turbomachines

Thème
Simulation du Comportement Transitoire
d’une Chaudière à Tubes de fumée - Démarrage à Froid

Présenté par :

Belkadi Achraf Taki Eddine

Encadreur: Dr. Ahmad Rahmani

Année universitaire: 2016 / 2017

‫بسم هللا الرحمن الرحيم‬
 Dédicaces

Je dédie ce modeste travail à

Mon père et ma mère, Mes frères et Ma sœur
A tous mes amis particulièrement : Mehdi, Yaccine,
Tarek, Zinou, Khaled, Amin, Houssem, Nedjmou
Tous mes proches, Tous mes collègues
A tous ceux qui ont contribués de près ou de loin dans ce travail
Et à tous mes enseignants et spécialement mon encadreur
Dr. Ahmad Rahmani

i
 Remerciement
Au nom de dieu, clément miséricordieux, que les grâces de dieu, sa paix, ses
salutations soit sur notre prophète MOHAMED.

Nous remercions en premier lieu, le dieu tout puissant, de nous avoir donné
courage, patience et violenté pour élaborer ce modeste travail.

La meilleure expression de remerciement et de gratitude pour l’élaboration d’un

projet à la hauteur des efforts fournis par notre encadrement –Mr. Ahmad
Rahmani

Nos promoteurs dont nous tenons à souligner leur implication et leurs conseils
pertinents qui ont éclairé un précieux apport à notre travail.

Nous remercions d’avantage le département de génie mécanique Qui nous a été

d’une aide précieuse.

ii
 ABSTRACT
In this work, we are interested to simulate the transient thermal-hydraulic behaviour of a 3-
pass fire-tube boiler during a cold start-up. A transient thermal model has been developed
based upon energy balance equations of the main components of the boiler and the relevant
correlations describing heat transfer phenomena taking place in the boiler. The variation of
the main thermal parameters describing the boiler dynamic behaviour, namely: temperatures
of combustion gases, tube walls and water as well as the heat transfer rates are also studied.
The validation of the proposed model has been performed by comparing the simulation results
against experimental data obtained during the boiler start-up. The comparative study shows a
reasonable agreement between the simulation and the experimental data of the boiler. Indeed,
the maximum error in predicting the main parameters of the boiler was found to be acceptable
and the deviations are mainly due to simplifications introduced in the model. The agreement
between the simulation and experimental data reveals the accuracy and suitability of the
proposed thermal model to estimate the dynamic behaviour of the studied boiler.

Keywords: Modeling and simulation, fire-tube boiler, transient simulation, heat transfer.

iii
 RESUME

Dans ce travail, nous sommes intéressés à la simulation du comportement thermo-hydraulique

en d'une chaudière à tubes de fumée à 3-passes lors d'un démarrage à froid. Un modèle
thermique transitoire a été développé en sur la base des équations du bilan énergétique des
principaux composants de la chaudière et des corrélations appropriées décrivant les
phénomènes de transfert de chaleur ayant lieu dans la chaudière. La variation des principaux
paramètres thermiques décrivant le comportement dynamique de la chaudière, à savoir: les
températures du gaz, de la paroi et celle de l'eau ainsi que les taux de transfert de chaleur sont
également étudiés. La validation du modèle proposé a été réalisée en comparant les résultats
de la simulation aux données expérimentales obtenues lors du démarrage de la chaudière.
L'étude comparative montre un accord raisonnable entre la simulation et les données
expérimentales de la chaudière. En effet, l'erreur maximale de prédiction des principaux
paramètres de la chaudière a été jugée acceptable et les écarts sont principalement dus aux
simplifications introduites dans le modèle. L'accord entre la simulation et les données
expérimentales révèle la précision et l'adéquation du modèle thermique proposé pour estimer
le comportement dynamique de la chaudière étudiée.

Mots-clés: Modélisation et simulation, chaudière à tubes de fumée, simulation transitoire,

transfert de chaleur.

iv
 SOMMAIRE
Introduction générale ……………………..……………………..…………….. …….. 01

Chapitre I: Introduction sur les générateurs de vapeur ……………………………. 05

I.1 Historique …………………………………………………………………....…….. 06

I.2 Classification des générateurs de vapeur ……………………………………………. 10
I.2.1 Chaudières à tubes de fumée …………………………………………………… 10
I.2.2 Chaudières à tubes d'eau ……………………………………………………….. 12
I.3 Problématiques des générateurs de vapeurs ………………………………………… 13
I.4 Facteurs conduisant à l’explosion des générateurs de vapeur ………………………. 16
I.4.1 Manque d’eau d’alimentation …………………………………………………... 16
I.4.2 Mauvaise combustion …………………………………………………………... 16
I.4.3 Fatigue et fragilisation des matériaux ………………………………………….. 16
I.4.4 Qualité d’eau …………………………………………………………………… 16

Chapitre II : Présentation de la chaudière à tube fumée …………………………… 18

II.1 Introduction ………………………………………………………………………... 19

II.2 Quelques exemple sur les chaudières à tubes de fumée ……………………………. 20
II.3 Calcul thermique dans la chaudière à tubes de fumée ……………………………... 22
II.3.1 Calcul de la combustion …………………………………………………….. 22
II.3.2 Calcul d’excès d’air …………………………………………………………. 23
II.3.3 Etude de transfert de chaleur ………………………………………………… 24
II.3.4 Echange de chaleur: Gaz → paroi interne …………………………………… 25
II.3.5 Echange de chaleur : paroi interne → paroi externe ………………………..... 25
II.3.6 Echange de chaleur : paroi externe → eau …………………………………… 26

Chapitre III : Modélisation et simulation de la chaudière à tubes de fumée ……… 27

III. Introduction ………………………………………………………………………… 28

III.2 Présentation de la chaudière à 03-passes………………………………………….. 28
III.3 Modélisation de la chaudière à tube de fumée……………………………………. 30
III.3.1 Stratégie de la modélisation ……………………………………………………. 31
III.3.2 Suppositions et simplifications …………………………………………………. 31
III.3.3 Modèle thermique ………………………………………………………………. 31
III.4 Calcul des corrélations de transfert de chaleur …………………………………… 34

v
 III.4.1 Convection côté gaz ………………………………………………………. 35
III.4.2 Convection côté eau ……………………………………………………….. 36
III.5 Propriétés thermos-physiques…………………………………………………….. 36
III.5.1 Propriétés thermos-physiques de l’eau…………………………………….. 36
III.5.2 Propriétés thermos-physiques de fumée…………………………………… 37
III.5.2 Propriétés thermo-physiques du métal ……………………………………. 37

Chapitre IV: Etude des performances thermiques de la chaudière lors du

démarrage ……………………………………………………………………………… 39

IV.1 Introduction………………………………………………………………………. 40
IV.2 Simulation du démarrage à froid………………………………………………….. 40
IV.3 Conditions initiales et aux limites……………………………………………….. 40
IV.4 Résultats et discussions……………………………………………………………. 42
IV.4.1 Validation en régime transitoire …………………………………………… 42
IV.4.2 Validation en régime stationnaire …………………………………………. 44
IV.4.3 Analyse des performances thermiques de la chaudière …………………… 44
IV.4.3.1 Température des parois des tubes …………………………………. 44
IV.4.32 Puissances thermiques échangées …………..……………………. 45
IV.5 Calcul d'erreur …………………………………………………………………….. 48
IV.6 Conclusion ………………………………………………………………………… 49

Conclusion générale …………………………………………………………………… 50

Bibliographie

vi
 LISTE DES FIGURES
Figure Titre Page
I.1 L’éolipile de Héron d'Alexandrie 1er siècle après J.C. 6
I.2 Chronologie des progrès de conception des principales machines à vapeur. 7
I.3 Historique de développement des générateurs de vapeur. 9
I.4 Présentation d’une chaudière à tubes de fumée. 11
I.5 Différents types des chaudières à tubes d’eau. 13
I.6 Explosion de la chaudière sur une locomotive ferroviaire vers 1850 14
Conséquences d’une explosion de la chaudière sur un Navire, 1905 à San
I.7 Diego California, l’explosion était due à l’éclatement d’une Soupape de 14
sûreté
I.8 Explosion de la chaudière sur une Locomotive, Texas, 1911 15
Explosion d’une chaudière de 70 bars et 375 t/h, le 19/01/2004 à 18h:40
I.9 Complexe GNL-SKIKDA-Algérie, l’accident a fait 27 morts et 74 blessés 15
avec la Destruction totale de la chaudière
II.1 Vue d’ensemble d’une chaudière à tubes de fumée 19
II.2 Coupe transversale de la chaudière à tubes de fumée à 3-passes. 20
II.3 Chaudière à tubes de fumée à simple passe. 21
II.4 Chaudière à tubes de fumée à 2-passes. 21
II.5 Chaudière à tubes de fumée à 3-passes. 22
II.6 Chaudière à tubes de fumée à 4-passes. 22
II.7 Transmission de chaleur à travers un cylindre. 24
III.1 Photographie de la chaudière à tubes de fumée 29
III.2 Coupe longitudinale de la chaudière à 3-passes 29
III.3 Modélisation de chaudière à tube de fumée à 3-passes 30
IV.1 Débit de gaz naturel (combustible) 41
IV.2 Débit de l'eau à l'entrée de la chaudière 41
IV.3 Température de l'eau à l'entrée de la chaudière 42
IV.4 Variation de la température de fumée sortie 2 ème passe lors du démarrage 43
IV.5 Variation de la température de fumée sortie 3 ème passe lors du démarrage 43
IV.6 Variation de la température de l'eau sortie chaudière lors du démarrage 43
IV.7 Variation de la température du métal des tubes lors du démarrage 45

vii
 Variation des puissances thermiques échangées vers l'eau durant le
IV.8 46
démarrage
IV.9 Variation des puissances convective et radiative au niveau du 1èr passe 47
ème
IV.10 Variation des puissances convective et radiative au niveau du 2 passe 47
IV.11 Variation des puissances convective et radiative au niveau du 3 ème passe 48
IV.12 Erreur relative entre la simulation et les données expérimentales 49

viii
 LISTE DES TABLEAUX

Tableau Titre Page

IV.1 Etude comparative en régime stationnaire 44

ix
 NOMENCLATURE

Symbole Désignation Unité

A Surface d’échange [𝑚 2 ]
Cp Chaleur massique [J/KgK]
E Excès d’air %
ei Erreur %
GAir Débit d’air [Kg/s]
GFuel Débit combustible [Kg/s]
Ggaz Débit des fumées [Kg/s]
g Facteur de transfert thermique
h Coefficient de transfert de chaleur [W/𝑚 2 K]
k Coefficient de conductivité thermique [W/𝑚 2 K]
M La masse molaire gramme
𝑚 Débit massique [Kg/S]
P Pression [bar]
PCI Pouvoir calorifique inférieur [J/Kg]
Q Puissance thermique [W]
Q Densité de flux de chaleur [W/m²]
Tad Température adiabatique des fumées [°C]
T Température [°C]
𝑇 Température moyenne [°C]
Ta Température ambiante [°C]
ΔT Différence de température [°C]
U Coefficient global de transfert de chaleur
V Coefficient de conductivité dynamique [Kg/m.s]
Re Nombre de Reynolds = G D/ μ
Pr Nombre de Prandtl = Cp μ/ k
Nu Nombre de Nusselt = h D/ k
α Coefficient de transfert de chaleur [W/m²K]
𝜀𝑔𝑖 Emissivité du gaz à la température de la paroi
𝜀𝑐𝑜2 Emissivité du CO2 à la température de la paroi
𝜀𝐻2𝑂 Emissivité du H2O à la température Tg

x
μ Viscosité dynamique [Kg / 𝑚 2 K]
ρ Masse volumique [Kg /𝑚 3 ]
σ Constante de Boltzmann = 5.67.10−8 [W/𝑚 2 𝐾 4 ]

Indices
conv Convection
ray Rayonnement
rc Rayonnement et convection
i,e Interne, Externe
entr Entrer
sort Sortie
ref Référence
g Gaz
p Paroi

xi
Introduction générale
 Introduction générale

La chaudière est l’une des composantes clés d’une installation thermique. Sa fonction
principale est de produire de l’eau chaude ou de la vapeur (à basse, moyenne ou haute
température) pour différentes raisons telles que le chauffage centralisé, la production d’eau
chaude sanitaire ou même pour la production de l’énergie électrique. Ce système est conçu
pour opérer à de longues périodes et à des variations de charge incessantes, rapides et parfois
importantes. Pendant l’opération, plusieurs problèmes peuvent être rencontrés et mener à
limiter l'opération normale de l'installation produisant des troubles sérieux dans le programme
de la production. Des accidents qui menacent l’intégrité structurale du système et provoqu ent
de sérieuses pertes économiques et parfois humaines peuvent en résulter.
Les chaudières à tubes de fumée sont exposées à un certain nombre de problèmes qui
résultent à des degrés plus ou moins élevés que pour le cas des chaudières à tubes d’eau. Les
défaillances les plus graves sont en relation avec les charges thermiques importantes
supportées par les tubes de la chaudière. En effet, le métal se trouve soumis à une variété de
contraintes qui accélèrent l’altération de cohésion des matériaux sous les effets du flux de
chaleur interne (haute température) et des phénomènes destructifs du côté externe (côté eau :
entartrage, corrosion, ébullition,…). Ces contraintes ont une influence sur le bon
fonctionnement de la chaudière, et parfois des conséquences graves, telles que des explosions.
La simulation numérique permet d’étudier les mécanismes de transfert de chaleur dans
chaque partie du système. Elle permet donc de prévoir les limites de sûreté de l’installation
pour une exploitation économique, sûre et durable. Actuellement, la simulation numérique est
devenue une voie incontournable dans l’étude de sûreté et de conception des installations
thermiques [1]. L’utilité grandissante des outils de simulation numérique repose
principalement sur le développement des méthodes numériques, l’avancement de la
programmation et l’apport de moyens informatiques puissants [2-3].
L’intérêt pratique de ce travail est de rassembler toutes les informations nécessaires sur
l’évolution des paramètres thermo-hydrauliques du système (perte de charge, coefficient de
transfert de chaleur, températures de fluide et de la paroi, pression, débit, …). Les résultats
obtenus intéressent beaucoup plus les utilités industrielles de production de la vapeur et
d’énergie électrique. En outre, cette étude permettra de mener correctement l’opération
d’expertise et de formuler des recommandations aux exploitants. Ce travail est un modèle de
calcul réaliste d’aide à la conception et à l’exploitation du système énergétique, en outre il
permettra de tester la réponse et la fiabilité des systèmes de contrôle et de protection lors des
accidents postulés.

2
 Introduction générale

Ce travail est structuré en quatre chapitres en plus d’une introduction et une conclusion
générale. Le premier chapitre présente une recherche bibliographique sur les chaudières
industrielles, leurs applications, les différents types ainsi que les facteurs conduisant à
l’explosion des chaudières à tube de fumée.
Le second chapitre présente une description détaillée d’une chaudière à tube de fumée de
03-passes de 90 kW. Cette chaudière qui est l’objet de notre étude est une chaufferie d’eau
fait partie d’une installation de chauffage centralisé. Une description succincte est présentée
pour les principales composantes de cette chaudière. Cette description concerne les différents
équipements qui le constituent ainsi que les systèmes de contrôle et de régulation associés.
Le troisième chapitre est consacré à la modélisation de la chaudière à tube de fumée de
03-passes. La modélisation de la chaudière consiste à subdiviser cette dernière en un
ensemble de volumes de gaz interconnectés vers le sens d’écoulement des fumées. À
l’intérieur des tubes, l’échange de chaleur sera considéré entre le volume gazeux et les
surfaces métalliques des tubes de la chaudière. Sur ces surfaces, on établira un bilan
énergétique qui tient compte le transfert de chaleur par rayonnement, convection et par
conduction. La chaleur produite sera transférée ensuite par convection vers l’eau qui assure le
refroidissement des surfaces externes des tubes. Le comportement thermique en régime
transitoire de la chaudière est donc défini par un système d’équations différentielles
ordinaires. La résolution des équations différentielles est effectuée par d’un code de calcul
informatique développé sous Matlab 6.1 à l'aide de la méthode Range-Kutta (ode45). Le
modèle thermique a été développé à partir des équations de conservation d’énergie pour les
principales composantes de la chaudière. Le code de calcul élaboré permet de prédire les
caractéristiques thermiques telles que les puissances thermiques échangées, les densités de
flux de chaleur et les températures de différentes composantes de la chaudière en régime
transitoire.
La simulation du comportement thermo-hydraulique de la chaudière lors d’un démarrage
à froid fera l’objet du quatrième chapitre. L’ensemble des résultats obtenus et leur discussion
sont présentés dans ce chapitre. La qualification de la modélisation est effectuée par la
confrontation des résultats de la simulation avec les relevés d’exploitation de la chaudière lors
de démarrage. L'étude comparative montre un accord raisonnable entre la simulation et les
données expérimentales de la chaudière. En effet, l'erreur maximale dans l’estimation des
principaux paramètres de la chaudière a été jugée acceptable et les écarts sont principalement
dus aux simplifications introduites dans le modèle. L'accord entre la simulation et les données
expérimentales révèle la précision et l'adéquation du modèle thermique proposé pour estimer

3
 Introduction générale

le comportement dynamique de la chaudière étudiée. En outre, les performances thermiques

de la chaudière étudiée ont été analysées de point de vue transfère thermique.

4
 Chapitre I
Introduction sur les générateurs
De vapeur
Chapitre I : Introduction sur les générateurs de vapeur

I.1 Historique

Depuis longtemps l’homme s’est intéressé aux générateurs de vapeur comme un moyen
efficace pour la production de l’énergie mécanique. L’idée d’utiliser la vapeur Comme force
motrice remonte au 1er siècle après le Christ avec l’éolipyle de Héron D’Alexandrie (Fig. I.1),
C'est une chaudière hermétique remplie en partie d’eau, placée sur le feu. Dans cet appareil
sortent deux tubes reliés à une sphère pouvant tourner autour d’un axe horizontal. De cette
sphère, deux autres tubes perpendiculaires à l’axe laissent sortir la vapeur qui, par propulsion,
fait tourner la sphère. Ce principe simple a servi pour l'ouverture et la fermeture des portes de
certains temples, l’éolipile est considéré comme un jouet. Malgré ça, cette petite Machine a
été considérée comme la première chaudière à vapeur dans l'histoire [1].

Fig. I.1: L’éolipile de Héron d'Alexandrie 1er siècle après J.C.

Depuis Héron d'Alexandrie, de nombreuses tentatives ont été réalisées, pour Mettre au point
une machine à vapeur efficace. Cependant, plusieurs machines à Vapeur ont été mises en
services depuis la formulation des lois fondamentales de la Thermodynamique. La figure I.2
décrit l’histoire d’évolution des machines à vapeur Depuis Denis Papin (1697) jusqu’à nos
jours [2].

6
Chapitre I : Introduction sur les générateurs de vapeur

Fig. I.2: Chronologie des progrès de conception des principales machines à vapeur.

Le développement de meilleures machines à vapeur est resté un objectif à atteindre Depuis

Denis Papin (1697), qui est le découvreur d’une nouvelle manière pour lever l’eau Par la force
du feu. Mais, ceci est vraiment qu’à partir de la fin du 17ème siècle Que les ingénieurs ont
développés les machines à vapeur modernes. La machine idéale De Carnot, réversible et sans
frottement, a permis de comprendre la conversion de la Chaleur en travail et de nombreux
perfectionnements l’ont amélioré tout au long du XIXème et au début du XXème siècle avec
l’introduction des machines à piston.
James Watt présente en 1768 l’usage de la vapeur à haute pression et l’action À double effet
sur les deux faces du piston [3]. En 1800, l’ingénieur américain Evans Mis au point la
première chaudière à tubes de fumée (Fire-Tube Boiler) qui servit dans les premières
locomotives.
La nécessité d’avoir des débits et pressions de vapeur importants trouve en 1867 à La mise au
point de la chaudière à tubes d’eau par les ingénieurs américains Babcok et Wilcox. Depuis,
celles-ci se sont sans continuer d’améliorer permettant d’avoir Notamment des rendements
élevés. En 1867, George Herman Babcock et Stephen Wilcox, ont développés un nouveau
type de générateurs de vapeur appelés les Chaudières à tubes d’eau (Water-Tube Boiler).

7
Chapitre I : Introduction sur les générateurs de vapeur

Nouvelle manière pour lever Machine a vapeur pour puiser l’eau des mines
L’eau par la force (thomas Savery en 1698)
De vapeur, (Denis papin en 1681).

Premier modèle de la pompe a eau Machine de thomas Newcomen améliorée

(Thomas Newcomen 1712) en 1763

8
Chapitre I : Introduction sur les générateurs de vapeur

Locomotive à vapeur, (james watt 1784)

Chauddière a tube d’eau Chauddière a tube fumée

( stirling boiler 1902)

Fig. 1.3: Historique de développement des générateurs de vapeur.

9
Chapitre I : Introduction sur les générateurs de vapeur

I.2. Classification des générateurs de vapeur

Les générateurs de la vapeur peuvent être classifiés de plusieurs façons: d'après la forme et la
disposition des tubes, nombre de réservoirs, position des réservoirs, Méthode de circulation de
l'eau et des fumées (naturelle ou forcée), service (marin ou Stationnaire), capacité (charge) et
les conditions thermiques. Généralement, les Chaudières à vapeur industrielles sont de deux
(02) catégories:
a) Chaudière à tubes de fumée, (Fire-Tube Boiler)
b) Chaudière à tubes d’eau, (Water-Tube Boiler)

I.2.1 Chaudières à tubes de fumée

Historiquement, ce sont les premières machines à vapeur utilisées dans la vieille Industrie,
comme c'est le cas pour le chemin de fer ou la navigation. L’appellation "Chaudière à tubes
de fumée" est relative aux gaz de fumée qui circulent à l’intérieur Des tubes de la chaudière
contrairement à l’autre type de chaudières à tubes d’eau. Les Gaz, résultants de la combustion
dans le foyer, circulent ensuite vers la zone de Convection qui est une série d’échangeurs de
chaleur formés par des faisceaux de tubes De faibles diamètres. Les tubes sont immergés dans
un volume d’eau stagnant avec lequel se transmet la chaleur. Donc, la température de fumées
diminue progressivement avant d’arriver à la cheminée. Les chaudières à tubes de fumées
sont habituellement limitées par une pression de 17bars et une charge de 25 tonnes/h [10].
Une pression de 10 bars est très commune pour ces unités. Leur principal usage est dans les
petites installations industrielles et pour les applications de chauffage. La figure I.4, montre
une chaudière à tubes de fumée à 3 passes.

10
Chapitre I : Introduction sur les générateurs de vapeur

Fig. I.4: Présentation d’une chaudière à tubes de fumée.

Les chaudières à tubes de fumée sont de plusieurs types :

- La chaudière Cornouaillaise (Cornish boiler)
- La chaudière Lancashire (Lancashire boiler)
- La chaudière Locomotive (Locomotive boiler)
- La chaudière à tubes de fumées verticale (Vertical Fire-Tube Boiler VFT)
- La chaudière à tubes de fumées horizontale (Horizontal Return Tubular Boiler HRT).

Les chaudières à tubes de fumées sont choisies pour les raisons suivantes:
- Perte minimum par rayonnement.
- N'exige aucune maçonnerie ou arrangement réfractaire.
- La conception empêche les fuites d'air frais et de combustibles.
- Grande capacité de production de vapeur pour un volume fini.
- N’exige pas de grandes espaces pour son installation.
- Robustesse et stabilité.
- Facilité d’exploitation.
- La teneur en eau est excellente pour des applications de chauffage.
- N'exige pas un établissement spécial.

11
Chapitre I : Introduction sur les générateurs de vapeur

I.2.2 Chaudières à tubes d’eau

Les chaudières à tubes d'eau ont été développées entre 1958 et 1978, elles sont Destinées à
travailler à des puissances et pressions trop élevées (660 MW et 220 bars). Ce type de
chaudière fournit un débit de vapeur supérieur à 20 tonnes/h, en moyenne et haute pression. Il
possède deux réservoirs appelés ballon distributeur (en partie Inférieure) et ballon collecteur
(ou encore ballon de vaporisation, en partie supérieure), reliés par un faisceau de tubes
vaporisateurs.
Dans cet ensemble, circule l’eau qui se Transforme en vapeur. Les gaz chauds produits par la
combustion sont directement en Contact avec les tubes vaporisateurs où se produit la
vaporisation. La vapeur ainsi Générée est collectée dans le ballon supérieur, l’eau
excédentaire est ramenée vers le Ballon inférieur par des tubes de chute non soumis à la
chaleur.
Il existe plusieurs types des chaudières à tubes d’eau, distingués selon leur design ou leur
utilisation. A titre d’exemple, on peut citer les chaudières (Fig. I.5): type-A, type-O, type-D et
ainsi que les chaudières nucléaires.

Chaudière type –A Chaudière type -O

12
Chapitre I : Introduction sur les générateurs de vapeur

Chaudière type –D Chaudière Nucléaire

Fig. I.5: Différents types des chaudières à tubes d’eau.

I.3 Problématique des générateurs de vapeur

Les générateurs de vapeur sont soumis à des conditions sévères de haute Température, haute
pression, environnement corrosif et un fonctionnement continu. Toutes ces conditions
menacent l’intégrité du système et peuvent dans certains cas, Provoquées l’explosion du
générateur de vapeur. Généralement, les explosions des Générateurs de vapeur sont de deux
(02) sortes :
1. Explosion des parties sous pression (côté eau).
2. Explosion de la chambre de combustion (côté feu)
Les explosions de chaudière dues à la surpression interne sont généralement Désastreuses,
souvent catastrophiques avec des dommages souvent étendus au Personnel. Voici quelques
exemples sur les conséquences d’explosion des générateurs de vapeur à travers l’histoire.

13
Chapitre I : Introduction sur les générateurs de vapeur

Fig. I.6: Explosion de la chaudière sur une locomotive ferroviaire vers 1850 [31].

Fig. I.7: Conséquences d’une explosion de la chaudière sur un Navire, 1905 à San Diego
California, l’explosion était due à l’éclatement d’une Soupape de sûreté [31].

14
Chapitre I : Introduction sur les générateurs de vapeur

Fig. I.8: Explosion de la chaudière sur une Locomotive, Texas, 1911[32].

Fig. I.9: Explosion d’une chaudière de 70 bars et 375 t/h, le 19/01/2004 à 18h:40 Complexe
GNL-SKIKDA-Algérie, l’accident a fait 27 morts et 74 blessés avec la Destruction totale de
la chaudière [33].

15
Chapitre I : Introduction sur les générateurs de vapeur

I.4 Facteurs conduisant à l’explosion des générateurs de vapeur

I.4.1 Manque d’eau d’alimentation

Le manque d’eau d’alimentation est considéré comme l’accident de référence pour les
installations thermiques équipées avec des générateurs de vapeur [4]. Il peut être provoqué par
l’arrêt des pompes alimentaires, la rupture de la tuyauterie principale d’alimentation, le
blocage de la vanne principale d’alimentation, par négligence ou encore par la défaillance de
la régulation [2]. Ces accidents menacent l'intégrité structurale du système et provoque
l’explosion des parties sous pression du générateur de vapeur [5].

I.4.2 Mauvaise combustion

Si la combustion est élevée sur la moyenne de fonctionnement, le niveau d’eau diminue

rapidement dans la chaudière et la pompe d’alimentation travaille beaucoup ce qui peut causer
sa surchauffe et par la suite son endommagement.

I.4.3 Fatigue et fragilisation du métal

Les chocs thermiques provoquent des cycles de dilatations restreintes, engendrent des
gradients de contraintes qui s’initieront généralement sur les défauts géométriques des
surfaces [6]. La fatigue à chaud est un autre facteur qui traduit une évolution en fonction du
temps et sous l’effet des contraintes de service qui ont la tendance à accélérer des phénomènes
de précipitation [7].

I.4.4 Qualité d’eau :

La qualité d’eau influe forcement sur les facteurs qui accélèrent la destruction des surfaces
internes des tubes vaporisateurs [8]. Il est obliger de traiter l’eau de la chaudière afin d’éviter
le dépôt de tartre à l’extérieur du tube du foyer et des tubes de fumée. En effet, le tartre
provoquerait un mauvais échange thermique, un temps de mis en pression-température plus
long, un risque de surchauffe au niveau du tube du foyer, une surconsommation du
combustible, une augmentation de la température des fumées à la sortie du générateur de
vapeur. L’expérience a montré qu’il est nécessaire de maintenir à une valeur convenable des
Caractéristiques chimiques de l’eau des différents circuits afin d’éviter :
a/ La corrosion : La corrosion peut définir l’ensemble des processus destructifs que Subit un
corps solide sous l’action de certains paramètres extérieurs, agents Atmosphériques ou
produits chimiques. Il existe plusieurs formes de corrosions [9] :

16
Chapitre I : Introduction sur les générateurs de vapeur

- Corrosion généralisée : cette forme de corrosion est celle que l’on rencontre le plus souvent,
l’attaque a lieu sur la totalité des surfaces métalliques sera plus au moins uniforme.
- Corrosion galvanique : le couplage de deux métaux à potentiels différents dans une solution
apporte une différence de potentiel et un passage d’électron. Le métal à faible résistance
devient donc une anode l’autre une cathode.
- Corrosion par piques (pitting): cette corrosion provoque des piques dans l’interface
métal/fluide, il existe par la combinaison entre les ions CL- et le métal lui-même. Les facteurs
influençant sur la corrosion sont :
- Le PH : plus que le PH d’eau évolue, la vitesse de la corrosion diminue et vice versa.
- Gaz carbonique : la bonne solubilité occasionne la formation d’acide carbonique
Qui s’ionise et active la corrosion correspondante.
- L’ammoniac (NH3) : L’ammoniac est principalement présenté par excès d’hydrazine ou
résultat de déchets organiques.
- L’oxygène dissous : l’oxygène est très corrosif notamment à chaud. Il favorise la réaction
chimique en fonction de la température, de l'humidité et de la pression. Il se produit une
attaque du métal et la formation d’un produit de corrosion à sa surface.
b/ L’entartrage : ou dépôt localisé, c’est la précipitation sous la forme cristalline directement
sur les parois externes des tubes de la chaudière. Le tartre joue le rôle d’un isolant thermique
au transfert de chaleur. La chaleur s’accumule dans la paroi ce qui engendre des points
chauds. Les principaux types de tartre sont :
- Carbonate de calcium(CaCO3)
- Carbonate de magnésium(MgCO3)
- Sulfate de calcium Ca (SO4)2
- Hydroxyde de magnésium Mg(OH)2
- Phosphate de calcium Ca3(PO4)2.

17
 Chapitre II
Présentation de la chaudière
À tube fumée
 Chapitre II : Présentation de la chaudière à tube fumée

II.1 Introduction

Une chaudière à tubes de fumées est constituée d'un grand réservoir d'eau traversé par des
tubes dans lesquels circulent les fumées. Le premier tube du parcours de fumées est un tube
de plus gros diamètre appelé le foyer. Ce type de construction est aujourd'hui utilisé presque
exclusivement pour les combustibles gazeux et liquides. En effet, la forme du foyer des
chaudières à tubes de fumées rend difficile l'extraction des cendres. Lorsqu'elles sont utilisées
avec des combustibles solides le foyer est placé à l'extérieur de la chaudière proprement dite.
Dans une chaudière à tubes de fumée, les gaz de combustion circulent à l'intérieur des tubes
tandis que l’eau entoure ceux-ci. L'ensemble est confiné dans une grande virole qui constitue
le corps cylindrique de la chaudière. On voit qu'un premier "tube" de très gros diamètre
constitue la chambre ou le foyer dans lequel se développe la flamme. Ensuite les fumées sont
canalisées vers des faisceaux tubulaires pour terminer leur cheminement dans un économiseur
éventuel qui assure le préchauffage des fumées par récupération sur l'enthalpie résiduelle des
fumées. L'ensemble du "tube foyer" et des faisceaux tubulaires est noyé dans le grand volume
d'eau contenu dans la virole externe qui constitue le corps de la chaudière. La figure II.1
montre une vue d’ensemble d’une chaudière à tubes de fumée de 03-passe.

Figure II.1: Vue d’ensemble d’une chaudière à tubes de fumée

La figure II.2 montre une coupe transversale de la chaudière à 3-passes. La chambre de

combustion représente le premier tube du parcours de fumées est un tube de plus gros
diamètre, qui se trouve dans le ballon même de la chaudière, sous le plan d’eau, collecte les

19
 Chapitre II : Présentation de la chaudière à tube fumée

gaz chauds en sortie de brûleur. Les gaz chauds, accumulés dans un premier caisson (1er passe
la chambre de combustion) à l’arrière de la chaudière, sont véhiculés par un groupe de tubes
immergés dans l’eau du ballon vers un second caisson (2ème passe) à l’avant de la chaudière.
Un second groupe de tubes immergés emmène les gaz vers un troisième caisson (3ème passe) à
l’arrière de la chaudière, ce troisième caisson débouche sur la cheminée pour évacuation des
fumées vers l’extérieur [1].

Figure II.2: Coupe transversale de la chaudière à tubes de fumée à 3-passes.

Ce type de chaudières est généralement réservé à des puissances inférieures à 30 MW. La

principale caractéristique des chaudières de ce type est donc de comporter un grand volume
d'eau sous pression. Donc elle peut être très destructive en cas d'explosion (un gros volume
d'eau se transforme instantanément en vapeur).

II.2 Quelques exemple sur les chaudières à tubes de fumée

Il est difficile de trouver une classification standard pour les chaudières à tubes de fumée,
celle-ci faisant intervenir trop de critères non significatifs lorsqu’ils sont pris isolément, trop
de variantes et trop de combinaisons possibles. La classification la plus courante est celle
basée sur la forme géométrique de la chaudière, précisément le nombre de passe des fumées.

20
Chapitre II : Présentation de la chaudière à tube fumée

Figure II.3 : Chaudière à tubes de fumée à simple passe.

Figure II.4: Chaudière à tubes de fumée à 2-passes.

21
 Chapitre II : Présentation de la chaudière à tube fumée

Figure II.5 : Chaudière à tubes de fumée à 3-passes.

Figure II.6 : Chaudière à tubes de fumée à 4-passes.

II.3 Calcul thermique dans la chaudière à tubes de fumée

II.3.1 Calcul de la combustion

La combustion est une réaction chimique entre un combustible et l'oxygène de l'air. Ses

principales caractéristiques sont les quantités respectives de combustible et d'air à mettre en

œuvre ainsi que la quantité de chaleur dégagée [12]. Cette réaction est accompagnée d'un

dégagement de chaleur et elle apparait sous forme de flamme.

Combustible (gasoil) + Comburant (Air) → Produits de combustion + Energie

La réaction chimique de combustion du gasoil CmH n s’écrit [25]:

22
 Chapitre II : Présentation de la chaudière à tube fumée

𝑛
𝐶𝑚𝐻𝑛 + 𝑚 + 𝑂2 → 𝑚𝐶𝑂2 + 𝑛/2𝐻2 𝑂 (𝐼𝐼.1)
4

On parlera de « combustion en excès d’air » chaque fois que l’on détectera la présence

d’oxygène dans les produits de combustion [12] . L’excès d’air peut résulter soit d’un réglage

du volet d’air au niveau du brûleur (brûleur à air soufflé), soit d’une impossibilité à régler

l’arrivée d’air (brûleur atmosphérique). Cet excès d’air est indispensable pour une

Combustion complète L’air en excès impliquera :

 une augmentation des pertes par les fumées,

 une diminution du rendement de la chaudière.

Ainsi cette augmentation entraînera un accroissement des dépenses énergétiques sur une

saison de chauffe (consommation de combustible plus importante), qui restent non

négligeables même pour de petites puissances de chaudières.

On caractérise l’air en excès par le facteur d’air Fa ou taux d’aération n en utilisant la relation

Volume d' air réellement admis au brûleur 𝑉𝑎 +𝑉𝐸𝐴

suivante : n  soit 𝑛 =
Volume d ' air théorique 𝑉𝑎

Avec : Va : Pouvoir comburivore [m3 (n) d’air / unité de combustible]

VEA : Volume d’air en excès [m3 (n) d’air en excès / unité de combustible]

𝑉𝐸𝐴 = 𝑉𝑎 𝑛 − 1

Pourcentage d’excès d’air EA% : 𝐸𝐴% = 𝑛 − 1 ∗ 100

𝑄𝑐𝑜𝑚 =𝑀𝑓𝑢𝑒𝑙 PCI (𝐼𝐼.2)

Qcom, est l’énergie dégagée par combustion et est évaluée en multipliant la masse du

combustible par le pouvoir calorifique inférieur PCI du combustible.

Selon le principe de la conservation d’énergie, la chaleur produite par combustion est égale à

la somme des énergies des réactifs (Air et fuel).

𝑀𝑔 𝐶𝑝𝑔 𝑇𝑎𝑑 − 𝑇𝑎 = 𝑀𝑓𝑢𝑒𝑙 PCI +𝑀𝑎𝑖𝑟 𝐶𝑝𝑎𝑖𝑟 𝑇𝑝𝑟 é − 𝑇𝑎 (𝐼𝐼.3)

23
 Chapitre II : Présentation de la chaudière à tube fumée

Pour une combustion complète et adiabatique, l’énergie dégagée porterait la température des

gaz à une température notée 𝑇𝑎𝑑 donnée par :

𝑀𝑓𝑢𝑒𝑙 𝑃𝐶𝐼 + 𝑀𝑎𝑖𝑟 𝐶𝑝𝑎𝑖𝑟 𝑇𝑝𝑟 é − 𝑇𝑎

𝑇𝑎𝑑 = 𝑇𝑎 + (𝐼𝐼. 4)
𝑀𝑔 𝐶𝑝 𝑔

On peut utiliser la température de la flamme comme la température adiabatique [23] :

𝑞 𝑓𝑢𝑟 𝑞 𝑓𝑢𝑟 𝑞 𝑓𝑢𝑟 𝑞 𝑓𝑢𝑟 𝑞 𝑓𝑢𝑟 2

𝑇𝑓𝑙𝑎𝑚𝑒 = 1048.15 + 0.18 + 0.048 − 20 − 0.0017 + 0.2
𝑚𝑔 𝑚𝑔 1000 𝑆1 𝑚𝑔 1000 𝑆1

(𝐼𝐼. 5)

𝑞𝑓𝑢𝑟 = 𝑄𝑐𝑜𝑚 1000

𝑆1 : La surface de 1er passe (chambre de combustion)

𝑚𝑔 : Débit massique pour le combustible

Le débit d’air nécessaire pour la combustion d’un Kilogramme du gasoil est :

𝑀𝑎𝑠𝑠𝑒 𝑑′𝑎𝑖𝑟 32 1 + 3,76 𝑀𝑎𝑖𝑟
𝐺𝑎𝑖𝑟 = 𝐺𝑓𝑢𝑒𝑙 𝐺𝑓𝑢𝑒𝑙 (𝐼𝐼. 6)
𝑀𝑎𝑠𝑠𝑒 𝑑𝑒 𝑓𝑢𝑒𝑙 𝑀𝑓𝑢𝑒𝑙
Selon la loi de conservation de la masse, le débit total des gaz de combustion est égal à La
somme des débits d’air et du combustible.
𝐺𝑔𝑎𝑧 = 𝐺𝑎𝑖𝑟 + 𝐺𝑓𝑢𝑒𝑙 (𝐼𝐼.7)
II.3.2Etude de transfert de chaleur
Durant leur passage à travers le tube cylindrique, le flux de chaleur traverse trois résistances

thermiques (Fig. II.7). La première est représentée par la couche limite des gaz de combustion

près de la paroi interne, la deuxième est la résistance de la paroi métallique, alors que la

troisième résistance est la couche limite côté eau au voisinage de la paroi externe du cylindre.

24
 Chapitre II : Présentation de la chaudière à tube fumée

Fig. II.7: Transmission de chaleur à travers un cylindre.

II.3.3 Echange de chaleur: Gaz → paroi interne :

La chaleur se transmet du gaz vers la paroi interne par convection et par rayonnement lorsque

la température du gaz est supérieure à 900°C. Ceci est valable dans la zone de radiation

(chambre de combustion) où la température des gaz de combustion est trop élevée.

Cependant, on néglige le rayonnement dans la zone de convection (2𝑒𝑚𝑒 et 3𝑒𝑚𝑒 passes).

𝑄𝑔
𝑇𝑔 -𝑇𝑝𝑖 = (𝐼𝐼.8)
𝑕 𝑖 𝐴𝑖

𝑕𝑖 est la somme des coefficients de transfert de chaleur par convection et par rayonnement

exprimé par l’équation (9), Ai la surface interne du cylindre et Tpi, la Température interne de

la paroi.

hi=hconv+hrad (𝐼𝐼.9)

25
 Chapitre II : Présentation de la chaudière à tube fumée

II.3.4 Echange de chaleur: paroi interne → paroi externe :

Le transfert de chaleur de la paroi interne vers la paroi externe s’effectue par conduction selon

la loi de Fourier [13] :

𝑄𝑔 𝑟𝑒
𝑇𝑝𝑖 -𝑇𝑝𝑒 = 𝑙𝑜𝑔 (𝐼𝐼.10)
2𝜋𝐿𝐾𝑝 𝑟𝑖

𝑇𝑝𝑖 Et 𝑇𝑝𝑒 sont respectivement les températures interne et externe de la paroi, L, est la longueur

du cylindre et kp la conductivité thermique de la paroi définie à la température moyenne de la

paroi :

𝑻𝒑𝟏+𝑻𝒑𝟐
𝑻𝒑 = (𝐼𝐼. 11)
𝟐

II.3.5 Echange de chaleur: Paroi externe → eau :

Au niveau des surfaces externes des tubes, le transfert de chaleur s’effectue parle mécanisme

d’ébullition nucléé dû à l’écart de température ∆𝑇𝑒𝑎𝑢 entre la paroi externe et la température

de l’eau . La chaleur extraite par ébullition est donnée par l’équation suivante :

𝑸𝒈
𝑻𝒑𝒆 − 𝑻𝒆𝒂𝒖 = (𝐼𝐼. 12)
𝒉𝒆 𝑨𝒆

La résolution du système d’équations 8, 10 et 12 permet d’obtenir la puissance thermique

échangée entre le volume gazeux et l’eau saturée.

𝑻𝒈− 𝑻𝒆𝒂𝒖
𝑸𝒈 = 𝒓𝒆 (𝐼𝐼. 13)
𝒍𝒐𝒈
𝟏 𝒓𝒊 𝟏
+ +𝒉
𝒉𝒊 𝑨 𝟐𝝅𝑳𝑲𝒑 𝒆 𝑨𝒆
𝒊

26
 Chapitre III
Modélisation et simulation de la
chaudière à tubes de fumée
Chapitre III : Modélisation et simulation de la chaudière à tubes de fumée

III.1 Introduction

La capacité thermique d’un générateur de vapeur est définie par ses limites mécaniques et
thermiques définissant les marges de sûreté d’exploitation. En effet, il est important de
déterminer la distribution des densités de flux de chaleur et les températures des parois des
tubes pour différents modes d’exploitation du générateur de vapeur.
Dans ce chapitre, on s’intéresse à la simulation du comportement thermique d’une chaudière à
tubes de fumée à 3 passes lors d'un démarrage à froid. La modélisation de la chaudière
consiste à subdiviser cette dernière en un ensemble de volumes gazeux interconnectés vers le
sens d’écoulement des gaz de combustion. L’échange de chaleur sera considéré entre le
volume gazeux et les surfaces métalliques des tubes de la chaudière. Sur ces surfaces, on
établira un bilan énergétique qui tient compte le transfert de chaleur par rayonnement,
convection et conduction. La chaleur produite sera transférée ensuite par convection vers
l’eau qui assure le refroidissement des surfaces externes des tubes.
L’approche mathématique est basée sur les lois de la combustion et de transfert de chaleur par
rayonnement, convection et par ébullition. Le modèle thermique transitoire a été développé a
partir des équations de conservation d’énergie pour les principales composantes de la
chaudière. Le comportement dynamique de la chaudière est donc défini par un système
d’équations différentielles ordinaires. La résolution des équations différentielles est effectuée
par un programme de calcul développé sous Matlab 6.1 à l'aide de méthode Range-Kutta. Le
code de calcul élaboré permet de prédire les caractéristiques thermiques telles que les
puissances thermiques échangées, les densités de flux de chaleur et les températures de
différentes composantes de la chaudière en régime transitoire.

III.2 Présentation de la chaudière à 03-passes

La Figure III.1 montre une photographie de la chaudière à tube de fumée objet de notre
étude. C'est une chaudière à 03-passes de 90 kW qui fait partie d’une installation de chauffage
centralisé et utilisée pour produire de l’eau chaude sanitaire. Les principaux paramètres
techniques de la chaudière sont :

- Puissance maximale : 90KW

- Combustible : Gaz Naturel
- Volume d’eau : 182 litres
- Longueur/Diamètre chambre de combustion : 836mm/334mm.
- Longueur/Diamètre 2ème passe : 652 mm/146.4 mm.
- Longueur/Diamètre 3ème passe : 826 mm/36.4 mm.
28
Chapitre III : Modélisation et simulation de la chaudière à tubes de fumée

Fig. III.1: Photographie de la chaudière à tubes de fumée [23].

La figure III.2 montre une coupe transversale de la chaudière.

Fig. III.2: Coupe longitudinale de la chaudière à 3-passes [23].

29
Chapitre III : Modélisation et simulation de la chaudière à tubes de fumée

III.3 Modélisation de la chaudière à tube de fumée

III.3.1 Stratégie de la modélisation

La modélisation d’une chaudière à tubes de fumée consiste à subdiviser le system en un

ensemble de volumes de contrôles élémentaires chacun représente une partie physique dans la
chaudière (liquide, solide ou gaz). La modélisation nécessite, d’une part une meilleure
compréhension des phénomènes physiques ayant lieu dans le système et d’autre part la
connaissance des données géométriques et physiques relatives à l’installation. Le schéma de
modélisation dépond du type de la chaudière (conception). L’approche suivie dans le cadre de
cette étude pour la modélisation de chaudière est similaire à celle proposé par Gutiérrez et al.
[26,24]. L’approche consiste à subdiviser le parcours des fumées en trois volumes de contrôle
isothermes, chacun représentant un volume gazeux entouré par une paroi métallique et le tous
immergé dans l'eau comme illustré sur la figure III.3.

Figure III.3: Modélisation de chaudière à tube de fumée à 3-passes.

Pour la chaudière à tube de fumée, les entrées sont des paramètres de conception (géométrie
de la chaudière et éléments éventuellement incorporés au système et des paramètres
opérationnels tels que le débit d’eau d'alimentation, le débit d'air et du combustible ainsi que
les températures correspondantes sans oublié la température ambiante.
Les paramètres de sorties sont les températures (eau de la chaudière, gaz de combustion,
vapeur saturée et surchauffée), les pressions (chaudière à eau chaude et chaudière à phase
vapeur), débit de vapeur, niveau d’eau, composition des gaz de combustion. En outre, il peut

30
Chapitre III : Modélisation et simulation de la chaudière à tubes de fumée

y avoir d'autres sorties en fonction des éléments complémentaires disposés comme

récupération de chaleur par clignotement de purge, de dégazeur, échangeur de chaleur,
conditionnement de la vapeur à travers une surchauffeur ou l'eau d'alimentation par un
économiseur.

III.3.2 Suppositions et simplifications

Afin de simplifier les calculs et de réduire le nombre de variables, un ensemble d'hypothèses

simplificatrices ont été considéré:
- la température des gaz de combustion est uniforme dans le volume de contrôle,
- la température de paroi métallique de tube est uniforme,
- les parois internes des tubes sont des surfaces grises,
- le transfert de chaleur au niveau des zones de renversement (boîte de transfert) est négligé,
- l'émissivité effective totale des parois est considérée 𝜀𝑝 = 0,8.

III.3.3 Modèle thermique

Les équations du bilan énergétique décrivant les différentes composantes de la chaudière sont
obtenues selon le premier principe de la thermodynamique.

𝒅𝑻
𝝆𝑽𝑪𝒑 = 𝑸𝒆𝒏𝒕𝒓. − 𝑸𝒔𝒐𝒓𝒕. (𝑰𝑰𝑰. 𝟏)
𝒅𝒕

a) Chambre de combustion

La chambre de combustion est modélisée par l’équation suivante :

𝒅𝑻𝒈𝟏
𝝆𝑽𝑪𝒑 = 𝑸𝒄𝒐𝒎 + 𝑸𝒂𝒊𝒓 − 𝑸𝒈𝟏𝒐 − 𝑸𝒓𝒂𝒅(𝒈𝟏−𝒑𝟏) − 𝑸𝒄𝒐𝒏𝒗 (𝒈𝟏−𝒑𝟏) (𝑰𝑰𝑰. 𝟐)
𝒈𝟏 𝒅𝒕

𝑄𝑐𝑜𝑚 : est l’énergie dégagée par combustion

𝑄𝑎𝑖𝑟 : est l’énergie dégagée par air

𝑄𝑟𝑎𝑑 (𝑔1−𝑝1) : Puissance échangée par rayonnement vers le métal (paroi)

𝑄𝑐𝑜𝑛𝑣 (𝑔1−𝑝1) : Puissance échangée par convection vers le matériau (paroi)

𝑸𝒄𝒐𝒎 = 𝒎𝒄𝒐𝒎 𝑷𝑪𝑰 + 𝑪𝒑,𝒄𝒐𝒎 𝑻𝒄𝒐𝒎 − 𝑻𝒓𝒆𝒇 (𝑰𝑰𝑰.3)

PCI : Pouvoir calorifique inferieur du combustible

31
Chapitre III : Modélisation et simulation de la chaudière à tubes de fumée

𝑸𝒂𝒊𝒓 = 𝒎𝒂𝒊𝒓 . 𝑪𝒑,𝒂𝒊𝒓 𝑻𝒂𝒊𝒓 − 𝑻𝒓𝒆𝒇 (𝑰𝑰𝑰.4)

𝑸𝒈𝟏𝒐 = 𝒎𝒈𝟏 . 𝑯𝒈𝟏 = 𝒎𝒈 . 𝑪𝒑,𝒈𝟏 . 𝑻𝒈𝟏 − 𝑻𝒓𝒆𝒇 (𝑰𝑰𝑰.5)

𝑄𝑔1𝑜 , L’énergie dégagée de la chambre de combustion (1 er passe) vers le 2ème passe. L’énergie
thermique échangée par rayonnement vers la paroi interne du tube de la chambre de
combustion est donnée par :
𝑸𝒓𝒂𝒅(𝒈𝟏−𝒑𝟏) = 𝒉𝒓𝒂𝒅 . 𝑨𝒊𝟏. 𝑻𝟒𝒈𝟏 − 𝑻𝟒𝒑𝟏 (𝑰𝑰𝑰.6)

hrad : Le coefficient de transfert de chaleur par rayonnement exprimé par l’équation (II.9). Le
débit massique du gaz de combustion est égale la somme des débits d’air et du combustible.

𝒎𝒈𝟏 = 𝒎𝒄𝒐𝒎 + 𝒎𝒂𝒊𝒓 (𝑰𝑰𝑰. 𝟕)

La quantité de chaleur transférée du gaz vers la paroi interne par convection est calculée par :

𝑸𝒄𝒐𝒏𝒗(𝒈𝟏−𝒑𝟏) = 𝒉𝒄𝒐𝒏𝒗. 𝑨𝒊𝟏. 𝑻𝒈𝟏 − 𝑻𝒑𝟏 (𝑰𝑰𝑰. 𝟖)

Avec, h conv est le coefficient de transfert thermique par convection exprimée par l’Eq. (II.9).

La température moyenne de la chambre de combustion est exprimée par la moyenne

algébrique entre la température à l’entrée du volume (T ad) et à la sortie (T g1o )

𝑻𝒂𝒅+𝑻𝒈𝟏𝒐
𝑻𝒈𝟏 = (𝑰𝑰𝑰. 𝟗)
𝟐

On suppose ici que la température du gaz à l’entrée de la chambre de combustion est égale à
celle de la température adiabatique Tad exprimée par l’expression (II.4).

Pour la paroi du tube de la chambre de combustion, le bilan énergétique est exprimé par
l’équation suivante :
𝒅𝑻𝒑𝟏
𝝆𝑽𝑪𝒑 = 𝑸𝒓𝒂𝒅(𝒈𝟏−𝒑𝟏) + 𝑸𝒄𝒐𝒏𝒗 (𝒈𝟏−𝒑𝟏) − 𝑸𝒄𝒐𝒏𝒗 (𝒑𝟏−𝒆𝒂𝒖) (𝑰𝑰𝑰. 𝟏𝟎)
𝒑𝟏 𝒅𝒕

Qconv(p1-eau) : La puissance thermique échangée avec l’eau en circulation au niveau de la paroi

extérieure du tube, exprimée par :

𝑸𝒄𝒐𝒏𝒗(𝒑𝟏−𝒆𝒂𝒖) = 𝒉𝒄(𝒑𝟏−𝒘) . 𝑨𝒐𝟏 . 𝑻𝒑𝟏 − 𝑻𝒘 (𝑰𝑰𝑰. 𝟏𝟏)

hC(p1-eau) : le coefficient de transfert thermique par convection coté eau 1 er passe.

32
Chapitre III : Modélisation et simulation de la chaudière à tubes de fumée

b) Zone de convection 2 (2ème passe)

Pour le deuxième passe l’équation qui traduit l’évolution de la température du gaz (Tg2) lors
du démarrage est exprimé par :

𝒅𝑻𝒈𝟐
𝝆𝑽𝑪𝒑 = 𝑸𝒈𝟏𝒐 − 𝑸𝒈𝟐𝒐 − 𝑸𝒓𝒂𝒅(𝒈𝟐−𝒑𝟐) − 𝑸𝒄𝒐𝒏𝒗 (𝒈𝟐−𝒑𝟐) (𝑰𝑰𝑰. 𝟏𝟐)
𝒈𝟐 𝒅𝒕

Les puissances échangées ont la même formulation sauf que l’indice correspond au 2 ème
passe.

𝑸𝒈𝟐𝒐 = 𝒎𝒈 . 𝒉𝒈𝟐 = 𝒎𝒈 . 𝑪𝒑,𝒈𝟐 . 𝑻𝒈𝟐 − 𝑻𝒓𝒆𝒇 (𝑰𝑰𝑰. 𝟏𝟑)

𝑸𝒓𝒂𝒅(𝒈𝟐−𝒑𝟐) = 𝒉𝒓𝒂𝒅 . 𝑨𝒊𝟐. 𝑻𝟒𝒈𝟐 − 𝑻𝟒𝒑𝟐 (𝑰𝑰𝑰. 𝟏𝟒)

𝑸𝒄𝒐𝒏𝒗(𝒈𝟐−𝒑𝟐) = 𝒉𝒄𝒐𝒏𝒗. 𝑨𝒊𝟐. 𝑻𝒈𝟐 − 𝑻𝒑𝟐 (𝑰𝑰𝑰. 𝟏𝟓)

Avec, h conv est le coefficient de transfert thermique par convection exprimée par l’Eq.(II.9).

La température moyenne du volume de gaz dans le 2 ème passe 𝑇𝑔2 est exprimée par la
moyenne algébrique entre la température à l’entrée du volume (Tg1o) et à la sortie (T g2o)

𝑻𝒈𝟏𝒐+𝑻𝒈𝟐𝒐
𝑻𝒈𝟐 = (𝑰𝑰𝑰. 𝟏𝟔)
𝟐

Pour la paroi du tube du 2 ème passe, le bilan énergétique est exprimé par:

𝒅𝑻𝒑𝟐
𝝆𝑽𝑪𝒑 = 𝑸𝒓𝒂𝒅(𝒈𝟐−𝒑𝟐) + 𝑸𝒄𝒐𝒏𝒗 (𝒈𝟐−𝒑𝟐) − 𝑸𝒄𝒐𝒏𝒗 (𝒑𝟐−𝒆𝒂𝒖) (𝑰𝑰𝑰. 𝟏𝟕)
𝒑𝟐 𝒅𝒕

𝑸𝒄𝒐𝒏𝒗(𝒑𝟐−𝒆𝒂𝒖) = 𝒉𝒄(𝒑𝟐−𝒆𝒂𝒖) . 𝑨𝒐𝟐. 𝑻𝒑𝟐 − 𝑻𝒆𝒂𝒖 (𝑰𝑰𝑰. 𝟏𝟖)

hC(p2-eau) : le coefficient de transfert thermique par convection coté eau 2 ème passe.

c) Zone de convection 3 (3ème passe)

Pour le 3ème passe l’équation qui traduit l’évolution de la température du gaz (Tg 3) lors du
démarrage est exprimé par :

𝒅𝑻𝒈𝟑
𝝆𝑽𝑪𝒑 = 𝑸𝒈𝟐𝒐 − 𝑸𝒈𝟑𝒐 − 𝑸𝒓𝒂𝒅(𝒈𝟑−𝒑𝟑) − 𝑸𝒄𝒐𝒏𝒗 (𝒈𝟑−𝒑𝟑) (𝑰𝑰𝑰. 𝟏𝟗)
𝒈𝟑 𝒅𝒕

Les puissances échangées ont la même formulation sauf que l’indice correspond au 3 ème
passe.

33
Chapitre III : Modélisation et simulation de la chaudière à tubes de fumée

𝑸𝒈𝟑𝒐 = 𝒎𝒈 . 𝒉𝒈𝟑 = 𝒎𝒈 . 𝑪𝒑,𝒈𝟑 . 𝑻𝒈𝟑 − 𝑻𝒓𝒆𝒇 (𝐈𝐈𝐈. 𝟐𝟎)

𝑸𝒓𝒂𝒅(𝒈𝟑−𝒑𝟑) = 𝒉𝒓𝒂𝒅 . 𝑨𝒊𝟑. 𝑻𝟒𝒈𝟑 − 𝑻𝟒𝒑𝟑 (𝐈𝐈𝐈. 𝟐𝟏)

𝑸𝒄𝒐𝒏𝒗(𝒈𝟑−𝒑𝟑) = 𝒉𝒄𝒐𝒏𝒗. 𝑨𝒊𝟑. 𝑻𝒈𝟐 − 𝑻𝒑𝟑 (𝐈𝐈𝐈. 𝟐𝟐)

Avec, h conv est le coefficient de transfert thermique par convection exprimée par l’Eq.(II.9).

La température moyenne du volume de gaz dans le 2 ème passe 𝑇𝑔3 est exprimée par la
moyenne algébrique entre la température à l’entrée du volume (T g2o) et à la sortie (T g3o).

𝑻𝒈𝟐𝒐+ 𝑻𝒈𝟑𝒐
𝑻𝒈𝟑 = (𝑰𝑰𝑰. 𝟐𝟑)
𝟐

Pour la paroi du tube du 3 ème passe, le bilan énergétique est exprimé par:

𝒅𝑻𝒑𝟑
𝝆𝑽𝑪𝒑 = 𝑸𝒓𝒂𝒅(𝒈𝟑−𝒑𝟑) + 𝑸𝒄𝒐𝒏𝒗 (𝒈𝟑−𝒑𝟑) − 𝑸𝒄𝒐𝒏𝒗 (𝒑𝟑−𝒆𝒂𝒖) (𝑰𝑰𝑰. 𝟐𝟒)
𝒑𝟑 𝒅𝒕

𝑸𝒄𝒐𝒏𝒗(𝒑𝟑−𝒆𝒂𝒖) = 𝒉𝒄(𝒑𝟑−𝒆𝒂𝒖) . 𝑨𝒐𝟑. 𝑻𝒑𝟑 − 𝑻𝒆𝒂𝒖 (𝑰𝑰𝑰. 𝟐𝟓)

hC(p3-eau) : le coefficient de transfert thermique par convection coté eau 3 er passe.

d) Volume d’eau

La variation de la température de la masse d’eau contenue dans la cuve de la chaudière est

exprimée par l’équation suivante. Le volume d’eau reçoit la chaleur par convection par les
tubes des trois passes et échange cette énergie avec l’eau qui traverse la chaudière.

𝒅𝑻𝒆𝒂𝒖
𝝆𝑽𝑪𝒑 = 𝑸𝒄𝒐𝒏𝒗 (𝒑𝟏−𝒆𝒂𝒖) + 𝑸𝒄𝒐𝒏𝒗 (𝒑𝟐−𝒆𝒂𝒖) + 𝑸𝒄𝒐𝒏𝒗 (𝒑𝟑−𝒆𝒂𝒖) − 𝑸é𝒄𝒉𝒂𝒏𝒈é𝒆 (𝑰𝑰𝑰. 𝟐𝟔)
𝒆𝒂𝒖 𝒅𝒕

La puissance thermique échangée vers l’eau est exprimée en fonction des températures
d’entré et de sortie d’eau.

𝑸é𝒄𝒉𝒂𝒏𝒈é𝒆 = 𝒎𝑪𝒑 𝑻𝒆𝒂𝒖_𝒔𝒐𝒓𝒕−𝑻𝒆𝒂𝒖_𝒆𝒏𝒕𝒓 (𝑰𝑰𝑰. 𝟐𝟕)

𝒆𝒂𝒖

III.4 Calcul des corrélations de transfert de chaleur

Le coefficient d’échange thermique (hg) dépend de la géométrie de la paroi, les propriétés

thermodynamiques du fluide (Viscosité, densité, conductivité,…), ainsi que le régime
d’écoulement laminaire ou turbulent.

34
Chapitre III : Modélisation et simulation de la chaudière à tubes de fumée

𝑻𝟒𝒈 −𝑻𝟒𝒆𝒂𝒖
𝒉𝒓𝒂𝒚 = 𝝈𝒈𝒓𝒂𝒚 (𝑰𝑰𝑰. 𝟐𝟖)
𝑻𝒈 −𝑻𝒆𝒂𝒖

𝑨𝑻
𝒈𝒓𝒂𝒚 = 𝟏 𝟏 (𝑰𝑰𝑰. 𝟐𝟗)
+𝑪 −𝟏
𝜺𝒑 𝜺𝒈

𝝈 : La constante de Boltzmann

𝒈𝒓𝒂𝒚 : gray est le facteur de transfert thermique total par rayonnement du gaz vers la paroi

L’émissivité totale de la paroi, 𝜀𝑝 =0.85, et l’émissivité totale du gaz, 𝜀𝑔 , C le rapport des

surfaces.

III.4.1 Convection côté gaz

Le coefficient de transfert de chaleur par convection en régime d'écoulement laminaire
(Re ≤ 2300) est déterminé par l'équation de B. Spang [27] :

𝑫 𝟏.𝟑𝟑
𝟎.𝟎𝟖𝟔 𝑹𝒆𝑷𝒓
𝑳
𝑵𝒖𝒈 = 𝟒. 𝟑𝟔𝟒 𝑫 𝟎.𝟖𝟑
(𝑰𝑰𝑰. 𝟑𝟎)
𝟏+𝟎.𝟏𝑷𝒓 𝑹𝒆
𝑳

𝒉𝒈 𝑫
Avec 𝑵𝒖𝒈 =
𝛌𝒈

Pour les régimes d'écoulement de transition et turbulent (Re > 2300), on s’est basé sur la
Corrélation de Petukhov [27] où toutes les propriétés du fluide sont évaluées à la température
Moyenne (Tg+Tp1)/2.

𝒇𝒈
𝑹𝒆𝑷𝒓
𝟖
𝑵𝒖𝒈 = 𝒇𝒈
(𝑰𝑰𝑰. 𝟑𝟏)
𝟐/𝟑
𝟏+𝟏𝟐.𝟕 𝟖
𝑷𝒓 −𝟏

Pour pouvoir appliquer la corrélation de Petukhov dans la région d’entrée et la zone de

transition, Gnielinski (1976) a modifié l’équation (III.31) en remplaçant le nombre de
Reynolds par (Re - 1000) :
𝒇𝒈
𝟖
(𝑹𝒆−𝟏𝟎𝟎𝟎)𝑷𝒓
𝑵𝒖𝒈 = 𝒇𝒈
(𝑰𝑰𝑰. 𝟑𝟐)
𝟐/𝟑
𝟏+𝟏𝟐.𝟕 𝑷𝒓 −𝟏
𝟖

fg, c’est le facteur de frottement.

Cette corrélation est utilisée pour l’estimation du coefficient de transfert de chaleur dans la

35
Chapitre III : Modélisation et simulation de la chaudière à tubes de fumée

Chambre de combustion, elle est multipliée par un coefficient de correction qui tient en
compte le rapport D/L.
𝒇𝒈
𝟖
(𝑹𝒆−𝟏𝟎𝟎𝟎)𝑷𝒓 𝑫 𝟐/𝟑
𝑵𝒖𝒈 = 𝟏+ (𝑰𝑰𝑰. 𝟑𝟑)
𝒇𝒈 𝟐/𝟑 𝑳
𝟏+𝟏𝟐.𝟕 𝟖
𝑷𝒓 −𝟏

III.4.2 Convection côté eau

La relation qui présente le transfert de chaleur par convection côté eau est donnée par :

𝟐
𝟏𝟔/𝟗 𝟏/𝟔
𝑲𝒇 𝑹𝒂
𝑵𝒖𝒆𝒂𝒖 = 𝟎. 𝟔 + 𝟎. 𝟑𝟖𝟕 𝟎.𝟓𝟓𝟗 𝟗/𝟏𝟔
(𝑰𝑰𝑰. 𝟑𝟒)
𝑫𝒐
𝟏+
𝑷𝒓

𝐾𝑓 : Coefficient de transfert de chaleur par conduction

𝐷𝑜 : Le diamètre extérieur

𝑅𝑎 , 𝑃𝑟 Sont les nombres de Rayleigh et Prandtl pour l’eau

III.5 Propriétés thermos-physiques

Les propriétés thermo-physiques (des produits de combustion, parois métalliques et l’eau),

tels que : émissivité, capacité calorifique, conductivité thermique, sont formulées sous forme
des corrélations polynomiales en fonction de la température.

III.5.1 Propriétés thermos-physiques de l’eau

Physical properties of sub-cooled water versus temperature are given according to [28].
𝝆𝒘 = 𝟗𝟗𝟗. 𝟗 + 𝟐. 𝟎𝟑𝟒𝐱𝟏𝟎−𝟐𝑻 − 𝟔. 𝟏𝟔𝟐𝐱𝟏𝟎−𝟑 𝑻𝟐 + 𝟐. 𝟐𝟔𝟏𝐱𝟏𝟎−𝟓𝑻𝟑 − 𝟒. 𝟔𝟓𝟕𝐱𝟏𝟎−𝟖𝑻𝟒
𝑪𝒑𝒘 = 𝟏𝟏𝟎𝟕𝟏. 𝟎𝟔𝟒𝟖𝟓 − 𝟔𝟎. 𝟏𝟗𝟖𝟐𝟔𝑻 + 𝟎. 𝟏𝟕𝟑𝟐𝟕𝑻𝟐 − 𝟏. 𝟔𝟒𝐱𝟏𝟎−𝟒𝑻𝟑
𝒌𝒘 = −𝟎. 𝟓𝟑𝟏𝟐 + 𝟎. 𝟎𝟎𝟓𝟒𝟑𝑻 − 𝟑. 𝟓𝟓𝟓𝟕𝟒𝐱𝟏𝟎−𝟔𝑻𝟐 − 𝟔. 𝟏𝟔𝟓𝟏𝐱𝟏𝟎−𝟗𝑻 𝟑
𝟏
𝝁𝒘 = 𝟒. 𝟐𝟖𝟒𝟒𝐱𝟏𝟎−𝟓 + 𝟐
𝟎. 𝟏𝟓𝟕 𝑻 + 𝟔𝟒. 𝟗𝟗𝟑 − 𝟗𝟏. 𝟐𝟗𝟔

36
Chapitre III : Modélisation et simulation de la chaudière à tubes de fumée

III.5.2 Propriétés thermos-physiques de fumée [25]

𝟏𝟏.𝟒𝟖𝟒𝒏+𝟎.𝟐𝟐𝟕𝟐𝒏𝑯
𝝆𝟎 = 𝟏. 𝟐𝟗𝟑
𝟏𝟏.𝟒𝟖𝟒𝒏+ 𝟎.𝟐𝟐𝟕𝟐𝒏+𝟎.𝟕𝟏𝟖𝟔𝟓 𝑯

𝟐
𝒕 𝒕
𝑪𝒑𝒈 = 𝟗𝟕𝟏. 𝟕 + 𝟏𝟎. 𝟒𝟗𝒎 + 𝟑𝟐𝟓. 𝟑𝟓 − 𝟒. 𝟗𝟕𝒎 − 𝟕𝟔. 𝟓𝟗 − 𝟔. 𝟎𝟕𝒎
𝟏𝟎𝟎𝟎 𝟏𝟎𝟎𝟎
𝟖.𝟗𝟑𝟔𝑯
𝒎=
𝟏𝟏.𝟒𝟖𝟒𝒏+𝟏+𝟎.𝟐𝟐𝟕𝟐𝒏𝑯

La chaleur massique totale 𝐶𝑝𝑔 des fumées est évaluée entre𝑇𝑎 et 𝑇𝑎𝑑 , elle est donnée en

fonction des chaleurs massiques des constituants du mélange gazeux.

La chaleur massique de chaque élément gazeux en fonction de la température des gaz𝑇𝑔 .

H est le pourcentage hydrogène dans le combustible utilisé

m: pourcentage d'humidité en masse des gaz de combustion

La conductivité thermique pour le gaz est :

𝑲𝒈 =
𝒕 𝒕 𝟐
𝟐𝟏. 𝟗𝟐𝟒 − 𝟎. 𝟎𝟑𝟑𝟕𝒎 + 𝟖𝟔. 𝟒𝟔𝟕 + 𝟎. 𝟎𝟗𝟔𝟔𝒎 − (𝟏𝟐. 𝟗𝟗𝟏𝟎. 𝟔𝟐𝟐𝟗𝒎) 𝟏𝟎−𝟑
𝟏𝟎𝟎𝟎 𝟏𝟎𝟎𝟎

La viscosité dynamique du gaz et présenté par :

𝝁𝒈 =
𝒕 𝒕 𝟐
𝟏𝟔. 𝟖𝟔 − 𝟎. 𝟏𝟏𝟎𝟔𝒎 + 𝟒𝟑. 𝟒𝟒𝟗 − 𝟎. 𝟏𝟏𝟏𝒎 − (𝟏𝟏. 𝟏𝟗 + 𝟎. 𝟎𝟗𝟖𝟓𝒎) 𝟏𝟎−𝟔
𝟏𝟎𝟎𝟎 𝟏𝟎𝟎𝟎

L’émissivité totale du gaz est donnée en fonction de la température par la formule :

𝜺𝒈 = 𝟎. 𝟑𝟎𝟓𝟒𝟏 − 𝟏. 𝟓𝟔𝟗𝟖. 𝟏𝟎−𝟒 𝑻𝒈 + 𝟔. 𝟓𝟕𝟏𝟖. 𝟏𝟎−𝟖 𝑻𝒈𝟐 − 𝟑. 𝟗𝟏𝟎𝟕. 𝟏𝟎−𝟏𝟏 𝑻𝒈𝟑

+ (𝟖. 𝟒𝟐𝟐𝟔. 𝟏𝟎−𝟏𝟓)𝑻𝒈𝟒

III.5.3 Propriétés thermos-physiques du métal

La capacité de chaleur et la conductivité thermique du matériau du tube de chaudière en
fonction de la température sont données par: comme:

37
Chapitre III : Modélisation et simulation de la chaudière à tubes de fumée

𝐶𝑝𝑝 = 322.4 + 0.4155 𝑇𝑝 + 273.15

𝑘𝑝 = 770.89 + 0.44525 𝑇𝑠 + 273.15

38
 Chapitre IV:
Etude des performances thermiques de la
chaudière lors du démarrage
Chapitre IV : Etude des performances thermiques de la chaudière lors démarrage

IV.1 Introduction

Ce chapitre concerne l’étude et la simulation du comportement dynamique de la chaudière à

tube de fumée lors d’un démarrage à froid. La qualification de la modélisation est effectuée
par la confrontation des résultats de la simulation avec les relevés d’exploitation de la
chaudière lors de démarrage. Les résultats du programme sont confrontés avec les relevées
d’exploitation de la chaudière à disponible dans le référence [23].

Une importance particulière est portée à l’étude des performances thermiques de la chaudière.
L’étude montre que la conception et la disposition des tubes dans la chaudière joue un rôle
déterminant sur le comportement thermique de la chaudière.

En plus, une étude qualitative est faite pour examiner l’influence des paramètres opératoires
(température ambiante, débit de carburant, l’air ambiant et l'excès d'air sur les performances
thermiques de la chaudière.

IV.2 Simulation du démarrage à froid

Selon les hypothèses données précédemment (chapitre III), le comportement thermo-

hydraulique de la chaudière est décrit à travers sept (07) équations différentielles
correspondant aux principales composantes du système, à savoir: volume de gaz, métal des
tubes des trois passes ainsi que du volume d'eau. Pour cette raison, un programme de calcul
informatique a été élaboré sous le logiciel Matlab pour résoudre les équations différentielles et
prédire le comportement thermique transitoire de la chaudière lors du démarrage.

IV.3 Conditions initiales et aux limites

Les conditions initiales pour le programme de calcul représentent l’état initial de la chaudière.
On a supposé que les principales températures sont égales à celles de la température ambiante
juste à l’instant de démarrage. La variation temporelle des paramètres d'entrée de la chaudière
(débit d'eau et débit d'air et température de l'eau d'entrée) est introduite dans le modèle comme
conditions aux limites.
Les paramètres thermo-hydrauliques (débits, températures) imposés lors du démarrage sont
donnés dans les figures 1 à 3 [23]. Ces paramètres sont imposés au modèle comme conditions
aux limites. Vue l'allure complexe des courbes, ils ont introduite au programme par la
méthode d'interpolation linaire offert par le logiciel Matlab. La méthode d'interpolation
linéaire est utilisée en raison de sa simplicité dans la construction de nouveaux points de

40
Chapitre IV : Etude des performances thermiques de la chaudière lors démarrage

données dans la gamme d'un ensemble discret de points de données expérimentaux. Autres
paramètres nécessaires à l'exécution du programme de calcul sont:

 Température ambiante : 23°C

 Température du gaz : 19.4 °C
 Composition chimique du gaz naturel : 95% CH 4 , 2.5% C2H6 et gaz inertes

Pouvoir calorifique inférieur du combustible : 39.7 MJ/Nm3
 Masse volumique du combustible : 0.81 kg/Nm3
 Excès d'air : 21.6 %

Figure IV.1: Débit de gaz naturel (combustible)

Figure IV.2: Débit de l'eau à l'entrée de la chaudière.

41
Chapitre IV : Etude des performances thermiques de la chaudière lors démarrage

Figure IV.3: Température de l'eau à l'entrée de la chaudière.

La méthode du 4ème/5ème ordre Runge-Kutta a été utilisée pour résoudre numériquement les
équations de régulation pour un temps de calcul total de 1300 s. La résolution est faite à l'aide
de la fonction Ode45.

IV.4 Résultats et discussions

Les résultats de simulation sont donnés sous forme des courbes décrivant le comportement
dynamique de la chaudière lors du démarrage. Dans notre cas, on a considéré, pour notre
comparaison, la température des fumées à la sortie de 2ème passe, la température des fumées à
la sortie de 3ème passe et la température de l’eau lors de démarrage. La confrontation des
résultats de simulation avec les relevés d’exploitation du générateur de vapeur sont données
dans les figures IV.4-6 pour le régime transitoire et dans le tableau IV.1 pour le régime
stationnaire.

IV.4.1 Validation en régime transitoire

La variation de la température de fumée sortie 2 ème et 3ème passe est illustrée sur les figures
IV.4 et IV.5. La figure IV.6 montre la variation temporelle de la température de l'eau sortie
chaudière. La comparaison montre qu’il y a une concordance acceptable entre l’expérience et
la simulation. L'erreur relative maximale dans la prédiction des températures (fumées et de
l’eau) varie entre 1.8% et 7.45%. Les résultats obtenus comportent quelques déviations qui
sont due principalement aux simplifications et suppositions introduites lors de la modélisation
du générateur de vapeur. Précisément la négligence du transfert de chaleur au niveau des
zones de transition (boîte de transition) entre les passages de la chaudière.

42
Chapitre IV : Etude des performances thermiques de la chaudière lors démarrage

Figure IV.4: Variation de la température de fumée sortie 2ème passe lors du démarrage.

Figure IV.5: Variation de la température de fumée sortie 3ème passe lors du démarrage.

Figure IV.6: Variation de la température de l'eau sortie chaudière lors du démarrage.

43
Chapitre IV : Etude des performances thermiques de la chaudière lors démarrage

IV.4.2 Validation en régime stationnaire

A la fin du transitoire (au-delà de1300 s), la chaudière atteigne un régime quasi-stationnaire

vue le maintiens des paramètres d'entrés de la chaudière à savoir: débit de gaz, débit d'eau et
la température de l'eau à l'entrée de la chaudière (Fig. 1-3).

Tableau IV.1: Etude comparative en régime stationnaire

Paramètre Mesure Simulation Erreur

Température d'eau sortie chaudière, (°C) 77.7 76.3 1.8 %

Température de fumée sortie 2 ème passe, (°C) 712 739 3.8 %

Température de fumée sortie 3 ème passe, (°C) 535 575 7.47 %

IV.4.3 Analyse des performances thermiques de la chaudière

IV.4.3.1 Température des parois des tubes

La prédiction des températures des parois et précisément celle de la chambre de combustion

est d’une importance majeur dans l’analyse de sûreté de chaudière. Sachant que le métal se
trouve soumis à des conditions de part et d’autre. Au côté gaz, la température de la flamme est
de l’ordre de 1500-1800°C et cela juste à quelque seconde du démarrage. Par contre, au côté
eau, le fluide est la température ambiante.

La figure IV.7 illustre l’évolution des températures des parois des tubes de chaque passe dans
la chaudière. Il est clair que la température de la paroi dépond de celle des fumées en
circulation. La température de tube de la chambre de combustion est supérieure à celle du 2ème
passe et celle de ce dernier est plus élevée que la température des tubes du 3 ème passe. L’écart
maximum de température entre la paroi de tube de la chambre de combustion et l’eau Tp1-T w
est de 20.6°C. Dans le 2ème passe, l’écart de température Tp2-T w est 15.5°C. Alors que pour le
3ème passe, la différence de température Tp3 -Tw est 5.47°C.

44
Chapitre IV : Etude des performances thermiques de la chaudière lors démarrage

Fig. IV.7: Variation de la température du métal des tubes lors du démarrage.

L’étude montre que la différence maximale de température entre la paroi de la chambre de

combustion et l’eau est en dessous la limite admissible et que la chaudière fonctionne sous des
conditions normale.

IV.4.3.2 Puissances thermiques échangées

La figure IV.8, montre la variation des puissances thermiques échangées vers l’eau durant le
démarrage de la chaudière. Les calculs montrent que l’échange de chaleur est plus important
au 1èr passe (chambre de combustion) avec un taux de 51.74%, suivie par le 3 ème passe par
34.68%. Alors que, le 2 ème passe fournie 12.58% de la puissance totale échangée.

Il est clair que l’échange thermique dans la chambre de combustion (0.9 m²) est le plus
important dans la chaudière. Ceci est dû à la température trop élevés des gaz de combustion à
cet endroit est qui atteint 1030°C.

On constate que le 3ème passe est plus efficace dans l’échange thermique par rapport au 2ème
passe. Ce qui n’est pas le cas dans d’autres types de chaudières à 3-passe où la puissance
échangée décroit tout au long du parcours des fumées [29 et 30]. Ceci est grâce à la surface
d’échange du 3ème passe (A3o=2.3m²) qui représente la surface d’un faisceau de 22 tube est qui
est plus importante que celle du 2 ème passe (A2o 0.32m²) formée par un seul tube. Donc, cette
conception (01 tube) ne favorise pas le transfert thermique à cet endroit, ce qui traduit par une
température relativement élevée, de 712°C, à la sortie du 2ème passe (Fig. IV.5).

45
Chapitre IV : Etude des performances thermiques de la chaudière lors démarrage

Ceci indique que le rayonnement thermique est présent dans le 2 ème passe et aussi dans le 3 ème
passe (Fig. IV.10 et IV.11).

Figure IV.8: Variation des puissances thermiques échangées vers l'eau durant le démarrage.

L’analyse du transfert de chaleur pour chaque passe montre que le transfert de chaleur entre
les fumées et les parois des tubes tout au long de la chaudière s’effectue mutuellement par
rayonnement et par convection. Dans la chambre de combustion (Fig. IV.9), le rayonnement
thermique est le mode de transfert de chaleur le plus dominant avec un taux de 90.6%. Alors
que la convection forcée gaz-paroi représente seulement 9.4%. Ceci est expliqué par l’effet de
la température élevé de la flamme. Sachant que le rayonnement thermique est proportionnel à
la 4ème puissance de la température (Q=f(T 4)).

Dans le 2ème passe (Fig. IV.10) et suite au déficit du transfert thermique, la température des
fumées reste élevée et le rayonnement reste toujours dominant avec un taux d’échange de
60.45% par rapport à 39.5% à la convection.

La figure IV.11 montre la variation des puissances échangée par rayonnement et par
convection au niveau du 3 ème passe. Le rayonnement existe toujours en raison de la
température élevée avec laquelle les fumées entrent le 3 ème passe. Le rayonnement thermique
dépasse légèrement la convection par 53% contre 47%.

46
Chapitre IV : Etude des performances thermiques de la chaudière lors démarrage

D’habitude, dans les chaudières à 3-passes avec faisceau tubulaire dans le 2ème et le 3ème
passe, le rayonnement thermique est négligeable dans le 3 ème passe.

Figure IV.9: Variation des puissances convective et radiative au niveau du 1èr passe.

Figure IV.10: Variation des puissances convective et radiative au niveau du 2ème passe.

47
Chapitre IV : Etude des performances thermiques de la chaudière lors démarrage

Figure IV.11: Variation des puissances convective et radiative au niveau du 3ème passe.

IV.5 Calcul d'erreur

La qualification du modèle mathématique est faite par l’estimation de l'écart entre les données
théoriques de la simulation et les mesures expérimentales de la chaudière lors du démarrage.
Ce pourcentage d'erreur peut être exprimé comme suit:

𝑿 𝒊 − 𝒀𝒊
𝒆𝒊 = . 𝟏𝟎𝟎 (𝑰𝑽. 𝟏)
𝑿𝒊
Avec, Xi est la valeur théorique obtenue par la simulation et Yi est la valeur expérimentale de
paramètres comparés.

La Figue IV.12 montre une comparaison qualitative entre les données théoriques et
expérimentales pour les principales températures décrivant le comportement thermique de la
chaudière, à savoir la température de l'eau sortant de la chaudière, la température du gaz à la
sortie de la 2ème passe et la température du gaz à la sortie de la 3 ème passe.

Généralement, un bon accord a été trouvé entre le modèle thermique proposé et les données
expérimentales. L'erreur maximale dans l'estimation des paramètres de la chaudière a été
jugée acceptable et les écarts sont principalement dus aux simplifications introduites au
modèle.

48
Chapitre IV : Etude des performances thermiques de la chaudière lors démarrage

Figure IV.12: Erreur relative entre la simulation et les données expérimentales.

IV.6 Conclusion
Le comportement dynamique de la chaudière à tube de fumée à 3-passe a été numériquement
étudié dans ce chapitre. La validation du modèle proposé a été faite en comparant les résultats
théoriques du code avec les mesures expérimentales de la chaudière lors d’un démarrage.
L'étude comparative montre un accord raisonnable entre la simulation et les données
expérimentales de la chaudière. En effet, l'erreur maximale dans l’estimation des principaux
paramètres de la chaudière a été jugée acceptable et les écarts sont principalement dus aux
simplifications introduites dans le modèle. L'accord entre la simulation et les données
expérimentales révèle la précision et l'adéquation du modèle thermique proposé pour prévoir
le comportement dynamique de la chaudière étudiée.
Les performances thermiques de la chaudière étudiée ont été analysées. L’étude montre que la
conception et la disposition des tubes dans la chaudière joue un rôle déterminant sur le
comportement thermique de la chaudière.

49
Conclusion générale
 Conclusion générale

Dans le cadre de cette étude on s’est intéressé à la modélisation et la simulation du

comportement transitoire d’une chaudière à tube de fumée de 03-passes lors d’un démarrage à
froid. L’objectif principal est l’étude du comportement thermo-hydraulique de la chaudière à
travers l’évolution de ses performances thermiques. Un code de calcul est élaboré pour la
résolution des équations de conservations différentielles régissant l’installation et l’estimation
des caractéristiques thermiques de la chaudière telle que les températures (eau, parois, fumée),
les puissances thermiques échangées. L’estimation des taux de transfert de chaleur par
rayonnement, convection et conduction dans chaque partie du système a été faite par une
combinaison des corrélations empiriques et semi-empiriques traitants Ces corrélations ont été
choisies en raison de leurs étendues d'utilisation dans une large gamme d’application. La
préparation des données d’entrée du problème nécessite des efforts considérables vue la
quantité importante d’informations requises pour la simulation. Un effort supplémentaire a été
déployé pour remédier au manque de certaines données souvent inévitables pour ce type de
travail. Cette étape, nécessite une connaissance approfondie du fonctionnement de
l’installation ainsi que l’ensemble des phénomènes physiques ayant lieu dans la chaudière.
Une validation de la modélisation est faite par la confrontation des résultats de simulation
obtenus par le code de calcul élaboré avec les données d’exploitation de la chaudière. L'étude
comparative montre que les résultats de la simulation concordent avec les relevés
expérimentaux.
En outre, une étude qualitative est faite pour analyser l’effet des paramètres Opératoires; débit
du combustible, excès d’air, surface d’échange sur les performances thermiques du générateur
de vapeur. L’étude a montré que le comportement de ces paramètres est influencé par les
conditions d’exploitation.
Ce travail nous a permis d’élaborer un modèle de base pouvant simuler non seulement le
fonctionnement de la chaudière mais également certains types d’accidents qui pourraient avoir
lieu lors du fonctionnement. La modélisation adoptée dans le cadre de cette étude est relative
à une chaudière à 3-passes et applicable pour d’autres types similaires. Notre modèle de calcul
reste toujours sujet à des améliorations futures.
L’intérêt pratique de ce travail est de rassembler toutes les informations nécessaires sur
l’évolution des paramètres thermo-hydrauliques du système. Les résultats obtenus intéressent
beaucoup plus les utilités industrielles de production de la vapeur et d’énergie électrique. En
outre, cette étude permettra de mener correctement l’opération d’expertise et de formuler des
recommandations aux exploitants. Ce travail est un modèle de calcul réaliste d’aide à la

51
 Conclusion générale

conception et à l’exploitation du système énergétique, en outre il permettra de tester la

réponse et la fiabilité des systèmes de contrôle et de protection lors des accidents postulés.

52
Bibliographie

[1] Christophe Alleau, "Chaudières à vapeur à combustible", fichier PDF.

[2] A. Rahmani, "Modélisation et simulation thermohydraulique d'un générateur de vapeur
industriel en régime transitoire par le code RELAP5/Mod3.2", Thèse de Doctorat, Décembre
2007.
[3] Encyclopédie scientifique de Yahoo, (www.yahoo.fr.10-05-2001)
[4] J. P. Ylén, M. Paljakka, T. Karhela, J. Savolainen & K. Juslin, "Experiences on Utilising
Plant Scale Dynamic Simulation in Process Industry", Proceedings 19 th European Conference
on Modelling and Simulation, 2005.
[5] Boiler Design, "Maintenance and Safety Study Guide, Boiler Design and Construction",
American Trainco Inc., SGBO001 306, 2006.
[6] Michel F. Asahby, "Matériaux-Microstructure et mise en oeuvre". Vol. 02, Edition Dunod,
1991.
[7] "Choix des Alliages Métalliques Réfractaires pour Four de Traitement Thermique", centre
des industries mécaniques, (CETIM-1997).
[8] A. Rahmani, "Etude du phénomène de la crise d’ébullition dans un générateur de vapeur" ,
Thèse de magistère, 2002.
[9] F.Armanet et J.H. Davidson, " Resistance A La Corrosion A Haute Température des
Aciers Inoxydables Et Alliages Réfractaires", Edition de paris.
[10] Philip J. Potter, "Power Plant Theory and design", 2nd Edition, 1959, pp. 208- 212.
[11] Trabelsie Soumia et Telailia Rafika, PFE "modélisation et simulation d’une
Chambre de combustion de la chaudière FCB par un code de calcul", 2005.
[12] Diagramme d’Ostwald du gaz naturel (Combustion) Fichier Doc
[13] Saccadura J.F 1982. "Initiation au transfert thermique", édition paris, pp87-181
[14] Truelove, J.S., "Thermal and Hydraulic Design of Heat Exchangers, Furnaces and
Combustion Chambers, Heat Exchanger Design Handbook", Hemispherical
Publishing Corporation. 1983
[15] D.C Groeneveled, S.C.Chengm L.K.Leung and C.N.Guyen, "Compulation of
single and two phase heat transfer rates suitable for water-cooled tubes and sub
channels", Novembre 1988, pp61-75.
[16] K. J. Astrôm, R.D. Bell, "Drum-boiler dynamics Automatica", Vol. 36, 2000, pp.
363-378.
[17] G.F. Hewitt, G.L. Shires & T.R. Bott, "Process Heat Transfer", CRC-Press, 1994,
pp. 777-798.
[18] Hans Dieter Baehr & Karl Stephan, "Heat and Mass Transfer", pp 370-371.
[19] R. Siegel & J. R. Howell, "Thermal Hadiation Heat Transfer", 2nd Edition, 1980
[20] Bernhard, Spang, "Correlations for convective heat transfer", presented at the
Chemical Engineers’Resource page, www.cheresources.com (2004).
[21] John H. Lienhard IV & John H. Lienhard V, "A Heat Transfer Textbook", 3rd
Edition, pp. 341-362.
[22]Pierre Wuithier, "Le pétrole, Rafinage et Génie Chimique", Tome II, 1972.2nd
revised edition, Mars 2006
[23] "Fire tube heat generators: Experimental analysis and modeling" Alberto Bisetto a,
Davide Del Col a, *, Marco Schievano b .
[24] K. Sørensen, C.M.S. Karstensen, T. Condra, N. Houbak, Modelling and
simulating fire tube boiler performance, in: Proceedings from SIMS 2003e44th
Conference on Simulation and Modeling, Session 2b, Lecture 7, 2003 in Västerås,
Sweden (ISBN 91-631-4716-5).
[25] D. Annaratone, Steam Generators, Chapter 7: Combustion, Springer-Verlag
Berlin Heidelberg 2008. DOI: 10.1007/978-3-540-77715-1 7, 129
[26] F.J Gutiérrez Ortiz, Modeling of fire-tube boilers, Applied Thermal Engineering
31 (2011) 3463-3478.
[27] Bernhard, Spang, "Correlations for convective heat transfer", presented at the
Chemical Engineers’Resource page, www.cheresources.com (2004).
[28] Lienhard-IV, J.H., Lienhard-V, J.H., 2003. Heat transfer textbook. 3rd Edition,
Phlogiston Press, Cambridge, USA.

[29] A. Rahmani, A. Dahia, "Thermal-hydraulic Modeling of the Steady-state

Operating Conditions of a Fire-tube Boiler", International Journal of Nuclear
Technology & Radiation Protection, pp. 29-37,Vol. 24, n. 1, 2009.

[30] A. Rahmani, S. Trabelsi, "Numerical investigation of heat transfer in 4-pass fire-

tube boiler", American Journal of Chemical Engineering, 2014; 2(5): 65-70.
Sites d'internet
[31] www.steamlocomotive.com
[32] http://en.wikipedia.org/wiki/Boiler_explosion
[33] www.BIBCochran.com