[cfMesh] improveMeshQuality

student666 · May 30, 2017, 01:49

Hi all,

currently I'm using cfMesh as meshing tool.
I found out the

Code:

improveMeshQuality

command, but I use it with default settings.

Code:

Default number of loops is 10
Default number of iterations is 50
Default number of surface iterations is 2
Using default quality threshold 0.1
No constraints applied on the smoothing procedure

With default settings, my feeling is that cfMesh tends to reduce the skewness of the faces, but generally non orthogonality increases, is it right?

Now I would like to ask if it is possibile to set it up to work on specific errors that arise with the

Code:

checkMesh

utility and if someone can share her/his experience with this option.

e.g. If checkMesh shows a huge non-ortogonality value (> 80) I would like to make

Code:

improveMeshQuality

work to reduce mainly this parameter.

This is a checkMesh after a single run with cartesianMesh

Code:

    Overall domain bounding box (-4.43393e-012 0.3235 -0.225969) (0.25 1.012 0.318)
    Mesh (non-empty, non-wedge) directions (1 1 1)
    Mesh (non-empty) directions (1 1 1)
    Boundary openness (-9.98583e-015 2.23961e-015 9.44861e-016) OK.
    Max cell openness = 8.63145e-016 OK.
    Max aspect ratio = 139.552 OK.
    Minimum face area = 2.08953e-009. Maximum face area = 2.05036e-005.  Face area magnitudes OK.
    Min volume = 4.11405e-013. Max volume = 5.75174e-008.  Total volume = 0.0601205.  Cell volumes OK.
    Mesh non-orthogonality Max: 88.6164 average: 6.16081
   *Number of severely non-orthogonal (> 70 degrees) faces: 407.
    Non-orthogonality check OK.
  <<Writing 407 non-orthogonal faces to set nonOrthoFaces
    Face pyramids OK.
Average skewness 0.0461554
 ***Max skewness = 5.73023, 266 highly skew faces detected which may impair the quality of the results
  <<Writing 266 skew faces to set skewFaces
    Coupled point location match (average 0) OK.

this is a first checkMesh output after improveMeshQuality has been run

Code:

Checking geometry...
    Overall domain bounding box (-4.43393e-012 0.3235 -0.225969) (0.25 1.012 0.318)
    Mesh (non-empty, non-wedge) directions (1 1 1)
    Mesh (non-empty) directions (1 1 1)
    Boundary openness (-1.03912e-014 2.23071e-015 9.15085e-016) OK.
    Max cell openness = 1.44747e-015 OK.
    Max aspect ratio = 55.296 OK.
    Minimum face area = 1.8092e-009. Maximum face area = 2.05036e-005.  Face area magnitudes OK.
    Min volume = 7.42452e-014. Max volume = 5.12649e-008.  Total volume = 0.0601206.  Cell volumes OK.
    Mesh non-orthogonality Max: 86.0757 average: 8.91875
   *Number of severely non-orthogonal (> 70 degrees) faces: 9567.
    Non-orthogonality check OK.
  <<Writing 9567 non-orthogonal faces to set nonOrthoFaces
    Face pyramids OK.
Average skewness 0.0406472
 ***Max skewness = 4.91016, 2 highly skew faces detected which may impair the quality of the results
  <<Writing 2 skew faces to set skewFaces
    Coupled point location match (average 0) OK.

Thanks in advance.

Regards.