<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=Windows-1252">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body dir="ltr">
<div style="font-family: Calibri,Arial,Helvetica,sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Dear Christophe,<br>
<br>
Thank you! That looks like it has done the trick!<br>
<br>
Kind Regards,<br>
</div>
<div style="font-family: Calibri,Arial,Helvetica,sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div id="signature">
<div id="divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:#000000; font-family:Calibri,Arial,Helvetica,sans-serif">
<a alt="BURN" style="text-decoration:none; outline:none; display:inline-block; color:8A8EB9">
<table class="ms-rteTable-clear" style="border-collapse:collapse; table-layout:fixed; border:1px solid transparent" cellspacing="0" height="133" width="445">
<tbody>
<tr class="ms-rteTableEvenRow-clear" style="border-collapse:collapse; border:1px solid transparent">
<td class="ms-rteTableEvenCol-clear" style="border-collapse:collapse; width:168px; border:1px solid transparent">
<a alt="BURN" style="text-decoration:none; outline:none; display:inline-block; color:8A8EB9"></a><a alt="FLOW" style="text-decoration:none; outline:none; display:inline-block; color:8A8EB9"><br>
<img title="Go with the FLOW" alt="FLOW" style="color:rgb(224,52,18)" height="49" width="141" src="https://dl.dropboxusercontent.com/s/6l3hn54rmrkkf2a/FLOW_logo_baseline.png?dl=0"></a><br>
<br>
<img title="“VUBâ€" alt="VUB" style="color:rgb(224,52,18)" height="64" width="150" src="https://dl.dropboxusercontent.com/s/3gvh7fp6ldc5otg/VUB_logo.jpg?dl=0"><br>
</td>
<td class="ms-rteTableOddCol-clear" style="border-collapse:collapse; width:270px; border:1px solid transparent">
<span style="font-size:12pt; font-weight:bold">Nathan Ricks</span><span style="font-size:12pt"> | </span><span style="font-size:12pt">PhD Student</span><br>
<span style="font-size:10pt"></span><br>
<span style="font-size:10pt">Department of Mechanical Engineering</span><br>
<span style="font-size:10pt">Thermo and Fluid dynamics (FLOW)</span><br>
<span style="font-size:8pt"><br>
<span style="font-size:10pt">Brussels Faculty of Engineering (</span><a href="http://bruface.eu/" style="color:rgb(224,52,18)"><font color="E03412"><span class="s3" style="font-size:10pt">BRUFACE</span></font></a><span style="font-size:10pt">)</span></span><br>
<span style="font-size:8pt"><span style="font-size:10pt">Vrije Universiteit Brussel (</span><a href="http://www.vub.ac.be" style="color:rgb(224,52,18)"><font color="E03412"><span class="s3" style="font-size:10pt">VUB</span></font></a><span style="font-size:10pt">)</span></span><br>
</td>
</tr>
</tbody>
</table>
<br>
</a><br>
<br>
<span style="font-size:8pt"></span></div>
</div>
<div id="appendonsend"></div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> Christophe Geuzaine <cgeuzaine@uliege.be><br>
<b>Sent:</b> Thursday, 11 April 2019 9:37 PM<br>
<b>To:</b> Nathan RICKS<br>
<b>Cc:</b> gmsh@onelab.info<br>
<b>Subject:</b> Re: [Gmsh] Errors when increasing NNodesByEdge on a small edge</font>
<div> </div>
</div>
<div class="BodyFragment"><font size="2"><span style="font-size:11pt;">
<div class="PlainText"><br>
<br>
> On 11 Apr 2019, at 18:22, Nathan RICKS <Nathan.James.Ricks@vub.be> wrote:<br>
> <br>
> Sorry for the repeated messages, just trying to sort out the issue and coming across new things/information that might help resolve the issue. In the .geo file attached, it demonstrates what "I think" is the cause of my errors - if you generate the 2D mesh
and look at the trailing edge, there are cells/nodes crossing over the lines/edges where they obviously shouldn't be. This occurs on both Mac and Linux for me at least. Have you seen this issue before by any chance?<br>
> <br>
<br>
It seems that you ran into a case where the small random perturbation we use in our 2D initial triangulation is too large. A new version of the algorithm without any such perturbation will be available in a future release; but in the meantime you can set e.g.
"gmsh -rand 1e-12" (or set Mesh.RandomFactor = 1e-12 in your script), which seems to solve the issue.<br>
<br>
Christophe<br>
<br>
<br>
> Thanks again,<br>
> <br>
> Nathan Ricks<br>
> VUB<br>
> <br>
> <br>
> <br>
> <br>
> <br>
> From: Nathan RICKS<br>
> Sent: Thursday, 11 April 2019 11:41 AM<br>
> To: Christophe Geuzaine<br>
> Cc: gmsh@onelab.info<br>
> Subject: Re: [Gmsh] Errors when increasing NNodesByEdge on a small edge<br>
> <br>
> I'm not sure if this helps determine the issue, but it seems if I generate the mesh using the Delaunay algorithm, there are cells that break through boundaries of the surfaces, and if I use Frontal-Delaunay (which I was) these correspond to the broken cells
shown in the first screenshots. I've attached two screenshots showing this, if you compare between the two you can see "correspondence of errors".<br>
> <br>
> From: Nathan RICKS<br>
> Sent: Thursday, 11 April 2019 9:00 AM<br>
> To: Christophe Geuzaine<br>
> Cc: gmsh@onelab.info<br>
> Subject: Re: [Gmsh] Errors when increasing NNodesByEdge on a small edge<br>
> <br>
> Dear Prof. Geuzaine,<br>
> <br>
> Thanks for the response. Indeed, I just tested this using Gmsh on my Mac, and the resultant mesh runs fine.<br>
> <br>
> So it seems the problem is restricted to the workstation, which is running Linux 64-bit (red hat). I just tried instead compiling from the source code and the error persists there also unfortunately. Are you aware of anything that may cause this? It's for
an optimisation routine, so thousands of meshes are generated, so using my Mac isn't feasible unfortunately.<br>
> <br>
> In regards to the boundary layer, I was originally using the BL field directly, however ran into issues as I couldn't have a growth rate "around" the airfoil spline (to have narrower cells at leading and trailing edges), as well as perpendicular to the airfoil.
Also, at the trailing edge using the fan node resulted in too small cells causing numerical errors, and without it the cells were too bad quality, also causing errors. Using the transfinite surfaces gave me flexibility in the boundary layers, allowing two
different growth rates perpendicular to the surface. With that it was running very well, until now on trying to include the blunt TE.<br>
> <br>
> Thank you again,<br>
> Kind Regards,<br>
> <br>
> Nathan Ricks.<br>
> <br>
> <br>
> <br>
> <br>
> From: Christophe Geuzaine <cgeuzaine@uliege.be><br>
> Sent: Thursday, 11 April 2019 12:09 AM<br>
> To: Nathan RICKS<br>
> Cc: gmsh@onelab.info<br>
> Subject: Re: [Gmsh] Errors when increasing NNodesByEdge on a small edge<br>
> <br>
> <br>
> <br>
> > On 10 Apr 2019, at 21:57, Nathan RICKS <Nathan.James.Ricks@vub.be> wrote:<br>
> > <br>
> > Hi all,<br>
> > <br>
> > I have quite a particular problem, but I'm hoping someone can help. I'm meshing an airfoil using transfinite surfaces in the boundary layer.<br>
> > <br>
> > I've just had to add a blunt trailing edge (mesh was working well with sharp trailing edge) to the airfoil, and now I have this one error stopping me, and I just can't figure it out...<br>
> > <br>
> > When I increase the NNodesByEdge on Field[11] from 6 to 15, the mesh ends up bad quality, with cell corners not matching other cells at the boundary of surfaces, and it causes SU2 to diverge instantly. It runs fine with 6, but crashes at 15. I've attached
two images showing that it is a mesh issue causing this (nothing else is being changed in the setup).<br>
> > <br>
> <br>
> I cannot seem to reproduce the issue on my Mac. Attached is the mesh generated on your .geo file:<br>
> <br>
> <br>
> > Other things such as reducing lcMin on Field[12] can also make it crash for example, however the above has visible errors in the mesh.<br>
> > <br>
> <br>
> Idem... What OS are you running on?<br>
> <br>
> PS: instead of generating the BL by hand, why not using the boundary layer field directly? Cf. e.g.
<br>
> <br>
> <a href="https://gitlab.onelab.info/gmsh/gmsh/blob/master/benchmarks/2d/naca12BL.geo">
https://gitlab.onelab.info/gmsh/gmsh/blob/master/benchmarks/2d/naca12BL.geo</a><br>
> <a href="https://gitlab.onelab.info/gmsh/gmsh/blob/master/benchmarks/2d/naca12_2d_trailing.geo">
https://gitlab.onelab.info/gmsh/gmsh/blob/master/benchmarks/2d/naca12_2d_trailing.geo</a><br>
> <br>
> > Does anyone have an idea on why I can't use the Distance/Threshold field as normal on this edge?<br>
> > <br>
> > If you can provide any guidance that would be greatly appreciated, it's so close...! I've already tested with the latest release with the same results.<br>
> > <br>
> > Thanks a lot,<br>
> > Kind Regards,<br>
> > <br>
> > Nathan Ricks<br>
> > <br>
> > <br>
> > <br>
> > <br>
> > <br>
> > <NNodes_15.png><NNodes_6.png><airfoil_level1.geo>_______________________________________________<br>
> > gmsh mailing list<br>
> > gmsh@onelab.info<br>
> > <a href="http://onelab.info/mailman/listinfo/gmsh">http://onelab.info/mailman/listinfo/gmsh</a><br>
> <br>
> — <br>
> Prof. Christophe Geuzaine<br>
> University of Liege, Electrical Engineering and Computer Science <br>
> <a href="http://www.montefiore.ulg.ac.be/~geuzaine">http://www.montefiore.ulg.ac.be/~geuzaine</a><br>
> <br>
> <br>
> <br>
> <airfoil_level1_edit.geo>_______________________________________________<br>
> gmsh mailing list<br>
> gmsh@onelab.info<br>
> <a href="http://onelab.info/mailman/listinfo/gmsh">http://onelab.info/mailman/listinfo/gmsh</a><br>
<br>
— <br>
Prof. Christophe Geuzaine<br>
University of Liege, Electrical Engineering and Computer Science <br>
<a href="http://www.montefiore.ulg.ac.be/~geuzaine">http://www.montefiore.ulg.ac.be/~geuzaine</a><br>
<br>
<br>
<br>
</div>
</span></font></div>
</body>
</html>