[fra_mat]

Hi all,

I am a student quite new to OpenFoam. I have read the User's guide and studies some tutorials, so now I am starting with a very simple case based on the WindAroundBuildings tutorial. The mesh is very coars as I just want to understand the snappyHexMesh process.

The blockMeshDict file is as following (blue text):

backgroundMesh
{
xMin -20; // L = 350
xMax 50;
yMin -20; // L = 280
yMax 30;
zMin 0;
zMax 10;
xCells 20;
yCells 10;
zCells 5;
}

convertToMeters 1;

vertices
(
($:backgroundMesh.xMin $:backgroundMesh.yMin $:backgroundMesh.zMin)
($:backgroundMesh.xMax $:backgroundMesh.yMin $:backgroundMesh.zMin)
($:backgroundMesh.xMax $:backgroundMesh.yMax $:backgroundMesh.zMin)
($:backgroundMesh.xMin $:backgroundMesh.yMax $:backgroundMesh.zMin)

($:backgroundMesh.xMin $:backgroundMesh.yMin $:backgroundMesh.zMax)
($:backgroundMesh.xMax $:backgroundMesh.yMin $:backgroundMesh.zMax)
($:backgroundMesh.xMax $:backgroundMesh.yMax $:backgroundMesh.zMax)
($:backgroundMesh.xMin $:backgroundMesh.yMax $:backgroundMesh.zMax)
);

blocks
(
hex (0 1 2 3 4 5 6 7)
(
$:backgroundMesh.xCells
$:backgroundMesh.yCells
$:backgroundMesh.zCells
)
simpleGrading (1 1 1)
);

edges
(
);

boundary
(
inlet
{
type patch;
faces
(
(0 3 7 4)
);
}

outlet
{
type patch;
faces
(
(1 5 6 2)
);
}

ground
{
type wall;
faces
(
(0 1 2 3)
);
}

frontAndBack
{
type symmetry;
faces
(
(0 4 5 1)
(3 2 6 7)
(4 7 6 5)
);
}

);

mergePatchPairs
(
);

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

My snappyHexMesh is as following (red text):

castellatedMesh on;
snap on;
addLayers off;

geometry
{
try.stl
{
type triSurfaceMesh;
name try;
}

refinementBox
{
type searchableBox;
min ( -5 -10 0);
max (20 10 5);
}
};

castellatedMeshControls
{
features
(
{ file "try.eMesh"; level 3; }
);

refinementSurfaces
{
try
{
level (3 3);
patchInfo { type wall; }
}
}

refinementRegions
{
refinementBox
{
mode inside;
levels ((1E15 2));
}
}

locationInMesh (1 1 1);
}

snapControls
{
explicitFeatureSnap true;
implicitFeatureSnap false;
}

addLayersControls
{
layers
{
"CAD.*"
{
nSurfaceLayers 2;
}
}

relativeSizes true;
expansionRatio 1.2;
finalLayerThickness 0.5;
minThickness 1e-3;
}

meshQualityControls
{
}

writeFlags
(
scalarLevels
layerSets
layerFields
);

mergeTolerance 1e-6;

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

So I want to refine the mesh around my Geometry (the .stl file), defining a box inside the blockMesh domain.

The code runs, but when I analyse the Mesh with Paraview, there is no refinement.

May someone help me with this?

Thank you!

[fra_mat]

I apologize, the error was in the surfaceFeatureExtractDict file.

This thread can be eliminated