[Technical] How does checkMesh calculate wrongOrientedFaces

yjhp1016 · February 27, 2024, 19:02

Hello, I often get error:" Error in face pyramids: 2 faces are incorrectly oriented." when I do checkMesh. I'm wondering how does checkMesh calculate this? I tried hard in reading source code, but didn't get the idea. Thanks for your help

cibanez · March 7, 2024, 04:24

Hi,

To check how those faces look like I recommend you to run checkMesh -writeAllFields to write the faceSet with wrongOrientedFaces and then foamToVTK -faceSet wrongOrientedFaces to generate a file in VTK/faceSet that you can open in Paraview.

A face is categorized as wrongly oriented when the facePyramidVolume becomes negative

Here you can see how a face is categorized as wrongly oriented when its face pyramid volume becames negative:
https://develop.openfoam.com/Develop...shCheck.C#L468

Here you can see how the pyramid volumes are being called. The signs make the volumes positive because the normal face vector points away from the owner cell:
https://develop.openfoam.com/Develop...shTools.C#L582

Here you can check the equation that is used to calculate the pyramid volume. This volume becomes negative when the angle between the height and the face normal is more than 90 degrees. The height vector goes from the face center to the cell center (apex) of either the neighbour or the owner cell.
https://develop.openfoam.com/Develop...pyramidI.H#L80

This means that you will get a wrongly oriented face when the normal of a face points in the opposite direction is supposed to. For a boundary face this means the face is so warped that the normal of the face points inwards to the computational domain instead of outwards.

February 27, 2024, 19:02	How does checkMesh calculate wrongOrientedFaces	#1
yjhp1016 New Member Paul Yang Join Date: Jun 2022 Posts: 2 Rep Power: 0	Hello, I often get error:" Error in face pyramids: 2 faces are incorrectly oriented." when I do checkMesh. I'm wondering how does checkMesh calculate this? I tried hard in reading source code, but didn't get the idea. Thanks for your help

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
How to calculate centerbody thrust of an aerospike nozzle in Fluent	Roh	FLUENT	5	February 21, 2022 12:30
Is there a way to the same aspect ratio from checkMesh using Paraview?	MazenDraw	OpenFOAM Pre-Processing	1	January 30, 2020 10:03
How to calculate Turbulent Intensity and Hydraulic Diameter for my problem?	Roh	FLUENT	1	June 17, 2018 04:12
calculate friction factor & nusselt number	soheil1991	FLUENT	3	March 11, 2017 10:30
How to calculate Torque for francis turbine	manish	CFX	4	March 15, 2007 03:57

March 7, 2024, 04:24		#2
cibanez New Member Cristóbal Join Date: Jan 2022 Location: Sweden Posts: 15 Rep Power: 4	Hi, To check how those faces look like I recommend you to run checkMesh -writeAllFields to write the faceSet with wrongOrientedFaces and then foamToVTK -faceSet wrongOrientedFaces to generate a file in VTK/faceSet that you can open in Paraview. A face is categorized as wrongly oriented when the facePyramidVolume becomes negative Here you can see how a face is categorized as wrongly oriented when its face pyramid volume becames negative: https://develop.openfoam.com/Develop...shCheck.C#L468 Here you can see how the pyramid volumes are being called. The signs make the volumes positive because the normal face vector points away from the owner cell: https://develop.openfoam.com/Develop...shTools.C#L582 Here you can check the equation that is used to calculate the pyramid volume. This volume becomes negative when the angle between the height and the face normal is more than 90 degrees. The height vector goes from the face center to the cell center (apex) of either the neighbour or the owner cell. https://develop.openfoam.com/Develop...pyramidI.H#L80 This means that you will get a wrongly oriented face when the normal of a face points in the opposite direction is supposed to. For a boundary face this means the face is so warped that the normal of the face points inwards to the computational domain instead of outwards.