61 inline Foam::scalar Foam::LandauTellerReactionRate::operator()
70 if (
mag(beta_) > vSmall)
77 if (
mag(Ta_) > vSmall)
89 expArg += C_/
pow(T, 2.0/3.0);
92 if (
mag(expArg) > vSmall)
110 if (
mag(beta_) > vSmall)
112 lta *=
pow(T, beta_);
118 if (
mag(Ta_) > vSmall)
125 if (
mag(B_) > vSmall)
127 scalar BT = B_/
cbrt(T);
132 if (
mag(C_) > vSmall)
134 scalar CT = C_/
pow(T, 2.0/3.0);
139 if (
mag(expArg) > vSmall)
144 return lta*(beta_+deriv)/T;
151 return NullObjectRef<List<Tuple2<label, scalar>>>();
LandauTellerReactionRate(const scalar A, const scalar beta, const scalar Ta, const scalar B, const scalar C)
Construct from components.
A list of keyword definitions, which are a keyword followed by any number of values (e...
A 1D array of objects of type <T>, where the size of the vector is known and used for subscript bound...
void dcidc(const scalar p, const scalar T, const scalarField &c, scalarField &dcidc) const
Species concentration derivative of the pressure dependent term.
scalar ddT(const scalar p, const scalar T, const scalarField &c) const
const List< Tuple2< label, scalar > > & beta() const
Third-body efficiencies (beta = 1-alpha)
dimensionedScalar exp(const dimensionedScalar &ds)
dimensionedScalar cbrt(const dimensionedScalar &ds)
An Ostream is an abstract base class for all output systems (streams, files, token lists...
void writeEntry(Ostream &os, const HashTable< T, Key, Hash > &ht)
scalar dcidT(const scalar p, const scalar T, const scalarField &c) const
Temperature derivative of the pressure dependent term.
void write(Ostream &os) const
Write to stream.
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
Landau-Teller reaction rate.
A wordList with hashed indices for faster lookup by name.
dimensioned< scalar > mag(const dimensioned< Type > &)
virtual Ostream & write(const token &)=0
Write next token to stream.
ITstream & lookup(const word &, bool recursive=false, bool patternMatch=true) const
Find and return an entry data stream.