Merge branch 'master' of /home/noisy3/OpenFOAM/OpenFOAM-dev

ReleaseNotes-dev

@@ -240,11 +240,14 @@
     + =snappyHexMesh=:
       + extrude across multi-processor boundaries
       + preserve faceZones during layering
+      + combining patch faces after snapping
       + *Warning*:
             - minMedianAxisAngle angle fixed w.r.t. 17x. Set to
               90 to get same behaviour as 130 in 17x.
             - nGrow did not work in 17x. Set to 0
               to get same behaviour as 130 in 17x.
+            - nSolveIter reinstated to smooth internal mesh. Set to 0
+              to get same behaviour at 17x.
             - new minTetQuality setting (see tracking). Disable for 17x
               compatibility.
 * Post-processing

applications/utilities/postProcessing/sampling/sample/sampleDict

@@ -150,29 +150,29 @@ surfaces
         interpolate     true;
     }
 
-    movingWall_constant
+    walls_constant
     {
         type            patch;
-        patchName       movingWall;
+        patches         ( ".*Wall.*" );
         // Optional: whether to leave as faces (=default) or triangulate
         // triangulate     false;
     }
 
-    movingWall_interpolated
+    walls_interpolated
     {
         type            patch;
-        patchName       movingWall;
+        patches         ( ".*Wall.*" );
         interpolate     true;
         // Optional: whether to leave as faces (=default) or triangulate
         // triangulate     false;
     }
 
-    movingNearWall_interpolated
+    nearWalls_interpolated
     {
         // Sample cell values off patch. Does not need to be the near-wall
         // cell, can be arbitrarily far away.
         type            patchInternalField;
-        patchName       movingWall;
+        patches         ( ".*Wall.*" );
         distance        0.0001;
         interpolate     true;
         // Optional: whether to leave as faces (=default) or triangulate

src/lagrangian/basic/particle/particleI.H

@@ -97,12 +97,19 @@ inline Foam::scalar Foam::particle::tetLambda
     scalar lambdaNumerator = (base - from) & n;
     scalar lambdaDenominator = (to - from) & n;
 
-    if (mag(lambdaDenominator) < SMALL)
+    // n carries the area of the tet faces, so the dot product with a
+    // delta-length has the units of volume.  Comparing the component of each
+    // delta-length in the direction of n times the face area to a fraction of
+    // the cell volume.
+
+    scalar tol = 1e3*SMALL*mesh_.cellVolumes()[cellI];
+
+    if (mag(lambdaDenominator) < tol)
     {
-        if (mag(lambdaNumerator) < SMALL)
+        if (mag(lambdaNumerator) < tol)
         {
             // Track starts on the face, and is potentially
-            // parallel to it.  +-SMALL/+-SMALL is not a good
+            // parallel to it.  +-tol/+-tol is not a good
             // comparison, return 0.0, in anticipation of tet
             // centre correction.
 
@@ -110,10 +117,11 @@ inline Foam::scalar Foam::particle::tetLambda
         }
         else
         {
-            if (mag((to - from)) < SMALL)
+            if (mag((to - from)) < tol/mag(n))
             {
-                // Zero length track, not along the face, face
+                // 'Zero' length track (compared to the tolerance, which is
-                // cannot be crossed.
+                // based on the cell volume, divided by the tet face area), not
+                // along the face, face cannot be crossed.
 
                 return GREAT;
             }
@@ -291,12 +299,14 @@ inline Foam::scalar Foam::particle::movingTetLambda
 
     }
 
-    if (mag(lambdaDenominator) < SMALL)
+    scalar tol = 1e3*SMALL*mesh_.cellVolumes()[cellI];
+
+    if (mag(lambdaDenominator) < tol)
     {
-        if (mag(lambdaNumerator) < SMALL)
+        if (mag(lambdaNumerator) < tol)
         {
             // Track starts on the face, and is potentially
-            // parallel to it.  +-SMALL)/+-SMALL is not a good
+            // parallel to it.  +-tol)/+-tol is not a good
             // comparison, return 0.0, in anticipation of tet
             // centre correction.
 
@@ -304,7 +314,7 @@ inline Foam::scalar Foam::particle::movingTetLambda
         }
         else
         {
-            if (mag((to - from)) < SMALL)
+            if (mag((to - from)) < tol/mag(n))
             {
                 // Zero length track, not along the face, face
                 // cannot be crossed.

src/lagrangian/basic/particle/particleTemplates.C

@@ -30,6 +30,7 @@ License
 #include "symmetryPolyPatch.H"
 #include "wallPolyPatch.H"
 #include "wedgePolyPatch.H"
+#include "meshTools.H"
 
 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
 
@@ -463,7 +464,7 @@ Foam::scalar Foam::particle::trackToFace
         //     << origId_ << " " << origProc_<< nl
         //     << "# face: " << tetFaceI_ << nl
         //     << "# tetPtI: " << tetPtI_ << nl
-        //     << "# tetBasePtI: " << mesh.tetBasePtIs()[tetFaceI_] << nl
+        //     << "# tetBasePtI: " << mesh_.tetBasePtIs()[tetFaceI_] << nl
         //     << "# tet.mag(): " << tet.mag() << nl
         //     << "# tet.quality(): " << tet.quality()
         //     << endl;

src/sampling/sampledSurface/sampledPatch/sampledPatch.C

@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2004-2010 OpenCFD Ltd.
+    \\  /    A nd           | Copyright (C) 2004-2011 OpenCFD Ltd.
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -46,15 +46,14 @@ Foam::sampledPatch::sampledPatch
 (
     const word& name,
     const polyMesh& mesh,
-    const word& patchName,
+    const wordReList& patchNames,
     const bool triangulate
 )
 :
     sampledSurface(name, mesh),
-    patchName_(patchName),
+    patchNames_(patchNames),
     triangulate_(triangulate),
-    needsUpdate_(true),
+    needsUpdate_(true)
-    patchFaceLabels_(0)
 {}
 
 
@@ -66,10 +65,9 @@ Foam::sampledPatch::sampledPatch
 )
 :
     sampledSurface(name, mesh, dict),
-    patchName_(dict.lookup("patchName")),
+    patchNames_(dict.lookup("patches")),
     triangulate_(dict.lookupOrDefault("triangulate", false)),
-    needsUpdate_(true),
+    needsUpdate_(true)
-    patchFaceLabels_(0)
 {}
 
 
@@ -81,6 +79,20 @@ Foam::sampledPatch::~sampledPatch()
 
 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 
+const Foam::labelList& Foam::sampledPatch::patchIDs() const
+{
+    if (patchIDs_.empty())
+    {
+        patchIDs_ = mesh().boundaryMesh().patchSet
+        (
+            patchNames_,
+            false
+        ).sortedToc();
+    }
+    return patchIDs_;
+}
+
+
 bool Foam::sampledPatch::needsUpdate() const
 {
     return needsUpdate_;
@@ -97,7 +109,10 @@ bool Foam::sampledPatch::expire()
 
     sampledSurface::clearGeom();
     MeshStorage::clear();
+    patchIDs_.clear();
+    patchIndex_.clear();
     patchFaceLabels_.clear();
+    patchStart_.clear();
 
     needsUpdate_ = true;
     return true;
@@ -111,21 +126,58 @@ bool Foam::sampledPatch::update()
         return false;
     }
 
-    const label patchI = mesh().boundaryMesh().findPatchID(patchName_);
+    label sz = 0;
-
+    forAll(patchIDs(), i)
-    if (patchI != -1)
     {
-        const polyPatch& p = mesh().boundaryMesh()[patchI];
+        label patchI = patchIDs()[i];
-        this->storedPoints() = p.localPoints();
+        const polyPatch& pp = mesh().boundaryMesh()[patchI];
-        this->storedFaces()  = p.localFaces();
 
-        // an identity map
+        if (isA<emptyPolyPatch>(pp))
-        patchFaceLabels_.setSize(faces().size());
-        forAll(patchFaceLabels_, i)
         {
-            patchFaceLabels_[i] = i;
+            FatalErrorIn("sampledPatch::update()")
+                << "Cannot sample an empty patch. Patch " << pp.name()
+                << exit(FatalError);
         }
 
+        sz += pp.size();
+    }
+
+    // For every face (or triangle) the originating patch and local face in the
+    // patch.
+    patchIndex_.setSize(sz);
+    patchFaceLabels_.setSize(sz);
+    patchStart_.setSize(patchIDs().size());
+    labelList meshFaceLabels(sz);
+
+    sz = 0;
+
+    forAll(patchIDs(), i)
+    {
+        label patchI = patchIDs()[i];
+
+        patchStart_[i] = sz;
+
+        const polyPatch& pp = mesh().boundaryMesh()[patchI];
+
+        forAll(pp, j)
+        {
+            patchIndex_[sz] = i;
+            patchFaceLabels_[sz] = j;
+            meshFaceLabels[sz] = pp.start()+j;
+            sz++;
+        }
+    }
+
+    indirectPrimitivePatch allPatches
+    (
+        IndirectList<face>(mesh().faces(), meshFaceLabels),
+        mesh().points()
+    );
+
+    this->storedPoints() = allPatches.localPoints();
+    this->storedFaces()  = allPatches.localFaces();
+
+
     // triangulate uses remapFaces()
     // - this is somewhat less efficient since it recopies the faces
     // that we just created, but we probably don't want to do this
@@ -134,7 +186,6 @@ bool Foam::sampledPatch::update()
     {
         MeshStorage::triangulate();
     }
-    }
 
     if (debug)
     {
@@ -148,10 +199,7 @@ bool Foam::sampledPatch::update()
 
 
 // remap action on triangulation
-void Foam::sampledPatch::remapFaces
+void Foam::sampledPatch::remapFaces(const labelUList& faceMap)
-(
-    const labelUList& faceMap
-)
 {
     // recalculate the cells cut
     if (&faceMap && faceMap.size())
@@ -161,9 +209,25 @@ void Foam::sampledPatch::remapFaces
         (
             UIndirectList<label>(patchFaceLabels_, faceMap)
         );
-    }
+        patchIndex_ = labelList
-}
+        (
+            UIndirectList<label>(patchIndex_, faceMap)
+        );
 
+        // Redo patchStart.
+        if (patchIndex_.size() > 0)
+        {
+            patchStart_[patchIndex_[0]] = 0;
+            for (label i = 1; i < patchIndex_.size(); i++)
+            {
+                if (patchIndex_[i] != patchIndex_[i-1])
+                {
+                    patchStart_[patchIndex_[i]] = i;
+                }
+            }
+        }
+    }
+}
 
 
 Foam::tmp<Foam::scalarField> Foam::sampledPatch::sample
@@ -257,7 +321,7 @@ Foam::tmp<Foam::tensorField> Foam::sampledPatch::interpolate
 void Foam::sampledPatch::print(Ostream& os) const
 {
     os  << "sampledPatch: " << name() << " :"
-        << "  patch:" << patchName()
+        << "  patches:" << patchNames()
         << "  faces:" << faces().size()
         << "  points:" << points().size();
 }

src/sampling/sampledSurface/sampledPatch/sampledPatch.H

@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2004-2010 OpenCFD Ltd.
+    \\  /    A nd           | Copyright (C) 2004-2011 OpenCFD Ltd.
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -25,7 +25,7 @@ Class
     Foam::sampledPatch
 
 Description
-    A sampledSurface on a patch. Non-triangulated by default.
+    A sampledSurface on patches. Non-triangulated by default.
 
 SourceFiles
     sampledPatch.C
@@ -57,8 +57,11 @@ class sampledPatch
 
     // Private data
 
-        //- Name of patch
+        //- Name of patches
-        const word patchName_;
+        const wordReList patchNames_;
+
+        //- Corresponding patchIDs
+        mutable labelList patchIDs_;
 
         //- Triangulated faces or keep faces as is
         bool triangulate_;
@@ -66,11 +69,16 @@ class sampledPatch
         //- Track if the surface needs an update
         mutable bool needsUpdate_;
 
-        //- Local patch face labels
+        //- For every face (or triangle) the originating patch
+        labelList patchIndex_;
+
+        //- For every face (or triangle) the index in the originating patch
         labelList patchFaceLabels_;
 
-    // Private Member Functions
+        //- Start indices (in patchFaceLabels_) of patches
+        labelList patchStart_;
 
+    // Private Member Functions
 
         //- sample field on faces
         template <class Type>
@@ -79,15 +87,34 @@ class sampledPatch
             const GeometricField<Type, fvPatchField, volMesh>& vField
         ) const;
 
-
         template <class Type>
         tmp<Field<Type> >
         interpolateField(const interpolation<Type>&) const;
 
-
         //- remap action on triangulation or cleanup
         virtual void remapFaces(const labelUList& faceMap);
 
+
+protected:
+
+        const wordReList& patchNames() const
+        {
+            return patchNames_;
+        }
+
+        const labelList& patchIDs() const;
+
+        const labelList& patchStart() const
+        {
+            return patchStart_;
+        }
+
+        const labelList& patchFaceLabels() const
+        {
+            return patchFaceLabels_;
+        }
+
+
 public:
 
     //- Runtime type information
@@ -101,7 +128,7 @@ public:
         (
             const word& name,
             const polyMesh& mesh,
-            const word& patchName,
+            const wordReList& patchNames,
             const bool triangulate = false
         );
 
@@ -133,21 +160,6 @@ public:
         virtual bool update();
 
 
-        const word patchName() const
-        {
-            return patchName_;
-        }
-
-        label patchIndex() const
-        {
-            return mesh().boundaryMesh().findPatchID(patchName_);
-        }
-
-        const labelList& patchFaceLabels() const
-        {
-            return patchFaceLabels_;
-        }
-
         //- Points of surface
         virtual const pointField& points() const
         {

src/sampling/sampledSurface/sampledPatch/sampledPatchTemplates.C

@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2004-2010 OpenCFD Ltd.
+    \\  /    A nd           | Copyright (C) 2004-2011 OpenCFD Ltd.
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -37,15 +37,11 @@ Foam::sampledPatch::sampleField
     // One value per face
     tmp<Field<Type> > tvalues(new Field<Type>(patchFaceLabels_.size()));
     Field<Type>& values = tvalues();
-
+    forAll(patchFaceLabels_, i)
-    if (patchIndex() != -1)
     {
-        const Field<Type>& bField = vField.boundaryField()[patchIndex()];
+        label patchI = patchIDs_[patchIndex_[i]];
-
+        const Field<Type>& bField = vField.boundaryField()[patchI];
-        forAll(patchFaceLabels_, elemI)
+        values[i] = bField[patchFaceLabels_[i]];
-        {
-            values[elemI] = bField[patchFaceLabels_[elemI]];
-        }
     }
 
     return tvalues;
@@ -63,15 +59,15 @@ Foam::sampledPatch::interpolateField
     tmp<Field<Type> > tvalues(new Field<Type>(points().size()));
     Field<Type>& values = tvalues();
 
-    if (patchIndex() != -1)
-    {
-        const polyPatch& patch = mesh().boundaryMesh()[patchIndex()];
     const labelList& own = mesh().faceOwner();
 
     boolList pointDone(points().size(), false);
 
     forAll(faces(), cutFaceI)
     {
+        label patchI = patchIDs_[patchIndex_[cutFaceI]];
+        const polyPatch& pp = mesh().boundaryMesh()[patchI];
+        label patchFaceI = patchFaceLabels()[cutFaceI];
         const face& f = faces()[cutFaceI];
 
         forAll(f, faceVertI)
@@ -80,7 +76,7 @@ Foam::sampledPatch::interpolateField
 
             if (!pointDone[pointI])
             {
-                    label faceI = patchFaceLabels()[cutFaceI] + patch.start();
+                label faceI = patchFaceI + pp.start();
                 label cellI = own[faceI];
 
                 values[pointI] = interpolator.interpolate
@@ -93,7 +89,6 @@ Foam::sampledPatch::interpolateField
             }
         }
     }
-    }
 
     return tvalues;
 }

src/sampling/sampledSurface/sampledPatchInternalField/sampledPatchInternalField.C

@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2004-2010 OpenCFD Ltd.
+    \\  /    A nd           | Copyright (C) 2004-2011 OpenCFD Ltd.
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -56,15 +56,27 @@ Foam::sampledPatchInternalField::sampledPatchInternalField
 )
 :
     sampledPatch(name, mesh, dict),
-    directMappedPatchBase
+    mappers_(patchIDs().size())
+{
+    const scalar distance = readScalar(dict.lookup("distance"));
+
+    forAll(patchIDs(), i)
+    {
+        label patchI = patchIDs()[i];
+        mappers_.set
         (
-        mesh.boundaryMesh()[sampledPatch::patchIndex()],
+            i,
+            new directMappedPatchBase
+            (
+                mesh.boundaryMesh()[patchI],
                 mesh.name(),                        // sampleRegion
                 directMappedPatchBase::NEARESTCELL, // sampleMode
                 word::null,                         // samplePatch
-        -readScalar(dict.lookup("distance"))
+                -distance                           // sample inside my domain
             )
-{}
+        );
+    }
+}
 
 
 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
@@ -167,7 +179,7 @@ Foam::tmp<Foam::tensorField> Foam::sampledPatchInternalField::interpolate
 void Foam::sampledPatchInternalField::print(Ostream& os) const
 {
     os  << "sampledPatchInternalField: " << name() << " :"
-        << "  patch:" << patchName()
+        << "  patches:" << patchNames()
         << "  faces:" << faces().size()
         << "  points:" << points().size();
 }

src/sampling/sampledSurface/sampledPatchInternalField/sampledPatchInternalField.H

@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2004-2010 OpenCFD Ltd.
+    \\  /    A nd           | Copyright (C) 2004-2011 OpenCFD Ltd.
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -54,9 +54,13 @@ namespace Foam
 
 class sampledPatchInternalField
 :
-    public sampledPatch,
+    public sampledPatch
-    public directMappedPatchBase
 {
+    // Private data
+
+        //- Mapping engines
+        PtrList<directMappedPatchBase> mappers_;
+
 
     // Private Member Functions
 

src/sampling/sampledSurface/sampledPatchInternalField/sampledPatchInternalFieldTemplates.C

@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2004-2010 OpenCFD Ltd.
+    \\  /    A nd           | Copyright (C) 2004-2011 OpenCFD Ltd.
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -36,20 +36,27 @@ Foam::sampledPatchInternalField::sampleField
     const GeometricField<Type, fvPatchField, volMesh>& vField
 ) const
 {
-    const mapDistribute& distMap = map();
-
     // One value per face
     tmp<Field<Type> > tvalues(new Field<Type>(patchFaceLabels().size()));
     Field<Type>& values = tvalues();
 
-    if (patchIndex() != -1)
+    forAll(patchStart(), i)
     {
+        // Get patchface wise data by sampling internal field
         Field<Type> interpVals = vField.internalField();
-        distMap.distribute(interpVals);
+        mappers_[i].map().distribute(interpVals);
 
-        forAll(patchFaceLabels(), elemI)
+        // Store at correct position in values
+        label end =
+        (
+            i < patchStart().size()-1
+          ? patchStart()[i+1]
+          : patchFaceLabels().size()
+        );
+
+        for (label triI = patchStart()[i]; triI < end; triI++)
         {
-            values[elemI] = interpVals[patchFaceLabels()[elemI]];
+            values[triI] = interpVals[patchFaceLabels()[triI]];
         }
     }
 
@@ -64,19 +71,24 @@ Foam::sampledPatchInternalField::interpolateField
     const interpolation<Type>& interpolator
 ) const
 {
-    // One value per vertex
+    label sz = 0;
+    forAll(patchIDs(), i)
+    {
+        sz += mesh().boundaryMesh()[patchIDs()[i]].size();
+    }
 
-    if (patchIndex() != -1)
+    Field<Type> allPatchVals(sz);
+    sz = 0;
+
+    forAll(patchIDs(), i)
     {
         // See directMappedFixedValueFvPatchField
-        const mapDistribute& distMap = map();
+        const mapDistribute& distMap = mappers_[i].map();
-
-        const polyPatch& pp = mesh().boundaryMesh()[patchIndex()];
 
         // Send back sample points to processor that holds the cell.
         // Mark cells with point::max so we know which ones we need
         // to interpolate (since expensive).
-        vectorField samples(samplePoints());
+        vectorField samples(mappers_[i].samplePoints());
         distMap.reverseDistribute(mesh().nCells(), point::max, samples);
 
         Field<Type> patchVals(mesh().nCells());
@@ -96,19 +108,36 @@ Foam::sampledPatchInternalField::interpolateField
         distMap.distribute(patchVals);
 
         // Now patchVals holds the interpolated data in patch face order.
-        // Interpolate to points. Note: points are patch.localPoints() so
+        // Collect.
-        // can use standard interpolation
+        SubList<Type>(allPatchVals, patchVals.size(), sz).assign(patchVals);
+        sz += patchVals.size();
+    }
 
-        return PrimitivePatchInterpolation<primitivePatch>
+    // Interpolate to points. Reconstruct the patch of all faces to aid
-        (
+    // interpolation.
-           pp
+
-        ).faceToPointInterpolate(patchVals);
+    labelList meshFaceLabels(allPatchVals.size());
-    }
+    sz = 0;
-    else
+    forAll(patchIDs(), i)
     {
-        return tmp<Field<Type> >(new Field<Type>(points().size()));
+        const polyPatch& pp = mesh().boundaryMesh()[patchIDs()[i]];
+        forAll(pp, i)
+        {
+            meshFaceLabels[sz++] = pp.start()+i;
         }
     }
 
+    indirectPrimitivePatch allPatches
+    (
+        IndirectList<face>(mesh().faces(), meshFaceLabels),
+        mesh().points()
+    );
+
+    return PrimitivePatchInterpolation<indirectPrimitivePatch>
+    (
+        allPatches
+    ).faceToPointInterpolate(allPatchVals);
+}
+
 
 // ************************************************************************* //

tutorials/incompressible/pisoFoam/les/pitzDaily/system/controlDict

@@ -124,7 +124,7 @@ functions
             nearWall
             {
                 type            patchInternalField;
-                patchName       lowerWall;
+                patches         ( lowerWall );
                 distance        1E-6;
                 interpolate     true;
                 triangulate     false;

tutorials/incompressible/pisoFoam/les/pitzDailyDirectMapped/system/controlDict

@@ -124,7 +124,7 @@ functions
             nearWall
             {
                 type            patchInternalField;
-                patchName       lowerWall;
+                patches         ( lowerWall );
                 distance        1E-6;
                 interpolate     true;
                 triangulate     false;

tutorials/mesh/snappyHexMesh/flange/Allclean

@@ -4,6 +4,7 @@ cd ${0%/*} || exit 1    # run from this directory
 # Source tutorial clean functions
 . $WM_PROJECT_DIR/bin/tools/CleanFunctions
 
+rm -rf 0 > /dev/null 2>&1
 rm -f ./*.obj > /dev/null 2>&1
 rm -rf constant/extendedFeatureEdgeMesh > /dev/null 2>&1
 rm -f constant/triSurface/flange.eMesh > /dev/null 2>&1

tutorials/mesh/snappyHexMesh/flange/Allrun

@@ -3,5 +3,5 @@
 . $WM_PROJECT_DIR/bin/tools/RunFunctions
 
 runApplication blockMesh
-surfaceFeatureExtract -includedAngle 135 -writeObj constant/triSurface/flange.stl flange
+runApplication surfaceFeatureExtract -includedAngle 135 -writeObj constant/triSurface/flange.stl flange
 runApplication snappyHexMesh -overwrite

tutorials/mesh/snappyHexMesh/flange/constant/polyMesh/blockMeshDict

@@ -1,8 +1,8 @@
 /*--------------------------------*- C++ -*----------------------------------*\
 | =========                 |                                                 |
 | \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
-|  \\    /   O peration     | Version:  1.6                                   |
+|  \\    /   O peration     | Version:  dev                                   |
-|   \\  /    A nd           | Web:      http://www.OpenFOAM.org               |
+|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
 |    \\/     M anipulation  |                                                 |
 \*---------------------------------------------------------------------------*/
 FoamFile
@@ -18,19 +18,19 @@ convertToMeters 1;
 
 vertices
 (
-    ( -0.04 -0.043 -0.04)
+    ( -0.0260093 -0.0275 -0.02375)
-    (  0.04 -0.043 -0.04)
+    (  0.0260093 -0.0275 -0.02375)
-    (  0.04  0.04 -0.04)
+    (  0.0260093  0.0225 -0.02375)
-    ( -0.04  0.04 -0.04)
+    ( -0.0260093  0.0225 -0.02375)
-    ( -0.04 -0.043 0.01)
+    ( -0.0260093 -0.0275 0.00225014)
-    (  0.04 -0.043 0.01)
+    (  0.0260093 -0.0275 0.00225014)
-    (  0.04  0.04 0.01)
+    (  0.0260093  0.0225 0.00225014)
-    ( -0.04  0.04 0.01)
+    ( -0.0260093  0.0225 0.00225014)
 );
 
 blocks
 (
-    hex (0 1 2 3 4 5 6 7) (20 20 12) simpleGrading (1 1 1)
+    hex (0 1 2 3 4 5 6 7) (15 15 9) simpleGrading (1 1 1)
 );
 
 edges

tutorials/mesh/snappyHexMesh/flange/constant/transportProperties

@@ -1,21 +0,0 @@
-/*--------------------------------*- C++ -*----------------------------------*\
-| =========                 |                                                 |
-| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
-|  \\    /   O peration     | Version:  1.6                                   |
-|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
-|    \\/     M anipulation  |                                                 |
-\*---------------------------------------------------------------------------*/
-FoamFile
-{
-    version     2.0;
-    format      ascii;
-    class       dictionary;
-    location    "constant";
-    object      transportProperties;
-}
-// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
-
-nu              nu [ 0 2 -1 0 0 0 0 ] 0.01;
-
-
-// ************************************************************************* //

tutorials/mesh/snappyHexMesh/flange/system/controlDict

@@ -1,8 +1,8 @@
 /*--------------------------------*- C++ -*----------------------------------*\
 | =========                 |                                                 |
 | \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
-|  \\    /   O peration     | Version:  1.6                                   |
+|  \\    /   O peration     | Version:  dev                                   |
-|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
 |    \\/     M anipulation  |                                                 |
 \*---------------------------------------------------------------------------*/
 FoamFile

tutorials/mesh/snappyHexMesh/flange/system/decomposeParDict

@@ -1,8 +1,8 @@
 /*--------------------------------*- C++ -*----------------------------------*\
 | =========                 |                                                 |
 | \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
-|  \\    /   O peration     | Version:  1.6                                   |
+|  \\    /   O peration     | Version:  dev                                   |
-|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
 |    \\/     M anipulation  |                                                 |
 \*---------------------------------------------------------------------------*/
 FoamFile

tutorials/mesh/snappyHexMesh/flange/system/fvSchemes

@@ -1,8 +1,8 @@
 /*--------------------------------*- C++ -*----------------------------------*\
 | =========                 |                                                 |
 | \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
-|  \\    /   O peration     | Version:  1.6                                   |
+|  \\    /   O peration     | Version:  dev                                   |
-|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
 |    \\/     M anipulation  |                                                 |
 \*---------------------------------------------------------------------------*/
 FoamFile
@@ -15,46 +15,16 @@ FoamFile
 }
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
-ddtSchemes
-{
-    default         Euler;
-}
-
 gradSchemes
 {
-    default         Gauss linear;
-    grad(p)         Gauss linear;
 }
 
 divSchemes
 {
-    default         none;
-    div(phi,U)      Gauss linear;
 }
 
 laplacianSchemes
 {
-    default         none;
-    laplacian(nu,U) Gauss linear corrected;
-    laplacian((1|A(U)),p) Gauss linear corrected;
 }
 
-interpolationSchemes
-{
-    default         linear;
-    interpolate(HbyA) linear;
-}
-
-snGradSchemes
-{
-    default         corrected;
-}
-
-fluxRequired
-{
-    default         no;
-    p               ;
-}
-
-
 // ************************************************************************* //

tutorials/mesh/snappyHexMesh/flange/system/fvSolution

@@ -1,8 +1,8 @@
 /*--------------------------------*- C++ -*----------------------------------*\
 | =========                 |                                                 |
 | \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
-|  \\    /   O peration     | Version:  1.6                                   |
+|  \\    /   O peration     | Version:  dev                                   |
-|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
 |    \\/     M anipulation  |                                                 |
 \*---------------------------------------------------------------------------*/
 FoamFile
@@ -15,32 +15,4 @@ FoamFile
 }
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
-solvers
-{
-    p
-    {
-        solver          PCG;
-        preconditioner  DIC;
-        tolerance       1e-06;
-        relTol          0;
-    }
-
-    U
-    {
-        solver          PBiCG;
-        preconditioner  DILU;
-        tolerance       1e-05;
-        relTol          0;
-    }
-}
-
-PISO
-{
-    nCorrectors     2;
-    nNonOrthogonalCorrectors 0;
-    pRefCell        0;
-    pRefValue       0;
-}
-
-
 // ************************************************************************* //

tutorials/mesh/snappyHexMesh/flange/system/snappyHexMeshDict

@@ -2,7 +2,7 @@
 | =========                 |                                                 |
 | \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
 |  \\    /   O peration     | Version:  dev                                   |
-|   \\  /    A nd           | Web:      http://www.OpenFOAM.org               |
+|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
 |    \\/     M anipulation  |                                                 |
 \*---------------------------------------------------------------------------*/
 FoamFile
@@ -144,8 +144,8 @@ castellatedMeshControls
     // section reachable from the locationInMesh is kept.
     // NOTE: This point should never be on a face, always inside a cell, even
     // after refinement.
-    locationInMesh (-0.033 -0.033 0.0033);
+    // This is an outside point locationInMesh (-0.033 -0.033 0.0033);
-
+    locationInMesh (-9.23149e-05 -0.0025 -0.0025); // Inside point
 
     // Whether any faceZones (as specified in the refinementSurfaces)
     // are only on the boundary of corresponding cellZones or also allow

tutorials/multiphase/interDyMFoam/ras/sloshingTank2D/system/controlDict

@@ -89,7 +89,7 @@ functions
             walls
             {
                 type        patch;
-                patchName   walls;
+                patches     (walls);
                 triangulate false;
             }
         );

tutorials/multiphase/interDyMFoam/ras/sloshingTank2D3DoF/system/controlDict

@@ -88,7 +88,7 @@ functions
             walls
             {
                 type        patch;
-                patchName   walls;
+                patches     (walls);
                 triangulate false;
             }
         );