Made ddtPhiCorr optional in moving-mesh solvers.

Controlled via the "ddtPhiCorr" switch in the PISO sub-dict,
defaults to false.

applications/solvers/incompressible/icoDyMFoam/icoDyMFoam.C

@@ -75,6 +75,9 @@ int main(int argc, char *argv[])
 #           include "correctPhi.H"
         }
 
+        // Keep the absolute fluxes for use in ddtPhiCorr
+        surfaceScalarField phiAbs("phiAbs", phi);
+
         // Make the fluxes relative to the mesh motion
         fvc::makeRelative(phi, U);
 
@@ -95,7 +98,11 @@ int main(int argc, char *argv[])
 
                 U = rAU*UEqn.H();
                 phi = (fvc::interpolate(U) & mesh.Sf());
-                    //+ fvc::ddtPhiCorr(rAU, U, phi);
+
+                if (ddtPhiCorr)
+                {
+                    phi += fvc::ddtPhiCorr(rAU, U, phiAbs);
+                }
 
                 adjustPhi(phi, U, p);
 
@@ -116,7 +123,7 @@ int main(int argc, char *argv[])
                     {
                         pEqn.solve(mesh.solver(p.name()));
                     }
-                
+
                     if (nonOrth == nNonOrthCorr)
                     {
                         phi -= pEqn.flux();

applications/solvers/multiphase/interDyMFoam/pEqn.H

@@ -5,12 +5,12 @@
     volVectorField HU = UEqn.H();
     U = rAU*HU;
 
-    surfaceScalarField phiU
+    surfaceScalarField phiU("phiU", (fvc::interpolate(U) & mesh.Sf()));
-    (
+
-        "phiU",
+    if (ddtPhiCorr)
-        (fvc::interpolate(U) & mesh.Sf())
+    {
-      + fvc::ddtPhiCorr(rAU, rho, U, phiAbs)
+        phiU += fvc::ddtPhiCorr(rAU, rho, U, phiAbs);
-    );
+    }
 
     phi = phiU +
         (

src/finiteVolume/cfdTools/general/include/readPISOControls.H

@@ -25,3 +25,9 @@
     {
         nOuterCorr = readInt(piso.lookup("nOuterCorrectors"));
     }
+
+    bool ddtPhiCorr = false;
+    if (piso.found("ddtPhiCorr"))
+    {
+        ddtPhiCorr = Switch(piso.lookup("ddtPhiCorr"));
+    }