[Lligres]

Howdy all,
This is my first post, so I hope I am doing everithing wright .
I'm designing a 3D mesh of a wing using gmsh in order to run it on Openfoam. I made a hibrid mesh, so the elements around the wing are hexahedra and the elements far from it are tetrahedra; between those hexahedra and tetrahedra there is a layer of pyramids in order to combine the quadrangle faces of the structured mesh with the tetrahedra nodes (It's done by default with the Delaunay algo).
When running the algorithm I have problems, I have two error messages which are the following: "Warning : Strange edge cavity (tet is deleted)" and "Warning : 12 ill-shaped tets are still in the mesh". I managed to reduce the number of ill-shapped tets from 2000 to 12 by reducing the aspect ratio of the quadrangles in contact with the un-structured volumes (this way the pyramids and the tetrahedra sorrounding it seem to be simpler and the algorithm seems to work better), but there are still cells that have 0 or negative volume according to OpenFOAM.

On the other hand, when runing the checkMesh code on openFoam i have several errors : " ***Zero or negative cell volume detected. Minimum negative volume: -0.00053318274, Number of negative volume cells: 83 ", "***Open cells found, max cell openness: 1, number of open cells 432", " ***Number of non-orthogonality errors: 373." , " ***Error in face pyramids: 660 faces are incorrectly oriented." and "***Max skewness = 45.903898, 27 highly skew faces detected which may impair the quality of the results".

I have been working on this for several days but now I'm at a point that I don't know how to advance further. I have attached the .geo file splitted in two .txt files, as it was to heavy to upload in one.

Any help regarding these issues would be very appreciated,
Kind regards,

Sergi Giménez Herrero

[Lligres]

Next I will attach some images of the mesh, the wing and the computational domain so it may be easier for you to understand the case.

If you pay atention to the wing, you can see that it has a "c" shaped end, it's because it is a c-wing and I am studying the effects of this wingtip.

Hope it can help a bit,
Kind regards to anyone reading this,

Sergi Giménez Herrero