[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};
> 
> 
>