<div dir="ltr"><div><div><div><div><div><div><div><div><div>Dear all,<br><br></div>I am using Gmsh 3.0.5 to generate a complex geometry: imagine a big tree, formed by many branches of cylinders. This tree is a hole in a cubic volume. I generate all the cylinders using <i>Cylinder</i> from the OpenCASCADE kernel. I use <i>BooleanUnion</i> to add all cylinders and then <i>BooleanDifference</i> to subtract this union of cylinders from the cubic volume. Sometimes, a couple of cylinders can overlap and I thought that was not going to be an issue, since it has worked for simpler examples. <br></div>There can be a lot of cylinders (up to 400) and their radii is originally decreasing as we go up in the tree but at first I am keeping all radii equal. <br><br></div>I have two problems when generating the mesh (<i>gmsh file.geo -3</i>):<br>1. I have managed to successfully generate the mesh for only 45 cylinders with constant radii. For 46 cylinders (and higher), Gmsh crashes. <br><br></div>2. For the original decreasing radii, for 45 cylinders, I get many errors. At first, I get<br><i>Info: Cannot bind existing OpenCASCADE solid 168 to second tag 2<br></i></div><i>Info: Could not preserve tag of 3D object 2 <->168></i><br></div>Then I get a series of errors like <i>Impossible to recover the edge..., Impossible to mesh periodic face</i>.... Then:<br></div><i>Hint: You may use -p/# to decrease the dihedral angle tolerance 0.1 <degree>.<br></i></div><i>Error: Self intersecting surface mesh, computing intersections <this could take a while><br></i></div><div>The resulting mesh has some portions of the cylinder surfaces missing. <br><br></div><div>I attached the following files: <br></div><div>- Tree with constant radii for 45 cylinders (<i>45_constant_radii.geo</i>)<br></div><div>- Tree with constant radii for 50 cylinders (<i>50_constant_radii.geo</i>)</div><div>- Tree with decreasing radii for 45 cylinders (<i>45_decreasing_radii.geo</i>)<br></div><div>I have commented the Physical Surface(0) of the cubic volume (line 10) and the Physical Volume(0) of the resulting volume (last line), in order to see the mesh on the tree surface. <br></div><div><i></i></div><i></i><div><i><br></i></div><div>Do you have any idea on how to avoid these errors? I would like to have all the cylinders, with the decreasing radii. <br><br></div><div>Thank you very much, any help is more than welcome!<br><br></div><div>Best regards,<br><br></div><div>Isabela<br></div><div><i></i></div><div><div><div><div><div><div><div><br></div></div></div></div></div></div></div></div>