Harmonic triangulations

We introduce the notion of harmonic triangulations: a harmonic triangulation simultaneously minimizes the Dirichlet energy of all piecewise linear functions. By a famous result of Rippa, Delaunay triangulations are the harmonic triangulations of planar point sets. We prove by explicit counterexample that in 3D a harmonic triangulation does not exist in general. However, we show that bistellar flips are harmonic: if they decrease Dirichlet energy for one set of function values, they do so for all. This observation gives rise to the notion of locally harmonic triangulations. We demonstrate that locally harmonic triangulations can be efficiently computed, and efficiently reduce sliver tetrahedra. The notion of harmonic triangulation also gives rise to a scalar measure of the quality of a triangulation, which can be used to prioritize flips and optimize the position of vertices. Tetrahedral meshes generated by optimizing this function generally show better quality than Delaunay-based optimization techniques.

[1]  Ulrich Pinkall,et al.  Computing Discrete Minimal Surfaces and Their Conjugates , 1993, Exp. Math..

[2]  Herbert Edelsbrunner,et al.  Sliver exudation , 1999, SCG '99.

[3]  LongChen,et al.  OPTIMAL DELAUNAY TRIANGULATIONS , 2004 .

[4]  David Eppstein,et al.  Dihedral bounds for mesh generation in high dimensions , 1995, SODA '95.

[5]  Jonathan Richard Shewchuk,et al.  Delaunay refinement algorithms for triangular mesh generation , 2002, Comput. Geom..

[6]  Eitan Grinspun,et al.  Discrete laplace operators: no free lunch , 2007, Symposium on Geometry Processing.

[7]  Samuel Rippa,et al.  Minimal roughness property of the Delaunay triangulation , 1990, Comput. Aided Geom. Des..

[8]  Ligang Liu,et al.  Revisiting Optimal Delaunay Triangulation for 3D Graded Mesh Generation , 2014, SIAM J. Sci. Comput..

[9]  Ronald Fedkiw,et al.  A Crystalline, Red Green Strategy for Meshing Highly Deformable Objects with Tetrahedra , 2003, IMR.

[10]  Herbert Edelsbrunner,et al.  Incremental topological flipping works for regular triangulations , 1992, SCG '92.

[11]  J. Shewchuk What Is a Good Linear Finite Element? Interpolation, Conditioning, Anisotropy, and Quality Measures , 2002 .

[12]  R. K. Smith,et al.  Mesh Smoothing Using A Posteriori Error Estimates , 1997 .

[13]  Joab R Winkler,et al.  Numerical recipes in C: The art of scientific computing, second edition , 1993 .

[14]  A. Connes,et al.  ON THE SPECTRAL CHARACTERIZATION , 2008 .

[15]  Herbert Edelsbrunner,et al.  Sliver exudation , 2000, J. ACM.

[16]  Ludmil T. Zikatanov,et al.  A monotone finite element scheme for convection-diffusion equations , 1999, Math. Comput..

[17]  Robin Sibson,et al.  Locally Equiangular Triangulations , 1978, Comput. J..

[18]  Alexander I. Bobenko,et al.  A Discrete Laplace–Beltrami Operator for Simplicial Surfaces , 2005, Discret. Comput. Geom..

[19]  R. Duffin Distributed and Lumped Networks , 1959 .

[20]  B. Joe Three-dimensional triangulations from local transformations , 1989 .

[21]  M. Yvinec,et al.  Variational tetrahedral meshing , 2005, SIGGRAPH 2005.

[22]  D. Chand,et al.  On Convex Polyhedra , 1970 .

[23]  G. Dziuk Finite Elements for the Beltrami operator on arbitrary surfaces , 1988 .

[24]  Mark Meyer,et al.  Discrete Differential-Geometry Operators for Triangulated 2-Manifolds , 2002, VisMath.

[25]  I. Babuska,et al.  ON THE ANGLE CONDITION IN THE FINITE ELEMENT METHOD , 1976 .

[26]  Hervé Delingette,et al.  Curved optimal delaunay triangulation , 2018, ACM Trans. Graph..

[27]  Jonathan Richard Shewchuk,et al.  Isosurface stuffing: fast tetrahedral meshes with good dihedral angles , 2007, ACM Trans. Graph..

[28]  J. D. Loera,et al.  Triangulations: Structures for Algorithms and Applications , 2010 .

[29]  Pierre Alliez,et al.  Interleaving Delaunay refinement and optimization for practical isotropic tetrahedron mesh generation , 2009, ACM Trans. Graph..

[30]  Long Chen,et al.  Mesh Smoothing Schemes Based on Optimal Delaunay Triangulations , 2004, IMR.

[31]  Marc Alexa,et al.  Discrete Laplacians on general polygonal meshes , 2011, ACM Trans. Graph..

[32]  Franz Aurenhammer,et al.  Power Diagrams: Properties, Algorithms and Applications , 1987, SIAM J. Comput..

[33]  Daniele Panozzo,et al.  Tetrahedral meshing in the wild , 2018, ACM Trans. Graph..

[34]  Oleg R. Musin Properties of the Delaunay triangulation , 1997, SCG '97.

[35]  Mathieu Desbrun,et al.  Weighted Triangulations for Geometry Processing , 2014, ACM Trans. Graph..

[36]  Charles L. Lawson,et al.  Transforming triangulations , 1972, Discret. Math..

[37]  M. Krízek,et al.  On the maximum angle condition for linear tetrahedral elements , 1992 .

[38]  Ligang Liu,et al.  A spectral characterization of the Delaunay triangulation , 2010, Comput. Aided Geom. Des..

[39]  Mariette Yvinec,et al.  Variational tetrahedral meshing , 2005, ACM Trans. Graph..