<div dir="ltr">Hi,<div><br></div><div>I try to create the mesh in a tube with Gmsh, but cannot refine the boundary layer.</div><div>I first try to use "Attractor" and "Threshold", but seems like I can only refine the inner mesh. If I set Field[2].LcMin = s/10; and Field[2].LcMax = s; the inner mesh is refined. However if I set Field[2].LcMin = s; and Field[2].LcMax = s/10;, mesh is homogeneous and is refined everywhere. Attached is the code. I also tried to use "BoundaryLayer". But gmsh always tell me that "FacesList" is an unknown option. Not sure what happens. Thank you for your help!!</div><div><br></div><div>Regards,</div><div>Hongtao</div><div><br></div><div><br></div><div><br></div><div><div>// Gmsh project created on Sun Nov 13 20:00:53 2016</div><div>s = 0.5e-5;</div><div>Len = 40e-6;</div><div>rp = 3e-6;</div><div><br></div><div>// Tube</div><div>Point(8) = {0,0,0,s};</div><div><br></div><div>Point(9) = {1e-5,0,0,s};</div><div>Point(10) = {0,1e-5,0,s};</div><div>Point(11) = {-1e-5,0,0,s};</div><div>Point(12) = {0,-1e-5,0,s};</div><div><br></div><div>Point(13) = {0,0,Len,s};</div><div><br></div><div>Point(14) = {1e-5,0,Len,s};</div><div>Point(15) = {0,1e-5,Len,s};</div><div>Point(16) = {-1e-5,0,Len,s};</div><div>Point(17) = {0,-1e-5,Len,s};</div><div><br></div><div>Circle(31) = {9,8,10};</div><div>Circle(32) = {10,8,11};</div><div>Circle(33) = {11,8,12};</div><div>Circle(34) = {12,8,9};</div><div><br></div><div>Circle(35) = {14,13,15};</div><div>Circle(36) = {15,13,16};</div><div>Circle(37) = {16,13,17};</div><div>Circle(38) = {17,13,14};</div><div><br></div><div>Line(39) = {9, 14};</div><div>Line(40) = {10, 15};</div><div>Line(41) = {11, 16};</div><div>Line(42) = {12, 17};</div><div>Line(300) = {13,8};</div><div><br></div><div>Line Loop(43) = {36, 37, 38, 35};</div><div>Plane Surface(44) = {43};</div><div>Physical Surface("inlet") = {44};</div><div><br></div><div>Line Loop(45) = {35, -40, -31, 39};</div><div>Ruled Surface(46) = {45};</div><div><br></div><div>Line Loop(47) = {38, -39, -34, 42};</div><div>Ruled Surface(48) = {47};</div><div><br></div><div>Line Loop(49) = {37, -42, -33, 41};</div><div>Ruled Surface(50) = {49};</div><div><br></div><div>Line Loop(51) = {36, -41, -32, 40};</div><div>Ruled Surface(52) = {51};</div><div>Physical Surface("boundary") = {46, 48, 50, 52};</div><div><br></div><div>Line Loop(53) = {34, 31, 32, 33};</div><div>Plane Surface(54) = {53};</div><div>Physical Surface("outlet") = {54};</div><div><br></div><div>Field[1] = Attractor;</div><div>Field[1].NNodesByEdge = 100;</div><div>Field[1].EdgesList = {300};</div><div>Field[2] = Threshold;</div><div>Field[2].IField = 1;</div><div>Field[2].LcMin = s;</div><div>Field[2].LcMax = s/10;</div><div>Field[2].DistMin = 0.6e-5;</div><div>Field[2].DistMax = 0.8e-5;</div><div><br></div><div>Background Field = 2;</div><div><br></div><div>Field[4] = BoundaryLayer;</div><div>Field[4].hfar = .5;</div><div>Field[4].hwall_n = .5;</div><div>Field[4].thickness = 0.2e-5;</div><div>Field[4].ratio = 1.1;</div><div>Field[4].FacesList = {46,48,50,52};</div><div>//BoundaryLayer Field = 4;</div><div><br></div><div><br></div><div>Surface Loop(55) = {44, 52, 50, 48, 46, 54};</div><div>Volume(56) = {55};</div><div>Physical Volume("Pore Throat") = {56};</div></div><div><br></div><div><br clear="all"><div><div class="gmail_signature"><div dir="ltr"><div>------------------</div>
<div>
<p style="margin:0in 0in 0pt"><font size="3"><font face="Times New Roman">Hongtao Yang</font></font></p><p style="margin:0in 0in 0pt"><font size="3"><font face="Times New Roman">Graduate Research Assistant</font></font></p>
<p style="margin:0in 0in 0pt"><font size="3"><font face="Times New Roman">Dept. of Petroleum & Geosystems Engineering</font></font></p>
<p style="margin:0in 0in 0pt"><font size="3"><font face="Times New Roman">The University of Texas at Austin</font></font></p></div></div></div></div>
</div></div>