sablier

testcases/sablier-2d/mesh.geo

+d=1e-4;
+L = 15*d;
+H = 50*d;
+lc = d;
+R=3*d;
+
+Point(1) = {0., 0., 0, lc};
+Point(2) = {L, 0., 0, lc};
+Point(3) = {L, H, 0, lc};
+Point(4) = {0.5*(L+R), H, 0, lc};
+Point(5) = {0.5*(L+R), H+d, 0, lc};
+Point(6) = {L, H+d, 0, lc};
+Point(7) = {L, H+d+H, 0, lc};
+Point(8) = {0., H+d+H, 0, lc};
+Point(9) = {0., H+d, 0, lc};
+Point(10) = {0.5*(L-R), H+d, 0, lc};
+Point(11) = {0.5*(L-R), H, 0, lc};
+Point(12) = {0., H, 0, lc};
+Line(1) = {1, 2};
+Line(2) = {2, 3};
+Line(3) = {3, 4};
+Line(4) = {4, 5};
+Line(5) = {5, 6};
+Line(6) = {6, 7};
+Line(7) = {7, 8};
+Line(8) = {8, 9};
+Line(9) = {9, 10};
+Line(10) = {10,11};
+Line(11) = {11,12};
+Line(12) = {12,1};
+
+Line Loop(1) = {1:12};
+
+Plane Surface(1) = {1};
+Physical Line("Box") = {1,2,3,4,5,6,8,9,10,11,12};
+Physical Line("Top") = {7};
+Physical Surface("Domain") = {1};
+Physical Point("PtFix") = {1};