[Gmsh] AttractorAnisoCurve

[Sathyanarayan Rao]

Hello gmsh team,


I am using attractor field to create finer elements around a line. However at end points of the curve (point 5 and 47) , there is no refinement around the line:


How to include these points to have a refinement around the line even at end points ? Please suggest me to achieve this ! There is no option of Nodeslist in


AnisoCurve attractor.



// mesh size description

cl_1   =  0.5;
cl_2   =  0.5;

// boundary points that forms Rhizotron
Point(1) = {-15, -40, 0, cl_1};
Point(2) = {15, -40, 0, cl_1};
Point(3) = {15, 5.000000e-01, 0, cl_1};
Point(4) = {-15, 5.000000e-01, 0, cl_1};


// lines that connect boundary
Line(1) = {1, 2};
Line(2) = {3, 2};
Line(3) = {4, 3};
Line(4) = {1, 4};
Line Loop(4) = {1, -2, -3, -4};
Plane Surface(1) = {4};


// no. of spline surfaces = 59

X1 ={0.0083964,-0.093494,-0.17967,-0.32248,-0.41083,-0.49924,-0.58584,-0.90283,-1.0525,-1.1341,-0.89573,-0.7397,-0.56147,-0.51866,-0.50047,-0.44851,-0.40311,-0.37085,-0.35881,-0.17094,-0.09302,-0.089766,-0.15064,-0.1671,-0.17991,-0.51334,-0.79397,-0.99539,-1.0278,-1.0591,-1.0959,-0.87662,-0.66386,-0.41572,-0.15319,0.064832,0.40769,0.73987,1.0809,1.3972,1.6833,2.0053,2.1812,2.1812};
Y1 ={0.0083964,-3.985,-4.9787,-5.7556,-6.5471,-7.3393,-8.1315,-8.866,-9.6516,-10.4436,-11.2023,-11.9829,-12.7599,-13.5586,-14.358,-15.1562,-15.9547,-16.7522,-17.5515,-18.3217,-19.0864,-19.8853,-20.6824,-21.4818,-22.2814,-22.9989,-23.7322,-24.5061,-25.3042,-26.1024,-26.897,-27.6281,-28.3989,-29.1507,-29.8971,-30.6644,-31.3839,-32.1038,-32.8173,-33.5493,-34.2945,-35.0217,-35.8015,-35.8015};


// Define spline surfaces

LN = 91;
nR = #X1[ ];
p0 = newp; Point(p0) =  {X1[0],Y1[0],0,cl_2};
Point{p0} In Surface {1};
For i In {1:nR-2}
p0 = newp; Point(p0)= {X1[i], Y1[i], 0, cl_2};
l0 = newl; Line(l0)  =    {newp-2,newp-1};
Line{l0} In Surface {1};
Point{p0} In Surface {1};
EndFor


NL = newl-1; NP = newp-1;
Line{5:NL} In Surface {1};
Point{5:NP} In Surface {1};
Physical Curve(90) = {5:NL};
Physical Line(10000) = {1}; Physical Line(10001) = {3}; Physical Line(10002) = {2}; Physical Line(10003) = {4}; Physical Surface(100) = {1};

Mesh.CharacteristicLengthExtendFromBoundary = 0;
Field[1] = AttractorAnisoCurve;
Field[1].EdgesList = {5:NL};
Field[1].NNodesByEdge = 20;
Field[1].dMax = 15*cl_2;
Field[1].dMin = cl_2;
Field[1].lMaxNormal = cl_2;
Field[1].lMaxTangent = cl_2;
Field[1].lMinNormal = cl_2;
Field[1].lMinTangent = cl_2;
Background Field = 1;
Mesh.CharacteristicLengthMax = cl_1;



________________________________
Best Regards,

Sathyanarayan Rao, PhD student
Earth and Life Institute/Environmental Sciences (ELI-e)
Université catholique de Louvain
c.037, Croix du Sud 2, Louvain-la-Neuve, Belgium
Phone: 010473827 ( intercom 73827)

________________________________
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://onelab.info/pipermail/gmsh/attachments/20180313/fd53b495/attachment.html>