diff --git a/applications/solvers/incompressible/icoDyMFoam/icoDyMFoam.C b/applications/solvers/incompressible/icoDyMFoam/icoDyMFoam.C
index 814e755326fe84ddd10355e73bb2a123bb400294..5b8d15c544ed49c036b92d41aa1059902ebf42cb 100644
--- a/applications/solvers/incompressible/icoDyMFoam/icoDyMFoam.C
+++ b/applications/solvers/incompressible/icoDyMFoam/icoDyMFoam.C
@@ -57,15 +57,15 @@ int main(int argc, char *argv[])
 #       include "readControls.H"
 #       include "CourantNo.H"
 
+        // Make the fluxes absolute
+        fvc::makeAbsolute(phi, U);
+
 #       include "setDeltaT.H"
 
         runTime++;
 
         Info<< "Time = " << runTime.timeName() << nl << endl;
 
-        // Make the fluxes absolute
-        fvc::makeAbsolute(phi, U);
-
         mesh.update();
 
         if (mesh.changing() && correctPhi)
@@ -96,11 +96,6 @@ int main(int argc, char *argv[])
                 U = rAU*UEqn.H();
                 phi = (fvc::interpolate(U) & mesh.Sf());
 
-                if (ddtPhiCorr)
-                {
-                    phi += fvc::ddtPhiCorr(rAU, U, phi);
-                }
-
                 if (p.needReference())
                 {
                     fvc::makeRelative(phi, U);
diff --git a/applications/solvers/multiphase/interDyMFoam/interDyMFoam.C b/applications/solvers/multiphase/interDyMFoam/interDyMFoam.C
index 9d6b9c4cd9e986d412e1340fdb6e1bb85c21965f..053a9fa813a3ec26642d666cea1a801fef70f7d2 100644
--- a/applications/solvers/multiphase/interDyMFoam/interDyMFoam.C
+++ b/applications/solvers/multiphase/interDyMFoam/interDyMFoam.C
@@ -67,15 +67,16 @@ int main(int argc, char *argv[])
     {
         #include "readControls.H"
         #include "CourantNo.H"
+
+        // Make the fluxes absolute
+        fvc::makeAbsolute(phi, U);
+
         #include "setDeltaT.H"
 
         runTime++;
 
         Info<< "Time = " << runTime.timeName() << nl << endl;
 
-        // Make the fluxes absolute
-        fvc::makeAbsolute(phi, U);
-
         scalar timeBeforeMeshUpdate = runTime.elapsedCpuTime();
 
         // Do any mesh changes
@@ -96,9 +97,6 @@ int main(int argc, char *argv[])
             #include "correctPhi.H"
         }
 
-        // Keep the absolute fluxes for use in ddtPhiCorr
-        surfaceScalarField phiAbs0("phiAbs0", phi);
-
         // Make the fluxes relative to the mesh motion
         fvc::makeRelative(phi, U);
 
diff --git a/applications/solvers/multiphase/interDyMFoam/pEqn.H b/applications/solvers/multiphase/interDyMFoam/pEqn.H
index 0fdba99cd696f12bfec9b17340dab5f75388b893..c6961f6c22c3e2be738c1ae4de1226fa431deb0b 100644
--- a/applications/solvers/multiphase/interDyMFoam/pEqn.H
+++ b/applications/solvers/multiphase/interDyMFoam/pEqn.H
@@ -5,11 +5,6 @@
     U = rAU*UEqn.H();
     surfaceScalarField phiU("phiU", (fvc::interpolate(U) & mesh.Sf()));
 
-    if (ddtPhiCorr)
-    {
-        phiU += fvc::ddtPhiCorr(rAU, rho, U, phiAbs0);
-    }
-
     phi = phiU +
     (
         fvc::interpolate(interface.sigmaK())*fvc::snGrad(gamma)