U.S. Environmental Protection Agency Industrial Environmental Research
Office of Research and Development Laboratory
Research Triangle Park, N.C. 27711
MECHANISM AND KINETICS OF
THE FORMATION OF NOX
AND OTHER
COMBUSTION POLLUTANTS:
Phase I. Unmodified Combustion
Interagency
Energy-Environment
Research and Development
Program Report
-------
RESEARCH REPORTING SERIES
Research reports of the Office of Research and Development, U.S.
Environmental Protection Agency, have been grouped into seven series.
These seven broad categories were established to facilitate further
development and application of environmental technology. Elimination
of traditional grouping was consciously planned to foster technology
transfer and a maximum interface in related fields. The seven series
are:
1. Environmental Health Effects Research
2. Environmental Protection Technology
3. Ecological Research
4. Environmental Monitoring
5. Socioeconomic Environmental Studies
6. Scientific and Technical Assessment Reports (STAR)
7. Interagency Energy-Environment Research and Development
This report has been assigned to the INTERAGENCY ENERGY-ENVIRONMENT
RESEARCH AND DEVELOPMENT series. Reports in this series result from
the effort funded under the 17-agency Federal Energy/Environment
Research and Development Program. These studies relate to EPA's
mission to protect the public health and welfare from adverse effects
of pollutants associated with energy systems. The goal of the Program
is to assure the rapid development of domestic energy supplies in an
environmentally—compatible manner by providing the necessary
environmental data and control technology. Investigations include
analyses of the transport of energy-related pollutants and their health
and ecological effects; assessments of, and development of, control
technologies for energy systems; and integrated assessments of a wide
range of energy-related environmental issues.
REVIEW NOTICE
This report has been reviewed by the participating Federal
Agencies, and approved for publication. Approval does not
signify that the contents necessarily reflect the views and
policies of the Government, nor does mention of trade names
or commercial products constitute endorsement or recommen-
dation for use.
This document is available to the public through the National Technical
Information Service, Springfield, Virginia 22161.
-------
EPA- 600/7- 76- 009a
August 1976
MECHANISM AND KINETICS
OF THE FORMATION OF NO
X
AND OTHER COMBUSTION POLLUTANTS
PHASE I. UNMODIFIED COMBUSTION
by
V.S. Engleman and W. Bartok
Exxon Research and Engineering Company
P. O. Box 8
Linden, New Jersey 07036
Contract No. 68-02-0224
ROAPNo. 21BCC-013
Program Element No. 1AB014
EPA Project Officer: W. Steven Lanier
Industrial Environmental Research Laboratory
Office of Energy, Minerals, and Industry
Research Triangle Park, NC 27711
Prepared for
U.S. ENVIRONMENTAL PROTECTION AGENCY
Office of Research and Development
Washington, DC 20460
-------
TABLE OF CONTENTS
Page
FOREWORD i
SUMMARY 1
1. INTRODUCTION 2
2. EXPERIMENTAL PROCEDURES AND RESULTS 3
2.1 Jet-Stirred Combustor 3
2.2 The Multiburner 16
2.2.1 The Multiburner Furnace 16
2.2.2 Burner Design 18
2.2.3 Sampling and Gas Analysis 22
2.2.4 Data Reduction 22
2.2.5 Results for Premixed Flames . 27
2.2.6 Results for Diffusion Flames ,- . . . . 31
3. THEORETICAL CALCULATIONS . 37
3.1 Equilibrium Calculations on Premixed Flames . 37
3.2 Survey of Coupled NO /Combustion Kinetics 37
X
3.3 Kinetics Calculations 44
3.4 Results of Theoretical Calculations for Methane/Air .... 61
3.4.1 80% Stoichiometric Air 61
3.4.2 100% Stoichiometric Air 62
3.4.3 120% Stoichiometric Air 66
3.5 Discussion 66
3.6 Reactions of Importance at the Interface Between the
Branching Zone and the Relaxation Zone 69
3.7 Comparison of Theoretical Calculations with
Experimental Results 72
4. CONCLUSIONS 74
APPENDIX A - MASTER LIST OF REACTIONS FOR 25 ALLOWED SPECIES .... A-l
APPENDIX B - CROSS INDEX OF REACTIONS FOR 25 SPECIES B-l
APPENDIX C - KINETIC DATA ANALYSIS FOR JET-STIRRED COMBUSTOR .... C-l
APPENDIX D - NO/NO PLOTS D-l
X
APPENDIX E - CO/C02/HC/02 DATA PLOTS E-l
APPENDIX F - DATA LISTINGS F-l
-------
- 1 -
FOREWORD
This report summarizes the results of Phase I of a study on
the "Definition of the Mechanism and Kinetics of the Formation of NOX
and Other Pollutants under Normal and Combustion Modification Conditions".
This study was conducted by Exxon Research and Engineering Company under
Contract 68-02-0224 funded by the Environmental Protection Agency.
The helpful comments and suggestions of Messrs. W. S. Lanier,
G. B. Martin, and D. W. Pershing, project officers for this contract,
are gratefully acknowledged. The skillful assistance of Messrs. F. D. Remyn
and R. M. Buono in conducting the laboratory portion of this program is
also acknowledged.
-------
- 1 -
SUMMARY
The objective of this study was to investigate the mechanism
and kinetics of the formation of NOX and other pollutants in combustion.
A combined experimental and theoretical study was undertaken for this
purpose. This report covers Phase I, the study of combustion under
unmodified conditions and calculations using the best available kinetic
data from the literature.
A jet-stirred combustor was used to extend the range and accuracy
of data taken in previous studies. Since precise measurements of
concentrations and temperatures are critical to allow meaningful results
in combustion calculations for a well-stirred reactor, additional studies
on the combustion of hydrogen, carbon monoxide, methane, and propane
with air were made. These studies of coupled combustion/pollutant
formation indicated substantial heat loss from a conventional stirred
reactor and the need for the development of an adiabatic stirred combustor.
A furnace capable of studying combustion under adiabatic
conditions (called the multiburner because it was designed for multifuel
firing capability in the same combustion zone) was further refined for
these studies. This unit is an electrically-heated furnace with a
zirconia muffle tube that is capable of attaining temperatures up to
about 2500°K. Premixed flames of the flat-flame and focused-flame type
as well as laminar and turbulent diffusion flames were studies using
methane and propane as fuels. These studies included both adiabatic
and heat-loss conditions. In addition, a few runs were made with wall
temperatures above the adiabatic flame temperatures.
Both stirred reactor and plug flow calculations were made.
Stirred reactor calculations indicated the need for more detailed
kinetics in hydrocarbon-air combustion for the prediction of NOX formation.
Plug flow calculations (with kinetic data available from the literature
at the time) indicated strong coupling between combustion reactions and
NOX formation in the flame zone. However, since no reaction rate data
were available to test the direct coupling between hydrocarbon fragments
and nitrogenous species, further elucidation of the kinetics is required
to determine the importance of such reactions for NOX formation.
-------
- 2 -
1. INTRODUCTION
Under Contract No. 68-02-0224, sponsored by the Environmental
Protection Agency, Exxon Research and Engineering Company conducted a
study to relate the kinetics and mechanism of pollutant formation reactions
and those of hydrocarbon reactions. The emphasis in this program was
placed on NOX formation and destruction reactions as they relate to
hydrocarbon combustion. This program was performed in two phases. This
report covers the work in Phase I which included studies of unmodified
combustion in a non-adiabatic stirred combustor and a plug flow combustor
capable of adiabatic operation. Preliminary kinetics calculations based
on the most complete kinetic data available from the literature at the
time were also conducted in Phase I. Phase II, which is covered in a
companion report, included experimental studies of modified combustion in
a plug flow system as well as studies in a newly developed adiabatic
stirred combustor. Kinetics calculations, using updated literature
information as well as estimates for potentially important reactions not
available from the literature, were also performed.
The purpose of these studies was to provide further understanding
of the coupling between combustion reactions and pollutant formation (with
emphasis on NOX) leading to the capability to perform accurate predictive
calculations on complex combustion systems. The ultimate aim of this and
other fundamental combustion research studies being conducted for the
Combustion Research Branch of the Industrial Environmental Research
Laboratory of EPA is developed control technology for stationary sources.
Clearly, in a practical system, there will be interactions
among chemistry, fluid mechanics and heat transfer. However, without
an understanding of the basic chemistry of combustion, agreement between
theory and experiment in a particular case may be fortuitous and will not
allow application of the same theory to a new system. Thus, the theory
becomes merely an empirical correlation. These studies on the chemistry
of combustion/pollutant formation were therefore undertaken, in order to
progress towards a theory of general applicability to a wide variety of
combustion systems. These studies were oriented toward studying combustion
under idealized conditions, either controlled by chemical kinetics or
with minimum, or well-defined fluid mechanic and heat transfer interactions.
-------
- 3 -
2. EXPERIMENTAL PROCEDURES AND RESULTS
Two experimental devices were used In the experimental portion
of this study: (1) the jet-stirred combustor—non-adiabatic stirred
reactor capable of operating under kinetically-limited combustion conditions
and (2) the multiburner—plug flow reactor, with interchangeable burners,
capable of operating under adiabatic conditions. This section covers the
experimental portion of the study, while the theoretical considerations
will be discussed in the following section.
2.1 Jet-Stirred Combustor
The jet-stirred combustor used in this program is a modification
of the Longwell-Weiss reactor (2-1) with hemispherical geometry. This
device was selected because it has been used extensively in fluid mechanic
and combustion modeling and because the combustion rates are limited by
chemical kinetics as opposed to transport effects (2-2) . The apparatus
used was described in the Final Report for the preceding program (conducted
under Contract No. CPA 70-90, Reference 2-3) whose purpose was to investigate
the basic factors affecting nitric oxide formation in the combustion of
fossil fuels. The reactor (Figure 2-1) consists of an outer shell of
fire-brick shaped as two halves of a sphere three inches in diameter.
The upper hemisphere is solid with the exception of the hole through which
the reactants are brought to the injector. The lower hemisphere is
hollowed out to a reaction zone of 1.5 inch diameter. The insulating
shell has twenty-five holes of 0.125 inch diameter, through which the
burned mixture exits.
Fuel and air are metered separately through calibrated rotameters,
preheated to the desired inlet temperature and then mixed before entering
the combustor. The temperature of the fuel/air stream is measured
immediately before injection. The fuel-air mixture enters the reaction
zone through a stainless steel injector which is a hemisphere into which
are drilled forty radial holes of 0.020 inch diameter. The reactants
enter the reaction zone as small sonic jets which stir the reactor
contents and produce a mixture of essentially uniform temperature and
composition in a characteristic time which is short compared with the
average residence time. The combustion experiments were conducted at
atmospheric pressure with a range of residence times from 1-1/2 to 4
milliseconds.
2-1 Longwell, J. P. and Weiss, M. A., Ind. Eng. Chem. 47, 1634 (1955).
2-2 Hottel, H. C., Williams, G. C., and Miles, G- A., Eleventh Symposium
(International) on Combustion, p. 771, The Combustion Institute, 1967.
2-3 Bartok, W., Engleman, V. S., and del Valle, E. G., Laboratory Studies
and Mathematical Modeling of NOX Formation in Combustion Processes,
Final Report EPA Contract CPA 70-90, December 1971.
-------
- 4 -
FIGURE 2-1
SCHEMATIC OF THE JET-STIRRED COMBUSTOR
INCONEL SPHERE
FIRE BRICK
INJECTOR
THERMOCOUPLE
PRE-MIXED AIR AND
FUEL INLET
FIRE-BRICK
WATER-COOLED PROBE
-------
- 5 -
Gas temperatures were measured up to temperatures slightly in
excess of 2000°K with a Pt/Pt-10% Rh thermocouple 0.010 inches diameter.
Radiation corrections were found to be negligible because of the intense
stirring in the combustion zone (see Appendix C). Duplicate runs were
made with silica coated and uncoated thermocouples in selected cases to
check for catalytic effects on temperature measurements. No measurable
differences were found. The thermocouple was movable so that traverses
could be taken during a run. It was found that the temperatures were
quite uniform throughout the reaction zone although slightly cooler in
the immediate vicinity of the injector sphere and in the immediate
vicinity of the outer wall when the thermocouple was drawn back into its
port.
The combustion gases were sampled through a water-cooled
stainless steel probe, 0.125 inches outside diameter and 0.033 inches
inside diameter similar to that used by Longwell and Weiss (2-1). The
probe was placed through a port in the shell of the combustor. The
quenched gases were pumped through a diaphragm pump to the combustion
gas analyzers. The most important species analysis in this study was
the measurement of NOX concentrations. During the preceding program
(CPA 70-90, Reference 2-3) these analyses were performed primarily with
the Envirometrics Multigas Analyzer although cross-checks were made with
the DuPont Photometric NOX analyzer. During the present program
(68-02-0224) a Thermo Electron Chemiluminescence Analyzer was used and
the results of both programs served as cross-checks of the data. Care
had to be taken when using the Envirometrics analyzer in the presence
of CO, although since the response time to CO is slower than the response
time to NOX, CO concentrations of up to 5-10% could be handled with
minimum interference by cyclic sampling/purge techniques which took
advantage of the different response times for NOX and for CO.
The hydrogen was Linde extra dry grade with a stated minimum
purity of 99.95%, the carbon monoxide was Matheson C.P. grade with a
stated minimum purity of 99.5% and the propane was Matheson C.P. grade
with a stated minimum purity of 99.0%. The air was Baker dry grade
manufactured by mixing nitrogen and oxygen with a typical argon concentration
of 450 ppm. The oxygen concentration of each cylinder was measured and
was generally found to be 21.0%.
Precise experimental measurements of concentrations and temperature
are extremely critical for achieving meaningful results in combustion
calculations for a well-stirred reactor. To assure the precise nature
of the measurements made in the preceding program (reference 2-3) and
to extend those measurements to include those parameters necessary for
kinetic calculations for the assessment of the importance of combustion
intermediates in coupled combustion/pollutant formation, further studies
were made with the jet-stirred combustor for hydrogen, carbon monoxide,
methane, and propane.
Temperature measurements for methane/air combustion are shown
in Figures 2-2 and 2-3. Figure 2-2 shows temperature vs. mixture ratio
with the solid line being the same one indicated in reference 2-3. The
dotted line is the extension of this line to fuel-rich conditions with
the peak indicated at about 90% stoichiometric air. In reference 2-3
temperature measurements were taken only under fuel lean conditions at
these flow rates. An indication of the temperature uniformity within the
-------
- 6 -
FIGURE 2-2
TEMPERATURE AS A FUNCTION OF MIXTURE RATIO
JET-STIRRED COMBUSTOR / METHANE - AIR
H
W
W
2000
1900
1800
1700
1600
50
O
FROM REFERENCE 2-3
J L
100
PERCENT STOICHIOMETRIC AIR
150
FIGURE 2-3
TEMPERATURE PROFILES AT REPRESENTATIVE STOICHIOMETRIES
JET-STIRRED COMBUSTOR / METHANE - AIR
EH
<
tf
2000
1900
1800
1700
1600
.0--° O O
0
a
•a
.A.. A.
0.25
0.50
M
Vo,
0.75
% ST AIR
73
123
n
'[]
65
141
163
1.00
1.25
-------
- 7 -
jet-stirred combustor with methane is given in Figure 2-3. Temperature
is plotted against (r/Ro)3} a normalized volumetric parameter, indicating
the portion of the reactor volume at a given temperature. It can be
seen that the temperatures are quite uniform throughout the combustor.
Near the injector sphere the temperature apparently drops off but this
is probably caused by a combination of convective cooling by the cool
inlet jets in the vicinity of the injector sphere and radiative cooling
of the thermocouple bead by the injector sphere. The temperatures used
for kinetics calculations are indicated by the horizontal dashed lines.
The NOX emissions from the jet-stirred combustor with methane/air
are indicated in Figure 2-4. It should be noted that peak NQx formation
occurs under slightly fuel-rich conditions. These data were measured
with a Thermo-electron NO/NCx chemiluminescence analyzer and are in very
good agreement with the results obtained in reference 2-3 using Envirometrics
electrochemical analyzers and a DuPont photometric analyzer.
The oxygen results were obtained with a polarographic oxygen
analyzer. They are shown in Figure 2-5 and compare favorably with those
reported in reference 2-3 which were taken with a paramagnetic oxygen
analyzer. The results are in good agreement with calculated values of
residual oxygen in a well-stirred combustor.
The effect of the water-cooled sampling probe on the temperature
of the jet-stirred combustor is dramatic if species probe traverses are
attempted. As long as the probe is outside the volume of the reactor
(up to a position flush with the wall) the temperature in the reactor
remains fairly constant. As the probe enters the reaction zone, the
temperature begins to drop, and it drops more than 100 K by the time
the probe is three fourths of the way to the center of the reactor. This
is illustrated in Figure 2-6 for which the uncooled thermocouple probe
was held in a fixed position while the cooled species sampling probe was
inserted into the reactor.
Temperature measurements for carbon monoxide/air combustion are
shown in Figures 2-7 and 2-8. Figure 2-7 shows temperature vs. mixture
ratio with the solid line being the one calculated on a theoretical
basis as indicated in the paper presented at the Fourteenth Combustion
Symposium and shown in Appendix C. The data have been extended on the
fuel rich side beyond the data available at the time of writing the
Combustion Symposium papert The temperature uniformity in the combustion
zone for carbon monoxide air is also good as indicated in Figure 2-8.
The temperatures used for kinetics calculations are indicated by the
horizontal dashed lines.
The NOX emissions for CO/air are shown in Figure 2-9 and were taken
with a chemiluminescence analyzer. The agreement with the results reported
in reference 2-3 is quite good on the lean side. There is some disagreement
under fuel rich conditions between 50 and 65% stoichiometric air but
the present results should be considered more reliable because the CO
content in the exhaust exceeds 10%. That high a concentration is difficult
to eliminate as an interference with the NOX Paristor (even with the
special cyclic sampling techniques developed (2-3)) when the NOX level
is below 50 ppm; the chemiluminescent analyzer does not have such a CO
interference.
-------
FIGURE 2-4
w
<
W
ra
§
x
O
OS
W
O
O
75
50
25
10
50
NO AS A FUNCTION OF MIXTURE RATIO
JET-lTIRRED COMBUSTOR / METHANE - AIR
O
O
O
6
50 100 150
PERCENT STOICHIOMETRIC AIR
FIGURE 2-5
OXYGEN AS A FUNCTION OF MIXTURE RATIO
JET-STIRRED COMBUSTOR / METHANE - AIR
u
O
o .o Q"°. I
1 - 1 _ 1 _ 1
100 150
PERCENT STOICHIOMETRIC AIR
-------
- 9 -
FIGURE 2-6
EFFECT OF WATER-COOLED SAMPLING PROBE POSITION
ON JET-STIRRED COMBUSTOR TEMPERATURE
o
EH
S
O
o
w
EH
W
w
H
1800
1700
1600
EXHAUST PORTS
COMBUSTOR—*-
_-o o
THERMOCOUPLE MAINTAINED
IN FIXED POSITION AT
r/fe0 = l.,00 ,
'NORMAL
[SAMPLING
POSITION
INJECTOR »
250
200
150
100
50
0
COOLED PROBE POSITION
(ARBITRARY UNITS)
-------
- 10 -
FIGURE 2-7
TEMPERATURE AS A FUNCTION OF MIXTURE RATIO
JET-STIRRED COMBUSTOR / CARBON MONOXIDE - AIR
w
EH
2100
200(_
1900
180C-
1700-
CALCULATED TEMPERA
TURE (SEE APPENDIX C)
1 1 1 1
1 1 1 1
I 1
50 100 150
PERCENT STOICHIOMETRIC AIR
FIGURE 2-8
TEMPERATURE PROFILES IN JET-STIRRED AT
REPRESENTATIVE STOICHIOMETRIES COMBUSTOR / CARBON MONOXIDE - AIR
O
O
m
tf
W
2100
2000
1900
1800
1700
si-
j
•3
A
D i
n ~
/•o--0 ° o c
o'
0'
1 1 1
PERCENT
STOICHIO-
METRIC
AIR
» 117
J 53
1^.
)^ 142
0 0.25 0.50 0.75 1.00 1.25
, ^ 3
-------
AS MEASURED
PH
&
NO
X
OXYGEN (PERCENT
co
01
o
o
-^
01
^
O
M
01
- 11 -
FIGURE 2-9
NOX AS A FUNCTION OF MIXTURE RATIO
JET-STIRRED COMBUSTOR / CARBON MONOXIDE - AIR
f
o
1
50
100 150
PERCENT STOICHIOMETRIC AIR
FIGURE 2-10
OXYGEN AS A FUNCTION OF MIXTURE RATIO
JET-STIRRED COMBUSTOR / CARBON MONOXIDE - AIR
10
- O EXPERIMENTAL
THEORETICAL (SEE APPENDIX
50
100 150
PERCENT STOICHIOMETRIC AIR
-------
- 12 -
The oxygen concentration is shown In Figure 2-10 and is in
good agreement with theoretical calculations indicated by the solid line.
Careful temperature probing of the reaction zone was accomplished
and a study of the effect of preheat on NOX formation was undertaken for
propane. Data were also taken with hydrogen fuel for better characterization
of the flame temperature vs. mixture ratio for this system in the jet-stirred
combustor, and to obtain more accurate NOX measurements under extremely lean
and extremely rich conditions.
The NOX measurements for the propane-air system are presented
in Figure 2-11. The agreement with data taken previously is good in the
case of no preheat. The NOX peak appears to occur on the slightly
fuel-rich size (at about 90% stoichiometric air).
The results for preheat in Figure 2-11 indicate only a slight
effect on NOX formation. However, temperature measurements show that
only at 150% stoichiometric air conditions did the combustion temperature
reflect fully the nominal amount of preheat. At 135% stoichiometric
air the difference in combustion temperatures was only half of the
preheat and there was no difference in combustion temperatures under fuel
rich conditions. This lack of recovery of prehat temperature in the
combustion zone is directionally consistent with adiabatic equilibrium
calculations in which higher preheat causes dissociation at high temperatures
and some of the sensible heat is converted to latent heat. At 150%
stoichiometric air the difference of about 100°C in combustion temperature
results in a NOX increase by more than a factor of two. This behavior
is in line with theoretical predictions of well-stirred reactor performance.
Temperature and oxygen measurements for the propane/air system
are shown in Figure 2-12. The difference between theoretical adiabatic
well-stirred reactor temperatures and experimentally increased temperatures
indicates a heat loss on the order of 30-40 cal./sec. over the conditions
studied. The measured oxygen concentrations match the calculated oxygen
concentrations for a well-stirred reactor operating at the measured
temperatures.
Probing of combustor temperature profiles showed that temperatures
were quite uniform in the jet-stirred combustor for propane-air. The
temperatures are plotted in Figure 2-13 against (r/Ro)3. As discussed
before, this parameter represents the fraction of total volume in a
given zone of the combustor. The temperatures were found to be quite
uniform throughout the combustor from the vicinity of the injector to
within a short distance of the wall. Temperature probing was performed
both with bare metal thermocouples and with thermocouples coated with
silica. Because of the intense stirring and the apparent lack of catalytic
effects on the thermocouples under the conditions studied, no difference
was observed between the measurements obtained with coated and uncoated
thermocouples.
Measurements were also made in the hydrogen-air system. The
measurements were made under fuel-rich conditions around 40% stoichiometric
air and under fuel lean conditions between 140% and 250% stoichiometric
air. The chemiluminescence analyzer permits more accurate analysis of
combustion gases at low NOX levels than we were able to achieve during
the previous investigation (2-3) which employed Envirometrics electro-
chemical analyzers.
-------
- 13 -
FIGURE 2-11
N<> FORMATION IN JET-STIRRED COMBUSTOR (PROPANE/AIR)
'••
X
Q
W
B
03
<
W
X
150
100
50''li.
O1,
-\Q
I- NEAR
BLOWOUT
k/
- NO PREHEAT
D - "WITH PRE-
HEAT
(SEE TEXT).
50
100
150
PERCENT STOICHIOMETRIC AIR
-------
- 14 -
FIGURE 2-12
TEMPERATURE AND OXYGEN LEVELS
PROPANE - AIR
w
E-i
W
W
H
2200 L
2100 L-
2000 I—
1900 (-
ADIABATI
WELL-ST
REACTOR
EXPERIMENTAL
100 150
PERCENT STOICHIOMETRIC AIR
-------
15 -
FIGURE 2-13
TEMPERATURE PROFILES IN JET-STIRRED COMBUSTOR
PROPANE - AIR
U
o
tt
o
H
— «-n — -D-— --Q ~ -n
j
B ^
-o
^. ••._-•_„-.•.. —
* 1 1 1
% ST AIR
nc
95
^ 120
1
60
150
0
0.25
0.50
0.75
1.00
1.25
-------
- 16 -
Calculations were performed to show that the jet-stirred
combustor operates as a well-stirred reactor for "clean" combustion
systems (hydrogen and carbon monoxide) and calculations have been performed
on propane to test a quasi-global combustion model for calculations of
NOX formation. The results (2-4) of this study are reported in Appendix C.
The calculations have been valuable to our current program to point out
critical areas for experimentation and further kinetic studies. As
discussed in detail in Appendix C, the theoretical analysis based on a
simplified quasi-global approach failed to produce predicted NOX
concentrations in agreement with the experimental results for the
propane/air system. This lack of agreement was particularly pronounced
under fuel-rich conditions. Therefore, it was concluded that the
investigation should be broadened to a detailed treatment of the kinetic
data and also, that further experimental results should be obtained under
conditions which eliminated heat loss as well as mixing limitations for
comparison in theoretical predictions.
2.2 The Multiburner
Studies in the multiburner were undertaken to allow investigation
of combustion reactions with well-controlled heat loss. The multiburner
was operated as a flow reactor for this portion of this study with both
hot and cold walls. The multiburner studies conducted indicated the
need for the development of a new combustor capable of behavior as a
well-stirred reactor under adiabatic conditions. These complementary
studies with stirred reactors and flow reactors provide data on a range
of conditions from well-mixed combustion controlled by chemical kinetics
to diffusion flames controlled by fluid mechanics.
2.2.1 The Multiburner Furnace
The multiburner (so-called because it was designed to burn gas,
oil or coal in the same idealized combustion zone) is an electrically-
heated furnace with a zirconia muffle tube that is capable of attaining
temperatures up to about 2500°K. The furnace is shown schematically in
Figure 2-14. The furnace heating zone has tungsten mesh heating elements
and molybdenum radiation shields and is purged with gettered argon to
extend furnace component lifetimes. High temperature gasketing for the
furnace is composed of zirconia felt. Temperatures inside the furnace
are monitored with tungsten-rhenium thermocouples which also serve as
the input to the temperature controller for the furnace when operating
in the automatic mode. Wall temperatures can also be monitored by optical
pyrometry through three sight ports provided in the side of the furnace.
The power supply is a 25 KVA AC unit which is also equipped with a thermal
watt converter for input to the controller in cases where it is desirable
to control power input rather than strictly the temperature of the furnace.
The controller on the power supply is also capable of operating in a
"manual" mode, in which case current input to the furnace is controlled.
2-4 Engleman, V. S., Bartok, W., Longwell, J. P., and Edelman, R. B.,
Fourteenth Symposium (International) on Combustion, p. 755, The
Combustion Institute, 1973.
-------
FIGURE" 2-14
MULTIBURNER FURNACE
SEAL
SEAL
MOLYBDENUM
RADIATION
SHIELDS
WATER-
COOLED
JACKET
TUNGSTEN-
HEATING
ELEMENT
STABILIZED
ZIRCONIA
MUFFLE TUBE
-------
- 18 -
The muffle tube in the furnace has a heated length of 15
inches and the "constant temperature section" of the heated zone is
uniform within + 10°C under combustion operation conditions. The muffle
tube is quite stable at high temperatures and can be exposed to varying
thermal environments in the temperature range above 1300°C. However,
extreme care must be exercised in heating up and cooling down the furnace
since the zirconia is most sensitive to thermal shock in the temperature
range around 900°C.
The multiburner can be operated with the combustion zone wall
temperature at any desired level and at the extremes it can either be
used to add heat to the combustion zone or, if desired, to remove heat.
Thus, the multiburner can operate- with all of the flexibility of a
laboratory combustor without the high heat loss generally found in
small scale devices. This permits the use of the multiburner under
combustion conditions that would not be possible with the typical
small scale combustor.
2.2.2 Burner Design
The burners used in this part of the study were designed to
be simple examples of the following types of burners:
• premixed flat flame
• premixed focused flame
• diffusion flame - with interchangeable fuel nozzles
The simple design of each of these burners (shown in Figures 2-15, 16
and 17) was aimed at ease of comparison with analytical calculations.
Each of the burners was readily interchangeable with the others because
they were all built into the same size water jacket, which served to
maintain constant burner temperatures during combustion experiments.
Even though the burners produced three different basic flame types,
the total burner face area, through which combustion gases pass, was
maintained constant for all three types.
The flat flame burner has a porous stainless steel disc across
its face which provides a flat velocity profile for the combustion gases
entering the furnace. The gases are mixed immediately before entering
the burner body and pass through a pair of mixing/calming plates before
exiting through the porous disc. A diagram of the flat flame burner
is given in Figure 2-15.
The premixed focused flame burner used in this study is char-
acteristic of a commercial gas furnace burner; the burner has a single
center hole and eight pilot holes to stabilize the flame. As in the
flat flame burner, the fuel and air are mixed immediately before entering
the burner body and pass through mixing plates before exiting from the
burner face. A diagram of the premixed focused flame is given in
Figure 2-16.
-------
_ 1 Q _
FIGURE 2-15
FLAT FLAME BURNER
POROUS S.S. DISC
WATER-COOLED
TUBE
FLOW STRAIGHTENERS
WATER-COOLING JACKET
GAS FLOW CHA
WATER INLET
PREMEXED
GAS INLET
WATER OUTLET
WATER INLET
-------
- 20 -
FIGURE 2-16
PREMKED FOCUSED FLAME BURNER
FLAME ARRESTOR
WATER COOLING JACKET
GAS FLOW CHANNEL
WATER INLET
FLAME ARRESTOR
PREMIXED GAS
UNIT
-------
- 21 -
FIGURE 2-17
DIFFUSION FLAME BURNER
FLOW STRAIGHTENERS
WATER COOLING JACKET
WATER INLET
AIR INLET
\
FUEL INLET
-------
— 22 —
The diffusion flame burner is of the concentric type with the
fuel flowing through the inner tube and the air flowing through the outer
one. Two types of fuel injector were tested; one with a single fuel
jet and the other with an orifice plate containing six 0.020 inch
diameter holes covering the fuel jet to produce turbulence and increase
mixing. The latter, called the "stabilized" diffusion burner was found
to give better performance. The single fuel jet burner was tested in
two different configurations; one in which the fuel jet diameter was
kept constant when changing fuels from methane to propane allowing the
velocities to change between the fuels, and one in which the fuel velocity
was kept constant for the two fuels by changing the fuel jet diameter.
A diagram of the diffusion flame burner is given in Figure 2-17.
2.2.3 Sampling and Gas Analysis
Combustion gas samples were extracted through a cooled quartz
probe drawn to a nozzle tip. The sampling probe was cooled to 65 C to
prevent water condensation in the probe. The Teflon sampling lines were
also heated to prevent water condensation before the water knockout. The
sample touched only quartz and Teflon lines until it was cooled sufficiently
to prevent further reaction at which point stainless steel lines were used.
The combustion gases were analyzed for NO and NOx with a
Thermoelectron Model 10A Chemiluminescence analyzer. The analysis for
NO is accomplished directly while the analysis for NOx is accomplished
by converting N0£ (and under certain conditions other nitrogenous
compounds) to NO and then analyzing for NO.
Analyses for CO and C02 were performed with two MSA 303 Analyzers,
both of the NDIR type. Stacked cells each with a dual range switch allowed
full scale ranges from 0.1% to 15% for CO and from 0.5% to 25% for C02-
Oxygen analysis was performed with a Beckman Model 742
polarographic oxygen analyzer with full scale ranges from 1% to 25%.
Hydrocarbon analyses were performed with a Beckman Model 400
flame ionization analyzer with full scale ranges calibrated from 50 ppm
to 5% hydrocarbon as methane.
2.2.4 Data Reduction
Experimental data as taken were coded onto computer input forms
and were reduced to measured species concentrations by a data reduction
program. The data produced in this phase of the program are presented
in the Appendix. A further subroutine can also be called to plot the
results on a printer plot. The data plots are also found in the Appendix.
A summary of runs completed in this phase is presented as
Table 2-1. The table is divided into three sections; section A covers
the cold wall tests, section B the hot wall tests on premixed flames,
and section C the hot wall tests on diffusion flames. For each of the
runs in the table fuel type, burner type, firing rate, mixture ratio,
and type of probing are given. Wall temperature is shown for the hot
wall tests. The experimental findings on premixed flames are discussed
in Section 2.2.5 and the findings on diffusion flames are discussed in
Section 2.2.6.
-------
TABLE 2-1
A. Cold
Run No.
101
102
103
104
105
106
107
108
109
110
111
112
113
114
116
119
141
* I/O =
AX =
RD =
Wall Tests
Fuel Type
Propane
Propane
Propane
Propane
Methane
Methane
Methane
Methane
Methane
Methane
Methane
Methane
Methane
Methane
Methane
Propane
Methane
Input/Output
Axial
Radial
SUMMARY OF RUNS
Burner Type
Flat Flame
Flat Flame
Flat Flame
Flat Flame
S. Diff**
S. Diff
S. Diff
S. Diff
S. Diff
S. Diff
Flat Flame
Focused Flame
Focused Flame
Focused Flame
Flat Flame
Flat Flame***
Flat Flame
MADE WITH MULTIBURNER (PHASE I)
Firing Rate
SCFM Air
0.854
0.546
0.546
0.546
2.172
2.172
2.172
2.172
0.264
0.264
0.546
2.067
2.067
2.067
0.481
2.094
0.481
Mixture Ratio
% Stoich Air
83
98
120
101
106
103
122
94
83
101
101
120
120
107
128
145
86
104
96
115
76
119
80
110
132
88
101
Probing*
I/O
I/O
I/O
AX
AX/RD
AX/RD
AX/RD
AX/RD
AX/RD
AX/RD
AX/RD
RD
AX/RD
AX/RD
AX/RD
AX/RD
AX/RD
AX/RD
AX/RD
AX/RD
AX/RD
AX/RD
AX/RD
AX/RD
AX/RD
AX/RD
AX/RD
** Stabilized Diffusion
*** High
Velocity Flame
(Lifted)
1-0
OJ
-------
TABLE 2-1 (CONTINUED)
SUMMARY OF RUNS MADE WITH MULTIBURNER (PHASE I)
B. Hot Wall Tests - Premixed Flames
Run No. Fuel Type Burner Type Firing Rate Mixture Ratio
SCFM Air % Stolen Air
115 Methane Flat Flame 0.481 161
141
60
117 Methane Flat Flame 0.481 161
119
80
118 Propane Flat Flame 0.481 161
61
120 Propane Focused Flame 0.481 156
121 Propane Focused Flame 0.481 125
101
101
82
117
78
95
117
95
78
101
130 Methane Flat Flame 0.481 140
131 Methane Flat Flame 0.481 120
78
111
101
Probing
RD
RD
RD
AX/RD
AX/RD
AX/RD
AX/RD
AX/RD
AX/RD
AX/RD
I/O
I/O
I/O
I/O
I/O
I/O
I/O
I/O
I/O
I/O
AX/RD
AX/RD
AX/RD
AX/RD
AX/RD
Wall
Temp.
°C
1703
1853
1783
1703
1703
2098
1723
1838
1768
2068
2068
2138
2138
2068
2068
2068
2138
2138
2138
2138
1858
2043
2043
2138
2218
I
N3
-C-
-------
TABLE 2-1 (CONTINUED)
SUMMARY OF RUNS MADE WITH MULTIBURNER (PHASE I)
B. Hot Wall
Run No.
131(cont)
132
133
134
135
136
137
138
139
140
Tests - Premixed Flames (cont'd)
Fuel Type
Me thane
Methane
Methane
Me thane
Propane
Propane
Propane
Propane
Methane
Me thane
Burner Type
Flat Flame
Focused Flame
Focused Flame
Focused Flame
Focused Flame
Focused Flame
Focused Flame
Focused Flame
Focused Flame
Focused Flame
Firing Rate
SCFM Air
0.481
0.481
2.094
0.481
0.481
2.094
2.094
0.481
0.481
2.094
Mixture Ratio
% Stolen Air
111
120
131
140
151
91
80
160
160
140
140
65
139
141
121
121
120
120
Probing
I/O
I/O
I/O
I/O
I/O
I/O
I/O
AX/RD
AX/RD
AX/RD
AX/RD
AX/RD
AX/RD
AX/RD
AX/RD
AX/RD
AX/RD
AX/RD
Wall
Temp.
°C
2218
2218
2218
2218
2218
2218
2218
1713
1713
1863
1858
1873
1898
1883
2068
2068
2043
2048
-------
TABLE 2-1 (CONTINUED)
SUMMARY OF RUNS MADE WITH MULTIBURNER (PHASE I)
C. Hot Wall Tests - Diffusion Flames
Run No.
122
123
124
125
126
127
128
129
Fuel Type
Methane
Methane
Me thane
Propane
Propane
Propane
Methane
Methane
Burner Type
S.
s.
S.
s.
s.
s.
Large
S.
Diff
Diff
Diff
Diff
Diff
Diff
Single Hole
Diff
Firing Rate
SCFM Air
2
2
2
2
2
2
2
2
.094
.094
.094
.094
.094
.094
.094
.094
Mixture Ratio Probing
% Stoich Air
160
140
120
80
160
140
120
100
80
100
160
141
114
80
160
140
140
AX/RD
AX/RD
AX/RD
AX/RD
I/O
I/O
I/O
I/O
I/O
AX
AX/RD
AX/RD
AX/RD
AX
AX/RD
AX/RD
RD
Wall
Temp.
°C
1723
1863
2048
2098
2228
2228
2228
2228
2228
2228
1698
1883
2143
2143
1713
1863
1863
I
N>
-------
- 27 -
2-2.5 Results for Premixed Flames
Species profiles have been obtained for premixed flames of
methane-air and propane-air in the multiburner with hot and cold walls.
Two types of premixed flames were studied: (1) the flat flame in which
the flame zone is flat across the burner face and the flame is essentially
one-dimensional, and (2) the shaped flame which is more representative of
premixed furnace burners although much more complex than the flat flame
in a fluid mechanical sense. The cold-wall runs are listed in Table 2-1(B).
The data listings for these runs are presented in Appendix F. The NOX
plots for the individual runs are presented in Appendix D and the plots
for CO/C02/HC/02 are presented in Appendix E. A summary of the types of
data obtained on premixed flames will be presented in this section while
the comparisons with theoretical calculations will be discussed in
Section 3.
Centerline species profiles (CO, C02, 02, NO and NOX) for a -
premixed flat flame of methane/air at 104% stoichiometric air are presented
in Figure 2-18. While the hydrocarbons are not plotted in this figure,
about 20,000 ppm were measured within 0.05 inches of the burner face,
but the concentration had dropped off to 70 ppm at 0.10 inches, and to
essentially zero at 0.50 inches. The oxygen profile drops rapidly from
about 7% within 0.05 inches to about 2.4% at 0.10 inches, and to 1.5% at
1 inch, settling at about 0.8% at greater distances. It was not possible
to probe close enough to the burner face to permit observation of the
increasing portion of the CO profile, but the decreasing portion is
plotted from a distance of 0.1 inches outward with only minor concentrations
of CO measurable after 5 inches. The C02 concentration increases rapidly,
reaching 90% of its ultimate value within 0.5 inches of the burner face.
The NOx plot exhibits substantial curvature in the region up to about
two inches from the burner face and increases only gradually thereafter.
It is interesting to note that the NO measurements are only
about 2/3 of the NOX measurements up to one inch and about 3/4 of the NOX
measurements up to two inches. At 10 inches the NO is about 90% of the
NOX= There are a number of possible interpretations for this early N02
peak.* One possible interpretation is that during the quenching process
in the probe, a portion of the oxygen atoms combine with the NO to form
N02- Calculations indicate that, with instantaneous quenching to 300°K
and assuming no heterogeneous reactions, 100 ppm of NO and 1000 ppm of
oxygen atoms react to form 38 ppm of N02 and 62 ppm of NO (allowing all
possible reactions between species containing nitrogen and oxygen but
ignoring species containing carbon or hydrogen). Such a calculation,
while oversimplified by the idealized quenching and limitation of species,
should still be somewhat conservative because of the probable underestimation
of oxygen atom concentration in the flame zone.
Cernansky, N. P., and Sawyer, R. F., "NO and N02 Formation in a Turbulent
Hydrocarbon/Air Diffusion Flame", Fifteenth Symposium (International)
on Combustion, The Combustion Institute, Pittsburgh, 1974; Merryman, E. L.
and Levy, A., "Nitrogen Oxide Formation in Flames: The Roles of N0£ and
Fuel Nitrogen", Fifteenth Symposium; Fenimore, C. P., "Ratio of Nitrogen
Dioxide to Nitric Oxide in Fuel-Lean Flames", Combustion & Flame,
volume 1, page 85, 1975.
-------
FIGURE 2-18
PREMDCED FLAT FLAME (METHANE / Am)
104% STOICHIOMETRIC AIR
W
u
tf
w
-------
- 29 -
To give an idea of the results of NO probing at other
stoichiometries, premixed flat flame data for methane/air are plotted
in Figure 2-19. At each stoichiometry on this plot, curvature can be
observed in the first inch or so from the burner face. As the stoichiometry
changes from 111% stoichiometric air to 160% stoichiometric air, the
temperature drop is sufficient to decrease post-flame NOX formation
despite increasing oxygen concentrations. This factor is accentuated
by the unheated walls. At 104% stoichiometric air, NOX is formed more
rapidly near the burner face, but the concentration flattens out after
about two inches so that the concentration at seven inches from the
burner face is lower than that for 111% stoichiometric air. The
combination of lower oxygen concentration in the post-flame zone and the
temperature drop caused by heat losses result in this lower post-flame
concentration.
-------
- 30 -
FIGURE 2-19
PREMISED FLAT FLAME (METHANE/AIR)
250 r
111% STOIC H. AIR
0
0 5 10
DISTANCE FROM BURNER (IN)
-------
- 31 -
Heated walls change the post-flame picture, as would be expected.
With the walls maintained at the adiabatic flame temperature at each
stoichiometry, peak NOX concentrations at 500-msec residence time are
observed at 100% stoichiometric air as shown in Figure 2-20. NOX concentra-
tions appear to approach zero at about 70% and at 150% stoichiometric air
when the walls are maintained adiabatic. However, if the walls are
maintained at a constant temperature of 2220°K, only minor differences
are noted on the fuel rich side (very little post-flame NOX formation) while
major differences appear under fuel lean conditions. Because of the sub-
stantial concentrations of oxygen present under fuel lean conditions, the
high temperature encourages the formation of post-flame NOX and rather than
dropping off, the NO levels increase on the fuel lean side, approaching
2000 ppm at 1507o stoichiometric air and 500 msec residence time. This
indicates that if substantial air preheat were used, high NO levels would
result on the fuel lean side. However, high preheat on the fuel rich side
does not have the same effect in a premixed flame. With interstage heat
removal to allow burnout at lower temperatures, after a fuel-rich first
stage, the NO levels resulting from high-temperature, high-oxygen operation
would be avoided.
2.2.6 Results for Diffusion Flames
Limited data were taken with diffusion flames to investigate
their behavior under a variety of heat loss conditions. Some investigations
were also made into the effect of fuel jet velocity on NOX formation.
The detailed data and centerline species plots are presented
in the Appendix. Only illustrative highlights will be covered in this
section. Under adiabatic wall conditions (wall temperature maintained
at the adiabatic flame temperature) the NO concentrations as a function
of stoichiometry are similar to those for premixed flames far downstream
of the burner face. In Figure 2-21 it can be seen that methane/air exhibits
a peak NO level of 1000 ppm at 100% stoichiometric air for 150 msec residence
time for adiabatic wall condition. Adiabatic flame temperatures vary with
stoichiometry. At 100% stoichiometric air, the adiabatic flame temperature
is about 2220°K, while at 60% stoichiometric air it is 1784°K and at 160%
stoichiometric air it is 1709°K for methane/air combustion. With the walls
heated to 2220K, the higher oxygen level at higher excess air causes more
rapid post-flame NO formation at a constant temperature. For a residence
time of 150 msec, peak NOX is shifted from 100% stoichiometric air to 1201
stoichiometric air when the walls are maintained at 2200K. As the residence
time is increased further towards the 500 msec residence time used with the
flat flame (with walls at 2220°K), the peak would shift toward leaner
mixtures, and the diffusion flame NO curve in Figure 2-21 would look more
and more like the flat-flame NO curve in Figure 2-20.
A comparison between the NOX levels for turbulent diffusion
flames of methane/air and propane/air with adiabatic walls is illustrated
in Figure 2-22. The NOX levels for propane/air are somewhat higher which
could be attributed in part to the higher adiabatic equilibrium temperature
for propane/air and in part to differences in flame chemistry. The curves
are similar in shape and behavior.
-------
- 32 -
FIGURE 2-20
PREMIXED FLAT FLAME (METHANE/AIR)
2500
2000
Q
W
CQ
H
1500
CM
1000
500
ADIABATIC
_OR 2220°K
WALL TE
Q
ADIABATIC WALL
60 80 100 120 140" 160
PERCENT STOICHOMETRIC AIR
180
-------
- 33 -
FIGURE 2-21
TURBULENT DIFFUSION FLAME (METHANE/AIR)
1500
g 1000
CQ
<
W
PH
Q 500
0
WALL = 2220°K
ADIABATIC
WALL
60 80 100 120 140 160
PERCENT STOICHIOMETRIC AIR
180
-------
- 34 -
FIGURE 2-22
TURBULENT DIFFUSION FLAME
1500
0
80 100 120 140 160
PERCENT STOICHIOMETRIC AIR
180
-------
- 35 -
A comparison of results for different fuel jet velocity and
induced turbulence yields some interesting results. Three types of fuel
jets were used in the same basic burner, All fuel jets had strictly
axial injection with no induced swirl. The fuel jet of the stabilized
burner had a 0.250 inch inside diameter which was covered with an end-plate
with six 0.020 inch holes drilled in a circular pattern. The fuel jet
of the large single hole burner had a 0.250 inch inside diameter fuel
jet with no end plate and the small single hole burner had 0.156 inch
inside diameter fuel jet with no end plate. The latter two fuel jets
were sized so that the velocity of the methane from the large jet matched
the velocity of the propane from the small jet at the same stoichiometry.
Table 2-2 shows that for methane/air under excess air conditions
the stabilized jet yields lower NO results than either the large single
hole or small single hole burners and the results drop off more rapidly as
the excess air level is increased. This is probably caused by the rapid
breakup of the multiple individual jets, resulting in behavior similar to
that found in a premixed flame. For the small single hole fuel nozzle the
measured NOX levels are essentially independent of excess air levels on an
"as measured" basis. This behavior seems to indicate that the fuel jet from
this burner maintains its integrity for a long distance downstream so that
stoichiometry plays a smaller role in NO level.
X
Propane/air gives higher NOX levels than methane/air for the
stabilized jet and for the small single hole, as expected from the higher
flame temperature. But the large single hole burner gives much lower NO
levels than either of the other two burners. This fuel/jet combination
gives the lowest jet velocity and may provide an explanation for the
results since in this case, the fuel velocity is smaller than the air
velocity. The low velocity fuel jet would cause lack of penetration of
the fuel into the air stream and recirculation of the surrounding air jet
to mix with the fuel stream near the burner face. This could cause
locally fuel rich combustion to yield results similar to that of staged
combustion. This subject was not pursued extensively at this time since
the major objective of this study was the understanding of the chemistry
rather than the fluid mechanics. However, it is noted as a possible area
for future investigation.
-------
TABLE 2-2
DIFFUSION FLAME NO/NO MEASUREMENTS
NO/NOX, ppm
Methane Propane
Percent Stoichlometric Air Percent Stoichiometric Air
Burner 120 130 140 110 125 150
Stabilized 30/35 16/18 5/10 55/60 25/30
(6 Hole)
Large 40/45 25/30 15/20 35/37 10/15
Single Hole
Small 45/55 45/50 40/45 68 60/65 48/50
i
u>
ON
Measurements taken after 100 msec residence time. ,
Velocity of methane from large single hole burner matches velocity of propane
small single hole burner at same stoichiometry. Air velocity constant.
-------
- 37 -
3. THEORETICAL CALCULATIONS
Theoretical calculations have been undertaken to obtain com-
parisons between experimental data obtained under well-defined conditions
and detailed chemical kinetics calculations. In addition, consideration
has been given to the significance of chemical reactions for which data
are available and the potential significance of reactions for which data
are not available. This section contains a discussion of the thermochemical
and kinetic calculations that have been performed thus far.
3.1 Equilibrium Calculations on Premixed Flames
Equilibrium calculations have been performed using the NASA
CEC71 computer program with the thermochemical information from the
JANAF Thermochemical Tables, Second edition (NSRDS-NBS 37). Equilibrium
compositions for the adiabatic flame temperature at one atmosphere for fuel/
air compositions with initial conditions of 298°K at one atmosphere have been
calculated for methane/air, propane/air, hydrogen/air, and carbon monoxide/
air. The information is presented in Tables 3-1 through 3-4. Only species
for which concentrations exceeded 5 ppm are included in the tables; species
considered in the calculations that did not exceed 5 ppm are listed at the
bottom of the table. Input air is considered to be dry and to contain
appropriate percentages of argon and carbon dioxide.
3.2 Survey of Coupled N0x/Combustion Kinetics
A survey of kinetic data for methane/air combustion was initiated
during this phase of the study. The status of the study at the end of this
phase is described below. An updated version of the survey will be found in
Reference 3-1.
The 25 species considered along with those species eliminated from
primary consideration are listed in Table 3-5 along with the factors that
went into making the species selections. The reactions included in this
portion of the study were those involving the 25 primary species for which
experimental or theoretical information was available in the literature.
A prescreening of the literature indicated that such information existed for
142 reactions. Thus, these reactions were numbered from 1 to 142* on an
interim basis. During the detailed reaction data accumulation it was
determined that there was not sufficient information in the literature to
yield a numerical estimate for some of these reactions (although the informa-
tion on some of them was sufficient to ascertain that the rates were slow).
It was also decided not to maintain separate headings for unimolecular
reactions of differing order. Thus, numerical data were available only on
127 of the 142 reactions of the interim list.
3-1 Engleman, V. S., Survey and Evaluation of Kinetic Data on Reactions
in Methane/Air Combustion, EPA-600/2-76-003, January, 1976.
* This interim numbering differs from the final numbering system in
Reference 3-1 which runs from 1 to 322. A cross-reference will be
found in Appendix A.
-------
TABLE 3-1
EQUILIBRIUM COMPOSITION OF PRODUCTS IN METHANE/AIR COMBUSTION
P=l ATM
INITIAL TEMP. = 298°K
PCT STOICH AIR
TEMP, DEC K
60
1784
80
2097
100
2226
120
2045
140
1861
160
1709
PRODUCTS
AR
CO
C02
H
H2
H20
NO
N02
N2
0
OH
02
PRODUCT MOLE FRACTIONS
0.00709
0.09824
0.03488
0.00011
0.11459
0.15114
0.00000
0.00000
0.59395
0.00000
0.00001
0.00000
0.00787
0.05328
0.05757
0.00057
0.03529
0.18548
0.00006
0.00000
0.65948
0.00000
0.00036
0.00001
0.00837
0.00898
0.08541
0.00039
0.00364
0.18296
0.00197
0.00000
0.70061
0.00021
0.00286
0.00459
0.00856
0.00086
0.07955
0.00004
0.00037
0.15890
0.00333
0.00000
0.71569
0.00016
0.00197
0.03057
0.00867
0.00011
0.06969
0.00000
0.00006
0.13859
0.00264
0.00000
0.72502
0.00005
0.00082
0.05434
0.00875
0.00002
0.06165
0.00000
0.00001
0.12261
0.00183
0.00000
0.73195
0.00001
0.00033
0.07285
MOLE FRACTION LESS THAN 0.000005
CN
CN2 C2
C(S) C CH CH2 CH20 CH3 CH4
C2H C2H2 C2H4 C2N C2N2 C20 C3 C302 C4 C5
HCM HCO HNO H02 H20(S) H20(L) H202 N NH NH2
NH3
N2C N2H4 N20 N204
U)
00
-------
TABLE 3-2
EQUILIBRIUM COMPOSITION OF PRODUCTS IN PROPANE/AIR COMBUSTION
P=l ATM
INITIAL TEMP. = 298°K
PCT STOICH AIR
TEMP DEC K
60
1822
80
2140
100
2254
120
2076
140
1888
160
1732
PRODUCTS
AR
CO
C02
H
H2
H20
NO
N02
N2
0
OH
02
PRODUCT MOLE FRACTIONS
0.00728
0.12308
0.04098
0.00013
0.09554
0.12281
0.00000
0.00000
0.61016
0.00000
0.00001
0.00000
0.00805
0.06366
0.07233
0.00065
0.02699
0.15342
0.00010
0.00000
0.67431
0.00001
0.00047
0.00002
0.00852
0.01183
0.10341
0.00042
0.00316
0.14841
0.00232
0.00000
0.71307
0.00028
0.00301
0.00555
0.00870
0.00132
0.09667
0.00005
0.00037
0.12285
0.00366
0.00000
0.72702
0.00020
0.00205
0.03111
0.00879
0.00017
0.08471
0.00001
0.00037
0.11232
0.00291
0 .00000
0.73510
0.00006
0.00086
0.05501
0.00886
0.00003
0.07484
0.00000
0.00001
0.09926
0.00201
0.00000
0.74096
0.00002
0.00034
0.07367
MOLE FRACTION LESS THAN 0.000005
CN
CN2 C2
C(S) C CH CH2 CH20 CH3 CH4
C2H C2H2 C2H4 C2N C2N2 C20 C3 C302 C4 C5
HCN HCO HNO H02 H20(S) H20(L) H202 N NH NH2
NH3 N2C N2H4 N20 N204
-------
TABLE 3-3
EQUILIBRIUM COMPOSITION OF PRODUCTS IN HYDROGEN/AIR COMBUSTION
P=l ATM
INITIAL TEMP. = 298°K
PCT STOICH AIR
TEMP, DEC K
60
2183
80
2349
100
2383
120
2216
140
2037
160
1884
PRODUCTS
AR
CO
C02
H
H2
H20
NO
N02
N2
0
OH
02
MOLE
C(S)
C2
C5
NH2
0.00625
0.00016
0.00005
0.00219
0.1866,
0.28061
0.00004
0.00000
0.52359
0.00000
0.00045
0.00000
FRACTIONS LESS THAN 0.000005
C CH CH2 CH20 CH3
C2H C2H2 C2H4 C2N C2N2
HCN HCO HNO H02 H20(S)
NH3 N2C N2H4 N20 N204
PRODUCT MOLE
0.00707
0.00013
0.00009
0.00349
0.07954
0.31478
0.00038
0.00000
0.59194
0.00007
0.00239
0.00012
CH4 CN
C20 C3
H20(L) H202
FRACTIONS
0.00762
0.00005
0.00019
0.00181
0.01534
0.32991
0.00260
0.00000
0.63717
0.00055
0.00687
0.00483
CN2 CO
C302 C4
N NH
0.00792
0.00001
0.00024
0.00029
0.00223
0.29115
0.00453
0.00000
0.66067
0.00049
0.00539
0.02707
0.00810
0.00000
0.00026
0.00004
0.00044
0.25801
0.00404
0.00000
0.67661
0.00020
0.00272
0.04957
0.00824
0.00000
0.00026
0.00001
0.00010
0.23054
0.00308
0.00000
0.68872
0.00007
0.00127
0.06770
o
I
-------
TABLE 3-4
EQUILIBRIUM COMPOSITION OF PRODUCTS IN CARBON MONOXIDE/AIR COMBUSTION
P=l ATM
INITIAL TEMP. = 298°K
PCT STOICH AIR
TEMP, DEC K
PRODUCTS
AR
CO
C02
NO
N02
N2
0
02
60
2306
0.00626
0.18884
0.28016
0.00052
0.00000
0.52382
0.00009
0.00031
80
2417
PRODUCT MOLE
0.00704
0.09379
0.30186
0.00292
0.00000
0.58815
0.00070
0.00555
100
2384
FRACTIONS
0.00754
0.04226
0.29690
0.00511
0.00000
0.62909
0.00106
0.01804
120
2285
0.00787
0.01582
0.27914
0.00589
0.00000
0 . 65 637
0.00083
0.03407
140
2151
0.00809
0.00472
0.25497
0.00545
0.00000
0.67485
0.00044
0.05148
160
2009
0.00824
0.00124
0.23023
0.00442
0.00000
0.68785
0.00018
0.06785
MOLE FRACTION LESS THAN 0.000005
C(S) C CN CN2 C2 C2N C2N2 C20 C3 C302
C4 C5 N N2C N20 N204
-------
- 42 -
TABLE 3-5
SPECIES CONSIDERED FOR KINETICS SURVEY ON METHANE/AIR COMBUSTION
SPECIES
C
CH
CHN
CHO
CH2
CH20
CH3
CH30
CH4
CN
CNO
CO
C02
C2
C2H
C2H2
C2H3
C2H4
C2H5
C2H6
C20
C4
C5
PRIMARY
X
X
X
X
X
X
X
X
X
X
X
SECONDARY
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
LOG K
AT 2000
-10
-10
-2
+3
-8
+1
-5
-3
-6
+7
+10
-12
-6
-3
-5
-1
-10
-15
-14
FACTORS INCLUDED IN
CONSIDERATION
possible role in soot formation
hydrocarbon radical
possible role in prompt NO
stable radical
hydrocarbon radical
combustion intermediate
hydrocarbon radical
possible role in ignition
starting material
possible role in prompt NO
possible role in NO/HC interaction
combustion product
combustion product
possible role in soot formation
C2 intermediates possibly
important for
fuel rich combustion of
methane or higher
hydrocarbons
J
difficult to form from CO
"^possible role in soot formation
/more likely from higher
J hydrocarbons
-------
- 43 -
TABLE 3-5 (CONTINUED)
SPECIES
H
HN
HNO
HO
H02
H2
H2N
H2N2
H20
H202
H3N
H3N2
H4N2
N
NO
N02
N03
N2
N20
N204
N205
0
02
Oq
PRIMARY
X
X
X
X
X
X
X
X
X
X
X
X
X
X
SECONDARY
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
LOG K
AT 2000K
-3
-8
-5
0
-3
0
-6
-12
44
-5
-14
-9
-2
-4
-10
0
-6
-15
-18
-3
0
-7
FACTORS INCLUDED IN
CONSIDERATION
combustion intermediate
possible role in prompt NO
possible role in prompt NO
combustion intermediate
combustion intermediate
combustion product
possible intermediate
not likely for methane/air
combustion product
possible role in ignition
possible role fuel rich
not likely for methane/air
not likely for methane/air
important role NO formation
of prime interest
oxidation of NO
higher oxidation state of NO
starting material
possible role in prompt NO
not likely for methane/air
not likely for methane/air
combustion intermediate
starting material
possible role in ignition
-------
- 44 -
The interim recommended rates for the reactions used in this
portion of the study are included here as Table 3-6. More up-to-date
recommendations will be found in Reference 3-1. All reactions included
in the table are covered by rate constants and/or comments. The rates
are also calculated for temperatures of 300°, 1500°, and 2500°K in the
direction indicated (F for forward and R for reverse). Where five asterisks
appear in the comments section, only limited estimates are available and
should be regarded with caution. It should be realized that recommended
rates are evolutionary. As further kinetic and thermochemical data become
available, recommendations require updating. Thus, any recommended rates
should be considered as subject to further evaluation.
Additional reactions can be written involving these species in
unimolecular and bimolecular reactions, some of which could prove to be
significant in the mechanism of methane combustion with pollutant
formation. Thus, a listing has been made of all the possible unimolecular
and bimolecular reactions among 25 species in the methane/air combustion
system. The listing (Appendix A) contains notations indicating the type of
reaction for those included in the survey. An evaluation of the potential
importance of each of these reactions in methane combustion will be found
in Reference 3-1. The 322 reactions in this "master list" are numbered
both with their interim reaction numbers, where applicable, and sequentially,
with "X" numbers which cross-reference this list to the complete survey in
Reference 3-1.
3.3 Kinetics Calculations
Kinetics calculations during this phase of the program were
performed using the version of the CKAP* program modified for use on an
IBM 370 computer. Difficulties were encountered in attempting to
prevent negative concentrations in plug-flow cases. The rate of disappearance
of hydrocarbons is so rapid as their concentrations approach zero that
concentrations of these substances frequently go negative instead of
leveling off near zero. By taking special care with step size control and
relative error criterion, the problem could be alleviated in most cases.
The initial set of calculations for reaction screening and data
comparisons have been run under isothermal conditions at the adiabatic flame
temperature in the region between 80 and 120% stoichiometric air at atmo-
spheric pressure. The isothermal approach offers a number of advantages:
computer time is shortened, screening is facilitated since the rate para-
meters are held constant, and ignition is easily accomplished. Additional
runs have been made for the case of adiabatic ignition.** These runs
indicate the same qualitative behavior as the isothermal runs and basically
similar proportions of reaction intermediates are computed, even though the
absolute concentrations are not identical for the two cases.
* Supplied by EPA for use in this program. Developed for detailed kinetic
analysis by Ultrasystems under EPA contract 68-02-0220.
** Such runs were accomplished by allowing a small amount of initial
reaction to occur at high temperature. After a temperature rise of
about 200°K, the system was dropped back to the lower temperature
(initial temperature + temperature rise) and ignition was allowed to
proceed.
-------
RECOMMENDED RATES
REACTION
IF. CH + CRN = CH2 + CN
2F. CH + CHO = CH2 + CO
3Ft CH + CH20 = CH2 + CHO
4F. CH + CH3 * CH2 + CH2
5F. CH + CH4 - CH2 + CH3
6F. CH + C02 = CHO + CO
7F. CH + H = CH2
8F. CH + H + M = CH2 + M
9R. CH2 + N a CH + HN
10F. CH + HN = CN + H2
11F. CH + HNO = CH2 + NO
TABLE 3-6
INTERIM RECOMMENDED RATES
1500-2500K. 9/30/73
C COMMENTS
8. #**** NOTE B.
BASED ON Tl (J-P)
1. ***** NOTE B.
BASED ON Tl (J-P)
4. ***** NOTE B.
BASED ON Tl (J-P)
5. ***** NOTE B.
BASED ON Tl (J-P)
6, ***** NOTE 8.
BASED ON Tl (J-PI
6. ***** NOTE A*
BASED ON Tl (EST)
0. BASED ON Tl (EST)
USE THIRD ORDER
SEE 8F.
NO DATA
WILL ESTIMATE
11.8 0.67 4.0.5 ***** NOTE A.
BASED ON MS (J-P)
NO DATA
PROBABLY SLOW
9R FASTER THAN 10F
31R FASTER THAN 10R
LOG A B
11.5 0.6
10.5 0.7
11.0 0.7
11.1 0.7
11.4 0.7
10.0 0.5
11.7 0.5
11.8 0.5
***** NOTE A.
BASED ON Tl (EST)
TEMP* KELVIN
300 1500 2000 2500
LOG K
7.2 12.2 12.6 12.8
11.5 12.6 12.7 12.8
9.8 12.6 12.9 13.0
9.2 12.6 12.9 13.0
8.8 12.7 13.1 13.3
6.9 10.7 11.0 11.2
12.9 13.3 13.4 13.4
-16.0 8.0 9.6 10.5
13.0 13.4 13.5 13.5
Ln
I
NOTE A - BEST ESTIMATE AVAILABLE, USE WITH CAUTION
NOTE B - ONLY ESTIMATE AVAILABLE, USE WITH CAUTION
NOTE C - LIMITED EXPERIMENTAL DATA AVAILABLE, USE WITH CAUTION
-------
RECOMMENDED RATES
REACTION LOG A
12F. CH -t- HO = CHO + H 11.7
13R. CH2 + 0 = CH + HO
14F. CH + HO = CO + H2
15F. CH + H02 = CHO + HO 11.7
16F. CH + H02 = CH2 + 02 10.0
17R. CH2 + H = CH + H2
18R. CH2 + HO = CH + H20
19F. CH + N = CN + H
20F. CH -t- NO = CHO + N
21F. CH -i- N2 = CHN + N
22F. CH + 0 = CO + H 11.7
23F. CH + 0 + M « CHO + M 16.
TABLE 3-6 (CONTINUED)
1500-2500K 9/30/73
B C COMMENTS
0.5 10. ***** NOTE A.
BASED ON Tl (EST)
11.3 0.7 26.
***** NOTE B.
BASED ON M8 (J-P)
NO DATA
PROBABLY SLOW
13R FASTER THAN 14F
14R UNFAVORABLE
0.5 6. ***** NOTE A.
BASED ON Tl !EST)
0.5 15. ***** NOTE B.
BASED ON Tl (EST)
11.5 0.7 26,
11.7 0.5
***** NOTE B.
BASED ON M8»T1 (J-P)
***** NOTE A.
BASED ON Tl (EST)
NO DATA
WILL ESTIMATE
NO DATA
PROBABLY SLOW
NO DATA
SPIN FORBIDDEN
0.5 4. ***** NOTE A.
BASED ON Tl (EST)
•0.5 0. ***** NOTE A.
BASED ON Tl (EST)
TEMPt KELVIN
300 1500 2000 2500
LOG K
5.7 11.8 12.3 12.5
-5.9 9.7 10.8 11.4
8.6 12.4 12.7 12.9
0.3 9.4 10.0 10.4
•5.7 9.9 11.0 11.6
8.6 12.4 12.7 12.9
10.0 12.7 12.9 13.0
14.8 14.4 14.3 14.3
i
j^
-------
28Rt CH3 + CN = CHN + CH2
29R. CH4 + CN = CHN + CH3
30R. CHO + CN = CHN + CO
31F. CHN + H = CN + H2
32F. CHN + HO = CN + H20
33R. CN + HN = CHN + N
34R. CN + HNO = CHN + NO
LOG A B C
11.7 0.5 6.
RECOMMENDED RATES
REACTION
24F. -CH + 02 = CHO + 0
25F. CH + 02 = CO + HO
26R. CN + H + M = CHN + M 16.5 -0.5
27R. CH20 + CN = CHN + CHO 11.1 0.7
TABLE 3-6 (CONTINUED)
1500-2500K 9/30/73
11.0 0.7
11.5 0.7
11.3 0.5
0.7 18.
11.3 0.6
11.0 0.5
11.6 0.5
COMMENTS
***** NOTE A.
BASED ON Tl (EST)
NO DATA
PROBABLY SLOW
***** NOTE A.
BASED ON Tl (EST)
***** NOTE A.
BASED ON Tl U-P)
***** NOTE A.
BASED ON Tl (J-P)
*****
BASED ON Tl (J-P)
XPT 100X J-P 300K
REQUIRE A=13.5 OR
REQUIRE E=0.
***** NOTE A.
BASED ON Tl (EST)
***** NOTE B.
BASED ON Tl (J-P)
***** NOTE A.
BASED ON Tl (J-P)
***** NOTE A.
BASED ON Tl (EST!
***** NOTE A.
BASED ON Tl (EST)
TEMP* KELVIN
300 1500 2000 2500
LOG K
8.6 12.4 12.7 12.9
15.3 14.9 14.8 14.8
10.6 12.9 13.1 13.2
10.5 12.8 13.0 13.1
9.6 13.0 13.3 13.4
12.5 12.9 13.0 13.0
0.0 11.0 11.7 12.2
9.1 12.5 12.7 12.9
10.8 12.3 12.4 12.5
12.8 13.2 13.3 13.3
-j
i
-------
TABLE 3-6 (CONTINUED)
RECOMMENDED RATES
REACTION
35F. CHN + 0 = CN + HO
LOG A
11.5
38F. CHO + CH2 = CH3 + CO
40F. CHO + CH3 = CH4 + CO
42R. CH2 + 0 = CHO + H
B
0.7
36F. CHO +M=CO+H+M 20.4 -1.5
37F. CHO + CHO = CH20 + CO 11.2 0.5
10.5 0.7
39F. CHO + CH3 = CH2 + CH20 11.2 0.7
11.5 0.5
41F. CHO + CH4 = CH20 + CH3 11.9
11.7
0.6
0.5
1500-2500K 9/30/73
C COMMENTS
17. ***** NOTE B.
BASED ON Tl (J-P)
16.8 SRI ESTIMATE 4/11/73
BASED ON HINSHELWOOO
-LINDEMANN THEORY
XPTL REF B65 GIVES
A=12.3»B=0.5»E=29.
NOT SENSITIVE XPT
0. ***** NOTE A«
BASED ON Tl (EST)
1. ***** NOTE A.
BASED ON Tl (J-P)
A, ***** NOTE B.
BASED ON Tl (J-P)
0. ***** NOTE A.
BASED ON Tl (EST)
9. ***** NOTE B.
BASED ON Tl (J-P)
44F. CHO + H = CO + H2
10.5 1.
4t
43R. CH20 + M = CHO + H + M 33.9 -4.5 87.
***** NOTE A.
BASED ON Tl (EST)
SRI ESTIMATE 4/11/73
BASED ON HINSHELWOOD
-LINDEMANN THEORY
SRI ESTIMATE 4/6/73
REF Tl ESTIMATES
A=12.2»B=0.5»E=0.
TEMPt KELVIN
300 1500 2000 2500
LOG K
0.8 11.2 12.0 12«4
4.4 13.2 13.6 13.8
12.4 12.8 12.9 12.9
11.5 12.6 12.7 12,
10.0 12.8 13.1 13.2
12.7 13.1 13.2 13.2
6.8 12.5 12.9 13.2
10.0 12.7 12.9 13.0
-40.6 6.9 9.5 11.0
13.0 13.7 13.8 13.9
00
I
-------
RECOMMENDED RATES
REACTION LOG A
45F. CHO + NO = CH20 + NO 11.5
46R. CH20 + 0 = CHO + HO
TABLE 3-6 (CONTINUED)
1500-2500K 9/30/73
C COMMENTS
0. ***** NOTE A.
BASED ON Tl (EST)
47F. CHO + HO « CO + H20
48R. CH20 + H = CHO + H2
49R. CH20 + HO = CHO + H20
50F. CHO + N = CO + HN
51F. CHO + NO = CO -f HNO
52F. CHO + 0 = CO + HO
53F. CH2 + CH4 = CHS + CH3 12.1 0.7
54F. CH2 + HNO = CH3 + NO
11.3 1. 4.4 BASED ON D22»N7»M27
USING B = 1
FIT MODIF ARRHENIUS
D22 BEST HI T
N7 & M27 BEST LO T
10.5 1. 0. SRI ESTIMATE 4/6/73
XPT 14. AT 2000K
OTHER ESTS 12.5
AT 2000K
10.1 1. 3.2 BASED ON W32
USING 8=1
NO DATA ABOVE 1000K
9.5 1. 1. SRI ESTIMATE 4/6/73
AGREES WITH REF Tl
11.3 0.5 2. ***** NOTE A.
BASED ON Tl (EST)
11.3 0.5 2. ***** NOTE A.
BASED ON Tl (EST)
11.5 1. 0.5 SRI ESTIMATE 3/28/73
OTHER ESTS 12.8
20. ***** NOTE A.
BASED ON Tl (J-P)
11.8 0.5 0. ***** NOTE A.
BASED ON Tl (EST)
TEMP, KELVIN
300 1500 2000 2500
LOG K
12.7 13.1 13.2 13.2
10.6 13.8 14.1 14.3
13.0 13.7 13.8 13.9
10.2 12.8 13.1 13.2
11.2 12.5 12.7 12.8
11.1 12.6 12.7 12.8
11.1 12.6 12.7 12.8
13.6 14.6 14.7 14.9
-0.7 11.4 12.2 12.7
13.0 13.4 13.5 13.5
-------
TABLE 3-6 (CONTINUED)
RECOMMENDED RATES
REACTION
55F. CH2 + HO = CH3 + 0
56R« CH3 + H = CH2 + H2
57R. CH3 + HO = CH2 + H20
58F. CH2 + 02 = CH20 + 0
59R. CH3 + 0 = CH20 + H
LOG
A
1500-2500K 9/30/73
B C COMMENTS
300
TEMP.
1500
KELVIN
2000
2500
LOG K
11.
11.
10.
11.
7
3
8
7
0.5
0.7
0.7
0.5
6.
3.
2.
7.
*****
BASED ON
*****
BASED ON
*****
BASED ON
*****
NOTE A.
Tl (EST)
NOTE B.
Tl (J-P)
NOTE B.
Tl (J-P)
NOTE A.
8
10
11
7
.6
.8
.1
.8
12
13
12
12
.4
.1
.7
• 3
12.7
13.3
12.9
12.6
12.9
13.4
13.0
12.8
12.3 0.5
991R. CH30 + M a CH20 + H + M 40.6 -7.5
BASED ON Tl (EST)
-0.3 BASED ON M25»D22
USING B - 0.5
M25 LO T» 022 HI T
22.6 SRI ESTIMATE 4/11/73
BASED ON HINSHELWOOD
-LINDEMANN THEORY
NOTE A.
60R. CH3 + 02 = CH20 + HO
13.5
30,
BASED ON D22
SRI RECOMMENDS
C = 30
D22 GIVES C=10
FOUR CENTER
13.8 13.9 14.0 14.0
5.6 13.5 13.4 13.1
-8.4 9.1 10.2 10.9
61R. CH4 + M = CHS + H + M 15.
17.3
0.
0.
104
62F. CH3 + HNO = CH4 + NO
11.7 0.5
BASED ON H32
FIRST ORDER
87.5 BASED ON H32
SECOND ORDER
RECENT XPT
0. ***** NOTE A.
BASED ON Tl (EST)
•60.8 -0.2 3.6 5.9
-46.4 4.6 7.7 9.7
12.9 13.3 13.4 13.4
Ol
o
-------
RECOMMENDED RATES
REACTION
63R. CH4 + 0 = CH3 + HO
TABLE 3-6 (CONTINUED)
1500-25QOK 9/30/73
LOG A B C COMMENTS
10. 1. 8. BASED ON W25»H25
USING B = 1.
64F. CHS + H02 = CH4 + 02 11.0 0.5
65R. CH4 + H = CH3 * H2 10.7 1.
66R. CH4 + HO = CH3 + H20 13.5 0.
992F. CHS + 02 = CH30 + 0 9.4 1,
67F. CN + NO = CO + N2 11.5 0.
993F. CN + 0 = CO + N 12. 0,
68F. CN + 02 = CO + NO 11.5 0.
69F. CO + HNO = C02 + HN 11.0 0.5
70F. CO + HO = C02 + H 9.6 0.5
6. ***** NOTE A.
BASED ON Tl (EST)
10. BASED ON W30
USING B = 1.
AGREES W/ OTHER
XPT AND EVAL
5. BASED ON W31
EVAL 300-2000K
AGREES G18» 300-500K,
28.5 SRI ESTIMATE 4/11/73
BASED ON EARLIER EST
USING B = 1.
0. ***** NOTE C.
BASED ON B78 (XPT)
0. BASED ON B78» R5
ASSUME C = 0
DATA 300-700K
0. ***** NOTE C.
BASED ON B86 (XPT)
UPPER LIMIT
7. ***** NOTE A.
BASED ON Tl (EST)
0. BASED ON TRANS STATE
PARAMETERS H = -1»
S = -30» CP = -3
LOG K = 11. AT 300K
BETTER FIT TO DATA
WITH A=5.9fB=2.»C=-l
TEMP* KELVIN
300 1500 2000 2500
LOG K
6.6 12.0 12.4 12.7
7.9 11.7 12.0 12.2
5.9 12.4 12.9 13.2
9.9 12.8 13.0 13.1
-8.9 8.4 9.6 10.3
11.5 11.5 11.5 11.5
12.0 12.0 12.0 12.0
11.5 11.5 11.5 11.5
7.1 11.6 11.9 12.1
10.8 11.2 11.3 11.3
Ol
I-1
-------
TABLE 3-6 (CONTINUED)
RECOMMENDED RATES
REACTION
71F. CO + H02 = C02 + HO
72R. C02 + H2 = CO + H20
73R. C02 + N = CO + NO
74F. CO -f N02 = C02 + NO
75F. CO + N20 = C02 + N2
76R. C02 = CO + 0
77R. C02 +M=CO+0+M
78F. CO + 02 = C02 + 0
79R. H2+M=H+H+M
80F. H + HN = H2 + N
LOG A
11.
9.0
11.3
12.5
11.0
1500-2500K. 9/30/73
B C COMMENTS 300
0. 10.5 SRI ESTIMATE 9/11/73 3.4
REF Lll RECOMMENDS
A=12. »B=0.»C=16.5
0.5 15. ***** NOTE A. -0.7
BASED ON Tl (EST)
0.5 25. BASED ON Tl -5.7
WITH HIGHER C
PROBABLY SLOW
SPIN FORBIDDEN
0. 30. BASED ON B47» J6 -9.4
LIMITED OLD DATA
0. 23. ***** NOTE C. -5.8
TEMP» KELVIN
1500 2000 2500
LOG K,
9.5 9.9 10,1
8.4 9.0 9.4
9.2 10.2 10.8
8.1 9.2 9.9
Ul
7.6 8.5 9.0 i
BASED ON L9 (XPT)
USE THIRD BODY
SEE 77R.
15. 0. 100. BASED ON C26» 01
COMPLEX REACTION
SLOW DISSOCIATION
RECENT DATA
13. 0. 60. BASED ON D22
MINIMUM E=50
HI T DATA FAIR AGRMT
14.3 0. 96. BASED ON B88
M = AR
BEST RECENT EVAL
11.8 0.5 8. SRI ESTIMATE
QPR NO. 2. 2/15/73
BASED ON TRANS STATE
PARAMETERS H = 7.
S = -22» CP = -3
-57.9 0.4 4.1 6.3
-30.7 4.3 6.4 7.8
-55.6 0.3 3.8 5.9
7.2 12.2 12.6 12.8
-------
TABLE 3-6 (CONTINUED)
RECOMMENDED RATES
REACTION
81F. H + HNO = HN + HO
82F. H + HNO = H2 + NO
83F. H + HO = H2 + 0
84R. H20 +M=HO+H+M
85F. H + H02 = HO + HO
86F. H + H02 = H2 + 02
87F. H + H02 = H20 + 0
88R. HO + H2 = H + H20
89F.H+N+M=HN+M
90R. HN + 0 = H + NO
LOG A
11.3
13.
9.9
15.5
14.4
13.4
13.0
B
0.5
0*
0.
0.
1500-2500K 9/30/73
C COMMENTS
13. ***** NOTE A.
BASED ON Tl (EST)
2.5 BASED ON H3 1
FIT DATA 200-2000K
7. BASED ON B88
CRIT EVAL 400-2000K
105. BASED ON B88
M = N2
CRIT EVAL 2000-6000K
DIFFICULT TO STUDY
1.9 BASED ON B88
CRIT EVAL 290-800K
0.7 BASED ON 888
CRIT EVAL 290-800K
13.4 0.
1. BASED ON Lll
EVAL 300-1000K.
LOWER THAN REC B88
12.9 AT 293K
5.2 BASED ON B88
EVAL, AGREES XPT
16.5 -0.5 0. ***** NOTE A.
BASED ON Tl
-------
RECOMMENDED RATES
REACTION
91R. HO + N = H + NO
92F.H+NO+M= HNO + M
93R. HNO + 0 = H + N02
94F. H + N02 = HO + NO
95R. HN + N * H + N2
96F. H + N20 = HN + NO
97R. HNO + N = H + N20
TABLE 3-6 (CONTINUED)
1500-2500K 9/30/73
LOG A B
11.8 0.5
16.3 Oi
10.7 0.5
14.5 0.
11.8 0.5
11.0 0.5
10.7 0.5
C COMMENTS
8. SRI ESTIMATE
QPR NO. 2» 2/15/73
BASED ON TRANS STATE
PARAMETERS H = 7
S = -22» CP = -3
K = 13.5 AT 320K
BASED ON REF C8
WOULD REQUIRE C = 0
ALSO USE C=0 FOR SCR
GIVES 13.45 AT 2000K
0. BASED ON C20. HlO
AT 300K
USING B = 0» C * 0
M = H2
DATA 200-700K INDIC
SLIGHT NEGATIVE E
0. ***** NOTE A.
BASED ON Tl (EST)
1.5 BASED ON B90
300-600K
NO DATA ON REVERSE
0. SRI ESTIMATE
QPR NO. 2* 2/15/73
BASED ON TRANS STATE
PARAMETERS H = -1
S = -22» CP = -3
30. ***** NOTE A.
BASED ON Tl (EST)
3. ***** NOTE A.
BASED ON Tl (EST)
TEMP. KELVIN
300 1500 2000 2500
LOG K
7.2 12.2 12.6 12.8
16.3 16.3 16.3 16.3
11.9 12.3 12.4 12.4
13.4 14.3 14.3 14.4
13.0 13.4 13.5 13.5
-9.6 8.2 9.4 10.1
9.8 11.9 12.0 12.1
Ul
I
-------
TABLE 3-6 (CONTINUED)
RECOMMENDED RATES
REACTION LOG A B
98F. H + N20 = HO + N2 13.9 0.
99F. H+0+M=HO+M 15.9 0.
100R. HO + 0 - H + 02 13.4 0.
1500-2500K
9/30/73
10IF. H + 02 + M - H02 + M 15.2 0,
102F. HN -i- HN a H2 + N2
103F. HN + HO = H2 + NO
104F. HN + HO = H20 •*• N 11.7 0.5
105R. HNO + N = HN + NO 11.0 0.5
106F. HN + N02 = HNO + NO 11.3 0.5
107F. HN + N20 = HNO + N2 11.0 0.5
C COMMENTS
15. BASED ON 890
CRIT EVAL 700-2500K
0. BASED ON S30
LIMITED DATA
REF B88 GIVES NO REC
0. BASED ON B88» W31
NOTE REF 688 REC'S
FOR 100F. A= 14.4.
B = 0.» C = 16.8
CRIT EVAL 300-2000K
1. BASED ON B88
M = AR
FOUR CENTER
PROBABLY SLOW
D8 GIVES 14 AT 2000K,
BAHN CITES 13 FROM
REF M8 (ERROR)
FOUR CENTER
PROBABLY SLOW
2. ***** NOTE A.
BASED ON Tl (EST)
2. ***** NOTE A.
BASED ON Tl (EST)
5. **** NOTE A.
BASED ON Tl (EST)
3. ***** NOTE A*
BASED ON Tl (EST)
TEMP* KELVIN
300 1500 2000 2500
LOG K
3.0 11.7 12.3 12.6
15.9 15.9 15.9 15.9
13.4 13.4 13.4 13.4
14.5 15.1 15.1 15.1
11.5 13.0 13.1 13.2
10.8 12.3 12.4 12.5
8.9 12.2 12.4 12.6
10.1 12.2 12.3 12.4
Ul
Ui
-------
RECOMMENDED RATES
TABLE 3-6 (CONTINUED)
1500-2500K 9/30/73
REACTION
108F. HN + 0 = HO + N
109F. HN + 0 + M = HNO + M
110R. HNO + 0 = HN + 02
111R. H02 + N = HN + 02
113F. HNO + HO = H20 + NO
114F. HNO + NO = HO + N20
115F. HNO + 0 = HO + NO
ll6Rt H2 + 02 = HO + HO
LOG A B
11.8 0.5
16.0 -0.5
11.0 0.5
11,
Oi
112R. HO + N + M = HNO + M 15.0 -0.5
14.0
12.
0.
0.
c
8.
0.
1.
12.3 0. 26<
11.7 0.5 0,
12.4 0. 39.
COMMENTS
SRI ESTIMATE
QPR NO. 2, 2/15/73
BASED ON TRANS STATE
PARAMETERS H = 7»
S = -22» CP * -3
***** NOTE A.
BASED ON Tl (EST)
***** NOTE A.
BASED ON Tl (EST)
***** NOTE C.
BASED ON K25 (XPT)
LOWER LIMIT
***** NOTE A.
BASED ON Tl (EST)
BASED ON H31
SRI ESTIMATE 3/28/73
BASED ON W23
REEXAMINATION OF
EXISTING H2-NO DATA
1000-1300K
***** NOTE A.
BASED ON Tl (EST)
BASED ON R2
1400-2500K
REF B88 INDICATES
R2 PROBABLY VALID
BUT MAKES NO REC
TEMP, KELVIN
300 1500 2000 2500
LOG K,
7.2 12.2 12.6 12.8
14.8 14.4 14.3 14.3
7.1 11,6 11.9 12.1
11.0 11.0 11.0 11.0
13.8 13.4 13.3 13.3
14.0 14.0 14.0 14.0
11.3 11,9 11.9 11.9
-6.6 8.5 9.5 10.0
12.9 13.3 13.4 13.4
-16.0 6.7 8.1 9.0
I
Ul
-------
RECOMMENDED RATES
REACTION
117F. HO + HO = H20 + 0
118F. HO + H02 = H20 + 02
119F. HO + NO = H02 + N
120R. H02 + NO = HO + N02
121R. H02 + M = HO + 0 + M
122R. H02 + 0 = HO + 02
123F. H2 + N02 = H20 + NO
124R. N2+M=N+N+M
125F. N + NO = N2 + 0
TABLE 3-6 (CONTINUED)
1500-2500K 9/30/73
LOG A B C
12.8 0. 1.
13. 0. 1.
12,
15.6 -1.5 64.
13.7 0,
21.6 -1.6 225,
13.2 0.
11.8 0.5
0.
0.
COMMENTS 300
BASED ON B88 12.1
EVALUATION 300-2000K
BASED ON Lll
H37 INDIC 14 AT 300K
B88 SUGGESTS 12.8 AS
LOWER LIMIT AT 300K
NO DATA
WILL ESTIMATE
SRI ESTIMATE 3/28/73
LOWER THAN J12.N6
WITHIN LIMITS Lll
SRI ESTIMATE 4/6/73
NO VALID DATA (888)
BASED ON Lll
300-1000K
LIMITED INFO
SLOW
K=3.5 AT 700K
BASED ON B90
M = N2
TEMP* KELVIN
1500 2000 2500
LOG K
12.7 12.7 12.7
12.3 12.9 12.9 12.9
9.8 11.6 11.7 11.7
-34.7 1.5 3.7 4.9
13.0 13.6 13.6 13.6
-146.3 -16.3 -8.3 -3.5
BASED ON B90 13.2
CRITICAL EVALUATION
SRI ESTIMATE 13.0
QPR NO. 2» 2/15/73
BASED ON TRANS STATE
PARAMETERS H = -1»
S = -22« CP = -3
13.2 13.2 13.2
13.4 13.5 13.5
I
Ul
-------
TABLE 3-6 (CONTINUED)
RECOMMENDED RATES
REACTION
126Ft N 4- NO 4- M = N20 + M
127F. N 4- N02 = NO + NO
1500-2500K
128F. N 4- N02 = N2 4- 0 + 0
129F. N 4- N02 = N2 + 02
130F. N 4- N02 = N20 + 0
131F. N 4- N20 «= NO + N2
132R
133F. N 4- 02 = NO 4- 0
LOG A
12.6
12.2
12.0
12.7
8.7
10.
20.6 -Ii5 150.
148,
9/30/73
COMMENTS
NOT LIKELY
GOES TO
STABILIZ
PROBABLY
N2 4- 0
S30-NO EVIDENCE
BASED ON P10
USING B = 0» C =
B90 INDICATES
P10 UPPER LIMIT
PROB 2.5X LOWER
BASED ON P10
USING B = 0» C *(
UPPER LIMIT
BASED ON P10
USING B - 0» C
PROBABLY UPPER
= 0
LIMIT
BASED
USING
ON
B
P10
= 0»
9*8
11.8
1.
0.5
6.3
8.
*#*##
BASED ON B19 (SEL)
BASED ON B90
M = AR» N2» 02
CAUTION RECOMMENDED
RASED ON M26» M « AR
(0) FORMATION MEAS.
BASED ON B90
EVALUATION 300-3000K
SRI ESTIMATE
QPR NO. 2» 2/15/73
BASED ON TRANS STATE
PARAMETERS H = 7,
S = -22. CP * -3
TEMP* KELVIN
300 1500 2000 2500
LOG K,
12.6 12.6 12.6 12.6
12.2 12.2 12.2 12»2
12.0 12.0 12.0 12.0
12.7 12.7 12.7 12.7
1.4 7.2 7.6
7.8
-92.4 -6.0 -0.7 2.4
-93.7 -7.5 -2.1 1.2
7.7 12.1 12.4 12.6
7.2 12.2 12.6 12.8
i
00
-------
TABLE 3-6 (CONTINUED)
RECOMMENDED RATES
REACTION
134F. NO + NO = N2 + 02
135R. N20 + 0 * NO + NO
136F. NO + N02 = N20 + 02
137F. NO + N20 = N02 + N2
138R. N02 +M=NO+0+M
139R. N02 + 0 = NO + 02
140R. N20
141R. N20 + 0 * N2 + 02
142Rt 02+M=0+0+M
LOG A B
14.0 0*
12.0
12.3
16.0
13.
14.
14.
18.4
19.4
0*
0.
0*
-1.
-1.
1500-2500K 9/30/73
C COMMENTS
SPIN FORBIDDEN
PROBABLY SLOW
OLD DATA DID NOT
MEASURE THIS RX
ACTUALLY MEAS 135F
28. BASED ON B90
EVAL 1200-2000K
SAME RATE AS 141R
NO TRANSITION STATE
NO EXPTL EVIDENCE
60. BASED ON B19 (EST)
40. BASED ON K30
USING A = 12.3
NO RECENT DATA
65. BASED ON B90» T4
M = AR* 02
1400-2400K
1. BASED ON B53
CRIT EVAL 300-600K
NO HI T DATA
50. APPROXIMATE
COVERS LIMITED DATA
IN 1500-2500K RANGE
28. BASED ON B90
REC 1200-2000K
SAME RATE AS 135R
118.7 BASED ON J4» M = AR
118.7 BASED ON J4i M = 02
CRITICAL EVALUATION
WIDE TEMP RANGE
1000-8000K FOR 02
300-15000K FOR AR
TEMP» KELVIN
300 1500 2000 2500
LOG K
-6.4 9.9 10.9 11.6
-31»7 3.3 5.4 6.8
-16.8 6.5 7.9 8.8
-31.4 6.5 8.9 10.3
12.3 12.9 12.9 12.9
-22.4 6.7 8.5 9.6
-6.4
-70.6
-69.6
9.9
•2.1
-1.1
10.9
2.1
3.1
11.6
4.6
5.6
I
Ul
-------
- 60 -
In the screening studies, reactions could be eliminated for one
of three reasons:
1. The reaction itself is too slow to play a major
role for any species.
2. The reaction is part of a loop which merely returns
its starting material.
3. A species involved in the reaction is of no
importance in either combustion or pollutant
formation.
These practical considerations must be handled with some caution since
reaction elimination based on species could vary with the goals of the
screening procedure. For example, calculations oriented primarily toward
heat release and combustion efficiency could easily ignore NO forming
reactions since these reactions play a very small role in that respect.
At peak NO formation rate, the NO forming reactions are involved in only
about 0.017= of the 0 atom appearance or disappearance reactions. Since NO is
of prime interest in these calculations, such reactions must be carefully
retained. On the other hand, since soot is not likely to form in methane/air
combustion unless substantial concentrations of hydrocarbons remain at the
end of the combustion zone, reactions precursor to soot formation can probably
be ignored with a fair degree of safety as long as calculations are restricted
to the region of 807o-1207o stoichiometric air at atmospheric pressure.
While progress toward the definition of a "minimum reaction set"
is an ultimate goal it is more important initially to appreciate the chemistry
of the coupled reactions. With 322 possible unimolecular and bimolecular
reactions which can be written for the 25 species in the kinetics survey,
it is entirely possible that some reaction beyond those 142 reactions for
which rate information was reported in the literature might be of importance.
While many of the reactions not reported in the literature are not elementary
reactions an evaluation of the potential importance of each of the remaining
reactions will be undertaken as a step towards the complete evaluation of
the combustion/pollutant formation mechanism for the 25 species considered.
-------
- 61 -
3.4 Results of Theoretical Calculations for Methane/Air
Screening calculations have been performed, using the interim
(September 1973) recommended rates for the kinetics survey, to gain insight
into the major reaction paths and major reactions in the coupled combustion/
pollutant formation for methane/air. The interim rates have since been
superseded by the recommended rates given in Reference 3-1 but the results
obtained from these calculations are instructive for the overall features
uncovered. However, the reader should be aware that the quantitative details
would be different if the revised rates were used and that the reaction set
used for these calculations was not complete. Therefore potentially impor-
tant paths for NO formation have been omitted from these calculations.
x
Most of the calculations were accomplished under isothermal con-
ditions at the adiabatic flame temperature as discussed earlier. These iso-
thermal calculations result in temperatures and reaction rates in the early
combustion zone which are higher than those observed experimentally. There-
fore, the calculations indicate a hydrocarbon burnout rate which is more
rapid than for the case of an adiabatic flame. However, the essential
features of the combustion are preserved as indicated by comparison with
calculations of adiabatic ignition from room temperature. The concentrations
calculated for the intermediates remain basically the same. The concentra-
tion levels are sufficiently close for screening purposes although the over-
shoots tend to reach higher concentrations (but for shorter times) in the
isothermal case.
Concentrations of selected species as a function of time for plug
flow reaction at stoichiometries of 80%, 100% and 120% stoichiometric air
are given in Tables 3-7, 3-8 and 3-9 respectively. It should be emphasized
that these are only preliminary calculations based on the limited number
of reactions reported in the literature. Such calculations are highly
informative but must be followed by more complete calculations for further
elucidation. The reaction times selected for these tables are illustrative
of some of the zone divisions to be discussed in Section 3.5. At this point
a discussion of the species concentrations themselves will be presented.
3.4.1 80% Stoichiometric Air
In Table 3-7 for 80% stoichiometric air, the first column of
concentrations indicates the initial conditions corresponding to zero time.
After 35 psec (assuming no ignition delay) the methane is essentially gone.
The CO has built up to essentially its maximum concentration and some of
the CO has already reacted to form C02. Before CO reaches its peak,
substantial concentrations of hydrocarbon intermediates as well as CHO and
CH20 are observed but these have dropped off by the time CO has peaked.
Concentrations of the radical intermediates H, HO, and 0 reach their peaks
at about the same time as CO, with H being 1.2 mole %, HO being 1.1% and
0 being 0.4%. Molecular hydrogen has also reached a high level (2.3%)
although water is much higher at 18.0%. NO has not yet formed to a
substantial extent, being on the order of 0.2 ppm but its rate of formation,
-------
- 62 -
as indicated in Table 3-7, is quite high. Two species (N and N02) present
in very low concentrations have quite different orders of importance.
Atomic nitrogen plays a major role in NO formation* through N + 02 and
N + OH and it is formed almost exclusively by the reaction N2 + 0 —> NO + N.
Thus N concentration is essentially determined by the 0 concentration.
On the other hand, N02 never builds up to any substantial extent in the
combustion zone and, almost as quickly as it is formed, it returns to NO.
After 600 psec CO and C02 have almost reached their ultimate
concentrations. Species H, HO and 0, while having decreased in concentration,
are still substantially above their equilibrium values. The concentrations
of H2 and H20 have essentially leveled off while NO is still increasing
but at a much slower rate.
After 4.6 msec, species concentrations remain fairly constant,
and although 0 and N atom concentrations are still somewhat above equilibrium,
they are not high enough to cause substantial incremental NO formation.
3.4.2 100% Stoichiometric Air
Table 3-8 summarizes the kinetics calculations for 100%
Stoichiometric air. After 12 psec (no ignition delay) methane has
disappeared and CO has reached its peak. Some of the CO has already reacted
to C02- Concentration of H2 has reached 2% and H20 has reached 13.6%.
Atomic nitrogen is over 0.1 ppm, a high mole fraction for this intermediate,
while NO has just started to form, being slightly in excess of 1 ppm.
Molecular oxygen has dropped from its initial concentration of 19% to
about 4% on its way to 0.5%. The rate of NO formation, at this point, is
extremely high (320,000 ppm/sec) but is dropping off quite rapidly as the
concentration of oxygen atom decreases. By the time 300 ysec have elapsed
the rate of NO formation has dropped to 90,000 ppm/sec and about 50 ppm of
NO has formed.
After 10 msec, CO and C02 concentrations have leveled off at
about 1% and 8% respectively. The concentrations of radical intermediates
have dropped off to within 150% to 200% of their equilibrated values and
the NO formation rate has dropped off to 7,000 ppm/sec. The calculated
concentration of 200 ppm represents an average formation rate of 20,000
ppm/sec, compared to an average formation rate of 150,000 ppm/sec up to
0.3 msec, and an average formation rate of 5,500 ppm/sec up to 100 msec.
At the 100 msec point the instantaneous formation rate has dropped to less
than 4,000 ppm/sec. The rapid initial formation of NO, even though probably
overpredicted because of the high initxai temperature, bears many resemblances
to what might be called "prompt NO" (3-2) .
* It should be noted that recommended rates were not available at the
time of these calculations for additional reactions, which are potentially
important for NO formation. Therefore, the inferences made here should
be regarded as preliminary and subject to further refinement.
(3^2) Fenimore, C. P., 13th Symposium (International) on Combustion,
p. 373, The Combustion Institute, 1971.
-------
- 63 -
TABLE 3-7
KINETIC CALCULATIONS FOR PLUG FLOW CH./AIR COMBUSTION AT 807» STOICH AIR
—4
T = 2097°K
Species
CH4
CO
C02
H
HO
H2
H20
N
NO
N02
0
°2
d(NO)
dt
Branching Zone
t = 0
11.6%
—
—
—
—
—
—
—
—
--
—
18.6%
—
35 y sec
14 ppm
9.4%
1.3%
1.2%
1.1%
2.3%
18.0%
_2
2.7 x 10 ppm
0.2 ppm
-4
10 ppm
0.4%
1.4%
32,000 ppm/sec
Relaxation Zone
600 \i sec
—
5.4%
5.7%
0.4%
0.2%
3.5%
18.3%
_2
1.0 x 10 ppm
4 . 8 ppm
10 ppm
0.02%
0.04%
1,600 ppm/sec
4.6 msec
—
5.3%
5.8%
0.07%
0.05%
3.5%
18.6%
1.6 x 10 ppm
5 . 7 ppm
10 ppm
0.0007%
0.002%
59 ppm/sec
Post-Flame Zone
10 msec
—
5.3%
5.8%
0.06%
0.04%
3.5%
18.6%
1.2 x 10 ppm
5 . 9 ppm
5 x 10 ppm
0.0004%
0.001%
36 ppm/sec
Note: See Section 3-5 for discussion of "zones".
-------
- 64 -
TABLE 3-8
KINETIC CALCULATIONS FOR PLUG FLOW CH./AIR COMBUSTION AT 100% STOICH AIR
T = 2222°K
Species
CH4
CO
C°2
H
HO
H2
H20
N
NO
N02
0
°2
d(NO)
dt
Branching Zone
t = 0
9.5%
—
—
—
—
—
—
—
—
—
—
19.0%
—
12 y sec
5 x 10~ ppm
7.7%
1.1%
2.0%
2.0%
2.0%
13.6%
0 . 11 ppm
1 . 4 ppm
0.0004 ppm
1.4%
4.1%
320,000 ppm/sec
Relaxation Zone
300 y sec
—
3.1%
6.1%
0.7%
1.1%
1.4%
16.1%
0.06 ppm
49 ppm
0.02 ppm
0.4%
1.9%
90,000 ppm/sec
.10 msec
—
1.0%
8.4%
0.06%
0.4%
0.4%
18.2%
0.02 ppm
200 ppm
0.04 ppm
0.03%
0.6%
7,000 ppm/sec
Post-Flame Zone
100 msec
—
0.8%
8.6%
0.04%
0.3%
0.3%
18.4%
0.01 ppm
550 ppm
0.09 ppm
0.02%
0.5%
3,700 ppm/sec
Note: See Section 3.5 for discussion of "zones".
-------
- 65 -
TABLE 3-9
KINETIC CALCULATIONS FOR PLUG FLOW CH,/AIR COMBUSTION AT 1207, STOICH AIR
—— — 4
T = 2044°K
Species
CH4
CO
co2
H
HO
H2
H20
N
NO
N02
0
°2
d(NO)
dt
Branching Zone
t = 0
8.0%
—
—
—
—
—
—
—
—
—
—
19.3%
—
12 y sec
8 ppm
6.7%
0.8%
1.6%
1.6%
1.3%
12.2%
0.02 ppm
0 . 3 ppm
—
1.5%
6.3%
85,000 ppm/sec
Relaxation Zone
900 u sec
. —
0.9%
7.1%
0.2%
0.7%
0.4%
15.1%
0.008 ppm
22 ppm
—
0.2%
3.6%
12,000 ppm/sec
40 msec
—
0.09%
7.9%
0.004%
0.2%
0.04%
15.9%
0.0009 ppm
86 ppm
—
0.02%
3.2%
1,000 ppm/sec
Post-Flame Zone
500 msec
—
0.08%
7.9%
0.004%
0.2%
0.04%
15.9%
0.0009
483 ppm
—
0.02%
3.2%
900 ppm/sec
Note: See Section 3.5 for discussion of "zones".
-------
- 66 -
3.4.3 120% Stoichiometric Air
At 120% Stoichiometric air (Table 3-9) the basic features of
the concentration profiles are the same as with lower air supply, except
that the N atom concentrations are lower and NO is produced more slowly.
At the 12 ysec point it can be seen that the concentrations of H, HO, and 0
are almost the same as in the 100% Stoichiometric air case, yet the N atom
concentration is lower by a factor of 5. As will be discussed later,
N atoms are produced primarily by the reaction N£ + 0 ^ NO + N which
has a 75 kcal activation energy and, therefore, proceeds more slowly at
2044°K than it does at 2222°K. By the same token, NO is produced more
slowly, and again its concentration is lower by a factor of 4 to 5 throughout.
If the temperature were increased to 2222°K, one might expect more rapid NO
formation at 120% Stoichiometric air than at 100% Stoichiometric air.
Indeed, this is observed experimentally (see Section 3.5).
3.5 Discussion
Theoretical calculations based on the recommended rates from
the kinetics survey of methane combustion indicate that, at all mixture
ratios studied, the combustion process under plug flow conditions can be
divided into three zones. For identification and convenience, these zones
will be called the branching zone, the relaxation zone and the post-flame
zone. The basic charapteristics of each zone are as follows: (it should
be noted that the reaction times, provided for reference would vary with the
temperature profile. The features discussed are valid for the general
combustion case.)
(1) The Branching Zone (10-100 ysec duration following ignition
delay, if any):
- is the zone in which the hydrocarbon is consumed.
- the hydrocarbon goes primarily to CO and H20 although some
C02 and H2 are formed in this zone.
- the most important feature of this zone is that the radical
intermediates such as 0, H, HO, etc., build up to extremely
high concentrations and N atoms are substantially in excess
of equilibrium although at low absolute concentrations.
- the individual intermediates do not appear to reach maximum
concentrations at precisely the same time, but the point at
which methane concentration has reached a low level (on the
order of 1 ppm) coincides roughly with the peak concentrations
of the intermediates.
- because of the short times involved, the NO concentration at
the end of the branching zone is quite low (on the order of
1 ppm) but the rate of NO formation is at its peak because
of the high concentrations of the intermediates.
- while very little NO is produced in this zone, it is the high
concentrations of intermediates which set the stage for rapid
NO formation and, therefore, the kinetic details of this zone
are extremely important.
-------
- 67 -
(2) The Relaxation Zone (5-50 msec duration):
- is the zone in which the concentrations of intermediates
approach their ultimate values at the specified temperature.
- the end of this zone is more difficult to define because the
approach to equilibrium is asymptotic. However, the precise
definition is not critical since it is an arbitrary
separation to distinguish this zone from the post-flame
zone in which intermediates have essentially reached
equilibrium. One possible definition of the end of the
zone is the point at which intermediate concentrations
(0, OH, H, N) have reached, e.g., 200% (or 150% or 110%)
of their equilibrium values. Another definition might
involve the rate of change of concentration of these
intermediates .
- in this zone, the rapid formation of NO takes place.
However, the stage has already been set by the high
concentrations of intermediates produced in the branching
zone.
the primary NO forming reactions in this zone, using the
preliminary rate survey, appear to be N£ + 0, N + 02 and
N + HO. However, the concentrations of these species
are not yet equilibrated in this zone. Thus, the rate of
formation of NO can be substantially higher than that
predicted by classical Zeldovich kinetics.
- even without "prompt NO" reactions, which may involve
nitrogeneous species other than atomic and molecular
nitrogen, the high concentrations of intermediates do
result in very rapid or "prompt" NO formation.
- it is still worth searching for "prompt NO" reactions.
The master set of 322 reactions includes several reactions
for which rate data are currently unavailable. An evaluation
of the potential importance of each of the 322 reactions would
provide a starting point for a further screening effort.
With reasonable estimates for the resulting set, the screening
framework would provide a basis for the evaluation of reactions
of importance.
- the formation of N atoms in these calculation is found to
result almost exclusively from the reaction N2 + 0 = NO + N.
The destruction of N atoms is primarily through NO forming
reactions N + 02 = NO + 0 and N + HO = NO + H. Other N-forming
reactions should be sought in the extension of this kinetics
survey.
- the formation of 0 atoms is found to result almost exclusively
from the reaction H + 02 = HO + 0.
-------
- 68 -
the formation of H atoms is found to result, in good part,
from the reaction CO + HO = C02 + H. The destruction of H
atoms occurs, in good part, through the reaction H + 02 =
HO + 0.
- HO is formed in the H + 0? reaction and destroyed in CO + HO
(although not exclusively).
- to summarize the nature of the rapid NO formation observed
in kinetic calculations using the preliminary rate survey:
+ the NO formed depends primarily on the concentrations of
N and 0 atoms.
4- N atom formation depends on the 0 atom concentration
through N£ + 0 = NO + N. N atom destruction depends
on HO + N = NO + H and N + 02 = NO + 0.
+ 0 atom formation depends on the H atom concentration
through H + 02 = HO + 0. 0 atom destruction is
primarily through reaction with CHO at the branching/
equilibration interface and through reaction with
^0 at the equilibration/post flame interface.
+ the H atom formation depends on the CO and HO concen-
trations through CO + HO = CO, + H. Destruction of
H atoms occurs through recombination with HO and CHO
and by the H + C>2 reaction.
+ the initial concentrations of CO, HO, H and 0 in the
equilibration zone are determined by reactions occurring
in the branching zone.
(3) The Post-Flame Zone (50 usec and beyond):
- the concentrations of intermediates are either equilibrated,
or nearly so, and the modified Zeldovich treatment (including
steady state assumptions) is a good approximation of the
behavior in this zone.
The above observations lead one back to questions about a quasi-
global approach in which the hydrocarbon can partially oxidize in one step
and the intermediates can burn out in a detailed manner. The behavior of
the hydrocarbon in the branching zone, at first glance, appears to be
quasi-global in nature. Hydrocarbon reacts very rapidly to produce CO and
H20 before a substantial concentration of NO is produced. However, in
addition to CO and ^0, the hydrocarbon disappearance results in the
appearance of high concentrations of other reactive intermediates.
In the Fourteenth Symposium (International) on Combustion, a
paper by Engleman _et_ al. (3-3) demonstrated that while CO and H2 combustion
in a jet-stirred combustor could be modeled by detailed kinetics, a
quasi-global model for propane was not adequate to predict NO formation.
The global hydrocarbon step in that mechanism was
3-3 V. S. Engleman, W. Bartok, J. P. Longwell, and R. B. Edelman,
Fourteenth Symposium (International) on Combustion, p. 755,
The Combustion Institute (1973).
-------
- 69 -
(I) HC + 02 > CO + H2
The CO and H2 were then allowed to complete combustion by a detailed mech-
anxsm It appears from detailed modeling that the global hydrocarbon step
is actually
(II) HC + 02 > CO + H20 + (reactive intermediates)
(H, HO, H2, 0)
In global step I an approximation to the true picture is made by assuming
all the hydrogen in the hydrocarbon goes through H2, a species which can
participate actively in branching reactions.
Thus, the quasi-global mechanism for hydrocarbon combustion
achieves an approximation of the overshoot of equilibrium by making the
assumption that the hydrocarbon is converted in one step to CO and H£
(see Appendix C). The H2 then participates in branching reactions which
produce high concentrations of intermediates similar to those found in
detailed calculations. However, detailed calculations indicate that the
actual mechanism produces high concentrations of active intermediates
while the hydrocarbon is converted to CO and ^0 (see Tables 3-7, 3-8,
and 3-9 for species concentrations at the end of the branching zone).
These calculations indicate that the quasi-global mechanism does not
adequately predict the overshoot.
3.6 Reactions of Importance at the Interface
Between the Branching Zone and the Relaxation Zone
The interface between the branching zone and the relaxation
zone is an extremely important region from the standpoint of rapid formation
of NO in the flame zone. This interface has been defined in a previous
section as the location where certain combustion intermediates have reached
their peak concentrations. Reactions of importance for formation and
destruction of these intermediates and for NO will be discussed in this
section. It should be reiterated here (and will be discussed further in
the next section) that the events which occur within the branching zone
are very important for the establishment of the high concentrations of
intermediates. The screening that has been done thus far has been more
concerned with determining controlling reactions and less concerned with
the absolute magnitude of rate constants or the concentrations predicted
by the calculations. However, the balance among the reactants, intermediates,
products has been of major concern, since the relative concentrations help
determine the magnitude and direction of the individual reactions.
At the end of the branching zone, NO is just barely starting
to form in ppm level concentrations. At this point the concentrations
of intermediates are high and the rate of NO formation is at its peak.
It is instructive to consider some of the important reactions at this point.
It should be emphasized that while the comments in this section are based
on detailed calculation using a large number of reactions, not all possible
reactions have been included because rate information was unavailable in
-------
- 70 -
some cases. The inclusion of "prompt NO" reactions in these calculations
could result in changes in the NO mechanism under specific conditions.*
Most of the observations discussed below, however, are expected to be valid.
At the point where carbon monoxide starts to be destroyed more
rapidly than it is being produced, the primary destruction reaction is
CO + HO > C02 + H
which accounts for 70-80% of its destruction. This reaction also produces
most of the H at this point although within the branching zone most of the
H comes from hydrocarbon fragments.
The H atoms, at this point, are being destroyed by the reaction
H + 0 • » HO + 0
which accounts for 40-50% of its destruction and also produces essentially
all of the 0 atoms. H is also being destroyed (25-30%) by
CO + H + M *• CHO + M
and (10-15%) by
H + H20 —>• HO +
The 0 atoms, as indicated above are formed entirely* * by
H + 02 - » HO + 0
at the interface, and are destroyed (45-50%) by
CHO + 0 •*•••> CO + HO
and (25-35%) by***
KLO + 0 —*- HO -1- HO
* "Prompt NO" reactions, in this context, are those involving hydrocarbon
fragments which result in production of N atoms; e.g. CH + N~ ->N.
** In the branching zone itself, some 0 atoms are formed by HO + HO ->H-0 4- 0.
*** In part of the branching zone this reaction is a producer of oxygen atoms.
-------
- 71 -
Based on calculations which do not consider the role of hydrocarbon
fragments in the production of N atoms, the following observations are noted.
While 0 atoms are extremely important for the production of NO, the NO
reactions are relatively unimportant in determining 0 atom concentrations.
NO reactions account for less than 1% of the total 0 atom production and
destruction. The main NO-producing reaction at the interface is
N2 + 0 •- >- NO + N
which produces about half of the NO and essentially all of the N atoms.
The reaction
N + 0 >• NO + 0
produces 25-30% of the NO and the reaction
HO + N -—». H + NO
produces about 20% of the NO at the interface. It will be noted that all
of the NO reactions also involve N atoms. In fact the N atoms are pro-
duced and destroyed by these reactions almost exclusively. More than 99%
of the N atoms are produced by
N2 + 0 = N + NO
while about 60% of the N is destroyed by
N + 02 ——^ NO + 0
and 40% is destroyed by
HO + N "' »• H + NO
Using only the reactions for which rate data are available, at one point
within the branching zone as much as 25% of N atoms are destroyed by
reactions with hydrocarbon fragments to produce CRN and CN. These species
do not react directly to produce NO, but produce either N atoms or N2
molecules through subsequent reactions. However, the main N atom reactions
for production and destruction in these calculations were those involving NO.
-------
- 72 -
3.7 Comparison of Theoretical Calculations
with Experimental Results
The comparisons between theoretical calculations discussed
in Section 3.4 and the experimental results discussed in Section 3 are
encouraging. As expected, because of the artificial constraints placed
on the temperature in the theoretical calculations, the NO levels for
the plug-flow simulation of the flat flame are somewhat overpredicted.
However, at this point it is not certain whether the overpredictions are
more a result of the high temperatures early in the calculations or because
the rate constants used for NO formation reactions may be too high. The
latter interpretation is favored by the fact that the theoretical
overprediction continues to be observed in the post-flame zone. The
temperature profile used for reaction screening may have a less dramatic
effect than initially expected because the major portion of the reactions
involving heat release are completed before substantial NO has formed.
The theoretical-experimental comparisons under stoichiometric or
excess air conditions are the closest. Under fuel rich conditions, NO con-
centrations are underpredicted. These underpredictions could be caused by
the lack of reactions between hydrocarbon fragments and molecular nitrogen,
or the high rate of disappearance of active oxygen-containing intermediates.
The high rate of disappearance of both hydrocarbon species and oxygen-
containing intermediates would be caused by the use of isothermal conditions
in the combustion zone where experimentally the temperature is increasing
from room temperature to flame temperature.
Under stoichiometric air conditions, plug flow calculations
overpredict experimentally measured concentrations of NO in the post-flame
zone by about a factor of two. Theoretical calculations indicate 700 ppm
at 150 msec while only 350 ppm are measured; 1200 ppm are calculated at
300 msec while 600 ppm are measured. In the post-flame zone the rate of
NO formation between 200 and 500 msec is measured to be 1500 ppm/sec while
the calculated rate is about 3000 ppm/sec. The comparison between theory
and experiment, on a preliminary basis, indicate that a reduction of the
NO formation rate constants by a factor of two could bring theory in line
with experiment. However, additional theoretical analysis is required
before such reductions may be considered permissible.
A similar picture is obtained at 120 percent stoichiometric air.
Experimental measurements at 40 msec indicate concentrations in the 40-50
ppm range while calculations indicate about 90 ppm. Measurements indicate
concentrations of about 170 ppm at 500 msec while calculations indicate
about 480 ppm. The rate of NO formation is measured to be about 400 ppm/sec
between 300 and 500 msec while calculations indicate about 900 ppm/sec.
Under fuel rich conditions experimental measurements indicate
a very rapid formation of NO in the flame zone and essentially a zero
net formation rate thereafter. Calculations indicate a similar behavior,
with essentially all the NO formed in the relaxation zone, and little or
none in the branching or post-flame zones. However, the levels of NO formed
experimentally are much higher than those calculated in the screening
-------
- 73 -
studies. This is likely to be caused, in part, by the overly rapid calculated
depletion of methane and oxygen in the branching zone because of the high
temperature assumed. However, another possible explanation is because of
the lack of reactions involving hydrocarbons in NO formation in these
calculations. This possibility will be pursued further in Phase II of this
study. Measured NO formation in the flame zone is on the order of 60 ppm,
while the calculated value is on the order of 6 ppm.
-------
- 74 -
4. CONCLUSIONS
1. The survey of 142 reactions appearing in the literature for 25
species in the Clfy/air system resulted in recommended rates for
127 reactions to be used for reaction screening. Since 322 reactions
could be postulated for these 25 species, additional attention is
needed to determine if any of these additional reactions could be
important.
2. Reaction screening using only these 127 reactions or smaller subsets
indicates that between 80% and 120% stoichiometric air for methane
combustion:
- NO formation rate is controlled primarily by oxygen atom
concentration.
- Oxygen atom concentrations are substantially above equilibrium
in the branching zone before significant NO has formed.
- Oxygen atom concentrations have peaked before significant NO
has formed.
- Oxygen atom concentration controls formation of N atoms which
participate in NO reactions. N atom destruction results in NO
formation by N + OH and N + 0 .
- While reaction screening for the 127 reaction set has not uncovered
reactions other than N£ + 0 which might produce N atoms in the
branching zone, it is still worth pursuing the search for other
N-forming reactions that might occur in the branching zone. Such
reactions might produce more substantial overshoots of N atom
concentrations which could provide more rapid NO formation early
in the flame zone.
3. Reaction screening from the 127 reaction set has resulted in narrowing
down the number of reactions that need be considered for NO formation
under the conditions studied. Caution must be exercised not to
extrapolate these conclusions outside the range studied without
rechecking the calculations under the new conditions. Further
screening will be required when other reactions from the 322 reaction
set (25 species) are tested for importance.
4. As additional reactions are studied or estimated the minimum set
of reactions may change:
- The new reaction may be added to the set by itself if it
provides a new path between two important intermediates
(can affect any zone).
- The new reaction may result in the rapid destruction of an
intermediate as it is formed, preventing overshoot from occurring.
In that case the minimum set would be reduced by eliminating
reactions of that intermediate (branching zone).
The new reaction may result in rapid production of an intermediate
in the branching zone and corresponding overshoot in its concentration.
Thus, a new intermediate may play a role in NO formation - possibly
HN, CHN or CN (branching zone).
-------
75 -
5. Premixed flat and focused flames as well as laminar and turbulent
diffusion flames have been studied in the multiburner under both
cold-and hot-wall conditons:
- The burners appeared to behave as idealized examples of the
types they represented. Concentrations and temperature profiles
have been measured in these flames.
- Under certain conditions the curvature in NO formation
(non-linear extrapolation back to the origin) could be measured.
6. Comparisons of isothermal plug-flow calculations at the adiabatic
flame temperature with experimental data indicate:
- NO formation is overpredicted by a factor of two in both the
flame zone and post-flame zone for fuel-lean combustion and in
the post-flame zone for stoichiometric air conditions. NO formation
is underpredicted for fuel rich combustion by an order of magnitude.
+ Overprediction in the post flame zone under fuel lean conditions,
where intermediates other than NO have reached essentially
equilibrium concentrations, indicate that the controlling NO
formation rate, N£ + 0 = NO + N, should be reduced by a factor
of two. (This reaction has a high activation energy and
therefore temperature uncertainty must be taken into account.)
+ The possible interaction between hydrocarbon species and
nitrog'e'rious species should be investigated to determine their
possible role under fuel-rich conditions.
- Extremely high concentration gradients are predicted by the plug
flow calculations which do not include diffusion effects. In a
real system, such concentration gradients would be smoothed
somewhat by diffusion.
- For the purpose of the present study that concerns itself with the
chemistry of combustion and pollutant formation reactions, it would
be desirable to eliminate diffusion effects from the experimental
system. To achieve that goal, the development of an adiabatic
sti-rr«d-combust-or is indicated. Development of this combustor type
which incorporates the principles of mixing in the jet-stirred
combustor and adiabatic operation in the multiburner will be
discussed in the report on Phase II of this contract study.
-------
A-l
APPENDIX A
1.
2*
3.
•
4.
5.
6.
8.
9l
10.
.
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A *
A 11.
A .
A .
A 12.
A 13.
A .
A 14.
A 15.
A 16.
A .
A .
A .
A 17.
A .
A
A
A
A
A
A
A 19.
A .
A 20.
A
A
A
A
A
A
A
A 21.
A .
A .
A .
A .
A 23.
A 22.
.
18.
.
.
.
.
CH
CH
CH
CH
CH
CH
CH
CH
CH
CH
CH
CH
CH
CH
CH
CH
CH
CH
CH
CH
CH
CH
CH
CH
CH
CH
CH
CH
CH
CH
CH
CH
CH
CH
CH
CH
CH
CH
CH
CH
CH
CH
CH
CH
CH
CH
CH
CH
CH
CH
CH
CH
CH
4-
4-
4-
4-
4-
4-
+
4-
•f
4-
+
4-
+
4-
•f
+
+
+
+
+
•f
+
•«•
4-
+
-t-
+
4-
-f
4-
4-
+
4-
4-
4-
4-
4-
4-
4-
4-
4-
4-
4-
4-
+
+
4-
4-
4-
4-
4-
4-
4-
CHN
CHO
CH20
CH20
CH3
CH30
CH30
CH30
CH4
C02
H
HN
HN
HN
HNO
HNO
HNO
HNO
HNO
HO
HO
HO
HO
H02
H02
H02
H02
H02
H2
H2
H20
H20
H20
H20
H20
N
N
NO
NO
NO
NO
N02
N02
N02
N02
NO?
N2
N2
N20
N20
N20
0
0
4- M
4- M
4- M
4- M
4- M
4- M
REACTIONS FOR 2!
CH2
CH2
CHO
CH3
CH2
CHO
CH2
CH4
CH2
CHO
CH2
CHN
CH2
CN
CHN
CHO
CH2
CH20
CN
CHO
CH2
CH20
CO
CHO
CH2
CH20
CO
C02
CH2
CH3
CHO
CH2
CH20
CH3
CH30
CHN
CN
CHN
CHO
CN
CO
CHN
CHO
CN
CO
C02
CHN
CN
CHN
CHO
CN
CHO
CO
4-
+
4-
4-
•f
•f
4-
+
+
"f
4-
4-
•f
4-
4-
•f
4-
4-
+
•f
4-
4-
4-
-f
4-
4-
4-
4-
+
4-
4-
4-
4-
•f
4-
-f
4-
4-
4-
•f
4-
4-
4-
4-
+
4-
+
CN
CO
CH2
CO
CH2
CH3
CH20
CO
CH3
CO
H
N
H2
HO
HN
NO
N
H20
H
0
H2
HO
02
0
H20
H2
H
H2
HO
H
0
H
0
N
HO
HN
02
NO
H02
HNO
HN
N
HN
NO
N2
HNO
H
+ M
+ M
•f M
-t- M
+ M
4- M
H-0 EXCHANGE
0 TRANSFER
H TRANSFER
H-0 EXCHANGE
N TRANSFER
N TRANSFER
0 TRANSFER
(HO) TRANSFER
H-N EXCHANGE
(HO) ADDITION
(HO) TRANSFER
H-0 EXCHANGE
C-H EXCHANGE
(H2) ADDITION
0 TRANSFER
(HO) TRANSFER
(H2) TRANSFER
(H20) ADDITION
N ADDITION
N TRANSFER
H-N EXCHANGE
H-0 EXCHANGE
N TRANSFER
0 TRANSFER
H-N EXCHANGE
H-0 EXCHANGE
C-N EXCHANGE
H-N EXCHANGE
N TRANSFER
0 TRANSFER
H-N EXCHANGE
X
X
X
X
X
X
X
X
X
1
2
3
4
5
6
7
8
9
X 10
X 11
X 12
X 13
X 14
X 15
X 16
X 17
X 18
X 19
X 20
X 21
X 22
X 23
X 24
X 25
X 26
X 27
X 28
X 29
X 30
X 31
X 32
X 33
X 34
X 35
X 36
X 37
X 38
X 39
X 40
X 41
X 42
X 43
X 44
X 45
X 46
X 47
X 48
X 49
X 50
X 51
X 52
X 53
-------
A-2
MASTER LIST OF REACTIONS FOR 25 ALLOWED SPECIES
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
24.
25.
.
26.
27.
28.
.
29.
30.
.
31.
.
.
,
.
,
.
32.
.
.
*
.
.
.
.
»
.
33.
.
34.
.
.
35.
.
.
.
.
»
36.
,
37.
38.
.
.
39.
40.
.
.
.
41.
42.
43.
44.
CH
CH
CH
CHN
CHN
CHN
CHN
CHN
CHN
CHN
CHN
CHN
CHN
CHN
CHN
CHN
CHN
CHN
CHN
CHN
CHN
CHN
CHN
CHN
CHN
CHN
CHN
CHN
CHN
CHN
CHN
CHN
CHN
CHN
CHN
CHN
CHN
CHN
CHO
CHO
CHO
CHO
CHO
CHO
CHO
CHO
CHO
CHO
CHO
CHO
CHO
CHO
CHO
+ 02
+ 02
-t- 02
+ CHO
+ CH2
+ CH20
+ CH3
+ CO
+ H
+ H
+ HN
+ HNO
+ HNO
+ HO
+ HO
+ HO
+ HO
+ H02
+ H02
+ H02
+ H2
+ H2
+ H20
-»• H20
+ H20
+ H20
+ N
+ NO
+ NO
"»- N02
+ 0
4- 0
+ 0
+ 02
+ 02
+ 02
4- 02
-t- CHO
+ CHO
+ CH2
+ CH20
+ CH20
+ CH3
+ CH3
+ CH30
+ CH30
+ CH4
+ CH4
+ H
+ H
+ H
+ M
M
=
=
3
=
=
=
S
E
=
=
=
=
=
=
=
3
=
S
=
=
=
=
=
£
=
=
=
s
=
=
=
=
s
=
=
=
=
=
=
=
=
=
s
=
=
=
=
s
=
s
=
=
=
CHO
CO
C02
CN
CH20
CH3
CH30
CH4
CHO
CH2
CN
CH2
CH2
CH20
CHO
CH2
CH20
CN
CHO
CH2
CH20
CH2
CH3
CH2
CH20
CH3
CH30
CN
CHO
CN
CHO
CHO
CN
CO
CHO
CN
CO
C02
CO
CH2
CH20
CH3
CH3
CH30
CH2
CH4
CH20
CH4
CH2
CH20
CH2
CH20
CO
4- 0
+ HO
+ H
+ H
+ CN
+ CN
+ CN
+ CN
+ CN
4- N
+ H2
-t-. ,\2
+ N20
+ N2
+ HN
+ NO
4- N
4- H20
+ HNO
4- N02
+ NO
+ HN
4- N
+ HNO
4- HN
+ NO
4- N
+ HN
+ N2
+ HNO
+ N20
•f N
+ HO
+ HN
+ NO
4- H02
+ HNO
-t- HN
+ H
+ C02
+ CO
+ CO
+ C02
+ CO
4- CH20
+ CO
+ CH20
4- C02
+ CH30
+ CH3
-t- 0
-t- H2
4- M
+ M
C TRANSFER
H TRANSFER
H-N EXCHANGE
H-N EXCHANGE
H-N EXCHANGE
N-(HO) EXCHANGE
N-0 EXCHANGE
H-N EXCHANGE
N-(HO) EXCHANGE
N-0 EXCHANGE
H-N EXCHANGE
N-(HO) EXCHANGE
H-N EXCHANGE
N-(H2) EXCHANGE
H-N EXCHANGE
N-(HO) EXCHANGE
N-(H2) EXCHANGE
CH TRANSFER
N-0 EXCHANGE
N-0 EXCHANGE
N-0 EXCHANGE
O-(HN) EXCHANGE
N-0 EXCHANGE
H TRANSFER
O-(HN) EXCHANGE
C TRANSFER
H-0 EXCHANGE
H-0 EXCHANGE
H TRANSFER
H TRANSFER
H-0 EXCHANGE
H-0 EXCHANGE
4- M
X 54
X 55
X 56
X 57
X 58
X 59
X 60
X 61
X 62
X 63
X 64
X 65
X 66
X 67
X 68
X 69
X 70
X 71
X 72
X 73
X 74
X 75
X 76
X 77
X 78
X 79
X 80
X 81
X 82
X 83
X 84
X 85
X 86
X 87
X 88
X 89
X 90
X 91
X 92
X 93
X 94
X 95
X 96
X 97
X 98
X 99
X100
X101
X102
X103
X104
X105
X106
-------
A-3
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
*
*
*
*
45.
•
46*
47.
.
.
.
.
48.
.
.
*
49.
*
.
.
50.
.
51.
•
.
*
52.
.
.
.
.
.
53.
.
.
.
.
54.
.
.
55.
.
.
.
.
56.
.
.
57.
.
.
.
*
CHO
CHO
CHO
CHO
CHO
CHO
CHO
CHO
CHO
CHO
CHO
CHO
CHO
CHO
CHO
CHO
CHO
CHO
CHO
CHO
CHO
CHO
CHO
CHO
CHO
CHO
CHO
CHO
CHO
CHO
CH2
CH2
CH2
CH2
CH2
CH2
CH2
CH2
CH2
CH2
CH2
CH2
CH2
CH2
CH2
CH2
CH2
CH2
CH2
CH2
CH2
CH2
CH2
+
+
+
+
•f
+
•+•
+
•*•
+
+
+
+
•f
+
+
•f
+
+
+
•f
+
+
+
+
•f
•f
+
+
4-
+
•f
•f
+
+
•f
•f
•f
+
+
+
•f
+
•f
+
•f
•f
+
•f
+
+
•f
•f
HN
HN
HN
HNO
HNO
HO
HO
HO
HO
H02
H02
H2
H2
H2
H2
H20
H20
H20
H20
N
N
NO
NO
NO
N02
N2
0
0
02
02
CH20
CH30
CH4
C02
H
HN
HNO
HNO
HNO
HO
HO
HO
H02
H02
H02
H2
H2
H20
H20
H20
H20
NO
N
MASTER LIST OF REACTIONS FOR 25 ALLOWED SPECIES
H-0 EXCHANGE
H TRANSFER
N-(HO) EXCHANGE
H-O EXCHANGE
H-0 EXCHANGE
H-O EXCHANGE
H TRANSFER
H-0 EXCHANGE
H-0 EXCHANGE
+ M
+ M
M
M
=
=
s
=
=
s:
a
=
E
S
=
r
=
=
s
=
s
ss
=
=
=
=
=
=
s
=
sz
=
=
=
s
=
=
=
X
s
=
s
=
=
=
s
=
=
S5
s
s
r
=
s
c
3
=
CH?
CH20
CN
CH2
CH20
CH2
CH20
CO
C02
CH20
C02
CH2
CH20
CHS
CH30
CH2
CH20
CH3
CH30
CN
CO
CN
CO
C02
C02
CN
CO
C02
CO
C02
CH4
CH20
CH3
CH20
CHS
CHS
CH20
CHS
CH30
CH20
CHS
CH30
CH20
CHS
CH30
CH3
CH4
CH20
CH3
CH30
CH4
CH20
CN
+ NO
+ N
+ H20
+ N02
+ NO
+ 02
+ 0
+ H20
+ H2
+ 02
+ H20
+ HO
+ H
+ 0
+ H02
+ HO
+ 02
+ 0
+ HO
+ HN
+ H02
+ HNO
+ HN
+ HNO
+ HNO
+ HO
+ H
+ H02
+ HO
+ CO
+ CHS
+ CHS
+ CO
+ N
+ HN
+ NO
+ N
+ H
+ 0
+ HO
+ 02
+ 0
+ H
+ H2
+ HO
+ H
+ 0
+ N
+ H2
+ M
0-(H2) EXCHANGE
(H2) ADDITION
H-0 EXCHANGE
0-(H2) EXCHANGE
(H2) TRANSFER
N-(HO) EXCHANGE
N-(HO) EXCHANGE
H-0 EXCHANGE
H-0 EXCHANGE
N-(HO) EXCHANGE
H-0 EXCHANGE
H TRANSFER
H-0 EXCHANGE
|