Hi SPIS community,
I’m trying to generate a quite simple geometry in GMSH with a thin surface but I keep getting error messages when I start the simulation. More specifically I want to generate an aluminum sphere with a thin patch near it. The original idea was to apply the patch directly on the sphere (for anyone who has used it, I’m trying to reproduce the SPENVIS-EQUIPOT geometry) but I couldn’t do that without generating mesh error messages so I opted for this other configuration and to connect the two nodes via a circuit.
I’ve read the manual but I can’t figure out where I’m going wrong. I’m leaving the .geo file below:
//Variables
sph_r = 0.5;
dist = 0.05;
thick = 0.02;
hs = 0.05;
b_r = 12;
sph_m = 0.05;
sph_m_c = 0.02;
patch_m = 0.01;
b_m = 2;
/*
spr_h: sphere radius
dist: distance sphere-patch
thick: patch thickness
hs: patch half side
b_r: boundary radius
sph_m: sphere mesh dimension
sph_m_2: sphere mesh dimension of the point closest to the patch
patch_m: patch mesh dimension
b_m: boundary mesh dimension
*/
//ORIGIN
Point(1) = {0, 0, 0, 0.5};
//SPHERE - Aluminum
Point(2) = {sph_r, 0, 0, sph_m};
Point(3) = {-sph_r, 0, 0, sph_m};
Point(4) = {0, sph_r, 0, sph_m};
Point(5) = {0, -sph_r, 0, sph_m};
Point(6) = {0, 0, sph_r, sph_m_c};
Point(7) = {0, 0, -sph_r, sph_m};
Circle(1) = {7, 1, 4};
Circle(2) = {4, 1, 3};
Circle(3) = {4, 1, 6};
Circle(4) = {4, 1, 2};
Circle(5) = {2, 1, 7};
Circle(6) = {7, 1, 3};
Circle(7) = {3, 1, 6};
Circle(8) = {6, 1, 2};
Circle(9) = {3, 1, 5};
Circle(10) = {7, 1, 5};
Circle(11) = {2, 1, 5};
Circle(12) = {6, 1, 5};
Line Loop(13) = {3, 8, -4};
Ruled Surface(14) = {13};
Line Loop(15) = {4, 5, 1};
Ruled Surface(16) = {15};
Line Loop(17) = {1, 2, -6};
Ruled Surface(18) = {17};
Line Loop(19) = {2, 7, -3};
Ruled Surface(20) = {19};
Line Loop(21) = {6, 9, -10};
Ruled Surface(22) = {21};
Line Loop(23) = {7, 12, -9};
Ruled Surface(24) = {23};
Line Loop(25) = {8, 11, -12};
Ruled Surface(26) = {25};
Line Loop(27) = {11, -10, -5};
Ruled Surface(28) = {27};
Physical Surface(29) = {16, 22, 18, 20, 14, 26, 24, 28};
Surface Loop(30) = {20, 18, 16, 14, 26, 28, 22, 24};
Volume(31) = {30};
//PATCH
Point(8) = {hs, hs, sph_r+dist, patch_m};
Point(9) = {-hs, hs, sph_r+dist, patch_m};
Point(10) = {hs, -hs, sph_r+dist, patch_m};
Point(11) = {-hs, -hs, sph_r+dist, patch_m};
Point(12) = {hs, hs, sph_r+thick+dist, patch_m};
Point(13) = {-hs, hs, sph_r+thick+dist, patch_m};
Point(14) = {hs, -hs, sph_r+thick+dist, patch_m};
Point(15) = {-hs, -hs, sph_r+thick+dist, patch_m};
Point(16) = {hs, hs, sph_r+thick+thick+dist, patch_m};
Point(17) = {-hs, hs, sph_r+thick+thick+dist, patch_m};
Point(18) = {hs, -hs, sph_r+thick+thick+dist, patch_m};
Point(19) = {-hs, -hs, sph_r+thick+thick+dist, patch_m};
Line(32) = {9, 8};
Line(33) = {8, 10};
Line(34) = {10, 11};
Line(35) = {11, 9};
Line(36) = {9, 13};
Line(37) = {13, 17};
Line(38) = {8, 12};
Line(39) = {12, 16};
Line(40) = {12, 14};
Line(41) = {14, 10};
Line(42) = {15, 14};
Line(43) = {11, 15};
Line(44) = {15, 13};
Line(45) = {15, 19};
Line(46) = {19, 17};
Line(47) = {13, 12};
Line(48) = {17, 16};
Line(49) = {16, 18};
Line(50) = {18, 19};
Line(51) = {14, 18};
Line Loop(52) = {32, 33, 34, 35};
Plane Surface(53) = {52};
Line Loop(54) = {36, 47, -38, -32};
Plane Surface(55) = {54};
Line Loop(56) = {36, -44, -43, 35};
Plane Surface(57) = {56};
Line Loop(58) = {34, 43, 42, 41};
Plane Surface(59) = {58};
Line Loop(60) = {33, -41, -40, -38};
Plane Surface(61) = {60};
Line Loop(62) = {37, 48, -39, -47};
Plane Surface(63) = {62};
Line Loop(64) = {39, 49, -51, -40};
Plane Surface(65) = {64};
Line Loop(66) = {44, 37, -46, -45};
Plane Surface(67) = {66};
Line Loop(68) = {42, 51, 50, -45};
Plane Surface(69) = {68};
Line Loop(70) = {46, 48, 49, 50};
Plane Surface(71) = {70};
Coherence;
Coherence;
Line Loop(72) = {47, 40, -42, 44};
Plane Surface(73) = {72};
Physical Surface(74) = {73};
Surface Loop(75) = {53, 55, 57, 59, 61, 73};
Volume(76) = {75};
Surface Loop(77) = {71, 67, 63, 65, 69, 73};
Volume(78) = {77};
Physical Volume(79) = {76};
Physical Volume(80) = {78};
//BOUNDARY
Point(20) = {b_r, 0, 0, b_m};
Point(21) = {-b_r, 0, 0, b_m};
Point(22) = {0, b_r, 0, b_m};
Point(23) = {0, -b_r, 0, b_m};
Point(24) = {0, 0, b_r, b_m};
Point(25) = {0, 0, -b_r, b_m};
Circle(81) = {22, 1, 20};
Circle(82) = {20, 1, 24};
Circle(83) = {24, 1, 22};
Circle(84) = {20, 1, 25};
Circle(85) = {25, 1, 22};
Circle(86) = {22, 1, 21};
Circle(87) = {21, 1, 25};
Circle(88) = {21, 1, 23};
Circle(89) = {23, 1, 24};
Circle(90) = {24, 1, 21};
Circle(91) = {25, 1, 23};
Circle(92) = {23, 1, 20};
Line Loop(93) = {85, 86, 87};
Ruled Surface(94) = {93};
Line Loop(95) = {86, -90, 83};
Ruled Surface(96) = {95};
Line Loop(97) = {82, 83, 81};
Ruled Surface(98) = {97};
Line Loop(99) = {84, 85, 81};
Ruled Surface(100) = {99};
Line Loop(101) = {91, -88, 87};
Ruled Surface(102) = {101};
Line Loop(103) = {88, 89, 90};
Ruled Surface(104) = {103};
Line Loop(105) = {82, -89, 92};
Ruled Surface(106) = {105};
Line Loop(107) = {92, 84, 91};
Ruled Surface(108) = {107};
Physical Surface(109) = {102, 94, 100, 96, 104, 108, 98, 106};
Surface Loop(110) = {94, 100, 108, 106, 98, 96, 104, 102};
Surface Loop(111) = {53, 55, 57, 59, 61, 69, 65, 63, 67, 71};
Volume(112) = {30, 110, 111};
Physical Volume(113) = {112};
Thanks so much to anyone who will be able to help me.