zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

specieThermo: Add guard in K(T) to protect against 0 moles in composition

Browse Source
This commit is contained in:
Henry
2012-10-04 18:10:26 +01:00
parent e9e09f11dc
commit c6eea6990f
1 changed files with 12 additions and 5 deletions
Download Patch File Download Diff File Expand all files Collapse all files

17
src/thermophysicalModels/specie/thermo/thermo/thermoI.H
View File

@ -298,15 +298,22 @@ 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
{ {
scalar arg = -this->nMoles()*this->g(p, T)/(this->RR*T); if (equal(this->nMoles(), SMALL))
if (arg < 600.0)
{ {
return ::exp(arg); return 1.0;
} }
else else
{ {
return VGREAT; scalar arg = -this->nMoles()*this->g(p, T)/(this->RR*T);
if (arg < 600.0)
{
return exp(arg);
}
else
{
return VGREAT;
}
} }
} }