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
2016-03-17 08:05:08 +00:00 · 2016-03-17 08:05:08 +00:00 · 926c29e95a
commit 926c29e95a
parent 156b2d9791
2 changed files with 4 additions and 3 deletions
--- 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)
 {
--- 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()