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

Commenting out sanity check in `RedlichKwongMFTP::pressure' and removing m_Pcurrent

Browse source
This commit is contained in:
Steven DeCaluwe 2017-02-17 11:59:01 -07:00 • committed by Ray Speth
parent 11a0727d5c
commit ecbd819e91
1 changed files with 5 additions and 5 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

10
src/thermo/RedlichKwongMFTP.cpp
View file

@ -171,11 +171,11 @@ doublereal RedlichKwongMFTP::pressure() const
doublereal T = temperature();
double molarV = meanMolecularWeight() / density();
double pp = GasConstant * T/(molarV - m_b_current) - m_a_current/(sqrt(T) * molarV * (molarV + m_b_current));
if (fabs(pp -m_Pcurrent) > 1.0E-5 * fabs(m_Pcurrent)) {
throw CanteraError(" RedlichKwongMFTP::pressure()", "setState broken down, maybe");
}
// if (fabs(pp -m_Pcurrent) > 1.0E-5 * fabs(m_Pcurrent)) {
// throw CanteraError(" RedlichKwongMFTP::pressure()", "setState broken down, maybe");
// }
return m_Pcurrent;
return pp;
}
void RedlichKwongMFTP::calcDensity()
@ -426,7 +426,7 @@ void RedlichKwongMFTP::getPartialMolarVolumes(doublereal* vbar) const
- 2.0 * m_pp[k] / (sqt * vpb)
+ m_a_current * b_vec_Curr_[k] / (sqt * vpb * vpb)
);
doublereal denom = (m_Pcurrent + RT() * m_b_current/(vmb * vmb) - m_a_current / (sqt * vpb * vpb)
doublereal denom = (pressure() + RT() * m_b_current/(vmb * vmb) - m_a_current / (sqt * vpb * vpb)
);
vbar[k] = num / denom;
}