Dear All<div><br></div><div>The geo file below</div><div><br></div><div><br></div><div><br></div><div><div>====================================================</div><div>lcar = 0.05;</div><div><br></div><div>Function Mysphere;</div>
<div><br></div><div> // In the following commands we use the reserved variable name</div><div> // `newp', which automatically selects a new point number. This</div><div> // number is chosen as the highest current point number, plus</div>
<div> // one. (Note that, analogously to `newp', the variables `newc',</div><div> // `news', `newv' and `newreg' select the highest number amongst</div><div> // currently defined curves, surfaces, volumes and `any entities</div>
<div> // other than points', respectively.)</div><div><br></div><div> p1 = newp; Point(p1) = {x, y, z, lcar} ;</div><div> p2 = newp; Point(p2) = {x+r,y, z, lcar} ;</div><div> p3 = newp; Point(p3) = {x, y+r,z, lcar} ;</div>
<div> p4 = newp; Point(p4) = {x, y, z+r,lcar} ;</div><div> p5 = newp; Point(p5) = {x-r,y, z, lcar} ;</div><div> p6 = newp; Point(p6) = {x, y-r,z, lcar} ;</div><div> p7 = newp; Point(p7) = {x, y, z-r,lcar} ;</div>
<div><br></div><div> c1 = newreg; Circle(c1) = {p2,p1,p7};</div><div> c2 = newreg; Circle(c2) = {p7,p1,p5};</div><div> c3 = newreg; Circle(c3) = {p5,p1,p4};</div><div> c4 = newreg; Circle(c4) = {p4,p1,p2};</div><div> c5 = newreg; Circle(c5) = {p2,p1,p3};</div>
<div> c6 = newreg; Circle(c6) = {p3,p1,p5};</div><div> c7 = newreg; Circle(c7) = {p5,p1,p6};</div><div> c8 = newreg; Circle(c8) = {p6,p1,p2};</div><div> c9 = newreg; Circle(c9) = {p7,p1,p3};</div><div> c10 = newreg; Circle(c10) = {p3,p1,p4};</div>
<div> c11 = newreg; Circle(c11) = {p4,p1,p6};</div><div> c12 = newreg; Circle(c12) = {p6,p1,p7};</div><div><br></div><div> // We need non-plane surfaces to define the spherical holes. Here we</div><div> // use ruled surfaces, which can have 3 or 4 sides:</div>
<div><br></div><div> l1 = newreg; Line Loop(l1) = {c5,c10,c4}; Ruled Surface(newreg) = {l1};</div><div> l2 = newreg; Line Loop(l2) = {c9,-c5,c1}; Ruled Surface(newreg) = {l2};</div><div> l3 = newreg; Line Loop(l3) = {c12,-c8,-c1}; Ruled Surface(newreg) = {l3};</div>
<div> l4 = newreg; Line Loop(l4) = {c8,-c4,c11}; Ruled Surface(newreg) = {l4};</div><div> l5 = newreg; Line Loop(l5) = {-c10,c6,c3}; Ruled Surface(newreg) = {l5};</div><div> l6 = newreg; Line Loop(l6) = {-c11,-c3,c7}; Ruled Surface(newreg) = {l6};</div>
<div> l7 = newreg; Line Loop(l7) = {-c2,-c7,-c12};Ruled Surface(newreg) = {l7};</div><div> l8 = newreg; Line Loop(l8) = {-c6,-c9,c2}; Ruled Surface(newreg) = {l8};</div><div> theloops[t] = newreg;</div><div> Surface Loop(theloops[t]) = {l8+1,l5+1,l1+1,l2+1,l3+1,l7+1,l6+1,l4+1};</div>
<div> </div><div>Return</div><div><br></div><div>x = 0.0 ; y = 0.0 ; z = 0 ; r = 2.0 ; t=0; </div><div>Call Mysphere ;</div><div>x = 0 ; y = 0.0 ; z = 0 ; r = 1.0 ; t=1; </div><div>Call Mysphere ;</div><div><br></div><div>
Volume (186)={theloops[]};</div><div>Physical Volume(999)= 186;</div><div>===============================</div><div><br></div></div><div>creates a mesh with constant lc=0.05</div><div>between r=1 and r=2 </div><div>I would like to have the lc value increase with growing r.</div>
<div>However, if I use different lc values for the two calls to Mysphere</div><div>the result is a sphere of radius 2 instead of the mesh for the desired shell!</div><div>Any ideas how to fix this would be appreciated!</div>
<div><br></div><div><br></div><div>regards</div><div><br></div><div>Moritz <br>-- <br>Prof M Braun Tel.:27-12-4298006/8027<br>Physics Department Fax.: 27-12-4293643<br>University of South Africa (UNISA) <br>
<a href="mailto:moritz.braun@gmail.com" target="_blank">moritz.braun@gmail.com</a><br>P.O. Box 392 <br>0003 <br>UNISA<br> South Africa <br><a href="http://moritz-braun.blogspot.com" target="_blank">http://moritz-braun.blogspot.com</a><br>
</div>