On Thu, Jan 16, 2020 at 1:21 AM Ramakrishnan Thirumalaisamy <
thirumalaisamyr...@gmail.com> wrote:

> Dear libmesh users,
>
> I want to create sphere with tetrahedral elements. I came across this
> function MeshTools::Generation::build_square. But this function supports
> only HEX8/27 for 3D.  Is there any other functions in libmesh that I could
> use to generate sphere with tetrahedral elements?
>

Hi Ramakrishnan,

Your email mentions both spheres and build_square, but I'm assuming what
you're really after is spheres. The implementation doesn't currently exist,
but you could do something similar to what's currently done in the
build_cube() case when Tets are selected. That is, first a mesh of
hexahedra is created and then each hex is split into IIRC 24 tets. This
isn't the most economical splitting but it is very "symmetric" and
therefore easy to make the split elements line up consistently. This
splitting also introduces new mid-face nodes, which could be "snapped" to
the sphere surface as necessary.

Another option would be to configure/build libmesh with the
--disable-strict-lgpl option and then use the libmesh Tetgen interfaces to
create a Delaunay mesh of the sphere. If you end up implementing either of
those, it would be great if you could contribute them back as either
library code or possibly an augmented miscellaneous_ex6.

Note that libmesh also supports reading several different mesh file
formats, so if you generate a sphere mesh in e.g. Gmsh, you would be able
to read it in without doing any major programming yourself...

-- 
John

_______________________________________________
Libmesh-users mailing list
Libmesh-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-users

Reply via email to