diff --git a/applications/solvers/basic/potentialFoam/createFields.H b/applications/solvers/basic/potentialFoam/createFields.H
index 97801c0221..40db86dacf 100644
--- a/applications/solvers/basic/potentialFoam/createFields.H
+++ b/applications/solvers/basic/potentialFoam/createFields.H
@@ -40,7 +40,7 @@ word pName("p");
 // Update name of the pressure field from the command-line option
 args.optionReadIfPresent("pName", pName);
 
-// Infer the pressure BCs from the velocity BCs
+// Infer the pressure BCs from the velocity
 wordList pBCTypes
 (
     U.boundaryField().size(),
@@ -71,6 +71,21 @@ volScalarField p
     pBCTypes
 );
 
+// Infer the velocity potential BCs from the pressure
+wordList PhiBCTypes
+(
+    p.boundaryField().size(),
+    zeroGradientFvPatchScalarField::typeName
+);
+
+forAll(p.boundaryField(), patchi)
+{
+    if (p.boundaryField()[patchi].fixesValue())
+    {
+        PhiBCTypes[patchi] = fixedValueFvPatchScalarField::typeName;
+    }
+}
+
 Info<< "Constructing velocity potential field Phi\n" << endl;
 volScalarField Phi
 (
@@ -84,7 +99,7 @@ volScalarField Phi
     ),
     mesh,
     dimensionedScalar("Phi", dimLength*dimVelocity, 0),
-    p.boundaryField().types()
+    PhiBCTypes
 );
 
 label PhiRefCell = 0;