Université LARBI BEN M’HIDI (ULBM)

Faculté des Sciences et Sciences Appliquées (FSSA)
Département de Génie Mécanique (DGM)
Option : Master Energétique/Installations énergétique et
turbomachines
Année/Semestre : Première année/S2

Cours de Combustion
Partie 1 : Thermochimie (Max 6 semaines)
Dr. MAMERI Abdelbaki

1 
 
  
 

Table des matières/ Partie 1 (thermochimie)

Chapitre 1 (1.5 semaines)
Rappels de la thermodynamique
1.1 Propriétés extensives et intensives

1.2 Equation d’état des gaz parfaits

1.3 Equations calorifiques d’état

1.4 Mélange des gaz parfaits

1.5 Chaleurs latentes de vaporisation

1.6 Première loi de la thermodynamique pour un système à masse constante
1.7 Première loi pour un volume de contrôle constant
1.8 Chimie des combustibles

Chapitre 2 (1.5 semaines)

Mélanges de réactifs et produits de combustion
2.1 Stœchiométrie
2.1 Rapport stœchiométrique, richesse et excédant d’air
2.2 Enthalpie et enthalpie de formation absolue (standard)
2.3 Enthalpie de combustion (chaleur de réaction)
2.4 Pouvoirs calorifiques

Chapitre 3 (1.5 semaines)

Température adiabatique de la flamme
3.1 Température adiabatique de la flamme à pression constante
3.2 Température adiabatique de la flamme à volume constant

Chapitre 4 (1.5 semaines)

Equilibre chimique
4.1 Deuxième loi de la thermodynamique
4.2 Fonction de Gibbs ou énergie libre

2 
 
  
 

Chapitre 1

Combustion et thermochimie

I Rappels de la thermodynamique

1.1 Propriétés extensives et intensives

La valeur numérique d’une propriété extensive dépend de la quantité, masse où nombre

de moles, d’une substance considérée. Les propriétés extensives sont notées par des majuscules ;

par exemple, V [m3] pour le volume, U [J] pour l’énergie interne, H[J](=U+PV) pour l’enthalpie,

etc.

D’autre part une propriété intensive (spécifique) est exprimée par unité de masse où de

moles, sa valeur numérique est indépendante da la quantité ou de la substance. Les propriétés

intensives sont notées par des minuscules ; par exemple, u[J/kg] pour l’énergie spécifique

interne, h[J/kg](=u+Pv) pour l’enthalpie spécifique, etc.

Il faut noter que la température et la pression sont des variables intensives, elles sont notées par

T et P. Les variables intensives pondérées par des moles sont notées par des barres, par exemple

et (J/mole).

Les propriétés extensives sont obtenues par multiplication de celles intensives

correspondantes, par unité de masse ou de moles, par la quantité de la masse ou de moles ;

c.à.d. :

(1.1)

, etc

3 
 
  
 

1.2 Equation d’état des gaz parfaits

Une équation d’état donne une relation entre la pression P, la température T et le volume

V (ou volume spécifique v) pour une substance. Pour un gaz parfait, dans lequel les forces

intermoléculaires et le volume des molécules sont négligeable, l’équation d’état s’écrit :

(1.2)

La constante spécifique du gaz R est reliée à la constante universelle des gaz Ru (=8315

[J/(kmolK]) et à la masse moléculaire du gaz M par :

(1.3)

La masse et la masse volumique, densité, 1/ .

1.3 Equations calorifiques d’état

Les expressions qui relient l’énergie interne (ou enthalpie) à la pression et la température

sont dites équations calorifiques d’état, c.à.d. :

, (1.4)

h , (1.5)

Pour une variation différentielle dans u ou h, on écrit :

(1.6)

(1.7)

4 
 
  
 

On introduit les capacités calorifiques à volume et pression constantes :

(1.8)

(1.9)

Pour un gaz parfait, l’énergie interne et l’enthalpie ne dépendent que de la température

cela se traduit par :

(1.10)

(1.11)

Si en plus la température est constante T = cste alors dT=0 ce qui donne du =0 et dh=0

c.a.d u et h sont constantes.

On peut calculer l’enthalpie et l’énergie interne pour un gaz parfait si on définit un état de

référence ‘ref’, on intègre les relations (1.10) et (1.11), on obtient :

→ (1.12)

→ (1.13)

Pour les gaz parfaits et réels, les capacités calorifiques, chaleurs spécifiques, Cp et Cv

sont fonction de la température. C’est une conséquence de l’énergie interne d’un gaz qui est

générée par le mouvement de translation, de rotation et vibration de ses molécules. Selon la

théorie quantique, le stockage en vibration et translation devient actif lorsque la température

augmente.

5 
 
  
 

La figure 1 montre schématiquement ces modes de stockage d’énergie pour une espèce

monoatomique où l’énergie interne est constituée par l’énergie cinétique de translation. Par

contre, pour une molécule diatomique on a les trois modes, translation, vibration et rotation.

Ces modèles très simples nous permettent de conclure que les chaleurs spécifiques des

molécules triatomiques sont plus importantes par rapport à celles des molécules diatomiques et

qui sont à leurs tours plus grandes que celles des molécules monoatomiques.

Translation 

Rotation  Vibration 

Translation 

Figure 1 Mouvements des molécules Mono et Di atomiques Figure 2 Chaleurs spécifiques des espèces

6 
 
  
 

1.4 Mélange des gaz parfaits

Soit un mélange composé de n gaz (espèces), on utilise deux concepts pour caractériser sa

composition: fractions molaires et fractions massiques. Si on a N1 moles de l’espèce « 1 », N2

moles de l’espèce « 2 », …….etc. La fraction molaire de l’espèce « i » dans le mélange est

définie par :

(1.14)
⋯

Similairement, la fraction massique de l’espèce « i » est définie par :

(1.15)
⋯

On peut vérifier facilement que :

∑ ⋯ 1
⋯ ⋯ ⋯

∑ ⋯ 1
⋯ ⋯ ⋯

On peut calculer la masse molaire du mélange gazeux en fonction des fractions molaires par :

é ∑ (1.16)

Cela donne une relation entre les fractions molaires et massiques :

∑ ∑
(1.17)
∑ é

é
Ou bien : (1.18)

Calculons la masse molaire du mélange en fonction des fractions massiques, on a :

é ∑ (1.19)

7 
 
  
 

Les fractions molaires des espèces sont aussi utilisées pour calculer les pressions

partielles des espèces correspondantes. La pression partielle, de l’espèce « i », Pi est la pression

qui règne si l’espèce « i » est isolée à la même température dans le même volume du mélange.

Pour les gaz parfaits, la pression du mélange est la somme des pressions partielles.

∑ (1.20)

La pression partielle est reliée à celle du mélange par :

(1.21)

Pour les mélanges des gaz parfaits, les propriétés spécifiques du mélange sont calculées

par la pondération en fractions massiques où molaires des propriétés spécifiques individuelles.

Par exemple, les enthalpies du mélange sont calculées par :

é ∑ (1.22)

é ∑ (1.23)

Il faut noter que l’enthalpie et l’énergie interne pour un gaz parfait ne dépendent pas de la

pression. L’entropie du mélange est aussi calculée par la même façon,

é , ∑ , (1.24)

̅ é , ∑ ̅ , (1.25)

Dans ce cas les entropies des espèces ( et ̅ ) dépendent des pressions partielles des
espèces. Les entropies peuvent être évaluées pour l’état standard (Pref=P0=1atm), puisque le
mélange des gaz a une seule température T, toutes les espèces ont cette température. L’entropie
s’écrit donc :

Δ ln ⇒ , , (1.26)

Δ ̅ ̅ ln ⇒ ̅ , ̅ , (1.27)

8 
 
  
 

Les entropies standard des différentes espèces sont données sous forme de tableaux

ou de polynômes.

1.5 Chaleurs latentes de vaporisation

Dans plusieurs procédés de combustions, le changement de phase liquide-vapeur est

important. On définit la chaleur latente de vaporisation hfg, par la chaleur nécessaire pour

vaporiser complètement une unité de masse de liquide à pression constante et à une température

donnée.

, , , (1.28)

Où T et P sont la température et la pression correspondante de saturation. La chaleur latente de

vaporisation est aussi connue par enthalpie de vaporisation. Les chaleurs latentes de vaporisation

pour différents combustibles à leurs points d’ébullition, à une pression de 1 atm, sont données

sous formes de tableaux.

Exemple

1. Soit un mélange gazeux dans une enceinte où la pression est 3 atm, la température 1500K
et la composition est donnée par la table.

N° Espèce Nombre de moles

1 0.095
2 6
3 7
4 34
5 0.005

1. Calculer la masse molaire du mélange,

2. Calculer les fractions molaires des espèces,
3. Calculer les fractions massiques des espèces,
4. Calculer la pression partielle de chaque espèce,
5. Calculer l’enthalpie du mélange.

9 
 
  
 

Solution
1. La masse molaire du mélange est donnée par :

é ∑ ∑ ∑

0.095 ∗ 28 6 ∗ 44 7 ∗ 18 34 ∗ 28 0.005 ∗ 30 28.55 Kg/Kmole

.

2. Les fractions molaires des espèces sont données par :

.
→ 2.017 10 , 1.274 10
⋯ . .

.
1.486 10 , 7.219 10 , 1.062 10
. . .

3. Les fractions massiques des espèces sont données par :

→ 2.017 10 9.974 10 ,
é .

1.274 10 1.963 10 , 1.486 10 9.369 10

. .

7.219 10 7.080 10 , 1.062 10 1.116 10

. .

4. Pour les pressions partielles peuvent être calculées par :

→ 2.017 10 ∗3 6.051 10 , 3.822 10 ,

4.458 10 , 2
2.166 , 3.186 10 .

L’enthalpie du mélange est :

é ∑ ∑ ∆ avec l’enthalpie de formation de l’espèce i et

∆ l’enthalpie sensible de la même espèce. Ces quantités sont données par les tables .
N° Espèce Nombre de moles (kj/kmol) ∆ 1500 (kj/kmol)
1 0.095 -110 541 38 847
2 6 -393 546 61 681
3 7 -241 845 48 181
4 34 0 38 404
5 0.005 90 297 39 732
é 2.017 10 110 541 38 847 1.274 10 393 546 61 681
1.486 10 241 845 48 181 7.219 10 0 38 404
1.062 10 90 297 39 732 28411.31 /

10 
 
  
 

1.6. Première loi de la thermodynamique pour un système à masse constante (système

fermé) : La première loi traduit la conservation d’énergie pour un système à masse constante
entre deux états « 1 » et « 2 » par :

∆ (1.29)
Chaleur ajoutée au  Travail exercé par le  Variation de l’énergie 
système pour aller de  système sur les limites en  totale en allant de l’état      
l’état «1» à «2».  allant de l’état « 1 » à « 2 ». «1» à «2». 

et sont des fonctions qui dépendent du chemin suivi et sont ajoutées aux

limites ; ∆ est la variation de l’énergie totale du système, qui est la somme de

l’énergie interne, cinétique et potentielle, c.à.d. :

(1.30)

L’équation (1.30) peut être écrite pour des quantités intensives ou spécifiques, cela donne :

∆ (1.31)
m=cte

Q
1
Cylindre/Piston 
2
W

Si l’évolution du système est instationnaire, l’équation de la conservation d’énergie s’écrit :

(1.32)
Taux instantané de chaleur  Taux instantané de travail réalisé  Taux de variation instantanée 
transférée au système   par le système où puissance  de l’énergie totale  du système 

Reprenons l’équation (1.31), si on suppose qu’il n’y a pas de travail externe et

d’échange de chaleur en phase de combustion (par exemple le cas des MCI où la combustion se

produit dans un laps de temps très court par rapport au mouvement du piston et au temps

caractéristique de l’échange de chaleur, aussi le volume de la chambre peut être considéré

constant durant cette phase) , dans ce cas l’équation se simplifie à :

11 
 
  
 

∆ 0 (1.33)

Ou bien :

0 (1.34)

Généralement les vitesses à l’entrée et à la sortie de la chambre de combustion sont

comparables ou même négligeable. Aussi les dimensions de la chambre ne sont pas importante

ce qui permet d’éliminer l’effet de la pesanteur. On écrit donc :

0 (1.35)

Sachant que (avec v est le volume massique ), cette équation s’écrit :

0 (1.36)

En réarrangeant les termes, on écrit :

(1.37)

En termes de grandeurs extensives, cette équation s’écrit (multipliant par la masse m):

Notons que le volume de la chambre de combustion est constant au moins le long de la phase de

combustion , cela donne :

(1.38)

1.7 Première loi pour un volume de contrôle fixe :

Prenons le cas d’un volume de contrôle où le fluide peut traverser les limites de l’entrée

« e » vers la sortie « s », le travail des forces de pression (travail d’aspiration et d’éjection) sur

les frontières est :

12 
 
  
 

La première loi s’écrit donc :

(1.39)
Taux d’énergie transférée au  Taux de travail réalisé par le  Taux net d’énergie (me)et de travail 
volume de contrôle à travers la  système en excluant le travail du  associé à la pression du fluide (mPv) à 
surface des limites   fluide et incluant celui des aubes  travers le volume de contrôle. 

Surface de contrôle

Ecoulement    0  

Qcv Wcv

Dans l’écriture de l’équation (1.39) on a supposé que :

1- Le volume de contrôle est fixe par rapport au système de coordonnées. Cela élimine toute

interaction du travail associé aux déplacements des limites et aussi la nécessité de

considérer les variations dans l’énergie cinétique et potentielle du volume de contrôle lui-

même.

2- Les propriétés du fluide sont partout indépendantes du temps.

3- Les propriétés du fluide sont uniformes à l’entrée et à la sortie du volume de contrôle,

cela permet de prendre des valeurs constantes au lieu de faire des intégrales de surface.

4- On suppose qu’il y’a une seule entrée et une seule sortie pour simplifier la forme finale

du résultat tout en gardant la possibilité de faire une extrapolation vers plusieurs entrées

et sorties.

L’énergie spécifique totale du fluide « e » est composée de l’énergie interne, cinétique et

potentielle c.a.d :

(1.40)

Le terme de pression dans l’équation (1.39) peut être combiné avec l’énergie spécifique

interne de (1.40) pour donner l’enthalpie :

(1.41)

13 
 
  
 

La combinaison de (1.39) et (1.41) et après réarrangement, on aura :

(1.42)

Pour un kg de fluide on aura :

(1.43)

Dans le cas des chambres de combustion des turbomachines, on néglige le travail externe
et les pertes de chaleur, aussi les vitesses à l’entrée et la sortie sont comparables et la chambre est
horizontale. L’équation (1.43) se simplifie à :

0 (1.44)

Ou bien en termes de quantités extensives :

0 (1.45)

1.8. Définitions et propriétés physico-chimiques des combustibles et carburants.

Un combustible est quelque chose qui peut brûler. Un carburant est un combustible qui alimente
un moteur pour produire de l’énergie mécanique.

1.8.1. Définitions des combustibles

Dans le langage courant, le terme « combustible » est souvent réservé aux produits utilisés pour le
chauffage (bois, charbon, produits pétroliers…). Dans les faits, ce qualificatif s’applique à tout composé
susceptible de s’unir à un oxydant (presque toujours l’oxygène de l’air) et capable de se consumer.

Dans la première acception, « les combustibles sont des produits dont la combustion en présence d'air
dans les brûleurs, foyers, fours ou chaudières fournit de l'énergie thermique. Celle-ci est utilisée dans le
secteur domestique et commercial (chauffage, cuisine, parfois éclairage), dans l'industrie (apport de
chaleur dans les réactions chimiques endothermiques, production de vapeur et d'électricité), ou encore
dans l'agriculture (séchage des récoltes).

Les combustibles se distinguent des carburants qui sont destinés à la production d'énergie mécanique dans
les moteurs. On notera que certains produits comme les gaz de pétrole liquéfiés (GPL), le fioul
domestique (FOD) et le fioul lourd peuvent être utilisés à la fois comme combustibles et carburants. »

Les combustibles se répartissent en trois grandes catégories selon leur état physique habituel. On
distingue ainsi:

 les combustibles solides (charbon, bois, paille…) ;

14 
 
  
 

 les combustibles liquides (GPL, FOD, fioul lourd et produits assimilés), tous issus du pétrole brut
et parfois, pour une très faible part, de la biomasse (dérivés d'huiles végétales) ;
 les gaz combustibles, le gaz naturel constituant, à lui seul, une source d'énergie primaire
importante.
La seconde acception plus large englobe évidemment un plus grand nombre de substances, dont la
vocation dans l'industrie, n'est pas forcément d'être brûlées pour produire de l'énergie. Dans ce cas
l’appellation combustible renvoie aux propriétés chimiques, aux précautions de manipulation et de
stockage.
Certains produits combustibles ont la propriété de s’enflammer vivement et de brûler avec production de
flammes ; ils sont qualifiés de « produit inflammable ». Pour caractériser l’inflammabilité des liquides, on
utilise la notion de point d’éclair. Le « point d'éclair » est la température minimale à laquelle, dans des
conditions d’essais spécifiés, un produit émet suffisamment de gaz inflammables capables de
s’enflammer momentanément en présence d’une source d’inflammation.
Par ailleurs, les gaz, vapeurs, brouillards de produits combustibles, mélangés à l’air, sont explosifs dans le
domaine de concentration compris entre la limite inférieure d'explosivité (LIE) et la limite supérieure
d’explosivité (LSE). En dessous de la LIE, le mélange est trop pauvre en combustible pour donner lieu au
phénomène explosif. Au-dessus de la LSE, le mélange est trop riche en combustible et ne contient pas
suffisamment d’oxygène.
Le triangle du feu montrant les trois éléments indispensables pour une combustion chimique

La combustion est une réaction chimique d’oxydation d’un combustible par un comburant, qui nécessite
une source d’énergie. Elle met en jeu trois éléments nécessaires :
 le combustible, la matière susceptible de brûler (solide : charbon… ; liquide : essence, solvants…
; gazeuse : propane, butane…), etc.
 le comburant, en se combinant avec un autre corps, permet la combustion (oxygène, peroxydes,
chlorates…)
 l’énergie d’activation, l’énergie minimum nécessaire au démarrage de la réaction chimique de
combustion ; elle est apportée par une source de chaleur, une étincelle…
qui constituent le triangle du feu.

L’absence d’un des trois éléments empêche le déclenchement de la combustion.

La suppression d’un des trois éléments arrête le processus.

Le feu s’éteint de lui-même, s’il n’y a pas assez de comburant, si le combustible manque ou si le foyer est
refroidi.

1.8.2. Types de carburants

Comme son nom l'indique, un carburant contient du carbone. Dans le cas de l'hydrogène par exemple on
parlera de combustible. Le carburant est souvent un liquide et parfois un gaz.

15 
 
  
 
Les carburants peuvent être classés de différentes manières. On peut par exemple distinguer les carburants
fossiles, des biocarburants, des ergols.
 Les carburants fossiles sont ceux qui proviennent de la transformation des matières organiques
mortes mélangées à divers minéraux à de grandes profondeurs. Cette transformation nécessite
plus d'un million d'années et se déroule à des températures et pressions très élevées. Les
carburants fossiles sont en quantité limitée sur terre. La diminution future de la production de
carburants est traitée dans l'article sur le pic de production. Les carburants issus du pétrole sont
aussi appelés hydrocarbures.
 Les biocarburants eux, proviennent de plantes ou animaux (non fossilisés). Parmi les produits
développés commercialement, citons le bioéthanol et le biodiesel.
 Les ergols liquides voir solides sont utilisés dans les moteurs de fusées anaérobie, rare type de
moteur capable de fonctionner en dehors de l'atmosphère.
1.8.3. Le carburant, un élément stratégique

Les carburants, encore le plus souvent dérivés du pétrole, deviennent stratégiques dans l'économie de
nombreux secteurs industriels, agricoles, forestiers, halieutiques, tourisme, aviation, marine marchande.
L'augmentation des prix du pétrole et la perspective d'un épuisement de la ressource ont incité les gros
utilisateurs de carburants à demander à leurs états des aides financières ou une défiscalisation.

1.8.4. Indice de cétane et d’octane

L'indice de cétane évalue la capacité d'un carburant à s'enflammer sur une échelle de 0 à 100. Il est
particulièrement important pour les moteurs Diesel où le carburant doit s'auto-enflammer sous l'effet de la
compression. Un carburant à haut indice de cétane est caractérisé par sa facilité à s'auto-allumer.

L'indice de cétane est au moteur Diesel ce que l'indice d'octane est au moteur à essence. Toutefois, ils
décrivent des qualités de carburant radicalement opposées, adaptées au type de moteur. Si le premier
décrit une capacité d'auto-inflammation recherchée par le moteur Diesel, le second décrit une capacité de
résistance à la détonation recherchée par le moteur à essence.

La définition a été posée par Wilfrid de Nercy, ingénieur chez Shell, en 1958.

On dit qu'un carburant a un indice de cétane de x, lorsque celui-ci a le même pouvoir d'auto-inflammation
qu'un mélange constitué de x% en volume de n-cétane et (100-x)% d'alpha-méthylnaphtalène.

Un additif pour carburant, le Tétranitrométhane peut améliorer l'indice de cétane du fuel-diesel.

Échelle
 0 : alpha-méthylnaphtalène (C10H7CH3) ;

16 
 
  
 
 15 : isocétane (2,2,4,4,6,8,8-heptaméthylnonane) (iso-C16H34) utilisé comme nouvelle référence à
la place de l'alpha-méthylnaphtalène ;
 40 : fioul domestique ;
 51 : gazole ;
 85 - 96 : diéthyl éther C4H10O ;
 100 : n-cétane (C16H34).

1.9 Notions sur la chimie des hydrocarbures

Les différentes familles des hydrocarbures sont distinguées par :
- Les liaisons entre les atomes de carbone carbone-carbone, (C-C) en simple liaison dites
alcanes ou (C=C) en double dites alcènes ou (C≡C) en triple liaisons dites aussi alcynes.
- Aussi au fait que les molécules sont arrangées en chaines ouvertes (toutes les limites des
chaines sont déconnectées) ou fermées aussi dites anneaux ou cycliques.
Les alcanes, alcènes et alkyles sont des structures de chaines ouvertes, tandis que les cyclanes et
les aromatiques ont des structures annulaires.
Pour les familles des chaines ouvertes, la nomenclature suivante est utilisée pour démontrer le
nombre d’atomes de carbone contenu dans un membre d’une famille particulière.
Nb Dénomination Nb Dénomination Nb Dénomination
atomes atomes atomes
1 Méth 7 Pent 9 Non
2 Eth 6 Hex 10 Dec
3 Prop 7 Hep 11 Undec
4 But 8 Oct 12 Dodec

Nom de la famille Formule moléculaire Liaison C-C

Alcanes Liaison simple
Alcènes Une seule double
Alcynes Une seule triple

C1 C2 C3 C4
Alcanes Méthane Ethane Propane Butane
Alcènes Ethène Propène Butène
Ethylène Propylène Butylène
Alcynes Ethyne Propyne Butyne
Acetylene Methyl acetylene Ethyl acetylene

17 
 
  
 

Ci-dessous quelques exemples d’hydrocarbures :

Alcanes

Méthane  Ethane  Propane Butane 

Alcènes

Ethène  Propène  Butène 

                                                     

Alcynes

Ethyne  propyne  Butyne 

18 
 
  
 

Chapitre 2

Mélanges de réactifs et produits de combustion

2.1 Stœchiométrie

La quantité stœchiométrique d’un oxydant (l’air par exemple) est la quantité nécessaire

pour bruler complètement une quantité de combustible. Si on alimente par une quantité

supérieure d’oxydant, le mélange est dit pauvre. Si par contre on alimente par une quantité

inférieure, le mélange est dit riche. L’oxydant (où l’air) stœchiométrique est déterminé par

l’écriture d’un simple bilan atomique, en supposant que le combustible réagit pour donner les

produits. Pour un hydrocarbure , la relation stœchiométrique peut être exprimée par :

3.76 → (2.1)

Pour déterminer le coefficients stœchiométriques , et , on fait un bilan atomique

pour chaque atome à savoir C, H, N, O. Commençons par le carbone, à gauche de la réaction on

a atomes, à droite donc . En procède de la même façon on trouve :

 Bilan sur : 2 donc /2

 Bilan sur : 2 2 2 /2 onc /4

 Bilan sur N : 3.76

On trouve la réaction équilibré de l’hydrocarbure général avec l’air :

3.76 → 3.76 (2.2)

Avec (2.3)

On a supposé que l’air est composé de 21% de et 79% de en volume, c.à.d. pour chaque

mole de , il y’a 3.76 moles .

19 
 
  
 

2.1 Rapport stœchiométrique, richesse et excédant d’air

Le rapport stœchiométrique est défini par le rapport de la masse ou débit d’air à celui du fuel,

on écrit :

(2.4)

La masse d’air peut être déduite de la réaction (29)

4.76 (2.5)

Celle du fuel, puisqu’on une seule mole :

1∗ (2.6)

En remplaçant par on obtient :

4.76 (2.7)

Le tableau suivant donne les valeurs de ce rapport pour quelques réactions :

17.11
8.00
11.4

Il faut noter que dans la deuxième réaction l’oxydant est l’oxygène pur, sa masse molaire
est différente de celle de l’air.

La richesse est un facteur très important dans l’étude de la performance des systèmes. Elle est

définie par le rapport des rapports stœchiométriques à l’état stœchiométrique et réel. Ce

paramètre indique si le mélange est riche en combustible où pauvre :

(2.8)

20 
 
  
 

1→ é
Pour 1→ é
1→ é é

L’excédent d’air est l’inverse de la richesse, il est définit par : (2.9)

Deux autres paramètres sont utilisés pour définir la stœchiométrie relative et l’air

stœchiométrique, le pourcentage d’air stœchiométrique et le pourcentage d’excès d’air :

%
% (2.10)

% é ′ 100% (2.11)

2.2 Enthalpie et enthalpie de formation absolue (standard)

Dans les systèmes réactifs, le concept des enthalpies standard où absolues est

extrêmement important. Pour chaque espèce, on définit une enthalpie standard où absolue qui est

la somme de l’enthalpie de formation, liée à la nature de l’espèce, et l’enthalpie sensible, Δ

associée à la variation de la température. On écrit pour l’enthalpie molaire absolue de l’espèce

« i »:

, Δ , (2.12)

Pour utiliser cette relation en pratique, il est nécessaire de définir un état standard où de

référence indiquée par « ref ». On utilise une température, Tref=25°C (298.15K), et une pression,

Pref=P0=1 atm (101325 Pa). En plus, on adopte la convention d’enthalpie de formation nulle pour

les éléments existant dans leurs états naturels à l’état de référence de T et P.

Par exemple, à 25°C et 1 atm, l’oxygène existe sous forme de molécule diatomique,

alors , 0.

21 
 
  
 

Pour former des atomes d’oxygène monoatomique à l’état standard, il faut rompre une

liaison chimique très forte. L’énergie de dissociation de la liaison d’O2 à 298 K est 498390

Kj/kmole de O2. Cette rupture donne deux atomes O. Alors, l’enthalpie de formation de

l’oxygène atomique est la moitié de ce celle de dissociation de l’ O2, , 249195 /

Les enthalpies de formation ont une interprétation physique très claire, c’est la variation

d’enthalpie associée à une rupture des liaisons chimiques d’un élément aux conditions standard

pour former un nouveau composé.

  La représentation graphique de

l’enthalpie absolue permet de comprendre

et d’utiliser ce concept. A 298.15 K, on

note que , 0 (par définition) et

l’enthalpie absolue de O est égale à

l’enthalpie de formation, puisque celle

sensible est nulle à 289.15K. A 4000K, on

remarque une contribution additionnelle de

l’enthalpie sensible dans celle absolue.

Les enthalpies de formation à l’état de référence sont données. Celles sensibles pour

plusieurs espèces sont données en fonction de la température sous forme de tableaux ou de

polynômes. 

22 
 
  
 

2.3 Enthalpie de combustion (chaleur de réaction)

Considérons un réacteur stationnaire parfait dans

lequel en injecte un mélange stœchiométrique en l’existence

des produits de réaction à l’état standard (25°C et 1atm). Le Réactifs    Produits  

Stoechio.  De combustion  
procédé de combustion est complet, le fuel est complètement

consommé pour donner le CO2 et l’H2O. Pour que les produits qcv 
et le gaz injecté soient à la même température, il faut évacuer

la chaleur en permanence du réacteur.

Cette quantité de chaleur est reliée aux enthalpies absolues des réactifs et produits

par (28) en l’absence du travail et en négligeant les vitesses et les hauteurs:

(2.13)

La définition de l’enthalpie de réaction où de combustion Δ (par unité de masse du mélange),

est Δ (2.14)

En termes de propriétés extensives, on écrit :

Δ (2.15)

Prenons par exemple la combustion du

méthane avec l’air à 25°C et 1 atm. L’enthalpie

des réactifs stœchiométriques du CH4 et l’air pour 1

kmole du fuel qui réagit est -748831 KJ. Aux mêmes

conditions (25°C et 1 atm), les produits de

combustion ont une enthalpie absolue de -877236 J.

L’enthalpie absolue de combustion des produits se

trouve au dessous de celle des réactifs.

L’enthalpie de combustion peut être calculée par :

Δ 877236— 74831

802405 .
23 
 
  
 

2.4 Facteur de pondération de l’enthalpie

Si on veut pondérer la quantité d’énergie ou l’enthalpie par unité de masse du fuel ou de

mélange, on utilise la relation :

Δ é é Δ (2.16)
é

Avec : (2.17)
é

Dans le calcul précédent on a :

Δ 50016
.

Pour le méthane (A/F) = 17.11 ce qui donne Δ 2761.8

é .

On note que la valeur de l’enthalpie de combustion dépend de la température choisie pour

son évaluation car les enthalpies des réactifs et des produits varient avec la température. La

distance entre les courbes HProd et Hreac n’est pas constante.

Pour un hydrocarbure , l’enthalpie de réaction (combustion) est donnée par :

Δ , , , (2.18)

2.5 Pouvoirs calorifiques

Le pouvoir calorifique supérieur, Pcs, est la chaleur de combustion calculée, par

kilogramme de combustible, en supposant que l’eau dans les produits de combustion est

condensée en liquide. Ceci libère le maximum d’énergie. Le pouvoir calorifique inférieur, Pci,

correspond au cas où l’eau n’est pas condensée et quitte la chambre sous sa forme vapeur. Pour

le méthane CH4 le Pcs est de 11% supérieur au Pci. Les valeurs des pouvoirs calorifiques sont

données sous forme de tableaux pour différents combustibles.

24 
 
  
 

Généralement, les pouvoirs calorifiques sont donnés en Kj/Kg de fuel, ils sont exprimés

en fonction des chaleurs (enthalpies) de réaction par :

(2.19)

A son tour le pouvoir calorifique supérieure s’écrit :

(2.20)

Exemples

Considérons un mélange stœchiométrique d’iso-octane avec l’air, déterminer :

1. La fraction molaire du fuel,
2. Le rapport A/F,
3. La fraction de dans les produits,
4. L’enthalpie de combustion,
5. Le pouvoir calorifique inférieur,
6. Le pouvoir calorifique supérieur.

Solutions

D’après la réaction globale stœchiométrique et complète, celle de l’iso-octane avec l’air s’écrit :

12.5 3.76 →8 9 3.76 ∗ 12.5

1. La fraction molaire du fuel est donnée par :

0.0165
. . ∗ .
2. Le rapport A/F est donné par :

4.76 ∗ 12.5 15.16

3. La fraction molaire de dans les produits est :

0.141
. ∗ .

4. La chaleur de réaction ∆ est donnée par :

∆ , , ,
8 393546 9 241845 1 208447 5116526
5. Le pouvoir calorifique inférieur est :

44882
6. Le pouvoir calorifique supérieur est :
44882 9∗ 48356.5

25 
 
  
 

Chapitre 3

Température adiabatique de la flamme

On définit deux températures adiabatiques de la flamme, une à pression constante et

l’autre à volume constant.

3.1 Température adiabatique de la flamme à pression constante

Commençons par les systèmes ou la combustion se fait à pression constante. Le premier

principe dicte que l’enthalpie absolue des réactifs à l’état initial (Ti=298K, P=1atm) est égale à

celle des produits à l’état final (T=Tad, P=1atm), c.à.d. l’application de (28) donne :

, , (3.1)

Où par unité de masse du mélange , , (3.2)

Cette définition est illustrée graphiquement

par la figure ci-contre. L’évaluation de la

température adiabatique de la flamme nécessite la

connaissance des produits de combustion. A la

température de la flamme, les produits se dissocient

et le mélange se compose de plusieurs espèces. La

température de la flamme est de l’ordre de quelques

milliers de Kelvin.

Pour un combustible , la réaction de combustion complète stœchiométrique pour un

mélange pauvre s’écrit :

3.76 ⇒ 3.76 (3.3)

26 
 
  
 

Si on injecte les réactifs à la température , l’enthalpie des réactifs s’écrit :

, é ∑ ̅ , (3 .4)

L’enthalpie des produits s’écrit à son tour :

, , ∑ ̅ , (3.5)
é

Ici les indices et parcourent les réactifs et les produit simultanément.

Le problème donc réside dans la détermination des intégrales car la température

adiabatique est le paramètre cherché. Pour cela il faut soit connaitre la dépendance des chaleurs

spécifiques de la température soit de supposer que ces dernières sont constante et de les estimer à

une température prescrite. Les chaleurs spécifiques sont approximées par des fonctions

polynômiales pour des intervalles bien définis, si on utilise ces fonctions on doit faire recours

aux méthodes numériques pour déterminer la température adiabatique de la flamme. Pour

simplifier les calculs, on suppose que les ̅ , sont constantes et on prescrit une température é

pour les évaluer, on peut donc écrire :

, , é ∑ ̅ , é (3.6)

En égalisant et , on peut calculer la température adiabatique de la flamme :

, , , é ∑ ̅ ,
é ∑ ̅ ,
(3.7)
é

On voit que le premier terme entre ( ) correspond à l’enthalpie de réaction du fuel

considéré Δ , on peut donc écrire la température adiabatique de la flamme en fonction du Pci

du combustible :

é ∑ ̅ ,
é ∑ ̅ ,
(3.8)
é

27 
 
  
 

Si on injecte les réactifs à la température de référence, cette relation se simplifie à :

é ∑ ̅ ,
(3.9)
é

Cette formule peut être généralisée pour l’injection d’un mélange de plusieurs fuels avec

l’air. Si les fuels ont les pouvoirs calorifiques et les masses molaires , on écrit :

∑ ,
é ∑ ̅ ,
(3.10)
é

3.2 Température adiabatique de la flamme à volume constant

Pour la température adiabatique de la flamme à volume constant, ce qui est le cas des

moteurs à combustion interne car la combustion est très rapide ce qui permet de supposer que le

volume est constant au moment de la combustion. La première loi de la thermodynamique (1.35)

s’écrit :

, , (3.11)

Puisque les tables thermodynamiques donnes les enthalpies et non pas les énergies internes on

utilise la relation U=H-PV pour écrire :

0 (3.12)

On utilise ensuite l’équation d’état du gaz pour éliminer le terme PV.

∑ (3.13)

∑ (3.14)

Ce qui donne 0 (3.15)

é
On peut écrire cette équation par unité de masse, reconnaissons que et
é
on obtient :

0 (3.16)

28 
 
  
 

Puisque la composition du mélange de produits à l’équilibre dépend de T et P ,

l’utilisation de (39) où (40) avec une loi pour le gaz parfait et une équation d’état appropriée,

, pour une gaz parfait, donne directement Tad.

L’équation (40) s’écrit : (3.17)

Le terme à gauche est similaire de celui déjà développé, on a donc :

, é ∑ ̅ , (3.18)

, , é ∑ ̅ , é (3.19)

Après réarrangement, on trouve :

é ∑ ̅ , ∑ ̅ , ∑ ̅ ,
∑ ̅ ,
(3.20)

Si é la température adiabatique s’écrit :

é ∑ ̅ ,
∑ ̅ ,
(3.21)

En général, le nombre de moles des réactifs et des produits sont très proches, on peut faire

l’hypothèse on écrit :

é ∑ ̅ ,
(3.22)

En comparant cette relation à celle déjà déduite pour ce cas de la pression constante, on
remarque que :

(3.23)

29 
 
  
 

Exemple :

On reprend la combustion de l’iso-octane où on injecte les réactifs à la température de référence.

Aussi les chaleurs spécifiques des produits sont évaluées à 2000K. Calculer la température

adiabatique de la flamme si :

1. La pression est constante,

2. Le volume est constant, calculer la pression en fin de combustion si celle initiale est 1
atm.

Solution :

La réaction de l’iso-octane avec l’air s’écrit :

12.5 3.76 →8 9 3.76 ∗ 12.5

1. Pour une pression constante, la température adiabatique de la flamme est donnée par :

é ∑ ̅ , é

Le est déjà calculé dans l’exemple précédent, 44882 , dans ce cas

é 2000 . Les ̅ , é sont obtenues des tables :

̅ , 2000 60.433, ̅ , 2000 51.143, ̅ , 2000 35.988 /

Le remplacement donne :

44882 ∗ 114
298 2240
8 ∗ 60.433 9 ∗ 51.143 3.76 ∗ 12.5 ∗ 35.988

2. Pour un volume constant, la température adiabatique de la flamme est :

é ∑ ̅ ,
∑ ̅ ,

30 
 
  
 

On calcule 60.5 et 64, cela donne

44882 ∗ 114 298 8 ∗ 60.433 9 ∗ 51.143 3.76 ∗ 12.5 ∗ 35.988 8.314 ∗ 60.5
2735
8 ∗ 60.433 9 ∗ 51.143 3.76 ∗ 12.5 ∗ 35.988 8.314 ∗ 64

Si on suppose que , on obtient 2730 , cette hypothèse est valable pour ces calculs

car l’écart absolue ne dépasse pas 2 pour mille.

Pour la pression en fin de combustion on applique l’équation d’état du gaz parfait au début et à la
fin de la combustion,

La division des deux équations membre par membre donne,

L’application numérique donne

2735
1 9.2
298

31 
 
  
 

Chapitre 4

Equilibre chimique

Dans les procédés de combustion à hautes températures, les produits de combustion ne

sont pas des mélanges simples de gaz idéals comme on a déjà vu dans les bilans atomiques. En

réalité, les espèces majoritaires se dissocient produisant d’autres espèces minoritaires. Sous

certaines conditions, ces espèces minoritaires sont présentes en grandes quantités. Par exemple,

la combustion idéale d’un hydrocarbure dans l’air donne les espèces suivantes : CO2, H2O, O2 et

N2. La dissociation de ces espèces et la réaction entre les nouvelles espèces donne : H2, OH, CO,

H, O, N, NO, et d’autres. Dans cette partie, on va calculer les fractions molaires de ces espèces à

une température et pression donnée. Cela est possible, en considérant la conservation des

éléments c.à.d. le nombre des atomes de C, H, O et N indépendamment de la façon avec laquelle

ils se combinent. Ils y’a plusieurs méthode de calculer l’équilibre chimique, on va considérer

l’approche basée sur la constante d’équilibre.

4.1 Deuxième loi de la thermodynamique

L’équilibre chimique trouve ses racines dans la deuxième loi de la thermodynamique.

Prenons un volume constant dans lequel une masse constante de réactifs en réaction adiabatique

forme des produits. Au cours de la réaction, la température et la pression augmentent jusqu'à

l’équilibre final. La première loi de la thermodynamique ne peut pas prédire toute seule l’état

final (température, pressure et composition), cela nécessite la deuxième loi. Pour illustrer le

problème considérons l’exemple de la combustion du CO.

→ (4.1)

32 
 
  
 

Si la température est suffisamment importante, le CO2 va se dissocier. Supposons que les

produits sont seulement le CO2, CO, et O2, on peut écrire :

→ 1 (4.2)
é

Où α est la fraction de CO2 dissociée. On peut calculer la température adiabatique de la

flamme en fonction de la fraction dissociée α. Par exemple, pour α=1, aucune chaleur n’est

dégagée et la température, la pression et la composition restent inchangée. Lorsque α=0, le

maximum de chaleur est dégagée et la pression et la température doivent avoir leurs maximum

permis par la première loi de la thermodynamique, cette variation est montrée par la figure.

Calculons maintenant l’entropie du mélange par sommer les entropies des différentes

espèces :

é , ∑ ̅ , 1 ̅ ̅ ̅ (4.3)

Les entropies individuelles des espèces sont calculées par :

̅ ̅ ̅ , (4.4)
é

33 
 
  
 

En traçant l’entropie du mélange en fonction de la fraction dissociée α, on voit qu’un

maximum est atteint pour une valeur intermédiaire de (1- α = 0.5) pour la réaction prise en

compte.

Dans notre exemple, on supposé que U,V et m sont constantes ce qui implique qu’il n’y a

pas d’échange de chaleur où de travail. Selon la deuxième loi la variation de l’entropie du

système s’écrit :

0 (4.5)

D’où la composition du système va spontanément tendre vers le point maximal d’entropie

de part et d’autre car dS est positive. Lorsque le maximum est atteint, aucune variation n’est

permise, cela rentre en contradiction avec le deuxième principe (dS<0). Formellement, la

condition de l’équilibre peut s’écrire :

, , 0 (4.6)

En résumé, si on fixe l’énergie interne U, le volume V et la masse m d’un système isolé,

l’application de la première et deuxième loi et l’équation d’état vont définir la température

d’équilibre, la pression et la composition chimique.

4.2 Fonction de Gibbs ou énergie libre

En réalité l’utilisation d’un système isolé d’une masse constante et un volume fixe n’a

aucun intérêt pratique. Par exemple, il est fréquent de calculer la composition d’un mélange à

une température, une pression et une stœchiométrie donnée. Pour ce problème, l’énergie libre

de Gibbs, G, remplace l’entropie comme la propriété thermodynamique la plus importante.

L’énergie libre de Gibbs est définie en terme d’autre propriétés thermodynamique par :

(4.7)

34 
 
  
 

La deuxième loi est exprimée en fonction de l’énergie libre de Gibbs par :

, , 0 (4.8)

Cette fonction diminue toujours pour une évolution spontanée d’un système à masse

constante isotherme, isobare et sans travail sauf celui des limites (P-dV). Ce principe permet de

calculer la composition à l’équilibre d’un mélange à une température et pression données. La

fonction de Gibbs atteint un minimum à l’équilibre, en contraste au maximum atteint pour

l’entropie.

Alors, à l’équilibre on a : , , 0 (4.9)

Pour un mélange de gaz idéals, la fonction de Gibbs pour l’espèce « i » est donnée :

̅, ̅, ⁄ (4.10)

Où ̅ , est la fonction de Gibbs de l’espèce pure à la pression de l’état standard (c.à.d. Pi=P0) et

Pi est la pression partielle. La pression à l’état standard, P0 dans le dénominateur, est prise par

convention 1 atm. Dans les systèmes réactifs, il est nécessaire d’utiliser la fonction de formation

de Gibbs,

̅ , ̅ ∑é é ̅ (4.11)

Où sont les coefficients stœchiométriques des éléments nécessaires pour former une mole de

l’espèce en question.

La fonction de Gibbs pour un mélange pur de gaz idéals est exprimée par :

é ∑ ̅, ∑ ̅, ⁄ (4.12)

Pour une température et pression constante, les conditions d’équilibre deviennent

é 0 (4.13)

35 
 
  
 

Où ∑ ̅, ⁄ ∑ ̅, ⁄ 0 (4.14)

Le second terme dans l’équation (4.14) est nul car :

̅, et / et ∑ ∑ 0

Cela donne :

é 0 ∑ ̅, ⁄ (4.15)

Pour une réaction qui s’écrit sous la forme générale :

⋯⇔ ⋯ (4.16)

Ou bien 0 0 ⋯⇔ 0 0 ⋯

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

coefficient stœchiométrique, c.à.d.

é 0 ,

é 0 ,

é 0 ,

é 0 ,…… (4.17)

En substituant dans (4.15) et on éliminant la constante de proportionnalité , on aura :

̅ , ⁄ ̅ , ⁄ ⋯

̅ , ⁄ ̅ , ⁄ ⋯ 0

36 
 
  
 

Après réarrangement en regroupant les termes ensemble

⁄ ⁄ . .
̅ , ̅ , ⋯ ̅ , ̅ , ⋯ ⁄ ⁄ . .

Le terme entre parenthèses à gauche de l’équation est dit la variation de la fonction de Gibbs de

l’état standard Δ c.à.d.

Δ ̅ , ̅ , ⋯ ̅ , ̅ , ⋯

Où alternativement

Δ ̅ , ̅ , ⋯ ̅ , ̅ , ⋯ (4.18)

L’argument du logarithme est définie par constante d’équilibre KP pour la réaction, elle est notée
par :

⁄ ⁄ . .
⁄ ⁄
(4.19)
. .

De ces définitions de Δ et , la condition d’équilibre chimique à une température et pression

constante est donnée par :

Δ (4.20)

Où bien exp Δ / (4.21)

Si Δ est positive, les réactifs vont être favorisés car est négative, ce qui implique

que elle-même est inférieure à l’unité c.à.d.

⁄ ⁄ . . ⁄ ⁄ . ..

Similairement, si Δ est négative la réaction favorise les produits. Une interprétation physique

peut être faite si on écrit la définition de Δ en fonction de l’enthalpie et l’entropie :

37 
 
  
 

Δ Δ Δ (4.22)

Qui peut être substituée dans l’équation de .

exp Δ / exp Δ /

Pour que soit supérieure à l’unité, ce qui favorise les produits, la variation de

l’enthalpie de la réaction, Δ , doit être négative, c.à.d., la réaction est exothermique et

l’énergie du système diminue. Aussi, une variation positive de l’entropie, qui indique un chaos

moléculaire plus important, donne des valeurs de 1.

Exemple :

On veut étudier l’effet de la pression sur la réaction d’équilibre ↔ à T = 2000K.

Pour cela on prend les pressions P =1, 3, 5, 8 et 10 atm.

1. Ecrire l’expression de la constante d’équilibre,

2. Calculer la constante Kp,

3. Ecrire les équations qui donnent la composition du mélange si on suppose que le

rapport atomique C/O = ½.

4. Calculer la composition du mélange pour les pressions P =1, 3, 5, 8 et 10 atm.

Solution :

1. La constante d’équilibre s’écrit :

⁄
⁄ ⁄ /

2. La constante d’équilibre peut être calculée par exp Δ /

Avec Δ 1 ̅ , 1 ̅ , ̅ , 396410 285948 0 110462 Kj/Kmole

38 
 
  
 

D’où exp 110462/ 8.314 ∗ 2000 767.5

3. Pour le calcul de la composition du mélange, on doit avoir trois équations car on a trois

inconnues à savoir : , et .

La première équation est

=1 (1)

De l’hypothèse C/O = ½, on sait que CO donne une mole de O et de C, CO2 donne deux moles

de O et une de C et finalement O2 donne deux moles de O. Ecrivant le rapport

1
2 2 2 2 2

On obtient 0 (2)

La troisième est obtenue de la constante d’équilibre

⁄
⁄ ⁄ / / /

D’où 0 (3)

Les équations (1), (2) et (3) donnent la composition du mélange à la température et pressions

voulues. C’est un système non linéaire et couplé, dans l’équation (3) on doit éliminer et

pour chercher par une méthode numérique.

De (2) on a 2 , le remplacement dans (1) donne 1 3 .

Si on remplace dans l’équation (3), on obtient 1 3 2 0, en réarrangeant cela

/
donne

39 
 
  
 

La convergence de cette forme d’écriture est vérifiée pour la méthode des approximations

successives. Si on part d’un estimé initial 0.000 en écrivant :

/
n =0,1,2,… (n est le numéro de l’itération)

P = 1 atm

Pour n=0 7.405 10 , n=1 7.407 10 , n=2 7.407 10

Donc 7.407 10 , 9.778 10 et 1.481 10

P = 3 atm

Pour n=0 5.208 10 , n=1 5.154 10 , n=2 5.155 10

Donc 5.155 10 , 9.845 10 et 1.031 10

P = 5 atm

Pour n=0 4.348 10 , n=1 4.348 10

Donc 4.348 10 , 9.870 10 et 8.696 10

P = 8 atm

Pour n=0 3.759 10 , n=1 3.731 10 , n=2 3.731 10

Donc 3.731 10 , 9.888 10 et 7.462 10

P = 10 atm

Pour n=0 3.497 10 , n=1 3.472 10 , n=2 3.472 10

Donc 3.472 10 , 9.990 10 et 6.944 10

40