ENH: Modification of hEqn in reactingOneDim.C introducing laplacian(kappa, T)

instead of laplacian(alpha, h). This is necessary when Cp is not uniform.

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_)