51 A_(
dict.lookup<scalar>(
"A")),
52 beta_(
dict.lookup<scalar>(
"beta")),
53 Ta_(
dict.lookup<scalar>(
"Ta")),
54 B_(
dict.lookup<scalar>(
"B")),
55 C_(
dict.lookup<scalar>(
"C"))
69 inline Foam::scalar Foam::LandauTellerReactionRate::operator()
79 if (
mag(beta_) > vSmall)
86 if (
mag(Ta_) > vSmall)
98 expArg += C_/
pow(
T, 2.0/3.0);
101 if (
mag(expArg) > vSmall)
120 if (
mag(beta_) > vSmall)
122 lta *=
pow(
T, beta_);
128 if (
mag(Ta_) > vSmall)
135 if (
mag(B_) > vSmall)
137 scalar BT = B_/
cbrt(
T);
142 if (
mag(C_) > vSmall)
144 scalar CT = C_/
pow(
T, 2.0/3.0);
149 if (
mag(expArg) > vSmall)
154 return lta*(beta_+deriv)/
T;
static const Foam::dimensionedScalar A("A", Foam::dimPressure, 611.21)
static const Foam::dimensionedScalar B("B", Foam::dimless, 18.678)
Graphite solid properties.
Landau-Teller reaction rate.
void preEvaluate() const
Pre-evaluation hook.
void postEvaluate() const
Post-evaluation hook.
void write(Ostream &os) const
Write to stream.
void ddc(const scalar p, const scalar T, const scalarField &c, const label li, scalarField &ddc) const
The derivative of the rate w.r.t. concentration.
LandauTellerReactionRate(const scalar A, const scalar beta, const scalar Ta, const scalar B, const scalar C)
Construct from components.
scalar ddT(const scalar p, const scalar T, const scalarField &c, const label li) const
The derivative of the rate w.r.t. temperature.
bool hasDdc() const
Is the rate a function of concentration?
An Ostream is an abstract base class for all output systems (streams, files, token lists,...
virtual Ostream & write(const char)=0
Write character.
A list of keyword definitions, which are a keyword followed by any number of values (e....
A wordList with hashed indices for faster lookup by name.
const dimensionedScalar c
Speed of light in a vacuum.
dimensionedScalar exp(const dimensionedScalar &ds)
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
void writeEntry(Ostream &os, const HashTable< T, Key, Hash > &ht)
dimensioned< scalar > mag(const dimensioned< Type > &)
dimensionedScalar cbrt(const dimensionedScalar &ds)
void T(FieldField< Field, Type > &f1, const FieldField< Field, Type > &f2)