[Gmsh] Question about mesh refinement with stating of characteristic lengths

Christophe Geuzaine cgeuzaine at ulg.ac.be
Thu Jul 10 14:54:35 CEST 2008


Moritz Nadler wrote:
> Thanks for you fast reply but do not get what I have to do.
> 
> I tried to make a dummy surface out of the dummy point but I failed
> I would be great if somebody could give me an example of a dummy geometry for 
> changing the local mesh refinement so I can understand how I have to do it.
> 

You can use Fields to impose characteristic lengths outside the geometry.

See demos/fields.geo for an example.


> kind regards
> 
> 
> Moritz
> 
> 
> //my current .geo file with a non working dummy geometry
> highPr = 0.000015;
> midPr = 0.0002;
> lowPr = 0.0006;
> vLowPr = 0.001;
> 
> Point(1) = {0,0,0,highPr};
> Point(2) = {0.25E-2,0,0,midPr};
> Point(3) = {0.5E-2,0,0,highPr};
> Point(4) = {0.5E-2,0.025E-2,0,highPr};
> Point(5) = {0.75E-2,0.025E-2,0,midPr};
> Point(6) = {2.25E-2,0.025E-2,0,midPr};
> Point(7) = {2.5E-2,0.025E-2,0,highPr};
> Point(8) = {2.5E-2,0,0,highPr};
> Point(9) = {2.75E-2,0,0,midPr};
> Point(10) = {4E-2,0,0,lowPr};
> Point(11) = {9.5E-2,0,0,vLowPr};
> Point(13) = {0,9.525E-2,0,vLowPr};
> Point(14) = {0,4.025E-2,0,lowPr};
> Point(15) = {0,2.775E-2,0,midPr};
> Point(16) = {0,2.525E-2,0,highPr};
> Point(17) = {0,2.275E-2,0,midPr};
> Point(18) = {0,0.775E-2,0,midPr};
> Point(19) = {0,0.525E-2,0,highPr};
> Point(20) = {0,0.275E-2,0,midPr};
> Point(21) = {0,0.025E-2,0,highPr};
> 
> // dummy gemetrey to lower the mesh density in this area (does not work)
> Point(100) = {0.25E-2,0.275E-2,0,midPr};
> Line(100) = {100,100};
> Line Loop(101) = {100, 100};
> Plane Surface(102) = {101};
> 
> Line(1) = {1,2};
> Line(2) = {2,3};
> Line(3) = {3,4};
> Line(4) = {4,5};
> Line(5) = {5,6};
> Line(6) = {6,7};
> Line(7) = {7,8};
> Line(8) = {8,9};
> Line(9) = {9,10};
> Line(10) = {10,11};
> //Line(11) = {11,12};
> //Line(12) = {12,13};
> Circle(11) = {11,1,13};
> Line(13) = {13,14};
> Line(14) = {14,15};
> Line(15) = {15,16};
> Line(16) = {16,17};
> Line(17) = {17,18};
> Line(18) = {18,19};
> Line(19) = {19,20};
> Line(20) = {20,21};
> Line(21) = {21,1};
> 
> Line Loop(22) = {1,2,3,4,5,6,7,8,9,10,11,13,14,15,16,17,18,19,20,21};
> Plane Surface(23) = {22};
> 
> //end file
> 
> 
> avid Colignon schrieb:
>> Hi,
>>
>> it doesn't work because a dummy point is not enough, the dummy point 
>> must be part of a dummy line, which must be part of a dummy surface.
>>
>> Or you can have a look at:
>>
>> http://www.geuz.org/pipermail/gmsh/2008/003141.html
>>
>> Cheers,
>>
>> Dave
>>
> 
> 
> _______________________________________________
> gmsh mailing list
> gmsh at geuz.org
> http://www.geuz.org/mailman/listinfo/gmsh
> 
> 


-- 
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science
http://www.montefiore.ulg.ac.be/~geuzaine