diff --git a/src/thermophysicalModels/chemistryModel/chemistryModel/chemistryModel/chemistryModel.C b/src/thermophysicalModels/chemistryModel/chemistryModel/chemistryModel/chemistryModel.C index 91f2143d3..52c31eb26 100644 --- a/src/thermophysicalModels/chemistryModel/chemistryModel/chemistryModel/chemistryModel.C +++ b/src/thermophysicalModels/chemistryModel/chemistryModel/chemistryModel/chemistryModel.C @@ -808,7 +808,7 @@ void Foam::chemistryModel::edm forAll(R.lhs(), s) { 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(); minYR = min(minYR, Yi[k]/Wk/vk);