diff --git a/applications/solvers/basic/potentialDyMFoam/potentialDyMFoam.C b/applications/solvers/basic/potentialDyMFoam/potentialDyMFoam.C
index d68eb5680..44945b07b 100644
--- a/applications/solvers/basic/potentialDyMFoam/potentialDyMFoam.C
+++ b/applications/solvers/basic/potentialDyMFoam/potentialDyMFoam.C
@@ -72,6 +72,11 @@ int main(int argc, char *argv[])
             U.internalField() = vector::zero;
         }
 
+        p.correctBoundaryConditions();
+        U.correctBoundaryConditions();
+
+        phi == (fvc::interpolate(U) & mesh.Sf());
+
 #       include "volContinuity.H"
 #       include "meshCourantNo.H"