specie/thermo: Use Pstd in the evaluation of equilibrium coefficients

This commit is contained in:
Henry Weller
2015-11-21 18:29:20 +00:00
parent 37430daa19
commit 6335be9356

View File

@ -281,7 +281,7 @@ template<class Thermo, template<class> class Type>
inline Foam::scalar
Foam::species::thermo<Thermo, Type>::K(const scalar p, const scalar T) const
{
scalar arg = -this->nMoles()*this->g(p, T)/(RR*T);
scalar arg = -this->nMoles()*this->g(Pstd, T)/(RR*T);
if (arg < 600.0)
{