Thank you again Jeremy, You’re solution (see below as I understand it) works. Yet, I have no clue why yours does and mine doesn’t.
Bart. ---- KanR = DefineNumber[ 0.1, Name "Parameters/KanR" ]; WanD = DefineNumber[ 0.2, Name "Parameters/WanD" ]; KanMeshsize = DefineNumber[ 0.01, Name "Parameters/KanMeshsize" ]; SetFactory("OpenCASCADE"); Rectangle(1) = {-2.5, -2.5, 0, 5, 5, 0}; Circle(5) = {0, 0, 0, KanR, 0, 2*Pi}; Curve Loop(2) = {5}; Plane Surface(2) = {2}; BooleanDifference{ Surface{1}; Delete; }{ Surface{2}; Delete; } Extrude {0, 0, WanD } { Surface{1}; } Point(15) = {0.05, 0, 0}; Point(16) = {0.05, 0, WanD}; Line(20) = {5, 15}; Line(21) = {15, 16}; Line(22) = {16, 14}; Curve Loop(8) = {21, 22, -18, 20}; Plane Surface(8) = {8}; Extrude {{0, 0, 1}, {0, 0, 0}, Pi} { Surface{8}; Layers{25}; } Field[1] = Distance; Field[1].EdgesList = {5}; Field[1].NNodesByEdge = 50; Field[2] = Threshold; Field[2].DistMax = 2.5; Field[2].DistMin = 0; Field[2].IField = 1; Field[2].LcMax = 0.25; Field[2].LcMin = KanMeshsize ; Background Field = 2; Mesh.CharacteristicLengthExtendFromBoundary = 0; Mesh.CharacteristicLengthFromPoints = 0; Mesh.CharacteristicLengthFromCurvature = 0;
_______________________________________________ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh