[Gmsh] BoundaryLayerSurface working

[Andreas Ruopp]

Dear all,

 
I’m wondering if/how the 3d prism layer functionality works in gmsh?

Attached, you will find a simple cubic channel with the boundary layer functionality. It gives no error, but unfortunately also no boundary layer L

 
Any help is appreciated!!!

 
Best regards,

 
Andy

 
 
Test.geo content:

 
cl__1 = 1;

Point(1) = {0, 0, 0, 1};

Point(2) = {1, 0, 0, 1};

Point(3) = {1, 1, 0, 1};

Point(4) = {0, 1, 0, 1};

Point(5) = {0, 0, 10, 1};

Point(6) = {1, 0, 10, 1};

Point(7) = {1, 1, 10, 1};

Point(8) = {0, 1, 10, 1};

Line(1) = {5, 6};

Line(2) = {7, 6};

Line(3) = {8, 5};

Line(4) = {8, 7};

Line(5) = {1, 2};

Line(6) = {3, 2};

Line(7) = {4, 1};

Line(8) = {4, 3};

Line(9) = {7, 3};

Line(10) = {8, 4};

Line(11) = {5, 1};

Line(12) = {6, 2};

Line Loop(14) = {1, -2, -4, 3};

Plane Surface(14) = {14};

Line Loop(16) = {5, -6, -8, 7};

Plane Surface(16) = {16};

Line Loop(18) = {2, 12, -6, -9};

Plane Surface(18) = {18};

Line Loop(20) = {7, -11, -3, 10};

Plane Surface(20) = {20};

Line Loop(22) = {1, 12, -5, -11};

Plane Surface(22) = {22};

Line Loop(24) = {4, 9, -8, -10};

Plane Surface(24) = {24};

Mesh.CharacteristicLengthFactor = 0.04;

Surface Loop(25)={14,16,18,20,22,24};

Volume(26)={25};

Field[1] = BoundaryLayer;

Field[1].FacesList = {-14,-16,-18,-20};

Field[1].hfar = 5.;

Field[1].hwall_n = 0.0000005;

Field[1].hwall_t = 0.05;

Field[1].thickness = 0.04;

Field[1].ratio = 1.5;

Field[1].fan_angle = 60;

Field[1].Quads = 1;

BoundaryLayer Field = 1;

Transfinite Line {1,2,3,4,5,6,7,8,9,10,11,12} = 100 Using Progression 1;

 
 
 
--

 
Andreas Ruopp

 
+49-711-685-60234

aruopp at ihs.uni-stuttgart.de

 
Universitaet Stuttgart

Institut fuer Stroemungsmechanik

und Hydraulische Stroemungsmaschinen (IHS)

Pfaffenwaldring 10

70550 Stuttgart

http://www.ihs.uni-stuttgart.de

 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20140102/b193410e/attachment.html>