Added a utility for generating boundary layers

utilities/generateBoundaryLayers/Make/files

+generateBoundaryLayers.C
+
+EXE = $(FOAM_USER_APPBIN)/generateBoundaryLayers

utilities/generateBoundaryLayers/Make/options

+EXE_INC = \
+    -I../../meshLibrary/lnInclude \
+    -I$(LIB_SRC)/triSurface/lnInclude \
+    -I$(LIB_SRC)/meshTools/lnInclude
+
+EXE_LIBS = \
+    -lmeshTools \
+    -ltriSurface \
+    -L$(FOAM_USER_LIBBIN) \
+    -lmeshLibrary

utilities/generateBoundaryLayers/generateBoundaryLayers.C

+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | cfMesh: A library for mesh generation
+   \\    /   O peration     |
+    \\  /    A nd           | Author: Franjo Juretic (franjo.juretic@c-fields.com)
+     \\/     M anipulation  | Copyright (C) Creative Fields, Ltd.
+-------------------------------------------------------------------------------
+License
+    This file is part of cfMesh.
+
+    cfMesh 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.
+
+    cfMesh 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 cfMesh.  If not, see <http://www.gnu.org/licenses/>.
+
+Description
+    Performs point relocations in the mesh (smoothing) in order to
+    improve quality measures. It does not make the mesh invalied.
+
+\*---------------------------------------------------------------------------*/
+
+#include "argList.H"
+#include "Time.H"
+#include "polyMeshGenModifier.H"
+#include "meshOptimizer.H"
+
+using namespace Foam;
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+int main(int argc, char *argv[])
+{
+    argList::validArgs.clear();
+
+    argList::validOptions.insert("nLoops", "int");
+    argList::validOptions.insert("nIterations", "int");
+    argList::validOptions.insert("nSurfaceIterations", "int");
+    argList::validOptions.insert("qualityThreshold", "scalar");
+    argList::validOptions.insert("constrainedCellsSet", "word");
+
+#   include "setRootCase.H"
+#   include "createTime.H"
+
+    //- read the settings
+    label nIterations(50);
+    label nLoops(10);
+    label nSurfaceIterations(2);
+    scalar qualityThreshold(0.1);
+
+    if( args.options().found("nLoops") )
+    {
+        nLoops = readLabel(IStringStream(args.options()["nLoops"])());
+    }
+    else
+    {
+        Info << "Default number of loops is 10" << endl;
+    }
+
+    if( args.options().found("nIterations") )
+    {
+        nIterations =
+            readLabel(IStringStream(args.options()["nIterations"])());
+    }
+    else
+    {
+        Info << "Default number of iterations is 50" << endl;
+    }
+
+    if( args.options().found("nSurfaceIterations") )
+    {
+        nSurfaceIterations =
+            readLabel(IStringStream(args.options()["nSurfaceIterations"])());
+    }
+    else
+    {
+        Info << "Default number of surface iterations is 2" << endl;
+    }
+
+    if( args.options().found("qualityThreshold") )
+    {
+        qualityThreshold =
+            readScalar(IStringStream(args.options()["qualityThreshold"])());
+    }
+    else
+    {
+        Info << "Using default quality threshold 0.1" << endl;
+    }
+
+    word constrainedCellSet;
+
+    if( args.options().found("constrainedCellsSet") )
+    {
+        constrainedCellSet = args.options()["constrainedCellsSet"];
+    }
+    else
+    {
+        Info << "No constraints applied on the smoothing procedure" << endl;
+    }
+
+    //- load the mesh from disk
+    polyMeshGen pmg(runTime);
+    pmg.read();
+
+    //- construct the smoother
+    meshOptimizer mOpt(pmg);
+
+    if( !constrainedCellSet.empty() )
+    {
+        //- lock cells in constrainedCellSet
+        mOpt.lockCellsInSubset(constrainedCellSet);
+
+        //- find boundary faces which shall be locked
+        labelLongList lockedBndFaces, selectedCells;
+
+        const label sId = pmg.cellSubsetIndex(constrainedCellSet);
+        pmg.cellsInSubset(sId, selectedCells);
+
+        boolList activeCell(pmg.cells().size(), false);
+        forAll(selectedCells, i)
+            activeCell[selectedCells[i]] = true;
+    }
+
+    //- clear geometry information before volume smoothing
+    pmg.clearAddressingData();
+
+    //- perform optimisation using the laplace smoother and
+    mOpt.optimizeMeshFV
+    (
+        nLoops,
+        nLoops,
+        nIterations,
+        nSurfaceIterations
+    );
+
+    //- perform optimisation of worst quality faces
+    mOpt.optimizeMeshFVBestQuality(nLoops, qualityThreshold);
+
+    //- check the mesh again and untangl bad regions if any of them exist
+    mOpt.untangleMeshFV(nLoops, nIterations, nSurfaceIterations);
+
+    Info << "Writing mesh" << endl;
+    pmg.write();
+
+    Info << "End\n" << endl;
+    return 0;
+}
+
+// ************************************************************************* //