Hi Netgen team
Using ZRefine or a 3D boundary layer and then generating the second-order mesh, the triangles are expanded as a 12-node prism.
When saving the mesh and trying to open it in the Netgen GUI, the mesh doesn’t render properly.
For example:
from netgen import meshing, occ
cyl = occ.Cylinder((0,0,0), occ.Z, r=1, h=1, bottom="bot", top="top")
cyl.faces.Min(occ.Z).Identify(cyl.faces.Max(occ.Z), "bot-top", type=meshing.IdentificationType.CLOSESURFACES)
mesh = occ.OCCGeometry(cyl).GenerateMesh(maxh=0.3, quad_dominated=True)
mesh.ZRefine("bot-top", [0.25,0.5,0.75])
mesh.SecondOrder()
mesh.Save("extrude.vol")
That said, would it be possible to add an option to generate a 15-node prism instead of the 12-node one? In FreeCAD, we don’t have the 12-node prism, so if it were possible to use the 15-node one, we could add ZRefine and 3D boundary layers to our Netgen command, not only for first-order elements but also for second-order elements.
