[Gmsh] 3D mesh crash. Coplanar points warning

[Benjamin SAUTHIER]

Thx for your answer. I modified the geo file to suit
the requirement.

gmsh seems to mesh (attached the end of the log,
debian gmsh 1.60-1) with no error and no warnings
anymore, but it quits silently or crash in gui mode
without saving the mesh, after having swapped edges.

Either there is something wrong in my mesh or the mesh
is too big (I expect 1-2 millions points, and I mesh
with 512Mb Ram on AMD 2.4GHz) but I can't figure out.

I attached the geo files (only the upper part of the
domain is meshed at the moment). This is for CFD
calculation using OpenFoam.

Thx

Benjamin


Info    : Nod=13100/26209 Elm=70689
Info    : Vol=20957.5 (8312 4707 81)
Info    : Nod=14410/26209 Elm=76421
Info    : Vol=20957.5 (8758 5481 171)
Info    : Nod=15720/26209 Elm=81932
Info    : Vol=20957.5 (9226 6215 279)
Info    : Nod=17030/26209 Elm=87501
Info    : Vol=20957.5 (9698 6922 410)
Info    : Nod=18340/26209 Elm=94418
Info    : Vol=20957.5 (10201 7596 543)
Info    : Nod=19650/26209 Elm=102515
Info    : Vol=20957.5 (10513 8491 646)
Info    : Nod=20960/26209 Elm=109927
Info    : Vol=20957.5 (10771 9451 738)
Info    : Nod=22270/26209 Elm=117038
Info    : Vol=20957.5 (11031 10425 814)
Info    : Nod=23580/26209 Elm=124157
Info    : Vol=20957.5 (11292 11376 912)
Info    : Nod=24890/26209 Elm=131237
Info    : Vol=20957.5 (11584 12289 1017)
Info    : Nod=26200/26209 Elm=138127
Info    : Vol=20957.5 (12077 13082 1041)
Info    : Swapped 2311 edges
Info    : Swapped 1624 edges
Info    : Swapped 1091 edges
Info    : Swapped 677 edges
Info    : Swapped 299 edges
Info    : Swapped 182 edges
Info    : Swapped 104 edges
Info    : Swapped 53 edges
Info    : Swapped 19 edges
Complété

Benjamin SAUTHIER wrote:
> Hello,
> 
> I want to mesh a domain bounded by a 2m long
cylinder
> inside an approximately 20*20*25m box.
> 
> The surface mesh of the cylinder is generated by
> extrusion of a circle (Transfinite mesh) to obtain
> quadrangles.
> 
> The characteristic length is max 0.5m at the
boundary
> of the box and min 0.06m on the cylinder.
> 
> When generating the 3D mesh I got the  warnings of
> coplanar points and "Unable to add points"
statements.
> Ultimately gmsh crashes.

Benjamin - Currently gmsh cannot generate 3D
tetrahedral meshes if the 
surface mesh contains quadrangles: we would need to
introduce pyramids 
to connect the tets with the quads, and that's not
implemented.

Christophe


	
	
		
___________________________________________________________________________
Appel audio GRATUIT partout dans le monde avec le nouveau Yahoo! Messenger 
Téléchargez cette version sur http://fr.messenger.yahoo.com
-------------- next part --------------
A non-text attachment was scrubbed...
Name: diffuseur5_norecombine.geo
Type: application/octet-stream
Size: 14994 bytes
Desc: 3225451133-diffuseur5_norecombine.geo
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20051001/7fe41cf9/attachment.geo>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: ring2_ring.geo
Type: application/octet-stream
Size: 4414 bytes
Desc: 1536246974-ring2_ring.geo
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20051001/7fe41cf9/attachment-0001.geo>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: ring2_ringwrap.geo
Type: application/octet-stream
Size: 12803 bytes
Desc: 1820655102-ring2_ringwrap.geo
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20051001/7fe41cf9/attachment-0002.geo>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: functions.geo
Type: application/octet-stream
Size: 3079 bytes
Desc: 256550975-functions.geo
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20051001/7fe41cf9/attachment-0003.geo>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: units.geo
Type: application/octet-stream
Size: 37 bytes
Desc: 3944362996-units.geo
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20051001/7fe41cf9/attachment-0004.geo>