interDyMFoam: Update the phi boundary field for consistency for the p_rgh BC

applications/solvers/multiphase/interFoam/interDyMFoam/pEqn.H

@@ -4,6 +4,8 @@
 
     volVectorField HbyA("HbyA", U);
     HbyA = rAU*UEqn.H();
+    phi.boundaryField() =
+        fvc::interpolate(U.boundaryField()) & mesh.Sf().boundaryField();
 
     surfaceScalarField phiHbyA
     (