[Gmsh] Hole in the geometry after meshing

Vadivelan C vc at tetrahedrix.com
Mon Jun 15 08:17:52 CEST 2009


Hi

I created a geometry using salome (did surface mesh and exported in stl
format). Then i imported to gmsh and seen the surface mesh which looks fine.
I generated 3d mesh and converted into parafoam format. When i seen the
sufrace in parafoam i found there was some holes. During 3d mesh generation
in gmsh, i got a warning like

Warning:  Polygon 1 has two identical verts in facet 14512.
Warning:  Polygon 1 has two identical verts in facet 14612.
Warning:  Polygon 1 has two identical verts in facet 15192.
Warning:  Polygon 1 has two identical verts in facet 15212.


Info    : Meshing volume 31 (Tetgen+Delaunay)
Info    : 6653 points 6491 edges and 13298 faces in the final mesh
Info    : 6653 points created -- Worst tet radius is 44.4765
Info    : 6663 points created -- Worst tet radius is 10.9067
Info    : 6683 points created -- Worst tet radius is 12.0333


Simillarly i got some warnings while converting gmsh to openfoam format .

--> FOAM Warning :
    From function polyMesh::polyMesh(... construct from shapes...)
    in file meshes/polyMesh/polyMeshFromShapeMesh.C at line 576
    Found 13222 undefined faces in mesh; adding to default patch.
Finding faces of patch 0
--> FOAM Warning :
    From function gmshToFoam
    in file gmshToFoam.C at line 891
    Could not match gmsh face 3(999 2356 535) to any of the interior or
exterior faces that share the same 0th point
--> FOAM Warning :
    From function gmshToFoam
    in file gmshToFoam.C at line 891
    Could not match gmsh face 3(467 2356 1687) to any of the interior or
exterior faces that share the same 0th point
--> FOAM Warning : Not using gmsh face 3(2395 2095 1415) since zero vertex
is not on boundary of polyMesh
--> FOAM Warning :
    From function gmshToFoam
    in file gmshToFoam.C at line 891
    Could not match gmsh face 3(2395 2095 1415) to any of the interior or
exterior faces that share the same 0th point
--> FOAM Warning : Not using gmsh face 3(2395 1562 2466) since zero vertex
is not on boundary of polyMesh
--> FOAM Warning :
    From function gmshToFoam
    in file gmshToFoam.C at line 891
    Could not match gmsh face 3(2395 1562 2466) to any of the interior or
exterior faces that share the same 0th point


I have attached those files and snapshots of the geometry. I dont know from
where these holes comes from.

Is it problem in salome or gmsh or openfoam ?

Thanks
 - Ranjith
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20090615/4852e5a2/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: gmsh.png
Type: image/png
Size: 86507 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20090615/4852e5a2/attachment.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: parafoam.png
Type: image/png
Size: 78527 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20090615/4852e5a2/attachment-0001.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: cylinder.geo
Type: application/octet-stream
Size: 2617 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20090615/4852e5a2/attachment.geo>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: cylindercone.tgz
Type: application/x-gzip
Size: 235837 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20090615/4852e5a2/attachment.tgz>