diff --git a/applications/solvers/incompressible/potentialFreeSurfaceFoam/resetPhiPatches.H b/applications/solvers/incompressible/potentialFreeSurfaceFoam/resetPhiPatches.H
new file mode 100644
index 0000000000000000000000000000000000000000..fe428b6e3bf606c25bacd523196415f6ff618a45
--- /dev/null
+++ b/applications/solvers/incompressible/potentialFreeSurfaceFoam/resetPhiPatches.H
@@ -0,0 +1,8 @@
+FieldField<fvsPatchField, scalar>& phibf = phi.boundaryField();
+const FieldField<fvPatchField, vector>& Ubf = U.boundaryField();
+const FieldField<fvsPatchField, vector>& Sfbf = mesh.Sf().boundaryField();
+
+forAll(phibf, patchI)
+{
+    phibf[patchI] = (Ubf[patchI] & Sfbf[patchI]);
+}