[Gmsh] 2D unstructured mesh problem
David Colignon
David.Colignon at ulg.ac.be
Mon Nov 10 12:42:57 CET 2008
Hi,
your Line and Line Loop definitions are wrong.
You should have a look at the tutorial t1.geo
Cheers,
Dave
--
David Colignon, Ph.D.
Collaborateur Logistique F.R.S.-FNRS
CÉCI - Consortium des Équipements de Calcul Intensif
ACE - Applied & Computational Electromagnetics
Institut Montefiore B28
Université de Liège
4000 Liège - BELGIQUE
Tél: +32 (0)4 366 37 32
Fax: +32 (0)4 366 29 10
WWW: http://hpc.montefiore.ulg.ac.be/
Agenda: http://www.google.com/calendar/embed?src=david.colignon%40gmail.com
Umberto Paoletti wrote:
> Hello,
>
> I found a problem in this very simple example.
> Some corners of the internal hole seem to be skipped by the mesh algorithm.
> Why?
>
> Thanks
> Umberto
>
>
> unit = 0.001000;
>
> lc = 0.1*unit;
>
> edge = 2*unit;
> holedge = 0.5*unit;
>
> z = 0*unit;
>
> //external edge
> x1 = 0*unit;
> y1 = 0*unit;
> x2 = x1+edge;
> y2 = y1;
> x3 = x2;
> y3 = y2+edge;
> x4 = x3-edge;
> y4 = y3;
> Point(1) = {x1,y1,z, 1*lc};
> Point(2) = {x2,y2,z, 1*lc};
> Point(3) = {x3,y3,z, 1*lc};
> Point(4) = {x4,y4,z, 1*lc};
> Line(1) = {1:4,1};
> Line Loop(2) = {1};
>
> x5 = 1*unit;
> y5 = 1*unit;
> x6 = x5+holedge;
> y6 = y5;
> x7 = x6;
> y7 = y6+holedge;
> x8 = x7-holedge;
> y8 = y7;
>
> //hole edge
> Point(5) = {x5,y5,z, 1*lc};
> Point(6) = {x6,y6,z, 1*lc};
> Point(7) = {x7,y7,z, 1*lc};
> Point(8) = {x8,y8,z, 1*lc};
> Line(3) = {5:8,5};
> Line Loop(4) = {3};
>
>
> Plane Surface(1) = {2,4};
>
>
>