Consider elementary reversible (or irreversible) reactions involving K chemical species that can be represented in the general form
(3–1) |
The stoichiometric coefficients are integer numbers and
is the chemical symbol for the k th species. The
superscript
indicates forward stoichiometric coefficients, while
indicates reverse stoichiometric coefficients. Normally, an elementary
reaction involves only three or four species; hence the
matrix is quite sparse for a large set of reactions. For non-elementary
reactions, Equation 3–1
also represents the reaction
expression, but the stoichiometric coefficients may be non-integers.
Note: Global reactions are sometimes stated with non-integer stoichiometric coefficients. Ansys Chemkin can accommodate non-integer stoichiometric coefficients.
The production rate of the k th species can be written as a summation of
the rate-of-progress variables for all reactions involving the k th
species
(3–2) |
where
(3–3) |
The rate of progress variable for the i th reaction is given by the difference of
the forward and reverse rates as
(3–4) |
where is the molar concentration of the k th species and
and
are the forward and reverse rate constants of the i th
reaction. As indicated in Equation 3–4
, the rate-of-progress of a
reaction is evaluated, by default, using the concentration of each reactant or product species
raised to the power of its stoichiometric coefficient. Thus, the rate-of-progress of a
reaction that includes species
with a coefficient of 2 will be second-order with respect to the
concentration of
. Equation 3–4
is always valid when mass-action
kinetics are obeyed, and when the mechanism is written in terms of elementary reactions. As it
is often useful to work with reduced chemistry, Gas-phase
Kinetics includes an option allowing the user to define an arbitrary
reaction order for a species in place of the coefficients used in Equation 3–4
. This option is described further
below.
The forward rate constants for the reactions are generally assumed to have the following Arrhenius temperature
dependence:
(3–5) |
where the pre-exponential factor , the temperature exponent
, and the activation energy
are specified. These three parameters are required input to the
Gas-phase Kinetics package for each reaction.
Note: Two gas constants, and
are used throughout this chapter and within Ansys Chemkin programs.
is used only in conjunction with the activation energy
and has compatible units. The reason for the duality is that many users
would rather use units of cal/mole for the activation energies even though other energy
units are used elsewhere.
In Equation 3–5
through Equation 3–10
, refers to the gas temperature, unless auxiliary reaction information is
provided to indicate that the reaction depends on a temperature associated with a particular
species. Such information would be specified using the auxiliary keyword,
TDEP
, which is described further in Table 3.7: Alphabetical Listing of Gas-phase Reaction Auxiliary Keywords
of the
Chemkin Input Manual
Input Manual. In the
case where the TDEP
keyword is included for reaction ,
represents the temperature of the species whose name follows the
TDEP
keyword.
In thermal systems, the reverse rate constants are related to the forward rate constants through the equilibrium constants
by
(3–6) |
Although is given in concentration units, the equilibrium constants are more easily
determined from the thermodynamic properties in pressure units; they are related by
(3–7) |
The equilibrium constants are obtained with the relationship
(3–8) |
The refers to the change that occurs in passing completely from reactants to
products in the i th reaction; specifically,
(3–9) |
(3–10) |
For reactions involving electrons, the use of equilibrium constants to determine reverse rates
is usually not appropriate. In some cases, detailed balancing on electron-driven reactions can
be applied using the Saha equation (see, for example, Mitchner and Kruger [14]) that relates the ionization and electron-third-body recombination reactions to the
species partition functions. While such relations can be used to calculate explicitly reverse
rates from forward rates, they are not part of the built-in features of
Gas-phase Kinetics. To avoid erroneous results, it is therefore required
that all reactions involving electron species must either be specified as forward reactions
only, or must include the reverse rate parameters explicitly stated using auxiliary keywords.
The specification of reverse-rate parameters is described in more detail in the
REV
entry in Table 3.7: Alphabetical Listing of Gas-phase Reaction Auxiliary Keywords
of the
Chemkin Input Manual
Input
Manual).