BUG: Fixing solid specie index in omega function

74430112 · sergio · Andrew Heather · 78524117 · 74430112
Commit 74430112 authored Aug 21, 2019 by sergio Committed by Andrew Heather Aug 21, 2019
Hide whitespace changes
Inline Side-by-side

Showing with 2 additions and 2 deletions

src/thermophysicalModels/solidChemistryModel/pyrolysisChemistryModel/pyrolysisChemistryModel.C ...ryModel/pyrolysisChemistryModel/pyrolysisChemistryModel.C +2 -2

No files found.
--- a/src/thermophysicalModels/solidChemistryModel/pyrolysisChemistryModel/pyrolysisChemistryModel.C
+++ b/src/thermophysicalModels/solidChemistryModel/pyrolysisChemistryModel/pyrolysisChemistryModel.C
@@ -2,7 +2,7 @@
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     |
-    \\  /    A nd           |
+    \\  /    A nd           | Copyright (C) 2019 OpenCFD Ltd.
     \\/     M anipulation  |
 -------------------------------------------------------------------------------
                            | Copyright (C) 2013-2017 OpenFOAM Foundation
@@ -264,7 +264,7 @@ Foam::pyrolysisChemistryModel<CompType, SolidThermo, GasThermo>::omega
    for (label s=0; s<Nl; s++)
    {
        label si = R.lhs()[s].index;
-        const scalar exp = R.lhs()[si].exponent;
+        const scalar exp = R.lhs()[s].exponent;

        kf *=
            pow(c1[si]/Ys0_[si][celli], exp)