diff --git a/applications/solvers/multiphase/reactingEulerFoam/reactingMultiphaseEulerFoam/pU/pEqn.H b/applications/solvers/multiphase/reactingEulerFoam/reactingMultiphaseEulerFoam/pU/pEqn.H
index 3fe95ee8b2..372916fbc0 100644
--- a/applications/solvers/multiphase/reactingEulerFoam/reactingMultiphaseEulerFoam/pU/pEqn.H
+++ b/applications/solvers/multiphase/reactingEulerFoam/reactingMultiphaseEulerFoam/pU/pEqn.H
@@ -141,9 +141,12 @@ while (pimple.correct())
                 new surfaceScalarField
                 (
                     IOobject::groupName("phiHbyA", phase.name()),
-                    fvc::flux(HbyAs[phase.index()])
-                  - phigFs[phase.index()]
-                  - ddtCorrByAs[phase.index()]
+                    MRF.relative
+                    (
+                        fvc::flux(HbyAs[phase.index()])
+                      - phigFs[phase.index()]
+                      - ddtCorrByAs[phase.index()]
+                    )
                 )
             );
         }
@@ -204,8 +207,6 @@ while (pimple.correct())
         }
     }
 
-    MRF.makeRelative(phiHbyA);
-
     // Construct pressure "diffusivity"
     surfaceScalarField rAUf
     (
diff --git a/applications/solvers/multiphase/reactingEulerFoam/reactingMultiphaseEulerFoam/pUf/pEqn.H b/applications/solvers/multiphase/reactingEulerFoam/reactingMultiphaseEulerFoam/pUf/pEqn.H
index 4999134331..4dab2a37af 100644
--- a/applications/solvers/multiphase/reactingEulerFoam/reactingMultiphaseEulerFoam/pUf/pEqn.H
+++ b/applications/solvers/multiphase/reactingEulerFoam/reactingMultiphaseEulerFoam/pUf/pEqn.H
@@ -128,13 +128,16 @@ while (pimple.correct())
             new surfaceScalarField
             (
                 IOobject::groupName("phiHbyA", phase.name()),
-                rAUfs[phase.index()]
-               *(
-                    fvc::flux(UEqns[phase.index()].H())
-                  + alphaRho0fs[phase.index()]
-                   *byDt(MRF.absolute(phase.phi()().oldTime()))
+                MRF.relative
+                (
+                    rAUfs[phase.index()]
+                   *(
+                        fvc::flux(UEqns[phase.index()].H())
+                      + alphaRho0fs[phase.index()]
+                       *byDt(MRF.absolute(phase.phi()().oldTime()))
+                    )
+                  - phigFs[phase.index()]
                 )
-              - phigFs[phase.index()]
             )
         );
     }
@@ -188,8 +191,6 @@ while (pimple.correct())
         }
     }
 
-    MRF.makeRelative(phiHbyA);
-
     // Construct pressure "diffusivity"
     surfaceScalarField rAUf
     (