diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/pU/UEqns.H b/applications/solvers/multiphase/twoPhaseEulerFoam/pU/UEqns.H
index 8e4af0811297df8dcc23749da67af70706a0f9d8..bc8ad0f556eea8f8da0d8bdd9f0415f7c7ced230 100644
--- a/applications/solvers/multiphase/twoPhaseEulerFoam/pU/UEqns.H
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/pU/UEqns.H
@@ -1,3 +1,5 @@
+Info<< "Constructing momentum equations" << endl;
+
 MRF.correctBoundaryVelocity(U1);
 MRF.correctBoundaryVelocity(U2);
 MRF.correctBoundaryVelocity(U);
diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/pUf/UEqns.H b/applications/solvers/multiphase/twoPhaseEulerFoam/pUf/UEqns.H
index 7fccb4098e48aafc9d55fd984e665dc0655df3de..257008f7ad8d6cabe3cb6b6485ecec912eeca1ea 100644
--- a/applications/solvers/multiphase/twoPhaseEulerFoam/pUf/UEqns.H
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/pUf/UEqns.H
@@ -1,3 +1,5 @@
+Info<< "Constructing face momentum equations" << endl;
+
 MRF.correctBoundaryVelocity(U1);
 MRF.correctBoundaryVelocity(U2);
 MRF.correctBoundaryVelocity(U);
diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/pUf/pEqn.H b/applications/solvers/multiphase/twoPhaseEulerFoam/pUf/pEqn.H
index d269577591286a068c8b5191cedc9bd5f2f0ee04..ab9a23197cd0b4268845e1ee11de5942231c94f6 100644
--- a/applications/solvers/multiphase/twoPhaseEulerFoam/pUf/pEqn.H
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/pUf/pEqn.H
@@ -358,5 +358,3 @@ if (thermo1.dpdt() || thermo2.dpdt())
 {
     dpdt = fvc::ddt(p);
 }
-
-#include "pUf/DDtU.H"
diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/twoPhaseEulerFoam.C b/applications/solvers/multiphase/twoPhaseEulerFoam/twoPhaseEulerFoam.C
index ffdbd2c0f92198245216b2eb1230e4e90de99697..f2bd6d9a9c7e1951e1c309d000c69dc18aec8f21 100644
--- a/applications/solvers/multiphase/twoPhaseEulerFoam/twoPhaseEulerFoam.C
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/twoPhaseEulerFoam.C
@@ -93,18 +93,18 @@ int main(int argc, char *argv[])
             fluid.correct();
 
             #include "contErrs.H"
-            #include "EEqns.H"
 
             if (faceMomentum)
             {
-                Info<< "Constructing face momentum equations" << endl;
                 #include "pUf/UEqns.H"
+                #include "EEqns.H"
                 #include "pUf/pEqn.H"
+                #include "pUf/DDtU.H"
             }
             else
             {
-                Info<< "Constructing momentum equations" << endl;
                 #include "pU/UEqns.H"
+                #include "EEqns.H"
                 #include "pU/pEqn.H"
                 #include "pU/DDtU.H"
             }