research/cantera
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

Removed unimplemented TST and LandauTeller kinetic rates

Browse source
This commit is contained in:
Ray Speth 2012-07-20 18:26:50 +00:00
parent 521b2957c9
commit d7463cc226
1 changed files with 0 additions and 73 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

73
include/cantera/kinetics/RxnRates.h
View file

@ -246,52 +246,6 @@ protected:
};
#ifdef INCL_TST
class TST
{
public:
static int type() {
return TSTRATE;
}
TST() {}
TST(const vector_fp& c) {
m_b.resize(10);
copy(c.begin(), c.begin() + 10, m_b.begin());
m_k = int(c[10]);
}
void update_C(const vector_fp& c) {
doublereal ck = c[m_k];
delta_s0 = m_b[0] + m_b[1]*ck + m_b[2]*ck*ck;
delta_e0 = m_b[5] + m_b[6]*ck + m_b[7]*ck*ck;
}
doublereal update(doublereal logT, doublereal recipT) const {
doublereal delta_s = delta_s0*(1.0 + m_b[3]*logT + m_b[4]*recipT);
doublereal delta_E = delta_e0*(1.0 + m_b[8]*logT + m_b[9]*recipT);
return logBoltz_Planck + logT + delta_s - delta_E*recipT;
}
doublereal updateRC(doublereal logT, doublereal recipT) const {
doublereal lres = update(logT, recipT);
return exp(lres);
}
void writeUpdateRHS(std::ostream& s) const {}
protected:
doublereal delta_s0, delta_e0;
int m_k;
vector_fp m_b;
};
#endif
//! Arrhenius reaction rate type depends only on temperature
/**
* A reaction rate coefficient of the following form.
@ -675,33 +629,6 @@ protected:
vector_fp dotProd_; //!< dot product of chebCoeffs with the reduced pressure polynomial
};
// class LandauTeller {
// public:
// static int type(){ return LANDAUTELLER; }
// LandauTeller(){}
// LandauTeller( const vector_fp& c ) : m_c(c) { m_c[0] = log(c[0]); }
// doublereal update(doublereal logT, doublereal recipT) const {
// return m_c[0] + m_c[1]*tt[1] - m_c[2]*tt[2]
// + m_c[3]*tt[3] + m_c[4]*tt[4];
// }
// //void writeUpdateRHS(ostream& s) const {
// // s << exp(m_logA);
// // s << " * exp(";
// // if (m_b != 0.0) s << m_b << " * tlog";
// // if (m_E != 0.0) s << " - " << m_E << " * rt";
// // if (m_E != 0.0) s << " - " << m_E << " * rt";
// // s << ");" << endl;
// // }
// //}
// protected:
// doublereal m_logA, m_b, m_E;
// };
//}
}
#endif