CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Open Source Meshers: Gmsh, Netgen, CGNS, ...

Boundary Layer Mesh in GMSH

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   November 20, 2013, 15:06
Default Boundary Layer Mesh in GMSH
  #1
New Member
 
Join Date: Nov 2013
Posts: 2
Rep Power: 0
Medu is on a distinguished road
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;
Attached Images
File Type: jpg 1.jpg (52.8 KB, 137 views)
Medu is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
How to set boundary layer of a moving body in ICEM to a mesh zone for dynamic mesh WassimCFD FLUENT 0 July 13, 2013 14:24
An error has occurred in cfx5solve: volo87 CFX 5 June 14, 2013 17:44
Add Mesh Layers doesnt work on the whole surface Kryo OpenFOAM Native Meshers: snappyHexMesh and Others 8 September 13, 2012 09:28
Boundary layer in a pipe Clementhuon OpenFOAM Native Meshers: snappyHexMesh and Others 6 March 12, 2012 13:41
Boundary layer generation problems ivan_cozza OpenFOAM Native Meshers: snappyHexMesh and Others 0 October 6, 2010 13:47


All times are GMT -4. The time now is 03:39.