<div dir="ltr"><div>Hi,</div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Jun 17, 2020 at 10:54 AM Christophe Geuzaine <<a href="mailto:cgeuzaine@uliege.be">cgeuzaine@uliege.be</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><br>
Hi Juan,<br>
<br>
This is expected: to get a 3D mesh that satisfies the Delaunay criterion we would need to modify the surface mesh - this is something we try hard to avoid.<br></blockquote><div><br></div><div><br></div><div><div>In TetGen, I started experimenting with using a background mesh from the initial mesh generated from the PLC.  Specifying smaller edge lengths on the nodes of the bad elements results in a greatly improved mesh.</div><div><br></div><div>For a simple cube experiment:</div><div>If I subdivided a tetrahedron into 4, using the the center of the circumsphere as the common vertex:</div><div>    The sum of the volumes of the subdivided elements was 50% higher than the actual volume of the cube.</div><div><br></div><div>Using the background mesh, where I specified smaller edge lengths at the vertices of the bad elements,</div><div>    The sum of the volumes of the subdivided elements matched the volume of the cube.</div><div><br></div><div>I will port this refinement strategy back to Gmsh when I have the opportunity. </div><div><br></div><div>Regards,</div><div><br></div><div>Juan</div></div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<br>
Christophe<br>
<br>
<br>
> On 14 Jun 2020, at 20:17, Juan Sanchez <<a href="mailto:juan.e.sanchez@gmail.com" target="_blank">juan.e.sanchez@gmail.com</a>> wrote:<br>
> <br>
> Hello,<br>
> <br>
> I am having a hard time getting a 3d delaunay mesh, so that the circumcenters of the tetrahedra are inside the element.<br>
> <br>
> In my simulator, I am seeing that about 50% of the elements are not delaunay if:<br>
> Mesh.CharacteristicLengthExtendFromBoundary = 1;<br>
> <br>
> and about 90% are not Delaunay if:<br>
> Mesh.CharacteristicLengthExtendFromBoundary = 0;<br>
> <br>
> Is there a way to improve this when running the tool from the command line or with an option in the geo file?<br>
> <br>
> I am using:<br>
> /Applications/Gmsh.app/Contents/MacOS/gmsh -format msh2 -3  3dblock.geo<br>
> <br>
> with version 4.5.6 on macOS 10.14.6 (Mojave).<br>
> <br>
> Regards,<br>
> <br>
> Juan<br>
> <3dblock.geo>_______________________________________________<br>
> gmsh mailing list<br>
> <a href="mailto:gmsh@onelab.info" target="_blank">gmsh@onelab.info</a><br>
> <a href="http://onelab.info/mailman/listinfo/gmsh" rel="noreferrer" target="_blank">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" rel="noreferrer" target="_blank">http://www.montefiore.ulg.ac.be/~geuzaine</a><br>
<br>
<br>
<br>
</blockquote></div></div>