ENH: Adding momemtum predictor to chtSimpleFoam

2025-11-28 03:28:01 +00:00 · 2017-11-06 09:03:42 -08:00
parent 0eeffc318b
commit 454e609fbc
2 changed files with 17 additions and 11 deletions
--- 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);
+    }
--- 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);