Merge branch 'master' of ssh://dm/home/dm4/OpenFOAM/OpenFOAM-dev

applications/utilities/mesh/generation/snappyHexMesh/snappyHexMeshDict

@@ -255,6 +255,8 @@ addLayersControls
     relativeSizes true;
 
     // Per final patch (so not geometry!) the layer information
+    // Note: This behaviour changed after 21x. Any non-mentioned patches
+    //       now slide unless nSurfaceLayers is explicitly mentioned to be 0.
     layers
     {
         sphere.stl_firstSolid
@@ -265,6 +267,17 @@ addLayersControls
         maxY
         {
             nSurfaceLayers 1;
+            // Per patch layer data
+            expansionRatio 1.3;
+            finalLayerThickness 0.3;
+            minThickness 0.1;
+        }
+
+        // Disable any mesh shrinking and layer addition on any point of
+        // a patch by setting nSurfaceLayers to 0
+        frozenPatches
+        {
+            nSurfaceLayers 0;
         }
     }
 
@@ -281,6 +294,7 @@ addLayersControls
 
     //- Minimum thickness of cell layer. If for any reason layer
     //  cannot be above minThickness do not add layer.
+    //  Relative to undistorted size of cell outside layer.
     //  See relativeSizes parameter.
     minThickness 0.25;
 
@@ -417,7 +431,7 @@ meshQualityControls
 // Flags for optional output
 // 0 : only write final meshes
 // 1 : write intermediate meshes
-// 2 : write volScalarFields with cellLevel and cell centres for postprocessing
+// 2 : write volScalarField with cellLevel for postprocessing
 // 4 : write current intersections as .obj files
 debug 0;
 

src/finiteVolume/fields/fvPatchFields/derived/codedMixed/codedMixedFvPatchField.H

@@ -38,7 +38,11 @@ Description
     myPatch
     {
         type            codedMixed;
-        value           uniform 0;
+
+        refValue        uniform (0 0 0);
+        refGradient     uniform (0 0 0);
+        valueFraction   uniform 1;
+
         redirectType    rampedMixed;   // name of generated BC
 
         code

src/finiteVolume/fields/fvPatchFields/derived/waveSurfacePressure/waveSurfacePressureFvPatchScalarField.C

@@ -38,15 +38,15 @@ License
 namespace Foam
 {
     template<>
-    const char* Foam::NamedEnum
+    const char* NamedEnum
     <
-        Foam::waveSurfacePressureFvPatchScalarField::timeSchemeType,
+        waveSurfacePressureFvPatchScalarField::timeSchemeType,
         3
     >::names[] =
     {
-        fv::EulerDdtScheme<scalar>::typeName.c_str(),
-        fv::CrankNicholsonDdtScheme<scalar>::typeName.c_str(),
-        fv::backwardDdtScheme<scalar>::typeName.c_str()
+        fv::EulerDdtScheme<scalar>::typeName_(),
+        fv::CrankNicholsonDdtScheme<scalar>::typeName_(),
+        fv::backwardDdtScheme<scalar>::typeName_()
     };
 }
 

src/mesh/autoMesh/Make/options

@@ -5,6 +5,7 @@ EXE_INC = \
     -I$(LIB_SRC)/lagrangian/basic/lnInclude \
     -I$(LIB_SRC)/meshTools/lnInclude \
     -I$(LIB_SRC)/edgeMesh/lnInclude \
+    -I$(LIB_SRC)/surfMesh/lnInclude \
     -I$(LIB_SRC)/triSurface/lnInclude
 
 LIB_LIBS = \
@@ -13,5 +14,6 @@ LIB_LIBS = \
     -llagrangian \
     -lmeshTools \
     -ledgeMesh \
+    -lsurfMesh \
     -ltriSurface \
     -ldistributed

src/regionModels/pyrolysisModels/Make/options

@@ -19,8 +19,7 @@ LIB_LIBS = \
     -lmeshTools \
     -lchemistryModel \
     -lspecie \
-    -lspecie \
-    -lsolid \
+    -lsolidSpecie \
     -lfluidThermophysicalModels \
     -lsolidChemistryModel \
     -lcompressibleTurbulenceModel \

src/regionModels/thermoBaffleModels/Make/options

@@ -17,4 +17,4 @@ LIB_LIBS = \
     -lfiniteVolume \
     -lmeshTools \
     -lOpenFOAM \
-    -lsolid
+    -lsolidSpecie

src/surfMesh/Make/files

@@ -26,6 +26,7 @@ $(surfaceFormats)/ftr/FTRsurfaceFormatRunTime.C
 $(surfaceFormats)/gts/GTSsurfaceFormatRunTime.C
 $(surfaceFormats)/nas/NASsurfaceFormatRunTime.C
 $(surfaceFormats)/obj/OBJsurfaceFormatRunTime.C
+$(surfaceFormats)/obj/OBJstream.C
 $(surfaceFormats)/off/OFFsurfaceFormatRunTime.C
 $(surfaceFormats)/ofs/OFSsurfaceFormatCore.C
 $(surfaceFormats)/ofs/OFSsurfaceFormatRunTime.C

src/surfMesh/surfaceFormats/obj/OBJstream.C

@@ -0,0 +1,323 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2012 OpenFOAM Foundation
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software: you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by
+    the Free Software Foundation, either version 3 of the License, or
+    (at your option) any later version.
+
+    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
+    for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
+
+\*---------------------------------------------------------------------------*/
+
+#include "OBJstream.H"
+//#include "token.H"
+#include "primitivePatch.H"
+
+// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
+
+defineTypeNameAndDebug(Foam::OBJstream, 0);
+
+
+// * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
+
+void Foam::OBJstream::writeAndCheck(const char c)
+{
+    if (c == '\n')
+    {
+        startOfLine_ = true;
+    }
+    else if (startOfLine_)
+    {
+        startOfLine_ = false;
+        if (c == 'v')
+        {
+            nVertices_++;
+        }
+    }
+    OFstream::write(c);
+}
+
+
+// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
+
+Foam::OBJstream::OBJstream
+(
+    const fileName& pathname,
+    streamFormat format,
+    versionNumber version,
+    compressionType compression
+)
+:
+    OFstream(pathname, format, version, compression),
+    startOfLine_(true),
+    nVertices_(0)
+{}
+
+
+// * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
+
+Foam::OBJstream::~OBJstream()
+{}
+
+
+// * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
+
+Foam::Ostream& Foam::OBJstream::write(const char c)
+{
+    writeAndCheck(c);
+    return *this;
+}
+
+
+Foam::Ostream& Foam::OBJstream::write(const char* str)
+{
+    for (const char* p = str; *p != '\0'; ++p)
+    {
+        writeAndCheck(*p);
+    }
+    return *this;
+}
+
+
+Foam::Ostream& Foam::OBJstream::write(const word& str)
+{
+    write(str.c_str());
+    return *this;
+}
+
+
+Foam::Ostream& Foam::OBJstream::write(const string& str)
+{
+    OFstream::write(token::BEGIN_STRING);
+
+    register int backslash = 0;
+    for (string::const_iterator iter = str.begin(); iter != str.end(); ++iter)
+    {
+        register char c = *iter;
+
+        if (c == '\\')
+        {
+            backslash++;
+            // suppress output until we know if other characters follow
+            continue;
+        }
+        else if (c == token::NL)
+        {
+            lineNumber_++;
+            backslash++;    // backslash escape for newline
+        }
+        else if (c == token::END_STRING)
+        {
+            backslash++;    // backslash escape for quote
+        }
+
+        // output pending backslashes
+        while (backslash)
+        {
+            OFstream::write('\\');
+            backslash--;
+        }
+
+        writeAndCheck(c);
+    }
+
+    // silently drop any trailing backslashes
+    // they would otherwise appear like an escaped end-quote
+
+    OFstream::write(token::END_STRING);
+    return *this;
+}
+
+
+Foam::Ostream& Foam::OBJstream::writeQuoted
+(
+    const std::string& str,
+    const bool quoted
+)
+{
+    if (quoted)
+    {
+        OFstream::write(token::BEGIN_STRING);
+
+        register int backslash = 0;
+        for
+        (
+            string::const_iterator iter = str.begin();
+            iter != str.end();
+            ++iter
+        )
+        {
+            register char c = *iter;
+
+            if (c == '\\')
+            {
+                backslash++;
+                // suppress output until we know if other characters follow
+                continue;
+            }
+            else if (c == token::NL)
+            {
+                lineNumber_++;
+                backslash++;    // backslash escape for newline
+            }
+            else if (c == token::END_STRING)
+            {
+                backslash++;    // backslash escape for quote
+            }
+
+            // output pending backslashes
+            while (backslash)
+            {
+                OFstream::write('\\');
+                backslash--;
+            }
+
+            writeAndCheck(c);
+        }
+
+        // silently drop any trailing backslashes
+        // they would otherwise appear like an escaped end-quote
+        OFstream::write(token::END_STRING);
+    }
+    else
+    {
+        // output unquoted string, only advance line number on newline
+        write(str.c_str());
+    }
+
+    return *this;
+}
+
+
+Foam::Ostream& Foam::OBJstream::write(const point& pt)
+{
+    write("v ") << pt.x() << ' ' << pt.y() << ' ' << pt.z()
+        << nl;
+    return *this;
+}
+
+
+Foam::Ostream& Foam::OBJstream::write(const point& pt, const vector& n)
+{
+    write(pt);
+    OFstream::write("vn ") << n.x() << ' ' << n.y()
+        << ' ' << n.z() << nl;
+    return *this;
+}
+
+
+Foam::Ostream& Foam::OBJstream::write(const edge& e, const pointField& points)
+{
+    write(points[e[0]]);
+    write(points[e[1]]);
+    write("l ") << nVertices_-1 << ' ' << nVertices_ << nl;
+    return *this;
+}
+
+
+Foam::Ostream& Foam::OBJstream::write(const linePointRef& ln)
+{
+    write(ln.start());
+    write(ln.end());
+    write("l ") << nVertices_-1 << ' ' << nVertices_ << nl;
+    return *this;
+}
+
+
+Foam::Ostream& Foam::OBJstream::write
+(
+    const face& f,
+    const pointField& points,
+    const bool lines
+)
+{
+    label start = nVertices_;
+    forAll(f, i)
+    {
+        write(points[f[i]]);
+    }
+    if (lines)
+    {
+        write('l');
+        forAll(f, i)
+        {
+            write(' ') << start+1+i;
+        }
+        write(' ') << start+1 << '\n';
+    }
+    else
+    {
+        write('f');
+        forAll(f, i)
+        {
+            write(' ') << start+1+i;
+        }
+        write('\n');
+    }
+    return *this;
+}
+
+
+Foam::Ostream& Foam::OBJstream::write
+(
+    const faceList& fcs,
+    const pointField& points,
+    const bool lines
+)
+{
+    SubList<face> allFcs(fcs, fcs.size());
+
+    primitivePatch pp(allFcs, points);
+
+    const pointField& localPoints = pp.localPoints();
+    const faceList& localFaces = pp.localFaces();
+
+    label start = nVertices_;
+
+    forAll(localPoints, i)
+    {
+        write(localPoints[i]);
+    }
+
+    if (lines)
+    {
+        const edgeList& edges = pp.edges();
+        forAll(edges, edgeI)
+        {
+            const edge& e = edges[edgeI];
+
+            write("l ") << start+e[0]+1 << ' ' << start+e[1]+1 << nl;
+        }
+    }
+    else
+    {
+        forAll(localFaces, faceI)
+        {
+            const face& f = localFaces[faceI];
+            write('f');
+            forAll(f, i)
+            {
+                write(' ') << start+f[i]+1;
+            }
+            write('\n');
+        }
+    }
+    return *this;
+}
+
+
+// ************************************************************************* //

src/surfMesh/surfaceFormats/obj/OBJstream.H

@@ -0,0 +1,164 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2012 OpenFOAM Foundation
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software: you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by
+    the Free Software Foundation, either version 3 of the License, or
+    (at your option) any later version.
+
+    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
+    for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
+
+Class
+    Foam::OBJstream
+
+Description
+    OFstream which keeps track of vertices
+
+SourceFiles
+    OBJstream.C
+
+\*---------------------------------------------------------------------------*/
+
+#ifndef OBJstream_H
+#define OBJstream_H
+
+#include "OFstream.H"
+#include "point.H"
+#include "edge.H"
+#include "face.H"
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+namespace Foam
+{
+
+
+/*---------------------------------------------------------------------------*\
+                          Class OBJstream Declaration
+\*---------------------------------------------------------------------------*/
+
+class OBJstream
+:
+    public OFstream
+{
+    // Private data
+
+        bool startOfLine_;
+
+        label nVertices_;
+
+
+    // Private Member Functions
+
+        void writeAndCheck(const char);
+
+public:
+
+    // Declare name of the class and its debug switch
+    ClassName("OBJstream");
+
+
+    // Constructors
+
+        //- Construct from pathname
+        OBJstream
+        (
+            const fileName& pathname,
+            streamFormat format=ASCII,
+            versionNumber version=currentVersion,
+            compressionType compression=UNCOMPRESSED
+        );
+
+
+    //- Destructor
+    ~OBJstream();
+
+
+    // Member functions
+
+        // Access
+
+            //- Return the name of the stream
+            label nVertices() const
+            {
+                return nVertices_;
+            }
+
+
+        // Ostream implementation
+
+            //- Write character
+            virtual Ostream& write(const char);
+
+            //- Write character string
+            virtual Ostream& write(const char*);
+
+            //- Write word
+            virtual Ostream& write(const word&);
+
+            virtual Ostream& write(const string&);
+
+            //- Write std::string surrounded by quotes.
+            //  Optional write without quotes.
+            virtual Ostream& writeQuoted
+            (
+                const std::string&,
+                const bool quoted=true
+            );
+
+
+        // Direct write functionality
+
+            //- Write point
+            Ostream& write(const point&);
+
+            //- Write point and vector normal ('vn')
+            Ostream& write(const point&, const vector&);
+
+            //- Write edge as points with line
+            Ostream& write(const edge&, const pointField&);
+
+            //- Write line
+            Ostream& write(const linePointRef&);
+
+            //- Write face as points with lines or filled polygon
+            Ostream& write
+            (
+                const face&,
+                const pointField&,
+                const bool lines = true
+            );
+
+            //- Write patch as points and faces with lines or filled polygons
+            Ostream& write
+            (
+                const faceList&,
+                const pointField&,
+                const bool lines = true
+            );
+
+};
+
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+} // End namespace Foam
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#endif
+
+// ************************************************************************* //

src/turbulenceModels/compressible/RAS/derivedFvPatchFields/wallFunctions/epsilonWallFunctions/epsilonWallFunction/epsilonWallFunctionFvPatchScalarField.C

@@ -77,10 +77,7 @@ epsilonWallFunctionFvPatchScalarField::epsilonWallFunctionFvPatchScalarField
     GName_("RASModel::G"),
     Cmu_(0.09),
     kappa_(0.41),
-    E_(9.8),
-    yPlusLam_(mutWallFunctionFvPatchScalarField::yPlusLam(kappa_, E_))
-
-
+    E_(9.8)
 {
     checkType();
 }
@@ -98,9 +95,7 @@ epsilonWallFunctionFvPatchScalarField::epsilonWallFunctionFvPatchScalarField
     GName_(ptf.GName_),
     Cmu_(ptf.Cmu_),
     kappa_(ptf.kappa_),
-    E_(ptf.E_),
-    yPlusLam_(ptf.yPlusLam_)
-
+    E_(ptf.E_)
 {
     checkType();
 }
@@ -117,9 +112,7 @@ epsilonWallFunctionFvPatchScalarField::epsilonWallFunctionFvPatchScalarField
     GName_(dict.lookupOrDefault<word>("G", "RASModel::G")),
     Cmu_(dict.lookupOrDefault<scalar>("Cmu", 0.09)),
     kappa_(dict.lookupOrDefault<scalar>("kappa", 0.41)),
-    E_(dict.lookupOrDefault<scalar>("E", 9.8)),
-    yPlusLam_(mutWallFunctionFvPatchScalarField::yPlusLam(kappa_, E_))
-
+    E_(dict.lookupOrDefault<scalar>("E", 9.8))
 {
     checkType();
 }
@@ -134,8 +127,7 @@ epsilonWallFunctionFvPatchScalarField::epsilonWallFunctionFvPatchScalarField
     GName_(ewfpsf.GName_),
     Cmu_(ewfpsf.Cmu_),
     kappa_(ewfpsf.kappa_),
-    E_(ewfpsf.E_),
-    yPlusLam_(ewfpsf.yPlusLam_)
+    E_(ewfpsf.E_)
 {
     checkType();
 }
@@ -151,8 +143,7 @@ epsilonWallFunctionFvPatchScalarField::epsilonWallFunctionFvPatchScalarField
     GName_(ewfpsf.GName_),
     Cmu_(ewfpsf.Cmu_),
     kappa_(ewfpsf.kappa_),
-    E_(ewfpsf.E_),
-    yPlusLam_(ewfpsf.yPlusLam_)
+    E_(ewfpsf.E_)
 {
     checkType();
 }

src/turbulenceModels/compressible/RAS/derivedFvPatchFields/wallFunctions/epsilonWallFunctions/epsilonWallFunction/epsilonWallFunctionFvPatchScalarField.H

@@ -107,9 +107,6 @@ protected:
         //- E coefficient
         scalar E_;
 
-        //- Y+ at the edge of the laminar sublayer
-        scalar yPlusLam_;
-
 
     // Protected Member Functions
 

src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/wallFunctions/epsilonWallFunctions/epsilonWallFunction/epsilonWallFunctionFvPatchScalarField.C

@@ -78,8 +78,7 @@ epsilonWallFunctionFvPatchScalarField::epsilonWallFunctionFvPatchScalarField
     GName_("RASModel::G"),
     Cmu_(0.09),
     kappa_(0.41),
-    E_(9.8),
-    yPlusLam_(nutkWallFunctionFvPatchScalarField::yPlusLam(kappa_, E_))
+    E_(9.8)
 {
     checkType();
 }
@@ -97,8 +96,7 @@ epsilonWallFunctionFvPatchScalarField::epsilonWallFunctionFvPatchScalarField
     GName_(ptf.GName_),
     Cmu_(ptf.Cmu_),
     kappa_(ptf.kappa_),
-    E_(ptf.E_),
-    yPlusLam_(ptf.yPlusLam_)
+    E_(ptf.E_)
 {
     checkType();
 }
@@ -115,8 +113,7 @@ epsilonWallFunctionFvPatchScalarField::epsilonWallFunctionFvPatchScalarField
     GName_(dict.lookupOrDefault<word>("G", "RASModel::G")),
     Cmu_(dict.lookupOrDefault<scalar>("Cmu", 0.09)),
     kappa_(dict.lookupOrDefault<scalar>("kappa", 0.41)),
-    E_(dict.lookupOrDefault<scalar>("E", 9.8)),
-    yPlusLam_(nutkWallFunctionFvPatchScalarField::yPlusLam(kappa_, E_))
+    E_(dict.lookupOrDefault<scalar>("E", 9.8))
 {
     checkType();
 }
@@ -131,8 +128,7 @@ epsilonWallFunctionFvPatchScalarField::epsilonWallFunctionFvPatchScalarField
     GName_(ewfpsf.GName_),
     Cmu_(ewfpsf.Cmu_),
     kappa_(ewfpsf.kappa_),
-    E_(ewfpsf.E_),
-    yPlusLam_(ewfpsf.yPlusLam_)
+    E_(ewfpsf.E_)
 {
     checkType();
 }
@@ -148,8 +144,7 @@ epsilonWallFunctionFvPatchScalarField::epsilonWallFunctionFvPatchScalarField
     GName_(ewfpsf.GName_),
     Cmu_(ewfpsf.Cmu_),
     kappa_(ewfpsf.kappa_),
-    E_(ewfpsf.E_),
-    yPlusLam_(ewfpsf.yPlusLam_)
+    E_(ewfpsf.E_)
 {
     checkType();
 }

src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/wallFunctions/epsilonWallFunctions/epsilonWallFunction/epsilonWallFunctionFvPatchScalarField.H

@@ -107,9 +107,6 @@ protected:
         //- E coefficient
         scalar E_;
 
-        //- Y+ at the edge of the laminar sublayer
-        scalar yPlusLam_;
-
 
     // Protected Member Functions