[kmurphy]

Hi,

I am trying to construct a mesh to represent a rectangular region with a single circular hole. I started with the mesh used in tutorials/potentalFoam/cylinder (which matches half of my required mesh) and tried to extend this. When I do not split the blocks, i.e., using

hex (# # # # # # # #) (1 1 1) simpleGrading(1 1 1)

for each block, everything seems OK. Program checkBlock does not generate any warnings and the resulting mesh looks like .

The output from checkMesh is in checkMesh_1.log and the blockMeshDict file is blockMeshDict_1 .

Adding grading to most of the blocks works as I would have expected, to produce the following mesh .

The output from checkMesh is in checkMesh_2.log and the blockMeshDict file is blockMeshDict_2 and everything appears OK.

However, when I try to add grading to any of the other blocks I do not obtain the expected (to me) mesh. In the following diagram and files I have only
added grading to two of the problem blocks to keep the mesh as simple as possible and the problem seem to be independent of the grading used in the other blocks.

.

The output from checkMesh is checkMesh_3.log and the blockMeshDict file is blockMeshDict_3 .

In the log of checkMesh there were reports of

High aspect ratio for cell 10: 4.44057e+199
High aspect ratio for cell 11: 4.44058e+199

and

Severe non-orthogonality detected for face ...

which doesn't appear to be too healthy.

I have tried to switch the order of the vertices without success. Any advice would be great,

Thanks, Kieran

Hi, Im currently working on such a geometry which I got thru the gmsh program and then converted it to OpenFOAM with help of the command gmshToFoam.

You could try this otherwise I dont know how to implement this.

/Erik