<div dir="ltr"><div><div>Greetings. Thanks to everyone for such a quality software project. I have a question I have been struggling with for some time, and am not quite sure it is even possible. I want to find a way to specify the internal structure somewhat explicitly for a mesh.<br>
<br></div>As an example, I would like to embed a cylinder into a cube. The code below does this, but the cylinder and larger bounding box are not topologically connected. I instead would like a single mesh structure, but internally nodes to be densely located on the cylinder. I realize I could sort of do this with attractors, but I want to be able to quite explicitly set the internal node location. <br>
<br></div><div>If there is an example showing this, I would greatly appreciate it. All I really need is some kind of combine function for the currently generated meshes. <br><br></div><div>I saw in the examples t11.geo sort of did this with fields. Is this my best bet? I would prefer a more explicit definition than that. <br>
<br></div><div>Thanks for any ideas. <br></div><div><br>///////////////////////////////////////////////<br></div>// My attempt .geo file<br><div> <br><div><div>lc = 1e-1; // Target element size<br>R = .25;   // Magnet Radius<br>
D0 = 10;   // Top of magnet<br>D1 = 11;   // Bottom of magnet<br><br>// Total Solution Space<br>X0 = -2.;<br>X1 =  2.;<br>Y0 = -2.;<br>Y1 =  2.;<br>Z0 =  9.;<br>Z1 = 12.;<br><br>////////////////////////////////////<br>// North Pole<br>
<br>Point(0) = { 0, 0, D0, lc};<br>Point(1) = { R, 0, D0, lc};<br>Point(2) = { 0, R, D0, lc};<br>Point(3) = {-R, 0, D0, lc};<br><br>// Connect up the points<br>Circle(4) = {1, 0, 3};<br>Line(5) = {3,1};<br>Line Loop(6) = {4,5};<br>
Plane Surface(7) = {6};<br><br>////////////////////////////////////<br>// South Pole<br><br>Point(8)  = { 0,  -0.0001, D0, lc};<br>Point(9)  = { R,  -0.0001, D0, lc};<br>Point(10) = { 0, -R, D0, lc};<br>Point(11) = {-R,  -0.0001, D0, lc};<br>
<br>// Connect up the points<br>Circle(12) = {11, 8, 9};<br>Line(13) = {9, 11};<br>Line Loop(14) = {12, 13};<br>Plane Surface(15) = {14};<br><br>//////////////////////////////////////<br>// Extrude<br>Extrude {0, 0, D1-D0} { Surface{7};  }<br>
Extrude {0, 0, D1-D0} { Surface{15}; }<br><br>/////////////////////////////////////<br>// Large Bounding box<br>Point(116)  = {X0, Y0, Z0};<br>Point(117)  = {X1, Y0, Z0};<br>Point(118)  = {X1, Y1, Z0};<br>Point(119)  = {X0, Y1, Z0};<br>
Line(120) = {116,117};<br>Line(121) = {117,118};<br>Line(122) = {118,119};<br>Line(123) = {119,116};<br>Line Loop(124) = {120, 121, 122, 123}; <br>Plane Surface(125) = {124};<br>Extrude {0, 0, Z1-Z0} { Surface{125}; }<br>
<br>//////////////////////////////////////<br>// Rings of sensitivity calculation<br>//Circle(222) = {2.*1, 0, 2.*2};<br><br><br>//////////////////////////////////////<br>// Volumes<br><br>Physical Volume("plus") = {2};<br>
Physical Volume("minus") = {1};<br>Physical Volume("background") = {3};<br>Coherence;<br>Coherence;<br><br clear="all"><div>-- Trevor<br></div>
</div></div></div></div>