Dear all,

I am having some trouble when trying to generate a mesh that has a threshold field around a line and a physical
field defined on any different line.

I am using gmsh 2.8.3

Here is my example:

Point(1) = {-1, -0.8, 0, 1.0};
Point(2) = {0.3, -0.8, 0, 1.0};
Point(3) = {0.3, 0.3, 0, 1.0};
Point(4) = {-1, 0.3, 0, 1.0};

Line(1) = {1, 2};
Line(2) = {2, 3};
Line(3) = {3, 4};
Line(4) = {4, 1};

Line Loop(5) = {4, 1, 2, 3};
Plane Surface(6) = {5};

Point(5) = {-0.5, -0.1, 0, 1.0};
Line(7) = {4, 5};

Field[1] = Attractor;
Field[1].EdgesList = {7};

Field[2] = Threshold;
Field[2].DistMax = 0.1;
Field[2].DistMin = 0.05;
Field[2].IField = 1;
Field[2].LcMax = 0.1;
Field[2].LcMin = 0.01;

Background Field = 2;
Physical Line(8) = {1};

Without the physical line, the mesh seems ok. With the physical line, something goes wrong. When I click "2d" in the gmsh gui, I get something that looks like the mesh I want. but after saving, the file only contains 12 or so nodes.

I tried detaching the point at the upper left corner (i.e. put another point at a slightly different position rather than using the same point for three lines). In that case, The .msh file contains more nodes, which seems ok at first, but then my FE code has trouble importing the mesh. Also, when opening the .msh file in gmsh, I don't see anything at all.

Thanks in advance for any help on this.
Best regards

_______________________________________________
gmsh mailing list
[email protected]
http://www.geuz.org/mailman/listinfo/gmsh

Reply via email to