[PetrosLazaris]

Hello to everyone.

I have used the checkMesh utility and it reported that I have some highly skewed faces and some severely non-orthogonal faces in my mesh.

This information is useful to me, however I want to know what is the face with the maximum skewness and the maximum non-orthogonality in my mesh. So is there any way to print the skewness/non-orthogonality field in my mesh or print the skewness/non-orthogonality of the problematic faces ?

Thank you very much in advance.

[unilord]

Originally Posted by PetrosLazaris

Hello to everyone.

I have used the checkMesh utility and it reported that I have some highly skewed faces and some severely non-orthogonal faces in my mesh.

This information is useful to me, however I want to know what is the face with the maximum skewness and the maximum non-orthogonality in my mesh. So is there any way to print the skewness/non-orthogonality field in my mesh or print the skewness/non-orthogonality of the problematic faces ?

Thank you very much in advance.

Hey

The command you should use is:

checkMesh -writeSets vtk

This will write the cell sets with nonOrthogonality or skewness higher than the defined in the meshQualityDict directory. It is written in a /postProcessing folder. You can then open paraview and open the .vtk files. Usually what I do to have a better visualisation is to have the whole geometry with a reduced opacity, and open the .vtk sets with "surface and edges".

Best Regards,
Pedro

P.S. - I recommend you to use the command

 checkMesh -allGeometry -allTopology

since it gives you more insights regarding the mesh.

[PetrosLazaris]

Thank you very much for your answer.

However, the utility you just told me to use does not print the skewness values, but only the id numbers of the highly skewed faces. I want the values of skewness of each face.

[unilord]

Originally Posted by PetrosLazaris

Thank you very much for your answer.

However, the utility you just told me to use does not print the skewness values, but only the id numbers of the highly skewed faces. I want the values of skewness of each face.

Ok I understand what you want now. I know that is possible to do, since I already did it. However, I believe only with ESI versions. Since I stopped using that function a long time ago since I do my pre processing in SALOME, I can't really help you. I remember I saw it in a video of József Nagy in youtube.

Let's hope someone else knows how to do this.

[Yann]

In the OpenCFD branch (openfoam.com) you can use the -writeAllFields and -writeAllSurfaceFields options when running checkMesh.
https://www.openfoam.com/news/main-n...sing#checkmesh

This is not available in the foundation branch (openfoam.org)

[PetrosLazaris]

Thak you very much for your answer.

Is there any way to do this in OpenFOAM.org too, by writing a custom c++ script ?