<div dir="ltr">Hi,<div><br></div><div>I've created the following mesh. For the boundary layer I've used extrude command. As I need to define physical volume for the farfield (bounded by the box and hemisphere). Here I considered two volumes created by a seperate Surface loops. One is volume is the extruded volume and the other volume is outside extruded surface (bounded by the block). </div><div><br></div><div>However, when I check the surface face, I see some overlapping in the boundary and around it (specially in volume surface).</div><div><br></div><div>Also I would like to have have qaud dominated mesh in the volume? Any method to get it ?</div><div><br></div><div>The critical task was to define the cut plane. </div><div><br></div><div><div>d = 0.0002; //partical diameter</div><div>r = d/2; //particle radius</div><div>L = 0.002; //length of channel</div><div>l = L/2; //half length of channel</div><div>H = 0.002; //height of channel</div><div>h = H/2; //half of channel height</div><div>W = 0.002; // width of channel</div><div>w = W/2; //half width of the channel</div><div>lc = 0.0001; // characteristic length of channel mesh</div><div>lc2 = 0.00001; //characteristic lengteh of sphere</div><div><br></div><div><br></div><div>//lc = H/10; charasteric length in box</div><div>//lc_s = r/40;charasteristic length in sphere</div><div><br></div><div><br></div><div>Point(1) = {0.0,0.0,0.0,lc2};</div><div>Point(2) = {r,0.0,0.0,lc2};</div><div>Point(3) = {0,r,0.0,lc2};</div><div>Circle(1) = {2,1,3};</div><div>Point(4) = {-r,0,0.0,lc2};</div><div>Point(5) = {0,-r,0.0,lc2};</div><div>Circle(2) = {3,1,4};</div><div>Circle(3) = {4,1,5};</div><div>Circle(4) = {5,1,2};</div><div>Point(6) = {0,0,-r,lc2};</div><div>//Point(7) = {0,0,0,lc2};</div><div>Circle(5) = {3,1,6};</div><div>Circle(6) = {6,1,5};</div><div>//Circle(7) = {5,1,7};</div><div>//Circle(8) = {7,1,3};</div><div>//Circle(9) = {2,1,7};</div><div>//Circle(10) = {7,1,4};</div><div>Circle(11) = {4,1,6};</div><div>Circle(12) = {6,1,2};</div><div>//Ellipse(013) = {5, 1, 1, 3};</div><div>//Ellipse(014) = {3, 1, 1, 2};</div><div>//Ellipse(015) = {2, 1, 1, 4};</div><div>//Ellipse(016) = {4, 1, 1, 5};</div><div><br></div><div><br></div><div><br></div><div>//Line Loop(13) = {2,8,-10};</div><div>//Ruled Surface(14) = {13};</div></div><div><div>//Line Loop(15) = {10,3,7};</div><div>//Ruled Surface(16) = {15};</div><div>//Line Loop(17) = {-8,-9,1};</div><div>//Ruled Surface(18) = {17};</div><div>Line Loop(19) = {-11,-2,5};</div><div>Ruled Surface(20) = {19};</div><div>Line Loop(21) = {-5,-12,-1};</div><div>Ruled Surface(22) = {21};</div><div>Line Loop(23) = {-3,11,6};</div><div>Ruled Surface(24) = {23};</div><div>//Line Loop(25) = {-7,4,9};</div><div>//Ruled Surface(26) = {25};</div><div>Line Loop(27) = {-4,12,-6};</div><div>Ruled Surface(28) = {27};</div><div>Line Loop(29) = {1,2,3,4}; //Circle created by the particle at the symmetry plane</div><div>//Recombine Surface{20,22,24,28};</div><div><br></div><div><br></div><div><br></div><div>//Geometry option ExtrudeReturnLateralEntities: Add lateral entities in lists returned by extrusion commands? Default value: 1,Saved in: General.OptionsFileName</div><div>//tmp is a list that all entities due to extrution are saved</div><div><br></div><div>old = Geometry.ExtrudeReturnLateralEntities;</div><div>Geometry.ExtrudeReturnLateralEntities = 0;</div><div>tmp[] = Extrude {</div><div>  Surface{20,22,24,28}; Layers{10, lc}; Recombine;</div><div>};</div><div>Geometry.ExtrudeReturnLateralEntities = old;</div><div><br></div><div>Point(100) = {-l,-h,-w, lc};</div><div>Point(101) = {-l,h,-w, lc};</div><div>Point(102) = {l,-h,-w, lc};</div><div>Point(103) = {l,h,-w, lc};</div><div>Point(104) = {-l,-h,0, lc};</div><div>Point(105) = {-l,h,0, lc};</div><div>Point(106) = {l,-h,0, lc};</div><div>Point(107) = {l,h,0, lc};</div><div>Line(165) = {100, 102};</div><div>Line(166) = {102, 103};</div><div>Line(167) = {103, 101};</div><div>Line(168) = {101, 100};</div></div><div><div>Line(169) = {100, 104};</div><div>Line(170) = {104, 106};</div><div>Line(171) = {106, 102};</div><div>Line(172) = {101, 105};</div><div>Line(173) = {105, 107};</div><div>Line(174) = {107, 103};</div><div>Line(175) = {105, 104};</div><div>Line(176) = {106, 107};</div><div>Line Loop(177) = {167, 168, 165, 166};</div><div>Plane Surface(178) = {177};</div><div>Line Loop(179) = {166, -174, -176, 171};</div><div>Plane Surface(180) = {179};<br></div></div><div><div>Line Loop(181) = {171, -165, 169, 170};</div><div>Plane Surface(182) = {181};</div><div>Line Loop(183) = {176, -173, 175, 170}; //symmetry plane exterior line loop</div><div>//Plane Surface(184) = {183,29};</div><div>Line Loop(185) = {172, 175, -169, -168};</div><div>Plane Surface(186) = {185};</div><div>Line Loop(187) = {172, 173, 174, 167};</div><div>Plane Surface(188) = {187};</div><div><br></div><div><br></div><div><br></div><div>Line Loop(10000) = {49,82,65,32}; //Line loop of boundary layer outer circle in symmetry plane</div><div>Plane Surface(10084) = {10000,183}; //plane surface bounded by the outer line loop of boundary layer and the symmetry plane exterior lines</div><div><br></div><div>//Recombine Surface{178,180,182,186,188,10084};</div><div><br></div><div>Physical Surface("<inlet>") = {186};</div><div>Physical Surface("<outlet>") = {180};</div><div>Physical Surface ("<symmetry_plane>") = {10084,41,58,88,71}; //41,58,88,71 are outer surfaces projected on the symmetry plane due to extrusion - examined</div><div>//Physical Surface ("<symmetry_plane>") = {183,29};</div><div><br></div><div>Physical Surface("<top>") = {188};</div><div>Physical Surface("<bottom>") = {182};</div><div>//Physical Surface("<front>") = {184};</div><div>Physical Surface("<wall>") = {178};</div><div><br></div><div><br></div><div>Physical Surface("<particle_wall>") = {22,20,28,24}; // curved surfaces of particle wall</div><div><br></div><div>Surface Loop(1000) = {22,20,28,24}; // the sphere (only need for internal mesh)</div><div>Surface Loop(1001) = {tmp[{0:6:2}]};  // The outside of the BL - index of 0,2,4,6 represent the curved surfaces of the outer boundary layer in tmp list</div><div>//Surface Loop(1001) = {46,63,80,97};  // The outside of the BL</div><div>Surface Loop(1002) = {188, 186, 10084, 180, 178, 182}; // the exterior surfaces bounded by the boundary layer exterior surface</div><div>//Surface Loop(1002) = {188, 186, 184, 180, 178, 182};</div><div><br></div><div>//Volume(1000) = {1000}; /// inside the sphere</div><div>//Volume(1001) = {1002, 1001}; /// FarField</div><div>Volume(1002) = {1001, 1002}; /// FarField outside the boundary layer</div><div>//Volume(1002) = {1000, 1002}; /// FarField outside the boundary layer</div><div><br></div><div>//Physical Volume("<flow domain>") = {1002,1,2,3,4}; //1,2,3,4 are boundary layer volumes + farfield outside the boundary layer</div><div><br></div><div><br></div><div>Mesh.Algorithm3D = 4;</div></div><div><br></div><div>Surface mesh with surface edges </div><div><img src="cid:ii_15ed2118e0e8417a" alt="Inline image 1" width="472" height="465"><br></div><div><br></div><div>Volume faces - Non uniformity around  boundary layer (Due to face overlapping?)</div><div><img src="cid:ii_15ed2204fef381a3" alt="Inline image 3" width="472" height="342"><br></div><div><img src="cid:ii_15ed21c4185978bb" alt="Inline image 2" width="472" height="319"><br></div><div><br></div><div>Cheers,</div><div>Nilanka</div><div><br></div></div>