[Gmsh] Unable to create 3D mesh

Christophe Geuzaine cgeuzaine at ulg.ac.be
Tue Feb 28 22:33:10 CET 2017


See attached corrected file and comments therein. Basically:

1) embedded lines/surfaces are not a substitute for correct topology! they are only provided to embed 0D/1D/2D meshes in 1D/2D/3D ones (not e.g. 1D in 1D, or 2D in 2D, or 3D in 3D)

2) to do the intersections automatically you will want to use the new SetFactory("OpenCASCADE") mode in recent nightly builds, and the "BooleanFragments" operation. See gmsh/demos/boolean for examples (there is no official documentation yet).


-------------- next part --------------
A non-text attachment was scrubbed...
Name: fixed.geo
Type: application/octet-stream
Size: 3519 bytes
Desc: not available
URL: <http://onelab.info/pipermail/gmsh/attachments/20170228/7326bc1e/attachment.geo>
-------------- next part --------------


> On 28 Feb 2017, at 18:54, adhitomo sulistyo <adhitomo1 at gmail.com> wrote:
> 
> Hi, I am trying to create a reservoir model with unstructured grids. The reservoir is box shaped, with two horizontal wells on one side, and one vertical well on the other side. The model is basically a box with three cylinders inside.
> 
> Initially I succeeded in creating the 3D mesh. However, later I found that the mesh from the cylinder and the box were different and intersecting with each other. Based on the information from the web, I need to embed the surfaces inside the volume.
> 
> After successfully embedding the surfaces, the 2D mesh looks good. There are no longer intersecting grids However, when I try to run the 3D mesh, it fails. I get the following error message:
> 
> Info    : Done meshing 2D (0.715054 s)
> Info    : Meshing 3D...
> Info    : Delaunay Meshing 4 volumes with 4 connected components
> Info    : Meshing volume 227 (Delaunay)
> Info    : 7174 points 36783 edges and 14324 faces in the initial mesh
> Info    : Found void region
> Info    : Found void region
> Info    : Found void region
> Info    : Found void region
> Error   : No tetrahedra in region 227 0
> Info    : 3D point insertion terminated (7204 points created):
> Info    :  - 0 Delaunay cavities modified for star shapeness
> Info    :  - 0 points could not be inserted
> Info    :  - 0 tetrahedra created in 0.010615 sec. (0 tets/sec.)
> Info    : Meshing volume 228 (Delaunay)
> Info    : 155 points 592 edges and 280 faces in the initial mesh
> Info    : Found region 228
> Info    : 0 points created - Worst tet radius is 1.10375 (PTS removed 0 0)
> Info    : 3D point insertion terminated (7185 points created):
> Info    :  - 0 Delaunay cavities modified for star shapeness
> Info    :  - 0 points could not be inserted
> Info    :  - 308 tetrahedra created in 0.004804 sec. (64113 tets/sec.)
> Info    : Meshing volume 229 (Delaunay)
> 
> I have also attached the gmsh file. 
> 
> Thank you for your help.
> 
> Kind regards,
> Adhitomo
> 
> 
> <test2>_______________________________________________
> gmsh mailing list
> gmsh at onelab.info
> http://onelab.info/mailman/listinfo/gmsh

-- 
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science 
http://www.montefiore.ulg.ac.be/~geuzaine

Free software: http://gmsh.info | http://getdp.info | http://onelab.info



More information about the gmsh mailing list