[dhruv]

Hi,

I am using OF 2.0 and sHM. I have a rectangular channel with some internal thickness. I want to have boundary layers on the inside walls of this channel. When I do the meshing with snappy using the featureEdge, the boundary layers are not extended through the entire length of the channel, especially near the corners of the geometry. I am attaching pictures of the generated Mesh. Can anyone tell me how to make a consistent layer? I am also attaching the sHMDict file.

*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM Extend Project: Open Source CFD |
| \\ / O peration | Version: 1.6-ext |
| \\ / A nd | Web: www.extend-project.de |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object snappyHexMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

castellatedMesh true;
snap true;
addLayers true;


// ************************************************** **************************
geometry
{
geo.stl
{
type triSurfaceMesh;
name geo;
}
refBox
{
type searchableBox;
min (-20 0 -20);
max (855 850 850);
}
}
// ************************************************** **************************
castellatedMeshControls
{
maxLocalCells 5000000;
maxGlobalCells 10000000;
minRefinementCells 5;
//maxLoadUnbalance 0.10;
nCellsBetweenLevels 3;

features
(

{
file "geo.eMesh";
level 2;
}

);

refinementRegions
{
refBox
{
mode inside;
levels ((0.1 1));
}
}
refinementSurfaces
{
geo
{
level (1 2);
}
};
resolveFeatureAngle 150;

locationInMesh (50 50 50);
allowFreeStandingZoneFaces true;
}

// ************************************************** **************************
snapControls
{
nSmoothPatch 5;
tolerance 4.0;
nSolveIter 30;
nRelaxIter 5;
nFeatureSnapIter 10;
}

// ************************************************** **************************
addLayersControls
{
layers
{
"geo_walls"
{
nSurfaceLayers 10;
}


}
relativeSizes true;
expansionRatio 1.05;
finalLayerThickness 0.1;
minThickness 0.1;
nGrow 1;
featureAngle 150;
nRelaxIter 5;
nSmoothSurfaceNormals 1;
nSmoothNormals 3;
nSmoothThickness 10;
maxFaceThicknessRatio 0.5;
maxThicknessToMedialRatio 0.3;
minMedianAxisAngle 130;
nBufferCellsNoExtrude 0;
nLayerIter 50;
nRelaxedIter 20;
}

// ************************************************** **************************
meshQualityControls
{
maxNonOrtho 65;
maxBoundarySkewness 20;
maxInternalSkewness 4;
maxConcave 80;
minFlatness 0.5;
minVol 1e-13;
minTetQuality 1e-30;
minArea -1;
minTwist 0.05;
minDeterminant 0.001;
minFaceWeight 0.05;
minVolRatio 0.01;
minTriangleTwist -1;
nSmoothScale 4;
errorReduction 0.75;
relaxed
{
}
}

debug 0;

mergeTolerance 1e-6;

// ************************************************** *********************** //