
Anisotropic geometryconforming dsimplicial meshing via isometric embeddings
Caplan, Philip Claude, Robert Haimes, David L. Caplan, Marshall C. Galbraith
Proceedings, 26th International Meshing Roundtable, Elsevier, Science Direct, September 1821 2017

INTERNATIONAL MESHING ROUNTABLE

26th International Meshing Roundtable
Barcelona, Spain
September 1821, 2017
Philip Claude Caplan, Massachusetts Institute of Technology, US, pcaplan@mit.edu
Robert Haimes, Massachusetts Institute of Technology, US, haimes@mit.edu
David L. Caplan, Massachusetts Institute of Technology, US, darmofal@mit.edu
Marshall C. Galbraith, Massachusetts Institute of Technology, US, galbramc@mit.edu
Abstract
We develop a dimensionindependent, Delaunaybased anisotropic mesh generation algorithm suitable for integration with adaptive numerical solvers. As such, the mesh produced by our algorithm conforms to an anisotropic metric prescribed by the solver as well as the domain geometry, given as a piecewise simplicial complex. Motivated by the work of L\'evy and Dassi \cite{Levy_2013_Vorpaline,Dassi_2014_Curvature_adapted_remeshing_of_CAD_surfaces,Dassi_2015_Anisotropic_Finite_Element_Mesh_Adaptation_via_Higher_Dimensional_Embedding,Dassi_2016_An_anisotropic_surface_remeshing_strategy_combining_higher_dimensional_embedding_with_radial_basis_functions}, we use a discrete manifold embedding algorithm to transform the anisotropic problem to a uniform one. This work differs from previous approaches in several ways. First, the embedding algorithm is driven by a Riemannian metric field instead of the Gauss map, lending itself to general anisotropic mesh generation problems. Second we describe our method for computing restricted Voronoi diagrams in a dimensionindependent manner which is used to compute constrained centroidal Voronoi tessellations. In particular, we compute restricted Voronoi simplices using exact arithmetic and support our data structures using convex polytope theory. Finally, since adaptive solvers require geometryconforming meshes, we offer a Steiner vertex insertion algorithm for ensuring the extracted dual Delaunay triangulation is homeomorphic to the input geometries.\newline The two major contributions of this paper are a method for isometrically embedding arbitrary meshmetric pairs in higher dimensional Euclidean spaces along with a dimensionindependent vertex insertion algorithm for producing geometryconforming Delaunay meshes. The former is demonstrated on a twodimensional anisotropic problem whereas the latter is demonstrated on both $3d$ and $4d$ problems.\newline
Download Full Paper (PDF)
Contact author(s) or publisher for availability and copyright information on above referenced article
