From 454e609fbca000a103c8405eda7109ea93bfedde Mon Sep 17 00:00:00 2001
From: sergio <sergio>
Date: Mon, 6 Nov 2017 09:03:42 -0800
Subject: [PATCH] ENH: Adding momemtum predictor to chtSimpleFoam

---
 .../chtMultiRegionSimpleFoam/fluid/UEqn.H     | 25 +++++++++++--------
 .../readFluidMultiRegionSIMPLEControls.H      |  3 +++
 2 files changed, 17 insertions(+), 11 deletions(-)

diff --git a/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/fluid/UEqn.H b/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/fluid/UEqn.H
index ed0402c2d0..f9712b548e 100644
--- a/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/fluid/UEqn.H
+++ b/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/fluid/UEqn.H
@@ -16,17 +16,20 @@
 
     fvOptions.constrain(UEqn);
 
-    solve
-    (
-        UEqn
-     ==
-        fvc::reconstruct
+    if (momentumPredictor)
+    {
+        solve
         (
+            UEqn
+        ==
+            fvc::reconstruct
             (
-              - ghf*fvc::snGrad(rho)
-              - fvc::snGrad(p_rgh)
-            )*mesh.magSf()
-        )
-    );
+                (
+                - ghf*fvc::snGrad(rho)
+                - fvc::snGrad(p_rgh)
+                )*mesh.magSf()
+            )
+        );
 
-    fvOptions.correct(U);
+        fvOptions.correct(U);
+    }
diff --git a/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/fluid/readFluidMultiRegionSIMPLEControls.H b/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/fluid/readFluidMultiRegionSIMPLEControls.H
index 66fd738d7a..3ce2a2b7b4 100644
--- a/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/fluid/readFluidMultiRegionSIMPLEControls.H
+++ b/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/fluid/readFluidMultiRegionSIMPLEControls.H
@@ -2,3 +2,6 @@
 
     const int nNonOrthCorr =
         simple.lookupOrDefault<int>("nNonOrthogonalCorrectors", 0);
+
+    const bool momentumPredictor =
+        simple.lookupOrDefault("momentumPredictor", true);