research/OpenFOAM-4.x
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

Handling exception of reactant with zero stoichiometric coefficient

Browse source
This commit is contained in:
ignis 2018-06-21 01:08:37 +09:00
parent cf0d9f904f
commit 2c655680a6
1 changed files with 1 additions and 1 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

2
src/thermophysicalModels/chemistryModel/chemistryModel/chemistryModel/chemistryModel.C
View file

@ -808,7 +808,7 @@ void Foam::chemistryModel<CompType, ThermoType>::edm
forAll(R.lhs(), s) forAll(R.lhs(), s)
{ {
const label k = R.lhs()[s].index; const label k = R.lhs()[s].index;
const scalar vk = R.lhs()[s].stoichCoeff; const scalar vk = R.lhs()[s].stoichCoeff + SMALL;
const scalar Wk = specieThermo_[k].W(); const scalar Wk = specieThermo_[k].W();
minYR = min(minYR, Yi[k]/Wk/vk); minYR = min(minYR, Yi[k]/Wk/vk);