Updates in pimpleDyMOversetFoam

Solving pEqn to a tight tolerance in all last PISO correctors.
Mopolino8 · Sep 22, 2016 · f719e7a · f719e7a
1 parent d6fad2e
commit f719e7a
Show file tree

Hide file tree

Showing 2 changed files with 19 additions and 5 deletions.
diff --git 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
@@ -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())
         {