[Gmsh] 2D unstructured mesh problem

Umberto Paoletti paoletti at circuit.kuee.kyoto-u.ac.jp
Mon Nov 10 11:49:21 CET 2008


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



-- 
---------------------------------------------------------------------------------------
Umberto Paoletti, Dr. Eng.,
Building A1, Room 402,
Graduate School of Engineering,
Department of Electrical Engineering,
Kyoto University, Nishikyo-ku,
Katsura, 615-8510, Kyoto, Japan
Tel.: +81-(0)75-383-2246
Fax.: +81-(0)75-383-2245
E-mail: paoletti at circuit.kuee.kyoto-u.ac.jp
---------------------------------------------------------------------------------------