Updates on behaviour of limiters. Oliver Borm.

src/finiteVolume/interpolation/surfaceInterpolation/limitedSchemes/LimitedScheme/NVDTVD.H

@@ -113,9 +113,14 @@ public:
             }
 
             // Stabilise for division
-            gradf = stabilise(gradf, VSMALL);
+            // Changed to SMALL to prevent FPE.  OB, 14/Jan/2011
+            gradf = stabilise(gradf, SMALL);
 
-            return 2*(gradcf/gradf) - 1;
+            // New formulation.  Oliver Borm and Aleks Jemcov
+            // HJ, 13/Jan/2011
+            return max(2*(gradcf/gradf) - 1, 0);
+
+//             return 2*(gradcf/gradf) - 1;
         }
 };
 

src/finiteVolume/interpolation/surfaceInterpolation/limitedSchemes/LimitedScheme/NVDVTVDV.H

@@ -114,9 +114,14 @@ public:
             }
 
             // Stabilise for division
-            gradf = stabilise(gradf, VSMALL);
+            //Changed to SMALL to prevent FPE.  OB, 14/Jan/2011
+            gradf = stabilise(gradf, SMALL);
 
-            return 2*(gradcf/gradf) - 1;
+            // New formulation.  Oliver Borm and Aleks Jemcov
+            // HJ, 13/Jan/2011
+            return max(2*(gradcf/gradf) - 1, 0);
+
+//             return 2*(gradcf/gradf) - 1;
         }
 };
 

src/finiteVolume/interpolation/surfaceInterpolation/limitedSchemes/upwind/upwind.H

@@ -121,7 +121,7 @@ public:
                         this->mesh()
                     ),
                     this->mesh(),
-                    dimensionedScalar("upwindLimiter", dimless, 0.0)
+                    dimless
                 )
             );
         }

src/finiteVolume/interpolation/surfaceInterpolation/limitedSchemes/vanAlbada/vanAlbada.H

@@ -78,7 +78,11 @@ public:
             faceFlux, phiP, phiN, gradcP, gradcN, d
         );
 
-        return r*(r + 1)/(sqr(r) + 1);
+        // New formulation.  Oliver Borm and Aleks Jemcov
+        // HJ, 13/Jan/2011
+        return (r + 1)/(r + 1/stabilise(r, VSMALL));
+
+//         return r*(r + 1)/(sqr(r) + 1);
     }
 };
 

src/finiteVolume/interpolation/surfaceInterpolation/schemes/linear/linear.H

@@ -120,7 +120,7 @@ template<class Type>
 tmp<GeometricField<Type, fvsPatchField, surfaceMesh> >
 linearInterpolate(const tmp<GeometricField<Type, fvPatchField, volMesh> >& tvf)
 {
-    tmp<GeometricField<Type, fvsPatchField, surfaceMesh> > tinterp = 
+    tmp<GeometricField<Type, fvsPatchField, surfaceMesh> > tinterp =
         linearInterpolate(tvf());
     tvf.clear();
     return tinterp;