aminem |
April 16, 2015 11:38 |
sHM quality of multi-region
3 Attachment(s)
Hi,
I have check my mesh with checkMesh -allGeometry -allTopology
I want to understand what does mean : 9122 illegal face (concave, zero area or negative cell pyramid mesh) and if it's affect my simulations?
How I can resolve the problem of concave mesh or negative cell pyramid mesh?
You can see in picture good mesh but when I check my mesh I found 2 failed mesh.
Thanks
Code:
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.3.0-f5222ca19ce6
Exec : checkMesh -allGeometry -allTopology
Date : Apr 16 2015
Time : 16:34:19
Host : "mesbahi-VirtualBox"
PID : 29101
Case : /home/mesbahi/OpenFOAM/mesbahi-2.3.0/run/travail/spotOFtest/spotOF
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create polyMesh for time = 0
Enabling all (cell, face, edge, point) topology checks.
Enabling all geometry checks.
Time = 0
Mesh stats
points: 3192315
faces: 9158547
internal faces: 9120147
cells: 2983210
faces per cell: 6.12719
boundary patches: 71
point zones: 0
face zones: 35
cell zones: 35
Overall number of cells of each type:
hexahedra: 2876652
prisms: 0
wedges: 0
pyramids: 0
tet wedges: 0
tetrahedra: 0
polyhedra: 106558
Breakdown of polyhedra by number of faces:
faces number of cells
6 30695
9 40889
12 21531
15 11246
18 2196
21 1
Checking topology...
Boundary definition OK.
Cell to face addressing OK.
Point usage OK.
Upper triangular ordering OK.
Face vertices OK.
Topological cell zip-up check OK.
Face-face connectivity OK.
Number of regions: 1 (OK).
Checking patch topology for multiply connected surfaces...
Patch Faces Points Surface topology Bounding box
boundaries 38400 38402 ok (closed singly connected) (-0.25 -0.25 -0.15) (0.25 0.25 0.35)
ALU 0 0 ok (empty)
ALU_slave 0 0 ok (empty)
Dissipateur 0 0 ok (empty)
Dissipateur_slave 0 0 ok (empty)
LED1 0 0 ok (empty)
LED1_slave 0 0 ok (empty)
LED10 0 0 ok (empty)
LED10_slave 0 0 ok (empty)
LED11 0 0 ok (empty)
LED11_slave 0 0 ok (empty)
LED12 0 0 ok (empty)
LED12_slave 0 0 ok (empty)
LED13 0 0 ok (empty)
LED13_slave 0 0 ok (empty)
LED14 0 0 ok (empty)
LED14_slave 0 0 ok (empty)
LED15 0 0 ok (empty)
LED15_slave 0 0 ok (empty)
LED16 0 0 ok (empty)
LED16_slave 0 0 ok (empty)
LED2 0 0 ok (empty)
LED2_slave 0 0 ok (empty)
LED3 0 0 ok (empty)
LED3_slave 0 0 ok (empty)
LED4 0 0 ok (empty)
LED4_slave 0 0 ok (empty)
LED5 0 0 ok (empty)
LED5_slave 0 0 ok (empty)
LED6 0 0 ok (empty)
LED6_slave 0 0 ok (empty)
LED7 0 0 ok (empty)
LED7_slave 0 0 ok (empty)
LED8 0 0 ok (empty)
LED8_slave 0 0 ok (empty)
LED9 0 0 ok (empty)
LED9_slave 0 0 ok (empty)
MCPCB 0 0 ok (empty)
MCPCB_slave 0 0 ok (empty)
source1 0 0 ok (empty)
source1_slave 0 0 ok (empty)
source10 0 0 ok (empty)
source10_slave 0 0 ok (empty)
source11 0 0 ok (empty)
source11_slave 0 0 ok (empty)
source12 0 0 ok (empty)
source12_slave 0 0 ok (empty)
source13 0 0 ok (empty)
source13_slave 0 0 ok (empty)
source14 0 0 ok (empty)
source14_slave 0 0 ok (empty)
source15 0 0 ok (empty)
source15_slave 0 0 ok (empty)
source16 0 0 ok (empty)
source16_slave 0 0 ok (empty)
source2 0 0 ok (empty)
source2_slave 0 0 ok (empty)
source3 0 0 ok (empty)
source3_slave 0 0 ok (empty)
source4 0 0 ok (empty)
source4_slave 0 0 ok (empty)
source5 0 0 ok (empty)
source5_slave 0 0 ok (empty)
source6 0 0 ok (empty)
source6_slave 0 0 ok (empty)
source7 0 0 ok (empty)
source7_slave 0 0 ok (empty)
source8 0 0 ok (empty)
source8_slave 0 0 ok (empty)
source9 0 0 ok (empty)
source9_slave 0 0 ok (empty)
Checking geometry...
Overall domain bounding box (-0.25 -0.25 -0.15) (0.25 0.25 0.35)
Mesh (non-empty, non-wedge) directions (1 1 1)
Mesh (non-empty) directions (1 1 1)
Boundary openness (-8.735056e-17 -1.252181e-15 -3.096387e-15) OK.
Max cell openness = 4.455474e-16 OK.
Max aspect ratio = 16.27134 OK.
Minimum face area = 4.189e-10. Maximum face area = 3.932657e-05. Face area magnitudes OK.
Min volume = 8.184723e-14. Max volume = 2.461116e-07. Total volume = 0.125. Cell volumes OK.
Mesh non-orthogonality Max: 75.63721 average: 9.745855
*Number of severely non-orthogonal (> 70 degrees) faces: 23.
Non-orthogonality check OK.
<<Writing 23 non-orthogonal faces to set nonOrthoFaces
Face pyramids OK.
Max skewness = 2.463072 OK.
Coupled point location match (average 0) OK.
***Error in face tets: 73 faces with low quality or negative volume decomposition tets.
<<Writing 73 faces with low quality or negative volume decomposition tets to set lowQualityTetFaces
*Edges too small, min/max edge length = 1.615488e-07 0.006278598, number too small: 6
<<Writing 6 points on short edges to set shortEdges
*There are 21518 faces with concave angles between consecutive edges. Max concave angle = 89.91077 degrees.
<<Writing 21518 faces with concave angles to set concaveFaces
Face flatness (1 = flat, 0 = butterfly) : min = 0.4107551 average = 0.9993399
*There are 82 faces with ratio between projected and actual area < 0.8
Minimum ratio (minimum flatness, maximum warpage) = 0.4107551
<<Writing 82 warped faces to set warpedFaces
Cell determinant (wellposedness) : minimum: 0.00615444 average: 11.70074
Cell determinant check OK.
***Concave cells (using face planes) found, number of cells: 87551
<<Writing 87551 concave cells to set concaveCells
Face interpolation weight : minimum: 0.05250315 average: 0.478492
Face interpolation weight check OK.
Face volume ratio : minimum: 0.01658261 average: 0.9001287
Face volume ratio check OK.
Failed 2 mesh checks.
End
|