s40571…

mesh_dcb_2d.geo
View Options

			dx = 0.0025;

			Point(1) = {0,0,0,dx};
			Point(2) = {0,0.05,0,dx};
			Point(3) = {0,-0.05,0,dx};
			Point(4) = {1,0,0,dx};
			Point(5) = {1,0.05,0,dx};
			Point(6) = {1,-0.05,0,dx};
			Line(1) = {1, 2};
			Line(2) = {2, 5};
			Line(3) = {5, 4};
			Line(4) = {1, 4};
			Line(5) = {1, 3};
			Line(6) = {6, 4};
			Line(7) = {3, 6};
			Line Loop(8) = {2, 3, -4, 1};
			Plane Surface(9) = {-8};
			Line Loop(10) = {5, 7, 6, -4};
			Plane Surface(11) = {10};
			Physical Surface("bulk") = {9,11};
			Physical Line("coh") = {4};

			Transfinite Surface "*";
			Recombine Surface "*";

			Mesh.SecondOrderIncomplete = 1;