diff --git a/src/regionModels/pyrolysisModels/reactingOneDim/reactingOneDim.C b/src/regionModels/pyrolysisModels/reactingOneDim/reactingOneDim.C
index db31437057b4f1c67210363d43f2e9cd2c4c8e53..10c8a4b4ac440af69bb2e30ead371a0a0c1bd4c4 100644
--- 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_)