twoPhaseEulerFoam: Correct continuity errors after temperature update

applications/solvers/multiphase/twoPhaseEulerFoam/contErrs.H

@@ -0,0 +1,13 @@
+volScalarField contErr1
+(
+    "contErr1",
+    fvc::ddt(alpha1, rho1) + fvc::div(alphaRhoPhi1)
+  - (fvOptions(alpha1, rho1)&rho1)
+);
+
+volScalarField contErr2
+(
+    "contErr2",
+    fvc::ddt(alpha2, rho2) + fvc::div(alphaRhoPhi2)
+  - (fvOptions(alpha2, rho2)&rho2)
+);

applications/solvers/multiphase/twoPhaseEulerFoam/correctContErrs.H

@@ -0,0 +1,7 @@
+contErr1 =
+    fvc::ddt(alpha1, rho1) + fvc::div(alphaRhoPhi1)
+  - (fvOptions(alpha1, rho1)&rho1);
+
+contErr2 =
+    fvc::ddt(alpha2, rho2) + fvc::div(alphaRhoPhi2)
+  - (fvOptions(alpha2, rho2)&rho2);

applications/solvers/multiphase/twoPhaseEulerFoam/pEqn.H

@@ -10,6 +10,9 @@ surfaceScalarField alpharAU2f(fvc::interpolate(alpha2*rAU2));
 // --- Pressure corrector loop
 while (pimple.correct())
 {
+    // Update continuity errors due to temperature changes
+    #include "correctContErrs.H"
+
     volVectorField HbyA1
     (
         IOobject::groupName("HbyA", phase1.name()),

applications/solvers/multiphase/twoPhaseEulerFoam/twoPhaseEulerFoam.C

@@ -76,20 +76,7 @@ int main(int argc, char *argv[])
             fluid.solve();
             fluid.correct();
 
-            volScalarField contErr1
+            #include "contErrs.H"
-            (
-                "contErr1",
-                fvc::ddt(alpha1, rho1) + fvc::div(alphaRhoPhi1)
-              - (fvOptions(alpha1, rho1)&rho1)
-            );
-
-            volScalarField contErr2
-            (
-                "contErr2",
-                fvc::ddt(alpha2, rho2) + fvc::div(alphaRhoPhi2)
-              - (fvOptions(alpha2, rho2)&rho2)
-            );
-
             #include "UEqns.H"
             #include "EEqns.H"
             #include "pEqn.H"

tutorials/multiphase/twoPhaseEulerFoam/RAS/fluidisedBed/0/T.particles

@@ -27,7 +27,10 @@ boundaryField
 
     outlet
     {
-        type               zeroGradient;
+        type               inletOutlet;
+        phi                phi.particles;
+        inletValue         uniform 300;
+        value              uniform 300;
     }
 
     walls

tutorials/multiphase/twoPhaseEulerFoam/laminar/fluidisedBed/0/T.particles

@@ -27,7 +27,10 @@ boundaryField
 
     outlet
     {
-        type               zeroGradient;
+        type               inletOutlet;
+        phi                phi.particles;
+        inletValue         uniform 300;
+        value              uniform 300;
     }
 
     walls