From 6faccd27d4bd1261fb3164482e1f2419e9d2ad63 Mon Sep 17 00:00:00 2001
From: Sergio Ferraris <sergio@alex.opencfd.co.uk>
Date: Tue, 6 Aug 2013 10:37:35 +0100
Subject: [PATCH] ENH: Adding heat transfer term to the hEq due to energy
 transfer to the gas

---
 src/regionModels/pyrolysisModels/reactingOneDim/reactingOneDim.C | 1 +
 1 file changed, 1 insertion(+)

diff --git a/src/regionModels/pyrolysisModels/reactingOneDim/reactingOneDim.C b/src/regionModels/pyrolysisModels/reactingOneDim/reactingOneDim.C
index 29325f602e..db31437057 100644
--- a/src/regionModels/pyrolysisModels/reactingOneDim/reactingOneDim.C
+++ b/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_)