thermophysicalModels/specie/thermo/thermo/thermoI.H: Reverting change to filter K based on the nMoles

This commit is contained in:
Henry
2012-11-01 13:17:02 +00:00
parent bfa8ffc25e
commit d0d6ef505e

View File

@ -298,22 +298,15 @@ template<class Thermo, template<class> class Type>
inline Foam::scalar inline Foam::scalar
Foam::species::thermo<Thermo, Type>::K(const scalar p, const scalar T) const Foam::species::thermo<Thermo, Type>::K(const scalar p, const scalar T) const
{ {
if (equal(this->nMoles(), SMALL)) scalar arg = -this->nMoles()*this->g(p, T)/(this->RR*T);
if (arg < 600.0)
{ {
return 1.0; return exp(arg);
} }
else else
{ {
scalar arg = -this->nMoles()*this->g(p, T)/(this->RR*T); return VGREAT;
if (arg < 600.0)
{
return exp(arg);
}
else
{
return VGREAT;
}
} }
} }