<div dir="ltr"><div>Dear all,</div><div><br></div><div>First of all, I would like to thanks the Gmsh developers and team for providing such a robust and versatile meshing tool.</div><div><br></div><div>I encountered a strange behaviour related to periodicities. When I run the attached file <b>not_working.geo</b>
with an alpha parameter (line 8 of the script) different to pi/4 (try
pi/8 if you want to see the error), the code correctly finds the
corresponding periodic surfaces <b>2-5 and 8-1</b> and performs the 2D meshing without any error. <br></div><div>Nevertheless, the 3D meshing outputs the following errors:</div><div><br></div><div><i>Info : Found two exactly self-intersecting facets.<br>Info : 1st: [297, 2694, 3786] #2<br>Info : 2nd: [297, 2694, 9] #2<br></i></div><div><i>Error : Invalid boundary mesh (overlapping facets) on surface 2 surface 2<br>Error : No elements in volume 35</i></div><div><br></div><div>Taking
a look more carefully at the 2D surface mesh, I found that 2 out of the
4 periodic faces have overlapping triangles (see attached image <b>not_working.png</b>).
This seems to be related to the orientation of the ellipse being
different on each periodic surface, but this just an hypothesis.<br></div><div>I tried with the Mesh.Algorithm options 1,5 and 6 without any success.<br></div><div><br></div><div>This
bug does not appear when using a pi/4 alpha parameter, which might be
due to the fact that the previously badly oriented edges are now
circles.<br></div><div><br></div><div>I managed to find a workaround that works for any alpha parameter (see <b>working.geo</b>) using ThruSections functions (see line 45 of <b>working.geo</b>), if that can help you debugging this issue.</div><div><br></div><div>My Gmsh building options are the following (printed by Gmsh):</div><div><br></div><div><i>Version : 4.5.3-git-5d81a7c49<br>License : GNU General Public License<br>Build OS : Linux64<br>Build date : 20200213<br>Build host : PC-Alain-BEROD<br>Build
options : 64Bit ALGLIB Ann Bamg Blas[generic] Blossom DIntegration
Dlopen DomHex Fltk Gmm Hxt Hxt3D Jpeg Kbipack Lapack[generic]
LinuxJoystick MathEx Med Mesh Metis Mmg3d Mpeg Netgen ONELAB
ONELABMetamodel OpenCASCADE OpenCASCADE-CAF OpenGL OptHom Parser Plugins
Png Post QuadTri Solver TetGen/BR Voro++ Zlib<br>FLTK version : 1.3.4<br>OCC version : 7.3.0<br>MED version : 4.0.0<br>Packaged by : alain<br>Web site : <a href="http://gmsh.info" target="_blank">http://gmsh.info</a><br>Issue tracker : <a href="https://gitlab.onelab.info/gmsh/gmsh/issues" target="_blank">https://gitlab.onelab.info/gmsh/gmsh/issues</a></i></div><div><br></div><div>Do not hesitate to contact me directly if something is not clear.</div><div><br></div><div>Kind regards,</div><div><br></div><div>Alain BEROD</div></div>