From 2c655680a60f9b07264fa5c0794695977dc46ffd Mon Sep 17 00:00:00 2001
From: ignis <ignis@postech.ac.kr>
Date: Thu, 21 Jun 2018 01:08:37 +0900
Subject: [PATCH] Handling exception of reactant with zero stoichiometric
 coefficient

---
 .../chemistryModel/chemistryModel/chemistryModel.C              | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

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<CompType, ThermoType>::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);