phasefield

acurnier_2D_3.geo
View Options

			cl1 = 2;
			cl2 = -0.15;
			cl3 = -0.5;
			Point(1) = {-1., cl3, 0, cl1};
			Point(2) = {-1., 0.5, 0, cl1};
			Point(3) = {1., 0.5, 0, cl1};
			Point(4) = {1., cl3, 0, cl1};
			Point(5) = {0.7, 0.7+cl2, 0, cl1};
			Point(6) = {0.7, 1.4+cl2, 0, cl1};
			Point(7) = {0, 1.4+cl2, 0, cl1};
			Point(8) = {0, 0.7+cl2, 0, cl1};
			Point(9) = {-0.7, 1.4+cl2, 0, cl1};
			Point(10) = {-0.7, 0.7+cl2, 0, cl1};
			Line(1) = {1, 4};
			Line(2) = {4, 3};
			Line(3) = {3, 2};
			Line(4) = {2, 1};
			Line(5) = {8, 5};
			Line(6) = {5, 6};
			Line(7) = {6, 8};
			Line(8) = {7, 8};
			Line(9) = {6, 7};
			Line(10) = {10, 8};
			Line(11) = {8, 9};
			Line(12) = {7, 9};
			Line(13) = {9, 10};
			Line Loop(10) = {5, 6, 7};
			Plane Surface(11) = {10};
			Line Loop(12) = {-7, 9, 8};
			Plane Surface(13) = {12};
			Line Loop(14) = {1, 2, 3, 4};
			Plane Surface(15) = {14};
			Line Loop(16) = {10, 11, 13};
			Plane Surface(17) = {16};
			Line Loop(18) = {-8, 12, -11};
			Plane Surface(19) = {18};
			Physical Line("Contact") = {3};
			Physical Line("Top") = {9, 12};
			Physical Line("Bottom") = {1};
			Physical Surface("Top_body") = {11, 13, 17, 19};
			Physical Surface("Bottom_body") = {15};