From f719e7acb26f277e95cb3c4adb31e411de7b0e47 Mon Sep 17 00:00:00 2001
From: Vuko Vukcevic <vuko.vukcevic@fsb.hr>
Date: Thu, 22 Sep 2016 14:46:12 +0200
Subject: [PATCH] Updates in pimpleDyMOversetFoam

Solving pEqn to a tight tolerance in all last PISO correctors.
---
 src/pimpleDyMOversetFoam/Make/options |  6 ++++--
 src/pimpleDyMOversetFoam/pEqn.H       | 18 +++++++++++++++---
 2 files changed, 19 insertions(+), 5 deletions(-)

diff --git a/src/pimpleDyMOversetFoam/Make/options b/src/pimpleDyMOversetFoam/Make/options
index 8b8daa1..aab2ae4 100644
--- a/src/pimpleDyMOversetFoam/Make/options
+++ b/src/pimpleDyMOversetFoam/Make/options
@@ -6,7 +6,8 @@ EXE_INC = \
     -I$(LIB_SRC)/transportModels \
     -I$(LIB_SRC)/transportModels/incompressible/singlePhaseTransportModel \
     -I$(LIB_SRC)/finiteVolume/lnInclude \
-    -I../oversetMesh/lnInclude
+    -I../oversetMesh/lnInclude \
+    -I../oversetDynamicFvMesh/lnInclude
 
 EXE_LIBS = \
     -ldynamicFvMesh \
@@ -20,4 +21,5 @@ EXE_LIBS = \
     -lfiniteVolume \
     -llduSolvers \
     -L$(MESQUITE_LIB_DIR) -lmesquite \
-    -L$(FOAM_USER_LIBBIN) -loversetMesh
+    -L$(FOAM_USER_LIBBIN) -loversetMesh \
+    -L$(FOAM_USER_LIBBIN) -loversetDynamicFvMesh
diff --git a/src/pimpleDyMOversetFoam/pEqn.H b/src/pimpleDyMOversetFoam/pEqn.H
index 904d8a2..866b714 100644
--- a/src/pimpleDyMOversetFoam/pEqn.H
+++ b/src/pimpleDyMOversetFoam/pEqn.H
@@ -31,10 +31,22 @@
         om.correctNonOrthoFluxes(pEqn, U);
 
         pEqn.setReference(pRefCell, pRefValue);
-        pEqn.solve
+
+        if
         (
-            mesh.solutionDict().solver(p.select(pimple.finalInnerIter()))
-        );
+            (
+                pimple.corrPISO() == 1
+             || pimple.corrPISO() == pimple.nCorrPISO()
+            )
+         && pimple.finalNonOrthogonalIter()
+        )
+        {
+            pEqn.solve(mesh.solutionDict().solver("pFinal"));
+        }
+        else
+        {
+            pEqn.solve();
+        }
 
         if (pimple.finalNonOrthogonalIter())
         {