diff --git a/applications/solvers/incompressible/pimpleFoam/UEqn.H b/applications/solvers/incompressible/pimpleFoam/UEqn.H
index 18b6beafe99d7023196061d1de49c9a1a0a3afb8..880c0066e363ad1dee9c52b4b86cf1d844fad65a 100644
--- a/applications/solvers/incompressible/pimpleFoam/UEqn.H
+++ b/applications/solvers/incompressible/pimpleFoam/UEqn.H
@@ -15,8 +15,6 @@ UEqn().relax();
 
 fvOptions.constrain(UEqn());
 
-volScalarField rAU(1.0/UEqn().A());
-
 if (pimple.momentumPredictor())
 {
     solve(UEqn() == -fvc::grad(p));
diff --git a/applications/solvers/incompressible/pimpleFoam/pEqn.H b/applications/solvers/incompressible/pimpleFoam/pEqn.H
index 13118296c5920f4eac34ca0a4bbda42b8b4821bf..18083aa810f548cb01e58108103e8a3ec46d9816 100644
--- a/applications/solvers/incompressible/pimpleFoam/pEqn.H
+++ b/applications/solvers/incompressible/pimpleFoam/pEqn.H
@@ -1,3 +1,4 @@
+volScalarField rAU(1.0/UEqn().A());
 volVectorField HbyA("HbyA", U);
 HbyA = rAU*UEqn().H();
 
diff --git a/applications/solvers/incompressible/pimpleFoam/pimpleDyMFoam/pEqn.H b/applications/solvers/incompressible/pimpleFoam/pimpleDyMFoam/pEqn.H
index 5ad9b80020493835da52896afd949e5237b28cc3..5b5fe720ad8fc9834d304f24a5c63928a6b1424c 100644
--- a/applications/solvers/incompressible/pimpleFoam/pimpleDyMFoam/pEqn.H
+++ b/applications/solvers/incompressible/pimpleFoam/pimpleDyMFoam/pEqn.H
@@ -1,3 +1,4 @@
+volScalarField rAU(1.0/UEqn().A());
 volVectorField HbyA("HbyA", U);
 HbyA = rAU*UEqn().H();
 
diff --git a/applications/solvers/incompressible/pisoFoam/pEqn.H b/applications/solvers/incompressible/pisoFoam/pEqn.H
index 5e4d1560ff71c2943a785dd738bf4aaaceae5df4..f8c5a4597805a6990d1280ce396b11e2d08be92a 100644
--- a/applications/solvers/incompressible/pisoFoam/pEqn.H
+++ b/applications/solvers/incompressible/pisoFoam/pEqn.H
@@ -1,5 +1,4 @@
 volScalarField rAU(1.0/UEqn.A());
-
 volVectorField HbyA("HbyA", U);
 HbyA = rAU*UEqn.H();
 surfaceScalarField phiHbyA