LTSInterFoam: Improved temporal damping

2013-09-26 22:39:39 +01:00 · 2013-09-26 22:39:39 +01:00 · ca6fb03260
commit ca6fb03260
parent 2e846f260d
1 changed files with 5 additions and 3 deletions
--- a/applications/solvers/multiphase/interFoam/LTSInterFoam/setrDeltaT.H
+++ b/applications/solvers/multiphase/interFoam/LTSInterFoam/setrDeltaT.H
@ -121,9 +121,11 @@
     && runTime.timeIndex() > runTime.startTimeIndex() + 1
    )
    {
-        rDeltaT =
-             rDeltaT0
-            *max(rDeltaT/rDeltaT0, scalar(1.0) - rDeltaTDampingCoeff);
+        rDeltaT = max
+        (
+            rDeltaT,
+            (scalar(1.0) - rDeltaTDampingCoeff)*rDeltaT0
+        );

        Info<< "Damped flow time scale min/max = "
            << gMin(1/rDeltaT.internalField())