<div dir="ltr">Dear all,<div>I hope you are well.</div><div>I am trying to generate a pure hexahedral mesh to be used with deal.II (where tets are not supported).<br></div><div><br></div><div>However, I regularly encounter negative jacobians in a very limited number of elements (<5), which prevents me from simulating my geometry.</div><div><br></div><div>I have tried to force the jacobian from being positive by adding the following parameters:</div><div>Mesh.ElementOrder = 2;<br>Mesh.HighOrderOptimize = 1;</div><div>This works very well for tet mesh and I can generate my mesh in a few seconds. However, as soon as I add the decomposition into hexes:</div><div>Mesh.SubdivisionAlgorithm = 2</div><div>GMSH slows down to a crawl and the mesh is unable to be generated even after 24h of meshing</div><div><br></div><div>I was wondering if there was a solution to my problem, since deal.II only support hex meshes. Is it possible to first fully generate the tet mesh and then subdivide it into hexahedron within GMSH? I would like to keep the mapping in order to ensure that the added points are correctly mapped to the geometry. I know there are tools outside (such as tethex) that can achieve this, but I would prefer not to go down this route.</div><div><br></div><div>Thanks!</div><div>Bruno</div><div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr">
<div><span><div><b>Bruno Blais<br></b></div><div><span>Ing., Ph.D.</span><br></div><div>Professeur Adjoint / Assistant Professor<br>Génie chimique / Chemical engineering<br>Bureau/Room: A-684.2.3<br>T: +1-514-340-4711, 4533<br></div><div><img src="https://docs.google.com/uc?export=download&id=11gpzGafLd7SxTmhwQSoeSvD8cHVlhck4&revid=0B7SC3o7bKzYqc1M4c3Vma0xhYW4rSTMzWmhTSDZVVTNRV2xVPQ" width="200" height="81">                                                     <br></div><div>C.P. 6079, succ. Centre-ville<br>Montréal, QC<br>CANADA H3C 3A7<div><div dir="ltr"><a href="mailto:bruno.blais@polymtl.ca" style="color:rgb(17,85,204)" target="_blank">E-mail: bruno.blais@polymtl.ca</a></div></div></div></span></div>
</div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div>