From 665f5243b1dd85f3dd2d9be55724e42a29235ff4 Mon Sep 17 00:00:00 2001
From: Henry Weller <http://cfd.direct>
Date: Thu, 17 Mar 2016 08:05:08 +0000
Subject: [PATCH] reactingEulerFoam: Moved fluid.correctThermo() into energy
 corrector loop to update the interfacial temperature. Resolves bug-report
 http://www.openfoam.org/mantisbt/view.php?id=2029

---
 .../reactingEulerFoam/reactingMultiphaseEulerFoam/EEqns.H     | 3 ++-
 .../reactingEulerFoam/reactingTwoPhaseEulerFoam/EEqns.H       | 4 ++--
 2 files changed, 4 insertions(+), 3 deletions(-)

diff --git a/applications/solvers/multiphase/reactingEulerFoam/reactingMultiphaseEulerFoam/EEqns.H b/applications/solvers/multiphase/reactingEulerFoam/reactingMultiphaseEulerFoam/EEqns.H
index c1d98bcb439..0130a872c90 100644
--- a/applications/solvers/multiphase/reactingEulerFoam/reactingMultiphaseEulerFoam/EEqns.H
+++ b/applications/solvers/multiphase/reactingEulerFoam/reactingMultiphaseEulerFoam/EEqns.H
@@ -33,9 +33,10 @@ for (int Ecorr=0; Ecorr<nEnergyCorrectors; Ecorr++)
             EEqn->solve();
         }
     }
+
+    fluid.correctThermo();
 }
 
-fluid.correctThermo();
 
 forAll(phases, phasei)
 {
diff --git a/applications/solvers/multiphase/reactingEulerFoam/reactingTwoPhaseEulerFoam/EEqns.H b/applications/solvers/multiphase/reactingEulerFoam/reactingTwoPhaseEulerFoam/EEqns.H
index 0ac8e841297..0cdbfba8b92 100644
--- a/applications/solvers/multiphase/reactingEulerFoam/reactingTwoPhaseEulerFoam/EEqns.H
+++ b/applications/solvers/multiphase/reactingEulerFoam/reactingTwoPhaseEulerFoam/EEqns.H
@@ -47,9 +47,9 @@ for (int Ecorr=0; Ecorr<nEnergyCorrectors; Ecorr++)
             E2eqn->solve();
         }
     }
-}
 
-fluid.correctThermo();
+    fluid.correctThermo();
+}
 
 Info<< phase1.name() << " min/max T "
     << min(phase1.thermo().T()).value()
-- 
GitLab