BUG: In these solvers adjustPhi in the pEqn was calculated "before"

fvOption.relativeFlux() statement. Therefore adjustPhi was using
the absolute phi instead if the relative when the MRF was active

applications/solvers/basic/potentialFoam/potentialFoam.C

@@ -59,9 +59,10 @@ int main(int argc, char *argv[])
     // function objects so do it ourselves
     runTime.functionObjects().start();
 
+    fvOptions.relativeFlux(phi);
+
     adjustPhi(phi, U, p);
 
-    fvOptions.relativeFlux(phi);
 
     for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
     {

applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/fluid/pEqn.H

@@ -19,12 +19,12 @@
         fvc::interpolate(rho)*(fvc::interpolate(HbyA) & mesh.Sf())
     );
 
+    fvOptions.relativeFlux(fvc::interpolate(rho), phiHbyA);
+
     bool closedVolume = adjustPhi(phiHbyA, U, p_rgh);
 
     phiHbyA += phig;
 
-    fvOptions.relativeFlux(fvc::interpolate(rho), phiHbyA);
-
     dimensionedScalar compressibility = fvc::domainIntegrate(psi);
     bool compressible = (compressibility.value() > SMALL);
 

applications/solvers/incompressible/pimpleFoam/pEqn.H

@@ -13,10 +13,10 @@ surfaceScalarField phiHbyA
   + fvc::ddtPhiCorr(rAU, U, phi)
 );
 
-adjustPhi(phiHbyA, U, p);
-
 fvOptions.relativeFlux(phiHbyA);
 
+adjustPhi(phiHbyA, U, p);
+
 // Non-orthogonal pressure corrector loop
 while (pimple.correctNonOrthogonal())
 {

applications/solvers/incompressible/simpleFoam/pEqn.H

@@ -5,10 +5,11 @@
     UEqn.clear();
 
     surfaceScalarField phiHbyA("phiHbyA", fvc::interpolate(HbyA) & mesh.Sf());
-    adjustPhi(phiHbyA, U, p);
 
     fvOptions.relativeFlux(phiHbyA);
 
+    adjustPhi(phiHbyA, U, p);
+
     // Non-orthogonal pressure corrector loop
     while (simple.correctNonOrthogonal())
     {

applications/solvers/incompressible/simpleFoam/porousSimpleFoam/pEqn.H

@@ -10,10 +10,11 @@ else
 
 UEqn.clear();
 surfaceScalarField phiHbyA("phiHbyA", fvc::interpolate(HbyA) & mesh.Sf());
-adjustPhi(phiHbyA, U, p);
 
 mrfZones.relativeFlux(phiHbyA);
 
+adjustPhi(phiHbyA, U, p);
+
 while (simple.correctNonOrthogonal())
 {
     tmp<fvScalarMatrix> tpEqn;