Go to the documentation of this file.00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028 namespace Foam
00029 {
00030
00031
00032
00033
00034 inline JanevReactionRate::JanevReactionRate
00035 (
00036 const scalar A,
00037 const scalar beta,
00038 const scalar Ta,
00039 const scalar b[]
00040 )
00041 :
00042 A_(A),
00043 beta_(beta),
00044 Ta_(Ta)
00045 {
00046 for (int n=0; n<nb_; n++)
00047 {
00048 b_[n] = b[n];
00049 }
00050 }
00051
00052
00053
00054 inline JanevReactionRate::JanevReactionRate
00055 (
00056 const speciesTable&,
00057 Istream& is
00058 )
00059 :
00060 A_(readScalar(is.readBegin("JanevReactionRate(Istream&)"))),
00061 beta_(readScalar(is)),
00062 Ta_(readScalar(is))
00063 {
00064 for (int n=0; n<nb_; n++)
00065 {
00066 is >> b_[n];
00067 }
00068
00069 is.readEnd("JanevReactionRate(Istream&)");
00070 }
00071
00072
00073
00074
00075 inline scalar JanevReactionRate::operator()
00076 (
00077 const scalar T,
00078 const scalar,
00079 const scalarField&
00080 ) const
00081 {
00082 scalar lta = A_;
00083
00084 if (mag(beta_) > VSMALL)
00085 {
00086 lta *= pow(T, beta_);
00087 }
00088
00089 scalar expArg = 0.0;
00090
00091 if (mag(Ta_) > VSMALL)
00092 {
00093 expArg -= Ta_/T;
00094 }
00095
00096 scalar lnT = log(T);
00097
00098 for (int n=0; n<nb_; n++)
00099 {
00100 expArg += b_[n]*pow(lnT, n);
00101 }
00102
00103 lta *= exp(expArg);
00104
00105 return lta;
00106 }
00107
00108
00109 inline Ostream& operator<<(Ostream& os, const JanevReactionRate& jrr)
00110 {
00111 os << token::BEGIN_LIST
00112 << jrr.A_ << token::SPACE << jrr.beta_ << token::SPACE << jrr.Ta_;
00113
00114 for (int n=0; n<JanevReactionRate::nb_; n++)
00115 {
00116 os << token::SPACE << jrr.b_[n];
00117 }
00118
00119 os << token::END_LIST;
00120
00121 return os;
00122 }
00123
00124
00125
00126
00127 }
00128
00129