37 reactantNames_(dict.
lookup(
"reactants")),
60 r_[i] = st[reactantNames_[i]];
67 inline Foam::scalar Foam::LangmuirHinshelwoodReactionRate::operator()
74 const scalar c0m =
pow(c[r_[0]], m_[0]);
75 const scalar c1m =
pow(c[r_[1]], m_[1]);
77 const scalar TaByT0 = Ta_[0]/
T;
78 const scalar TaByT1 = Ta_[1]/
T;
79 const scalar TaByT2 = Ta_[2]/
T;
81 const scalar k0 = A_[0]*
pow(T, beta_[0])*
exp(-TaByT0);
82 const scalar k1 = A_[1]*
pow(T, beta_[1])*
exp(-TaByT1);
83 const scalar k2 = A_[2]*
pow(T, beta_[2])*
exp(-TaByT2);
85 return k0/
pow(a_ + k1*c0m + k2*c1m, m_[2]);
96 const scalar c0m =
pow(c[r_[0]], m_[0]);
97 const scalar c1m =
pow(c[r_[1]], m_[1]);
99 const scalar TaByT0 = Ta_[0]/
T;
100 const scalar TaByT1 = Ta_[1]/
T;
101 const scalar TaByT2 = Ta_[2]/
T;
103 const scalar k0 = A_[0]*
pow(T, beta_[0])*
exp(-TaByT0);
104 const scalar k1 = A_[1]*
pow(T, beta_[1])*
exp(-TaByT1);
105 const scalar k2 = A_[2]*
pow(T, beta_[2])*
exp(-TaByT2);
109 (beta_[0] + TaByT0)*k0
110 - m_[2]*k0*((beta_[1] + TaByT1)*k1*c0m + (beta_[2] + TaByT2)*k2*c1m)
111 /(a_ + k1*c0m + k2*c1m)
112 )/(
pow(a_ + k1*c0m + k2*c1m, m_[2])*
T);
119 return NullObjectRef<List<Tuple2<label, scalar>>>();
#define forAll(list, i)
Loop across all elements in list.
Langmuir-Hinshelwood reaction rate for gaseous reactions on surfaces.
scalar ddT(const scalar p, const scalar T, const scalarField &c) const
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.
dimensionedScalar exp(const dimensionedScalar &ds)
LangmuirHinshelwoodReactionRate(const speciesTable &species, const dictionary &dict)
Construct from dictionary.
An Ostream is an abstract base class for all output systems (streams, files, token lists...
scalar dcidT(const scalar p, const scalar T, const scalarField &c) const
Temperature derivative of the pressure dependent term.
const List< Tuple2< label, scalar > > & beta() const
Third-body efficiencies (beta = 1-alpha)
void writeEntry(Ostream &os, const HashTable< T, Key, Hash > &ht)
void write(Ostream &os) const
Write to stream.
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
T lookupOrDefault(const word &, const T &, bool recursive=false, bool patternMatch=true) const
Find and return a T,.
A wordList with hashed indices for faster lookup by name.
const dimensionedScalar c
Speed of light in a vacuum.
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.