[Gmsh] GMSH doesn’t want to make grid inside sphere
Christophe Geuzaine
cgeuzaine at ulg.ac.be
Sat Jun 20 18:10:44 CEST 2009
Andrey wrote:
> Dear GMSH community!
>
>
> Let’s consider the following geometry:
>
> Input data to GMSH is a triangulated sphere of radius=1.
>
> Each point of triangulation has constant characteristic length.
>
>
> I pass triangulated sphere to GMSH in order to get meshed volume inside:
>
>
> gmsh-win.exe -3 boundary.geo -o mesh.mesh -format mesh -order 1 -algo
> iso -optimize
>
>
> Then I consider two different cases (input geometry is the same, the
> only difference is characteristic length):
>
>
> 1) If characteristic length=0.025 – everything is OK.
>
> See attached input: boundary_good.zip and produced: sphere_good.png
>
> Output of quality for this case:
>
> Info : Opti : 0.00 < QUAL < 0.10 : 0 elements
>
> Info : Opti : 0.10 < QUAL < 0.20 : 0 elements
>
> Info : Opti : 0.20 < QUAL < 0.30 : 8 elements
>
> Info : Opti : 0.30 < QUAL < 0.40 : 113 elements
>
> Info : Opti : 0.40 < QUAL < 0.50 : 3134 elements
>
> Info : Opti : 0.50 < QUAL < 0.60 : 112218 elements
>
> Info : Opti : 0.60 < QUAL < 0.70 : 285051 elements
>
> Info : Opti : 0.70 < QUAL < 0.80 : 389672 elements
>
> Info : Opti : 0.80 < QUAL < 0.90 : 184604 elements
>
> Info : Opti : 0.90 < QUAL < 1.00 : 15858 elements
>
> 2) If characteristic length=0.05 – produced mesh is not correct.
>
> See attached input: boundary_bad.zip and produced: sphere_bad.png
>
> Output of quality for this case:
>
> Info : Opti : 0.00 < QUAL < 0.10 : 0 elements
>
> Info : Opti : 0.10 < QUAL < 0.20 : 2707 elements
>
> Info : Opti : 0.20 < QUAL < 0.30 : 2033 elements
>
> Info : Opti : 0.30 < QUAL < 0.40 : 923 elements
>
> Info : Opti : 0.40 < QUAL < 0.50 : 649 elements
>
> Info : Opti : 0.50 < QUAL < 0.60 : 670 elements
>
> Info : Opti : 0.60 < QUAL < 0.70 : 606 elements
>
> Info : Opti : 0.70 < QUAL < 0.80 : 425 elements
>
> Info : Opti : 0.80 < QUAL < 0.90 : 174 elements
>
> Info : Opti : 0.90 < QUAL < 1.00 : 15 elements
>
>
> Summary:
>
> First case is more complex, because of larger amount of produced
> tetrahedra and because 2D-frames (of input data) are also meshed by GMSH
> (according to characteristic lengths). But first case is handled
> correctly, the second - not.
>
This is weird: I cannot reproduce the problem on my machine... Does it
still happen with recent nightly snapshots?
>
> Ant suggestions?
>
>
> Best regards,
>
> Andrey
>
>
>
> ------------------------------------------------------------------------
>
>
> ------------------------------------------------------------------------
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> gmsh mailing list
> gmsh at geuz.org
> http://www.geuz.org/mailman/listinfo/gmsh
--
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science
http://www.montefiore.ulg.ac.be/~geuzaine