[Gmsh] Problems to refine near a vertex of a L-shaped domain

Durufle marc marc_durufle at yahoo.fr
Tue Dec 5 13:13:04 CET 2006


I tried to mesh a L-shaped domain. And I want to
refine locally on the vertex where the solution will
be singular. I observed that all is fine if the
refinement is large enough. If the parameter of
refinement is small (0.001 here), the produced mesh is
ugly.
I attached the details of the mesh when you zoom on
the refined vertex.

The geo file is the following :

lc = 0.1;

a = 1.0;
b = 1.0;
h = 1.0;

refine = 0.001;

Point(1) = {0, 0, 0, lc*refine};
Point(2) = {a, 0, 0, lc} ;
Point(3) = {0, b, 0, lc} ;
Point(4) = {a, -h, 0, lc} ;
Point(5) = {-h, b, 0, lc} ;
Point(6) = {-h, -h, 0, lc} ;

Line(1) = {1,3} ;
Line(2) = {3,5} ;
Line(3) = {5,6} ;
Line(4) = {6,4} ;
Line(5) = {4,2} ;
Line(6) = {2,1} ;

Line Loop(1) = {1,2,3,4,5,6} ;

// Ruled Surface(1) = {1} ;
Plane Surface(1) = {1};

Physical Line(1) = {1,2,3,4,5,6} ;

MySurface = 1;
Physical Surface(MySurface) = {1} ;

Do you have an idea to fix this problem ?


	

	
		
___________________________________________________________________________ 
Découvrez une nouvelle façon d'obtenir des réponses à toutes vos questions ! 
Profitez des connaissances, des opinions et des expériences des internautes sur Yahoo! Questions/Réponses 
http://fr.answers.yahoo.com
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Lshaped1.jpg
Type: image/jpeg
Size: 95443 bytes
Desc: 479974726-Lshaped1.jpg
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20061205/a4e713bd/attachment.jpg>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Lshaped2.jpg
Type: image/jpeg
Size: 101569 bytes
Desc: 3298822053-Lshaped2.jpg
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20061205/a4e713bd/attachment-0001.jpg>