[Kinetics] Limit sticking species to gas phase

Closes #1210
Cantera · Mar 5, 2022 · 57db852 · 57db852
1 parent 367250d
commit 57db852
Show file tree

Hide file tree

Showing 2 changed files with 1 addition and 2 deletions.
diff --git a/src/kinetics/InterfaceRate.cpp b/src/kinetics/InterfaceRate.cpp
@@ -145,7 +145,7 @@ void StickingCoverage::setContext(const Reaction& rxn, const Kinetics& kin)
         bool foundStick = false;
         for (const auto& sp : rxn.reactants) {
             size_t iPhase = kin.speciesPhaseIndex(kin.kineticsSpeciesIndex(sp.first));
-            if (iPhase != iInterface) {
+            if (kin.thermo(iPhase).phaseOfMatter() == "gas" && iPhase != iInterface) {
                 // Non-interface species. There should be exactly one of these
                 if (foundStick) {
                     throw InputFileError("StickingCoverage::setContext",

diff --git a/test/data/haca2.yaml b/test/data/haca2.yaml
@@ -88,4 +88,3 @@ reactions:
   rate-constant: {A: 2.2e+12, b: 0.0, Ea: 7.5}
 - equation: OH + Csoot-H + CB-CB3 => Csoot-* + CO + H2  # Reaction 8
   sticking-coefficient: {A: 0.13, b: 0.0, Ea: 0.0}
-  sticking-species: OH