[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