<div dir="ltr"><div><div><div><div><div>Dear GMSH Users,<br><br></div><div>I answer my own question and have an additional one.<br></div>I looked into the debug file (tato.pos, toto.pos) to search for the reason behind my problem and what I found was that it's not that there are some additional cells but rather that there are cells missing in the other parts of the boundary layer. The tato.pos file is correct while the toto.pos file produced has some cells removed.<br>

</div>Looking deeper into the source code, I realised that it is due to the call to the function filterOverlappingElements at line 711 of the meshGFace.cpp.<br></div>Can someone tell me what is the purpose of this function ?<br>

</div>Furthermore, I found out that the only other places where it is called is at line 913 of the meshGRegion.cpp. But there that line is commented.<br></div>Is it an oversight and should the call to filterOverlappinElements be also commented in the meshGFace.cpp file ?<br>

<br>Thank you for any insight you may provide.<br><div><div><div><div><div><div><div><div class="gmail_extra"><br clear="all"><div><div dir="ltr">Nguyen Anh Khoa DOAN<div>Last year dual degree student - Université Libre de Bruxelles - SUPAERO</div>

<div>Aerospace Engineering - Major: Propulsion & Complex Systems and Simulation</div></div></div>
<br><br><div class="gmail_quote">On Wed, May 21, 2014 at 10:11 AM, Nguyen Anh Khoa Doan <span dir="ltr"><<a href="mailto:nguyen.doan@supaero.isae.fr" target="_blank">nguyen.doan@supaero.isae.fr</a>></span> wrote:<br>

<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div><div><div><div>Hello GMSH users,<br><br></div>I have a problem regarding the use of the BoundaryLayer feature in gmsh for a 2D mesh.<br>

I have tried to create a boundary layer mesh around a square and use the FanNodesLists option to get a fan around the  corners. However this creates additional quad-cells for some parts of the fan and thus, I get triangular cells with a really bad aspect ratio just beyond the boundary layer (see picture attached). The same problem arises when trying to mesh the Boundary Layer around an airfoil at the TE (see picture).<br>


<br></div>The code for the square:<br>// Gmsh project created on Fri May 16 14:29:34 2014<br>Point(1) = {-4, -4, 0, 1.0};<br>Point(2) = {-4, 4, 0, 1.0};<br>Point(3) = {4, 4, 0, 1.0};<br>Point(4) = {4, -4, 0, 1.0};<br><br>


Point(5) = {-1, -1, 0, 0.1};<br>Point(6) = {1, -1, 0, 0.1};<br>Point(7) = {1, 1, 0, 0.1};<br>Point(8) = {-1, 1, 0, 0.1};<br><br>Line(1) = {1, 2};<br>Line(2) = {2, 3};<br>Line(3) = {3, 4};<br>Line(4) = {4, 1};<br>Line(5) = {5, 6};<br>


Line(6) = {6, 7};<br>Line(7) = {7, 8};<br>Line(8) = {8, 5};<br>Line Loop(9) = {1, 2, 3, 4};<br>Line Loop(10) = {6, 7, 8, 5};<br><br>Plane Surface(11) = {9, 10};<br><br>Field[1] = BoundaryLayer;<br>Field[1].EdgesList = {5,6,7,8};<br>


//Field[1].NodesList = {5,6,7,8};<br>Field[1].FanNodesList = {5,6,7,8};<br>Field[1].hfar = 0.05;<br>Field[1].hwall_n = 0.001;<br>Field[1].hwall_t = 0.003;<br>Field[1].ratio = 1.1;<br>Field[1].thickness = 0.01;<br>Field[1].Quads = 1;<br>


BoundaryLayer Field = 1;<br><br></div>Any help would be extremely welcome.<br><br></div>Regards,<br><div><div><div><div><div><br clear="all"><div><div dir="ltr">Nguyen Anh Khoa DOAN<div>Last year dual degree student - Université Libre de Bruxelles - SUPAERO</div>


<div>Aerospace Engineering - Major: Propulsion & Complex Systems and Simulation</div></div></div>
</div></div></div></div></div></div>
</blockquote></div><br></div></div></div></div></div></div></div></div></div>