geo files

test_bem/doubleHole.geo

+// Outside
+Lx_out = 2.0;
+Ly_out = 1.5;
+nx_out = 40;
+ny_out = 30;
+phi_out = 100.;
+
+// Hole 1
+Lx1 = 0.25;
+Ly1 = 0.25;
+cx1 = -0.5;
+cy1 = 0.35;
+nx1 = 8;
+ny1 = 8;
+phi1 = 215.;
+
+// Hole 2
+Lx2 = 0.25;
+Ly2 = 0.25;
+cx2 = 0.5;
+cy2 = -0.35;
+nx2 = 8;
+ny2 = 8;
+phi2 = 250.;
+
+// Outside
+Point(1) = {-Lx_out/2, -Ly_out/2, 0, 1.0};
+Point(2) = {Lx_out/2, -Ly_out/2, 0, 1.0};
+Point(3) = {Lx_out/2, Ly_out/2, 0, 1.0};
+Point(4) = {-Lx_out/2, Ly_out/2, 0, 1.0};
+
+// Hole 1
+Point(5) = {cx1 - Lx1/2, cy1 - Ly1/2, 0, 1.0};
+Point(6) = {cx1 - Lx1/2, cy1 + Ly1/2, 0, 1.0};
+Point(7) = {cx1 + Lx1/2, cy1 + Ly1/2, 0, 1.0};
+Point(8) = {cx1 + Lx1/2, cy1 - Ly1/2, 0, 1.0};
+
+// Hole 2
+Point(9) = {cx2 - Lx2/2, cy2 - Ly2/2, 0, 1.0};
+Point(10) = {cx2 - Lx2/2, cy2 + Ly2/2, 0, 1.0};
+Point(11) = {cx2 + Lx2/2, cy2 + Ly2/2, 0, 1.0};
+Point(12) = {cx2 + Lx2/2, cy2 - Ly2/2, 0, 1.0};
+
+// Outside
+Line(1) = {1, 2};
+Line(2) = {2, 3};
+Line(3) = {3, 4};
+Line(4) = {4, 1};
+
+// Hole 1
+Line(5) = {5, 6};
+Line(6) = {6, 7};
+Line(7) = {7, 8};
+Line(8) = {8, 5};
+
+// Hole 2
+Line(9) = {9, 10};
+Line(10) = {10, 11};
+Line(11) = {11, 12};
+Line(12) = {12, 9};
+
+Curve Loop(1) = {4, 1, 2, 3};
+Curve Loop(2) = {8, 5, 6, 7};
+Curve Loop(3) = {12, 9, 10, 11};
+
+Plane Surface(1) = {1, 2, 3};
+
+Transfinite Curve {1} = nx_out + 1 Using Progression 1;
+Transfinite Curve {2} = ny_out + 1 Using Progression 1;
+Transfinite Curve {3} = nx_out + 1 Using Progression 1;
+Transfinite Curve {4} = ny_out + 1 Using Progression 1;
+
+Transfinite Curve {5} = ny1 + 1 Using Progression 1;
+Transfinite Curve {6} = nx1 + 1 Using Progression 1;
+Transfinite Curve {7} = ny1 + 1 Using Progression 1;
+Transfinite Curve {8} = nx1 + 1 Using Progression 1;
+
+Transfinite Curve {9} = ny2 + 1 Using Progression 1;
+Transfinite Curve {10} = nx2 + 1 Using Progression 1;
+Transfinite Curve {11} = ny2 + 1 Using Progression 1;
+Transfinite Curve {12} = nx2 + 1 Using Progression 1;
+
+//Recombine Surface {1};
+
+Physical Curve("outside", 5) = {1, 2, 3, 4};
+Physical Curve("hole1", 6) = {5, 6, 7, 8};
+Physical Curve("hole2", 7) = {9, 10, 11, 12};
+
+SetNumber("Boundary Conditions/outside/dirichlet", phi_out);
+SetNumber("Boundary Conditions/hole1/dirichlet", phi1);
+SetNumber("Boundary Conditions/hole2/dirichlet", phi2);
\ No newline at end of file

test_bem/hole.geo

+Lx_in = 1.0;
+Ly_in = 0.1;
+cx = 0.0;
+cy = -0.6;
+
+Lx_out = 2.0;
+Ly_out = 1.5;
+
+phi_in = 200.;
+phi_out = 100.;
+
+nx_in = 40;
+ny_in = 4;
+
+nx_out = 80;
+ny_out = 60;
+
+Point(1) = {-Lx_out/2, -Ly_out/2, 0, 1.0};
+Point(2) = {Lx_out/2, -Ly_out/2, 0, 1.0};
+Point(3) = {Lx_out/2, Ly_out/2, 0, 1.0};
+Point(4) = {-Lx_out/2, Ly_out/2, 0, 1.0};
+
+
+Point(5) = {-Lx_in/2 + cx, -Ly_in/2 + cy, 0, 1.0};
+Point(6) = {-Lx_in/2 + cx, Ly_in/2 + cy, 0, 1.0};
+Point(7) = {Lx_in/2 + cx, Ly_in/2 + cy, 0, 1.0};
+Point(8) = {Lx_in/2 + cx, -Ly_in/2 + cy, 0, 1.0};
+
+Line(1) = {1, 2};
+Line(2) = {2, 3};
+Line(3) = {3, 4};
+Line(4) = {4, 1};
+
+Line(5) = {5, 6};
+Line(6) = {6, 7};
+Line(7) = {7, 8};
+Line(8) = {8, 5};
+
+Curve Loop(1) = {4, 1, 2, 3};
+Curve Loop(2) = {8, 5, 6, 7};
+
+Plane Surface(1) = {1, 2};
+
+Transfinite Curve {1} = nx_out + 1 Using Progression 1;
+Transfinite Curve {2} = ny_out + 1 Using Progression 1;
+Transfinite Curve {3} = nx_out + 1 Using Progression 1;
+Transfinite Curve {4} = ny_out + 1 Using Progression 1;
+
+Transfinite Curve {5} = ny_in + 1 Using Progression 1;
+Transfinite Curve {6} = nx_in + 1 Using Progression 1;
+Transfinite Curve {7} = ny_in + 1 Using Progression 1;
+Transfinite Curve {8} = nx_in + 1 Using Progression 1;
+
+//Recombine Surface {1};
+
+Physical Curve("inside", 5) = {5, 6, 7, 8};
+Physical Curve("outside", 6) = {1, 2, 3, 4};
+
+SetNumber("Boundary Conditions/inside/dirichlet", phi_in);
+SetNumber("Boundary Conditions/outside/dirichlet", phi_out);
\ No newline at end of file

test_bem/impossible.geo

+L = 1.0;
+
+phi_bottom = 200.;
+phi_sides = 100.;
+
+n = 100;
+
+Point(1) = {0, 0, 0, 1.0};
+Point(2) = {L, 0, 0, 1.0};
+Point(3) = {L, L, 0, 1.0};
+Point(4) = {0, L, 0, 1.0};
+Line(1) = {1, 2};
+Line(2) = {2, 3};
+Line(3) = {3, 4};
+Line(4) = {4, 1};
+Curve Loop(1) = {4, 1, 2, 3};
+Plane Surface(1) = {1};
+Transfinite Curve {3, 1} = n+1 Using Progression 1;
+Transfinite Curve {4, 2} = n+1 Using Progression 1;
+Transfinite Surface {1};
+
+Recombine Surface {1}; // quads instead of triangles
+
+Physical Curve("left_edge", 5) = {4};
+Physical Curve("bottom_edge", 6) = {1};
+Physical Curve("right_edge", 7) = {2};
+Physical Curve("top_edge", 8) = {3};
+
+// additional parameters given to the solver
+SetNumber("Boundary Conditions/left_edge/dirichlet", phi_sides);
+SetNumber("Boundary Conditions/bottom_edge/dirichlet", phi_bottom);
+SetNumber("Boundary Conditions/right_edge/dirichlet", phi_sides);
+SetNumber("Boundary Conditions/top_edge/neumann", 0.);
\ No newline at end of file

test_bem/l_shape.geo

+w_inlet = 0.5;
+w_outlet = 1.0;
+L_inlet = 1.0;
+L_outlet = 1.5;
+
+phi_inlet = 200.;
+phi_outlet = 100.;
+
+n_trans_inlet = 10;
+n_long_inlet = 10;
+n_trans_outlet = 10;
+n_long_outlet = 15;
+
+Point(1) = {0, 0, 0, 1.0};
+Point(2) = {w_inlet, 0, 0, 1.0};
+Point(3) = {w_inlet + L_outlet, 0, 0, 1.0};
+Point(4) = {w_inlet + L_outlet, w_outlet, 0, 1.0};
+Point(5) = {w_inlet, w_outlet, 0, 1.0};
+Point(6) = {w_inlet, w_outlet + L_inlet, 0, 1.0};
+Point(7) = {0, w_outlet + L_inlet, 0, 1.0};
+Point(8) = {0, w_outlet, 0, 1.0};
+
+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, 1};
+
+Curve Loop(1) = {8, 1, 2, 3, 4, 5, 6, 7};
+
+Plane Surface(1) = {1};
+
+Transfinite Curve {1} = n_trans_inlet + 1 Using Progression 1;
+Transfinite Curve {2} = n_long_outlet + 1 Using Progression 1;
+Transfinite Curve {3} = n_trans_outlet + 1 Using Progression 1;
+Transfinite Curve {4} = n_long_outlet + 1 Using Progression 1;
+Transfinite Curve {5} = n_long_inlet + 1 Using Progression 1;
+Transfinite Curve {6} = n_trans_inlet + 1 Using Progression 1;
+Transfinite Curve {7} = n_long_inlet + 1 Using Progression 1;
+Transfinite Curve {8} = n_trans_outlet + 1 Using Progression 1;
+
+//Recombine Surface {1};
+
+Physical Curve("sides", 5) = {1, 2, 4, 5, 7, 8};
+Physical Curve("inlet", 6) = {3};
+Physical Curve("outlet", 7) = {6};
+
+SetNumber("Boundary Conditions/sides/neumann", 0.);
+SetNumber("Boundary Conditions/inlet/dirichlet", phi_inlet);
+SetNumber("Boundary Conditions/outlet/dirichlet", phi_outlet);
\ No newline at end of file

test_bem/rectangle.geo

+Lx = 1.5;
+Ly = 1.0;
+
+phi_left = 200.;
+phi_right = 100.;
+
+nx = 20;
+ny = 5;
+
+Point(1) = {0, 0, 0, 1.0};
+Point(2) = {Lx, 0, 0, 1.0};
+Point(3) = {Lx, Ly, 0, 1.0};
+Point(4) = {0, Ly, 0, 1.0};
+Line(1) = {1, 2};
+Line(2) = {2, 3};
+Line(3) = {3, 4};
+Line(4) = {4, 1};
+Curve Loop(1) = {4, 1, 2, 3};
+Plane Surface(1) = {1};
+Transfinite Curve {3, 1} = nx+1 Using Progression 1;
+Transfinite Curve {4, 2} = ny+1 Using Progression 1;
+Transfinite Surface {1};
+
+Recombine Surface {1}; // quads instead of triangles
+
+Physical Curve("left_edge", 5) = {4};
+Physical Curve("bottom_edge", 6) = {1};
+Physical Curve("right_edge", 7) = {2};
+Physical Curve("top_edge", 8) = {3};
+
+// additional parameters given to the solver
+SetNumber("Boundary Conditions/left_edge/dirichlet", phi_left);
+SetNumber("Boundary Conditions/bottom_edge/neumann", 0.);
+SetNumber("Boundary Conditions/right_edge/dirichlet", phi_right);
+SetNumber("Boundary Conditions/top_edge/neumann", 0.);