Maamri, Rachid
Maamri, Rachid
THÈSE PRÉSENTÉE À
L’UNIVERSITÉ DU QUÉBEC À TROIS-RIVIÈRES
COMME EXIGENCE
DU DOCTORAT EN INGÉNIERIE
Par
Rachid Maamri
MÉLANGES
Octobre 2014
ii
RÉSUMÉ
microprocesseur de ce même moteur. Une approche analytique basée sur le modèle de Wiebe
a été élaborée et des tests sur un moteur à petite cylindrée fonctionnant avec le mélange du
gaz naturel comprimé et d'hydrogène ont été réalisés.
X X
Yves Dubé Lotfi Toubal
Directeur de thèse Co-directeur de thèse
iv
ABSTRACT
The inventor of the compression ignition engine, German Rudolf Diesel was
motivated by the desire to replace the coal fueled steam engine by an engine running on fossil
fuels with a higher efficiency and also much cleaner. These considerations remain topical
today more than ever. The will of legislator to ensure environmental compliance has forced
the auto industry to focus on other alternate cleaner source of energy or other than
hydrocarbons and also on engine optimization of our vehicles.
While we await the generalization of even cleaner vehicles (hydrogen fueled engines,
fuel cell, electric vehicle ... etc.), engines fueled by substitute fuels represent an interesting
alternative to those operating on conventional fuels (gasoline and diesel).
The objective of this thesis is to evaluate the potential of natural gas and hydrogen
in pure form or mixed with other fuels, in internal combustion engines, as well as the potential
of the external combustion engine.
The research topic addressed in this paper marks the beginning of a collaboration
between three research teams from three universities: the Université du Québec à Trois-
Rivières, Canada, the National Technical University “Kharkov Polytechnic Institute”,
Ukraine and Kharkov National Automobile and Highway University, Ukraine. During my
first training course held at the Kharkov Polytechnic Institute, the research focused on the use
of pure hydrogen in internal combustion engines and the development of an optimization
model for an external combustion engine. A mathematical model of the combustion-
expansion phase in a hydrogen fueled engine was developed and tests on a new concept of
external combustion engine were carried out. In the second training course, which was held at
the Ukraine and Kharkov National Automobile and Highway University, engine experimental
aspect was the guideline parallel to the first training course. Thus, tests and calculations were
carried out for an engine running on pure natural gas and another fueled by a mixture natural
gas / hydrogen. A method of calculation of the equilibrium composition of the combustion
products of the gas engine using the Zeldovitch-Polarny method was developed and compared
with the experimental data. Experimental tests of a gas engine converted from a diesel engine
have been carried out and a mathematical model based on the Wiebe model with the accurate
calculation of the parameter of burning m was developed. We also carried out tests on a
natural gas engine with constant addition of hydrogen and developed maps for the
management system with microprocessor of the same engine. An analytical approach based
on the Wiebe model was developed and tests on a small engine running on a mixture of
compressed natural gas and hydrogen has been achieved.
v
This thesis is the fruit of works resulting from this partnership and it is the foundation
of a promising research axis which will certainly be continued in the future.
vi
Remerciements
Je tiens à remercier le professeur Yves Dubé ainsi que le professeur Loutfi Toubal du
département de Génie mécanique de l’UQTR, qui ont dirigé cette étude, de m’avoir assisté
financièrement et moralement, de m’avoir donné l’opportunité de réaliser ce travail de thèse
et de m’avoir soutenu au cours de cette thèse.
Je tiens à remercier également tous les membres de jury qui ont accepté de juger ma
thèse :
vii
Je dédie ce travail à :
-Toute ma famille et surtout mon unique fille Soumia, car sans leur support et sans leur
confiance je n’aurais jamais réussi à finaliser ce travail !
-Benaceur Benchohra
-Tribu de Larbaa
viii
Résumé ................................................................................................................................ ii
Abstract .............................................................................................................................. iv
Remerciements ................................................................................................................... vi
Nomenclature .................................................................................................................... xx
Chapitre 1........................................................................................................................... 1
Introduction ................................................................................................................. 1
1.2 Problématique..................................................................................................... 2
Chapitre 2........................................................................................................................... 7
Chapitre 3......................................................................................................................... 21
3.9 Références......................................................................................................... 86
Partie II. Carburants de substitution à l’état pur dans les MCI .............................. 101
4.5 Méthode de calcul des émissions des oxydes d'azote .................................... 111
4.6 Comparaison des résultats des calculs avec les données expérimentales ...... 113
5.4.1 Choice of the compression ratio and the maximum charging pressure ... 123
Partie III. Carburants de substitution en mélange dans les MCI ............................. 147
d’hydrogène.............................................................................................................. 147
9.1 Analyse des carburants de substitution à l’état pur dans les MCI.................. 177
9.2 Analyse des Carburants de substitution en mélange dans les MCI ................ 178
Article I:
Calculation method of combustion products equilibrium composition of gas
engine with spark ignition / Maamri R., Abramtshuk F.I., Kabanov A., Dubé
Y., Toubal L., Kodjo A. / / Eastern European Journal of Enterprise
Technologies, Kharkov Ukraine. - 2011- No. 5/6 (53) - 63 p - pp. 37-43.
Article II:
R. Maamri, F. I. Abramtshuk, A. N. Kabanov, M. S. Lipinsky, Y. Dubé, L.
Toubal, A. Kodjo. Investigation in Adjusting the Parameters of a Diesel
Engine Converted to Forced Aspiration Gas Engine, American Journal of
Vehicle Design, 2013, Vol. 1, No. 1, 9-15.
Article III:
R. Maamri, A. P. Martshenko, O. O. Osetrov, Y. Dubé, L. Toubal, A. Kodjo.
Analyze and Mathematical Modeling of the Combustion Process of One-
Cylinder Spark Ignited Hydrogen Fueled, EngineAmerican Journal of
Vehicle Design, 2013, Vol. 1, No. 2, 21-24.
Article IV:
Choosing of performance parameters of light-duty engine running on natural
gas and hydrogen mixture / Maamri R., Abramtshuk F.I., Kabanov A., Dubé
Y., Toubal L., Kodjo A. / / Automotive Transport, Kharkov Ukraine. - 2011.
– No. 29 - pp. 152-159.
Article V:
R. Maamri, W. G. Dyatshenko, A. I. Varonkov, O. U. Linkov, Y. N.
Nikitenko, Y. Dubé, L. Toubal, A. Kodjo, Development of external
combustion engines, American Journal of Vehicle Design, 2013, Vol. 1, No.
2, 25-29.
xv
Figure 2.1 Utilisation de l’hydrogène dans les moteur à combustion interne ............ 15
Figure 3.3 Les émissions spécifiques des HC, NOx et du CO avant le catalyseur, 30%
RGE barre 14 ............................................................................................. 34
Figure 3.4 Les émissions spécifiques des HC, NOx et du CO après le catalyseur, 30%
RGE barre 14 ............................................................................................. 35
Figure 3.10 Délais d’allumage mesurés pour mélange avec 50% d’H2 en fonction de la
pression avec température constante Tc = 995 ± 4 K et φ = 1.0, φ = 0.5 .. 57
Figure 3.12 Evolution des vitesses de propagation de flamme des différents mélanges
méthane/hydrogène/air en fonction de la richesse (Yu et al., [123]) ........ 59
Figure 3.15 Evolution des vitesses en fonction de la richesse pour des mélanges
CH4/H2 dans l’air. Croix: 35% H2, cercles: 25% H2, triangles: 15% H2,
carrés: 5% H2, diamants: 0% H2; traits: modélisation (Coppens, [12]) .... 62
Figure 3.16 Comparaison des vitesse de flammes expérimentales des mélanges CH4,
CH4/10%H2, CH4/20%H2, CH4/40%H2, Tini = 298 K, P = 1 atm, (1: Yu et
al., [123], 2: Halter et al., [130], 3: Coppens et al., [91]) ......................... 65
Figure 3.21 Schéma de l’écoulement à un seul jet dans un plan de stagnation ............ 71
Figure 3.28 Variation des émissions des HC en fonction de la vitesse du moteur ....... 77
Figure 5.7 Emissions variation of the engine 6GCHN13/14 by load (n = 1500 min-1):
– – – – naturally aspiration; –––– – forced aspiration ............................. 131
Figure 6.2 Selection of Wiebe exponent m (a) and burn duration φz (b) in formulas (1)
and (2) for a regime when n = 600 min-1 and Ne = 1.3 kW. By points the
measured diagram is shown. .................................................................. 142
Figure 6.3 Combustion speed (а) and pressure- crank angle diagram (b) with Wiebe
exponent m. Burn duration φz = 22 °crank angle. By points the CFR engine
experimental diagram is shown. n = 600 min-1, Ne =1.3 kW.. ................ 143
Figure 6.4 Influence of a spark advance θоп (a), compression ratio ε (b) and air–fuel
excess ratio α (c,d) on measured and calculated pressure-crank angle
diagrams. By points experimental diagrams are shown .......................... 144
Figure 7.1 Système d’alimentation en combustible sur banc d’essai. ...................... 150
Figure 7.3 Influence d’ajout d'hydrogène sur l'énergie d'inflammation du mélange gaz
naturel et air pour divers coefficients d’excès d'air ................................. 156
Figure 7.8 Variation de la concentration des composants toxiques selon CEVM.... 160
Figure 8.1 Schematic view of the external combustion heat engine ........................ 168
Figure 8.2 Variation in the charge parameters in the combustion chamber in the zone
above the piston with crank angle ........................................................... 172
Figure 8.3 Cylinder pressure external combustion engine (n = 5000 min-1, air pressure
in the zone-buffer ps = 5 MPa and the air to fuel ratio in the combustion
chamber α = 3) ........................................................................................ 173
xviii
Figure 8.4 Influence of the pressure of combustion products ps in the zone located above
the piston on the indicated parameters of the external combustion engine at
α = 3......................................................................................................... 174
Figure 8.5 Influence of the variation of the air to fuel ratio α in the combustion chamber, on
the indicated parameters of the external combustion engine at ps = 5 МPа ..... 175
xix
Tableau 3.3 Les résultats obtenus avec le nouveau concept du moteur contre ceux du
moteur à gaz prototype .............................................................................. 36
Tableau 3.5 Récapitulatif des conditions expérimentales des travaux relatifs aux
vitesses de flammes méthane/hydrogène .................................................. 63
Tableau 5.2 The results of tests with forced aspiration engine 6GCHN13/14 with the
ESC test cycle.......................................................................................... 134
Tableau 5.3 The results of tests naturally aspirated engine 6GCH13/14 with the ESC
test cycle .................................................................................................. 134
Tableau 7.1 Intervalle de variation des paramètres du moteur lors de la résolution des
problèmes d'optimisation des paramètres- régulateurs du moteur .......... 155
Tableau 7.2 Résultats des tests réalisés sur le moteur selon le cycle ESC.................. 161
xx
NOMENCLATURE
Abréviations françaises
Acronymes anglo-saxons
Chimie
Ar Argon
CH4 Méthane
C2H6 Éthane
C3H8 Propane
C4H10 Butane
CO2 Dioxyde de carbone
H2 Hydrogène
HC Hydrocarbures imbrûlés
HFC Hydrofluorocarbures
GN Gaz naturel
N2O Protoxyde d'azote
OH Radical Hydroxyle
Indices
ad Avance
av Advance
С Carbone
c Suralimentation (compresseur)
Car Carburant
d Détonation
e Effectif
H Hydrogène
НС Hydrogène et carbone
i Espèce chimique
ini Initiale
l Laminaire
N Azote
t Turbuente
th, frein Théorique au frein
tol Toléré
pap Papillon
p.с Produits de combustion
var Variable
Paramètres géométriques
Mécanique
Thermodynamique et combustion
Chapitre1
Introduction
1.2 Problématique
Le champ de recherche concernant les carburants de substitution est assez jeune, mais
l'intérêt qu'il suscite se fait sentir dans de plus en plus de pays touchés par le problème de la
pollution issue des gaz à effet de serre. Depuis environ une trentaine d'années, un nombre
important de laboratoires se sont intéressés à la problématique des carburants de
remplacement. Ces études ont permis d’identifier le gaz naturel, le méthane, le gaz du pétrole
liquéfié (GPL) et l’hydrogène comme étant des carburants susceptibles d’être utilisés dans le
domaine automobile et qui peut remplacer les carburants conventionnels telles que: l’essence
et le carburant Diesel, afin de réduire les émissions des gaz à effet de serre et améliorer la
qualité de l’air.
La modélisation mathématique est considérée comme l’un des instruments les plus
efficaces pour étudier les cycles dans les moteurs à combustion interne. Lors de l’élaboration
des modèles mathématiques, plusieurs questions surgissent. Ces questions ne sont pas encore
traitées dans la littérature technique. Par exemple, la dynamique de dégagement de chaleur
dans les moteurs à hydrogène, l’influence des paramètres constructifs et des paramètres de
réglage durant le processus de fonctionnement du moteur. Dans ce cadre, l’analyse de la
phase de combustion de l’hydrogène et l’élaboration du modèle mathématique, qui prennent
en considération les particularités de la combustion, l’estimation de l’influence des paramètres
constructifs et des paramètres de réglage du moteur, sont nécessaires.
En raison de la dégradation de la situation écologique dans le monde, beaucoup
d'attention est donnée aux questions de toxicité des moteurs à combustion interne. En
conséquence, le calcul des composants toxiques dans les produits de combustion s'avère ainsi
indispensable. La garantie d'exécution de cette tâche permet d'économiser énormément les
ressources humaines et matérielles à l'aide de la réduction du nombre des tests.
L’ajout d'hydrogène aux divers carburants connus permet non seulement à réduire la
quantité d'hydrogène nécessaire à titre de combustible, mais il permet également d’augmenter
d'une manière significative les propriétés écologiques de l'automobile [1, 2, 3, 4, 5, 6, 7, 8,10,
11,12,13,14].
L'analyse réalisée des sources bibliographiques [59-116] a montré l'opportunité
d'utilisation du mélange gaz naturel-hydrogène (mélange) à titre de combustible moteur (voir
3
§ 3.3-§ 3.5). C'est pourquoi le but de cette étude est le choix des paramètres optimaux du
fonctionnement d’un moteur fonctionnant avec mélange gaz naturel-hydrogène.
La conversion des moteurs Diesel des camions aux moteurs à gaz à allumage par
étincelles est actuellement un sujet d’actualité. L'application de ce concept entraine la chute de
la puissance du moteur, qui s’avère particulièrement considérable aux régimes de la charge
maximale. Pour la compensation des pertes de puissance, qui apparaissant lors de la conversion du
moteur diesel vers le gaz, il sera rationnel d'utiliser la suralimentation. Il en résulte que le choix et le
fondement des paramètres régulateurs du moteur à gaz suralimenté converti à partir du moteur
diesel sont absolument nécessaires.
L’une des solutions efficaces pour remédier aux problèmes d’épuisement imminent
des ressources fossiles et de la flambée de leur prix est le développement d’un nouveau
concept de moteur: le moteur à combustion externe, qui est une source prometteuse d'énergie
mécanique dans le domaine du transport. Il a de bonnes caractéristiques techniques: un
rendement très élevé, la capacité de fonctionner avec différents carburants, une structure
simple et les coûts d'exploitation faibles.
1.3 Objectifs
Dans ce travail, nous avons comme objectif d’étudier l’influence de l’utilisation des
carburants de substitution, tels que le gaz naturel et l’hydrogène à l’état pur et en mélange, sur
les performances des moteurs à combustion interne à allumage par étincelles et à allumage par
compression. De plus, nous étudions la phase de combustion des moteurs alimentés par ces
carburants. Un deuxième objectif de cette thèse concerne l’étude des performances d’un
concept original d’un moteur à combustion externe. Plus spécifiquement, il s'agit :
• d’analyser la phase de combustion des moteurs fonctionnant avec de tels mélanges,
• d'étudier l'influence d’utilisation des carburants de substitution sur les émissions des
composants toxiques dans les gaz d’échappement,
• d'étudier l'influence d’utilisation des carburants de substitution sur les performances
du moteur, telles que : la puissance, le rendement et la consommation du carburant,
• d’optimiser les moteurs alimentés par ces carburants,
• d’étudier le concept de moteurs à combustion externe, du point de vue des
performances, tel que: la pression et la température des produits de combustion à l’intérieur
du cylindre, le rendement, la consommation, l’utilisation de différents carburants...etc.
4
Les résultats ainsi obtenus contribueront à une meilleure compréhension des différents
phénomènes physico-chimiques complexes qui ont lieu dans les moteurs à combustion interne
alimentés par des carburants de substitution et à justifier l’utilisation future de ce type de
carburants et des moteurs à combustion externe.
1.4 Méthodologie
Première étape : utilisation des carburants de substitution à l’état pur dans les moteurs à
combustion interne
Dans cette étape, on a réalisé des calculs et des expérimentations sur des moteurs à
combustion interne fonctionnant seulement avec un seul carburant de substitution.
On a exposé les résultats d'étude d’un moteur à gaz suralimenté, converti à partir d’un
Diesel. Pour ce moteur, on a précisé le modèle mathématique du cycle de combustion basé sur le
modèle de Wiebe avec le calcul précis du paramètre de combustion m.
Deuxième étape : utilisation des carburants de substitution en mélange dans les moteurs à
combustion interne
Dans cette deuxième étape, on a réalisé des études analytiques et effectué des tests sur
des des moteurs à combustion interne alimentés cette fois par différents mélanges de
carburants de remplacement.
5
On a exposé les résultats d'étude d’un moteur à petite cylindrée à gaz fonctionnant
avec le mélange du gaz naturel et d'hydrogène. Pour ce moteur, on a précisé le modèle
mathématique de la combustion avec l’indice variable de combustion de Wiebe. Des tests ont
été réalisés afin de valider les résultats de calcul. En vertu des études réalisées, on a reçu les
cartographies pour les systèmes de gestion du moteur fonctionnant avec le mélange du gaz
naturel et d'hydrogène.
Ce moteur utilise de l'air dans des réservoirs à des pressions de 30 - 50 МPа et à une
température ambiante en tant que fluide de travail et l'hydrogène, les alcools ou les carburants
traditionnels issus de minéraux organiques, comme source d'énergie chimique. Un modèle
mathématique a été développé, testé et vérifié afin de simuler les cycles d'un moteur 4 temps à
allumage commandé fonctionnant à l'essence.
Chapitre 2
Notions fondamentales
Avant d’entamer toute discussion sur les carburants de substitution, il sera souhaitable
de donner des notions de base sur les différents types de moteurs existants fonctionnant avec
différents carburants.
Avantages et inconvénients :
Les moteurs à combustion interne, qui sont relativement petits et légers, compensent
un couple assez faible par une vitesse de rotation élevée. Les carburants conventionnels qui
font fonctionner ce type de moteur sont peu encombrants et faciles à se procurer, ce qui en fait
des moteurs à explosion tout à fait indiqués pour équiper de petits véhicules roulants, mais
aussi des appareils volants. Leur succès est dû aussi à leur facilité d'utilisation et de
maintenance. De plus, ils ne sont pas délicats et peuvent fonctionner sans problème avec
divers carburants, sans qu'il soit nécessaire de procéder à des modifications majeures.
L’alcool ou le gaz peut remplacer l’essence, quant au carburant diesel, il peut être remplacé
par des huiles végétales.
Mais comme tout mécanisme conçu par les humains, ces moteurs possèdent des
avantages et inconvénients. En plus du fait qu' ils ne sont pas efficaces à haute altitude, là où
la teneur en oxygène est faible, la combustion entraîne le rejet de gaz potentiellement
polluants et voire même toxiques . Ils sont, de ce fait, malgré les efforts engagés dans leur
amélioration, désignés comme une des principales sources de pollution atmosphérique.
8
C’est l’Allemand Rudolf Diesel, qui a conçu en 1893 le premier prototype du moteur
qui porte son nom jusqu’à nos jours. Il s’agissait d’un moteur à 4-temps dont le combustible
est injecté, par le bais d’un système d’injection, directement à forte pression dans la chambre
de combustion. Ce type de moteur, dont le taux de compression est élevé, a connu une
expansion rapide dans le domaine de l’automobile depuis sa création.
Avantages du diesel :
• Le rendement est supérieur à celui d’un moteur à essence: une proportion plus grande
de la chaleur est convertie en travail.
• Le carburant Diesel coûte environ 10 % moins cher que l’essence.
• La consommation moyenne est moins élevée que le moteur essence.
Inconvénients du diesel :
• De fortes contraintes thermiques et mécaniques sont exercées sur les composants
mécaniques du moteur, d’où la nécessité de les surdimensionner.
• Le coût d’entretien est environ 20 % plus élevé que celui d’un moteur à essence.
• L’étanchéité entre le cylindre et le piston est plus difficile a réaliser.
• Le moteur est plus bruyant.
Les moteurs diesel émettent un certain nombre de polluants parmi lesquels :
• Le dioxyde de carbone CO2.
• L’oxyde d’azote NOx.
• Les hydrocarbures HC.
• Les fines particules (dont le diamètre est inférieur à 2.5 micromètres).
• Le dioxyde de soufre SO2.
9
Le moteur à allumage commandé est relativement plus léger et plus vif dans les bas
régimes et lorsqu'il est froid, car l’allumage se fait par l’intermédiaire d’une bougie. Il est
relativement plus silencieux, puisque l’explosion du mélange air-carburant dans un moteur
Diesel, à la suite de l’auto-inflammation, provoque une onde de choc importante dont le bruit
ressemble à une sorte de claquement. Généralement, la mise en marche du moteur à allumage
commandé est plus rapide aux bas régimes, car l'allumage est réalisé par une étincelle
produite par la bougie. Les régimes moteurs étant de surcroît plus élevés, les véhicules de
sport et les voitures de luxe sont généralement munis de moteurs à allumage commandé.
Il est cependant plus facile d'obtenir une puissance plus importante, mais à haut
régime seulement à cause d’un faible couple à bas régime.
Néanmoins, le rendement thermodynamique théorique du moteur à essence est
sensiblement inférieur que celui du moteur Diesel, elle dépasse rarement 30 %. La différence
tient essentiellement aux taux de compression : plus le taux de compression est élevé,
meilleur est le rendement, car le carburant est brûlé plus complètement.
La consommation de carburant est plus élevée dans les moteurs à allumage
commandé, non seulement à cause de la question du rendement déjà évoquée, mais aussi à
cause de la valeur énergétique plus faible de l’essence. À puissance égale, un moteur à
essence aura un régime de rotation supérieur et, en contrepartie, un couple inférieur à bas
régime. Les voitures à essence, même dotées d'un pot catalytique produisent beaucoup plus de
CO2 que les diesels.
santé humaine: le NO2 et les particules fines entraînent des maladies respiratoires et de
l’asthme.
• Les véhicules alimentés au gaz naturel sont moins bruyants que ceux roulants au
carburant Diesel. Cet avantage est un aspect important dans les zones urbaines, ou
actuellement le transport en commun utilise généralement le diesel. De plus, étant donné
qu'un moteur au gaz naturel émet moins de vibrations de résonance qu'un moteur diesel, le
moteur s'use moins et les frais d'entretien sont moins élevés.
• Le gaz naturel étant plus léger que l'air, il se disperse et se dilue rapidement dans
l'air ambiant. Contrairement à l'essence et au gaz du pétrole liquéfié, Il s’enflamme
uniquement quand sa concentration dans l’air se situe entre 5 et 15 %. En dessous de 5 %, le
mélange est trop pauvre pour brûler et au-dessus de 15 %, le mélange est trop riche et ne
brûlera pas non plus. Le point d’inflammation du gaz naturel est d’approximativement
1076 °F, pour le diesel, il est de 500 °F, alors que l’essence s’enflamme à 428 °F. Le gaz
naturel prend donc moins facilement feu que le diesel et l’essence.
Le gaz naturel est transporté en toute sécurité par gazoducs souterrains. Ceci permet
aussi de distribuer du gaz naturel dans les stations-service des zones à forte densité de
population.
Un réservoir de gaz naturel résiste mieux à une collision qu’un réservoir à essence
traditionnel. Le gaz naturel est stocké dans des cylindres dont la paroi mesure 10 à 20 mm
d'épaisseur et qui sont réalisés en matériaux durables et comprennent des dispositifs de
sécurité intégrés.
• Les véhicules roulant au gaz naturel sont actuellement plus chers que ceux
fonctionnant avec les carburants traditionnels. Néanmoins, le gaz naturel véhicule est
significativement moins cher que les carburants classiques tels que l’essence ou le diesel: 1 kg
de gaz naturel compressé coûte à la pompe de 30 à 40 % moins cher qu’un litre de diesel.
L'écart de prix entre entre un véhicule au gaz et un véhicule essence ou diesel va diminuer au
fil du temps.
Grâce à la combustion plus propre du gaz naturel, il y aura moins de frais d’entretien
et la durée de vie du moteur est nettement plus longue que celle d’un moteur essence ou
diesel.
En plus, le gaz naturel possède d’autres avantages tels que:
• Couple disponible à bas régime plus élevé qu’à l’essence.
• Stabilité de fonctionnement accrue.
• Bruits de fonctionnement atténués.
11
• Son utilisation comme carburant ne connait pas les mêmes problèmes techniques que
ceux rencontrés lors de l’utilisation des carburants liquides, tels que: la vaporisation
insuffisante, le bouchon de vapeur (vapor lock), le mélange pauvre, etc.
• Le domaine d’inflammabilité de l’hydrogène dans l’air est très étendu ( de 4 à 75 % en
volume), ce qui rend possible son inflammation avec large gamme de mélanges air-carburant.
Ainsi, le fonctionnement en mélange pauvre est possible, ce qui contribue à une économie de
carburant, assure une réaction de combustion plus complète ainsi qu’une température de
combustion finale plus faible. La quantité de polluants émis lors de la combustion de
l’hydrogène avec l’air tels que les oxydes d’azote est ainsi réduite. Il a une énergie minimale
d’inflammation (énergie d’allumage) très faible, ce qui permet d’enflammer les mélanges
pauvres et assurer un allumage rapide.
Toutefois, l’utilisation de l’hydrogène comme combustible dans le domaine automobile a
des inconvénients :
• L'hydrogène peut s'enflammer ou exploser au contact de l'air. Il doit donc être utilisé
avec précaution, de ce fait, il faut éviter tout risque de fuite, et toute situation confinée peut
s'avérer dangereuse.
• Malgré que l'hydrogène est l'élément le plus abondant dans la nature, il n'en existe
pratiquement pas à l'état pur, c'est-à-dire sous forme de dihydrogène. Ce dernier n'est pas
disponible sous forme brute, mais a besoin d'être synthétisé. La synthèse d'hydrogène
nécessite plus d'énergie que celui-ci peut en produire. La production d'hydrogène nécessite de
l'énergie produite par les centrales thermiques fonctionnant avec du charbon, du pétrole ou
avec du gaz naturel, libérant ainsi du dioxyde de carbone ainsi que d'autres gaz à effet de
serre, ce qui ne ferait que déplacer le problème de la non-utilisation des énergies fossiles et
augmenterait leur consommation par rapport à une utilisation directe de ces énergies dans le
domaine automobile.
• Du fait de sa la légèreté, l'hydrogène, à poids égal, occupe un volume beaucoup plus
important qu'un autre gaz. Concrètement, même comprimés à 700 bars, 4.6 litres
d'hydrogène sont encore nécessaires pour produire autant d'énergie qu'avec 1 litre d'essence.
Ces volumes importants sont une contrainte pour le transport et le stockage sous forme
gazeuse. Pour ces raisons, il sera très difficile de concevoir un réservoir sûr, compact, léger et
moins dispendieux tout en gardant une autonomie relativement élevée.
• En plus d'être facilement explosif dans certaines concentrations à l'air confiné, mais
pas à l'air libre, sous forme de gaz peu comprimé, il occupe trop de place sur un véhicule;
sous forme de gaz très comprimé, il y aura une augmentation du risque d’auto-inflammation
13
qui est une technologie parfaitement maîtrisée et qui dispose d’une importante infrastructure
de distribution et de transport, compatible avec celle de l’hydrogène, permettrait à la
technologie du mélange des deux carburants d’en bénéficier. En plus, le GNC étant déjà plus
propre que les autres carburants conventionnels, il le sera proportionnellement plus si l’on
ajoute certains pourcentages d’hydrogène.
L’injection de l’hydrogène rend la combustion plus complète, ce qui diminue la
consommation et l’encrassement du moteur ainsi qu’une réduction d’émission de polluants
notamment les gaz (NOx, CO, HC) (Tab. 2.1).
2.8 Réferences
[1] NASA Technical Note, May 1977, "Emissions and Total Energy Consumption of a
Multicylinder Piston Engine Running on Gasoline and a Hydrogen-Gasoline Mixture "
(Accessed 2008-08-08).
[2] Idaho National Laboratory on Fuel Enhancement. INL News Release Archive 12/2009.
[4] Mathur H.B., Das L.M., « Performance characteristics of a Hydrogen Fueled SI Engine
using Timed Manifold Injection », dans Int. J. Hydrogen Energy, no vol 16, pp. 115-117,
1991.
[5] Per Tunestal, Magnus Christensen, Patrik Einewall, Tobias Andersson, and Bengt
Johansson, « Hydrogen Addition For Improved Lean Burn Capability of Slow and Fast
Natural Gas Combustion Chambers », dans SAE Technical Paper Series, no 2002-01-2686,
janvier 2002, p. 7–8.
[6] Tsolakis A, Megaritis A, Wyszynski ML, "Application of exhaust gas fuel reforming in
compression ignition engines fueled by diesel and biodiesel fuel mixtures" Energy & Fuels 17
(6): 1464-1473, Nov-Dec 2003.
[7] Yougen Kong, Sam Crane, Palak Patel and Bill Taylor, « NOx Trap Regeneration with an
On-Board Hydrogen Generation Device », dans SAE Technical Paper Series, no 2004-01-
0582, janvier 2004, p. 6–7.
[8] Thorsten Allgeier, Martin Klenk and Tilo Landenfeld, « Advanced Emissions and Fuel
Economy Control Using Combined Injection of Gasoline and Hydrogen in SI-Engines », dans
SAE Technical Paper Series, no 2004-01-1270, janvier 2004, p. 11–12.
[9] « Guidelines for Use of Hydrogen Fuel in Commercial Vehicles Final Report, November
2007 » sur le site de la Federal Motor Carrier Safety Administration, United States
Department of Transportation, novembre 2007.
[10] Hydrogen Fuel Injection System - Technology Fact Sheet for Canadian Hydrogen
Energy Company Ltd, 2005.
[11] Stephen Leahy, « Truckers Choose Hydrogen Power » [archive], dans Wired, novembre
2005.
20
[12] Changwei Ji and Shuofeng Wang, « Effect of hydrogen addition on combustion and
emissions performance of a spark ignition gasoline engine at lean conditions », International
Journal of Hydrogen Energy, vol. 34, septembre 2009.
[13] Marius J. Rauckis, William J. McLean, « The Effect of Hydrogen Addition on Ignition
Delays and Flame Propagation in Spark Ignition Engines », Combustion Science and
Technology, vol. 19, 1979.
[14] Fukutani and Kunioshi, « Fuel Mixing Effects on Propagation of Premixed Flames -
hydrogen plus carbon monoxide flames », Bulletin of the Chemical Society of Japan, 1992.
21
Chapitre 3
Revue de littérature
Les auteurs qui ont prêté une attention particulière au refroidissement accrue, au
"contrôle amélioré de l’huile" par le montage de différents segments, par une augmentation du
balayage des gaz, etc., attribuent le fonctionnement sans retour de flammes qui en résulte, à la
réduction des points chauds, mais ont en même temps (parfois peut-être sans s'en rendre
compte), aux mesures prises pour supprimer la combustion anormale.
Certains auteurs mentionnent qu’une diminution du taux de compression par
l'abaissement de la température de la chambre de combustion pourrait augmenter la résistance
au retour des flammes [9, 18]; d’autres disent que l'augmentation du taux de compression est
conseillée, afin d’augmenter la surface de la chambre de combustion par rapport au volume,
ce qui améliore le transfert de chaleur et le refroidissement des gaz résiduaires [5, 19]. Un
taux de compression élevé permet également de réduire la quantité de résidus. Les deux
suggestions sont valides et indiquent l'existence d'un taux de compression optimal: son
augmentation aura comme conséquence l’augmentation de la puissance grâce à un rendement
effectif élevé jusqu'à un certain point, où le mélange doit être allégé afin d'éviter le pré-
allumage et la chute de la puissance [6, 20].
Le phénomène du cliquetis (ou encore cognements) dans les moteurs à hydrogène, qui
est une combustion anormale entraînant une résonance de l'explosion sur les parois de
la chambre de combustion et du piston, a été mal étudié. Pour ce est qui du retour des
flammes, certaines causes ont été mélangées ou leurs effets surestimés, mais en ce qui
concerne le cliquetis, il y avait des revendications totalement contradictoires dans la
littérature.
Pour commencer, la plupart des articles ne mentionnent pas que la résistance au
cliquetis est une propriété du mélange carburant / air, c’est à dire, en précisant les indices
d'octane sans définir la richesse du mélange correspondante. Certains prétendent que l'indice
d'octane est très faible [21, 22], d'autres affirment qu'il est très élevé [14, 5].
Un auteur d’un article est même allé dire qu’à la fois «l'hydrogène a un indice d'octane
élevé effectif» et «l'indice d'octane effectif d'hydrogène est plutôt faible» [23]. Seuls de très
rares documents mentionnent les indices d'octane en fonction de la richesse [19, 24]. Des tests
ont rapporté que l'hydrogène peut agir comme un agent anti-cliquetis lorsqu'il est ajouté à un
isooctane sans plomb [5].
Il existe certaines preuves que les causes du cliquetis des moteurs à hydrogène
pourraient être différentes de celles des moteurs à essence, ou le cliquetis est dû
essentiellement aux vitesses de flamme excessives [5]. Ainsi, la réduction de la vitesse de
24
montée en pression peut être plus efficace pour contrôler les cliquetis que la limitation de la
durée de combustion [8].
Il est à noter que le travail expérimental et théorique de Karim et al. [22, 26] fait état
de très grandes zones de cliquetis, où des mélanges stoechiométriques ont tendance à
détonner, même à des taux de compression aussi faibles que 6:1. Comme ces résultats sont en
désaccord avec toutes les autres expériences rapportées dans la littérature, il semble qu’ils
sont probablement affectés par des causes inconnues aux auteurs.
L’examen de la littérature liée aux études expérimentales sur les moteurs à allumage à
hydrogène montre que le pré-allumage pour certains auteurs est un facteur limitant pour les
taux de compression, le moment d’allumage et les richesses du mélange, plutôt que le
cliquetis. Les mesures réalisées avec un taux de compression de 11:1 et une pression de
suralimentation de 0.85 bar (manométrique) sur des mélanges stoechiométriques ont été
rapportés [6], ainsi que des mesures sur des mélanges pauvres en utilisant des taux de
compression de 14:1 et plus [14, 25], le tout sans aucune apparence de cliquetis.
Plusieurs méthodes de formation du mélange ont été testées dans les moteurs à
hydrogène, afin d’éviter le phénomène du retour des flammes:
• formation externe du mélange avec un carburateur à gaz [4, 19],
• formation externe du mélange avec induction parallèle, c'est à dire: un moyen de
retarder l'introduction de l'hydrogène, par exemple une conduite de carburant fermée par une
vanne séparée sur le dessus de la soupape d'admission qui s'ouvre uniquement lorsque la
soupape d'admission est assez levée [27],
• formation externe du mélange de gaz avec un carburateur et injection d'eau [9, 21],
parfois avec recirculation des gaz d'échappement supplémentaire (RGE) [28],
• formation externe du mélange avec collecteur chronométré ou injection du carburant
dans l’orifice d’admission (ICOA) [14, 2,6, 11, 16, 20, 29], parfois aussi avec des moyens
d'induction parallèle [24],
• la formation du mélange interne par injection directe (ID) [30, 31, 32, 33].
L'injection d'eau, l'introduction tardive de l'hydrogène et l'injection directe sont toutes
principalement destinées à retarder ou à empêcher le retour des flammes soit par un
refroidissement supplémentaire ou en évitant l’inflammation du mélange pendant la phase
d'admission. Au cours de la dernière décennie, seulement l’injection du carburant dans
l’orifice d’admission et l’injection directe (lors de la course de compression ou plus tard) ont
25
été utilisées, comme les autres méthodes sont moins flexibles et incontrôlables. Il a été
démontré que la formation externe du mélange à l'aide de l’injection du carburant dans
l’orifice d’admission mène à l’augmentation des rendements du moteur, à la prolongation du
fonctionnement avec mélange pauvre et à la réduction des émissions des NOx par rapport à
l’injection directe [25, 34]. Ceci est la conséquence de l'homogénéité élevée du mélange en
raison de longues durées de formation du mélange dans les moteurs à injection dans l’orifice
d’admission (MIOA) ainsi que les durées de formation du mélange dans les moteurs à
injection directe (MID), étant donné que les turbulences générées de l’admission contribuent
moins à la formation du mélange. En outre, le coût et la complexité sont significativement
plus faibles pour les MIOA que pour les MID [7] et l'adaptation d'un moteur existant est
possible. D'autre part, la puissance d'un moteur à hydrogène avec formation du mélange
externe est limitée en raison de la diminution du rendement volumétrique: en raison de la
faible densité de l'hydrogène et l'exigence de petites quantités d’air des mélanges
stoechiométriques, le volume du cylindre occupé par de l'hydrogène dans un mélange
stoechiométrique atteint les 29.5%. Il en résulte une diminution de la teneur en énergie
volumique d'environ 18% pour l'hydrogène par rapport à l'essence. Si l'injection directe est
utilisée pour introduire l'hydrogène après que la soupape d'admission soit fermée, la puissance
maximale peut être supérieure de 17% par rapport à l'essence.
Un avantage important de l’injection directe par rapport à l’injection dans l’orifice
d’admission est l'impossibilité du retour de flamme. Cela augmente aussi la puissance
maximale de l’injection directe par rapport à l’injection dans l’orifice d’admission comme les
plus riches mélanges peuvent être utilisés sans crainte de retour de flamme. Le pré-allumage
peut encore se produire sauf si l'injection très tardive est utilisée. La formation externe du
mélange offre une plus grande liberté en ce qui concerne les méthodes de stockage: l’injection
directe pendant la course de compression a besoin d'hydrogène à haute pression et nécessite
donc le stockage de l'hydrogène liquide.
Ainsi, la formation du mélange externe et celle interne ont leurs avantages et leurs
inconvénients. L’injection directe est préférable pour des performances à pleine charge
(puissance maximale), l’injection dans l’orifice d’admission est préférable à charge partielle
(rendement maximal du moteur). Des modèles de moteurs ont été proposés en utilisant les
deux techniques de formation de mélanges [34, 35, 36].
Si un moteur à hydrogène est conçu pour fonctionner avec une seule vitesse de
rotation /puissance, par exemple, pour la production d'électricité ou pour un véhicule hybride,
le fonctionnement avec minimum d’émissions et avec un rendement très élevé, est possible
sans post-traitement des gaz d’échappement (dont le rendement pourrait se détériorer au fil du
temps).
Les émissions des NOx inférieures à 10 ppm, voire 1 ppm, avec des rendements
indiqués de près de 50% sont possibles [25, 38, 39]. L'hydrogène est le seul carburant avec
lequel cela est possible (avec les hydrocarbures, la diminution des émissions NOx à mélange
pauvre entraîne une augmentation des émissions d'hydrocarbures imbrûlés).
Dans cette section, une tentative est faite pour fournir un aperçu complet des
caractéristiques de conception des moteurs à hydrogène du point de vue de ses avantages et de
ses inconvénients.
• Les bougies d'allumage: on recommande d’utiliser une bougie d’allumage froide pour
éviter que les températures des électrodes de la bougie dépasseent la limite d'auto-
inflammation et provoquent ainsi le retour de flamme [3, 10]. Ne pas utiliser des bougies avec
des électrodes en platine, car cela peut être un catalyseur pour oxydation de l'hydrogène [5, 7]
(le platine a été utilisé dans l'échappement pour oxyder l'hydrogène non brûlé [27]).
• Le système d'allumage: on peut éviter l'inflammation incontrôlée due à l'énergie
d'inflammation résiduelle par une bonne mise à terre du système d'allumage ou la
modification de la résistance électrique du câble d'allumage [9, 10]; éviter l'induction
d'allumage dans le câble d'allumage à proximité, par exemple à l'aide d’un système de bobine
sur la bougie; fournir un système d'allumage à haute tension puisque l'allumage des mélanges
d'hydrogène demande une tension d'allumage secondaire accrue [9, 10, 40] (probablement en
raison de la plus faible concentration d'ions de la flamme d'hydrogène), le système de bobine
sur bougie peut aussi satisfaire à cette condition. Par ailleurs, l'écartement des électrodes peut
être diminué pour abaisser la tension d'allumage, ce n'est pas un problème pour les moteurs à
hydrogène, car il n'y aura presque pas de formation de dépôts. Des électrodes avec des
écartements aussi petits que 0.25 mm ont été utilisées [4] (bien que l'écart ait ensuite été porté
à 0.5 mm en raison des difficultés de démarrage à froid dues à la condensation de l'eau à la
pointe de la bougie).
• Le système d'injection: on peut fournir une injection synchronisée, soit en utilisant
l'injection dans l’orifice d’admission et la programmation du système d’injection de telle sorte
28
qu'une période initiale de refroidissement de l'air est crée dans la phase initiale de la course
d'admission et la fin de l'injection est telle que tout l'hydrogène est admis, sans laisser des
traces d'hydrogène dans le collecteur d'admission lorsque la soupape se ferme, ou par
injection directe lors de la course de compression. Des injecteurs à haut débit sont nécessaires
dans les deux cas, plusieurs injecteurs par cylindre peuvent atténuer cette exigence. La
synchronisation décrite ici peut ne pas être nécessaire puisque les travaux ont été rapportés
où aucune relation entre le calage de l'injection et le retour des flammes ou le rapport
d’équivalence limité du pré-allumage n’a été trouvée [14]. L’injection synchronisée diminue
aussi la quantité de carburant non brûlé dans le collecteur d'admission à un moment donné, en
limitant la gravité du retour des flammes s’il se produit.
• Les points chauds: on peut éviter les points chauds dans la chambre de combustion qui
pourraient provoquer le pré-allumage ou le retour des flammes, en utilisant les soupapes
d'échappement refroidies; utiliser les culasses à soupapes multiples pour abaisser encore la
température de la soupape d'échappement [7, 8, 9]; assurer un contrôle adéquat de l'huile;
fournir d'autres passages du liquide de refroidissement autour des soupapes et autres zones
présentant de fortes contraintes thermiques [11] (si possible), retarder l’injection du carburant
afin de créer une période de refroidissement de l'air (à l'aide du collecteur synchronisé ou
l’injection directe); assurer le piégeage approprié (par exemple en utilisant des calages
variables [14, 6]) afin de réduire la température des gaz résiduels.
• Les segments de piston et les volumes de crevasse: on peut diminuer l’espacement du
cordon du piston pour éviter que les flammes d'hydrogène se propagent dans le cordon, Swain
et al. [11], utilisent un espacement de 0.152 mm pour éteindre la flamme d'hydrogène;
modifier les volumes de crevasse et/ou les segments du piston dans le but de réduire la
redistribution de mélange non brûlé du second cordon vers le cordon supérieur [11, 12, 13]
(prévention «ravitaillement» de la flamme du cordon supérieur pendant l’échappement et
l’admission).
• Les sièges de soupapes et les injecteurs: le pouvoir lubrifiant très faible de l'hydrogène
doit être pris en compte, les matériaux de siège de soupape appropriés doivent être choisis [7,
9] et la conception de l'injecteur doit en tenir compte. C'est le cas avec n'importe quel
combustible gazeux sec (comme le gaz naturel), mais peut être plus critique pour l'hydrogène
(le gaz naturel comprimé contient de petites quantités de pétrole provenant de la vapeur
d'huile dans le compresseur alors que les compresseurs d'hydrogène ont normalement plus
d’espacements plus étroits pour limiter le taux de fuite).
29
peut être nécessaire pour limiter les émissions de NOx. Ceci ne peut être réalisé qu’avec un
système d'entraînement à commande électronique.
La littérature sur la simulation des moteurs à hydrogène est assez rare. Fagelson et al.
[43] utilisent un modèle à deux zones quasi dimensionnel pour calculer la puissance et les
émissions de NOx d'un moteur à combustion interne à hydrogène. Ils utilisent un modèle
semi-empirique de la combustion turbulente de la forme ut = A· B· Re· ul, où A et B sont des
constantes, Re est le nombre de Reynolds basé sur le diamètre du piston, la vitesse moyenne
du piston et les propriétés des gaz brûlés; ut et ul sont respectivement les vitesses de
combustion turbulente et laminaire. La propagation de la flamme sphérique est prise en
charge, le transfert de chaleur est négligé et la formation de NOx est calculée en utilisant 10
zones à masse constante dans le mécanisme des gaz brûlés et le mécanisme élargi de
Zeldovitch. La vitesse de combustion stratifiée est calculée à partir d'une réaction globale du
second ordre avec une énergie d'activation estimée.
Le modèle est validé par rapport aux mesures en faisant varier le rapport d'équivalence
et le calage de l'allumage seulement. Prabhu-Kumar et al. [44] utilisent ce modèle pour
prédire la performance d'un moteur à hydrogène suralimenté, mais aucune modification n'est
apportée au modèle d'origine. Ils rapportent une surestimation de la vitesse d’augmentation de
pression (et donc de la vitesse de combustion).
Keck [45] a réalisé des mesures optiques sur un moteur, fonctionnant au propane ainsi
qu’à l'hydrogène. Il a utilisé un modèle d'entraînement turbulent pour comparer les tendances
prédites avec les tendances observées expérimentalement.
Johnson [46] a utilisé le code de simulation du moteur Kiva-3V développé au
Laboratoire Nationale de Los Alamos selon la norme du modèle standard «eddy-turnover»
pour simuler un moteur à hydrogène avec rapport équivalence et efficacité volumétrique fixes.
Le modèle standard contient un paramètre libre qui est adapté pour l'hydrogène maintenu
constant pour faire varier la vitesse et la séquence d'allumage du moteur. Le modèle est validé
par rapport aux tests rapportés dans la référence [38].
Fontana G. et al. [47] ont modifié le code Kiva-3V, afin de simuler un moteur alimenté
par un mélange hydrogène/essence. Ils ont utilisé un modèle hybride où le taux de réaction
globale est donné soit par le modèle standard «eddy-turnover » ou par un taux de réaction
pondéré basé sur deux expressions de taux de réaction global, l'un pour la combustion
31
d'hydrogène et l'autre pour l'essence. Ils ont validé le modèle de fonctionnement à l'essence et
calculé les effets de l'addition de diverses concentrations d'hydrogène à l'essence.
Les modèles zéro-dimensionnels et multidimentionnels ont été utilisés pour la
simulation d'un moteur à hydrogène à l'Université Technique Tchèque [48, 49]. Un modèle
zéro-dimensionnel sur la base du code «GT-Power» est utilisé avec la loi de Wiebe monté sur
les taux mesurés de dégagement de la chaleur. Le mécanisme élargi de Zeldovitch est utilisé
pour calculer les émissions de NOx. Le «modèle Eulérien avancé multizone» a été développé
pour la simulation multi-dimensionnelle. Ce modèle est un mélange d'éléments zéro-
dimensionnels et multi-dimensionnels: la chambre de combustion semble être limitée à des
géométries simples en raison des limitations de la génération de la grille, et le transfert de
chaleur est modélisé pour le contenu du cylindre comme un volume en vrac. Le modèle de
combustion est un modèle semi-empirique « PDF-like», qui s'appuie sur un taux mesuré de
dégagement de chaleur et l'hypothèse d'un front de flamme hémisphérique pour suivre la
propagation des flammes. Une deuxième option pour le modèle de combustion est basée sur
un mécanisme de réaction détaillée, mais n'est apparemment pas intégrée dans le code
complet du moteur encore en raison de «difficultés numériques». Les auteurs utilisent le
premier modèle de combustion pour les optimisations, car il permet des temps de calcul
limités, mais leur but final est l'intégration de la cinétique chimique dans le code.
Enfin, Ma et al. [50] utilisent un modèle zéro-dimentionnel en utilisant la loi de
Wiebe. Le modèle est utilisé pour calculer les effets de la variation du taux de compression et
du calage de l'allumage et de déterminer un diamètre optimal du cylindre, pour un rapport
d'équivalence fixe. Aucune validation des données expérimentales n’a été réalisée et toute
autre justification de l'extrapolation en dehors des conditions pour lesquelles l'ajustement est
valable n’est donnée.
Les études sur l’utilisation du gaz naturel comme carburant de substitution se sont
développées ces dernières années. En effet le gaz naturel est un mélange gazeux qui possède
un fort potentiel de dépollution, mais aussi dans l’aspect performance quant à son utilisation
comme carburant. Il a été nécessaire d’en étudier les aspects les plus prometteurs dus à ses
propriétés gazeuses.
Goto [51] a analysé le processus de la formation du mélange à proximité de la bougie
d’allumage dans un moteur à injection directe, alors que Rubas et al [52] ont étudié dans un
moteur au gaz naturel à injection directe, la pénétration du gaz dans la chambre de
32
combustion. Hiltner et al [53] ont montré que l’injection séquentielle dans un moteur
fonctionnant au gaz naturel a une influence sur la distribution du carburant dans la chambre de
combustion. Aesoy et al [54] et Naber et al [55] ont cherché à connaître l’influence de la
composition du gaz naturel sur le délai d’allumage.
Quatre mélanges énumérés dans le tableau 3.1 ont été choisis pour représenter la variabilité
observée dans la composition du gaz naturel.
Tableau 3.1 Mélanges de gaz naturel sélectionnés pour la combustion à injection directe
(Naber et al [55]).
Constituant Méthane Moyenne Superethanol Écrêtement des
du mélange (% en vol.) (% en vol.) (% en vol.) pointes (% en vol.)
N2 0 1.37 7.76 4.34
CO2 0 1.21 1.43 1.13
CH4 100.00 93.12 80.56 87.15
C2H6 0 3.20 8.99 3.00
C3H8 0 0.70 1.00 1.00
C4H10 0 0.40 0.26 0.26
Mole, poids 16.04 17.36 19.02 18.58
Méthane N° 100 91 75 94
L'appareil expérimental utilisé est constitué de deux composants principaux: une cuve
de combustion à volume constant et un système d'injection de carburant. La température
minimale nécessaire au retard d’allumage raisonnable du moteur diesel (<2 ms) pour tous les
rapports examinés du gaz naturel ajouté est d'environ 1150K, Il correspond étroitement aux
conditions au PMH obtenues dans un moteur diesel à injection directe à gaz naturel à
allumage par compression. Aux hautes températures, les délais d’allumage mesurés sont
limités par le taux d'injection de carburant. L'augmentation élevée de la concentration des
hydrocarbures dans le gaz naturel diminue le délai à l'allumage d'environ 0.8 ms
comparativement à celui du méthane pur sous les conditions de moteurs diesel. Cependant,
d'autres augmentations de concentrations de propane et de butane pourraient conduire à leur
condensation dans le système d'injection, ce qui entraîne des problèmes de contrôle de
carburant. Une bonne convergence entre les tendances dans les effets mesurés
expérimentalement de la composition du gaz naturel sur le délai d'allumage et ceux prédits par
un modèle cinétique, indique que les variations du délai d'allumage produites par des
différences de composition sont principalement dues à des effets cinétiques. Enfin, l'indice de
33
méthane ne corrèle pas bien avec le délai d'allumage mesuré et prédit dans des conditions du
moteur diesel.
Patrik Einewall et al. [56] à LIT en Suède ont effectué des tests où ils ont comparé le
fonctionnement moteur avec mélange pauvre avec le fonctionnement aux conditions
stoechiométriques avec RGE et un catalyseur à trois voies pour un moteur alimenté au gaz
naturel. Le moteur (TG103/G10A) a été développé à l'origine pour fonctionner avec du
carburant diesel et remanié par Volvo pour le fonctionnement avec du gaz naturel (voir le
tableau 3.2). Dans ces tests, le carburant est injecté dans l’orifice d’admission (voir figure 3.1)
et le moteur est équipé d'un système RGE refroidi (voir figure 3.2).
Ils ont constaté que le fonctionnement aux conditions stoechiométriques avec RGE et
catalyseur (voir figure 3.3 et figure 3.4) a donné beaucoup moins de NOx et des HC
(respectivement 10-30 et 360-700 fois plus faible) avec seulement une légère diminution du
rendement par rapport au fonctionnement avec mélange pauvre. Les émissions du CO,
toutefois, ont été trouvées à environ 10 fois plus élevées que pour la solution RGE, mais il
convient de noter que le contrôle de λ n'a pas été optimisé pour l'utilisation d'un catalyseur à
trois voies.
Fig. 3.3 Les émissions spécifiques des HC, NOx et du CO avant le catalyseur,
30% RGE (PME = 14 bars) (Patrik Einewall et al. [56]).
35
Fig. 3.4 Les émissions spécifiques des HC, de NOx et du CO après le catalyseur, 30% RGE
(PME = 14 bars) (Patrik Einewall et al. [56]).
Nellen et al. [57] ont développé un concept de moteur à gaz naturel pour les
applications de cogénération stationnaire. Le moteur fonctionne aux conditions
stoechiométriques et muni d’un système de RGE refroidi et d'un catalyseur à trois voies. Ils
ont réalisé une pression moyenne effective maximale (PMEM) égale à 23 bars avec la charge
thermique du moteur et le niveau du stress mécanique similaires à celles du même moteur
diesel à injection directe. Sur toute la plage des charges du moteur étudié (PMEM = 12 à
23 bars), les niveaux des émissions des NOX sont inférieurs à 2 mg/Nm3 avec 5% de O2 et de
CO, les niveaux des émissions des hydrocarbures sont inférieurs à 75 mg/Nm3(mg de
contaminant par Nm3 de gaz sec).
Dans le tableau 3.3, on résume les données des performances du moteur obtenues, en
les comparant avec celles du même moteur optimisé fonctionnant avec mélange pauvre, sans
et avec catalyseur d’oxydation. L'amélioration de toute évidence est frappante.
36
Tableau 3.3 Les résultats obtenus avec le nouveau concept du moteur contre ceux du moteur à
gaz prototype (Nellen et al. [57]).
Rendement
NOX CO HC PMEM
du moteur Mode de
mg/Nm3 g/kwh mg/Nm3 g/kwh mg/Nm3 g/kwh bar % Fonctionnement
(5% O2) (5% O2) (5% O2)
RGE, λ=1 et
~
~1 ≤ 75 ≤ 0.21 ≤ 30 ≤ 0.08 12/20 ~ 40/42 catalyseur
0.003
3-voies
RGE, λ>1 et sans
< 250 < 0.79 < 650 < 2.05 ~1000 ~3.15 12 ~36.5
catalyseur
λ>1 avec
250 0.76 100 0.30 1500 4.54 12 38 catalyseur
d’oxydation
Caractéristiques du Moteur :
Les tests selon la procédure d'essai fédérale Américaine, ont révélé les émissions
suivantes : 0.049 g / bhp-h (g/kWh) · 0,7457 = g/bhp-h) des NOx, 4.153 g/bhp-h de CO, de
0.002 g/bhp-h de particules solides, 0.435 g/bhp-h d’hydrocarbures et 0.000 g/bhp-h
d’hydrocarbures non-méthaniques. La consommation spécifique du carburant au frein était de
2% au-dessus du moteur à mélange pauvre calibré à 2 g/bhp-h de NOx.
Zuohua Huang et al. [61] ont étudié les caractéristiques de la combustion à différents
moments d’injection d’un moteur à injection directe fonctionnant avec un mélange de gaz
naturel et d’hydrogène avec la durée d’injection et le moment d’allumage fixées. Cette étude a
montré que:
• L’injection tardive réduit le taux d’excès d’air et rend le mélange plus pauvre.
• La pression efficace augmente quand on avance le moment d’injection du carburant.
Elle atteint la valeur maximale au moment d’injection 190 °CA avant PMH.
• La durée de combustion diminue quand on avance le moment d’injection du
combustible.
• Les émissions des NOx et du CO2 augmentent quand on avance les moments
d’injection du carburant. Toujours selon l’article l’ajout d’hydrogène dans le gaz naturel peut
réduire la concentration du CO2.
I. Saanum et al. [62] ont réalisé des tests sur un moteur Volvo TD100 à allumage
commandé, du 9.6 litres, fonctionnant au gaz naturel et avec un mélange d’hydrogène/gaz
naturel en proportion 25/75 (en volume). On a testé deux stratégies pour de faibles émissions
de gaz toxiques: fonctionnement avec mélange pauvre et fonctionnement selon les conditions
stoechiométriques avec RGE et avec un catalyseur à trois voies. Les tests ont été réalisés dans
le laboratoire de moteur à combustion de l'Université de Lund en Suède. Le banc d'essai est le
même que celui utilisé en [56].
39
Avec la régulation lambda (la sonde de régulation située avant le catalyseur mesure
l'oxygène se trouvant dans les gaz d'échappement. Elle envoie un signal au calculateur
("ECU" : "Engine Control Unit") pour modifier la composition du mélange via le système
d’injection), les émissions du CO dans le cas de la RGE et du catalyseur sont les mêmes ou
inférieurs aux émissions des CO lors du fonctionnement en mélange pauvre sans catalyseur à
la fois pour le gaz naturel et le mélange avec 25% d'hydrogène. Les émissions du CO et des
HC pour le cas de la combustion avec mélange pauvre sont faibles pour le mélange avec 25%
d’hydrogène pour la même λ. La possibilité de réduire les NOx sans augmenter les émissions
des HC par l'ajout d'hydrogène pour prolonger la limite du mélange pauvre est confirmée. Le
rendement thermique maximal est plus élevé pour le fonctionnement aux mélanges pauvres
40
que pour le fonctionnement avec la RGE, tant pour le gaz naturel, que pour l’hydrogène
enrichi avec du gaz naturel.
P. Tunestål et al. [63] ont effectué des tests sur un moteur alimenté avec un mélange
de gaz naturel et d’hydrogène. Ils ont testé deux chambres de combustion avec différents
niveaux de turbulence. Ils ont constaté que l'ajout d'hydrogène conduit à une réduction des
NOx et des HC par rapport au gaz naturel pur. L'effet de l'addition d'hydrogène a été plus
prononcé pour la turbulence lente de la chambre de combustion que pour la chambre à haute
turbulence.
Des tests approfondis des mélanges hydrogène / gaz naturel comme combustible ont
été réalisés par Munshi et al. [64]. Les tests sur dynamomètre ont révélé une réduction de 50%
des émissions des NOx pour le fonctionnement stationnaire et la réduction de 56% dans un
cycle transitoire lorsque le moteur est alimenté avec un mélange de 20 % (en vol.) d’H2 par
rapport au gaz naturel. L'Agence Transit SunLine a testé deux bus avec ces moteurs
fonctionnant avec mélange pauvre.
Raman et al. [65] ont reporté que l’hydrogène est un additif profitable pour le gaz
naturel. Il donne la possibilité au moteur de fonctionner avec un mélange pauvre sous des
conditions de charge partielle et il améliore la pression moyenne efficace pour un papillon
complètement ouvert aux alentours de la limite inférieure de la combustion. De plus, les
auteurs ont observé une réduction de NOx avec 15% à 20% d’hydrogène dans le mélange,
avec une certaine augmentation des HC à la suite du mélange carburant/air ultra-pauvre.
Choudhuri et al. [66] ont réalisé une étude expérimentale des flammes de diffusion du
mélange gaz naturel/hydrogène. L’écoulement des gaz frais est laminaire. Les radicaux libres
sont mesurés par Fluorescence induite par Laser (FIL). Les profils radiaux de radicaux et la
température sont mesurés à 8 mm au-dessus du brûleur pour des mélanges gaz naturel/(0-
100%)hydrogène/air. Les fractions molaires maximales OH, CH, O et H augmentent avec la
quantité d’hydrogène et cette évolution converge avec les prédictions du mécanisme de Peters
[67].
Dagaut et Nicolle [68] et Dagaut et Dayma [69] ont examiné la combustion de
mélanges gaz naturel/(0, 40, 75%)H2/air en réacteur ouvert parfaitement agité. L’étude [68]
est réalisée à pression atmosphérique, pour des températures comprises entre 900 et 1450 K et
pour un temps de séjour fixé à 120 ms. Des mélanges pauvres (φ = 0.3) et stoechiométriques
de gaz naturel (89%CH4/9%C2H6 /2%C3H8) sont étudiés. Dagaut et al. [69] ont fixé la
pression à 10 atm, le domaine de température à 900-1200 K et le temps de séjour à 250 ms.
Le gaz naturel est un mélange (90%CH4/10%C2H6) et son oxydation est étudiée aux richesses
41
0.3, 0.6 et 1. Les profils d’évolution des espèces stables sont établis en fonction de la
température du réacteur. Comme on le voit bien sur Figure 3.7, les résultats expérimentaux
convergent avec les prédictions du mécanisme de Dagaut décrit dans [70].
La substitution partielle de l’hydrogène au méthane à pression atmosphérique aussi
bien qu’à 10 atm, conduit à une plus grande propension du mélange à s’oxyder,
particulièrement en mélange pauvre. La température à laquelle la combustion du combustible
produit des intermédiaires (CO, éthylène) au même niveau de concentration est réduite de 50
à 100 K en présence d’hydrogène.
Fig. 3.7 Oxydation de mélanges GN/40%H2 dans un réacteur parfaitement agité aux richesses
0.3 et 1 à une pression de 1 à 10 atm (Dagaut et Nicolle [68] et Dagaut et Dayma [69]).
Huang et al. [71] ont mesuré les vitesses de flamme laminaire gaz naturel/hydrogène/
/air dans une enceinte sphérique (fig. 3.8) à 298 K et à pression atmosphérique.
ux −u0 x
= 0.00737 exp �20.38� + 0.00334 (3.1)
u100 − u0
En 1993, un test a été réalisé par Swain [72] sur un moteur Toyota, 1.6 l, quatre
cylindres, afin de déterminer l’influence de l’utilisation de l’Hythane sur les performances du
moteur sus-mentionné. Le test a été exécuté avec une vitesse de rotation constante égale à
1000 tr/min, une avance à l’allumage (AVO) optimal et une charge minimale. L’Hythane
comme carburant (pour la même richesse), si on le compare au méthane, augmente le
rendement thermique et les émissions des NOx, cependant il diminue l’AVO, les
hydrocarbures imbrulés et les oxydes de carbone. En combinaison avec d’autres carburants,
43
Le gaz naturel comprimé enrichi par hydrogène (GNCH) est une alternative potentielle
aux carburants issus des hydrocarbures dans les moteurs à combustion interne. Des recherches
expérimentales sur l’utilisation de ce genre de carburant ont été conduites pendant plus de dix
années et sont toujours en cours.
M. Bysveen [73] a évalué le rendement et les émissions d'un moteur fonctionnant
respectivement au gaz naturel comprimé (GNC) et avec un mélange de gaz naturel et
d'hydrogène. Les mélanges de gaz naturel comprimé et d'hydrogène sont appelés mélanges
GNCH. Deux différents mélanges ont été étudiés; le GNC pur, et 29% H2 (en vol.) dans GNC
(nommé GNCH). Pour chaque mélange, quatre régimes ont été testés dans des conditions de
pleine charge. Pour chaque mélange et pour chaque régime moteur, on a choisi quatre
différents coefficients d’excès d'air (λ). La limite de mélange pauvre du GNC pur testé dans
ce moteur est d'environ λ = 1.8 et la limite du GNCH est encore plus basse, c'est à dire
d’environ λ = 2.0. Ceci est accepté et conforme à la littérature. Compte tenu des émissions
totales des hydrocarbures imbrûlés et des oxydes d'azote, il est facile de constater que l'effet
positif de l'ajout d'hydrogène au GNC augmente en fonction de lambda. En ajoutant de
l'hydrogène au GNC, une augmentation de la puissance permet d’atteindre la limite du
mélange pauvre étudiée. Le rendement du GNCH est supérieur à celui du GNC pour le même
λ et la différence de rendement théorique ηth, frein entre GNCH et GNC pour le même λ,
augmentera avec l’augmentation du rapport d’excès d’air.
Nagalingam et al. [74] ont mené des tests sur un moteur AVL monocylindrique à
allumage par étincelles, quatre temps, muni de quatre soupapes, d’un alésage de 125 mm,
d’une course de 130 mm et d’un taux de compression de 11.73, fonctionnant avec un mélange
de GNC/H2 en proportion (en %) de: 100/0, 80/20, 50/50, 0/100. Ils ont conclu que la perte
maximale de puissance pour le moteur fonctionnant à l’hydrogène pur est de l’ordre de 23%
et la réduction maximale du rendement thermique indiqué atteint les 12% pour une vitesse de
rotation de 1200 t/min. L’augmentation des dimensions des moteurs stationnaires
fonctionnant avec de l’hydrogène va compenser les pertes en puissance. La période
d’allumage optimale diminue au-dessus des 20° avant le point mort bas (PMB) pour le
mélange avec 100% d’hydrogène. Les émissions de NOx augmentent lors de l’utilisation de
44
1.3 (φ devrait être supérieure à 0.77), alors que pour avoir des faibles émissions des NOx, λ
devrait être égale au moins à 1.5. Il est constaté que la teneur en hydrogène égale à 80%,
garantit un fonctionnement moteur sans anomalies (sans retour des flammes et sans
détonation), quel que soit le facteur d'excès d'air. Il est démontré aussi que, pour obtenir un
rendement maximal du moteur pour toute la gamme de charges tout en prenant en
considération les basses émissions, la composition du mélange doit être modifiée en fonction
de la charge du moteur.
Das et al. [78] ont étudié le rendement thermique et la consommation spécifique du
carburant d’un moteur à combustion interne fonctionnant au gaz naturel comprimé et avec de
l’hydrogène. Ils ont constaté que la consommation est réduite et le rendement thermique est
amélioré avec l’hydrogène pur comme carburant par rapport au moteur fonctionnant avec du
GNC. Le rendement thermique était de l’ordre de 31.19% pour l’hydrogène et de 27.59%
pour le GNC.
Ma F. et al [79] ont réalisé une étude expérimentale se concentrant sur les effets
d'ajout d'hydrogène sur les performances et les émissions du moteur. Les données
expérimentales ont été compilées sur moteur suralimenté à allumage par étincelle
fonctionnant avec un mélange pauvre de gaz naturel dans diverses conditions, y compris le
ralenti. Des données ont été prises lors du fonctionnement avec du gaz naturel pur aussi bien
qu'avec du gaz naturel comprimé (GNCH) enrichi avec de l’hydrogène, y compris les
fractions volumétriques d'hydrogène de 20, 30 et 40%. Pour la comparaison, la capacité
calorifique inférieure pour chaque taux hydrogène/ gaz naturel est demeurée identique à
chaque condition de fonctionnement du moteur. Les résultats prouvent que, dans des
conditions de ralenti, l'ajout d'hydrogène peut réduire le rendement du moteur ainsi que les
émissions d'oxyde de carbone et des hydrocarbures. Ce sont les deux principaux grands
problèmes liés au régime de ralenti. On constate également que cette réduction est plus
évidente aux moments d’allumage relativement retardés. En conditions de fonctionnement
normal, l'addition d'hydrogène s'est avérée bénéfique simultanément pour l'économie du
combustible et pour les émissions des HC. Bien que l'addition d'hydrogène augmente les
émissions de NOx, ceci peut être compensé en retardant le moment d’allumage et en faisant
fonctionner le moteur avec un mélange pauvre, qui est réalisable par la vitesse de combustion
plus rapide et la limite plus large d'inflammation d'hydrogène.
Ma F. et al [80] ont réalisé un travail expérimental visant l’étude des effets d'ajout
d'hydrogène sur les comportements de la combustion et les variations cycle par cycle dans un
moteur 6 cylindres en ligne, 6.2 litres, suralimenté alimenté par gaz naturel et fonctionnant
46
avec un mélange pauvre. Sur la figure 3.9, on présente le schéma du système d’alimentation
en carburant. Les tests ont été effectués avec différents : fractions d'hydrogène, richesses,
moments d’allumage, pressions absolues dans la tubulure d’admission et vitesses de rotation
du moteur.
On a constaté que l'ajout d'hydrogène contribue beaucoup plus à réduire la durée de
formation de la flamme, dont la réduction exerce considérablement des effets positifs sur la
réduction des variations cycle par cycle que sur la réduction de la durée de propagation de
flamme. Basés sur ces résultats de combustion, les effets d’ajout d’hydrogène sur les
variations cycle par cycle pour différents : moments d’allumage, richesses, pressions absolues
à la tubulure d’admission et vitesses de rotation du moteur, ont été également examinés et on
a conclu que les variations cycle par cycle, qui sont d'une importance élevée dans des moteurs
à combustion interne fonctionnant avec mélange pauvre, peuvent en effet être abaissées par
l'ajout d'hydrogène. L'efficacité thermique du moteur ainsi que les émissions des gaz
d’échappement ont été également examinées lors des tests. On a observé qu'après l’ajout
d'hydrogène, le rendement thermique du moteur pourrait être gardé relativement élevé sur un
large intervalle de richesse moteur.
quasi-dimensionnel est utilisé pour simuler le fonctionnement du moteur pour douze régimes
et les résultats de la simulation sont comparés avec ceux obtenus pour le moteur en
exploitation. Il y a un très bon accord entre la simulation et l’expérimentation pour tous les
rapports du mélange hydrogène/gaz naturel utilisés, ce qui confirme l'efficacité du modèle
proposé.
Ma F. et al [82] ont étudié les caractéristiques de la combustion et les émissions d'un
moteur à allumage par étincelle à injection dans la tubulure d’admission avec différents
rapports de mélange H2/ GNC (0-50% en volume), pour différents moments d’allumage. On
constate alors qu’avec l'augmentation de la fraction d’hydrogène dans le mélange, la valeur du
couple maximal délivré par le moteur diminue et le rendement thermique indiqué augmente.
Le rendement thermique indiqué à faible charge est élevé. La durée de combustion est réduite
en augmentant la fraction d'hydrogène dans le mélange GNCH (gaz naturel carburant enrichi
en hydrogène). Les émissions de NOx, CO et des HC diminuent avec l'augmentation de l'angle
d’avance à l'allumage, et augmentent avec l'augmentation de la charge. Avec le même
moment d’allumage, les émissions de NOx et des CO augmentent avec l'augmentation du taux
d’enrichissement en hydrogène et les émissions des HC présentent une tendance inverse.
Ma F. et al [83] ont réalisé une étude expérimentale visant l’étude de l'augmentation
de la durée de fonctionnement en mélange pauvre d’un moteur à allumage commandé par
l'ajout d'hydrogène, sur un moteur six cylindres fonctionnant au gaz naturel. Quatre rapports
de mélange H2/ GNC (0%, 10%, 30% et 50% en volume). Les effets de diverses conditions de
fonctionnement du moteur sur l’aptitude au fonctionnement en mélange pauvre ont été
également examinés. Des résultats de tests ont été alors analysés d'un point de vue
combustion. Les résultats prouvent que la limite du fonctionnement du moteur avec mélange
pauvre pourrait être prolongée en ajoutant l’H2 et en augmentant la charge (pression dans le
collecteur d'admission). L'effet de la vitesse de rotation du moteur sur la limite du
fonctionnement du moteur avec mélange pauvre est minime. A petite charge, l'augmentation
de la vitesse de rotation du moteur est bénéfique pour prolonger la limite de fonctionnement
au mélange pauvre, qui n’est pas le cas pour une charge élevée. Les effets de la vitesse de
rotation du moteur sont encore plus faibles lors d’enrichissement du moteur avec de
l’hydrogène.
Ma F. et al [84] ont réalisé une étude expérimentale examinant les effets d’ajout
d'hydrogène sur la variation de cycle par cycle (VCC) d'un moteur à allumage par étincelle,
ayant 6 cylindres et fonctionnant au gaz naturel avec injection commandée par papillon. Deux
types de carburants ont été testés à titre de comparaison : le GNC et le mélange
48
hydrogène/GNC en rapport de 80/20 (en volume). L'ajout d'hydrogène s’est avéré efficace
pour le contrôle des émissions de NOx et des HC imbrulés et afin de maintenir des faibles
variations cycle par cycle dans les moteurs à allumage par étincelle fonctionnant en mélange
pauvre.
Andy Burke et al. [85] ont utilisé un modèle de moteur nouvellement développé muni
de réactions chimiques détaillées, pour prédire la formation de NOx à l’intérieur du cylindre
dans des conditions réelles du moteur. Le modèle se compose de quatre parties: (1) un modèle
de moteur simplifié permettant de prédire la pression maximale dans le cylindre pour les
mélanges CH4/H2 ; (2) un modèle semi-empirique permettant de prédire la pression et les
courbes de consommation de carburant au cours d'un cycle du moteur, qui est basé sur la
réalisation du couple maximal de freinage (CFM); (3) un modèle de la structure de flamme
détaillée afin de prédire la formation de NOx au cours du cycle du moteur avec l'utilisation de
mécanisme de réaction chimique GRI, et (4) un modèle chimique détaillé pour prédire
l'influence de la compression et la détente à l’intérieur du cylindre sur l'évolution des
différentes parcelles de gaz brûlé dans le moteur . Le modèle a été appliqué à l' ensemble du
cycle du moteur de l'admission à l'échappement.
Les principales conclusions sont:
1. Un nouveau modèle dynamique de formation de NOx dans les moteurs alimentés
par le mélange CNG/H2 a été développé et appliqué à un moteur.
2 . Le modèle utilise le mécanisme chimique cinétique GRI 3.0 pour le processus de
combustion, et il tient compte de la structure locale détaillée de la flamme, et la variation de
NOx due à la compression et de la détente des produits de combustion.
3. Le modèle néglige l'influence du mélange des produits de combustion après la
combustion, ainsi que l'influence de la turbulence sur la structure de la flamme.
4 . Les simulations de CH4 et CH4/H2 mélanges a donné des valeurs de NOx qui sont
typiques dans le domaine. Cela est vrai sur une gamme de rapports d'équivalence qui vont du
stoechiométrique au trop pauvre.
5 . Le modèle prédit que des mélanges pauvres de CH4/H2 peuvent produire de très
faibles émissions de NOx .
6. Les résultats indiquent qu'il existe une sensibilité élevée de NOx au rapport
d'équivalence. Par exemple, la variation des émissions de NOx pour les rapports d’équivalence
entre 0.7 et 0.6 peut varier d’un ordre de grandeur .
Les moteurs à allumage commandé alimentés avec du gaz naturel comprimé enrichi
par l’hydrogène (GNCH), ont de nombreux avantages par rapport à ceux fonctionnant avec de
49
l'essence, avec du carburant Diesel ou avec du gaz naturel, en particulier du point de vue
contrôle des émissions. Ma, F. et al. [86] ont mené des études expérimentales pour améliorer
en permanence la configuration du moteur à GNCH et la stratégie de contrôle visant à faire
pleinement usage de ce nouveau carburant. Dans le même but, ils ont présenté un modèle
pour simuler le cycle de fonctionnement des moteurs à GNCH, qui sera applicable pour
différents rapports hydrogène/gaz naturel. Les paramètres fondamentaux du modèle
thermodynamique, le modèle de propagation de la flamme turbulente et l'équation connexes
ont été mis en place. Si on considère que le facteur qui influe le plus la pertinence du modèle
sus-mentionné pour différents rapports hydrogène/gaz naturel, est la vitesse de flamme
laminaire, les méthodes de traitement de la vitesse de combustion laminaire dans le modèle
ont ensuite été décrites un peu plus en détail. Après avoir déterminé les constantes du modèle
par calibration, les résultats des simulations ont été comparés avec les données expérimentales
obtenues pour la pression du cylindre pour différents rapports du mélange hydrogène/gaz
naturel, différents moments d'allumage et différentes richesses du moteur. Les données ainsi
obtenues montrent que les résultats des simulations et des tests convergent assez bien, sauf
pour les mélanges trop pauvres où les problèmes d'une combustion incomplète deviennent
plus sévères.
Afin d'étudier l'influence d’ajout de 0-50 % (en volume) d'hydrogène sur les
performances du moteur au ralenti, Ma, F. et al. [87] ont conduit une étude expérimentale sur
un moteur six cylindres alimenté avec du gaz naturel avec injection commandée au moyen
d'un papillon. Des tests ont été effectués sous divers rapports d'excès d’air et moments
d'allumage. Les résultats prouvent que l'ajout d'hydrogène a diminué d’une façon
remarquable l'émission du méthane non brulé (CH4), tandis qu'elle n'a exercé aucun effet
significatif sur la réduction du CO. Les émissions de NOx sont relativement basses au ralenti
comparativement au reste des émissions. L'ajout d'hydrogène combiné avec le retard à
l’allumage s’est avéré efficace pour réduire les émissions au ralenti. Le coefficient de
variation de pression effective moyenne indiquée et le taux de combustion partielle
pourraient être simultanément réduits par l'ajout d'hydrogène, ce qui a mené à l'amélioration
de la stabilité du ralenti. En conséquence, l'amélioration de la stabilité du ralenti a diminué la
consommation de carburant. On conclut ainsi que l’ajout est une approche efficace et
pratique pour améliorer la stabilité du ralenti et pour diminuer les émissions.
Ma F. et al. [88] dans leur étude ont employé des méthodes d’analyse statistique pour
optimiser un moteur à allumage par étincelle alimenté par des mélanges de gaz naturel et
d'hydrogène. Les paramètres de contrôle optimaux sont déterminés pour être la fraction
50
Puisque le méthane est le principal constituant du gaz naturel, nous élargirons notre
étude bibliographique aux mélanges méthane/hydrogène.
Kahraman et al. [90] lors de leur étude expérimentale d’un moteur à allumage
commandé (4 cylindres) alimenté par des différents rapports CH4-H2 (0% H2, 10% H2, 20%
H2 et 30% H2) ont relevé qu’en augmentant la fraction d’hydrogène dans le carburant, le peak
de la pression du cylindre se rapproche de la position du point mort haut et que les émissions
HC, CO et CO2 sont réduites avec l’augmentation de la fraction d’hydrogène dans le mélange.
51
F.H.V. Coppens et al. [91] ont effectué des mesures de formation des oxydes d’azote
dans les flammes du mélange hydrogène + méthane + air sous pression atmosphérique et
température ambiante de 298 K. La concentration de NOx dans ces flammes a été mesurée
dans les gaz brûlés en utilisant un échantillonnage de la sonde. Dans les mélanges méthane +
hydrogène + air, les concentrations de NOx ont été mesurées à partir d’une distance fixe du
brûleur en fonction de la richesse. Deux maxima en concentration ont été observés: en
mélange stoechiométrique, dû au mécanisme thermique de formation de NOx de Zeldovich et
en mélange riche, avec richesse autour au 1.3, dû au mécanisme de formation de NOx de
Fenimore. En mélange pauvre, l’ajout d’hydrogène possède un effet mineur sur les émissions
de NOx; cependant, en mélange riche, ces émissions diminuent d’une manière significative.
Les calculs numériques concordent bien avec les tests.
Akansu et al. [92] ont aussi testé un moteur 4 cylindres à allumage commandé
fonctionnant avec un mélange méthane/hydrogène en proportion : 0, 10%, 20% et 30% d’H2
en volume, en faisant varier la richesse globale de 0.6 à 1.2. Ils ont conclu que :
• Les émissions des NOx augmentent avec l’ajout d’hydrogène. Les émissions des
CO et des CO2 diminuent avec l’augmentation de la fraction d’hydrogène dans le mélange
carburant.
• Le rendement thermique, dans le cas des mélanges pauvres, augmente avec l’ajout
d’hydrogène.
Fairweather et al. [93] ont étudié l’effet de l’ajout d’hydrogène sur la vitesse de
flamme turbulente dans une enceinte. Ils ont observé que dans les conditions
stoechiométriques, le rapport entre la vitesse turbulente et la vitesse laminaire de combustion
ne varie pas avec l’ajout d’hydrogène pour une intensité turbulente fixée.
Yusuf dans sa thèse de Maitrise a utilisé un rapport de 99% de méthane (CH4) et
80.08/19.92 du mélange CH4/ H2 [94]. Il a testé un moteur Nissan quatre cylindres de type
510, ayant : la vitesse de rotation maximale de 5200 tr/min, l’alésage 85 mm, la course 86
mm, le taux de compression de 8.5:1, 1952 cm3 et la puissance maximale de 92 ch. Il a
analysé le taux de croissance du front de la flamme (TCFF) et le moment d’allumage pour
différentes richesses. La richesse de 0.535 était choisie à cause de la combustion du méthane
avec mélange pauvre. L'analyse des résultats montre que l'ajout d'une quantité d'hydrogène au
méthane compris entre 0.4 et 0.6, avec un rapport d'équivalence proche du stoechiométrique
et une vitesse de rotation du moteur située entre 2000 et 3000 tr/min, produit des
améliorations importantes des performances du moteur et des émissions.
52
Yusuf [95] a testé un autre moteur Toyota quatre cylindres de type 2TC, ayant comme:
alésage de 85 mm, course de 70 mm, taux de compression de 9.0:1, 1.6 litres, 1588 cm3,
puissance maximale de 88 hp et vitesse de rotation maximale de 6000 tr/min. Le moteur
fonctionnant avec un mélange CH4/ H2 a été testé pour une vitesse de rotation de 1000 tr/min,
avec le meilleur rendement d’avance à l’allumage et avec la charge minimale. Quand le
mélange méthane/hydrogène est comparé avec du méthane pur pour les mêmes richesses, on
remarque une augmentation du rendement thermique et des émissions des NOx . Alors que le
rendement dimunue au détriment des concentrations des hydrocarbures imbrulés et du CO. De
plus, la limite inférieure de la combustion du gaz naturel était réduite de 0.61 jusqu’a 0.54.
Cette limite de combustion était définie comme étant opérationnelle au moins aux environs de
38% des cycles de la combustion non complétée. En ajoutant de l’hydrogène, les richesses
devraient être réduites de 15 % sans réduction de la durée de la combustion et de l’allumage.
Hoekstra et al. [96,97] ont testés un moteur à allumage par étincelle, 17CV, 1700
tr/min avec des rapports de mélanges CH4/ H2: 100/0, 89/11, 80/20, 72/28 et 64/36 en faisant
varier les richesses. Ils ont constaté que l’ajout de 20% (en volume) d’hydrogène dans le
carburant permet l’extension de la limite de stabilité de la richesse de 0.75 (sans hydrogène) à
0.67. En outre, les émissions des NOx augmentent et celles des hydrocarbures diminuent avec
l’augmentation da la richesse. Néanmoins, l’ajout de 28% ou de 36% (en volume)
d’hydrogène à des mélanges pauvres (φ<0.8), permet d’obtenir des valeurs d’émissions NOX
assez basses (<250ppm) avec des émissions HC modérées.
Liu et al. [98] ont étudié les caractéristiques de la détonation dans le mélange CH4/ H2
avec des proportions en % : 100/0, 90/10, 70/30, 50/50, 30/70 et 0/100. Ils ont simulé des
modèles mathématiques et comparé les données reçues avec l’étude ultérieure sur les
caractéristiques de la détonation. Ils ont constaté que si l’hydrogène et le méthane possèdent
une température d’admission élevée, les régions de la détonation correspondent aux valeurs
expérimentales. Cependant, l’hydrogène manifeste une grande déviation dans les données
expérimentales pour des températures d’admission basses, à cause de la haute énergie issue de
la propagation rapide de la flamme au cours de la combustion de l’hydrogène.
Karim [99] a étudié les limites de la détonation associées avec des mélanges contenant
différents pourcentages de H2 et de CH4, différentes richesses et différentes températures. Il a
constaté que lorsque l’hydrogène et le méthane sont mélangés relativement en petites
quantités, l’excellente résistance à la détonation du méthane n’est pas amoindrie.
Karim et al. [100] ont fait une étude expérimentale sur un moteur CFR avec un
mélange méthane/hydrogène (CH4/H2) en proportion (en %): 100/0, 90/10, 80/20, 70/30,
53
60/40, 50/50, 40/60, 30/70 et 20/80, en faisant varier les richesses et l’avance à l’allumage.
Pour 10°APMH et 20°APMH, avec l’augmentation de la concentration d’hydrogène dans le
moteur, la puissance augmente; mais pour 30°APMH, avec l’augmentation de la
concentration d’hydrogène dans le moteur, la puissance diminue. La puissance est maximale
aux environs de 20°APMH. Le moment d’allumage possède une influence défavorable sur le
rendement thermique. La région de détonation a été définie pour différents taux de
compression et différents rapports méthane-hydrogène. Si une certaine quantité d’hydrogène
est ajoutée au méthane comme carburant du moteur à allumage par étincelle, les performances
de ce moteur chutent d’une manière drastique.
Blarigan et Keller [101] ont étudié le mélange CH4/H2 en proportion (%): 100/0,
70/30 et 0/100. Ils ont réalisé des tests et des simulations avec un code trois-dimentionnel
CFD. Ils ont reporté qu’il fût possible de concevoir une unité de puissance auxiliaire à haut
rendement avec émissions équivalentes à zéro, pour un véhicule hybride, avec de l’hydrogène
ou avec un mélange de 30% d’hydrogène et 70% de gaz naturel.
Wong et al. [102] ont étudié les effets cinétiques de la présence de certains carburants
gazeux et des produits de la réaction du préallumage avec l’hydrogène dans les moteurs. Il ont
démontré que la présence de ces combustibles avec l'hydrogène peut influencer de manière
significative les vitesses des réactions responsables de l’auto-allumage, qui contrôlent la
tendance des moteurs à allumage commandé à s’auto-allumer. Il ont également montré que la
recirculation des gaz d'échappement peut modifier de manière significative le comportement
d'auto-inflammation de l'hydrogène et de ses mélanges avec du méthane ou du propane dans
les moteurs.
Shrestha et al. [103] ont étudié de leur coté le mélange CH4/H2 en proportions (en %):
100/0, 90/10, 80/20, 70/30, 20/80 et 10/90 avec différents taux de compression, en faisant
varier les richesses. Ils ont remarqué que l’ajout d’une certaine quantité d’hydrogène au
méthane dans le moteur à allumage commandé améliore les performances du moteur,
particulièrement lors du fonctionnement avec des mélanges à richesse relativement faible. La
concentration optimale d’hydrogène dans le mélange combustible qui peut mener à un gain en
puissance et éviter la détonation, est de l’ordre de 20–25% en volume, sur la plage des
conditions considérées.
Shudo et al. [104] ont analysé la combustion et les émissions dans un moteur alimenté
par un mélange pauvre de méthane et d’hydrogène dans un moteur à injection directe à charge
stratifiée prémélangée. Les tests ont été exécutés sur un moteur à allumage par étincelle
monocylindrique à quatre temps, avec D × L: 85 × 88mm et ε=13. Les résultats ont montré
54
Nagalingam et al. en 1983 [108]. Ils ont étudié l’ajout d'hydrogène au méthane dans les
proportions suivantes: 0, 20, 50 et 100%, pour une seule vitesse de rotation du moteur. Ils ont
conclu que l'enrichissement du méthane par l'hydrogène contribue à la prolongation de la
limite de la combustion du mélange pauvre due à sa nature même, mais il diminue la
puissance à cause du bas pouvoir calorifique volumétrique global. Le rendement thermique
indiqué a diminué avec l'ajout d’hydrogène, probablement dû à un rapport décroissant de la
puissance au frein, à la puissance de frottement et à la valeur croissante du pouvoir calorifique
du carburant. Le point d’allumage optimal a baissé jusqu'à 20° avant le point mort haut pour
des mélanges d’hydrogène pur indiquant une vitesse de la flamme élevée. Les émissions des
NOx ont augmenté pour l'hydrogène pur dû à la température élevée de combustion, avec des
valeurs maximales se produisant pour une richesse plus faible. Les émissions des HC
imbrûlés et de CO ont diminué. Toutes les tendances se sont comportées comme prévu, avec
une puissance maximale au alentour de la richesse stoechiométrique, un rendement et des
émissions de NOx maximaux pour des richesses de 0.8-0.9, et avec des émissions des HC qui
atteignent une valeur minimale pour des richesses allant de de 0.8 jusqu’à 0.9. Le CO a été
seulement produit par des mélanges ayant la richesse au-dessus de 0.95.
En 1992, une étude a été réalisée à l'université de Melbourne sur un moteur
monocylindrique pour une vitesse de rotation fixe de 1200 tr/mn [109]. Le but de l'étude était
de réduire la variabilité cycle à cycle et l’enrichissement du gaz naturel à l'allumage par de
l'hydrogène. Le processus se réfère à un système d’allumage à jet assisté d’hydrogène (AJAH)
et à l’emploi d’une préchambre pour contrôler ce jet. La quantité d'hydrogène utilisée varie de
2 à 10% en masse du méthane. L'ajout d'hydrogène augmente le rendement thermique, la
puissance et les émissions des NOx, tout en diminuant le point d'avance à l’allumage optimal
et les émissions des hydrocarbures imbrûlés. Cependant, il est difficile de dire si c'est dû
uniquement aux propriétés de l'hydrogène, ou à la variation dans la stoechiométrie de la
combustion.
En 1995 Cattelan et al. [110] ont réalisé une étude sur un moteur 3.1 litres, six
cylindres avec un mélange de 15% d'hydrogène et 85% de méthane, en faisant varier la
vitesse de rotation, la charge et le moment d’allumage, pour des richesses allant de 1 jusqu’à
1.03. Les résultats ont prouvé que l’ajout de 15% d'hydrogène contribue généralement à
l’amélioration du rendement, à la diminution des émissions des hydrocarbures, de CO et de
CO2 et à l’augmentation des émissions de NOx.
En 1996, une étude complète a été exécutée par Collier et al. [111] sur un moteur 8
cylindres, 4.6 litres. On a étudié l’ajout d'hydrogène dans le méthane en proportion de (en %):
56
Fig. 3.10 Délais d’allumage mesurés pour mélange avec 50% d’H2 en fonction de la
pression avec température constante Tc = 995 ± 4 K et φ = 1.0, φ = 0.5 (N.B. Anikin et al
[113]).
Park et al. [116] ont réalisé une étude numérique sur les effets d’ajout d’hydrogène et
de la vapeur sur une flamme de diffusion à contre-courant méthane-air (fig. 3.11) avec le
mécanisme GRI-Mech 3.0 (Smith et al., [117]).
Les auteurs ont constaté que l’ajout d’hydrogène conduisait à une augmentation de la
température de la flamme. Ils ont adopté les indices des émissions du CO, CO2 et de NOx
proposés par Takeno et Nishioka[118].
58
𝐿
∫0 𝑊𝑖 𝜔̇𝑖 𝑑𝑥
𝐸𝐼𝑖 = 𝐿 (3.4)
∫0 𝑊𝑐𝑎𝑟𝑏 𝜔̇𝑐𝑎𝑟𝑏 𝑑𝑥
Scholte et Vaags en 1959 [119] et Miller et al. en 1963 [120] étaient les premiers
chercheurs qui ont réalisé des travaux sur les vitesses de propagation de la flamme
méthane/hydrogène. Scholte et Vaags [119] ont mesuré les vitesses de flammes laminaires
méthane/hydrogène/air stabilisées sur un brûleur. Les auteurs ont constaté que les vitesses de
flamme augmentent considérablement avec l’ajout d’hydrogène. Miller et al. [120] ont étudié
l’influence d’ajout de 0.5% par volume de méthane sur la vitesse de flamme hydrogène/air par
imagerie de particules d’oxyde de magnésium ensemencées. L’effet d’ajout d’hydrogène
devient significatif seulement quand la fraction d’hydrogène dans le mélange atteint les 30%.
Liu et al. [121] ont mesuré les vitesses de flamme méthane (0-100%)/ hydrogène par la
méthode du tube pour une gamme de richesse allant de 0.6 jusqu’à 1.3. Tanoue et al. [122]
ont mesuré les vitesses de flamme laminaire étirée d’un mélange méthane /(0-40%)
hydrogène dans une bonbonne sphérique. Toutefois, ces auteurs n’ont pas pris en compte
l’étirement du front de flamme et leurs mesures ne peuvent pas donc être directement
comparées aux autres valeurs de la littérature.
Yu et al. [123] ont réalisé une étude expérimentale des vitesses des flammes
laminaires du mélange méthane/air et propane/air avec et sans addition de petites quantités
59
Fig. 3.12 Évolution des vitesses de propagation de flamme des différents mélanges
méthane/hydrogène/air en fonction de la richesse (Yu et al. [123]).
Une corrélation a été proposée pour décrire l’évolution des vitesses de propagation de
la flamme des mélanges méthane/hydrogène/air en fonction de la quantité d’hydrogène
substituée au méthane:
Dans cette équation, SL0 (φ = φF ) et SL0 (φF , R H ) sont respectivement les vitesses de flamme
sans et avec hydrogène, RH représente la fraction molaire d’hydrogène présent dans le
combustible.
Ren et al. [125] ont mesuré les vitesses de flamme de mélanges méthane/hydrogène à
pression atmosphérique dans les conditions de combustion en mélange pauvre (φ = 0.73, 0.68
et 0.63) pour une fraction molaire d’hydrogène limitée (0 à 8%) dans le combustible. Une
augmentation de 20% est obtenue entre une flamme méthane (22 cm.s-1) et une flamme
mélange méthane/hydrogène CH4/8%H2 (27 cm.s-1) à une richesse de 0.73. Ces valeurs
expérimentales convergent avec les prédictions du mécanisme GRI 3.0 (Smith et al.[117]).
Law et Kwon [124] ont réalisé une étude expérimentale de la propagation des flammes
sphériques prémélangées allumées par étincelle d’un mélange méthane/hydrogène/air dans
une chambre de combustion à pression constante dans les conditions TPN (température et
pression normales). Ce travail a été réalisé afin d’évaluer l’eficacité de la substitution partielle
d’une partie de l’hydrogène au méthane (5, 10 et 15%), d’améliorer la sécurité d’utilisation
d’hydrogène en général et les performances du moteur à combustion interne en particulier
ainsi que la stabilité et la propagation d’une flamme hydrogène/air. L'étude expérimentale a
été réalisée à l'aide de la chambre de combustion à pression constante décrite dans [126].
Trois différentes richesses ont été étudiées (φ = 0.6; 1 et 1.67). La séquence de propagation de
la flamme a été visionnée par strioscopie et enregistrée à l’aide d’une caméra numérique à
grande vitesse de mouvement, jusqu'à 8000 images/s. Les mesures ont été limitées aux
flammes ayant des rayons plus grands que 5mm pour contre le mauvais allumage et plus
petits que 20mm pour éviter les interférences du mur. La propagation de la flamme plane, non
étirée, stable et adiabatique dans le domaine infini a été simulée en utilisant le code PREMIX
de la flamme premélangée, laminaire, unidimensionnelle et stable (Kee et al. [127]). Tous les
calculs ont été effectués avec le mécanisme d'oxydation des hydrocarbures détaillée du Qin et
al. [128], ce qui permet la diffusion à plusieurs composants, la diffusion thermique et les
propriétés de transport variables [129]. Les auteurs ont constaté que les vitesses de
propagation de flamme diminuent respectivement de 17, 30 et 40% lorsqu’on remplace 5, 10
et 15% d’hydrogène par du méthane dans le cas des flammes pauvres et stoechiométriques.
L’effet du méthane est beaucoup plus important dans les conditions riches en combustible.
Par exemple, la vitesse de flamme est réduite de 66% lorsqu’on remplace 15% d’hydrogène
par du méthane. Les prédictions du mécanisme de Qin et al. [128] sont satisfaisantes pour les
mélanges stoechiométriques et riches, mais surestiment les valeurs expérimentales des
mélanges pauvres.
61
Halter et al. [130] ont réalisé une étude expérimentale sur les effets de la pression (1, 3
et 5 atm) et de la substitution de 10 et 20% d’hydrogène sur les vitesses de propagation de la
flamme méthane/air. Le dispositif expérimental est constitué d'une chambre de combustion
sphérique. La visualisation de la flamme est obtenue par la méthode d'ombroscopie classique.
La lumière parallèle a été créée par trois lentilles planes convexes (ayant des distances
focales : 104 mm (L1), 300 mm (L2) et 22 mm (L3)) et de deux miroirs sphériques (SM1 et
SM2,
f = 500 mm), comme indiqué sur la Fig. 3.13. Les images ont été enregistrées à l'aide d'une
caméra vidéo à grande vitesse (Phantom V5) fonctionnant à 11.200 images /sec avec un
temps d'exposition de 50 ms.
Les résultats expérimentaux ont été comparés aux calculs en utilisant un schéma
cinétique chimique détaillé (GRIMECH 3.0). Les profils d’évolution des vitesses
expérimentales de flammes prémélangées CH4/H2/air présentent un maximum obtenu au
voisinage d’une richesse égale à 1.1. Une vitesse maximale de 37.9 cm.s-1 est obtenue pour un
mélange CH4/air à pression atmosphérique; de 39.3 cm.s-1 pour un mélange CH4/10%H2/air et
de 42.4 cm.s-1 pour un mélange CH4/20%H2/air. À la stoechiométrie, la vitesse d’une flamme
CH4/10%H2 à 1 atm est de 40.81 cm.s-1, de 23.8 cm.s-1 à 3 atm et de 21.1 cm.s-1 à 5 atm.
Enfin, lorsque la pression augmente, la vitesse de combustion laminaire diminue pour tous les
mélanges. Le domaine de pression est limité à 0.5 MPa en raison de l'apparition des
instabilités à des pressions supérieures à cette valeur. . L’augmentation est prédite par le
mécanisme GRI 3.0 (Smith et al. [117]), mais les valeurs modélisées sont surestimées.
Coppens et al. [91] ont réalisé des mesures de la vitesse de flamme adiabatique dans
des mélanges (CH4 + H2)/O2/N2. On a fait varier la fraction d’hydrogène dans le combustible
de 0 à 35%. Des flammes non étirées ont été stabilisées sur un brûleur à plaque perforée à
pression atmosphérique (fig. 3.14).
62
Fig. 3.14 Le brûleur utilisé dans le procédé de flux thermique (Coppens et al. [91]).
Fig. 3.15 Évolution des vitesses en fonction de la richesse pour des mélanges CH4/H2 dans
l’air. Croix : 35% H2, cercles : 25% H2, triangles : 15% H2, carrés : 5% H2, diamants : 0% H2;
traits : modélisation (Coppens et al. [91])
Mandilas et al. [131] ont réalisé une étude expérimentale afin d’étudier l’influence de
l’ajout de 30% d’hydrogène au méthane sur la vitesse laminaire de flamme. Les mesures ont
été effectuées à une pression de 5 bars et la richesse a été variée de la limite d’allumage dans
le mélange pauvre (φ = 0.6) à la limite d’allumage dans le mélange riche (φ = 1.4). La vitesse
63
de flamme de CH4/H2 est supérieure à celle du méthane pur dans les conditions de mélange
pauvre et stoechiométrique allant même jusqu’à doubler pour un mélange très pauvre. Par
exemple, la vitesse maximale du méthane est de 28 cm.s-1 et celle du mélange CH4/H2 est de
38 cm.s-1. Par contre, la vitesse de flamme de méthane/hydrogène est proche de celle du
méthane pour des richesses supérieures à 2 d’après les auteurs.
L'augmentation de la vitesse de flamme laminaire est principalement responsable de
l’augmentation de vitesse de flamme turbulente lors d'ajout d'hydrogène. Pour les mélanges
pauvres iso-octane/air, l’addition d'hydrogène entraîne le rehaussement de la vitesse de
flamme turbulente au-delà de celle prévue sur la base de la vitesse de flamme stratifiée.
Dans le Tableau 3.5, on présente une synthèse des conditions expérimentales des
différentes études réalisées sur les vitesses de flammes de propagation de méthane en
présence d’hydrogène.
Tableau 3.5. Récapitulatif des conditions expérimentales des travaux relatifs aux vitesses de
flammes méthane/hydrogène
Miller et al.,
Brûleur 93 φ = 1.75 298 1
[120]
Halter et Enceinte
0.7 ≤ φ ≤ 1.2 298 1
al.[130] sphérique 0, 10, 20
1
Coppens et al., 5, 15, 25,
Brûleur 0.7 ≤ φ ≤ 1.5 298
[91] 35
Fig. 3.16 Comparaison des vitesses de flammes expérimentales des mélanges CH4,
CH4/10%H2, CH4/20%H2, CH4/40%H2, Tini = 298 K, P = 1 atm, (1 : Yu et al. [123], 2 : Halter
et al. [130], 3 : Coppens et al. [91]).
Cheng et al. [114] ont réalisé une étude expérimentale sur les délais d'auto-
inflammation de onze mélanges méthane/hydrogène/oxygène dilués dans de l’argon en
utilisant la technique de choc réfléchi (tube à choc). Les conditions expérimentales couvrent
des températures de 800 à 2400 K, des pressions de 1 à 3 atm et des richesses de 0.5 à 1.5.
Une corrélation de type Arrhenius généralisée est obtenue par l'ajustement des données
expérimentales avec une méthode de régression non linéaire pour le méthane et pour
l’hydrogène pur. On a observé que l’énergie d’activation globale du mélange est
proportionnelle au rapport entre la concentration du CH4 et du H2. Ils ont proposé une
expression pour le délai d'auto-inflammation du mélange méthane-hydrogène:
1−𝜉 𝜉
𝜏 = 𝜏𝐶𝐻4 · 𝜏𝐻2 , (3.6)
Fig. 3.17 Comparaison des délais d’allumage pour le mélange stoechiométrique méthane-air
en utilisant l’équation d’Oppenheim avec différentes concentrations d’hydrogène à 1 atm
(Cheng et al. [114])
67
Fig. 3.18 Délais d’auto-inflammation mesurés et calculés du mélange testé. Symboles qui
désignent les résultats expérimentaux: ( ), méthane pur (φ =1), ( ), mélange #1
(85%CH4/15%H2), ( ), mélange #2 (65%CH4/35%H2). Les lignes représentent le délai
d’auto-inflammation du modèle actuel. Une erreur horizontale ( ̶ ) représente l'incertitude
typique dans la température expérimentale. (Huang et al. [133]).
Hadj Ali [136] a mesuré en tube à choc (fig. 3.19) les délais d’auto-inflammation du
mélange CH4/10%H2 pour différentes richesses (φ = 0.5, 1 et 2) et pour une pression
𝑝 = 10 𝑎𝑡𝑚 et du mélange CH4/40%H2 pour des richesses (φ = 0.2, 0.5 et 1) et une pression
𝑝 = 14 𝑎𝑡𝑚 pour des températures comprises entre 1250 et 2000 K. Les résultats
expérimentaux obtenus avec les deux mélanges méthane / hydrogène (10 et 40% H2) ont servi
à l'établissement d’une corrélation exprimant le délai d'auto-inflammation en fonction de la
concentration de chaque composant du mélange (exprimée en mol/m3) et de l'énergie
d'activation (exprimée en J/mol). R étant la constante des gaz parfaits (8.314 J/mol.K) et T la
température en Kelvin.
Cette corrélation est valable dans l'intervalle de température [1340-2040 K], dans la gamme
de pression [8.7-15.1 bar] et la gamme de richesse [0.2-1]. Le coefficient de corrélation, R2,
est de 0.938. Hadj Ali a constaté que les délais mesurés à 1625 K augmentent
considérablement lorsque la richesse augmente, mais l'effet est moins marqué à 1777 K.
69
Fig. 3.19 Dispositif expérimental du tube à choc en acier inoxydable (Hadj Ali [136]).
Chaumeix et al. [137]) ont réalisé une série de mesures en tube à choc couplé à la
spectrométrie d'émission, sur les délais d’auto-inflammation de mélanges CH4/(0, 20% et
60%)H2/O2, pauvres (φ = 0.4 ÷ 0.75) et stoechiométriques ( φ = 1) fortement dilués dans
l’argon. La température varie de 1250 à 2000 K et la pression derrière les ondes de choc
réfléchies se situait entre 0.15 et 1.6 MPa. Les délais d’auto-inflammation varient de façon
exponentielle avec la température. Les longs délais d’auto-inflammation sont obtenus pour
des mélanges pauvres: le délai d’auto-inflammation d’un mélange CH4/60%H2 est de 300 μs à
1650 K et à une richesse du mélange de 0.4, de 610 μs à une richesse de 0.7 et à la
stoechiométrie du mélange, il est égal à 900 μs . Une corrélation de régression multiple est
appliquée sur les valeurs expérimentales pour obtenir une expression d’Arrehenius modifiée
du délai d’auto-inflammation τi, en fonction des conditions initiales de température, pression
et concentrations.
de fonctionnement des températures et des pressions utilisées dans cette étude, les différences
entre les valeurs mesurées et calculées du délai d’auto-inflammation sont inférieures à 10%
pour les combustibles purs et supérieures à 25% pour les mélanges hydrogène / méthane.
Très peu d’études expérimentales ont été publiées sur la structure de flammes
prémélangées méthane/hydrogène. Néanmoins, on peut citer les travaux de Ren et al. [125],
qui ont étudié la dynamique et la structure de mélanges pauvres CH4/air. Le dispositif
expérimental est une configuration d’écoulement à un seul jet dans un plan de stagnation,
variante d’une configuration à jets opposés. La figure 3.21 illustre le schéma de la
configuration, qui a également été utilisé dans des études précédentes [147]. Les flammes sont
pratiquement en conditions adiabatiques et la température est maintenue aux environs de 900
K. Les auteurs présentent une structure de flamme CH4/8%H2 de richesse ϕ = 0.64 peu
détaillée. Les profils de fraction molaire de H2O, O2, CH4, CO2, H2 sont présentés (Figure
3.22, référence [125]). En général, un bon accord entre les données expérimentales et les
simulations numériques a été trouvé, renforçant ainsi la validité du mécanisme GRI 3.0 décrit
dans [117] pour ce mélange combustible. (GRI-Mech 3.0 est un mécanisme optimisé conçu
pour modéliser la combustion du gaz naturel, y compris la formation de NOx et de la chimie
de rebrûlage).
Fig. 3.21 Schéma de l’écoulement à un seul jet dans un plan de stagnation (Ren et al.
[125]).
72
Fig. 3.22 Profils de fractions molaires expérimentales (symboles) et modélisées (traits) des
espèces H2O, O2, CH4, CO2, H2 dans une flamme CH4/8%H2, φ= 0.65, (Ren et al., [125]).
Bagdanavicius A. et al. [143] ont réalisé une étude sur l'effet d’ajout d'hydrogène dans
un moteur fonctionnant au méthane sur la flamme turbulente prémélangée. Les flammes
pauvres, stoechiométriques et riches ont été testées à l'aide d'un brûleur Bunsen à pressions
absolues de 3 et 7 bars, et des températures de 473, 573 et 673 K. Les mesures d'épaisseur de
la queue de la flamme ont montré que la pression, la température et l'ajout d'hydrogène ont
peu d'effet sur la structure du méthane et du mélange méthane/d'hydrogène pour le brûleur
Bunsen testé. Deux importantes caractéristiques de la structure de flamme turbulente, à savoir
l’épaisseur de flamme δT et la densité de surface de flamme Σ (air de flamme par unité de
volume), ont été quantifiées et analysées. La densité de surface de flamme et le coefficient du
modèle Bray-Moss-Libby g/(σyLy) pour le mélange étudié, ont été déduites. On n’a pas
observé un effet significatif de la température, la pression et de la richesse.
Gauducheau et al. [144] ont mené une étude numérique des structures de flammes
libres de prémélange méthane/hydrogène/air. Les profils des espèces ont été calculés pour les
flammes (90%CH4/10%H2)/O2/N2 et (80%CH4/20%H2)/O2/N2 pour une richesse de 0.6 et à
une pression de 30 atm et une température initiale de 700 K. Le mécanisme GRI 2.11
(Bowman et al., [135]) a été utilisé. Les profils des espèces et de la température sont restés
presque identiques à ceux de la flamme méthane/air. Par contre, les concentrations des
radicaux H, OH et O sont légèrement plus élevées (10%) dans le cas d’un mélange contenant
de l’hydrogène. Pour ce qui est des produits de combustion, on constate une diminution de
12% de CO2, de 15% de CO et une augmentation de 4% de H2O. Cependant, les auteurs n’ont
pas discuté ces résultats.
73
Tseng [145] a réalisé une étude numérique, en utilisant le mécanisme GRI 2.11
(Bowman et al., [135]), de l’effet de l’hydrogène sur la température maximale de flamme dans
le cas de mélanges CH4/H2/air pauvres (φ = 0.8) contenant 60% d’hydrogène, dans les
bruleurs à milieu poreux. Toutefois, l'ajout d'hydrogène au méthane ne semble pas avoir
d'importants effets sur la température maximale de la flamme. L'énergie libérée dans la zone
de flamme par la réaction chimique est équilibrée par la convection et la conduction. Bien que
l'hydrogène a des températures de flamme adiabatique supérieures au méthane, il a aussi une
conductivité plus élevée. De plus, la convection augmente également dans le cas de
l’hydrogène en raison de la vitesse de flamme plus élevée. Par conséquent, la température
maximale de la flamme varie peu quand on ajoute de l’hydrogène au méthane dans les
brûleurs à milieu poreux. Les profils de fraction molaire des radicaux H et OH sont aussi
présentés pour les flammes CH4/air et CH4/60%H2/air. Les concentrations des radicaux dans
la flamme du brûleur en milieu poreux sont plusieurs fois plus élevées que celles de la flamme
libre. Ces fortes concentrations de radicaux peuvent conduire à l'accélération de la vitesse de
réaction et à l’amincissement de la zone de réaction, et donc à l’augmentation de la vitesse de
flamme. Tseng n’a pas présenté que les profils de fraction molaire de H et de OH (Tableau
3.6).
Tableau 3.6 Compositions initiales des structures de flammes CH4/H2 et espèces analysées
(Tseng [145])
Profils de
Auteurs Flamme φ CH4 H2 O2 N2 fraction
molaire
Gauducheau 15 espèces
CH4/air 0.6 0.0591 0 0.1971 0.7437
et al., [144] ( à 30 atm)
Gauducheau 15 espèces
CH4/20%H2/air 0.6 0.0551 0.0138 0.1951 0.7361
et al., [144] ( à 30 atm)
Tseng,
CH4/air 0.8 0.0773 0 0.1933 0.7294 H, OH
[145]
Tseng,
CH4/60%H2/air 0.8 0.0349 0.0523 0.1918 0.7211 H, OH
[145]
Ren et al., O2, CH4, H2,
CH4/8%H2/air 0.62 0.0596 0.0052 0.1965 0.7387
[125] H2O, CO2
74
Uykur et al. [146] ont étudié les profils de température adiabatique de flammes
laminaires prémélangées méthane/air stoechiométriques et pauvres (φ = 0.6) avec 20%
d’hydrogène en utilisant les méthodes de simulation cinétique et chimique. Le logiciel de
simulation cinétique CHEMKIN a été utilisé avec le mécanisme cinétique GRI 3.0 (Smith et
al., [117]) à pression atmosphérique. Bien que tous les mélanges semblent avoir des profils de
température similaires, les températures finales (températures, à la fin du domaine de calcul)
varient pour les différents mélanges comme on le voit sur la Fig. 3.23. Bien que l'hydrogène
ait une température de flamme plus élevée par rapport au méthane, la température de flamme
du mélange a été augmentée de moins de 1%, même avec l'ajout de 20% de H2. L’évolution
de la fraction molaire finale de CO et de NOx a été calculée pour les valeurs suivantes de la
richesse: 0.6 ≤ ϕ ≤ 1.4. L'ajout de 10 à 20% d'hydrogène dans le mélange n’a qu’un simple
effet sur l'amélioration de la vitesse de la flamme. Cependant, l'ajout d'oxygène et
d'hydrogène dans le même rapport que l'on trouve dans l'eau s'est avéré bénéfique.
L'amélioration de la vitesse de flamme du mélange méthane / air par l'addition de 10%
d'hydrogène et de l'oxygène était équivalente aux améliorations obtenues par l'ajout de 20%
d'hydrogène seulement. Dans les mélanges proches de stoechiométriques, l'ajout d'oxygène
augmente sensiblement les concentrations de NOx, mais pour les mélanges pauvres, on ne
prédit pas d'augmentation de NOx.
On peut citer le travail expérimental effectué par İsmet Çelikten, Can Çinar et E. İshak
Yildiz [148] sur l’effet d’ajout d’hydrogène au gaz de pétrole liquéfié (GPL) (0%, 5% 10% et
15% en Volume), du point de vue performances et émissions dans un moteur quatre cylindres
à allumage par étincelles (voir fig.3.24). L’essence a été utilisée comme carburant de base
pour la comparaison. Les tests ont été réalisés pour une plage de vitesse de rotation allant de
1800 t/min jusqu’à 3200 t/min, avec papillon complètement ouvert.
Fig. 3.24 Schéma des équipements expérimentaux (İsmet Çelikten et al. [148])
Fig. 3.27 Variation des émissions du CO en fonction de la vitesse du moteur (İsmet Çelikten
et al. [148]).
Fig. 3.28 Variation des émissions des HC en fonction de la vitesse du moteur (İsmet Çelikten
et al. [148]).
Choi et al. [149] ont étudié les effets d’ajout d’hydrogène dans un moteur fonctionnant
au gaz de pétrole liquéfié (GPL), sur les émissions des gaz d’échappement, sur l’efficacité
thermique et sur les performances. L’augmentation du taux d’hydrogène diminue la pression
effective moyenne (PEM), l’efficacité thermique ainsi que les émissions des hydrocarbures.
78
L’ajout de 20% d’hydrogène favorisera la diminution des émissions des NOx de 20% en
comparaison avec la combustion avec pur GPL.
• Des températures suffisamment basses qui ne nécessitent pas l’utilisation des alliages
couteux.
La figure 3.29 montre le cycle Bourque sous forme de diagramme. L'écoulement d'air
est montré en vert et celui de la vapeur est en bleu. La figure 3.30, extraite du brevet original,
montre le cycle dans le plan température-entropie. Le côté gauche de la figure 3.30 montre le
cycle de la vapeur à pleine puissance. Le côté droit de la figure 3.30 montre le cycle de la
vapeur avec une vanne d'étranglement entre le générateur de vapeur et le détendeur (chemin
G-G ').
Fig. 3.30 Cycle de Bourque dans le plan température-entropie (Robert Bourque [150]).
Kang Ma et al. [151], ont étudié les effets des lois de transfert de chaleur sur
l’optimisation du moteur pour un travail maximal. La loi généralisée de transfert radiatif de
80
chaleur [𝑞 ∝ ∆(𝑇 𝑛 )] est introduit dans un modèle de moteur à combustion externe ayant un
piston mobile. Des exemples numériques d’optimisation avec les lois de transfert de chaleur
linéaires phénoménologiques (n = -1), de Newton (n = 1), en carré (n = 2), cubiques (n = 3) et
radiatives (n = 4) sont présentés et les résultats obtenus sont comparés. Ces résultats montrent
que le travail et le rendement du moteur en mode de fonctionnement optimal cyclique et en
mode de fonctionnement optimal semi-cyclique diminuent avec l'augmentation de la
conductivité thermique, le travail, le taux de compression et le rendement en mode
d’opération semi-cyclique sont plus grandes que ceux en mode d’opération cyclique. Bien que
toutes les courbes de volume en fonction du temps du mode de fonctionnement optimal
cyclique et des arcs d’Euler-Lagrange (EL) en mode de fonctionnement optimal semi-
cyclique sont à peu près sinusoïdales et se composent de trois étapes avec les cinq lois de
transfert de chaleur, les courbes avec diverses lois de transfert de chaleur sont différentes.
La modélisation d’un moteur Stirling est un problème complexe auquel de nombreux
chercheurs se sont attelés. Elle met en jeu la modélisation des flux massiques, du flux
thermique et des évolutions thermodynamiques.
La complexité de l’écoulement du fluide de travail et des échanges de chaleur associés
dans l’absorbeur et le régénérateur rend la modélisation et la prédiction des performances de
ce type de moteur, difficiles à réaliser. Dans le régénérateur, l’écoulement est turbulent,
tridimensionnel et surtout oscillant et dans le réchauffeur (l’absorbeur); il est équivalent à un
écoulement instationnaire dans une conduite cylindrique de petites dimensions soumise à un
flux thermique non uniforme [152].
A cet effet, il existe trois niveaux d’analyse thermodynamique des moteurs Stirling:
les analyses du premier ordre (dites approximatives), du deuxième ordre (découplées), et du
troisième ordre (couplée) [153]. Celles du premier ordre permettent une évaluation
approximative des performances du moteur. Les analyses du deuxième ordre permettent une
première optimisation. Les méthodes de troisième ordre servent à réaliser des simulations
fines de la distribution de masse, de pression et de température dans le moteur [154]. A ce
niveau, l’utilisation des codes de calcul est indispensable.
G. Schmidt [155] en 1871, a développé la première modélisation zéro-dimensionnelle
d’un moteur Stirling. Il a proposé une analyse, connue sous le nom d'analyse idéale, en se
basant sur les modèles développés ultérieurement. L'analyse idéale repose sur la division du
moteur en ses cinq régions élémentaires. On suppose qu'on peut définir une température
caractéristique de chacun des cinq espaces élémentaires.
Les principales hypothèses du cycle Schmidt sont:
81
Cependant, la surface d'échange dans les moteurs réels aura tendance à être plutôt
adiabatique. La première analyse de la non-isothermie a été proposée par Finkelstein [156]. Il
a développé une analyse dans laquelle il a supposé que les deux espaces de travail sont
parfaitement isolés. Il a introduit dans ce modèle le concept de température conditionnelle qui
dépend du sens du flux. Le moteur Stirling est considéré comme un ensemble de cinq
compartiments et les pertes thermiques et mécaniques sont négligées. Martini [157] a effectué
une simulation numérique sur un moteur Stirling, divisé en cinq espaces de travail en
supposant que les températures du fluide de travail dans l’espace de compression et celui de
détente sont uniformes et les températures du réchauffeur et du refroidisseur sont constantes,
tandis que la température du régénérateur a un profil linéaire.
Les moteurs Stirling réels sont confrontés aux différents types de déperditions, telles
que les pertes thermiques et mécaniques. Les études menées par Ureili et Berchowitz [158]
montrent nettement que le modèle adiabatique à transfert de chaleur imparfait ou modèle dit
‘à écoulement quasi permanent’ est plus performant que le modèle dit ‘à écoulement
permanent’.
Le prolongement de ce travail de modélisation adiabatique idéale consistera à
développer des modèles instationnaires. Parlak et al. [159] développent un modèle
adiabatique qui diffère des modèles idéaux présentés précédemment. Cette modélisation, en se
basant sur celle de Urieli et Berchowitz [158], considère les volumes de compression et de
détente adiabatique et les échangeurs imparfaits. De ce fait, si la température de chaque zone
est uniforme, chacune évolue au cours du cycle.
M. Abbas et al. [160] ont présenté une modélisation dynamique d’un moteur Stirling
de type gamma, basée une approche quasistationnaire. Le programme de calcul élaboré en
82
Fig. 3.31 Schema du micro moteur thermique (H. Bardaweel et al. [161]).
83
Les résultats montrent que le cycle de travail est à peu près de forme rectangulaire et
se compose de deux processus de température constante et deux processus de volume
constant. Les variations maximales de pression et de volume mesurées à l'intérieur de la
cavité du moteur sont respectivement 45 Pa et 0.55 mm3. Les résultats montrent que pour un
apport de chaleur de 1 mJ, le moteur fonctionne sur une différence de température très faible.
Cette faible différence de température est principalement attribuable au stockage thermique
important de la structure du moteur, des membranes et des mèches. Le rendement exergétique
(appelé aussi rendement selon la seconde loi ou rendement rationnel) mesuré du micromoteur
thermique est de 16%.
84
Pour ce qui est des moteurs fonctionnant à l’hydrogène pur, le retour des flammes a
été un obstacle particulièrement tenace au développement de ce type de moteurs. La plupart,
sinon la totalité, de la littérature mentionne que le retour des flammes se produit seulement
quand le mélange combustible est présent dans la tubulure d'admission (formation externe du
mélange gazeux), c’est à dire dans les moteurs à injection indirecte du mélange combustible.
Plusieurs méthodes de formation du mélange dans les moteurs à hydrogène ont été
testées, afin d’éviter le phénomène du retour des flammes. Dans ce contexte, on peut citer
l'injection d'eau, l'introduction tardive de l'hydrogène et l'injection directe, qui sont toutes
principalement destinées à retarder ou à empêcher le retour des flammes.
Quant au phénomène du cliquetis dans les moteurs alimentés par hydrogène, il a été
mal expliqué. Il y’avait des revendications totalement contradictoires dans la littérature.
Néanmoins, on peut mentionner que la réduction de la vitesse de montée en pression peut être
plus efficace pour contrôler les cliquetis que la limitation de la durée de combustion
L’examen de la littérature lié aux études expérimentales sur ce type de moteurs montre
que le pré-allumage est un facteur limitant pour les taux de compression, le moment
d’allumage et les richesses du mélange.
La littérature sur la simulation des moteurs à hydrogène est assez rare. Néanmoins,
certains auteurs ont utilisé le modèle à deux zones quasi-dimensionnel, les modèles zero-et
multi-dimensionnels et le code de simulation du moteur Kiva-3V, afin de simuler de tels
moteurs. La caractéristique importante de Kiva-3V est que c'est un code source ouverte, qui
peut être modifié selon les demandes de l’utilisateur. L’approche zéro-dimensionnelle ne
permet pas d’atteindre les résultats obtenus par des modèles dimensionnels ou par
l’expérimentation, mais présente un intérêt majeur pour les études paramétriques. La
turbulence est un phénomène tridimensionnel qui est bien simulé par les simulations
multidimensionnelles, qui permettent de déterminer les coefficients de transfert de chaleur
sans faire appel à des expérimentations en laboratoire.
Les études sur l’utilisation du gaz naturel à l’état pur comme carburant de substitution
en remplacement des énergies fossiles, se sont beaucoup développées ces dernières années. le
gaz naturel est un mélange gazeux qui possède un fort potentiel dans l’aspect dépollution,
mais aussi dans l’aspect performance quant à son utilisation comme carburant.
85
Par contre, le mélange du gaz naturel avec l’hydrogène reste un domaine très peu
étudié. On a commencé à s’intéresser à la combustion de ce genre de mélange que depuis peu
de temps. Quelques articles seulement ont été publiés sur ce sujet.
L’ajout d’hydrogène au gaz naturel donne la possibilité au moteur de fonctionner avec
un mélange pauvre sous des conditions de charge partielle.
L’ajout de l’hydrogène au gaz naturel améliore: le rendement du moteur et sa pression
moyenne efficace. Il réduit les durées de combustion, l’émission de CO2 (en raison de la
diminution du taux de C/H) et la consommation de carburant. Les émissions de NOx sont, au
contraire, plus élevées que celles du gaz naturel pur en raison de l’augmentation de la
température à l’intérieur du cylindre, pour un rapport d'équivalence donnée. Quoi qu'il en soit,
l'utilisation des mélanges pauvres ou l'EGR réduit vraiment les émissions de NOx et entraîne
une augmentation supplémentaire du rendement du moteur. Néanmoins, la possibilité de
réduire le NOx sans augmenter les émissions des HC par l'ajout d'hydrogène pour prolonger la
limite du mélange pauvre est confirmée. Une bonne combustion du mélange Hydrogene /Gaz
naturel contribue à maintenir les émissions des HC faibles. En plus, la vitesse de flamme
laminaire du gaz naturel est inférieure à celui de l'essence et l'addition de l'hydrogène, qui
présente une vitesse de flamme laminaire environ huit fois supérieure à celle du gaz naturel,
améliore de manière significative cette propriété principale de combustion.
Le gaz naturel comprimé enrichi par hydrogène (GNCH) est vraisemblablement une
alternative potentielle aux carburants issus des hydrocarbures dans les moteurs à combustion
interne. Des recherches expérimentales sur l’utilisation de ce genre de carburant ont été
conduites pendant plus de dix années et sont toujours en cours.
Bien que l'addition d'hydrogène augmente les émissions des NOx, elle s'est avérée
bénéfique simultanément pour l'économie du combustible, la pression moyenne effective, le
rendement thermique indiqué et pour les émissions des HC.
Plusieurs nouveaux moteurs dont le moteur de Strinling peuvent permettre de réduire
les problèmes associés à la combustion classique. Cependant, ils restent sujettes à la
variabilité cycle à cycle ou variabilité cyclique de combustion (VCC) qui limite leurs
performances. Ces variabilités, qui sont d'une grande importance dans les moteurs
fonctionnant avec un mélange pauvre, peuvent en effet être abaissées par l'ajout d'hydrogène.
L’ajout d’hydrogène peut aussi faire diminuer le couple maximal et le coefficient de variation
(CoV) de la pression maximale (Pmax) et de la pression effective moyenne indiquée (PEMI).
Peu d’études ont été réalisées sur l’effet d’ajout d’hydrogène au méthane. Cette
addition conduit à une augmentation des vitesses de propagation de flammes et une
86
3.9 Références
[1] Bardon M.F. and Haycock R.G. The hydrogen research of R.O. King. 14th World
Hydrogen Energy Conference, invited paper, Montreal, Canada, 2002.
87
[2] MacCarley C.A. A study of factors influencing thermally induced backfiring in hydrogen
fuelled engines, and methods for backfire control. 16th IECEC conference, Atlanta, USA,
1981.
[3] Das L.M. Near-term introduction of hydrogen engines for automotive and agricultural
application. International Journal Hydrogen Energy, 27:479–487, 2002.
[4] Lucas G.G. and Morris L.E. The backfire problem of the hydrogen engine. Symposium
organized by the university’s internal combustion engine group, King’s College, London, UK,
1980.
[5] Das L.M. Hydrogen-oxygen reaction mechanism and its implication to hydrogen engine
combustion. Int. J. Hydrogen Energy, 21:703–715, 1996.
[6] Berckmuller M. et al. Potentials of a charged SI-hydrogen engine. SAE, paper N°2003-01-
3210, 2003.
[7] Stockhausen W.F. et al. Ford P2000 hydrogen engine design and vehicle development
program. SAE, paper N°2002-01-0240, 2002.
[8] R.K. Rajput. Internal Combustion Engines. Laxmi Publications (P) LTD, 2005, 1027
pages, p. 59.
[9] Project Coordinator Motor Vehicles and Road Transport, TUV Rheinland e.V. for the
Federal Ministry for Research and Technology. Alternative energy sources for road transport
- hydrogen drive test. Technical report, TUV Rheinland, 1990.
[10] Kondo T., Iio S., and Hiruma M. A study on the mechanism of backfire in external
mixture formation hydrogen engines –about backfire occurred by the cause of the spark plug -
SAE, paper N°971704, 1997.
[11] Swain M.R., Schade G.J., and Swain M.N. Design and testing of dedicated hydrogen
fueled engine. SAE, paper N°961077, 1996.
[12] Koyanagi K., Hiruma M. and Furuhama S. Study on mechanism of backfire in hydrogen
engines. SAE, paper N°942035, 1994.
[13] Lee J.T., Kim Y.Y., Lee C.W., and Caton J.A. An investigation of a cause of backfire
and its control due to crevice volumes in a hydrogen fueled engine. ASME Spring Technical
Conference, paper N°2000-ICE-284, San Antonio, USA, 2000.
[14] Tang X. et al. Ford P2000 hydrogen engine dynamometer development. SAE, paper
N°2002-01-0242, 2002.
88
[15] Akagawa H. et al. Development of hydrogen injection clean engine. 15th World
Hydrogen Energy Conference, paper nr 28J-05, Yokohama, Japan, July 2004.
[16] Lee S.J., Yi H.S., and Kim E.S. Combustion characteristics of intake port injection type
hydrogen fueled engine. International Journal Hydrogen Energy, 20:317–322, 1995.
[18] Jorach R.W. Low nitrogen oxide emission combined with high power density by using
the hydrogen combustion method with advanced internal mixture formation. MTZ
Worldwide, 58-4:5–8, 1997.
[19] Jing-Ding L., Ying-Qing L., and Tian-Shen D. Improvement on the combustion of a
hydrogen fueled engine. Int. J. Hydrogen Energy, 11:661–668, 1986.
[20] Natkin R.J. et al. Hydrogen IC engine boosting performance and NOx study. SAE, paper
N°2003-01-0631, 2003.
[21] Binder K. and Withalm G. Mixture formation and combustion in hydrogen engine using
hydrogen storage technology. Int. J. Hydrogen Energy, 7:651–659, 1982.
[22] Li H. and Karim G.A. Knock in spark ignition hydrogen engines. International Journal
Hydrogen Energy, 29:859–865, 2004.
[23] Karim G.A. Hydrogen as a spark ignition engine fuel. Int. J. Hydrogen Energy, 28:569–
577, 2003.
[24] Heffel J.W., Johnson D.C., and Shelby C. Hydrogen powered Shelby Cobra: vehicle
conversion. SAE, paper N°2001-01-2530, 2001.
[25] Smith J.R., Aceves S., and Van Blarigan P. Series hybrid vehicle and optimized
hydrogen engine design. SAE, paper N°951955, 1995.
[26] Li H. and Karim G.A. Hydrogen fuelled spark ignition engines: predictive and
experimental performance. ASME Spring Technical Conference, paper nr ICES2003-548,
Salzburg, Austria, 2003.
[27] Olavson L.G., Baker N.R., Lynch F.E., and Meija L.C. Hydrogen fuel for underground
mining machinery. SAE, paper N°840233, 1984.
[28] Davidson D., Fairlie M., and Stuart A.E. Development of a hydrogen-fuelled farm
tractor. Int. J. Hydrogen Energy, 11:39–42, 1986.
89
[29] Heffel J.W., McClanahan M.N., and Norbeck J.M. Electronic fuel injection for hydrogen
fueled internal combustion engines. SAE, paper N°981924, 1998.
[30] Meier F. et al. Cycle-resolved hydrogen flame speed measurement with high-speed
Schlieren technique in a hydrogen direct injection SI engine. SAE, paper N°942036, 1994.
[31] Furuhama S. Problems of forecasting the future of advanced engines and engine
characteristics of the hydrogen injection with LH2 tank and pump. Journal of Engineering for
Gas Turbines and Power, 119:227–242, 1997.
[32] Guo L.S., Lu H.B., and Li J.D. A hydrogen injection system with solenoid valves for a
four-cylinder hydrogen-fuelled engine. Int. J. Hydrogen Energy, 24:377–382, 1999.
[33] Kim J.M., Kim Y.T., Lee J.T., and Lee S.Y. Performance characteristics of hydrogen
fueled engine with the direct injection and spark ignition system. SAE, paper N°952498,
1995.
[34] Yi H.S., Min K., and Kim E.S. The optimised mixture formation for hydrogen fuelled
engines. International Journal Hydrogen Energy, 25:685–690, 2000.
[35] Peschka W. and Escher W.J.D. Germany’s contribution to the demonstrated technical
feasibility of the liquid-hydrogen fueled passenger automobile. SAE, paper N°931812,
1993.
[36] Rottengruber H. et al. A high-efficient combustion concept for direct injection hydrogen
internal combustion engine. 15th World Hydrogen Energy Conference, paper N°28J-01,
Yokohama, Japan, 2004.
[37] Heffel J.W. NOx emission reduction in a hydrogen fuelled internal combustion engine at
3000 rpm using exhaust gas recirculation. International Journal Hydrogen Energy, 28:1285 -
1292, 2003.
[38] Van Blarigan P. Development of a hydrogen fueled internal combustion engine designed
for single speed/power operation. SAE, paper N°961690, 1996.
[39] Aceves S.M. and Smith J.R. Hybrid and conventional hydrogen engine vehicles that
meet EZEV emissions. SAE, paper N° 970290, 1997.
[40] Gerbig F. et al. Potentials of the hydrogen combustion engine with innovative hydrogen-
specific combustion process. Fisita World Automotive Congress, paper N°
F2004V113, Barcelona, Spain, 2004.
90
[41] DeLuchi M.A. Hydrogen vehicles: an evaluation of fuel storage, performance, safety,
environmental impacts, and cost. Int. J. Hydrogen Energy, 14:81–130, 1989.
[42] Strebig K.C. and Waytulonis R.W. The bureau of mines hydrogen powered mine vehicle.
SAE, paper N°871678, 1987.
[43] Fagelson J.J., McLean W.J., and de Boer P.C.T. Performance and NOx emissions of
spark-ignited combustion engines using alternative fuels - quasi one-dimensional modeling. I.
hydrogen fueled engines. Combustion Science and Technology, 18:47–57, 1978.
[44] Prabhu-Kumar G.P., Nagalingam B., and Gopalakrishnan K.V. Theoretical studies of a
spark-ignited supercharged hydrogen engine. Int. J. Hydrogen Energy, 10:389–397, 1985.
[45] Keck J.C. Turbulent flame structure and speed in spark-ignition engines. 19th Symp.
(Int.) on Combustion, pages 1451–1466, 1982.
[47] Fontana G., Galloni E., Jannelli E., and Minutillo M. Numerical modeling of a spark
ignition engine using premixed lean gasoline-hydrogen-air mixtures. 14th World Hydrogen
Energy Conference, Montreal, Canada, 2002.
[48] Takats M. and al. Hydrogen fueled reciprocating engine as an automotive prime mover,
Fisita World Automotive Congress, paper nr F98T/P693, Paris, France, 1998.
[49] Polasek M., Macek J., Takats M., and Vıtek O. Application of advanced simulation
methods and their combination with experiments to modeling of hydrogen fueled engine
emission potentials. SAE, paper N°2002-01-0373, 2002.
[50] Ma J., Su Y., Zhou Y., and Zhang Z. Simulation and prediction on the performance of a
vehicles hydrogen engine. Int. J. Hydrogen Energy, 28:77–83, 2003.
[51] Goto, Y. “Mixture Formation and Ignition in a Direct Injection Natural Gas Engine”,
JSME International Journal, Series B Vol. 42, N° 2, 1999.
[52] Rubas, P.J., Raul, M.A., Martin, G.C., Coverdill, R.E., Lucht, R.P., Peters, J.E.,
DelVecchio, K.A. “Methane Jet Penetration in a Direct-Injection Natural Gas Engine” SAE,
paper N° 980143, 1998.
91
[53] Hiltner, J., Samimy, M. “The Impact of Injection Timing on In-Cylinder Fuel
Distribution in a Natural Gas Powered Engine”, SAE, paper N° 971708, 1997.
[54] Aesoy, V., Valland, H. “The Influence of Natural Gas Composition on Ignition in a
Direct Injection Gas Engine Using Hot Surface Assisted Compression Ignition”, SAE, paper
N°961934, 1996.
[55] Naber, J. D., Siebers, D. L., Di Julio, S. S., Westbrook, C. K. “Effects of Natural Gas
Composition on Ignition Delay under Diesel Conditions”, Combustion and Flame, Vol. 99 pp
192-200, 1994.
[56] P. Einewall, P. Tunestål, B. Johansson, “Lean Burn Natural Gas vs. Stoichiometric
Operation with EGR and a Three Way Catalyst”, SAE 2005-01-0250
[57] C. Nellen, K. Boulouchos, “Natural Gas Engines for Cogeneration: Highest Efficiency
and Near-ZeroEmissions through Turbocharging, EGR and 3-Way Catalytic Converter”, SAE
paper N°2000-01-2825.
[58] T. Reppert, J. Chiu, “Heavy-Duty Waste Hauler with Chemically Correct Natural Gas
Engine Diluted with EGR and Using a Three-Way Catalyst”, Subcontract report NREL/SR-
540-38222, September 2005.
[61] Zuohua Huang, Jinhua Wang, Bing Liu, Ke Zeng, Jinrong Yu, and Deming Jiang,
Combustion Characteristics of a Direct-Injection Engine Fueled with Natural Gas-Hydrogen
Blends under Various Injection Timings, Energy & Fuels 2006, 20, 1498-1504.
[66] Choudhuri A.R. and Gollahalli S.R., ``Intermediate radical concentrations in hydrogen-
natural gas blended fuel jet flames``, International Journal of Hydrogen Energy, 29, 12, 1293-
1302, 2004.
[67] Peters, N. & Rogg, B. (Eds.). Reduced Kinetic Mechanisms for Application in
Combustion Systems, Lecture Notes in Physics, New Series, M 15, Springer Verlag. 1993.
[68] Dagaut P. et Nicolle A., ``Experimental and detailed kinetic modeling study of hydrogen-
enriched natural gas blend oxidation over extended temperature and equivalence ratio
ranges``, Proceedings of the Combustion Institute, 30, 2, 2631-2638, 2005.
[69] Dagaut P. and Dayma G., ``Hydrogen-enriched natural gas blend oxidation under high-
pressure conditions: Experimental and detailed chemical kinetic modeling``, International
Journal of Hydrogen Energy, 31, 4, 505-515, 2006.
[70] Dagaut P., ``On the kinetics of hydrocarbons oxidation from natural gas to kerosene and
diesel fuel``, Physical Chemistry Physics, 4, 11, 2079-2094, 2002.
[71] Huang Z., Zhang Y., Zeng K., Liu B., Wang Q. Et Jaing D., ’’Measurements of laminar
burning velocities for natural gas-hydrogen-air mixtures’’, Combustion and Flame, 146, 1-2,
302-311, 2006b.
[72] Swain M.R., Yusuf M.J., Dulger Z., Swain M.N. The effects of hydrogen addition on
natural gas engine operation. SAE, paper N°932775, 1993.
[74] Nagalim B., Duebel F., Schmillen K., Performance study using natural gas, hydrogen-
supplemented natural gas and hydrogen in AVL research engine. International Journal
Hydrogen Energy 1983; 9(8): 715–20.
[75] Das LM. Utilization of hydrogen-CNG blend in internal combustion engine, 11th World
Hydrogen Energy Conference, Stuttgart, Germany, June 23–28, 1996. p. 1513–35.
93
[76] Larsen J.F., Wallace J.S., Comparison of emissions and efficiency of a turbocharged
lean-burn natural gas and hythane-fueled engine. J Eng Gas Turbines Power 1997; 119:218–
26.
[77] Sierens R, Rosseel E. Variable composition hydrogen/natural gas mixtures for increased
engine efficiency and decreased emissions. Journal Eng Gas Turbines Power 2000; 122:135–
40.
[78] Das L.M., Gulati R., Gupta P.K. A comparative evaluation of the performance
characteristics of a spark ignition engine using hydrogen and compressed natural gas as
alternative fuels. International Journal Hydrogen Energy 2000; 25:783–93.
[79] Ma F., Ding S., Wang Y., Wang M., Jiang L., Naeve N., Zhao S. (2009a). Performance
and Emission Characteristics of a Spark-Ignition (SI) Hydrogen-Enriched Compressed
Natural Gas (HCNG) Engine Under Various Operating Conditions Including Idle Conditions.
Energy & Fuels, 23, 2113-3118.
[80] Ma, F., Ding, S., Wang, Y., Wang, Y., Wang, J., Zhao, S. (2008a). Study on combustion
behaviors and cycle-by-cycle variations in a turbocharged lean burn natural gas S.I. engine
with hydrogen enrichment. International Journal of Hydrogen Energy, 33, 7245-7255.
[81] Ma, F., Liu, H., Wang, Y., Wang, J., Ding, S., Zhao, S. (2008b). A Quasi-Dimensional
Combustion Model for SI Engines Fuelled by Hydrogen Enriched Compressed Natural Gas.
SAE Paper 2008-01-1633.
[82] Ma, F., Liu, H., Wang, Y., Li, Y., Wang, J., Zhao, S. (2008c). Combustion and emission
characteristics of a port-injection HCNG engine under various ignition timings. International
Journal of Hydrogen Energy, 33, 816-822.
[83] Ma, F., Wang, Y. (2008d). Study on the extension of lean operation limit through
hydrogen enrichment in a natural gas spark-ignition engine. International Journal of Hydrogen
Energy, 33, 1416-1424.
[84] Ma, F., Wang, Y., Liu, H., Li, Y., Wang, J., Ding, S. (2008e). Effects of hydrogen
addition on cycle-by-cycle variations in a lean burn natural gas spark-ignition engine.
International Journal of Hydrogen Energy, 33, 823-831.
[85] Burke, Andy; McCaffrey, Zach; Miller, Marshall; Collier, Kirk; Mulliga, Neal:
Hydrogen Bus Technology Validation Program, Institute of Transportation Studies 2005.
[86] Ma, F., Wang, Y., Wang, J., Zhao, S., Yin, Y., Cheng, W., Zhou, M. (2008g).
Development and validation of a quasi-dimensional combustion model for SI engine fuelled
by HCNG with variable hydrogen fractions. SAE, paper N°2008-01-1580.
94
[87] Ma, F., Wang, J., Wang, Y., Wang, Y., Li, Y., Liu, H., Ding, S. (2008h). Influence of
different volume percent hydrogen/natural gas mixtures on idle performance of a CNG
engine. Energy & Fuels, 22, 1880-1887.
[88] Ma, F., Wang, J., Wang, Y., Wang, Y., Zhong, Z., Ding, S., Zhao, S. (2008i). An
investigation of optimum control of a spark ignition engine fueled by NG and hydrogen
mixtures, International Journal of Hydrogen Energy, 33, 7592-7606.
[89] Ma, F., Wang, Y., Wang, J., Zhao, S., Yin, Y., Cheng, W., Zhou, M. (2008j).
Development and Validation of an On-line Hydrogen-Natural Gas Mixing System for Internal
Combustion Engine Testing. SAE, paper N°2008-02-1508.
[94] Yusuf MJ. In cylinder flame front growth rate measurement of methane and hydrogen
enriched methane fuel in a spark ignited internal combustion engine. Master thesis, University
of Miami, 1990.
[95] Yusuf MJ. Lean Burn natural gas fueled engines: engine modification versus hydrogen
blending. PhD thesis, University of Miami, 1993.
[96] Hoekstra RL, Collier K, Mulligan N, Demonstration of hydrogen mixed gas vehicles,
10th World hydrogen Energy Conference, Cocoa Beach, USA, June 20–24, 1994.
[97] Hoekstra RL, Collier K, Mulligan N, Chew L. Experimental study of clean burning
vehicle fuel. International Journal Hydrogen Energy 1995; 20:737–45.
[98] Liu Z, Karim GA. Knock characteristics of dual-fuel engines fueled with hydrogen fuel.
International Journal Hydrogen Energy 1995; 20: 919–24.
95
[99] Karim GA. Hydrogen as an additive to methane for engine applications. 11th World
Hydrogen Energy Conference, Stuttgart, Germany, June 23–28, 1996. p. 1921–526.
[101] Blarigan PV, Keller JO. A hydrogen fueled internal combustion engine designed for
single speed/power operation. International Journal Hydrogen Energy 1998; 23: 603–9.
[102] Wong YK, Karim GA. A kinetic examination of the effects of the presence of some
gaseous fuels and preignition reaction products with hydrogen in engines. International
Journal Hydrogen Energy 1999; 24:473–8.
[103] Shrestha SO, Karim GA. hydrogen as an additive to methane for spark ignition engine
applications. International Journal Hydrogen Energy 1999; 24:577–86.
[105] Bauer CG, Forest TW. Effect of hydrogen addition on performance of methane-fueled
vehicles. Part I: effect on S.I. engine performance. International Journal Hydrogen Energy
2001; 26:55–70.
[106] Bauer CG, Forest TW. Effect of hydrogen addition on performance of methane-fueled
vehicles. Part II: driving cycle simulation. International Journal Hydrogen Energy 2001;
26:71–90.
[107] Y.K. Wong, G.A. Karim. An analytical examination of the effects of hydrogen addition
on cyclic variations in homogeneously charged compression-ignition engines. International
Journal of Hydrogen Energy 25 (2000) 1217-1224.
[108] Nagalingam B, Duebel F, Schmillen K. Performance study using natural gas, hydrogen
supplemented natural gas and hydrogen in AVL research engine. International Journal of
Hydrogen Energy 1983; 8(9):715-20.
[109] Kyaw ZH. Chemical control of S.I. engine combustion. Ph.D. thesis, Department of
Mechanical and Manufacturing Engineering, The University of Melbourne, 1992. p. 38-45.
[110] Cattelan A, Wallace J. Exhaust emission and energy consumption effects from
hydrogen supplementation of natural gas. SAE, paper N°952497, 1995.
[111] Collier K, Hoekstra RL, Mulligan N, Jones C, Hahn D. Untreated exhaust emissions of
a hydrogen-enriched CNG production engine conversion. SAE, paper N°960858, 1996.
96
[112] Hoekstra RL, Van Blarigan P, Mulligan N. NOx emissions and efficiency of hydrogen,
natural gas, and hydrogen = natural gas blended fuels. SAE, paper N°961103, 1996.
[113] N.B. Anikin, S. Gersen, A.V. Mokhov, H.B. Levinsky. Autoignition of Methane/
Hydrogen Mixtures in a Rapid Compression Machine. Laboratory for Fuel and Combustion
Science, University of Groningen, The Netherlands Gas Engineering and Technology,
Groningen, The Netherlands. Proceedings of the European Combustion Meeting 2007.
[115] E.L. Petersen, D.M. Kalitan, S. Simmons, G. Bourgue, H.J. Curran and J.M. Simmie,
Proceedings of the Combustion Institute 31 (2007) doi:10.1016/j.proci.2006. 08.34.
[116] Park J., In Keel S., Han Yun J. and Kwon Kim T., ``Effects of addition of electrolysis
products in methane-air diffusion flames``, International Journal of Hydrogen Energy, 32, 16,
4059-4070, 2007.
[117] Smith, P., Golden, D., Frenklach, M., Moriarty, N., Eiteneer, B., Goldenberg, M.,
Bowman, C., Hanson, R., Song, S., Gardiner, W., Lissianski, V., and Qin, Z. (1999). http://
www.me.berkeley.edu/gri_ mech/.
[118] Takeno T. and Nishioka M., ``Species conservation and emission indices for flames
described by similary solutions``, Combustion and flame, 92, 4, 465-468, 1993.
[119] Scholte T.G. and Vaags P.B., ``Burning velocities of mixture of hydrogen, carbon
monoxide and methane with air’’, Combustion and Flame, 3, 511-524, 1959.
[120] Miller D.R., Evers R.L. and Skinner G.B., ``Effects of various inhibitors on hydrogen-
air flame speeds’’, Combustion and Flame, 7, 1, 137-142, 1963.
[121] Liu Y., Lenz B. and Leuckel W., ``Investigation on the laminar and Turbulent Burning
Velocities of Premixed Lean and Rich Flames of Methan-Hydrogen-Air Mixtures’’, Progress
in Astronautics and Aeronautics, 131, (Dyn. Deflagrations React. Syst.: Flames), 1991, pp.
259-274.
[122] Tanoue K., Goto S., Shimada F. and Hamatake T., ’’Effects of hydrogen addition on
stretched premixed laminar methane flames (1st report, effects on laminar burning velocity)’’,
Nippon Kikai Gakkai Ronbunshu, B Hen/Transactions of the Japan Society of Mechanical
Engineers, Part B, 69, 677, 162-168, 2003.
[123] Yu G., Law C.K. and Wu C.K., ’’Laminar flame speeds of hydrocarbon + air mixtures
with hydrogen addition’’, Combustion and flame, 63, 3, 339-347, 1986.
97
[124] Law C.K., Kwon O.C., ’’Effects of hydrocarbon substitution on atmospheric hydrogen-
air flame propagation’’, International Journal of Hydrogen Energy, 29, 8, 867-879, 2004.
[125] Ren J.Y., Qin W., Egolfopoulos F.N. and Tsotsis T.T., ’’Strain-rate effects on
hydrogen-enhanced lean premixed combustion’’, Combustion and Flame, 124, 4, 717-720,
2001.
[126] Tse S.D., Zhu D.L., Law C.K., ’’Morphology and burning rates of expanding spherical
flames in H2/O2/inert mixtures up to 60 atmospheres’’, Proc Combust Inst, 28 (2000), pp.
1793–1800.
[127] Kee R.J., Rupley F.M. and Miller J.A., ’’CHEMKIN II: A FORTRAN chemical
kinetics package for the analysis of gas phase chemical kinetics’’, Sandia National
Laboratories Report No SAND89-8009B, SAND89-8009B, 1991.
[128] Z. Qin, V. Lissianski, H. Yang, W.C. Gardiner Jr, S.G. Davis, H. Wang, Combustion
chemistry of propane a case study of detailed reaction mechanism optimization Proc Combust
Inst, 28 (2000), pp. 1663–1669.
[129] Kee RJ, Dixon-Lewis G, Warnatz J, Coltrin ME, Miller JA. A FORTRAN computer
code package for the evaluation of gas-phase, multi-component transport properties.
Albuquerque, NM: Sandia National Laboratories Report SAND86-8246; 1992.
[130] Halter F., Chauveau C., Djebaili-Chaumeix N. and Gokalp I., ’’Characterization of the
effects of pressure and hydrogen concentration on laminar burning velocities of methane-
hydrogen-air mixtures’’, Proceedings of the Combustion Institute, 30, 1, 201-208, 2005.
[131] Mandilas C., Ormsby M.P., Sheppard C.G.W. and Woolley R., ’’Effects of hydrogen
addition on laminar and turbulent premixed methane and iso-octane-air flames’’, Proceedings
of the Combustion Institute, 31, 1, 1443-1450, 2007.
[132] Lifshitz A., Scheller K., Burcat A. and Skinner G.B., ’’Shock-tube investigation of
ignition in methane-oxygen-argon mixtures’’, Combustion and Flame, 16, 3, 311-321, 1971.
[133] Huang J., Bushe W.K., Hill P.G. and Munshi S.R., ``Experimental and kinetic study of
shock initiated ignition in homogeneous methane-hydrogen-air mixtures at engine-relevant
conditions``. International Journal of Chemical Kinetics, 38, 4, 221-233, 2006a.
[134] Fotache C.G., Kreutz T.G. and Law C.K., ``Ignition of hydrogen-enriched methane by
heated air``, Combustion and Flame, 110, 4, 429-440, 1997.
98
[135] Bowman C.T., Hanson R.K., Davidson D.F., Gardiner W.C., Lissianki V., Smith G.P.,
Golden D.M., Frenklach M. and Goldenberg M., ``GRI 2.11``. http://www.me.
berkley.edu/gri_mech/, 1995.
[137] Chaumeix N., Pichon S., Lafosse F. and Paillard C.E., ``Role of chemical kinetics on
the detonation properties of hydrogen/natural gas/air mixtures``, International Journal of
Hydrogen Energy, 32, 13, 2216-2226, 2007.
[138] Marinov N.M., Pitz, W.J., Westbrook C.K., Vencitore C.K., Castaldi M.J., Senkan
S.M., ``Aromatic and polycyclic Aromatic Hydrocarbon Formation I a Laminar Premixed n-
Butane flame Combustion and flame``, 114 192-213, 1998.
[140] Konnov A.A., ``Detailed reaction mechanism for small hydrocarbons combustion``,
Release 0.5, http://homepages.vub.ac.be/~akonnov/, 2000.
[141] Gersen S., Anikin N.B., Mokhov A.V. and Levinsky H.B., ``Ignition properties of
methane/hydrogen mixtures in a rapid compression machine``, International Journal of
Hydrogen Energy, 33, 7, 1957-1964, 2008.
[142] Petersen E.L, Kalitan D.M., Simmons S., Bourque G., Curran H.J. and Simmie
J.M., ``Methane/propane oxydation at high pressures : Experimental and detailed chemical
kinetic modeling``, Proceedings of the Combustion Institute, 31, 1, 447-454, 2007.
[143] Bagdanavicius A, Bowen PJ, Syred N, Crayford AP, Turbulent Flame Structure of
Methane-Hydrogen Mixtures at Elevated Temperature and Pressure, Combustion Science and
Technology , 185 (2) (2013) 350-361 ISSN 0010-2202.
[144] Gauducheau J.L., Denet B., Searby G., ``A Numerical Study of Lean CH4/H2/Air
Premixed Flames at High Pressure``, Combustion Science and Technology, 137, 1-6, 81-99,
1998.
[145] Tseng C.-J., ``Effects of the Hydrogen Addition on Methane Combustion in a Porous
Medium Burner``, International Journal of Hydrogen Energy, 27, 6, 699-707, 2002.
99
[146] Uykur C., Henshaw P.F., Ting D.S.K. and Barron R.M., ``Effects of Addition of
Electrolysi Products on Methane/Air Premixed Laminar Combustion``, International Journal
of Hydrogen Energy, 26, 3, 265-273, 2001.
[147] Egolfopoulos, F. N., Zhang, H., and Zhang, Z., ``Wall Effects on the Propagation and
Extinction of Steady, Strained, Laminar Premixed Flames``, Combustion and Flame, 109, 1–
2, 237–252, 1997.
[148] İsmet Çelikten, Can Çinar, E. İshak Yildiz, The effects of hydrogen enriched LPG on
engine performance and exhaust emissions in spark ignition engine, Teknoloji, Volume 10,
(2007), Issue 3, 215-222.
[149] Choi, G.H., Chung, Y.J., Han, S.B., “Performance and Emissions of a Hydrogen
Enriched LPG Internal Combustion Engine 1400 rpm”, International Journal of Hydrogen
Energy, 30(1), 77-82, 2005.
[150] Robert Bourque. A Compact External Combustion Engine with High Part-Load
Efficiency. ASME 2010 4th International Conference on Energy Sustainability (ES2010)
Volume 1 May 17–22, 2010, pp. 735-746.
[151] Kang Ma, Lingen Chen, Fengrui Sun. Optimizations of a model external combustion
engine for maximum work output with generalized radioactive heat transfer law. International
Journal of Energy and Environment. Volume 2, Issue 4, 2011, pp.723-738.
[153] M. Abbas, N. Said and B. Boumeddane, ‘Thermal Analysis of Stirling Engine Solar
Driven’, Revue des Energies Renouvelables, Vol. 11, N°4, pp. 503 – 514, 2008.
[154] H. Lamrani, ‘Modélisation Energétique des Moteurs Stirling’, Thèse de Doctorat, Ecole
des Mines de Nantes, France, 1995.
[155] Schmidt G 1871 The Theory of Lehmann's Calorimetric Machine Z. ver. Dtsch. Ing. 15
part I
[156] Finkelstein T., Generalized thermodynamic analysis of Stirling engines. Paper 118B,
SAE Winter Annual Meeting 1960.
[157] W.R. Martini, ‘Stirling Engine Design Manual’, DOE/NASA/3152–3178/1 NASA CR-
135382; 1978.
100
[158] Urieli I. and Berchowitz D. M. , Stirling cycle engine analysis. Adam Hilger, Bristol
1984.
Chapitre 4
Cet article a fait l’objet d’une publication dans: Eastern European Journal of
Enterprise Technologies, Kharkov Ukraine - 2011- N° 5/6 (53) - 63 p - pp. 37-43 - ISSN
1729-3774. Il a été presenté aussi lors de la conférence scientifico-méthodologique et de
recherche scientifique, section des moteurs à combustion interne et sous-section de
l’amélioration des moteurs à piston à combustion interne, qui a eu lieu à l’Université
Technique de l’Automobile et des Ponts et Chaussées de Moscou Russie le 29 Janvier - 07
Février 2013. La version originale de l’article en langue russe se trouve dans l’annexe A.
Ce chapitre présente les travaux effectués pour l’article I, qui porte sur la méthode de
calcul de la composition d’équilibre des produits de combustion d’un moteur à gaz à allumage
par étincelle.
Pour déterminer la composition d'équilibre des produits de combustion dans le moteur
à gaz à allumage par étincelle, on suggère d'utiliser un système de dix (10) équations avec dix
(10) inconnues, fondé sur six (06) réactions chimiques, trois (03) équations de la balance
matérielle et l'équation de la loi de Dalton. Pour la résolution de ce système d'équations, il est
recommandé d'utiliser la méthode Zeldovitch-Polarny. Ensuite, on a comparé les résultats des
calculs avec ceux expérimentaux effectués au préalable sur le moteur à gaz, quatre (04)
temps, à allumage par etincelle 4GCH7.5/7.35.
4.1. Résumé
Mots clefs: Méthode de calcul, processus de combustion, toxicité, moteur à gaz, composition
d'équilibre.
4.2. Introduction
Actuellement pour calculer le taux des composants toxiques dans les gaz
d'échappement d'un moteur à allumage par étincelles, on utilise essentiellement la méthode de
calcul de la composition d'équilibre des produits combustibles, proposé par le professeur V.A.
Zvonov V. A [1]. Pour le calcul de la concentration des NOx, la méthode donnée est
complétée par l'équation cinétique de Zeldovitch [1].
Le but de cette étude est l'analyse des approches existantes de résolution du système
d'équations non linéaires dans la méthode de calcul du taux des composants toxiques dans les
gaz d'échappement du moteur à allumage par étincelle et l'élaboration des recommandations
quant au choix de la méthode de calcul.
CO + H 2 O ←
→ CO 2 + H 2
Cependant, dans les conditions des moteurs à combustion interne à gaz avec allumage
par étincelle avec α ≥ 1 le déroulement d'une telle réaction est pratiquement impossible.
L’approche la plus rationnelle est proposée dans [5], où on utilise un système de 10
équations avec 10 inconnues. Le système contient le minimum nécessaire d'équations pour
trouver les pressions partielles du CO et des NOx, de plus ces équations sont utiles pour
l'utilisation dans les conditions des moteurs à combustion interne à gaz avec allumage par
étincelle.
Il est supposé que dans les produits de combustion, les réactions suivantes se
déroulent:
104
1.CO 2 ← → CO+0.5O 2 ;
2.H 2 O ← →
2H +0.5O ;
2
3.H 2 O ←→
OH+0.5H 2 ;
(4.1)
4.H 2 ← →
2H;
5.O 2 ← → 2O;
6.0.5N 2 +0.5O 2 ←→ NO.
Il est considéré qu'à titre du combustible, on utilise le méthane et les 10 gaz suivants
entrent dans la composition des produits de combustion: CO2, CO, H2O, H2, O2, N2, NO, OH,
O et H.
Les équations d'équilibre chimique sont à la base du système d'équations:
1
PCO ⋅ PO2 2
= K1 , (4.2)
PCO2
1
PH2 ⋅ PO2 2
= K2 , (4.3)
PH2O
1
POH ⋅ PH2 2
= K3 , (4.4)
PH2O
PH 2
= K4 , (4.5)
PH2
PO 2
= K5 , (4.6)
P O2
PNO
1 1
= K6 , (4.7)
PO2 ⋅ PN2
2 2
Les valeurs des constantes d'équilibre K1 … K6 sont calculées à partir des tableaux
dans [5, 7].
Les équations (4.2)…(4.7) sont complétées par trois équations de la balance matérielle
(4.8)…(4.10)
105
(
OCar + χ ОС O AIR μ О 2 PCO2 + PCO + A
= ,
)
( )
(4.8)
СCar + χ ОС C AIR μ С PCO + PCO 2
où A = PH2O + 2 PO2 + POH + PNO + PO ; OCar - la part massique de l'atome de l'oxygène dans le
carburant; O AIR - la fraction massique de l'atome d'oxygène dans l'air; СCar - la fraction
carbone dans l'air; χ ОС - le rapport réel des masses des atomes O et C dans le mélange
H Car + χ НС ⋅ H AIR μН ⋅ ( B)
= , (4.9)
(
СCar + χ НС ⋅ CAIR μ С ⋅ PCO + PCO
2
)
carburant; H AIR - la fraction massique de l'atome d’hydrogène dans l'air; χ НС - le rapport réel
des masses des atomes H et C dans le mélange imbrulé; μ Н - masse moléculaire de l'atome
d’hydrogène, kg/kmol:
(
N Car + χ NC ⋅ N Air μ N ⋅ 2 ⋅ PN2 + PNO
= ,
)
( )
(4.10)
СCar + χ NC ⋅ CAir μ С ⋅ PCO2 + PCO
massique de l'atome d'azote dans l'air; χ NC - le rapport réel des masses des atomes N et C dans
Dans la première approximation, supposons que les pressions partielles des gaz
suivants sont égales à zéro,
PO=
2
= P=
PNO O P=
H = 0.
POH (4.12)
où,
μ C OCar + χ ОС ⋅ O Air
=
W ⋅ . (4.14)
μ O СCar + χ ОС ⋅ CAir
1
2 ⋅ PH2O + 2 ⋅ PH2
= W2 , (4.15)
PCO2 + PCO
où,
H Car + χ НС ⋅ H Air
W= μC ⋅ . (4.16)
СCar + χ НС ⋅ C Air
2
2 ⋅ PN2
= W3 , (4.17)
PCO2 + PCO
où,
μ C N Car + χ NC ⋅ N Air
=
W ⋅ . (4.18)
μ N СCar + χ NC ⋅ C Air
3
PCO ⋅ PH2O K1
= . (4.20)
PCO2 ⋅ PH2 K2
107
Et l’on obtient,
=
PH2O PCO2 (W1 − 2 ) + PCO (W1 − 1) ; (4.22)
𝑃𝐻2 = 0.5𝑊2 �𝑃𝐶𝑂2 + 𝑃𝐶𝑂 � − 𝑃𝐶𝑂2 (𝑊1 − 2) − 𝑃𝐶𝑂 (𝑊1 − 1). (4.23)
PN2 = (
0.5 ⋅ W3 PCO2 + PCO . ) (4.24)
𝑃𝐶𝑂2 + 𝑃𝐶𝑂 + 0.5𝑊2 𝑃𝐶𝑂2 + 0.5𝑊2 𝑃𝐶𝑂 + 0.5𝑊3 𝑃𝐶𝑂2 + 0.5𝑊3 𝑃𝐶𝑂 = 𝑃 (4.26)
( PCO 2
)
+ PCO (1 + 0.5W2 + 0.5W3 ) =
P, (4.27)
ou,
P
PCO2 + PCO = (4.28)
1 + 0.5W2 + 0.5W3
Ayant mis la partie droite de l'équation (4.28) dans l'égalité (4.24), on obtient l'expression
suivante pour le calcul de la pression partielle de l'azote:
0.5W3 P
PN 2 = , (4.29)
1 + 0.5W2 + 0.5W3
W3 P
PN 2 = . (4.30)
2 + W2 + W3
𝑃𝐶𝑂 + 0.5𝑊2 𝑃𝐶𝑂 + 0.5𝑊3 𝑃𝐶𝑂 = 𝑃 − 𝑃𝐶𝑂2 − 0.5𝑊2 𝑃𝐶𝑂2 − 0.5𝑊3 𝑃𝐶𝑂2 (4.31)
On divise l’équation obtenue par 0.5(2 + 𝑊2 +𝑊3 ), et on obtient la formule pour calculer la
pression partielle du CO dans les produits de combustion:
2P
=PCO − PCO2 . (4.33)
2 + W2 + W3
Afin d’obtenir la formule pour calculer la pression partielle de Н2О dans les produits de
combustion, on effectuera les actions suivantes.
En plaçant (4.33) dans (4.22), on obtient:
2∙𝑃
𝑃𝐻2 𝑂 = 𝑃𝐶𝑂2 ∙ (𝑊1 − 2) + �2+𝑊 +𝑊 − 𝑃𝐶𝑂2 � (𝑊1 − 1) (4.34)
2 3
2𝑊 𝑃 2𝑃
𝑃𝐻2 𝑂 = 𝑃𝐶𝑂2 𝑊1 − 2𝑃𝐶𝑂2 + 2+𝑊 1 𝑊 − 𝑊1 𝑃𝐶𝑂2 − 2+𝑊 + 𝑃𝐶𝑂2 . (4.35)
2+ 3 2+ 𝑊3
On obtient ultimement:
2 P (W1 − 1)
=PH 2O − PCO2 . (4.36)
2 + W2 + W3
Pour obtenir l'équation pour le calcul de la pression partielle Н2, nous accomplirons ainsi la
substitution de l'expression (4.33) dans l'équation (4.23):
où,
2P
=С − PCO2 .
2 + W2 + W3
2𝑃 𝑊 2𝑃
𝑃𝐻2 = 0.5𝑊2 2+𝑊 +𝑊 − 𝑊1 𝑃𝐶𝑂2 + 2𝑃𝐶𝑂2 − 2+𝑊1 +𝑊 + 𝑊1 𝑃𝐶𝑂2 + 𝐶. (4.38)
2 3 2 3
On obtient:
W2 P − 2W1P + 2 P
=PH 2 + PCO2 , (4.39)
2 + W2 + W3
P ( 2 − 2W1 + W2 )
=
PH 2 PCO2 + . (4.40)
2 + W2 + W3
Pour mettre au point l'équation du second degré relativement à la pression partielle du СО2,
nous nous servirons de l'égalité (4.20). En transférant tous les composants de l'équation (4.20)
à gauche, on obtient:
PCO ⋅ PH 2O K1
− =
0. (4.41)
PCO2 ⋅ PH 2 K2
K 2 PCO PH 2O − K1PCO2 PH 2 =
0. (4.42)
2P
= M1 . (4.43)
2 + W2 + W3
2 P (W1 − 1)
= M2 . (4.44)
2 + W2 + W3
P ( 2 − 2W1 + W2 )
= M3 . (4.45)
2 + W2 + W3
110
On réécrit l’équation (4.42) pour que toutes les valeurs des pressions partielles soient
exprimées par l’intermédiaire de la pression partielle du СО2,
2 2
𝐾2 𝑀1 𝑀2 − 𝐾2 𝑀1 𝑃𝐶𝑂2 − 𝐾2 𝑀2 𝑃𝐶𝑂2 + 𝐾2 𝑃𝐶𝑂2
− 𝐾1 𝑃𝐶𝑂2
− 𝐾1 𝑃𝐶𝑂2 𝑀3 = 0 (4.47)
En formattant les équations précedentes on obtient l'équation du second avec comme variable
𝑃𝐶𝑂2 :
2 (𝐾
𝑃𝐶𝑂2 2 − 𝐾1 ) + 𝐾2 𝑀1 𝑀2 − 𝑃𝐶𝑂2 (𝐾2 𝑀1 + 𝐾2 𝑀2 + 𝐾1 𝑀3 ) = 0 (4.48)
K 2 M1 + K 2 M 2 + K1M 3 − D
PCO2 = . (4.50)
2 ( K 2 − K1 )
Ayant trouvé selon la formule (4.30) la pression partielle de N2 et ayant résolu l'équation
(4.48) par rapport à CO2, on trouve les pressions partielles CO, H2O et H2 successivement à
l'aide des formules (4.33), (4.36) et (4.40).
Ensuite nous calculons les valeurs des pressions partielles des cinq composants restant des
produits de combustion.
La pression partielle О2 est égale à:
2
K1 ⋅ PCO2
PO2 = . (4.51)
P CO
K3 ⋅ PH 2O
POH = . (4.52)
PH 2
La pression partielle de l'hydrogène:
111
=
PH K 4 ⋅ PH 2 . (4.53)
=
PO K5 ⋅ PO2 . (4.54)
=
PNO K 6 PN2 ⋅ PO2 . (4.55)
D’après la référence [1] on adment que les émissions d’oxydes d’azote sont consitués
de monoxyde d’azote à 97 %.
0. O 2 ←→ O + O;
→ NO + N;
k1d
1. N 2 + O ←
k1i
k 2d (4.56)
2. N + O 2 ← →
NO + O;
k 2i
→ 2 ⋅ NO.
k 4d
3. N 2 + O 2 ←
k 4i
Dans le système de réactions chimiques (4.56) k1d, k1i, k2d. k2i. k4d. k4i - sont les constantes des
vitesses des réactions de (1), (2) et (3) respectivement.
L'équation de calcul cinétique pour la formation de NO pour le système des réactions (4.56) a
la forme suivante [1]:
drNO p ⋅ rN 2 2
rNO k1d ⋅ rO k 4d ⋅ rO2
=− 1 + , (4.57)
dϕ 254.4nTp.c. K 4(NO) rO2 rN 2
1 + k1i ⋅ rNO 2
k 2d ⋅ rO2
112
K 4( NO ) 21.12ехр(−181000 / ( RTp.c )) ;
= (4.61)
Début
i=0, j=0
Calcul de NOx
i=i+1
Non
j=1
Calcul d’état d’équilibre
1-ère approximation
j=j+1
Oui
Non
i = ip.com
Oui
FIN
Fig. 4.2 Les résultats des calculs de la toxicité du moteur 4GCH7.5/7.35 (n = 3650 min-1)
1 - Méthode Zeldovitch-Polarny; 2 - Méthode de Koutsenko (Newton); 3 - Expérimentale
Comme on peut voir sur la figure 4.2, la différence dans les résultats du calcul selon la
méthode numérique de Newton et selon la méthode Zeldovitch-Polarny est pratiquement
inexistante. L’erreur dans les calculs est inférieure à 10% par rapport aux résultats
expérimentaux. De plus, la méthode Zeldovitch-Polarny est beaucoup plus facile à réaliser et
ne nécessite pas des profondes connaissances en programmation et permet d’utiliser de
simples complexes mathématiques de programmation.
4.7. Conclusions
1. Pour déterminer la composition d'équilibre des produits de combustion dans le moteur
à gaz à allumage par étincelle, on suggère d'utiliser un système de 10 équations avec 10
inconnues, fondé sur 6 réactions chimiques, 3 équations de la balance matérielle et l'équation
de la loi de Dalton.
115
3. La comparaison des résultats des calculs avec ceux expérimentaux a démontré que la
différence dans les résultats de calcul selon la méthode numérique de Newton et selon la
méthode Zeldovitch-Polarny est pratiquement inexistente. Les résultats analytiques et
expérimentaux présentent une bonne corrélation et la différence est inférieure à 10%. Par
ailleurs, la méthode Zeldovitch-Polarny est nettement plus facile à appliquer.
4.8. Références
[1] V.A.Zvonov. La toxicité des moteurs à combustion interne. Manuel scolaire pour les
écoles supérieures / V.A.Zvonov - 2-ème édition avec modifications et suppléments. Moscou,
éditions « Machinostroyéniyé », 1981. - 154 p.
[6] Bade Shrestha S.O. A Predictive Model for Gas Fueled Spark Ignition Engine
Applications / Bade Shrestha S.O., Karim G.A. - Calgary, CA: University of Calgary, 1999. -
18 p. - (Preprint / University of Calgary: SAE 1999-01-3482).
Chapitre 5
Ce chapitre présente les travaux effectués pour l’article II, qui porte sur le choix et le
fondement des paramètres régulateurs du moteur à gaz suralimenté converti à partir du moteur
diesel, à l'exemple du moteur 6GCHN13/14. Pour ce moteur, on précise le modèle mathématique
du cycle de combustion avec l’indice variable de la combustion de Wiebe.
seulement de calculer les concentrations des oxydes d'azote NOx dans les gaz d'échappement
avec une précision suffisante, mais aussi d'obtenir les données nécessaires pour le calcul de la
probabilité de détonation et son estimation quantitative .
Le taux de compression du moteur diesel d'origine est beaucoup plus élevé que le
moteur diesel modifié, ce qui donne aux pièces de moteur à gaz a aspiration forcé une grande
fiabilité. Par conséquent, le choix du taux de compression est influencé seulement par la
probabilité de détonation se produisant aux charges maximales .
Un critère de détonation sans dimension 𝑘𝑐 peut ainsi être défini comme étant le
rapport de deux énergies relatives libérées au cours de la combustion .
La résolution du problème d'optimisation des paramètres de réglage du moteur à gaz a
été réalisée au moyen de la technique de l'espace des paramètres en utilisant les suites de
Sobol.
La concentration en NOx doit être mesurée par ailleurs en trois points de la zone de
contrôle, choisis au hasard par le service technique réalisant les essais. Les émissions de NOx
mesurées pour ces trois points doivent être proches des valeurs calculées pour les points
d’essai de la norme les plus proches. Cette vérification vise à assurer l’efficacité du contrôle
des émissions dans la plage des régimes d’exploitation types du moteur. Les trois points de
mesure additionnels sont choisis en utilisant des méthodes statistiques agréées de prélèvement
aléatoire (s’il y a écart entre les valeurs mesurées et les valeurs calculées, le moteur ne
satisfait pas l’essai de réception par type, voir l’article ci-dessous).
5.1 Abstract
Many cities and countries today have programs to convert older, polluting diesel
transit buses and trucks to run on clean, economical natural gas. Properly implemented, this
is an excellent way to quickly reduce fuel costs, clean up the air and reduce noise with
minimum capital costs. A poorly executed conversion program, however, can lead to higher
exhaust emissions, much higher fuel consumption, unacceptable power losses, poor
durability and high maintenance costs.
This article gives the results of our investigation of a forced aspiration gas engine
6GCHN13/14 converted from diesel. The mathematical model of the combustion process
specified for the gas engine uses a variable Wiebe combustion factor.
Use of the developed characteristic maps has shown improvement of engine power
indicators in comparison with that of a naturally aspirated engine, as well as a decrease
in average operational emissions with the use of adjustable forced aspiration.
5.2 Introduction
Converting diesel engines of trucks into gas spark ignited engines is now a subject of
actuality [1]. However the execution of this task faces a number of difficulties. For example,
the chamber of combustion of an original diesel engine can’t support the thermal load created
after conversion to gas spark ignition without further modifications. Also the stringent
regulation imposed on automobile internal combustion engines concerning gas emissions
requires the use in the converted engines of the concept of "lean burning" [1].
Application of this concept causes a decrease in the power of the engine, especially in
modes of maximum loads. To compensate these losses of engine power of the converted
diesel engine, it is expedient to use supercharging.
Projects to convert diesel engines to spark-ignited forced-aspiration gas engines using
this concept of “lean burning” are found all over the world. In Russia such work on a KAP-
740 diesel engine has been reported in [2]. In the USA, a similar program is pursued by
companies such as John Deere [3], Cummins [4], Detroit Diesel [5], Caterpillar [6] and
others. In Europe such gas fuel conversion activity takes place for public transport, as
exemplified by Iveco City Class buses [7]. In Ukraine similar work is carried out in LNTU
[8], KHNAHU [1] and other centres.
Analysis of the literature has shown that in the majority of studies on gas engines
converted from diesel engines, forced-aspiration is either used immediately, or planned for the
near future.
The test bench for carrying out the measurements, and also for verifying the validity
of the analytical models, is shown in Figure 7.1. This work was executed in the Gas
Laboratory of the Internal Ignition Engines Section in KHNAHU. The bench was created
around a gas spark ignited engine with supercharging 6GCHN13/14, converted from a diesel
engine YMP-236.
120
In order to permit a smooth change of drag torque on the engine shaft, during tests, the bench
is equipped with an electric loading device MPB-100 (1 in Figure 5.1). The power from the
engine is delivered to the balancing machine via a mechanical gearbox VAZ2103 2 and a
cardan joint 3. This gearbox can work in a wide frequency range of rotation of the engine
crankshaft. The bench is equipped with an additional air-fan 4, and an additional heat
exchanger 5, with which it is possible to model a temperature regime which approaches real
running conditions. For the same purpose, the bench also includes a noise muffler 6 and a resonator
7 from the regular exhaust system of a Daewoo Model Sens. A damper valve 8 is installed to
measure the exhaust system resistance.
In order to explore a variety of gases, the bench is equipped with a fourth-generation gas
system which includes: several cylinders with a variety of gases 9, a two-level gas pressure
reducer 10, high and low pressure lines, a rail with gas injectors 11 and a control panel 12. The
fuel gas is selected by a series of switches 13 on board 14.
Investigations in the field of operator workplace ergonomics [9] have shown that the
greatest informational content is provided by visual observations of low-frequency processes.
121
Hence a series of registering devices are included, such as damped analogical sensors (15, 16, 17, 18,
19), and a LED digital display (20 - 25). The bench also includes a gas flow meter (26) and an air
flow meter (27) to measure these input parameters. The stand is also equipped with an oil
temperature measuring device (21) and an engine load control system (28). Finally, an engine
electronic control unit 29 is mounted on the test bench. The stand is equipped with all the
necessary equipment to determine actual and indicated indicators of power and fuel economy
and toxicity.
θ + 18
A= , (5.2)
40
B = α+0.00025 , (5.3)
0.005 ⋅ ηv + 0.005
C= , (5.4)
0.01
n + 1100
E= , (5.5)
3000
122
𝜋𝑐 +0.001
𝐹= , (5.6)
0.95
where ϕi is the relative burn angle of the combustion phase i, ϕi = 0…1 (0-100% burn
duration); n is the frequency of rotation of the crankshaft of the engine, in min-1; α is the air-
fuel excess ratio; ηv is the volumetric efficiency; θ is the spark advance, crank angle before
the top dead center (TDC), and πc is the rate of charging pressure rise. Moreover, the burn
duration:
where, 𝜋𝑐 or CPR is total pressure ratio is the ratio of the compressor discharge air pressure to
the ambient air pressure (the ratio of the air total pressure exiting the compressor to the air
pressure entering the compressor).
𝛾
𝑝02 𝑁𝑡𝑐 𝑚̇𝑐 �𝑇01
𝜋𝑐 = = 𝑓𝜋𝑐 � , � = (𝑇02 ⁄𝑇01 )𝛾−1 , (5.9)
𝑝01 �𝑇01 𝑝01
The total temperature ratio 𝑇02 ⁄𝑇01 across the compressor is related to the pressure ratio by
the isentropic flow equations.
𝑝01 : air pressure entering the compressor.
𝑝02 : air total pressure exiting the compressor.
𝑇02 : total temperature exiting the compressor.
𝑇01 : temperature entering the compressor.
γ : ratio of specific heats.
As the concentration of toxic components in the exhaust gases cannot be calculated
with sufficient accuracy with the model of Wiebe because the impossibility of determining
the temperature distribution in the cylinder, a two-zone model as described in [12] has been
added in our model. The advantage of the two-zone model is that it allows not only the
calculation of the concentrations of the nitrogen oxides NOx in the exhaust gases with
sufficient accuracy, but also to obtain the data necessary for the calculation of the possibility
of detonation occurring, and its quantitative estimation.
123
5.4.1 Choice of the compression ratio and the maximum charging pressure
The compression ratio choice is determined by two factors: the mechanical durability
of the engine components and the probability of early detonation occurring at maximum
loads.
The compression ratio of the original diesel engine is much higher than the modified
diesel engine. Which gives the forced aspiration gas engine parts a high reliability. Therefore,
the choice of compression ratio is influenced only by the second factor.
Autoignition and knock (the sound of an auto-ignition initiated pressure wave inside
the engine cylinder) occurs when end gases release sufficient amount of energy above a
specific threshold.
An appropriate criterion for knock in spark ignition engines must be based on the total
energy released within the end gas due to auto-ignition reaction activity per unit of the
corresponding instantaneous volume. Such a specific energy release needs to be compared for
any engine, fuel and operating conditions to the corresponding total energy release that would
take place normally due to regular flame propagation. The latter needs to be related to the
total cylinder size and volume. A dimensionless knock criterion 𝑘𝑐 described in [13] can thus
be defined as the ratio of the following two relative energy releases during the course of the
combustion process:
where the numerator represents the total energy released by the unburned end gas pre-ignition
chemical reaction activity per unit of the instantaneous charge volume (Vt), while the
denominator is the total energy that can be released normally through flame propagation over
the whole cycle per unit of cylinder swept volume (Vo). Vt , the instantaneous volume at time
t, is the sum of the combustion chamber volume, the volume displaced by the piston
(depending on the crank angle) and mechanical deformations produced by the gas pressure
and inertial efforts; Vo , the cylinder volume swept by the piston displacement, is the full
stroke movement of the piston.
The energy released by the end gas self-reactions is defined as:
𝑡
𝐻𝑢𝑡 = ∫𝑡 𝑑𝐸𝑢 = (ℎ𝑠𝑡 − ℎ𝑡 )𝑚𝑢 , (5.11)
𝑠𝑡
124
where h is the enthalpy of the mixture per unit of mass (ℎ = ∑ 𝑐𝑖 ℎ𝑖 ) and 𝑚𝑢 is the mass of
the unburnt zone which is the end gas region at any instant of time t. Subscripts st and t
indicate the values at the spark discharge and at any instant of time, respectively.
The energy released due to normal combustion can be simplified as:
𝐻𝑜 = ℎ𝑜 𝑚𝑜 , (5.12)
where ℎ𝑜 is the effective heating value of the fuel and 𝑚𝑜 is the initial or total mass.
Since the spark timing is usually set such that much of the combustion process takes
place in the vicinity of the top dead center position so as to achieve optimum efficiency and
power [𝑉𝑡 ≈ 𝑉𝑐 (clearance volume)], then the volume ratio 𝑉𝑅, where 𝑉𝑅 = 𝑉0⁄𝑉𝑡 , can be
approximated to (ε – 1) [the ratio is calculated by the following formula: 𝜀 = 𝑉0⁄𝑉𝑐 + 1 ≈
≈ 𝑉0⁄𝑉𝑡 + 1]. Consequently, equation 7.9 may be written as:
The knock criterion of equation (5.13) can be used to test for the onset of knock and to
represent effectively engine knock intensity in gas-fueled spark ignition engines.
The example described in [13] and other similar results show that the operating
condition for the knock limit may be determined theoretically by assuming a value of 1.5 for
maximum knock criterion value at borderline knock operation. Excess of this value means
possible detonation presence at the end of the on the regime, and its intensity is directly
proportional to the value Kcmax. The detonation occurring is influenced by two major factors:
compression ratio ε and rate of charging pressure rise 𝜋𝑐 .
Supercharging is designed to increase the density of the air into the cylinders by
increasing its pressure in the intake manifold. Engine power thus increases as the work cycle
for a supercharged engine is much superior to that of a naturally aspirated engine. However
the compression ratio should not be excessive so as to avoid reaching maximum cycle
pressures which may be too high. Research [6] has shown that it is possible to obtain the best
engine power if the compression ratio is reduced, thus increasing charging pressure to a
maximal value. The design of the turbocharger and the mechanical durability of parts of the
combustion chamber can reach this value.
Calculations using the criterion Kc have shown that at the maximum rate of charging
pressure rise provided by a ТКР-9-12-07 turbo-compressor on regimes of maximum load
(πc_max = 1.45), a compression ratio ε = 11.8 allows to operate without premature detonations.
125
Table 5.1 Range of variation of engine parameters for optimizing engine performance
α θ n 𝜋𝑐 𝜑𝑡𝑟
° crank angle (deg) before
- min-1 - %
TDC
1 5 800 1 0
1.5 40 2100 1.45 100
The ranges of change of variable parameters for solving the optimization problems of
regulation parameters listed in Table 5.1 and there are selected from the following
considerations.
For α < 1 the fuel burns only partially and, as a consequence, the emissions of CO and
СН increase strongly. For α > 1.5, the NOx concentration in the exhaust gases decreases
significantly and CH emissions suddenly increase. However because of a sudden increase in
126
cyclic instability of the operating process and an increase in frequency of pre-ignition, the
power and the fuel economy of gas engines decrease considerably,.
The limits for change of spark advance θ, the crank angle (deg) before TDC, was
selected on the basis of the experimental data received as a result of preliminary tests on the
gas engine 6GCHN13/14.
The range for change of the frequency of rotation of the crankshaft n, min-1 was
decided from considerations of guarantee of stable working conditions of the engine at any
combination of other factors. The throttle positions go from fully closed (𝜑𝑡𝑟 = 0 %) to fully
open (𝜑𝑡𝑟 = 100 %).
The rate of aspiration pressure rise 𝜋𝑐 (due to transporting the fuel-air mixture through
the compressor) varies from total absence of forced aspiration (𝜋𝑐 = 1) to the maximum rate of
aspiration (𝜋𝑐 = 1.45).
The analytical and experimental study carried out has shown that operating the gas
engine without forced aspiration allowed it to reach the nominal power regime of the basic
unmodified diesel engine Ne = 130 kW only at α ≈ 1.
The strength margins of the engine allows an increased charging pressure in a nominal
regime to 𝜋𝑐 = 1.45, allowing it to reach a theoretical power of 180 kW at α = 1. However
this leads to a thermal overstrain of parts of the engine and pre-detonation occurring.
To reduce the influence of these last factors, the engine must be operated at a leaner
mixture. Research has shown that to assure the engine operating without pre-detonation in
nominal regime, it is necessary to make the mixture lean to α ≈ 1.3. This operational limit at
𝜋𝑐 = 1.45 will produce power Ne = 130 kW at n = 2100 min-1 and torque Me = 670 N·m at
n = 1500 min-1.
NOx emissions in regimes of maximum load at 𝜋𝑐 = 1.45 and α ≈ 1.3 show a decrease
of 35-45 % in comparison with operating at 𝜋𝑐 = 1.45 and α ≈ 1. So, reducing the speed at the
maximum relative load (Pe = 100%) is appropriate to carry through the increase of excess air
ratio to α ≈ 1.5. The effectiveness of ignition at high values of α provides high concentration
of fuel near the spark plug provided by the forced aspiration. So, load reduction to Pe ≈ 0.6 is
the value of α is 1.05. This approach can reduce NOx emissions on medium engine speed on
60…70 %, compared to the naturally aspirated version, operating at low values of α.
Thus, on the engine 6GCHN13/14 the mixed regulation of power is offered to use: by
means of change of throttle opening rate and α change. Besides, there is an appearance of
third operating factor - rise of pressure charging rate 𝜋𝑐 .
Figure 5.2 to 5.5 show the characteristic maps for the optimized engine performance,
namely the dependence of the respective parameters α, 𝜋𝑐 , 𝜑𝑡𝑟 and θ respectively at n and Pe .
These data will serve as reference values for a microprocessor control system of the engine.
From figure 5.6, it is visible that at use of forced aspiration on regimes of external
high-speed characteristics, Ne increases on 26…31 %, Me: on 35…41 %, although specific
effective fuel consumption ge increases on 7…9 %. Despite the rise of power, α increases
lightly on 3…4 %.
When converting the average power diesel engine to gas engine, adjustable forced
aspiration can be used not only for the preservation of the power at the same level as the base
engine, but also to reduce emissions of nitrogen oxides .
In Figures 5.7 and 5.8, we compare the emissions from the naturally aspirated engine
6GCH13/14 with those of the forced aspiration engine 6GCHN13/14.
131
The ESC test cycle [15] has been introduced, for emission certification of heavy-
duty diesel engines in Europe. The ESC is a 13-mode, steady-state test.
The ESC test is characterized by high average load factors and very high exhaust gas
temperatures.
The choice of speeds to determine the test modes 𝑛𝐴 , 𝑛𝐵 , 𝑛𝐶 is implemented with the
help of auxiliary speeds 𝑛ℎ𝑖 and 𝑛𝑙0 . The high speed 𝑛ℎ𝑖 is determined by calculating 70% of
the declared maximum net power. The highest engine speed where this power value occurs
(i.e. above the rated speed) on the power curve is defined as 𝑛ℎ𝑖 . Here 𝑛ℎ𝑖 = 2200 min-1. The
low speed 𝑛𝑙0 is determined by calculating 50% of the declared maximum net power. The
lowest engine speed where this power value occurs (i.e. below the rated speed) on the power
curve is defined as 𝑛𝑙0 . Here 𝑛𝑙0 = 1500 min-1.
The engine speeds to be used during the test are then calculated from the following
formulas:
In accordance with the dependencies 𝑛𝐴 =1675 min-1, 𝑛𝐵 =1850 min-1, 𝑛𝐶 =2025 min-1.
Each mode is characterized by the mode significance factor δз, demonstrating the
conditionnal relative proportion of time during which, the engine work at this mode.
Tests are conducted at low idle speed and at 4 values of the relative torque: 0.25, 0.5,
0.75 and 1.0. The results are shown in Table 5.2.
134
Table 5.2 The results of tests with forced aspiration engine 6GCHN13/14 with the ESC test
cycle.
№ δз М кр n NOx CO CH
– – % min-1 g/(kW·h) g/(kW·h) g/(kW·h)
1 0.15 0 XX 0.32 0.96 2.1
2 0.08 1.0 nA 3.05 1.22 1.24
3 0.1 50 nB 2.52 1.22 1.12
4 0.1 75 nB 2.74 1.19 1.16
5 0.05 50 nA 2.44 1.13 1.15
6 0.05 75 nA 2.66 1.17 1.16
7 0.05 25 nA 1.86 1.16 1.31
8 0.09 100 nB 3.8 1.15 1.19
9 0.1 25 nB 1.92 1.12 1.05
10 0.08 100 nC 3.55 1.15 1.1
11 0.05 25 nC 2 1.04 0.94
12 0.05 75 nC 2.82 1.11 1
13 0.05 50 nC 2.68 1.2 1.07
Table 5.3 The results of tests naturally aspirated engine 6GCH13/14 with the ESC test cycle.
№ δз М кр n NOx CO CH
-1
– – % min g/(kW·h) g/(kW·h) g/(kW·h)
1 0.15 0 XX 0.44 0.98 2.16
2 0.08 1.0 nA 4.24 1.24 1.28
3 0.1 50 nB 3.50 1.24 1.15
4 0.1 75 nB 3.81 1.21 1.19
5 0.05 50 nA 3.39 1.15 1.18
6 0.05 75 nA 3.70 1.19 1.19
7 0.05 25 nA 2.59 1.18 1.35
8 0.09 100 nB 5.28 1.17 1.23
9 0.1 25 nB 2.67 1.14 1.08
10 0.08 100 nC 4.93 1.17 1.13
11 0.05 25 nC 2.78 1.06 0.97
12 0.05 75 nC 3.92 1.13 1.03
13 0.05 50 nC 3.73 1.22 1.1
135
components the gas naturally aspired engine 6GCH13/14 have show, gr/(кW⋅h) [6]:
g NOx = 3.28; g CO = 1.15; g CH = 1.31.
5.8 Conclusions
Analysis of the literature showed that the use of adjustable forced aspiration to the gas
engine can eliminate deficiencies of the concept of a "poor combustion", retaining its
advantages.
Tests on the diesel-to-gas converted engine 6GCHN13/14 with forced aspiration and
using the characteristic maps developed in part 6 and summarized in Figures 5.2 to 5.5 have
shown a definite improvement of the engine power indicators in comparison with the
naturally-aspired configuration: that can be summarized as a 7-9% improvement in fuel
consumption, 26-31 % in Ne, 35-41 % in Me and a 10-28% reduction in noxious gas
emissions.
Tests of the gas engine 6GCHN13/14 on 13-regime cycle ESC have shown that
decrease in average operational emissions at use of adjustable charging on the gas engine
6GCHN13/14 in comparison with a naturally aspired variant has give: on NOx: 28 %, on
CO: 10 %, on СН: 18 %.
5.9 References
[3] D.L. McCaw and W.A. Horrell. On-Road Development of John Deere 6081 Natural Gas
Engine: Final Technical Report - Deere & Company - 2001. - 46p.
[4] M.M. Kamel. Development of a Cummins ISL Natural Gas Engine at 1.4g/bhp-hr NOx +
NMHC Using PLUS Technology: Final Report - Cummins Westport, Inc - 2005. - 18p.
[6] G.B. Cox, K.A. DelVecchio, W.J. Hays. Development of a Direct-Injected Natural Gas
Engine System for Heavy-Duty Vehicles: Final Report - Caterpillar Inc. - 2002. - 67 p.
[7] N-O Nylund, J. Laurikko, M. Ikonen. Pathways for Natural Gas into Advanced Vehicles -
Brussel: IANGV. - 2002. - 105 p.
[8] В. Zakartschuk, O.P.Sitovsky, І.С. Kozatschuk. Analytical and Experimental Study of Gas
engine, Converted from Diesel/Automobile Transport. - 2005- №16 - p. 276 - 278.
[10] Wiebe I.I. The new about working cycle of the engine. - Moscow: Mashgiz, 1962 - 271p.
[11] Kabanov Alexander Nikolaevich. Decrease of Emissions of Harmful Substances with the
Exhaust Gases of Transport Diesel Engines by their Converting in Gas engines: Dissertation
Candidate Technical Science: 05.05.03 - Х, 2007. - 206 pages.
[13] Attar A.A. Optimization and Knock Modelling of a Gas Fueled Spark Ignition Engine:
PhD thesis: Mechanical Engineering / Calgary, Canada, 1997. - 248 p.
Chapitre 6
Ce chapitre a fait l’objet d’une publication dans la revue scientifique American Journal of
Vehicle Design, 2013, Vol. 1, N° 2, 21-24:
Rachid Maamri, Alexandre Petrovitsh Martshenko, O. O. Osetrov, Yves Dubé, Lotfi Toubal,
Agbossou Kodjo, Analyze and mathematical modelling of the combustion process of one-
cylinder spark ignited hydrogen fueled engine.
Ce chapitre présente les travaux effectués pour l’article III, qui porte sur l’analyse de
la phase de combustion de l’hydrogène dans un moteur 4 temps monocylindrique à allumage
commandé, l’élaboration du modèle mathématique, qui prend en considération les
particularités de la phase de combustion et l’estimation de l’influence des paramètres
constructifs et de réglage du moteur. A l’aide d’une simulation, on va valider une étude
expérimentale effectuée au préalable sur un moteur à combustion interne à allumage
commandé CFR, 4 temps, monocylindrique fonctionnant à l’hydrogène et ayant : un alésage
de 82.5 mm, une course de 114.2 mm, une cylindrée de 612.5 cm3 et une vitesse de rotation
constante égale à 600 tours par minute (1CH8.255/11.42). L’étude sus-mentionnée a été
réalisée à l’Université de Ghent en Belgique à titre d’une thèse de doctorat [2].
𝜑
Ou, 𝜑� = �𝜑𝑧 , 𝜑, 𝜑𝑧 - la durée relative de la combustion, la durée atuelle à partir du
début de la combustion et la durée de l’inflammation, presentées en angles de rotation du
vilebrequin; C - constante (par exemple, la constante C = ln (1 − 0.999) = − 6.908 qui
correspond au rendement de la combustion égale à 99.9 %); m - exposant de Wiebe ou indice
du caractère de la combustion - constante qui définit la dynamique de la combustion. Les
constantes 𝜑𝑧 et m dans ce modèle sont choisies à l’aide d’une méthode empirique issue de la
condition qui assure la similitude entre les diagrammes pression vs angle de rotation du
vilebrequin, calculés et mesurés (la loi de la combustion).
Dans ce travail, on a défini l'influence des paramètres constructifs et ceux du régime
du processus de travail sur les lois de la combustion et, par conséquent, sur les paramètres 𝜑𝑧
et m des modèles de Wiebe.
6.1 Abstract
6.2 Introduction
essential material resources. To hydrogen lacks carry the high cost price of its production,
difficulty of storage of hydrogen onboard a vehicle, absence of an infrastructure, etc.
Advantages of hydrogen cause development of the works connected with use of given
fuel in internal combustion engines [1-7].
It is known that one of the most effective tools of research of processes in internal
combustion engine and optimisation of its design is mathematical modelling. By working out
of mathematical models there is some question, not shined widely enough in the technical
literature.
It concerns dynamics of a thermal emission in the hydrogen engine and the influence
on working process of constructive and adjusting parameters of the engine.
CFR engine intended originally for definition of combustible octane number has been
converted for hydrogen use. On the engine have established a hydrogen injection system in
the inlet manifold. The quantity of submitted fuel changed regulation of duration of the valve
lift upon the injector nozzle.
In a cylinder head the piezoelectric sensor cooled by water for pressure measurement
has been established. Also loads, frequency of rotation of a cranked shaft, TDC position, mass
consumption of air and hydrogen, mixture pressure in front of valves and exhaust gases
pressure were measured. In the exhaust gases the concentration of oxygen О2, hydrogen Н2,
nitrogen oxides NO and nitrogen dioxide NO2, carbon monoxide CO, carbon dioxide СО2 and
unburned hydrocarbons СН was defined. Frequency of rotation in all tests was supported by a
constant n = 600 min-1, compression ratio changed in a range from 7 to 9.5.
The experimental data received on engine CFR, have allowed to specify model of
combustion of hydrogen in the four-cycle engine with spark ignition.
Now are known single - two- and multi-zones combustion models. Laws of mass and
energy conservation are put in a basis single-zone models. Quantity of the burned fuel by time
141
and speed of the heat release are set by empirical or semi-empirical dependences (for
example, Wiebe's technique). The co-ordinate of time is introduced.
In multi-zones models of the equation of laws of mass and energy conservation in the
differential form are supplemented with the equations of Nave-Stocks describing turbulent
flow of gas, and, besides time, spatial co-ordinates are entered.
In two-zone models the spatial co-ordinate is the radius of front of a flame.
Combustion process is modelled on the basis of the data about speed of propagation of a
flame in laminar and turbulent environments.
The inconvenient of two- and multi-zone models is the necessity and the complexity
of the use of a considerable quantity of the boundary conditions received as a result of
experimental or theoretical researches. Speed of propagation of a flame, as a rule, is set by
empirical formulas on the basis of experimental data about speed of laminar and turbulent
burning in bombs [2, 4, 5, 8, and 9]. The account of specificity of combustion in internal
combustion engine is made by empirical factors that reduce accuracy of calculations at
researches for other types of engines. The field of application of the given models - detailed
theoretical study of combustion.
In the present work for research and modelling of process of combustion of hydrogen
the thermodynamic model used is Wiebe based. This single-zone semi-empirical model
calculates the rate of mass that is burned and the speed of combustion are put (See annexe
A.I):
m +1
x = 1 − exp(Cϕ ), (6.1)
dx m +1 m m +1
= −C ϕ exp(Cϕ ), (6.2)
dϕ ϕz
𝜑
Where 𝜑� = �𝜑𝑧 , 𝜑, 𝜑𝑧 - accordingly relative time of combustion, current time from the
beginning of combustion and burn duration presented in angles of rotation of a crankshaft;
C - a constant (for example, the constant C = ln (1 − 0.999) = − 6.908 corresponds to a
combustion efficiency of 99.9 %); m - Wiebe exponent or indicator of character of
combustion - a constant defining dynamics of process of combustion. Constants 𝜑𝑧 and m in
the given model are selected by an empirical method from a condition of providing the
similarity between calculated and measured pressure-crank angle diagrams (the combustion
law). As a rule, researchers define dependences 𝜑𝑧 and m from constructive and regime
parameters of working process, and then use the given dependences in optimisation and in
142
other works. In this work, we have to define the influence of constructive and regime
parameters of working process on the law of combustion and, accordingly, parameters 𝜑𝑧 𝑎nd
m Wiebe's models.
On figure 6.2, the influence of burn duration 𝜑𝑧 and the Wiebe exponent m on the form of the
pressure-crank angle diagram, are shown.
Fig. 6.2 Selection of Wiebe exponent m (a) and burn duration 𝜑𝑧 (b) in formulas (1) and (2)
for a regime when n = 600 min-1 and Ne = 1.3 kW. By points the measured diagram is shown.
It is visible that for a regime when n = 600 min-1 and Ne = 1.3 kW, the best similarity
between calculated and measured diagrams is reached at use of 𝜑𝑧 = 22°crank angle and
m = 2.75. However and in this case the inclination of calculated and measured curves of
pressure does not coincide. It is obvious that combustion occurs more intensively in an initial
phase and less intensively in a final phase, that is to say the combustion law, is close to the
triangular form (figure 6.3). We offer a variable Wiebe exponent, changing on linear
dependence from 𝑚1 = 1.5 to 𝑚2 = 4:
m= (m2 − m1 ) ⋅ ϕ + m1 (6.3)
143
dx 1 m
−C ϕ exp(Cϕ ) ( m + 1) ⋅ ϕ + ϕ ln(ϕ ) ⋅ ( m2 − m1 ) .
m +1 m m +1
= (6.4)
dϕ ϕz
Fig. 6.3 Combustion speed (а) and pressure- crank angle diagram (b) with Wiebe exponent m.
Burn duration 𝜑𝑧 = 22 °crank angle. By points the CFR engine experimental diagram is
shown. n = 600 min-1, Ne =1.3 kW.
This figure shows, that the use of a variable Wiebe exponent m provide a good
similarity between measured and calculated pressure-crank angle diagrams for the regime
with n = 600 min-1 and Ne = 1.3 kW.
144
Fig. 6.4 Influence of a spark advance θad (a), compression ratio ε (b) and air-fuel excess ratio
α (c,d) on measured and calculated pressure-crank angle diagrams. By points experimental
diagrams are shown
Figure 6.4 shows the influence of the spark advance θad, the compression ratio ε and
the air–fuel excess ratio α on measured and calculated pressure-crank angle diagrams. At the
time of the modelling of Wiebe exponent, we put linear dependence (3), where coefficients
𝑚1 and 𝑚2 don’t change and we take the combustion duration 𝜑𝑧 as a constant equal to 22°
crank angle (figure 6.4. а-c).
We observe that the variation of the spark advance and the compression ratio of the
engine, don’t have any influence on the dynamics and on the hydrogen burn duration. In the
same time, the influence of the air-fuel excess ratio on the burn duration is considerable
(figure 6.4. c). This fact confirms measurements. We found an empirical relation between 𝜑𝑧
and α: for the regime with n = 600 min-1 and Ne = 1.3 kW:
145
ϕz
= (12.652 ⋅ α 2
− 14.273 ⋅ α + 9.1265 ) ⋅ψ , (6.4)
where ψ - coefficient, defining the influence of the turbulence and other factors on
combustion duration.
The use of this relation allows a good similarity between measurements and
calculations (figure 6.4. d). The relation (6.4) allows an acceptable similarity between
calculated and measured results for CFR engine only for regimes with n = 600 min-1.
Obviously, for the regimes with wide crank angle, the burn duration decreases because the in-
cylinder mixture turbulence increases. Moreover, for engines with other dimension and a
different configuration of the inlet manifold, the in-cylinder turbulence will change.
Consequently, the coefficient ψ must be reduced or increased.
6.5 Conclusion
6.6 References
[1] Mishenko A.I. Application of hydrogen for automobile engines / A.I.Mishenko. - Kiev:
Sciences. Dumka, 1984. – 143 p.
[3] Kim J.M. Performance characteristics of hydrogen fueled engine with the direct
injection and spark ignition system / J.M. Kim, Y.T. Kim, J.T. Lee, S.Y. Lee//SAE. – paper
N° 95249. – 1995.
146
[4] Polasek M. Application of advanced simulation methods and their combination with
experiments to modelling of hydrogen fueled engine emission potentials / M. Polasek, J.
Macek, M.Takats, O. Vıtek//SAE, paper N° 2002-01-0373. – 2002.
[6] Sierens R. Experimental study of a hydrogen fuelled engine. Proceedings ASME Spring
Engine Technology Conference / R Sierens., S. Verhelst//ICE 34-3. – Paper N °2000-ICE-
285. – pp. 9-16. – San Antonio. – 2000.
[7] Kobayashi H. and Kawazoe H. Flame instability effects on the smallest wrinkling scale
and burning velocity of high pressure turbulent premixed flames. 28th Symposium
(International) on Combustion, pages 375–382, 2000.
[8] Lipatnikov A.N. Turbulent flame speed and thickness: phenomenology, evaluation, and
application in multidimensional simulations / A.N. Lipatnikov, J. Chomiak//Progress in
Energy and Combustion Science. – 28:1–74.–2002.
[9] Keck J.C. Early flame development and burning rates in spark ignition engines and their
cyclic variability/J.C. Keck, Heywood J.B., and Noske G. //SAE. – paper N°870164. – 1987.
147
Chapitre 7
Ce chapitre a fait l’objet d’une publication dans: Journal of Road Transport Scientific
Papers, Kharkov Ukraine - 2011- N°43 - 278 p- pp. 152-159 - ISSN 2219-8342. Il a été
presenté aussi en langue russe à la conférence scientifique internationale: Problems and
prospects of the automotive and road transport, qui a eu lieu à l’Université de l’Automobile et
des Ponts et Chaussées de Kharkiv Ukraine le 3-4 Novembre 2011. La version originale de
l’article en langue russe se trouve dans l’annexe A.
Ce chapitre présente les travaux effectués pour l’article IV, qui porte sur le choix des
paramètres de fonctionnement d’un moteur à petite cylindrée fonctionnant avec un mélange
de gaz naturel et d’hydrogène. Afin de réaliser cette étude, deux approches ont été utilisées:
une analytique et une autre expérimentale.
À titre du modèle de calcul, on a utilisé le modèle de Wiebe avec précisions qui
consistent à considérer que l’indice de la combustion de Wiebe m n’est pas une constante,
mais varie au cours de la combustion.
Puisque les matières toxiques se forment principalement dans la flamme et dans les
produits de combustion et l'apparition de la détonation dans le moteur dépend directement de
la température du mélange imbrûlé lors de la combustion, le modèle précisé de Wiebe était
complété avec un modèle à deux zones du calcul de la combustion.
Pour l’éstimation de la concentration des matières toxiques dans le cylindre du moteur
sur la base des températures du mélange brûlé, issues du calcul mené avec le modèle à deux
zones, on a utilisé la méthode du calcul de la composition d'équilibre des produits de
combustion avec précision du contenu de NOx selon l'équation cinétique, proposée par
Zvonov.
148
7.1 Résumé
On expose les résultats d'étude du moteur à petite cylindrée à gaz fonctionnant avec le
mélange du gaz naturel et d'hydrogène. Pour ce type de moteurs, on précise le modèle
mathématique de la combustion avec l’indice variable de combustion de Wiebe.
7.2 Introduction
Plusieurs auteurs indiquent qu'avec l'épuisement des stocks du pétrole, le gaz naturel
deviendra le combustible le plus répandu pour les moteurs à combustion interne [1, 5, 14
etc.]. Pour améliorer les propriétés écologiques de ce combustible, on doit le mélanger avec
l'hydrogène [2].
Par exemple, la voie la plus efficace pour la réduction de la toxicité du moteur à gaz
est l'utilisation du concept de «la combustion pauvre». De plus, la combinaison de l'hydrogène
avec le gaz naturel permet d’éliminer l’inconvénient principal - la combustion incomplète du
mélange gazeux et le préallumage, ainsi elle permet d'augmenter la limite inférieure de
l'inflammation du mélange selon [3]. En outre la conversion du moteur fonctionnant avec le
combustible [liquide] à base de pétrole, avec le gaz naturel comprimé permet d'augmenter la
durée de la combustion qui, à son tour, mène à l'augmentation des émissions des
hydrocarbures imbrûlés CH dans les gaz d’échappement. L’ajout de l'hydrogène dans le gaz
naturel augmente la vitesse de la flamme, en réduisant la durée de la combustion et en
réduisant ainsi les émissions des CH [4, 5].
Les études réalisées dans [6, 7], montrent qu'une petite augmentation de la
concentration de l'hydrogène dans le combustible à gaz contribue à l'augmentation de la
150
L'analyse des publications a démontré qu'il serait opportun d'utiliser le mélange du gaz
naturel avec l’hydrogène (mélange de carburants) en tant que carburant pour moteur. Ainsi, le
but de l'étude est la sélection des paramètres optimaux du fonctionnement du moteur à petite
cylindrée 4CH7.5/7.35, fonctionnant avec le mélange de carburants.
Pour exécuter les tests sur le moteur 4CH7.5/7.35, fonctionnant avec le mélange du
gaz naturel comprimé et d'hydrogène, on utilisait un banc d’essai décrit dans [10]. Afin
d’assurer le fonctionnement du banc d’essai avec le mélange du gaz naturel et d'hydrogène,
on utilisait le système d’alimentation en hydrogène décrit dans [11]. Ce système assure
l’alimentation indépendante en trois sortes de combustible (l’essence, le gaz naturel et
l'hydrogène), dans différentes proportions. En outre, pour l'alimentation en essence, on
utilisait le circuit électronique d’alimentation du combustible du moteur MeMz-307. La vue
générale du système d’alimentation en gaz naturel et en hydrogène, installé sur le banc d’essai
est présentée sur la figure 7.1.
Le taux de compression sur le banc d’essai pouvait être modifié à l'aide de trois ensembles de
pistons assurant assurant les valeurs suivantes de ε : ε = 9.8; ε = 12 et ε = 14.
θ + 18 ηv − 0.25 ε + 1 n + 500
m=
var i 10.639 ⋅ ϕi ⋅ (α + 0.00025 ) ⋅ × ⋅ ⋅ −
40 0.8 11 5000 , (7.1)
−2.3 ⋅ψ − 28.025 ⋅ ϕi2 + 98.045 ⋅ ϕi3 − 156.86 ⋅ ϕi4 + 86.88 ⋅ ϕi5
où, φ
� i - angle relatif de la combustion, φ
� i = 0 … 1; n - la fréquence de rotation du vilebrequin
du moteur, min-1; α - coefficient d’excès d'air; 𝜂𝜈 - coefficient de remplissage; θ - angle
d’avance à l’allumage, ° tours-vilebrequin jusqu'au point mort haut; ε - taux de compression;
ψ - volume d'hydrogène ajouté au gaz naturel, ψ = 0 … 1.
La relation entre la durée de la combustion et la part de l'hydrogène dans le mélange
combustible ψ sur le régime établi est exposée sur la figure 7.2.
152
Fig. 7.2 Relation entre la durée de la combustion et la part d'hydrogène dans le mélange
combustible (les paramètres du régime : n = 1000 min-1; ε = 12; α = 1; ηv = 0.7;
θ = 22° tours-vilebrequin jusqu'au point mort haut)
Pour le contrôle du modèle de Wiebe avec les équations obtenues, on a réalisé des
tests supplémentaires sur 10 régimes. En comparant les valeurs de pression indiquée pi
obtenues suite aux expérimentations avec celles calculées selon la méthode proposée, on
obtient l’erreur de calcul < 5%.
Puisque les matières toxiques se forment principalement dans la flamme et dans les
produits de combustion et l'apparition de la détonation dans le moteur dépend directement de
la température du mélange imbrûlé lors de la combustion, le modèle précisé de Wiebe était
complété avec le modèle à deux zones du calcul de la combustion, décrit dans [15].
Pour la modélisation de la concentration des matières toxiques dans le cylindre du
moteur sur la base des températures du mélange brûlé, issues du calcul mené avec le modèle à
deux zones, on a utilisé la méthode du calcul de la composition d'équilibre des produits de
combustion avec précision du contenu de NOx selon l'équation cinétique, proposée par
Zvonov. Sa description est exposée dans [12].
153
( hc′ − h ) ⋅ 1 − x ⋅ ε − 1
=
Kd ( )( )
Hu (7.3)
Tableau 7.1. Intervalle de variation des paramètres du moteur lors de la résolution des
problèmes d'optimisation des paramètres régulateurs du moteur
α θ n ψ Φpap
-1
– °RV jusqu’au PMH min - %
0
1 5 800 0
Fig. 7.3 Influence d’ajout d'hydrogène sur l'énergie d'inflammation du mélange gaz naturel et
air pour divers coefficients d’excès d'air
7.5.1 Cartographies
En vertu des études réalisées sont reçues les cartographies pour les systèmes de
gestion du moteur, exposées sur les figures 7.4, 7.5 et 7.6.
157
Fig. 7.8 Variation de la concentration des composants toxiques selon la CPCM 4CH7.5/7.35
_______
- essence;
_._._
- gaz naturel;
_ _ _ _
- mélange combustible.
La toxicité du moteur était estimée à l'aide du cycle à 13 modes ESC. Les résultats des
tests du moteur 4CH7.5/7.35, alimenté par un mélange combustible, avec l'utilisation de ce
cycle sont exposés dans le tableau 7.2.
Tableau 7.2 Résultats des tests réalisés sur le moteur selon le cycle ESC
№ δз M кр n g NOx g CO g CH
g/kWh
g/kWh
g/kWh
Min-1
%
–
1 0.15 0 nXX – – –
2 0.08 100 nA 6.1 0.56 0.16
3 0.1 50 nB 2.5 0.51 0.08
4 0.1 75 nB 4.9 0.48 0.11
5 0.05 50 nA 2.7 0.54 0.12
6 0.05 75 nA 4.2 0.56 0.14
7 0.05 25 nA 0.9 0.32 0.11
8 0.09 100 nB 6.8 0.53 0.13
9 0.1 25 nB 2.1 0.25 0.1
10 0.08 100 nC 7.2 0.5 0.14
11 0.05 25 nC 1.7 0.28 0.05
nC.
13
12 0.05 75 5.1 0.56 0.04
=g ∑(g ⋅δ ).
i =1
i зi
Les émissions moyennes des composants toxiques, g/(kWh), sont définis comme la
somme de la multiplication des émissions des composants toxiques sur treize modes du cycle
sur les coefficients de la valeur de ces régimes.
13
=g ∑ ( gi ⋅ δ зi ) (7.4)
i =1
7.6 Conclusions
3. À la base du moteur à petite cylindrée à gaz 4CH7.5/7.35, un banc d’essai a été conçu,
permettant de réaliser des tests sur trois sortes de combustible dans n'importe quelles
proportions – l’essence, l'hydrogène et le gaz naturel.
7.7 Références
[1] Kirillov N.G. Analyse de la perspective des divers combustibles alternatifs pour moteur: le
gaz naturel liquéfié - le combustible moteur du ХХІ-ème siècle / N.G.Kirill,
A.N.Lazarev//Dvigatelestroenie. - 2010. - № 1. – p. 26-33.
[2] Watson H., Milkins E. Hydrogen and methane - automotive fuels of the future? / SAE,
Australia 1975 3/4, p. 17-19.
[3] Jacob K. Performance of Lean-Burn Natural Gas Fueled Engine on Specific Fuel
Consumption, Power Capacity and Emissions - Calgary, CA: University of Calgary, 1990 - 18
p. - (Preprint / University of Calgary: SAE № 901495).
[4] Karim G.A., L'Al-Alusi Y.H. Some Considerations of Cyclic Variations in Spark Ignition
Engines Fueled with Gaseous Fuels - Calgary, CA: University of Calgary, 1984. - 12 p. -
(Preprint / University of Calgary: SAE № 840232).
[5] Michael R.S., Mirza J.Y., Zafer D., Matthew N.S. The Effect of Hydrogen Addition on
Natural Gas Engine Operation - Calgary, CA: University of Calgary, 1993. – 16 p. - (Preprint
/ University of Calgary: SAE № 932775).
[6] Karim G.A. Comparative Studies of Methane and Propane as Fuels for Spark Ignition and
Compression Ignition Engines / SAE Transactions. - 1984. - № 92. p. 3677-3688.
[8] Michenko A.I. Utilisation de l'hydrogène pour les moteurs automobiles - Kiev : Sciences,
doumka, 1984. - 143 p.
[12] Kabanov A.N. Réduction des émissions des substances nuisibles avec les gaz
d’échappement des moteurs Diesels par voie de leur conversion aux moteurs à gaz: Thèse
candidat des sciences techniques: 05.05.03 - KH, 2007. – 206 p.
[14] Attar A.A. Optimization and Knock Modeling of a Gas Fueled Spark Ignition Engine:
PhD thesis: Mechanical Engineering / Calgary, Canada, 1997. - 248 p.
Chapitre 8
Ce chapitre présente les travaux effectués pour l’article V, qui porte sur le
développement d'un nouveau concept de moteur à combustion externe. L’objectif de ce travail
était de déceler les avantages de ce type de moteur sur le plan des paramètres tels que: le
rendement, la puissance, la consommation spécifique du carburant, la pression à l’intérieur de
la chambre de combustion, le coefficient d’excès d’air et la température des produits de
combustion.
Pour ce faire, on a développé un modèle qui permet de simuler l’influence de la
pression, de la température des produits de combustion et de la variation du coefficient
d’excès d’air dans la chambre de combustion sur les performances du cycle réel du moteur à
combustion externe. On a simulé aussi la variation des paramètres de la charge en fonction de
l'angle de rotation du vilebrequin et la pression à l’intérieur de la chambre de combustion.
Cette simulation a été réalisée tout en s’appuyant sur certaines données résultant d'essais sur
un moteur à combustion interne.
L’objectif principal de cette simulation de moteur à combustion externe est de prédire les
performances de ce type de moteur afin de les comparer à celles d’un moteur à combustion
interne, et déduire par la suite s’il sera profitable de l’utiliser dans le domaine automobile.
Afin d’effectuer une telle simulation, il était nécessaire de développer un modèle pour estimer
l'influence de la pression 𝑝𝑠 et de la température 𝑇𝑠 des produits de combustion dans la
chambre de combustion sur les performances du cycle réel du moteur à combustion externe en
utilisant certaines des données résultant d'essais sur un moteur à combustion interne (variation
166
8.1 Abstract
8.2 Introduction
The principal schematic view of this heat engine on the basis of external combustion
engine is presented in Figure 8.1. The air from the tank 1, under a pressure of 30-50 MPa and
at ambient temperature, flows via main 2 to the reducer 3, then by main 4 to the surge tank 5
in which the pressure is kept at 2-5 MPa, as monitored by the pressure sensor 6. From the
surge tank, the air by drain 7 arrives in the combustion chamber 8 which is equipped with a
fuel supply system 9 (e.g.: an injector) and an ignition system 10 of the air-fuel mixture (e.g.:
a spark plug). The combustion chamber 8 is connected by channels 11 and 12 to the inlet
valves 13 and 14 of the cylinders 15 and 16. The valves 13 and 14 are actuated with a fast
action device (e.g.: solenoid valves). The products of combustion from the zone located above
the pistons of cylinders 15 and 16 are evacuated through the exhaust valves 19 and 20 and the
exhaust mains 21 and 22 to the exhaust manifold 23 and then the atmosphere (this happens
normally during the displacement of the piston towards top dead center (TDC)). The
temperature of the products of combustion at the end of the exhaust manifold is controlled by
the sensor 24. The engine power, for whatever external load, is established using the variation
168
of the quantity of fuel injected into the combustion chamber 8, the duration of opening of the
inlet valves 13 and 14, and the air pressure in the surge tank.
The estimate of the influence of the pressure 𝑝𝑠 and temperature 𝑇𝑠 of the products of
combustion in the combustion chamber on the performances of the real cycle of the external
combustion engine is carried out using analytical calculations together with the use of certain
data resulting from tests on an internal combustion engine (variation of the coefficients of
discharge of the valves according to the height of the lifting of the valves, of the coefficients
of heat transfer according to the dimensional parameters of the engine, of the parameters of
the operating cycles, of the values of the temperatures of the wall surface of the zone located
above the pistons, and of the temperatures of the combustion chamber). Since the volume of
the surge tank is much larger than the volume of the compression chamber, into which arrive
the products of combustion from the combustion chamber, the pressure in the combustion
chamber is taken as constant. The temperature of the products of combustion is given by the
equation for the thermal balance for the combustion chamber:
where, 𝑑𝑀 - variation of the mass of the charge in the zone located above the piston
conditioned by the inlet dMs and the exit 𝑑𝑀𝑏 of the charge in an interval of the calculation
time Δτ; 𝑑𝑄 - heat supplied to the load or removed from the load (in an interval of the
calculation time Δτ); 𝑝, 𝑉, 𝑇, 𝑀 - charge parameters at the beginning of the interval of the
calculation time.
Figure 8.2 shows the results of the calculation of the charge parameters variations in the zone
above the piston as a function of crank angle rotation φ. On figure 8. 3 - the variation of charge
pressure in the zone above the piston is presented depending on its volume (P-V diagram). The
calculations are performed for a 4-cylinder automobile engine with a bore of 88 mm, a piston stroke of
82 mm, a crankshaft rotational speed n = 5000 min-1, and with the ratio of the maximum value of the
effective square of flow areas of outlet valves to the piston square 2𝜇𝐵 𝑓𝐵 /𝐹𝑛 = 0.25 and the inlet
valve to piston 𝜇𝑠 𝑓𝑠 /𝐹𝑛 = 0.025. At the given value of relation 𝜇𝑠 𝑓𝑠 /𝐹𝑛 , the duration of opening of the
inlet valve is approximately 20°crankshaft rotation, pressure of products of combustion in the
combustion chamber 𝑝𝑠 = 5𝑀𝑃, temperature 𝑇𝑠 ≈ 900𝐾, the indicated power will make
approximately 𝑁𝑖 ≈ 95.6𝑘𝑊. The temperature of products of combustion in the zone above the
piston in the beginning of opening of outlet valves (φ = 30° crankshaft rotation before BDC) does not
exceed 320 K, i.e. the heat brought to a charge in the combustion chamber is used sufficiently effectively
(figure 8.2).
From the diagram Pi-Vi (figure 8.3), are defined:
𝐿𝑖 = 𝐿𝑠 + 𝐿𝑝 − (𝐿𝑏 + 𝐿𝐶 ), (8.5)
where, 𝐿𝑠 - work on the section of admission of combustion products in the zone located
above the piston; 𝐿𝑝 - work of expansion of combustion products in the zone located above
the piston; 𝐿𝑏 - work spent for removal of combustion products from the zone located above
171
the piston; 𝐿𝐶 - work of compression of combustion products remained in the zone located
above the piston;
• indicated power, kW:
𝑁𝑖 = 𝑘 ∙ 𝐿𝑖 ∙ 10−3 , (8.6)
• indicated efficiency of the engine at filling of tanks by compressed air in the fuel
station:
energy compressed air arriving in a combustion chamber, kJ/s; 𝐺𝑎 - air consumption through
the combustion chamber, kg/s; 𝑅 - gas constant for 1 kg of air, kJ/(kg·K); 𝑘𝑚 - mean value of
an adiabatic exponent at expansion of air from pressure 𝑝𝑐 (𝑇 = 𝑇0 ) to pressure 𝑝0 ;
where, 𝐵ℎ - hourly fuel consumption; indicated efficiency of the engine at tanks filling by the
compressed air directly onboard a vehicle:
𝑛−1
where, 𝑁𝑐 = 𝐺𝑎 ∙ 𝑚 ∙ 𝑛/(𝑛 − 1) ∙ 𝑅 ∙ 𝑇0 ∙ �(𝑝𝑘 ⁄𝑝0 ) 𝑛 − 1� ∙ 1/𝜂𝑐 - power of the compressor;
Fig. 8.2 Variation in the charge parameters in the combustion chamber in the zone above the
piston with crank angle
173
Fig. 8.3 Cylinder pressure external combustion engine (n = 5000 min-1, air pressure in the
zone-buffer 𝑝𝑠 = 5 MPa and the air to fuel ratio in the combustion chamber α = 3)
The principal results of the simulation of the dependency of the various engine parameters on
input conditions are shown in the next two figures: in Figures 8.4 as a function of piston pressure at a
fixed α, and in Figure 8.5 as a function of α at a fixed piston pressure.
From figure 8.4: At constant value of the coefficient of air excess in the combustion chamber
(α = 3), the temperature of the products of combustion also does not change (𝑇𝑠 ≈ 1000 K).
The temperature of the products of combustion 𝑇𝑒 at the moment of the beginning of opening
of outlet valves decreases markedly from 823 K at 𝑝𝑠 = 3 МPа to 387 K at pressure 𝑝𝑠 = 5 МPа in
the Figure 8.4). The specific indicated air rate gai decreases (to 5.7 kg/(kW·h)), and the specific
indicated fuel rate gi1 (to 0.102 kg/(kW·h)) decreases also by some 10%. The Indicated power
increases by some 167% from 50.5 kW at 𝑝𝑠 = 3 МPа to 119 kW at 𝑝𝑠 = 6 МPа, as do the two
efficiencies accordingly indicated efficiency 𝜂𝑖1 increases from 0.585 to 0.672.
174
Fig. 8.4 Influence of the pressure of combustion products ps in the zone located above the piston on
the indicated parameters of the external combustion engine at α = 3
175
Fig. 8.5 Influence of the variation of the air to fuel ratio α in the combustion chamber, on the
indicated parameters of the external combustion engine at 𝑝𝑠 = 5 МPа.
8.6 Conclusion
The main conclusions obtained from the present study of external combustion engine
are as follows:
176
In the external combustion engine with a filling of tanks by compressed air, the efficiency of
transformation of chemical energy of fuel into mechanical work of gases essentially (on 60 - 70 %)
surpasses the efficiency of traditional internal combustion engines with a spark ignition.
The external combustion engine, in consequence of the raised values of the air to fuel ratio and
low values of the maximum combustion temperature (𝑇𝑠 <1300K), possesses high ecological
characteristics without the use of additional devices for neutralizing the exhaust gases.
In the external combustion engine, it is possible to use any kinds of fuel, both gaseous and
liquid.
However practical use of the external combustion engine as power plant in vehicles, for
example on buses for service conditions in cities, demands: the conducting of an expensive program
of development of the engine and its systems; (among others: the manufacturing from polymeric
materials of the tanks for compressed air capable to withstand operating pressures up to 50 МPа), and
the creation of filling stations and car repair shops necessary to keep asll this stock rolling.
8.7 References
[1] World energy: Development Perspectives Until 2010/ Translation from English by
MIREK Commission; under Redaction of Y.N. Starshinov. – M.: Energia, 1980. -256 p.
[2] Guy Negre and Cyril Negre. Patent USA N°6305171B1 F01K 7/34. Method and Device
for Additional Thermal Heating for Motor Vehicle Equipped Compressed Air Injection, Publ.
23.10.2001.
[3] Dyatshenko V.G. Internal Combustion Engine Theory - KNADU, 2009. - 500 p.
[4] Kang Ma, Lingen Chen, Fengrui Sun. Optimizations of a Model External Combustion
Engine for Maximum Work Output with Generalized Radioactive Heat Transfer law.
International Journal of Energy and Environment. Volume 2, Issue 4, 2011, pp.723-738.
[5] Yehuda B. Band, Oded Kafri, and Peter Salamon. Optimization of a Model External
Combustion Engine. Journal of Applied Physics, Volume 53, Issue 1, Jan. 1982, pp. 29 - 33.
[6] Song, H. J.; Chen, L. G.; Sun, F. R. Optimisation of a Model External Combustion Engine
with Linear Phenomenological Heat Transfer Law. Journal of the Energy Institute, Volume
82, Number 3, September 2009, pp. 180-183(4).
177
Chapitre 9
Conclusions et Perspectives
Les résultats ainsi obtenus contribueront à une meilleure compréhension des différents
phénomènes physico-chimiques complexes qui ont lieu dans les moteurs à combustion interne
alimentés par des carburants de substitution à l’état pur et en mélange et à justifier l’utilisation
future de ce type de carburants et des moteurs à combustion externe.
Afin d’atteindre nos objectifs, nous avons divisé ce travail en trois parties principales :
9.1 Utilisation des carburants de substitution à l’état pur dans les MCI
Les principaux résultats de cette étude dans le moteur à combustion interne peuvent se
résumer comme suit:
Pour ce qui est des moteurs alimentés par l’hydrogène pur, la courbe de vitesse
d’inflammation est proche de la forme triangulaire. Dans le modèle de combustion de Wiebe,
l’utilisation de l’indice variable du caractère de la combustion m, qui varie selon une ligne,
permet de décrire d’une façon plus précise la dynamique de combustion d’hydrogène. La
variation du taux de compression et de l’angle d’avance à l’injection influence a un degré de
178
Dans cette partie, on a réalisé des tests sur un moteur à gaz naturel à petite cylindrée
avec ajout d’hydrogène, en élaborant des cartographies pour le système de gestion avec
microprocesseur de ce même moteur. On a élaboré aussi une approche analytique basée sur le
modèle de Wiebe.
Nous avons montré que :
Pour réaliser cette étude, on a réalisé des tests sur un nouveau concept de moteur à
combustion externe et on a élaboré un modèle d’estimation de l’influence de la pression et de
la température des produits de combustion à l’intérieur du cylindre sur les performances de ce
179
même moteur.
Les principales conclusions obtenues durant cette étude sont les suivantes :
Le rendement de transformation de l'énergie chimique du combustible en travail
mécanique des gaz dépasse essentiellement (de 60 - 70%) le rendement des moteurs à
combustion interne conventionnels à allumage par étincelle.
Le moteur à combustion externe, à cause des valeurs trop élevées de la richesse et des
faibles valeurs de la température de combustion maximale (𝑇𝑠 < 1300 𝐾), présente des
caractéristiques écologiques élevées sans l'utilisation de dispositifs supplémentaires pour
neutraliser les gaz d'échappement.
Dans le moteur à combustion externe, il est possible d'utiliser tous les types de
carburant, tant gazeux et liquides.
9.4 Perspectives
Les études effectuées au cours de cette recherche ont permis d'avoir une meilleure
compréhension de l'influence de l’utilisation des carburants de substitution à l’état pur et en
mélange, sur les performances des moteurs à combustion interne ainsi que l’utilisation du
concept des moteurs à combustion externe. Ces études sont donc un premier pas important
vers une aide précieuse à l’emploi de nouveaux types de carburants et de moteurs ainsi qu'à
l'amélioration des modèles mathématiques d’optimisation des moteurs fonctionnant avec de
tels combustibles, qui permettront bientôt de réduire le nombre de tests expérimentaux longs
et laborieux. Ces études auront aussi permis d'ouvrir de nouvelles voies de recherche plus
appliquée dans la compréhension de la phase de combustion des moteurs sus-mentionnés et
des paramètres physiques pouvant l’influencer. Cependant, avant d'aboutir à une solution
optimale, d'autres études expérimentales et analytiques devront être entreprises et de
préférence, en suivant les recommandations présentées ci-après.
Les résultats obtenus los de l’analyse de la phase de combustion de l’hydrogène en
utilisant le modèle de Wiebe, vont permettre une description plus précise de la dynamique de
combustion d’hydrogène. Il serait donc intéressant d’envisager le couplage de cette étude avec
d’autres études utilisant différents modèles de combustion.
Il a été aussi démontré que la variation du taux de compression et de l’angle d’avance
à l’injection à un degré moindre, influence la durée et la dynamique de la phase de
combustion plus que la variation du coefficient d’excès d’air. Dans ce contexte, il sera
souhaitable, d’étudier l’effet des autres paramètres sur la phase de combustion, tels que:
l’avance à l’allumage, la richesse, la vitesse de rotation et le couple moteur.
180
Pour ce qui est de la toxicité dans les moteurs à combustion interne, il serait donc
intéressant de comparer la méthode de calcul de la composition d'équilibre des produits de
combustion, proposé par le professeur V.A. Zvonov et complétée par l'équation cinétique de
Zeldovitch avec d’autres méthodes existantes afin de déceler les avantages et les
inconvénients de cette méthode.
Le choix et le fondement des paramètres régulateurs du moteur à gaz suralimenté
converti à partir d’un moteur diesel en utilisant le modèle mathématique du cycle de combustion de
Wiebe avec le calcul précis du paramètre de la combustion m et en prenant en considération la variation
de la pression de suralimentation, s’avèrent très utiles pour des futures études de ce type de
moteurs, du point de vue des émissions de gaz à effets de serre.
Lors de notre étude visant à choisir les paramètres optimaux du fonctionnement d’un
moteur à allumage par étincelle alimenté par un mélange gaz naturel/hydrogène, nous avons
constaté que l’ajout de l'hydrogène au gaz naturel permet d'éliminer les inconvénients
d’utilisation du concept de «la combustion pauvre», tout en gardant son avantage. Les
émissions des gaz polluants ont chuté d’une façon significative: NOx de 72 %; CO de 99 % et
CH de 96 %. On a eu un gain minime au niveau de la consommation du carburant et
malheureusement une perte de puissance. Il serait donc préférable et souhaitable même
d’envisager le couplage de cette étude avec une étude détaillée d’un autre moteur pour une
évaluation technico-économique plus détaillée de cette alternative intéressante pour la
réduction des émissions polluantes et le développement des moteurs propres.
Les moteurs à combustion externe s’avèrent une solution prometteuse dans le domaine
des transports. Avec un rendement et un rapport de mélange élevés, une basse température de
combustion, ce genre de moteurs peut fonctionner avec différents types de carburant, qu’il
soit liquide ou gazeux. La présente étude expérimentale et analytique a été réalisée sur un
moteur fonctionnant seulement à l’essence. On attend que dans un futur proche, d’autres
travaux liés à ce domaine voirent le jour, mais cette fois avec différents types de carburants.
181
ANNEXE A
The work of the Russian scientist Nikolai Semenov in the area of chemical kinetics
and chain chemical reactions was well established and known in the mid-1950s and engine
researchers had already started wondering whether better models, based on the theory of
radical chain reactions, could be constructed to characterize engine combustion processes and
thus contribute towards the design of better engines. The interest was further strengthened
following the award to Semenov, jointly with Sir Cyril Hinshelwood, of the Nobel prize for
chemistry in 1956 for their work in chemical reaction kinetics. Wiebe was one of the first to
attempt linking chain chemical reactions with the fuel reaction rate in internal combustion
engines and his approach was based on the premise that a simple one-step rate equation will
not be adequate to describe complex reacting systems such as those occurring in an internal
combustion engine. On the other hand, developing and solving rate equations which account
for the simultaneous and sequential interdependent chain and chain branching reactions would
be prohibitive. He reasoned that, for engineering purposes, the details of chemical kinetics of
all the reactions could be bypassed and a general macroscopic reaction rate expression could
be developed based on the concept of chain reactions.
In real combustion systems such as in an internal combustion engine, chain reactions
may proceed sequentially and in parallel with the formation of intermediate species
comprising free atoms and radicals. Wiebe called these intermediate species active centres.
Being highly reactive, active centres play a pivotal role in chemical kinetics reactions and are
crucial to the reaction path. For the initiation of the reaction, a certain concentration of active
centres (initial centres) is required and these centres can be generated by heating the air–fuel
mixture or providing electrical charge. During the combustion process a very large number of
active centres will exist in close proximity to the molecules of the main reacting substances.
Combustion produces molecules of final products and more active centres that are capable of
starting new reaction cycles. Excluded from this scheme are chain inhibiting centres resulting
from radicals or free atoms colliding with a third body (any other radical, atom, or molecule
in the system or the system walls). As combustion progresses, concentration of reactants
decreases progressively, causing corresponding decrease in the reaction rate.
In the combustion reaction of hydrogen forming water, a leading role is played by the
elementary process
182
OH+O2=H2O+H (Reaction I)
Reaction III provides the chain carriers H and OH for two more branching elementary
processes (I and II). The three elementary processes (I, II, and III) make-up and elementary
reaction cycle which is continuously repeated as the chain reactions propagate. The highly
reactive radicals OH produced by processes (II) and (III) are the active centres.
Similarly, the oxidation of carbon monoxide can be represented by the following
elementary reaction cycle
O+CO=CO2 (Reaction V)
Fig.1 Schematic diagram of elementary chemical processes, elementary cycles, and effective
events of chain reactions as a function of time (based on the scheme in [2]: Elementary
reaction cycle; Effective reaction event
On the basis of this concept, Wiebe postulated that the incremental change in the
number of molecules of the main reactants dN participating in the effective reaction events in
the time interval t to t + dt is directly proportional to the change in the number of effective
centres dNe – dN = ndNe where n is the constant of proportionality. In differential form
𝑑𝑁 𝑑𝑁𝑒
− 𝑑𝑡 = 𝑛 . (1)
𝑑𝑡
𝜌= 𝑑𝑡
,
𝑁
where N is the instantaneous number of molecules of the initial reactants. Equation (1) can
now be written as
𝑑𝑁
− 𝑁
= 𝑛𝜌𝑑𝑡 (2)
184
from which
𝑁 𝑡
𝑙𝑛 𝑁 = − ∫0 𝑛𝜌𝑑𝑡
0
𝑡
𝑁 = 𝑁0 𝑒 − ∫0 𝑛𝜌𝑑𝑡 (3)
where No is the number of moles of the main reactant in the combustion system at the start of
reactions.
If the burn fraction of the initial reactant at any moment in time t is defined as
𝑁0 −𝑁
𝑥= 𝑁0
(4)
then
𝑡
𝑥 = 1 − 𝑒 − ∫0 𝑛𝜌𝑑𝑡 (5)
𝑑𝑥 𝑡
𝑤= 𝑑𝑡
= 𝑛𝜌𝑒 − ∫0 𝑛𝜌𝑑𝑡 = 𝜌𝑛(1 − 𝑥) (6)
If the relative density of the effective centres can be expressed as 𝜌 = 𝑘𝑡 𝑚 , where k and m are
𝑑𝑥 𝑚+1
𝑑𝑡
= 𝐾𝑡 𝑚 𝑒 −[𝐾⁄(𝑚+1)]𝑡 (8)
𝑚+1
𝑥𝑑 = 1 − 𝑒 −[𝐾⁄(𝑚+1)]𝑡𝑧
from which
𝑚+1
𝑒 −[𝐾⁄(𝑚+1)]𝑡𝑧 = 1 − 𝑥𝑑 (9)
𝑚+1
𝑒 −[𝐾⁄(𝑚+1)]𝑡 = 1−𝑥 (10)
Taking the logarithm of equations (9) and (10) and dividing one by the other yields, we obtain
the Wiebe function for the non-dimensional burn fraction x as function of time t
𝑚+1
𝑥 = 1 − 𝑒 𝐶(𝑡⁄𝑡𝑧 ) (11)
The derivative w (burn rate) of Wiebe function as functions of time t can now be written as
𝑑𝑥 (𝑚+1) 𝑡 𝑚 𝑚+1
𝑤= 𝑑𝑡
= −𝐶 𝑡𝑧
�𝑡 � 𝑒 −𝐶(𝑡⁄𝑡𝑧 ) (12)
𝑧
𝑑𝑥 (𝑚+1) 𝑚+1
𝑤= 𝑑𝑡
= −𝐶 𝜑𝑧
𝜑� 𝑚 𝑒 −𝐶𝜑� (14)
The equation
dQ=dH-Vdp
is also
dq=dh-vdp.
Enthalpy h, is made up of quantities which are functions of state (u, p, v), enthalpy is
also a function of state. Just as internal energy, the enthalpy of a pure substance can be
represented as a function of any two properties, or parameters, of state, for instance, of
pressure p and temperature T:
h=f(p,t). (1)
The mathematical formulation of the first law of thermodynamics for the case where the only
dq=du+pdv, (3)
xdy=d(xy)-ydx, (4)
pdv=d(pv)-vdp, (5)
ANNEXE B
Voici la version originale de l’article qui a fait l’objet d’une publication dans la revue
scientifique American Journal of Vehicle Design, 2013, Vol. 1, N° 2, 21-24.
Cet article fait partie des articles qui figurent dans la présente thèse dans le chapitre 6.
Néanmoins, j’ai decidé de mettre sa version originale en langue russe dans l’annexe B.
А.П. Марченко, д-р техн. наук, А.А. Осетров, канд. техн. наук,
И. Дубей, д-р философии, Маамри Р., докторант
Актуальность проблемы
Цель исследования
Экспериментальное оборудование
m +1
x = 1 − exp(Cϕ ), (1)
dx m +1 m m +1
= −C ϕ exp(Cϕ ), (2)
dϕ ϕz
Рис. 3. Влияние изменения угла опережения зажигания θоп (а), степени сжатия ε (б) и
коэффициента избытка воздуха α (в, г) на экспериментальную и расчетную
индикаторные диаграммы. Точками показаны экспериментальные диаграммы
Видно, что для режима n = 600 мин-1, Ne = 1,3 кВт наилучшее совпадение расчетной и
экспериментальной диаграмм достигается при использовании 𝜑𝑧 = 22 гр. п.к.в. и
m = 2,75. Однако и в этом случае наклон расчетных и экспериментальных кривых
давления не совпадает. Очевидно, что сгорание происходит более интенсивно в
начальной фазе и менее интенсивно в конечной фазе, то есть закон сгорания близок к
треугольной форме (рис. 3). Нами предложен переменный показатель динамики
сгорания, изменяющийся по линейной зависимости от m1 = 1,5 до m2 = 4:
m= (m2 − m1 ) ⋅ ϕ + m1 . (3)
В этом случае скорость сгорания определяется по формуле:
dx 1 m
−C ϕ exp(Cϕ ) ( m + 1) ⋅ ϕ + ϕ ln(ϕ ) ⋅ ( m2 − m1 ) .
m +1 m m +1
=
dϕ ϕz
Выводы
Литература
3. Kim J.M. Performance characteristics of hydrogen fueled engine with the direct
injection and spark ignition system / J.M. Kim, Y.T. Kim, J.T. Lee, S.Y. Lee // SAE. - paper
nr 95249. - 1995.
9. Keck J.C. Early flame development and burning rates in spark ignition engines and
their cyclic variability /J.C. Keck, Heywood J.B., and Noske G. // SAE. - paper nr870164. -
1987.
195
Cet article a fait l’objet d’une publication dans: Eastern European Journal of
Enterprise Technologies, Kharkov Ukraine - 2011- N° 5/6 (53) - 63 p - pp. 37-43 - ISSN
1729-3774. Il a été presenté aussi lors de la conférence scientifico-méthodologique et de
recherche scientifique, section des moteurs à combustion interne et sous-section de
l’amélioration des moteurs à piston à combustion interne, qui a eu lieu à l’Université
Technique de l’Automobile et des Ponts et Chaussées de Moscou Russie le 29 Janvier - 07
Février 2013. Cet article fait partie de la présente thèse dans le chapitre 4 (Traduction en
français) et ici on trouve sa version originale.
УДК 621.43.052
Abstract. The toxic factors calculation method of gas engine has been shown. Comparison of
calculation results with experimental data has been performed.
Key words: calculation method, combustion process, toxic, gas engine, equilibrium
composition.
Введение
Анализ публикаций
Цель исследования
Так как механизм развития цепных реакций в процессе сгорания смесей воздуха и
углеводородов изучен не до конца и в основном представляет собой теоретические
допущения, у разных авторов предполагается разное количество уравнений. Так, в [6]
использовалась система из 155 элементарных реакций для 39 компонентов. Однако,
увеличение числа уравнений увеличивает сложность решения системы уравнений.
Кроме того, многие авторы указывают, что учитывать многие промежуточные
продукты реакций не обязательно [2, 5, 6 и др.].
→ CO 2 + H 2 .
CO + H 2 O ←
1.CO 2 ← → CO+0.5O 2 ;
2.H 2 O ← →
2H +0.5O ;
2
3.H 2 O ←→
OH+0.5H 2 ;
(1)
4.H 2 ← →
2H;
5.O 2 ← → 2O;
6.0.5N 2 +0.5O 2 ←→ NO.
1
PCO ⋅ PO2 2
= K1 , (2)
PCO2
1
PH2 ⋅ PO2 2
= K2 , (3)
PH2O
1
POH ⋅ PH2 2
= K3 , (4)
PH2O
PH 2
= K4 , (5)
PH2
PO 2
= K5 , (6)
P O2
PNO
1 1
= K6 , (7)
PO2 ⋅ PN2
2 2
(
O т + χ ОС Oвозд μ О 2 PCO2 + PCO + A
= ,
)
( )
(8)
С т + χ ОС Cвозд μ С PCO + PCO 2
где A = PH2O + 2 PO2 + POH + PNO + PO ; O т - массовая доля атома кислорода в топливе;
Oвозд – массовая доля атома кислорода в воздухе; С т - массовая доля атома углерода в
H т + χ НС ⋅ H возд μН ⋅ ( B)
= , (9)
(
С т + χ НС ⋅ Cвозд μ С ⋅ PCO + PCO
2
)
где B = 2 ⋅ PH2O + 2 ⋅ PH2 + POH + PH ; H т - массовая доля атома водорода в топливе; H возд -
- массовая доля атома водорода в воздухе; χ НС - действительное соотношение масс
(
N т + χ NC ⋅ N возд μ N ⋅ 2 ⋅ PN2 + PNO
= ,
)
( )
(10)
Ст + χ NC ⋅ Cвозд μ С ⋅ PCO2 + PCO
где N т - массовая доля атома азота в топливе; N возд - массовая доля атома азота в
PO=
2
= P=
PNO O P=
H = 0.
POH (12)
где
μ C O т + χ ОС ⋅ Oвозд
=
W ⋅ . (14)
μ O С т + χ ОС ⋅ Cвозд
1
2 ⋅ PH2O + 2 ⋅ PH2
= W2 , (15)
PCO2 + PCO
где
H т + χ НС ⋅ H возд
W= μC ⋅ . (16)
С т + χ НС ⋅ Cвозд
2
2 ⋅ PN2
= W3 , (17)
PCO2 + PCO
где
μ C N т + χ NC ⋅ N возд
=
W ⋅ . (18)
μ N С т + χ NC ⋅ Cвозд
3
PCO ⋅ PH2O K1
= . (20)
PCO2 ⋅ PH2 K2
Отсюда
=
PH2O PCO2 (W1 − 2 ) + PCO (W1 − 1) ; (22)
𝑃𝐻2 = 0.5𝑊2 �𝑃𝐶𝑂2 + 𝑃𝐶𝑂 � − 𝑃𝐶𝑂2 (𝑊1 − 2) − 𝑃𝐶𝑂 (𝑊1 − 1). (23)
PN2 = (
0.5 ⋅ W3 PCO2 + PCO . ) (24)
�𝑃𝐶𝑂2 + 𝑃𝐶𝑂 � + 𝑃𝐶𝑂2 (𝑊1 − 2) + 𝑃𝐶𝑂 (𝑊1 − 1) + 0.5𝑊2 �𝑃𝐶𝑂2 + 𝑃𝐶𝑂 � − 𝑃𝐶𝑂2 (𝑊1 − 2) −
−𝑃𝐶𝑂 (𝑊1 − 1) + 0.5𝑊2 �𝑃𝐶𝑂2 + 𝑃𝐶𝑂 � − 𝑃𝐶𝑂2 (𝑊1 − 2) − 𝑃𝐶𝑂 (𝑊1 − 1) + 0.5𝑊3 �𝑃𝐶𝑂2 +
+𝑃𝐶𝑂 ) = 𝑃 (25)
𝑃𝐶𝑂2 + 𝑃𝐶𝑂 + 0.5𝑊2 𝑃𝐶𝑂2 + 0.5𝑊2 𝑃𝐶𝑂 + 0.5𝑊3 𝑃𝐶𝑂2 + 0.5𝑊3 𝑃𝐶𝑂 = 𝑃 (26)
( PCO 2
)
+ PCO (1 + 0.5W2 + 0.5W3 ) =
P, (27)
или
P
PCO2 + PCO = . (28)
1 + 0.5W2 + 0.5W3
0.5W3 P
PN 2 = , (29)
1 + 0.5W2 + 0.5W3
𝑃𝐶𝑂 + 0.5𝑊2 𝑃𝐶𝑂 + 0.5𝑊3 𝑃𝐶𝑂 = 𝑃 − 𝑃𝐶𝑂2 − 0.5𝑊2 𝑃𝐶𝑂2 − 0.5𝑊3 𝑃𝐶𝑂2 (31)
2P
=PCO − PCO2 . (33)
2 + W2 + W3
2∙𝑃
𝑃𝐻2 𝑂 = 𝑃𝐶𝑂2 ∙ (𝑊1 − 2) + �2+𝑊 +𝑊 − 𝑃𝐶𝑂2 � (𝑊1 − 1) (34)
2 3
2𝑊 𝑃 2𝑃
𝑃𝐻2 𝑂 = 𝑃𝐶𝑂2 𝑊1 − 2𝑃𝐶𝑂2 + 2+𝑊 1 𝑊 − 𝑊1 𝑃𝐶𝑂2 − 2+𝑊 + 𝑃𝐶𝑂2 (35)
2+ 3 2+ 𝑊3
2 P (W1 − 1)
=PH 2O − PCO2 . (36)
2 + W2 + W3
203
2𝑃 𝑊𝑃
𝑃𝐻2 = 0.5𝑊2 2+𝑊 +𝑊 − 𝑊1 𝑃𝐶𝑂2 + 2𝑃𝐶𝑂2 − 2+𝑊1+𝑊 + 𝑊1 𝑃𝐶𝑂2 + 𝐶. (38)
2 3 2 3
W2 P − 2W1P + 2 P
=PH 2 + PCO2 , (39)
2 + W2 + W3
P ( 2 − 2W1 + W2 )
=
PH 2 PCO2 + . (40)
2 + W2 + W3
K 2 PCO PH 2O − K1PCO2 PH 2 =
0. (42)
2 P (W1 − 1)
= M2 . (44)
2 + W2 + W3
P ( 2 − 2W1 + W2 )
= M3 . (45)
2 + W2 + W3
Перепишем уравнение (42) так, чтобы все парциальные давления в нем были выражены
через парциальное давление СО2
2 2
𝐾2 𝑀1 𝑀2 − 𝐾2 𝑀1 𝑃𝐶𝑂2 − 𝐾2 𝑀2 𝑃𝐶𝑂2 + 𝐾2 𝑃𝐶𝑂2
− 𝐾1 𝑃𝐶𝑂2
− 𝐾1 𝑃𝐶𝑂2 𝑀3 = 0 (47)
.
2 (𝐾
𝑃𝐶𝑂2 2 − 𝐾1 ) + 𝐾2 𝑀1 𝑀2 − 𝑃𝐶𝑂2 (𝐾2 𝑀1 + 𝐾2 𝑀2 + 𝐾1 𝑀3 ) = 0 (48)
K 2 M1 + K 2 M 2 + K1M 3 − D
PCO2 = . (50)
2 ( K 2 − K1 )
Парциальное давление О2
2
K1 ⋅ PCO2
PO2 = . (51)
PCO
K3 ⋅ PH 2O
POH = . (52)
PH 2
=
PH K 4 ⋅ PH 2 . (53)
=
PO K5 ⋅ PO2 . (54)
Парциальное давление NO
=
PNO K 6 PN2 ⋅ PO2 . (55)
0. O 2 ←→ O + O;
→ NO + N;
k1п
1. N 2 + O ←
k1р
k 2п (56)
2. N + O 2 ← → NO + O;
k 2р
k 4п
→
3. N 2 + O 2 ← 2 ⋅ NO.
k 4р
В системе химических реакций (56) k1п, k1р, k2п. k2р. k4п. k4р - константы скоростей
реакций (1), (2) и (3) соответственно.
drNO p ⋅ rN 2 2
rNO k1п ⋅ rO k 4п ⋅ rO2
= 1 − × + , (57)
dϕ 254.4nTп.с. K 4(NO) rO2 rN 2
1 + k1p ⋅ rNO 2
k 2п ⋅ rO2
K 4( NO ) 21.12ехр(−181000 / ( RTпс )) ;
= (61)
207
НАЧАЛО
Да
Загрузка j=1
исходных данных
Нет
i = 0, j = 0
Расчёт оксидов
азота
i=i+1
Нет
j=1
Расчёт равновес.
j=j+1 Да
сост. 1-е прибл.
j=0
Расчёт равновес.
сост. 2-е прибл.
Вывод
Нет результатов
∆<∆доп
расчёта
Нет
i = isgor
Да
КОНЕЦ
WCO
%
0,04
0,03
0,02
0,01
0
WNOx
млн-1
1200
1000
800
600
400
200
Выводы
Литература
6. Bade Shrestha S.O. A Predictive Model for Gas Fueled Spark Ignition Engine
Applications / Bade Shrestha S.O., Karim G.A. - Calgary, CA: Univercity of Calgary, 1999. -
18 p. - (Preprint / Univercity of Calgary: SAE 1999-01-3482).
210
Cet article a fait l’objet d’une publication dans: Journal of The National Polytechnical
University (KPI), Kharkov, Ukraine - 2011- N°54: New Solutions in Modern Technologies -
156 p - pp. 18-26 - ISSN 2079-5459. Il a été presenté aussi lors de la conférence scientifico-
méthodologique et de recherche scientifique, section des moteurs à combustion interne et
sous-section de l’amélioration des moteurs à piston à combustion interne, qui a eu lieu à
l’Université Technique de l’Automobile et des Ponts et Chaussées de Moscou Russie le 29
Janvier - 07 Février 2013. Cet article ne fait pas partie de ceux qui figurent dans la présente
thèse. Néanmoins, j’ai decidé de mettre sa version originale dans l’annexe B.
Results of investigation of gas engine with supercharging converted from diesel have
been given. For this type of engines math model of combustion process with variable Viebe
combustion factor has been specified.
Key words: toxicity, supercharging, gas engine.
1. Введение
3. Цель исследования
4. Экспериментальный стенд
5. Математическая модель
=
mvar i 10.639 ⋅ ϕi ⋅ A ⋅ B ⋅ C ⋅ E ⋅ F + J , (1)
где
θ + 18
A= , (2)
40
B = α+0.00025 , (3)
0.005 ⋅ ηv + 0.005
C= , (4)
0.01
n + 1100
E= , (5)
3000
π k + 0.001
F= , (6)
0.95
( hc′ − h ) ⋅
=Kd (1 − x ) ⋅ ( ε − 1) , (9)
Hu
где hc′ - удельная энтальпия несгоревшей смеси в момент начала сгорания, кДж/кг; h -
удельная энтальпия несгоревшей смеси в любой момент времени, кДж/кг; Hu - низшая
теплота сгорания, кДж/кг, x - доля выгоревшего топлива; ε - степень сжатия.
216
α θ n πk φдр
- град. пкв до ВМТ хв-1 - %
1 5 800 1 0
1,5 40 2100 1,45 100
низком значении Pe .
Таким образом, на минимальных оборотах холостого хода n = 800 мин-1 значение α
составляет α = 1,05.
Такой подход позволяет снизить выбросы NOx на средних режимах на 60…70 % по
сравнению с безнаддувным вариантом, работающим при малых значениях α.
219
…10 % на всех режимах ВСХД. Это позволяет снизить выбросы оксидов азота NOx на
режимах ВСХД на 8…15 %. Выбросы СО и СН в обоих двигателях находятся
практически на одинаковом уровне.
На частичных режимах давление наддува снижается, параллельно увеличивается α.
Кроме того, осуществляется смешанное регулирование смеси на частичных режимах.
Благодаря этому на частичных режимах выбросы NOx снижаются на 19…36 %. На
режимах холостого хода разница в выбросах NOx практически отсутствует.
Результаты испытаний газового двигателя 6ГЧН13/14 с использованием 13-
режимного цикла ESC показали следующие средние эксплуатационные выбросы
нормируемых токсичных компонентов, г/(кВт⋅ч): g NOx = 2,36; g CO = 1,03; g CH = 1,07.
При этом средние эксплуатационные выбросы нормируемых токсичных компонентов
безнаддувного газового двигателя 6ГЧ13/14 составили, г/(кВт⋅ч) [6]: g NOx = 3,28;
g CO = 1,15; g CH = 1,31.
Снижение средних эксплуатационных выбросов при установке регулируемого
наддува на газовый двигатель 6ГЧ13/14 по сравнению с безнаддувным вариантом
составило: по NOx - 28 %, по СО - 10 %, по СН - 18 %.
Выводы
Литература
3. Cox G.B. Development of a Direct-Injected Natural Gas Engine System for Heavy-
Duty Vehicles: Final Report [Текст] / G.B. Cox, K.A. DelVecchio, W.J. Hays - Caterpillar
Inc. - 2002. - 67 p.
4. Nylund N-O. Pathways For Natural Gas Into Advanced Vehicles [Текст] / N-O
Nylund, J. Laurikko, M. Ikonen. - Brussel: IANGV. - 2002. - 105 p.
8. Attar A.A. Optimization and Knock Modelling of a Gas Fueled Spark Ignition Engine
[Текст]: PhD thesis: Mechanical Engineering / A.A. Attar. Calgary, Canada, 1997. - 248 p.
Cet article a fait l’objet d’une publication dans: Journal of The National Polytechnical
University (KPI), Kharkov, Ukraine - 2011- N°43 - 155 p - pp. 61-69 - ISSN 2079-5459.
Cet article ne fait pas partie de ceux qui figurent dans la présente thèse. Néanmoins, j’ai
decidé de mettre sa version originale dans l’annexe B.
Введение
Концепция «бедного горения» (то есть горения смеси в цилиндре при высоких α)
является одним из наиболее эффективных способов снижения выбросов NOx
двигателем с искровым зажиганием. При этом применениие добавки водорода к
газовому топливу позволяет устранить основные недостаток данной концепции –
неполное сгорание рабочей смеси и пропуски вспышек [3]. Кроме того, конвертация
двигателя, работающего на жидком нефтяном топливе, на сжатый природный газ
приводит к увеличению продолжительности сгорания, что, в свою очередь, приводит к
увеличению выбросов несгоревших углеводородов СН с отработавшими газами.
Добавка водорода к газовому топливу увеличивает скорость пламени, сокращая
продолжительность сгорания и снижая таким образом выбросы СН [4, 5].
Следует также отметить, что низкая плотность водорода при нормальных условиях
сводит на нет эту прибавку мощности при содержании водорода в газовом топливе
более 10 % по объёму [6, 7]. Кроме того, при α < 2 добавка водорода снижает
детонационную стойкость газового топлива [8].
227
Цель исследования
ϕZ
град. пкв
50
40
30
20
0 0,2 0,4 0,6 0,8 ψ
( hc′ − h )
=
Kd ⋅ (1 − x ) ⋅ ( ε − 1) , (1)
Hu
где hc′ - удельная энтальпия несгоревшей смеси в момент начала сгорания, кДж/кг; h -
удельная энтальпия несгоревшей смеси в любой момент времени, кДж/кг; Hu - низшая
теплота сгорания, кДж/кг, x - доля выгоревшего топлива; ε - степень сжатия.
Расчёты показали, что для ψ = 0,1 (при α = 1) предельной степенью сжатия является
ε = 12, что было подтверждено экспериментально. Как следствие, выбор оптимальных
параметров работы двигателя проводился для степени сжатия ε = 12.
230
Экспериментальный стенд
Степень сжатия на стенде могла изменяться при помощи трёх комплектов поршней,
обеспечивающих следующие значения ε: ε = 9,8; ε = 12 и ε = 14.
Характеристические карты
α
2,0
1,8
1,6
1,4
1,2
θ
град. пкв
до ВМТ
40
30
20
10
1,5
1,4
1,3
1,2
1,1
1,0
0,9
v
v
0,82
0,80
0,78
0,76
Me
Нм
110
Me
100
90
ge
г
кВт год
380
360
ge
340
320
300
280
Ne
кВт
50
Ne
40
30
20
10
1000 1500 2000 2500 3000 3500 4000 4500 5000 5500 n, хв-1
_______
- бензин;
__ . __ . __
- природный газ;
_____
- смесевое топливо
Рис. 5. Внешняя скоростная характеристика двигателя 4ГЧ7,5/7,35
234
NOx
млн-1
NOx
2500
2000
1500
1000
500
CH
млн-1
180
CH
160
120
CO
%
3,0
2,5
2,0
CO
1,5
1,0
0,5
1000 1500 2000 2500 3000 3500 4000 4500 5000 5500 n, хв-1
_______
- бензин;
__ . __ . __
- природный газ;
_____
- смесевое топливо.
Рис. 6. Внешняя скоростная характеристика двигателя 4ГЧ7,5/7,35
Как видно из рис. 5, при переводе двигателя на смесевое топливо с ψ = 0,1 практически
удалось избежать потерь мощности на номинальном режиме по сравнению с работой
при ψ = 0 (Ne = 50 кВт). На режиме максимального крутящего момента Memax
(n = 3500 мин-1) снижение мощности при переходе с природного газа на смесевое
топливо составляет 7 %.
235
№ δз M кр n g NOx g CO g CH
г/(кВт·ч)
г/(кВт·ч)
г/(кВт·ч)
мин-1
%
–
1 0,15 0 nXX – – –
2 0,08 100 nA 8,1 1,11 0,44
3 0,1 50 nB 4,6 1,32 0,32
4 0,1 75 nB 5,6 1,16 0,38
5 0,05 50 nA 4,1 1,34 0,4
6 0,05 75 nA 5,0 1,21 0,43
7 0,05 25 nA 2,9 1,44 0,51
8 0,09 100 nB 8,6 1,05 0,35
9 0,1 25 nB 3,5 1,29 0,38
10 0,08 100 nC 9,1 1,0 0,4
11 0,05 25 nC 4,2 1,18 0,24
12 0,05 75 nC 6,8 1,11 0,11
13 0,05 50 nC 5,5 1,24 0,18
13
=g ∑ ( gi ⋅ δ зi ) , (2)
i =1
gCH = 0,48.
Таким образом, снижение средних эксплуатационных выбросов двигателя 4ГЧ7,5/7,35
по сравнению с бензиновым двигателем-прототипом МеМЗ-307 составило: по NOx - 61
%; по СО - 96 %; по СН - 78 %.
Выводы
Литература
2. Watson H., Milkins E. Hydrogen and methane - automotive fuels of the future / SAE,
Australia, 1975, 3/4, p. 17-19.
5. Michael R.S. The Effect of Hydrogen Addition on Natural Gas Engine Operation /
Michael R.S., Mirza J.Y., Zafer D., Matthew N.S. – Calgary, CA: Univercity of Calgary,
1993. - 16 p. - (Preprint / Univercity of Calgary: SAE № 932775).
6. Karim G.A. Comparative Studies of Methane and Propane as Fuels for Spark Ignition
and Compression Ignition Engines / Karim G.A., Wierzba I. // SAE Transactions. - 1984. - №
92. - pp. 3677-3688.
Майстренко // Механика и машиностроение: сб. науч. тр. - Харьков: ХПИ. – 2010. - №1.
- С. 66-73.
11. Абрамчук Ф.И. Способ подачи водорода для питания автомобильного двигателя
/ Ф.И. Абрамчук, А.Н. Кабанов, А.А. Дзюбенко, А.П. Кузьменко // Автомобильный
транспорт: сб. науч. тр. - Харьков: ХНАДУ. - 2010. - № 26. - С. 63-66.
14. Attar A.A. Optimization and Knock Modelling of a Gas Fueled Spark Ignition Engine:
PhD thesis: Mechanical Engineering / A.A. Attar. Calgary, Canada, 1997. - 248 p.
Cet article a fait l’objet d’une publication dans: Journal of Road Transport Scientific Papers,
Kharkov Ukraine - 2011- N°43 - 278 p- pp. 152-159 - ISSN 2219-8342. Il a été presenté en
langue russe lors de la conférence scientifique internationale: Problems and prospects of the
automotive and road transport, qui a eu lieu à l’Université de l’Automobile et des Ponts et
Chaussées de Kharkiv Ukraine le 3-4 Novembre 2011. Cet article fait partie de la présente
thèse dans le chapitre 7 (Traduction en français) et ici on trouve sa version originale.
УДК 621.43.052
R. Maamri, PhD student, University of Quebec in Trois Riviéres Canada, F.I. Abramchuk,
Ph.D., Associate Professor, KhNAHU, Kharkov, A.N. Kabanov, Ph.D., Associate Professor,
KhNAHU, Kharkov, Y. Dubé, Ph.D., Professor, University of Quebec in Trois Riviéres
Canada, L. Toubal, Ph.D., Professor, University of Quebec in Trois Riviéres Canada, A.
Kodjo, Ph.D., Professor, University of Quebec in Trois Riviéres Canada
Abstract. Results of investigation of gas engine with supercharging converted from diesel
have been given. For this type of engines math model of combustion process with variable
Wiebe combustion factor has been specified.
Key words: toxicity, supercharging, gas engine.
Введение
Цель исследования
Экспериментальный стенд
Степень сжатия на стенде могла изменяться при помощи трёх комплектов поршней,
обеспечивающих следующие значения ε: ε = 9,8; ε = 12 и ε = 14.
243
θ + 18 ηv − 0.25 ε + 1 n + 500
mvar i 10.639 ⋅ ϕi ⋅ ( α + 0.00025 ) ⋅
= × ⋅ ⋅ −
40 0.8 11 5000 , (1)
−2.3 ⋅ ψ − 28.025 ⋅ ϕi2 + 98.045 ⋅ ϕi3 − 156.86 ⋅ ϕi4 + 86.88 ⋅ ϕi5
ϕZ
град. пкв
50
40
30
20
0 0,2 0,4 0,6 0,8 ψ
( hc′ − h )
=
Kd ⋅ (1 − x ) ⋅ ( ε − 1) , (3)
Hu
где hc′ - удельная энтальпия несгоревшей смеси в момент начала сгорания, кДж/кг; h -
удельная энтальпия несгоревшей смеси в любой момент времени, кДж/кг; Hu - низшая
теплота сгорания, кДж/кг, x - доля выгоревшего топлива; ε - степень сжатия.
Результаты исследования
1,6
Q
мДж
1,4
1,0
α = 1,3
0,8
0,6
α = 1,2
0,4
α = 1,1
0,2
α = 1,0
2 4 6 8 H2, %
Рис. 3. Влияние добавки водовода на энергию воспламенения смеси природного газа и
воздуха при различных коэффициентах избытка воздуха
249
Характеристические карты
ψ
%
60
40
20
0 1000
2000
0,2 3000 4000
5000 6000
0,4
0,6
0,8
n, мин-1
1,0
Pe
α
4
1 1000
2000
3000
0,2 4000
5000
0,4 6000
0,6
n, мин-1
0,8
1,0
Pe
θ
град. пкв
до ВМТ
40
30
20
10
0,2
0,4
0,6
0,8
1,0
Pe n, мин-1
В связи с тем, что дроссельная заслонка в двигателе открыта постоянно, под внешней
скоростной характеристикой понимается скоростная характеристика, снятая при
максимальной цикловой подаче топлива в цилиндре.
251
α
1,5
1,4
1,3
1,2
1,1
1,0
0,9
ηv
0,82
0,80
0,78
0,76
Me
Нм
110
100
90
ge
г
кВт год
380
360
340
320
300
280
Ne
кВт
50
40
30
20
10
1000 1500 2000 2500 3000 3500 4000 4500 5000 5500 n, хв-1
_______
- бензин;
_._._.
- природный газ;
_ _ _ _
- смесевое топливо.
Рис. 7. Внешняя скоростная характеристика двигателя 4ГЧ7,5/7,35
252
NOx
млн-1
2500
2000
1500
1000
500
CH
млн-1
180
160
120
CO
%
3,0
2,5
2,0
1,5
1,0
0,5
1000 1500 2000 2500 3000 3500 4000 4500 5000 5500 n, хв-1
_______
- бензин;
_._._._
- природный газ;
_ _ _ _ _
- смесевое топливо.
Рис. 8. Внешняя скоростная характеристика двигателя 4ГЧ7,5/7,35
Как видно из рис. 7, при переводе двигателя на смесевое топливо практически удалось
избежать потерь мощности на номинальном режиме по сравнению с работой при ψ = 0
(Ne = 50 кВт). На режиме максимального крутящего момента Memax (n = 3500 мин-1)
снижение мощности при переходе с природного газа на смесевое топливо составляет 5
%.
5 %, на режиме Memax - на 12 %.
г/(кВт·ч)
г/(кВт·ч)
г/(кВт·ч)
мин-1
%
–
1 0,15 0 XX – – –
2 0,08 100 nA 6,1 0,56 0,16
3 0,1 50 nB 2,5 0,51 0,08
4 0,1 75 nB 4,9 0,48 0,11
5 0,05 50 nA 2,7 0,54 0,12
6 0,05 75 nA 4,2 0,56 0,14
7 0,05 25 nA 0,9 0,32 0,11
8 0,09 100 nB 6,8 0,53 0,13
9 0,1 25 nB 2,1 0,25 0,1
10 0,08 100 nC 7,2 0,5 0,14
11 0,05 25 nC 1,7 0,28 0,05
12 0,05 75 nC 5,1 0,56 0,04
13 0,05 50 nC 3,7 0,47 0,06
13
=g ∑ ( gi ⋅ δ зi ) , (4)
i =1
gCH = 0,09.
Литература
2. Watson H., Milkins E. Hydrogen and methane - automotive fuels of the future? /
SAE, Australia, 1975, 3/4, p. 17-19.
5. Michael R.S. The Effect of Hydrogen Addition on Natural Gas Engine Operation /
Michael R.S., Mirza J.Y., Zafer D., Matthew N.S. - Calgary, CA: Univercity of Calgary,
1993. - 16 p. - (Preprint / Univercity of Calgary: SAE № 932775).
6. Karim G.A. Comparative Studies of Methane and Propane as Fuels for Spark Ignition
and Compression Ignition Engines / Karim G.A., Wierzba I. // SAE Transactions. - 1984. - №
92. - pp. 3677-3688.
11. Абрамчук Ф.И. Способ подачи водорода для питания автомобильного двигателя
/ Ф.И. Абрамчук, А.Н. Кабанов, А.А. Дзюбенко, А.П. Кузьменко // Автомобильный
транспорт: сб. науч. тр. - Харьков: ХНАДУ. - 2010. - № 26. - С. 63-66.
14. Attar A.A. Optimization and Knock Modelling of a Gas Fueled Spark Ignition
Engine: PhD thesis: Mechanical Engineering / A.A. Attar. Calgary, Canada, 1997. - 248 p.
ANNEXE C
Journal of The National Polytechnical University (KPI), Kharkov, Ukraine - 2011- N°54:
New Solutions in Modern Technologies - 156 p - pp. 18-26 - ISSN 2079-5459
1. http://sciencelit.net/author/Maamri_R.html
2. http://goo.gl/QP2YV
3. http://goo.gl/Fo0EK
4. http://goo.gl/6UvuE (9-conference)
5. http://www.pandia.ru/text/77/447/3862.php (conference)
6. http://zaz.gendocs.ru/docs/2800/index-688450.html?page=14 (conference)
7. http://www.tiatd.ru/science/conference/lukans_6_tezis/ (conference)
1. http://dspace.khadi.kharkov.ua/dspace/handle/123456789/206
2. http://cyberleninka.ru/article/n/vybor-parametrov-raboty-malolitrazhnogo-dvigatelya-
rabotayuschego-na-smesi-prirodnogo-gaza-i-vodoroda
3. http://elibrary.ru/contents.asp?issueid=1007180
4. http://elibrary.ru/item.asp?id=17305825
5. http://eco-madi.ru/sites/default/files/Sbornik_lch_2013.pdf