CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Meshing Format & General Technical

checkMesh states mesh is OK yet faces are bent pentagons

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   March 20, 2014, 05:40
Default checkMesh states mesh is OK yet faces are bent pentagons
  #1
New Member
 
Richel Bilderbeek
Join Date: Nov 2013
Posts: 14
Rep Power: 3
assert is on a distinguished road
My question is: how can this mesh be correct (and how can the original mesh be incorrect)?

This is a follow-up of checkMesh states that faces do not form a cell, mesh appears correct , where I had a AFAIKS correct mesh, yet checkMesh gave two errors.

On that mesh I did the following commands:
Code:
renumberMesh
zipUpMesh
This resulted in a 'correct' (that is: according to checkMesh) mesh, with six faces, two of which are bent polygons. The image below shows the point indices, the others show the face indices. Of the bent polygon faces, I number the vertex indices.



How can this mesh be correct (and how can the original mesh be incorrect)?

Below I put the content of the AFAIK most important files, you can view all at http://www.richelbilderbeek.nl/ToolO...dCorrected.htm .

checkMesh output (I removed some text, indicated with [...]):

Code:
[...]
Build  : 2.1-88b2f2ae3a0b
Exec   : checkmesh
Date   : Mar 20 2014
Time   : 08:03:54
Host   : "AIRBEAR-W7-I5"
PID    : 2764
Case   : D:/Projects/Test/ToolOpenFoamErrorOpenCellsFoundCorrected
nProcs : 1
SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create polyMesh for time = 0

[...]

Time = 1

Mesh stats
    points:           8
    internal points:  0
    faces:            8
    internal faces:   0
    cells:            2
    boundary patches: 6
    point zones:      0
    face zones:       0
    cell zones:       0

Overall number of cells of each type:
    hexahedra:     0
    prisms:        0
    wedges:        0
    pyramids:      0
    tet wedges:    0
    tetrahedra:    0
    polyhedra:     2

Checking topology...
    Boundary definition OK.
    Cell to face addressing OK.
    Point usage OK.
    Upper triangular ordering OK.
    Face vertices OK.
   *Number of regions: 2
    The mesh has multiple regions which are not connected by any face.
  <<Writing region information to "1/cellToRegion"

Checking patch topology for multiply connected surfaces ...
    Patch               Faces    Points   Surface topology                  
    back                1        4        ok (non-closed singly connected)  
    bottom              2        4        ok (non-closed singly connected)  
    front               1        4        ok (non-closed singly connected)  
    left                1        4        ok (non-closed singly connected)  
    right               1        4        ok (non-closed singly connected)  
    top                 2        6        ok (non-closed singly connected)  

Checking geometry...
    Overall domain bounding box (-1 -1 0) (1 1 1)
    Mesh (non-empty, non-wedge) directions (1 1 1)
    Mesh (non-empty) directions (1 1 1)
    Boundary openness (-2.28837e-018 -1.83069e-017 0) OK.
    Max cell openness = 5.55112e-017 OK.
    Max aspect ratio = 2 OK.
    Minimum face area = 1. Maximum face area = 2.23607.  Face area magnitudes OK.
    Min volume = 0.731939. Max volume = 0.731939.  Total volume = 1.46388.  Cell volumes OK.
    Non-orthogonality check OK.
    Face pyramids OK.
    Max skewness = 0.34936 OK.
    Coupled point location match (average 0) OK.

Mesh OK.

End
points, as shown at the top-left of the image:

Code:
FoamFile
{
    version     2.0;
    format      ascii;
    class       vectorField;
    location    "1/polyMesh";
    object      points;
}

8((0 1 1) (-1 0 1) (0 1 0) (0 -1 1) (1 0 0) (0 -1 0) (-1 0 0) (1 0 1))
faces, as shown in the image:

Code:
FoamFile
{
    version     2.0;
    format      ascii;
    class       faceList;
    location    "1/polyMesh";
    object      faces;
}

8
(
4(0 2 6 1)
3(5 6 4)
3(4 6 2)
4(1 6 5 3)
4(4 7 3 5)
4(0 7 4 2)
5(3 7 4 6 1)
5(1 6 4 7 0)
)
boundaries:

Code:
FoamFile
{
    version     2.0;
    format      ascii;
    class       polyBoundaryMesh;
    location    "1/polyMesh";
    object      boundary;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

6
(
    back
    {
        type            zeroGradient;
        nFaces          1;
        startFace       0;
    }
    bottom
    {
        type            zeroGradient;
        nFaces          2;
        startFace       1;
    }
    front
    {
        type            zeroGradient;
        nFaces          1;
        startFace       3;
    }
    left
    {
        type            zeroGradient;
        nFaces          1;
        startFace       4;
    }
    right
    {
        type            zeroGradient;
        nFaces          1;
        startFace       5;
    }
    top
    {
        type            zeroGradient;
        nFaces          2;
        startFace       6;
    }
)
owner:

Code:
FoamFile
{
    version     2.0;
    format      ascii;
    class       labelList;
    note        "nPoints:8  nCells:2  nFaces:8  nInternalFaces:0";
    location    "1/polyMesh";
    object      owner;
}

8(1 0 1 0 0 1 0 1)
neighbour:

Code:
FoamFile
{
    version     2.0;
    format      ascii;
    class       labelList;
    note        "nPoints:8  nCells:2  nFaces:8  nInternalFaces:0";
    location    "1/polyMesh";
    object      neighbour;
}

0()
Again, all and complete files are at http://www.richelbilderbeek.nl/ToolO...dCorrected.htm
assert is offline   Reply With Quote

Old   March 24, 2014, 10:10
Default Solved
  #2
New Member
 
Richel Bilderbeek
Join Date: Nov 2013
Posts: 14
Rep Power: 3
assert is on a distinguished road
The source error checkMesh states that faces do not form a cell, mesh appears correct has been solved and the answer can be found there.
__________________
I will never accept a 'Friend' request, I will always accept beer
assert is offline   Reply With Quote

Reply

Tags
bent face, checkmesh, mesh, pentagon, renumbermesh

Thread Tools
Display Modes

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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
snappyhexmesh remove blockmesh geometry philipp1 OpenFOAM Running, Solving & CFD 2 December 12, 2014 11:58
fluentMeshToFoam multidomain mesh conversion problem Attesz OpenFOAM Other Meshers: ICEM, Star, Ansys, Pointwise, GridPro, Ansa, ... 12 May 2, 2013 10:52
external flow with snappyHexMesh chelvistero OpenFOAM 11 January 15, 2010 20:43
BlockMeshmergePatchPairs hjasak OpenFOAM Native Meshers: blockMesh 11 August 15, 2008 07:36
fluent add additional zones for the mesh file SSL FLUENT 2 January 26, 2008 12:55


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