VoroCrust is the first provably correct algorithm for conforming Voronoi meshing of non-convex and non-manifold domains with guarantees on the quality of both surface and volume elements.

Polyhedral meshes are increasingly becoming an attractive option with particular advantages over traditional meshes for certain applications. What has been missing is a robust polyhedral meshing algorithm that can handle broad classes of domains exhibiting arbitrary curved boundaries and sharp features. In addition, the power of primal-dual mesh pairs, exemplified by Voronoi-Delaunay meshes, has been recognized as an important ingredient in numerous formulations.  VoroCrust is the first provably correct algorithm for conforming Voronoi meshing of non-convex and non-manifold domains with guarantees on the quality of both surface and volume elements.  A robust refinement process estimates a suitable sizing field that enables the careful placement of Voronoi seeds across the surface circumventing the need for clipping and avoiding its many drawbacks. The algorithm has the flexibility of filling the interior by either structured or random samples, while all sharp features are preserved in the output mesh. We demonstrate the capabilities of the algorithm on a variety of models and compare against state-of-the-art polyhedral meshing methods based on clipped Voronoi cells establishing the clear advantage of VoroCrust output. 

VoroCrust has attracted the attention of Earth Scientists at Sandia and Los Alamos National Laboratories due to its unique capabilities.  SNL is incorporating VoroCrust into a framework for simulation of deep geologic repositories for nuclear waste. A LANL geoscience webpage states: "Polyhedral meshing is important for flow and transport codes which include TOUGH2, FEHM, PFLOTRAN, and MODFLOW which are based on the two-point flux discretization. While the solution to flow/transport is stable without an orthogonal mesh, it is not accurate. Geometry conforming Voronoi tessellations are difficult to produce if the geometry is complex. So, if you want orthogonal discretization and you want complex geometry, you need conforming Voronoi tessellations."

The VoroCrust software is owned by Sandia National Laboratories.

To obtain the code, please use the form on the Contact Us page.