46 for (
int i=0; i<n_; i++)
60 co_(species.find(
"CO")),
61 c3h6_(species.find(
"C3H6")),
62 no_(species.find(
"NO"))
69 A_[i] = coeffs[i].
first();
70 Ta_[i] = coeffs[i].second();
84 return A_[0]*
exp(-Ta_[0]/
T)/
87 *
sqr(1 + A_[1]*
exp(-Ta_[1]/
T)*c[co_] + A_[2]*
exp(-Ta_[2]/
T)*c[c3h6_])
88 *(1 + A_[3]*
exp(-Ta_[3]/
T)*
sqr(c[co_])*
sqr(c[c3h6_]))
89 *(1 + A_[4]*
exp(-Ta_[4]/
T)*
pow(c[no_], 0.7))
100 coeffs[i].
first() = A_[i];
101 coeffs[i].second() = Ta_[i];
static const Foam::dimensionedScalar A("", Foam::dimPressure, 611.21)
A 1D vector of objects of type <T> with a fixed length <N>.
T & first() noexcept
The first element of the list, position [0].
Power series reaction rate.
An Ostream is an abstract base class for all output systems (streams, files, token lists,...
virtual Ostream & writeKeyword(const keyType &kw)
Write the keyword followed by an appropriate indentation.
virtual bool write(const token &tok)=0
Write token to stream or otherwise handle it.
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
ITstream & lookup(const word &keyword, enum keyType::option matchOpt=keyType::REGEX) const
friend Ostream & operator(Ostream &, const faMatrix< Type > &)
virtual bool write()
Write the output fields.
A wordList with hashed named lookup, which can be faster in some situations than using the normal lis...
OBJstream os(runTime.globalPath()/outputName)
dimensionedScalar exp(const dimensionedScalar &ds)
dimensionedSymmTensor sqr(const dimensionedVector &dv)
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
constexpr char nl
The newline '\n' character (0x0a)
#define forAll(list, i)
Loop across all elements in list.