Cannot create 2D mesh for simple geometry
Hey everyone, I'm having issues creating a 2D mesh on fairly simple geometry. I think the problem may have something to do with some circles I am using which are tangent to one another, but I'm not sure how to address the issue. I've isolated the geometry which is failing to mesh in the following .geo file:
Reproduced here for convenience:
cl__1 = 1e-06;
Point(1) = {0, 0, 0, cl__1};
Point(2) = {6e-05, 0, 0, cl__1};
Point(3) = {6e-05, -6e-05, 0, cl__1};
Point(4) = {0.00012, 0, 0, cl__1};
Point(5) = {0.00019, 0, 0, cl__1};
Point(6) = {0.00019, 7e-05, 0, cl__1};
Point(7) = {0.00026, 0, 0, cl__1};
Point(8) = {0, 0, 3e-05, cl__1};
Point(9) = {6e-05, 0, 3e-05, cl__1};
Point(10) = {6e-05, -6e-05, 3e-05, cl__1};
Point(11) = {0.00012, 0, 3e-05, cl__1};
Point(12) = {0.00019, 0, 3e-05, cl__1};
Point(13) = {0.00019, 7e-05, 3e-05, cl__1};
Point(14) = {0.00026, 0, 3e-05, cl__1};
Circle(1) = {1, 2, 3};
Circle(2) = {3, 2, 4};
Circle(3) = {4, 5, 6};
Circle(4) = {6, 5, 7};
Circle(5) = {8, 9, 10};
Circle(6) = {10, 9, 11};
Circle(7) = {11, 12, 13};
Circle(8) = {13, 12, 14};
Line(9) = {1, 8};
Line(10) = {14, 7};
Curve Loop(1) = {-8, -7, -6, -5, -9, 1, 2, 3, 4, -10};
Plane Surface(1) = {1};
Physical Surface("sidewall") = {1};
Here is the output of the gmsh console
-------------------------------------------------------
Version : 4.8.4
License : GNU General Public License
Build OS : Linux64-sdk
Build date : 20210428
Build host : gmsh.info
Build options : 64Bit ALGLIB ANN Bamg Blas[petsc] Blossom Cgns DIntegration Dlopen DomHex Eigen Fltk Gmm Hxt Jpeg[fltk] Kbipack Lapack[petsc] LinuxJoystick MathEx Med Mesh Metis Mmg Mpeg Netgen ONELAB ONELABMetamodel OpenCASCADE OpenCASCADE-CAF OpenGL OptHom PETSc Parser Plugins Png[fltk] Post QuadTri Solver TetGen/BR Voro++ Zlib[fltk]
FLTK version : 1.4.0
PETSc version : 3.14.4 (real arithmtic)
OCC version : 7.6.0
MED version : 4.1.0
Packaged by : geuzaine
Web site : https://gmsh.info
Issue tracker : https://gitlab.onelab.info/gmsh/gmsh/issues
-------------------------------------------------------
Info : Running '/home/nstone/py/lib/python3.8/site-packages/gmsh-4.8.4-Linux64-sdk/lib/../bin/gmsh' [Gmsh 4.8.4, 1 node, max. 1 thread]
Info : Started on Fri Aug 6 17:00:29 2021
Info : Reading '/home/nstone/Desktop/Tucker_Research/hdf/sidewall_edited.geo'...
Info : Done reading '/home/nstone/Desktop/Tucker_Research/hdf/sidewall_edited.geo'
Info : Meshing 1D...
Info : Meshing curve 1 (Circle)
Info : Meshing curve 2 (Circle)
Info : Meshing curve 3 (Circle)
Info : Meshing curve 4 (Circle)
Info : Meshing curve 5 (Circle)
Info : Meshing curve 6 (Circle)
Info : Meshing curve 7 (Circle)
Info : Meshing curve 8 (Circle)
Info : Meshing curve 9 (Line)
Info : Meshing curve 10 (Line)
Info : Done meshing 1D (Wall 0.0102664s, CPU 0.010042s)
Info : 884 nodes 894 elements
Info : Meshing 2D...
Info : Meshing surface 1 (Plane, Frontal-Delaunay)
Info : :-( There are 4 intersections in the 1D mesh (curves 1 4 5 8)
Info : 8-| Splitting those edges and trying again
Info : :-( There are 10 intersections in the 1D mesh (curves 1 1 4 4 4 5 5 5 8 8)
Info : 8-| Splitting those edges and trying again
Info : :-( There are 11 intersections in the 1D mesh (curves 1 4 4 5 5 8 8 1 4 5 8)
Info : 8-| Splitting those edges and trying again
Info : :-( There are 11 intersections in the 1D mesh (curves 1 4 4 5 8 8 1 4 5 5 8)
Info : 8-| Splitting those edges and trying again
Info : :-( There are 12 intersections in the 1D mesh (curves 1 4 4 5 8 8 1 5 1 5 4 8)
Info : 8-| Splitting those edges and trying again
Info : :-( There are 12 intersections in the 1D mesh (curves 1 5 1 5 1 4 5 8 4 4 8 8)
Info : 8-| Splitting those edges and trying again
Info : :-( There are 12 intersections in the 1D mesh (curves 1 4 5 8 1 4 4 5 5 8 8 1)
Info : 8-| Splitting those edges and trying again
Info : :-( There are 12 intersections in the 1D mesh (curves 4 5 5 8 4 8 4 8 1 1 1 5)
Info : 8-| Splitting those edges and trying again
Info : :-( There are 12 intersections in the 1D mesh (curves 1 5 4 8 4 5 8 1 4 8 1 5)
Info : 8-| Splitting those edges and trying again
Info : :-( There are 11 intersections in the 1D mesh (curves 4 4 8 4 8 1 5 1 5 1 5)
Info : 8-| Splitting those edges and trying again
Info : :-( There are 8 intersections in the 1D mesh (curves 1 1 4 5 5 8 1 5)
Info : 8-| Splitting those edges and trying again
Info : Meshing surface 1 (Plane, Frontal-Delaunay)
Info : :-( There are 7 intersections in the 1D mesh (curves 1 5 8 1 1 5 4)
Info : 8-| Splitting those edges and trying again
Info : :-( There are 11 intersections in the 1D mesh (curves 1 1 4 5 5 8 8 4 5 1 4)
Info : 8-| Splitting those edges and trying again
Info : :-( There are 9 intersections in the 1D mesh (curves 1 5 8 4 4 4 8 5 8)
Info : 8-| Splitting those edges and trying again
Info : :-( There are 9 intersections in the 1D mesh (curves 8 4 5 4 8 1 5 4 8)
Info : 8-| Splitting those edges and trying again
Info : :-( There are 8 intersections in the 1D mesh (curves 1 5 4 8 4 8 8 4)
Info : 8-| Splitting those edges and trying again
Info : :-( There are 8 intersections in the 1D mesh (curves 4 8 5 5 4 4 1 8)
Info : 8-| Splitting those edges and trying again
Info : :-( There are 5 intersections in the 1D mesh (curves 5 4 8 1 5)
Info : 8-| Splitting those edges and trying again
Info : :-( There are 5 intersections in the 1D mesh (curves 5 8 1 5 4)
Info : 8-| Splitting those edges and trying again
Info : :-( There are 7 intersections in the 1D mesh (curves 1 5 4 8 8 4 4)
Info : 8-| Splitting those edges and trying again
Info : :-( There are 8 intersections in the 1D mesh (curves 1 5 8 4 4 8 8 4)
Info : 8-| Splitting those edges and trying again
Info : :-( There are 7 intersections in the 1D mesh (curves 4 4 8 5 4 8 1)
Info : 8-| Splitting those edges and trying again
Info : Meshing surface 1 (Plane, Frontal-Delaunay)
Info : :-( There are 6 intersections in the 1D mesh (curves 5 1 4 8 1 5)
Info : 8-| Splitting those edges and trying again
Info : :-( There are 8 intersections in the 1D mesh (curves 1 5 5 1 8 1 5 4)
Info : 8-| Splitting those edges and trying again
Info : :-( There are 10 intersections in the 1D mesh (curves 4 8 4 8 1 4 1 5 5 1)
Info : 8-| Splitting those edges and trying again
Info : :-( There are 6 intersections in the 1D mesh (curves 1 5 4 1 5 8)
Info : 8-| Splitting those edges and trying again
Error : Unable to recover the edge 2211 (81/172) on curve 4 (on surface 1)
Warning : Aborted 2D meshing
Info : Done meshing 2D (Wall 0.515472s, CPU 0.507859s)
Info : 1113 nodes 1123 elements
Error : ------------------------------
Error : Mesh generation error summary
Error : 1 warning
Error : 1 error
Error : Check the full log for details
Error : ------------------------------
Please let me know if anyone has any suggestions, thanks so much for reading this far!