[Gmsh] BoundaryLayer Field broken?

Kilian Lackhove lackhove at ekt.tu-darmstadt.de
Thu Jan 16 11:26:01 CET 2014


Hello List,

I am trying to use the BoundaryLayer Field to create a prism boundary layer on 
one face of a simple cube and came across several issues. 

First, creating boundary layery by faces doesnt seem to work correctly, 
instead the faces edges need to be selected. Second, the mesh seems only be 
refined when hwall_n and hwall_t are set to the same values and third, quad 
recombination (the most important function) does not seem to have any effect. 
I attached a .geo file and three .png files to demonstrate the problems.

Am i doing something wrong or is there a problem in gmsh?

Kind Regards,

Kilian Lackhove



-- 
Kilian Lackhove, M.Sc.
Fachgebiet für Energie- und Kraftwerkstechnik
L1|08 114
Technische Universität Darmstadt
Jovanka-Bontschits-Straße 2
D-64287 Darmstadt
Germany

Tel.: +49 6151 16 - 5186
Fax: +49 6151 16 - 6555
e-mail: lackhove at ekt.tu-darmstadt.de
-------------- next part --------------
Mesh.LcIntegrationPrecision=1e-3;
Mesh.CharacteristicLengthMax = 0.1;
Mesh.CharacteristicLengthMin = 0.01;
Mesh.Algorithm3D = 1; // Delaunay
Mesh.Algorithm = 5; // Delaunay

lc = 0.1;

Point(1) = {0, 0, 0, lc};
Point(2) = {1, 0, 0, lc};
Point(3) = {1, 1, 0, lc};
Point(4) = {1, 1, 1, lc};
Point(5) = {0, 1, 1, lc};
Point(6) = {0, 0, 1, lc};
Point(7) = {0, 1, 0, lc};
Point(8) = {1, 0, 1, lc};

Line(1) = {7, 3};
Line(2) = {3, 4};
Line(3) = {4, 5};
Line(4) = {5, 7};
Line(5) = {7, 1};
Line(6) = {1, 2};
Line(7) = {2, 3};
Line(8) = {4, 8};
Line(9) = {8, 6};
Line(10) = {6, 5};
Line(11) = {1, 6};
Line(12) = {2, 8};

Line Loop(13) = {5, 11, 10, 4};
Plane Surface(14) = {13};
Line Loop(15) = {6, 7, -1, 5};
Plane Surface(16) = {15};
Line Loop(17) = {7, 2, 8, -12};
Plane Surface(18) = {17};
Line Loop(19) = {8, 9, 10, -3};
Plane Surface(20) = {19};
Line Loop(21) = {3, 4, 1, 2};
Plane Surface(22) = {21};
Line Loop(23) = {11, -9, -12, -6};
Plane Surface(24) = {23};

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

// specifying surfaces doesnt seem to work (See bl_test_1.png).
Field[1] = BoundaryLayer;
Field[1].FacesList = {14};
Field[1].hwall_n = 0.01;
Field[1].hwall_t = 0.01;
Field[1].hfar = 0.1;
Field[1].thickness = 0.1;
Field[1].ratio = 1.5;

// edges do work, but choosing different values for hwall_t and hwall_n doesnt seem to be recognized (See bl_test_2.png).
Field[2] = BoundaryLayer;
Field[2].EdgesList = {5, 11, 10, 4};
Field[2].hwall_n = 0.01;
Field[2].hwall_t = 0.05;
Field[2].hfar = 0.1;
Field[2].thickness = 0.1;
Field[2].ratio = 1.5;

// Quad recombination seems to be broken, too (See bl_test_3.png).
Field[3] = BoundaryLayer;
Field[3].EdgesList = {5, 11, 10, 4};
Field[3].hwall_n = 0.01;
Field[3].hwall_t = 0.01;
Field[3].hfar = 0.1;
Field[3].thickness = 0.1;
Field[3].ratio = 1.5;
Field[3].Quads = 1;

// Set to 1, 2 or 3
Background Field = 3;
-------------- next part --------------
A non-text attachment was scrubbed...
Name: bl_test_1.png
Type: image/png
Size: 268240 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20140116/fce4cd3a/attachment.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: bl_test_2.png
Type: image/png
Size: 256168 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20140116/fce4cd3a/attachment-0001.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: bl_test_3.png
Type: image/png
Size: 229862 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20140116/fce4cd3a/attachment-0002.png>