[Gmsh] boundary layer in a tube, code attached

[Hongtao Yang]

Hi,

I try to create the mesh in a tube with Gmsh, but cannot refine the
boundary layer.
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!!

Regards,
Hongtao



// Gmsh project created on Sun Nov 13 20:00:53 2016
s = 0.5e-5;
Len = 40e-6;
rp = 3e-6;

// Tube
Point(8) = {0,0,0,s};

Point(9) = {1e-5,0,0,s};
Point(10) = {0,1e-5,0,s};
Point(11) = {-1e-5,0,0,s};
Point(12) = {0,-1e-5,0,s};

Point(13) = {0,0,Len,s};

Point(14) = {1e-5,0,Len,s};
Point(15) = {0,1e-5,Len,s};
Point(16) = {-1e-5,0,Len,s};
Point(17) = {0,-1e-5,Len,s};

Circle(31) = {9,8,10};
Circle(32) = {10,8,11};
Circle(33) = {11,8,12};
Circle(34) = {12,8,9};

Circle(35) = {14,13,15};
Circle(36) = {15,13,16};
Circle(37) = {16,13,17};
Circle(38) = {17,13,14};

Line(39) = {9, 14};
Line(40) = {10, 15};
Line(41) = {11, 16};
Line(42) = {12, 17};
Line(300) = {13,8};

Line Loop(43) = {36, 37, 38, 35};
Plane Surface(44) = {43};
Physical Surface("inlet") = {44};

Line Loop(45) = {35, -40, -31, 39};
Ruled Surface(46) = {45};

Line Loop(47) = {38, -39, -34, 42};
Ruled Surface(48) = {47};

Line Loop(49) = {37, -42, -33, 41};
Ruled Surface(50) = {49};

Line Loop(51) = {36, -41, -32, 40};
Ruled Surface(52) = {51};
Physical Surface("boundary") = {46, 48, 50, 52};

Line Loop(53) = {34, 31, 32, 33};
Plane Surface(54) = {53};
Physical Surface("outlet") = {54};

Field[1] = Attractor;
Field[1].NNodesByEdge = 100;
Field[1].EdgesList = {300};
Field[2] = Threshold;
Field[2].IField = 1;
Field[2].LcMin = s;
Field[2].LcMax = s/10;
Field[2].DistMin = 0.6e-5;
Field[2].DistMax = 0.8e-5;

Background Field = 2;

Field[4] = BoundaryLayer;
Field[4].hfar = .5;
Field[4].hwall_n = .5;
Field[4].thickness = 0.2e-5;
Field[4].ratio = 1.1;
Field[4].FacesList = {46,48,50,52};
//BoundaryLayer Field = 4;


Surface Loop(55) = {44, 52, 50, 48, 46, 54};
Volume(56) = {55};
Physical Volume("Pore Throat") = {56};


------------------

Hongtao Yang

Graduate Research Assistant

Dept. of Petroleum & Geosystems Engineering

The University of Texas at Austin
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://onelab.info/pipermail/gmsh/attachments/20161121/cb9a93f0/attachment.html>