CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   Open Source Meshers: Gmsh, Netgen, CGNS, ... (http://www.cfd-online.com/Forums/openfoam-meshing-open/)
-   -   Boundary Layer Mesh in GMSH (http://www.cfd-online.com/Forums/openfoam-meshing-open/126563-boundary-layer-mesh-gmsh.html)

Medu November 20, 2013 14:06

Boundary Layer Mesh in GMSH
 
1 Attachment(s)
Hi,

I am trying to create boundary layers around some walls in my model but there is something screwing up.

1. When I run the code only with BoundaryLayer 1 it works. But if I run it with both BoundaryLayers it only creates BoundaryLayer 2.

2. BoundaryLayer 2 isn't working correct (see picture).


Does anyone know how the soultion of my problems could be.

Code:

capillary = 0.0001;
Jet_Disrupter = 0.0001;
Walls = 0.001;

Point(1) = {0, 0, 0, capillary};
Point(2) = {0.006, 0, 0, capillary};
Point(3) = {0.006, 0.000254, 0, capillary};
Point(4) = {0.006, 0.0008, 0, Walls};
Point(5) = {0.004, 0.0008, 0, Walls};
Point(6) = {0.004, 0.0064, 0, Walls};
Point(7) = {0, 0.0064, 0, Walls};
Point(8) = {0, 0.0123, 0, Walls};
Point(9) = {0.0045, 0.0123, 0, Walls};
Point(10) = {0.0595, 0.0123, 0, Walls};
Point(11) = {0.1043, 0.0011, 0, Walls};
Point(12) = {0.1045, 0.0011, 0, Walls};
Point(13) = {0.1045, 0.0123, 0, Walls};
Point(14) = {0.1300, 0.0123, 0, Walls};
Point(15) = {0.1300, 0, 0, Walls};
Point(16) = {0.02415, 0.0033, 0, Jet_Disrupter};
Point(17) = {0.02465, 0.0033, 0, Jet_Disrupter};
Point(18) = {0.02415, 0, 0, Jet_Disrupter};
Point(19) = {0.02465, 0, 0, Jet_Disrupter};
Point(20) = {-0.048, 0.000254, 0, capillary};
Point(21) = {-0.048, 0, 0, capillary};
Point(22) = {-0.046, 0.0015875, 0, Walls};
Point(23) = {-0.046, 0.005, 0, Walls};
Point(24) = {-0.062, 0, 0, Walls};
Point(25) = {-0.062, 0.005, 0, Walls};
Point(26) = {-0.048, 0.0015875, 0, Walls};



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





Line Loop(27) = {13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12};
Plane Surface(28) = {27};
Rotate {{1, 0, 0}, {0, 0, 0}, -0.08} {
  Surface{28};
}
Extrude {{1, 0, 0}, {0, 0, 0}, 0.04} {
  Surface{28};
Layers{1}; // 1 Element thick
  Recombine;
}


Surface Loop(142) = {141, 59, 28, 63, 67, 71, 75, 78, 82, 86, 89, 96, 100, 104, 108, 112, 116, 120, 124, 128, 132, 136, 140};
Volume(143) = {142};
Physical Volume(144) = {1};
Physical Surface(145) = {96, 100};
Physical Surface(146) = {104, 108, 112, 116, 120, 124, 128, 132, 136};
Physical Surface(147) = {140};
Physical Surface(148) = {59, 63, 71, 67, 89, 82, 86};
Physical Surface(149) = {75, 78};
Physical Surface(150) = {141};
Physical Surface(151) = {28};



//BoundaryLayer1

Field[1] = BoundaryLayer;
Field[1].EdgesList = {13, 14, 15, 16};
Field[1].NodesList = {9, 10, 11, 12, 13};
Field[1].hfar = 0.0005;
Field[1].hwall_n = 0.0000005;
Field[1].hwall_t = 0.001;
Field[1].thickness = 0.0001; //0.0004
Field[1].ratio = 1.5;
Field[1].fan_angle = 30;
Field[1].Quads = 1;
BoundaryLayer Field = 1;

//BoundaryLayer2

Field[2] = BoundaryLayer;
Field[2].EdgesList = {3, 4, 5, 6, 7, 8, 9, 10, 11};
Field[2].NodesList = {3, 4, 5, 6, 7, 8, 20, 22, 23, 26};
Field[2].hfar = 0.0005;
Field[2].hwall_n = 0.0000005;
Field[2].hwall_t = 0.001;
Field[2].thickness = 0.0001;
Field[2].ratio = 1.5;
Field[2].fan_angle = 30;
Field[2].Quads = 1;
BoundaryLayer Field = 2;



All times are GMT -4. The time now is 20:24.