Changing header in tutorial files and updating pEq.H

applications/solvers/multiphase/MPPICInterFoam/pEqn.H

@@ -1,11 +1,10 @@
 {
-    volVectorField HbyA("HbyA", U);
-    HbyA = rAUc*UEqn.H();
+    volVectorField HbyA(constrainHbyA(rAUc*UEqn.H(), U, p_rgh));
 
     surfaceScalarField phiHbyA
     (
         "phiHbyA",
-        (fvc::interpolate(HbyA) & mesh.Sf())
+        fvc::flux(HbyA)
       + alphacf*fvc::interpolate(rho*rAUc)*fvc::ddtCorr(U, phi)
     );
 
@@ -23,15 +22,8 @@
 
     phiHbyA += phig;
 
-    // Update the fixedFluxPressure BCs to ensure flux consistency
-    setSnGrad<fixedFluxPressureFvPatchScalarField>
-    (
-        p_rgh.boundaryField(),
-        (
-            phiHbyA.boundaryField()
-          - MRF.relative(mesh.Sf().boundaryField() & U.boundaryField())
-        )/(mesh.magSf().boundaryField()*rAUcf.boundaryField())
-    );
+    // Update the pressure BCs to ensure flux consistency
+    constrainPressure(p_rgh, U, phiHbyA, rAUcf, MRF);
 
     while (pimple.correctNonOrthogonal())
     {