[Gmsh] Problem with thresholdfield in combination with physical lines

Sebastian Otten sebi-otten at gmx.de
Tue Aug 11 03:02:49 CEST 2015


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