[chandra shekhar pant]

Hello All,


I am trying to simulate the blunt edge aerofoil with a slot at a distance of 60% chord length. I had created the stl file in gmsh and imported the stl into the snappyHexMesh. The mesh around the aerofoil looks steppy (in stairs, steep). The screen shot is attached herewith. I tried to check the stl file by using checkSurfaceMesh and it shows ok.
Complete report of the checkSurfaceMesh, snappyHexMesh, surfaceHexMesh dic along with the screenshot of mesh is attached.

Any help is highly welcomed.
Thanks!
Also snappyHexMesh gives this warning: Checking faces in error :
non-orthogonality > 65 degrees : 111012
faces with face pyramid volume < 1e-20 : 0
faces with face-decomposition tet quality < 1e-30 : 1
faces with concavity > 80 degrees : 827
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.02 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.02 : 2
faces on cells with determinant < 0.001 : 14
Finished meshing with 111856 illegal faces (concave, zero area or negative cell pyramid volume)

[pconen]

Hey,

I hope I understand your question in the right way.
If you want to get away the "stairs" you can try to turn

nCellsBetweenLevels 1;

into

nCellsBetweenLevels 2;

[chandra shekhar pant]

Hello Philipp Cone,



Many thanks for your reply, I will implement and update the results. Thanks again for the prompt and useful reply.

[chandra shekhar pant]

It haven't helped a lot. Please check the updated snap-shot herewith. I am just wondering if this issue may have been due to the meshing while creating the stl file in gmsh? Although not sure, if this is the issue of snappyHex or problem with the stl/gmsh. Any help or comment is highly welcomed.


Thanks!

[pconen]

Hey,

sorry I am really new to OpenFOAM.
Maybe this can help.
Can you tell me why you use

resolveFeatureAngle -1;

?
I never saw a negativ value for this.
I would try out

resolveFeatureAngle 30;

But remember to change the surfaceFeaturesDict and the surfaceFeaturesExtractDict.
With a value of 30 you have to use

 includedAngle   150;

[chandra shekhar pant]

Hello,

Thanks for your valuable inputs. I tried to changed whatever you have suggested like, changing the nCellsBetweenLevels, resolveFeatureAngle to 2 and 30 respectively. Also change the surfaceFeatureExtract dir. according to your inputs. The snappyHexMesh dic and the surfaceFeatureExtract dir along with the updated snapshot of the mesh is attached. Is that issue could be from the stl file ? The stl file is attached also herewith.



Again thanks for the help would highly appreciate any further inputs.


Thanks!

[pconen]

Hey,

of course, it could also be due to the STL. As already mentioned, I'm not a professional yet.
Unfortunately I have no experience with GMSH. I can only recommend the use of Salome. With it very solid surface nets can be created. For the meshing I recommend Mefisto for 2D. Under 1D you can choose Wire Distance. In a box next to it you can choose a local length of 0.002m for example. That depends on your component size.
Otherwise, I can't think of another quick tip for you. I am sorry.

[chandra shekhar pant]

Hello,
I am truly very happy and thankful with your inputs. These inputs have really helped me, I will try to switch these software. Thanks again!

[pconen]

Hey,

your welcome.
Fell free to ask further questions!