peculiar issue with opencascade : lines embedded in Volume conflict with non intersecting surface
Hi everyone, I have come up with a peculiar issue with a simple hexahedral mesh creation in GMSH GUI. Specifically I am including a surface and two lines in a box volume (model simplified from original to locate the issue), and I use OPENCASCADE to configure any intersecting components and the embedment of lower dimension entities. The problem is that the line behind the surface is eventually removed from the volume when I include the surface in the volume through the booleanfragment command (see image bellow). If I move the line so that its projection to the surface plane is outside the surface it works ok. In any case gmsh does not give any errors.
thank in advance for any help!
Here is the simple .geo file:
SetFactory("OpenCASCADE");
Point(1) = {30, 80, 0, 400}; Point(2) = {30, 80, -16.5, 400}; Point(3) = {30, 0, 0, 400}; Point(4) = {30, 0, -16.5, 400};
Point(5) = {-6.9, 1.5, -4, 400}; Point(6) = {-14.45, 1.5, -7, 400}; Point(7) = {36, 79.5, -4, 400}; Point(8) = {44, 79.5, -7, 400};
Line(1) = {1, 2}; Line(2) = {1, 3}; Line(3) = {3, 4}; Line(4) = {4, 2};
Line(5) = {5, 6}; Line(6) = {7, 8};
//This surface is the issue in combination with Line(6)!!!!!!! Line Loop(7) = {3,4,-1,2}; Plane Surface(1) = {7} ;
//Create box volume: Box(1) = {-50, -50, -30, 130, 180, 30};
v1() = BooleanFragments {Volume{:}; Delete; }{ Curve{:}; Delete; }; v2() = BooleanFragments {Volume{:}; Delete; }{ Surface{:}; Delete; };
Mesh.MshFileVersion = 2.2;