reactionTableRhs

Computes the right hand side of the reaction rate equation and the reaction energy change rate. This can then be used in a time integration scheme. Any number of reactions can be added for a given set of species.

\[\notag \begin{align} \frac{d n_{i}}{d t}=s_{i} \end{align}\]
\[\notag \begin{align} \frac{d e}{d t}=s_{re} \end{align}\]

Parameters

species (string vector)
List of species to include in the reactions.
fileName
Input file containing the reaction rate constants data (REACTIONS), specific heats data (CP) and energy of formation (EOF). Refer to SpeciesDataFile for the input data format.
maxRate
user specified value to limit the maximum rate of reactions.
outputEnergyRate
option to specify whether to compute reaction energy using specific heats. NOTE: if the option is false, then in vector requires only first and second variable as inputs. out vector requires only first variable as input.

Parent Updater Data

in (string vector, required)

1st Variable

number densities of species \(m^{-3}\)

2nd Variable

average temperature of the fluid \(K\)

3rd Variable

specific heat at constant pressure of the species \(\frac{J}{kg K}\)

out (string vector, required)

1st Variable

time rate of change of species density \(\frac{1}{m^{3} s}\)

2nd Variable

time rate of change of energy \(\frac{J}{s}\)

Example

<Updater sourceUpdater>
  kind = equation1d
  onGrid = domain

  in = [speciesDensity, temperature, specificHeat]
  out = [speciesDensitySource, reactionEnergySource]

  equations = [reactionSource]

   <Equation reactionSource>
     kind = reactionTableRhs
     species = [N2, N, O2, O, NO, NO_p1, e]
     fileName = air7Species.txt
     maxRate = 1.0e28
     outputEnergyRate = 1
  </Equation>

</Updater>