<div dir="ltr"><div><div>I had not seen tutorial t.15 which includes the 'surface in volume' command. This is almost exactly what I want! Except...<br> <br>After mesh generation, when you click 'surface edges' I do not want the internal surfaces recognized as surfaces. <br>
<br></div>The reason for this is that I am then exporting this to VTK format and then a FEM program I have written. And surfaces are recognized as outer solution space surfaces and have boundary conditions placed on them, which of course, I do not want applied internally. <br>
<br></div>Is there a way to scrub the internal `surface' but retain the mesh structure?<br> <br></div><div class="gmail_extra"><br clear="all"><div>-- Trevor<br></div>
<br><br><div class="gmail_quote">On 29 July 2014 14:34, Trevor Irons <span dir="ltr"><<a href="mailto:trevorirons@gmail.com" target="_blank">trevorirons@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<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>
</blockquote></div><br></div>