[Gmsh] Attractor around a point in 3D

Christophe Geuzaine cgeuzaine at ulg.ac.be
Sun Jan 18 13:07:58 CET 2009


Srinath Madhavan wrote:
> Hi all,
> 
> I have a simple cube and a point at the center. Using attractors, I'm trying to refine the mesh around the point and make it coarse as we move away from the point. The code I use is pasted below:
> 
> ************CODE BEGINS*************
> 
> lc = 0.15;
> Point(1) = {-0.5, -0.5, -0.5, lc};
> Point(2) = {0.5, -0.5, -0.5, lc};
> Point(3) = {0.0, 0.0, 0.0, lc};
> Line(3) = {1, 2};
> 
> e1[] = Extrude{0, 1, 0}{Line{3};};
> e2[] = Extrude{0, 0, 1}{Surface{e1[1]};};
> 
> // Attractor field on point 3. This field returns the
> // distance to point 3.
> Field[1] = Attractor;
> Field[1].NodesList = {3};
> 
> // Threshold field using the return value of the attractor field 1 in
> // order to define a simple change in element size around the
> // attractors (i.e., around point 3)
> //
> // LcMax -                         /------------------
> //                               /
> //                             /
> //                           /
> // LcMin -o----------------/
> //        |                |       |
> //     Attractor       DistMin   DistMax
> Field[2] = Threshold;
> Field[2].IField = 1;
> Field[2].LcMin = lc / 30;
> Field[2].LcMax = lc;
> Field[2].DistMin = 0.05;
> Field[2].DistMax = 0.5;
> 
> // Use minimum of all the fields as the background field
> Field[3] = Min;
> Field[3].FieldsList = {2};
> Background Field = 3;
> 
> // Don't extend the elements sizes from the boundary inside the domain
> Mesh.CharacteristicLengthExtendFromBoundary = 0;
> 
> Physical Volume("domain") = {e2[1]};
> Physical Surface("wall") = {29,28,24,20,16,7};
> 
> ************CODE ENDS*************
> 
> The mesh generation works but I don't see any refinement around point 3. If however, I change the line:
> 
> Field[1].NodesList = {3}; TO Field[1].NodesList = {1};
> 
> Then I see refinement around point 1. This works for point 2 as well.
> 
> Can someone point out what I'm doing wrong here?


Works OK over here. Make sure you don't use the "Netgen" 3D algorithm: 
it does not support Fields (cf. the documentation).



> 
> Thanks & Best Regards,
> Srinath Madhavan
> PhD Student
> University of Alberta
> 
> Free as in Freedom. May the source be with you!
> 
> 
>       
> 
> _______________________________________________
> 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