[Time4Tea]

Hi, I am trying to mesh a 2D geometry using gmsh, with a boundary layer mesh applied to two of the edges. However, it seems to be failing, with results such as as following:







The .geo file I am using is as follows:

element_size = 2.0000;
Point(1) = {0.7500, 0.0000, 0.0000, element_size};
Point(2) = {37.5000, 0.0000, 0.0000, element_size};
Point(3) = {37.5000, 40.0000, 0.0000, element_size};
Point(4) = {4.9, 40.0000, 0.0000, 0.05000};
Point(5) = {0.6, 25.448, 0.0000, 0.05000};
Point(6) = {0, 25.625, 0.0000, 0.05000};
Point(7) = {0.0000, 25.0000, 0.0000, 0.05000};
Point(8) = {0.0000, 24.7500, 0.0000, 0.2500};
Point(9) = {0.0000, 0.0000, 0.0000, 0.2500};

Line(10) = {1,2};
Line(11) = {2,3};
Line(12) = {3,4};
Line(13) = {4,5};
Circle(14) = {5,6,7};
Line(15) = {7,8};
Line(16) = {8,9};
Line(17) = {9,1};

Curve Loop(1) = {10,11,12,13,14,15,16,17};

Plane Surface(1) = {1};

//Boundary Layer
Field[1] = BoundaryLayer;
Field[1].AnisoMax = 1000;
Field[1].Quads = 1;
Field[1].Thickness = 0.2;
Field[1].CurvesList = {13,14};
Field[1].Ratio = 1.4;
Field[1].Size = 0.02;
Field[1].SizeFar = 0.2;
Field[1].IntersectMetrics = 1;

BoundaryLayer Field = 1;

Mesh.Algorithm = 6;
Recombine Surface {1};

Does anyone have any ideas why it is not working? I have tried playing with many of the boundary layer and sizing settings, but unfortunately nothing seems to work.