Modification of hEq in the pyrolysis model to account for gas enthalpy released

03c28e5d · Sergio Ferraris · a3bad4c2 · 03c28e5d
Commit 03c28e5d authored 11 years ago by Sergio Ferraris
--- a/src/regionModels/pyrolysisModels/reactingOneDim/reactingOneDim.C
+++ b/src/regionModels/pyrolysisModels/reactingOneDim/reactingOneDim.C
@@ -323,7 +323,7 @@ void reactingOneDim::solveEnergy()
      - fvm::laplacian(alpha, h_)
     ==
        chemistrySh_
-      + fvm::Sp(solidChemistry_->RRg(), h_)
+      - fvm::Sp(solidChemistry_->RRg(), h_)
    );

    if (gasHSource_)