ENH: Adding heat transfer term to the hEq due to energy transfer to the gas

src/regionModels/pyrolysisModels/reactingOneDim/reactingOneDim.C

@@ -323,6 +323,7 @@ void reactingOneDim::solveEnergy()
       - fvm::laplacian(alpha, h_)
      ==
         chemistrySh_
+      + fvm::Sp(solidChemistry_->RRg(), h_)
     );
 
     if (gasHSource_)