[Gmsh] structured and unstructured meshes bug

Christophe Geuzaine cag32 at case.edu
Wed Mar 22 02:11:12 CET 2006


Michel Aubourg wrote:
> Re-bonjour,
> 
> Avec le fichier joint cela ira mieux.
> 

Hi Michel - Unfortunately, that's inherent to the Delaunay technique as 
implemented in the default "isotropic" algorithm. You can use the 
combined Delaunay+advancing front technique from the "netgen" algorithm 
to fix this. (Add e.g. "Mesh.Algorithm3D = 4;" at the beginning of your 
file; beware that your volume boundary has to be defined with exterior 
pointing normals when using this algorithm.)

Best,

Christophe

> Cordialement
> 	Michel Aubourg
> 
> 
> ------------------------------------------------------------------------
> 
> wBoite = 8;
> lBoite = 12;
> hBoite = 4;
> 
> hSubstrat = 0.635;
> 
> wStrip = 0.63;
> lStrip = 7.2;
> 
> lPlot = 1.5;
> 
> lGap = ( lBoite - lStrip ) / 2 - lPlot;
> 
> mSurface = 0.999999 * wStrip / 2;
> mVolume = 0.999999 * 2;
> 
> Point(1) = {0,0,0,mSurface};
> Extrude {wStrip/2,0,0} {
>   Point{1};
>   Layers {{wStrip/2/mSurface},{1.0}};
> }
> Extrude {(wBoite-wStrip)/2,0,0} {
>   Point{2};
>   Layers {{(wBoite-wStrip)/2/mSurface},{1.0}};
> }
> Extrude {0,lPlot,0} {
>   Line{1,2};
>   Layers {{lPlot/mSurface},{1.0}};
> }
> Extrude {0,lGap,0} {
>   Line{3,7};
>   Layers {{lGap/mSurface},{1.0}};
> }
> Extrude {0,lStrip/2,0} {
>   Line{11,15};
>   Layers {{lStrip/2/mSurface},{1.0}};
> }
> 
> volSubstrat[] =
> Extrude {0,0,-hSubstrat} {
>   Surface{6,10,14,18,22,26};
>   Layers {{hSubstrat/mSurface},{1.0}};
> }
> ;
> 
> Point(64) = {0,0,hBoite,mVolume};
> 
> Extrude {wBoite/2,0,0} {
>   Point{64};
>   Layers {{wBoite/2/mVolume},{1.0}};
> }
> Extrude {0,lBoite/2,0} {
>   Line{159};
>   Layers {{lBoite/2/mVolume},{1.0}};
> }
> 
> Line(164) = {1,64};
> Line(165) = {65,3};
> Line(166) = {15,67};
> Line(167) = {66,12};
> 
> Line Loop(168) = {9,17,25,166,-162,165};
> Plane Surface(169) = {168};
> Line Loop(170) = {23,166,-160,167,19};
> Plane Surface(171) = {170};
> Line Loop(172) = {167,-20,-12,-4,164,161};
> Plane Surface(173) = {172};
> Line Loop(174) = {1,2,-165,-159,-164};
> Plane Surface(175) = {174};
> Surface Loop(176) = {163,175,-6,-10,169,-18,-26,-171,173,-22,-14};
> 
> Volume(177) = {176};
> 
> // ligne
> Physical Line(178) = {33};
> 
> // zMin
> Physical Surface(179) = {70,114,158,136,48,92};
> // plot
> Physical Surface(180) = {6};
> // gap
> Physical Surface(181) = {14};
> // strip
> Physical Surface(182) = {22};
> // substrat
> Physical Surface(183) = {10,18,26};
> // zMax
> Physical Surface(184) = {163};
> 
> // yMin
> Physical Surface(185) = {57,175,35};
> // yMax
> Physical Surface(186) = {153,171,131};
> 
> // xMin
> Physical Surface(187) = {47,91,135,173};
> // xMax
> Physical Surface(188) = {61,105,149,169};
> 
> // substrat
> Physical Volume(189) = {volSubstrat[]};
> // air
> Physical Volume(190) = {177};
> 
> 
> ------------------------------------------------------------------------
> 
> _______________________________________________
> gmsh mailing list
> gmsh at geuz.org
> http://www.geuz.org/mailman/listinfo/gmsh


-- 
Christophe Geuzaine
Assistant Professor, Case Western Reserve University, Mathematics
http://www.case.edu/artsci/math/geuzaine