From d730c49e86e6a5aab4e5a2ce3312dc465b311e93 Mon Sep 17 00:00:00 2001
From: Sergio Ferraris <sergio@alex.opencfd.co.uk>
Date: Mon, 2 Sep 2013 09:56:39 +0100
Subject: [PATCH] ENH: Modification of hEqn in reactingOneDim.C introducing
 laplacian(kappa, T) instead of laplacian(alpha, h). This is necessary when Cp
 is not uniform.

---
 .../pyrolysisModels/reactingOneDim/reactingOneDim.C             | 2 ++
 1 file changed, 2 insertions(+)

diff --git a/src/regionModels/pyrolysisModels/reactingOneDim/reactingOneDim.C b/src/regionModels/pyrolysisModels/reactingOneDim/reactingOneDim.C
index 6fa334158a..36ff9f4cad 100644
--- a/src/regionModels/pyrolysisModels/reactingOneDim/reactingOneDim.C
+++ b/src/regionModels/pyrolysisModels/reactingOneDim/reactingOneDim.C
@@ -324,6 +324,8 @@ void reactingOneDim::solveEnergy()
     (
         fvm::ddt(rho_, h_)
       - fvm::laplacian(alpha, h_)
+      + fvc::laplacian(alpha, h_)
+      - fvc::laplacian(kappa(), T())
      ==
         chemistrySh_
       - fvm::Sp(solidChemistry_->RRg(), h_)