Handling of U under-relaxation in last corrector

applications/solvers/compressible/rhoPimpleFoam/UEqn.H

@@ -8,9 +8,16 @@ tmp<fvVectorMatrix> UEqn
 );
 
 if (oCorr == nOuterCorr - 1)
+{
+    if (mesh.solutionDict().relax("UFinal"))
+    {
+        UEqn().relax(mesh.solutionDict().relaxationFactor("UFinal"));
+    }
+    else
     {
         UEqn().relax(1);
     }
+}
 else
 {
     UEqn().relax();

applications/solvers/compressible/rhoPorousMRFPimpleFoam/UEqn.H

@@ -8,9 +8,16 @@ tmp<fvVectorMatrix> UEqn
 );
 
 if (oCorr == nOuterCorr-1)
+{
+    if (mesh.solutionDict().relax("UFinal"))
+    {
+        UEqn().relax(mesh.solutionDict().relaxationFactor("UFinal"));
+    }
+    else
     {
         UEqn().relax(1);
     }
+}
 else
 {
     UEqn().relax();

applications/solvers/incompressible/pimpleDyMFoam/UEqn.H

@@ -5,13 +5,20 @@
       + turbulence->divDevReff(U)
     );
 
-    if (ocorr != nOuterCorr - 1)
+    if (ocorr == nOuterCorr - 1)
     {
-        UEqn.relax();
+        if (mesh.solutionDict().relax("UFinal"))
+        {
+            UEqn.relax(mesh.solutionDict().relaxationFactor("UFinal"));
         }
         else
         {
-        UEqn.relax(mesh.solutionDict().relaxationFactor("UFinal"));
+            UEqn.relax(1);
+        }
+    }
+    else
+    {
+        UEqn.relax();
     }
 
     if (momentumPredictor)

applications/solvers/incompressible/pimpleFoam/UEqn.H

@@ -8,9 +8,16 @@ tmp<fvVectorMatrix> UEqn
 );
 
 if (oCorr == nOuterCorr - 1)
+{
+    if (mesh.solutionDict().relax("UFinal"))
+    {
+        UEqn().relax(mesh.solutionDict().relaxationFactor("UFinal"));
+    }
+    else
     {
         UEqn().relax(1);
     }
+}
 else
 {
     UEqn().relax();

applications/solvers/multiphase/twoLiquidMixingFoam/UEqn.H

@@ -15,9 +15,16 @@
     );
 
     if (oCorr == nOuterCorr-1)
+    {
+        if (mesh.solutionDict().relax("UFinal"))
+        {
+            UEqn.relax(mesh.solutionDict().relaxationFactor("UFinal"));
+        }
+        else
         {
             UEqn.relax(1);
         }
+    }
     else
     {
         UEqn.relax();