From 45a4fbda2c607646a193cc530eca90e71b33693b Mon Sep 17 00:00:00 2001
From: andy <andy>
Date: Wed, 28 Sep 2011 13:58:08 +0100
Subject: [PATCH] ENH: Added ddtPhiCorr to pimpleDyMFoam on a switch

---
 .../solvers/incompressible/pimpleFoam/pimpleDyMFoam/pEqn.H   | 5 +++++
 .../incompressible/pimpleFoam/pimpleDyMFoam/readControls.H   | 2 ++
 2 files changed, 7 insertions(+)

diff --git a/applications/solvers/incompressible/pimpleFoam/pimpleDyMFoam/pEqn.H b/applications/solvers/incompressible/pimpleFoam/pimpleDyMFoam/pEqn.H
index 7e84ef50fe..3e8a587768 100644
--- a/applications/solvers/incompressible/pimpleFoam/pimpleDyMFoam/pEqn.H
+++ b/applications/solvers/incompressible/pimpleFoam/pimpleDyMFoam/pEqn.H
@@ -7,6 +7,11 @@ if (pimple.nCorr() <= 1)
 
 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/incompressible/pimpleFoam/pimpleDyMFoam/readControls.H b/applications/solvers/incompressible/pimpleFoam/pimpleDyMFoam/readControls.H
index 8f48f5d7d4..2ad7f610f0 100644
--- a/applications/solvers/incompressible/pimpleFoam/pimpleDyMFoam/readControls.H
+++ b/applications/solvers/incompressible/pimpleFoam/pimpleDyMFoam/readControls.H
@@ -8,3 +8,5 @@
     const bool checkMeshCourantNo =
         pimpleDict.lookupOrDefault("checkMeshCourantNo", false);
 
+    const bool ddtPhiCorr =
+        pimpleDict.lookupOrDefault("ddtPhiCorr", true);