From 7f8dbac2cc0f054936369231cb274768121b4165 Mon Sep 17 00:00:00 2001
From: Henry <Henry>
Date: Wed, 30 Mar 2011 13:50:09 +0100
Subject: [PATCH 1/2] LTSInterFoam: Improved the near-interface Courant number
 evaluation

---
 .../interFoam/LTSInterFoam/setrDeltaT.H       | 15 ++++----
 .../finiteVolume/fvc/fvcSmooth/fvcSmooth.C    | 35 +++++--------------
 2 files changed, 18 insertions(+), 32 deletions(-)

diff --git a/applications/solvers/multiphase/interFoam/LTSInterFoam/setrDeltaT.H b/applications/solvers/multiphase/interFoam/LTSInterFoam/setrDeltaT.H
index 3a4e6385ec..1e8d8da3e7 100644
--- a/applications/solvers/multiphase/interFoam/LTSInterFoam/setrDeltaT.H
+++ b/applications/solvers/multiphase/interFoam/LTSInterFoam/setrDeltaT.H
@@ -11,7 +11,7 @@
 
     scalar rDeltaTSmoothingCoeff
     (
-        piso.lookupOrDefault<scalar>("rDeltaTSmoothingCoeff", 0.02)
+        piso.lookupOrDefault<scalar>("rDeltaTSmoothingCoeff", 0.1)
     );
 
     label nAlphaSpreadIter
@@ -21,17 +21,17 @@
 
     scalar alphaSpreadDiff
     (
-        piso.lookupOrDefault<label>("alphaSpreadDiff", 0.2)
+        piso.lookupOrDefault<scalar>("alphaSpreadDiff", 0.2)
     );
 
     scalar alphaSpreadMax
     (
-        piso.lookupOrDefault<label>("alphaSpreadMax", 0.99)
+        piso.lookupOrDefault<scalar>("alphaSpreadMax", 0.99)
     );
 
     scalar alphaSpreadMin
     (
-        piso.lookupOrDefault<label>("alphaSpreadMin", 0.01)
+        piso.lookupOrDefault<scalar>("alphaSpreadMin", 0.01)
     );
 
     label nAlphaSweepIter
@@ -63,11 +63,14 @@
     {
         // Further limit the reciprocal time-step
         // in the vicinity of the interface
+
+        volScalarField alpha1Bar = fvc::average(alpha1);
+
         rDeltaT.dimensionedInternalField() = max
         (
             rDeltaT.dimensionedInternalField(),
-            pos(alpha1.dimensionedInternalField() - 0.01)
-           *pos(0.99 - alpha1.dimensionedInternalField())
+            pos(alpha1Bar.dimensionedInternalField() - alphaSpreadMin)
+           *pos(alphaSpreadMax - alpha1Bar.dimensionedInternalField())
            *fvc::surfaceSum(mag(phi))().dimensionedInternalField()
            /((2*maxAlphaCo)*mesh.V())
         );
diff --git a/src/finiteVolume/finiteVolume/fvc/fvcSmooth/fvcSmooth.C b/src/finiteVolume/finiteVolume/fvc/fvcSmooth/fvcSmooth.C
index 94a6a02598..ef22aeee8e 100644
--- a/src/finiteVolume/finiteVolume/fvc/fvcSmooth/fvcSmooth.C
+++ b/src/finiteVolume/finiteVolume/fvc/fvcSmooth/fvcSmooth.C
@@ -154,20 +154,13 @@ void Foam::fvc::spread
         const label own = owner[facei];
         const label nbr = neighbour[facei];
 
-        if
-        (
-            (alpha[own] > alphaMin && alpha[own] < alphaMax)
-         || (alpha[nbr] > alphaMin && alpha[nbr] < alphaMax)
-        )
+        if (mag(alpha[own] - alpha[nbr]) > alphaDiff)
         {
-            if (mag(alpha[own] - alpha[nbr]) > alphaDiff)
-            {
-                changedFaces.append(facei);
-                changedFacesInfo.append
-                (
-                    smoothData(max(field[own], field[nbr]))
-                );
-            }
+            changedFaces.append(facei);
+            changedFacesInfo.append
+            (
+                smoothData(max(field[own], field[nbr]))
+            );
         }
     }
 
@@ -188,20 +181,10 @@ void Foam::fvc::spread
                     alpha.boundaryField()[patchi].patchNeighbourField()
                 );
 
-                if
-                (
-                    (alpha[own] > alphaMin && alpha[own] < alphaMax)
-                 || (
-                        alphapn[patchFacei] > alphaMin
-                     && alphapn[patchFacei] < alphaMax
-                    )
-                )
+                if (mag(alpha[own] - alphapn[patchFacei]) > alphaDiff)
                 {
-                    if (mag(alpha[own] - alphapn[patchFacei]) > alphaDiff)
-                    {
-                        changedFaces.append(facei);
-                        changedFacesInfo.append(smoothData(field[own]));
-                    }
+                    changedFaces.append(facei);
+                    changedFacesInfo.append(smoothData(field[own]));
                 }
             }
         }

From 4802014f7a950297b826388a88331fec805866b2 Mon Sep 17 00:00:00 2001
From: Henry <Henry>
Date: Wed, 30 Mar 2011 13:50:29 +0100
Subject: [PATCH 2/2] LTSInterFoam/wigleyHull: further optimisation

---
 tutorials/multiphase/LTSInterFoam/wigleyHull/0/alpha1.org    | 4 +---
 .../multiphase/LTSInterFoam/wigleyHull/system/controlDict    | 2 +-
 .../multiphase/LTSInterFoam/wigleyHull/system/fvSolution     | 5 +++--
 3 files changed, 5 insertions(+), 6 deletions(-)

diff --git a/tutorials/multiphase/LTSInterFoam/wigleyHull/0/alpha1.org b/tutorials/multiphase/LTSInterFoam/wigleyHull/0/alpha1.org
index cdafeac48a..005c02f6a9 100644
--- a/tutorials/multiphase/LTSInterFoam/wigleyHull/0/alpha1.org
+++ b/tutorials/multiphase/LTSInterFoam/wigleyHull/0/alpha1.org
@@ -29,9 +29,7 @@ boundaryField
 
     outlet
     {
-        type            inletOutlet;
-        inletValue      uniform 0;
-        value           uniform 0;
+        type            zeroGradient;
     }
 
     sides
diff --git a/tutorials/multiphase/LTSInterFoam/wigleyHull/system/controlDict b/tutorials/multiphase/LTSInterFoam/wigleyHull/system/controlDict
index ec7bd439cf..0ebcfacfb3 100644
--- a/tutorials/multiphase/LTSInterFoam/wigleyHull/system/controlDict
+++ b/tutorials/multiphase/LTSInterFoam/wigleyHull/system/controlDict
@@ -23,7 +23,7 @@ startTime       0;
 
 stopAt          endTime;
 
-endTime         3000;
+endTime         2000;
 
 deltaT          1;
 
diff --git a/tutorials/multiphase/LTSInterFoam/wigleyHull/system/fvSolution b/tutorials/multiphase/LTSInterFoam/wigleyHull/system/fvSolution
index f1a31fb6cd..cacd579cd5 100644
--- a/tutorials/multiphase/LTSInterFoam/wigleyHull/system/fvSolution
+++ b/tutorials/multiphase/LTSInterFoam/wigleyHull/system/fvSolution
@@ -83,9 +83,10 @@ PISO
     nAlphaSubCycles 1;
     cAlpha          1;
 
-    maxCo           0.8;
+    maxCo           0.9;
     maxAlphaCo      0.2;
-    rDeltaTSmoothingCoeff 0.02;
+
+    rDeltaTSmoothingCoeff 0.1;
     rDeltaTDampingCoeff 1;
     maxDeltaT       1;
 }