CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Meshing & Mesh Conversion

[Commercial meshers] non-orthogonality errors, incorrectly oriented faces

Register Blogs Community New Posts Updated Threads Search

Like Tree17Likes
  • 2 Post By fertinaz
  • 1 Post By franjo_j
  • 6 Post By franjo_j
  • 1 Post By Mueye
  • 7 Post By Rasmusiwersen

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   July 26, 2010, 12:18
Default non-orthogonality errors, incorrectly oriented faces
  #1
Member
 
Markus Weinmann
Join Date: Mar 2009
Location: Stuttgart, Germany
Posts: 77
Rep Power: 17
cfdmarkus is on a distinguished road
Hi all,

I am trying to load a fairly complex ccm mesh.
The mesh has been tested in StarCCM and seems to work fine.

When I load the mesh in OpenFOAM I get a bunch off errors and warnings (see below). The case does run for a few iterations before the simulation crashes.

I am particulary concerned about those 2 errors:

Error in face pyramids: 143 faces are incorrectly oriented
Number of non-orthogonality errors: 15.

How important are they and do they need fixing in StarCCM?

Any help is appreciated.

Markus






Checking topology...
Boundary definition OK.
Point usage OK.
<<Found 1101 neighbouring cells with multiple inbetween faces.
Upper triangular ordering OK.
<<Writing 2297 unordered faces to set upperTriangularFace
Topological cell zip-up check OK.
Face vertices OK.
***Number of duplicate (not baffle) faces found: 306
***Number of faces with non-consecutive shared points: 540
<<Writing 1405 faces with incorrect edges to set edgeFaces
Number of regions: 1 (OK).



Checking geometry...
Overall domain bounding box (-2.18829e-07 -7.47931 -0.277615) (2.86 12.042 4.66777)
Mesh (non-empty, non-wedge) directions (1 1 1)
Mesh (non-empty) directions (1 1 1)
Boundary openness (-3.71575e-13 -8.51015e-14 -1.25127e-12) OK.
Max cell openness = 1.52626e-15 OK.
Max aspect ratio = 131.378 OK.
***Zero or negative face area detected. Minimum area: 0
<<Writing 1 zero area faces to set zeroAreaFaces
Min volume = 7.78059e-14. Max volume = 6.85685e-05. Total volume = 254.056. Cell volumes OK.
Mesh non-orthogonality Max: 121.889 average: 4.92652
*Number of severely non-orthogonal faces: 2880.
***Number of non-orthogonality errors: 15.
<<Writing 2895 non-orthogonal faces to set nonOrthoFaces
***Error in face pyramids: 143 faces are incorrectly oriented.
<<Writing 140 faces with incorrect orientation to set wrongOrientedFaces
***Max skewness = 88.8617, 233 highly skew faces detected which may impair the quality of the results
<<Writing 233 skew faces to set skewFaces
*Edges too small, min/max edge length = 0 0.0543461, number too small: 34
<<Writing 34 points on short edges to set shortEdges
*There are 5188 faces with concave angles between consecutive edges. Max concave angle = 89.9778 degrees.
<<Writing 5188 faces with concave angles to set concaveFaces
Face flatness (1 = flat, 0 = butterfly) : average = 0.999981 min = 0.0551237
*There are 447 faces with ratio between projected and actual area < 0.8
Minimum ratio (minimum flatness, maximum warpage) = 0.0551237
<<Writing 447 warped faces to set warpedFaces
Cell determinant (wellposedness) : minimum: 2.11744e-08 average: 8.24052
***Cells with small determinant found, number of cells: 189153
<<Writing 189153 under-determined cells to set underdeterminedCells

Failed 6 mesh checks.
cfdmarkus is offline   Reply With Quote

Old   July 27, 2010, 03:06
Default
  #2
Senior Member
 
sega's Avatar
 
Sebastian Gatzka
Join Date: Mar 2009
Location: Frankfurt, Germany
Posts: 729
Rep Power: 20
sega is on a distinguished road
Can you show the mesh how it is displayed in ParaView?
Anything strange on first glance?

Or can you even post the complete mesh?
__________________
Schrödingers wife: "What did you do to the cat? It's half dead!"
sega is offline   Reply With Quote

Old   July 27, 2010, 04:32
Default
  #3
Member
 
Markus Weinmann
Join Date: Mar 2009
Location: Stuttgart, Germany
Posts: 77
Rep Power: 17
cfdmarkus is on a distinguished road
Hi

the mesh has more than 30 mio. cells, so it is a bit hard to post the mesh here. Also, there is nothing very obvious wrong with the mesh when displayed in paraview.

The important question for me is which errors/warnings can be accepted/-ignored (skewed, non-orthogonal cells, multiple inbetween faces, unordered faces) and which need to be fixed in the grid generator (face orientation?, ...).

Thanks,
Markus
cfdmarkus is offline   Reply With Quote

Old   October 20, 2010, 06:58
Thumbs down
  #4
New Member
 
Mark Herskind
Join Date: Oct 2010
Posts: 12
Rep Power: 15
MarkHerskind is on a distinguished road
No comments to this?

I have a similar problem with a mesh.. I'm trying to use snappyHexMesh on a terrain (large scale CFD for wind farm purposes), and when I have the sufficient cell size, I get errors:

Code:
 
Checking geometry...
    Overall domain bounding box (404868 4758933 490) (420868 4774933 7500)
    Mesh (non-empty, non-wedge) directions (1 1 1)
    Mesh (non-empty) directions (1 1 1)
    Boundary openness (-4.989805e-18 -2.407661e-17 -3.776756e-14) OK.
 ***High aspect ratio cells found, Max aspect ratio: 3.113e+203, number of cells 19
  <<Writing 19 cells with high aspect ratio to set highAspectRatioCells
    Minumum face area = 1.61935. Maximum face area = 23599.98.  Face area magnitudes OK.
    Min volume = 2e-300. Max volume = 1737135.  Total volume = 1.317329e+12.  Cell volumes OK.
    Mesh non-orthogonality Max: 179.9979 average: 7.740474
   *Number of severely non-orthogonal faces: 141.
 ***Number of non-orthogonality errors: 151.
  <<Writing 292 non-orthogonal faces to set nonOrthoFaces
 ***Error in face pyramids: 280 faces are incorrectly oriented.
  <<Writing 164 faces with incorrect orientation to set wrongOrientedFaces
 ***Max skewness = 5.783191, 3 highly skew faces detected which may impair the quality of the results
  <<Writing 3 skew faces to set skewFaces
Failed 4 mesh checks.
Unfortunately, I can not share the mesh due to confidentiality reasons..

I'm wondering about the severity of these failed checks..
I my opinion the aspects ratio cells are not really of concern, as they are the remote cells in the very top of the domain..
The nonOrthogonality is neither of much concern I believe..
The skewness is only 3 cells, which I have located outside my area of interest

The one I'm wondering about is the incorrectly oriented faces.. How serious is this, and what would you recommend doing to avoid this error?
MarkHerskind is offline   Reply With Quote

Old   February 23, 2011, 11:32
Default
  #5
Member
 
Fatih Ertinaz
Join Date: Feb 2011
Location: Istanbul
Posts: 64
Rep Power: 15
fertinaz is on a distinguished road
Hi.

I don't know if you still have the same problem since it is bit outdated however I would like to say that your computations will not be able to convergence as long as wrong oriented faces exist in your mesh. You may skip the others, like non-orthogonal faces and skewness etc.

If you could solve this issue, I would like to know your snappy settings for the boundary layer addition part in the snappyHexMeshDict. I have a similar trouble and now I'm trying to figure out what the optimal settings are.

Best,
Fatih.
ashish.vinayak and TGS like this.
fertinaz is offline   Reply With Quote

Old   February 8, 2012, 11:35
Default
  #6
Member
 
Aqua
Join Date: Oct 2011
Posts: 96
Rep Power: 14
aqua is on a distinguished road
Quote:
Originally Posted by fertinaz View Post
Hi.

I don't know if you still have the same problem since it is bit outdated however I would like to say that your computations will not be able to convergence as long as wrong oriented faces exist in your mesh. You may skip the others, like non-orthogonal faces and skewness etc.

If you could solve this issue, I would like to know your snappy settings for the boundary layer addition part in the snappyHexMeshDict. I have a similar trouble and now I'm trying to figure out what the optimal settings are.

Best,
Fatih.
Hey, Fatih,
I have the same problem now, after snappyHexMesh, some error like "incorrectly oriented faces" showed up.
Could you please tell me how to solve the problem?
Thank you so much!
Aqua
aqua is offline   Reply With Quote

Old   August 9, 2014, 06:21
Default
  #7
Member
 
Ye Zhang
Join Date: Dec 2009
Location: Delft,Netherland
Posts: 92
Rep Power: 16
kiddmax is on a distinguished road
Hi all,

Do you have any solutions now? I also have the same problem "incorrectly oriented faces" in my mesh. I do not know why. I generated the mesh by using Pointwise.

Any suggestions are appreciated!

Ye
kiddmax is offline   Reply With Quote

Old   November 17, 2014, 16:00
Default
  #8
Senior Member
 
Franjo Juretic
Join Date: Aug 2011
Location: Velika Gorica, Croatia
Posts: 124
Rep Power: 16
franjo_j is on a distinguished road
Send a message via Skype™ to franjo_j
Hello,

What is the format (ascii or binary) is used for storing the OpenFOAM mesh after conversion? In general, it is best to use binary because it preserves all information about the vertex coordinates.
This is especially critical in bnd layers, which may be of the same thickness as the truncated part of the vertex coordinates. Therefore, it may happen that a valid mesh becomes invalid after conversion.
Incorrectly oriented faces mean that the mesh is invalid because there exist face normals pointing into the owner cell instead of outside. This causes conservation problems because the fluxes through such face contribute to the cell with a wrong sign.
I hope this helps a bit.

Kind Regards,

Franjo

Quote:
Originally Posted by kiddmax View Post
Hi all,

Do you have any solutions now? I also have the same problem "incorrectly oriented faces" in my mesh. I do not know why. I generated the mesh by using Pointwise.

Any suggestions are appreciated!

Ye
Stefanie.S.W. likes this.
franjo_j is offline   Reply With Quote

Old   June 8, 2015, 07:19
Default
  #9
Senior Member
 
Join Date: Mar 2015
Posts: 250
Rep Power: 12
KateEisenhower is on a distinguished road
Quote:
Originally Posted by cfdmarkus View Post
Hi all,

I am trying to load a fairly complex ccm mesh.
The mesh has been tested in StarCCM and seems to work fine.

When I load the mesh in OpenFOAM I get a bunch off errors and warnings (see below). The case does run for a few iterations before the simulation crashes.

I am particulary concerned about those 2 errors:

Error in face pyramids: 143 faces are incorrectly oriented
Number of non-orthogonality errors: 15.

How important are they and do they need fixing in StarCCM?

Any help is appreciated.

Markus






Checking topology...
Boundary definition OK.
Point usage OK.
<<Found 1101 neighbouring cells with multiple inbetween faces.
Upper triangular ordering OK.
<<Writing 2297 unordered faces to set upperTriangularFace
Topological cell zip-up check OK.
Face vertices OK.
***Number of duplicate (not baffle) faces found: 306
***Number of faces with non-consecutive shared points: 540
<<Writing 1405 faces with incorrect edges to set edgeFaces
Number of regions: 1 (OK).



Checking geometry...
Overall domain bounding box (-2.18829e-07 -7.47931 -0.277615) (2.86 12.042 4.66777)
Mesh (non-empty, non-wedge) directions (1 1 1)
Mesh (non-empty) directions (1 1 1)
Boundary openness (-3.71575e-13 -8.51015e-14 -1.25127e-12) OK.
Max cell openness = 1.52626e-15 OK.
Max aspect ratio = 131.378 OK.
***Zero or negative face area detected. Minimum area: 0
<<Writing 1 zero area faces to set zeroAreaFaces
Min volume = 7.78059e-14. Max volume = 6.85685e-05. Total volume = 254.056. Cell volumes OK.
Mesh non-orthogonality Max: 121.889 average: 4.92652
*Number of severely non-orthogonal faces: 2880.
***Number of non-orthogonality errors: 15.
<<Writing 2895 non-orthogonal faces to set nonOrthoFaces
***Error in face pyramids: 143 faces are incorrectly oriented.
<<Writing 140 faces with incorrect orientation to set wrongOrientedFaces
***Max skewness = 88.8617, 233 highly skew faces detected which may impair the quality of the results
<<Writing 233 skew faces to set skewFaces
*Edges too small, min/max edge length = 0 0.0543461, number too small: 34
<<Writing 34 points on short edges to set shortEdges
*There are 5188 faces with concave angles between consecutive edges. Max concave angle = 89.9778 degrees.
<<Writing 5188 faces with concave angles to set concaveFaces
Face flatness (1 = flat, 0 = butterfly) : average = 0.999981 min = 0.0551237
*There are 447 faces with ratio between projected and actual area < 0.8
Minimum ratio (minimum flatness, maximum warpage) = 0.0551237
<<Writing 447 warped faces to set warpedFaces
Cell determinant (wellposedness) : minimum: 2.11744e-08 average: 8.24052
***Cells with small determinant found, number of cells: 189153
<<Writing 189153 under-determined cells to set underdeterminedCells

Failed 6 mesh checks.
Hello all,

Does anyone know if the number of stars (*, ** or ***) before the error message indicates the severeness of the problem. Let's say if a complaint goes with ***, it should be fixed under any cirmumstances?

Best regards,

Kate
KateEisenhower is offline   Reply With Quote

Old   June 9, 2015, 10:50
Default
  #10
Senior Member
 
Franjo Juretic
Join Date: Aug 2011
Location: Velika Gorica, Croatia
Posts: 124
Rep Power: 16
franjo_j is on a distinguished road
Send a message via Skype™ to franjo_j
Hello,

Quote:
Originally Posted by KateEisenhower View Post
Hello all,

Does anyone know if the number of stars (*, ** or ***) before the error message indicates the severeness of the problem. Let's say if a complaint goes with ***, it should be fixed under any cirmumstances?
Your assumption is correct. The errors with three stars are the most important.
From my experience with OpenFOAM, the most important check is for face pyramids, and they must not be negative. If this check fails, it is almost impossible to get a meaningful solution because the discretisation is no longer conservative. Zero area faces make the mesh "disconnected" locally, and that can also cause problems with convergence.
Most other metrics like non-orthogonality, skewness, flatness can cause problems when they are located in regions of high gradient variation, like boundary layers, shear layers, etc.
Please bear in mind that it is extremely difficult to generate meshes satisfying all quality metrics, especially in complex geometries, and with small number of cells, so my suggestion is to focus on face pyramids, and zero are faces first.
I hope this helps a bit.

Regards,

Franjo
__________________
Principal Developer of cfMesh and CF-MESH+
www.cfmesh.com
Social media: LinkedIn, Twitter, YouTube, Facebook, Pinterest, Instagram
franjo_j is offline   Reply With Quote

Old   November 21, 2015, 03:37
Default
  #11
Member
 
Saurabh Tandon
Join Date: Nov 2015
Location: Austin
Posts: 43
Rep Power: 10
STutexas is on a distinguished road
Hi

I am a new user of Open foam. This is a very interesting post. Could you let me know if there is a way to change the mesh in Openfoam after I get the errors or would I have to remesh.

Saurabh
STutexas is offline   Reply With Quote

Old   December 6, 2016, 09:05
Default
  #12
New Member
 
Omur Icke
Join Date: Apr 2016
Posts: 6
Rep Power: 10
eromir is on a distinguished road
To avoid from this problem you can use command as below.

minVol 1e-13;
minTetQuality 1e-16;
minTwist 0.05;
eromir is offline   Reply With Quote

Old   April 13, 2017, 11:13
Default
  #13
New Member
 
Muye Ge
Join Date: Feb 2016
Location: Sweden
Posts: 2
Rep Power: 0
Mueye is on a distinguished road
Hi!

From your post I think it is a write precision problem.

When you converting the mesh from a .ccm file, try to change the write precision to 15 in system/controldict.

Regards
lourencosm likes this.
Mueye is offline   Reply With Quote

Old   September 20, 2019, 05:05
Default
  #14
Member
 
Rasmus Iwersen
Join Date: Jan 2019
Location: Denmark
Posts: 81
Rep Power: 8
Rasmusiwersen is on a distinguished road
Quote:
Originally Posted by cfdmarkus View Post
Hi all,

I am trying to load a fairly complex ccm mesh.
The mesh has been tested in StarCCM and seems to work fine.

When I load the mesh in OpenFOAM I get a bunch off errors and warnings (see below). The case does run for a few iterations before the simulation crashes.

I am particulary concerned about those 2 errors:

Error in face pyramids: 143 faces are incorrectly oriented
Number of non-orthogonality errors: 15.

How important are they and do they need fixing in StarCCM?

Any help is appreciated.

Markus






Checking topology...
Boundary definition OK.
Point usage OK.
<<Found 1101 neighbouring cells with multiple inbetween faces.
Upper triangular ordering OK.
<<Writing 2297 unordered faces to set upperTriangularFace
Topological cell zip-up check OK.
Face vertices OK.
***Number of duplicate (not baffle) faces found: 306
***Number of faces with non-consecutive shared points: 540
<<Writing 1405 faces with incorrect edges to set edgeFaces
Number of regions: 1 (OK).



Checking geometry...
Overall domain bounding box (-2.18829e-07 -7.47931 -0.277615) (2.86 12.042 4.66777)
Mesh (non-empty, non-wedge) directions (1 1 1)
Mesh (non-empty) directions (1 1 1)
Boundary openness (-3.71575e-13 -8.51015e-14 -1.25127e-12) OK.
Max cell openness = 1.52626e-15 OK.
Max aspect ratio = 131.378 OK.
***Zero or negative face area detected. Minimum area: 0
<<Writing 1 zero area faces to set zeroAreaFaces
Min volume = 7.78059e-14. Max volume = 6.85685e-05. Total volume = 254.056. Cell volumes OK.
Mesh non-orthogonality Max: 121.889 average: 4.92652
*Number of severely non-orthogonal faces: 2880.
***Number of non-orthogonality errors: 15.
<<Writing 2895 non-orthogonal faces to set nonOrthoFaces
***Error in face pyramids: 143 faces are incorrectly oriented.
<<Writing 140 faces with incorrect orientation to set wrongOrientedFaces
***Max skewness = 88.8617, 233 highly skew faces detected which may impair the quality of the results
<<Writing 233 skew faces to set skewFaces
*Edges too small, min/max edge length = 0 0.0543461, number too small: 34
<<Writing 34 points on short edges to set shortEdges
*There are 5188 faces with concave angles between consecutive edges. Max concave angle = 89.9778 degrees.
<<Writing 5188 faces with concave angles to set concaveFaces
Face flatness (1 = flat, 0 = butterfly) : average = 0.999981 min = 0.0551237
*There are 447 faces with ratio between projected and actual area < 0.8
Minimum ratio (minimum flatness, maximum warpage) = 0.0551237
<<Writing 447 warped faces to set warpedFaces
Cell determinant (wellposedness) : minimum: 2.11744e-08 average: 8.24052
***Cells with small determinant found, number of cells: 189153
<<Writing 189153 under-determined cells to set underdeterminedCells

Failed 6 mesh checks.
Hi Markus (and fellow readers)

I've found a solution to some similar errors as seen below:

Checking geometry...
Overall domain bounding box (-2.99663862296822e-95 0 0) (50 50 1)
Mesh has 3 geometric (non-empty/wedge) directions (1 1 1)
Mesh has 3 solution (non-empty) directions (1 1 1)
Boundary openness (2.92902493088356e-18 -2.15183472189403e-18 9.8378233986928e-16) OK.
Max cell openness = 3.76326914376673e-16 OK.
Max aspect ratio = 36.8327501227017 OK.
Minimum face area = 1.13891991650133e-07. Maximum face area = 0.781677819197455. Face area magnitudes OK.
Min volume = 4.15463340177619e-09. Max volume = 0.156337183992222. Total volume = 2496.89481044183. Cell volumes OK.
Mesh non-orthogonality Max: 133.637450973096 average: 8.3113546910368
*Number of severely non-orthogonal (> 70 degrees) faces: 15.
***Number of non-orthogonality errors: 4.
<<Writing 19 non-orthogonal faces to set nonOrthoFaces
***Error in face pyramids: 70 faces are incorrectly oriented.
<<Writing 67 faces with incorrect orientation to set wrongOrientedFaces
***Max skewness = 14.8751178328007, 27 highly skew faces detected which may impair the quality of the results
<<Writing 27 skew faces to set skewFaces
Coupled point location match (average 0) OK.

Failed 3 mesh checks.

"Number of non-orthogonlaity errors" was removed by increasing the number of cells between two refinement levels. In my case, i have two refinement regions, and the number of butter layers was set to 3. Increasing it to 5 removed this error.

"Error in face pyramids" was removed by lowering the explicit feature edge refinement from 6 to 5. I am not entirely convinced why this helped. Also, this removed the "Number of severely non-orthogonal" error.

"Max skewness" was removed by setting the surface refinement min and max to the same level, in this case going from level (5 6) to level (5 5).

As you can see from the following checkMesh, the mesh passes the check and is ready for computation. I hope this helps any who might encounter the same challenge.

Checking geometry...
Overall domain bounding box (-1.4981364335015e-95 0 0) (50 50 1)
Mesh has 3 geometric (non-empty/wedge) directions (1 1 1)
Mesh has 3 solution (non-empty) directions (1 1 1)
Boundary openness (1.42427331543913e-18 5.28677104697883e-18 2.49325550741594e-15) OK.
Max cell openness = 3.69995665943141e-16 OK.
Max aspect ratio = 9.31165615736311 OK.
Minimum face area = 6.11134067207085e-05. Maximum face area = 0.781674938655092. Face area magnitudes OK.
Min volume = 1.19560045411929e-06. Max volume = 0.156336729887031. Total volume = 2496.89481965278. Cell volumes OK.
Mesh non-orthogonality Max: 39.4652589611894 average: 8.0090802589223
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 3.55799244148157 OK.
Coupled point location match (average 0) OK.

Mesh OK.

Hope this helps!
Rasmusiwersen is offline   Reply With Quote

Reply


Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
[snappyHexMesh] Help with Snappy: no layers growing GianF OpenFOAM Meshing & Mesh Conversion 2 September 23, 2020 08:26
[Gmsh] Vertex numbering is dense KateEisenhower OpenFOAM Meshing & Mesh Conversion 7 August 3, 2015 10:49
dynamic Mesh is faster than MRF???? sharonyue OpenFOAM Running, Solving & CFD 14 August 26, 2013 07:47
[Commercial meshers] fluentMeshToFoam multidomain mesh conversion problem Attesz OpenFOAM Meshing & Mesh Conversion 12 May 2, 2013 10:52
[snappyHexMesh] external flow with snappyHexMesh chelvistero OpenFOAM Meshing & Mesh Conversion 11 January 15, 2010 19:43


All times are GMT -4. The time now is 06:48.