Skip to Main content Skip to Navigation
Conference papers

Geometric Meshing of Discrete Surfaces

Patrick Laug 1 Houman Borouchaki 2 
1 Gamma3 - Automatic mesh generation and advanced methods
Inria Paris-Rocquencourt, ICD - Institut Charles Delaunay
2 Gamma3 - Automatic mesh generation and advanced methods
ICD - Institut Charles Delaunay, Inria Saclay - Ile de France
Abstract : Most computer-aided design (CAD) systems include a boundary representation (B-rep) where a surface is modeled by a collection of patches, each patch being defined by a function sigma from a parametric domain in R2 to a face in R3. To generate a mesh of such a surface, an indirect approach can be used [1]. In this case, the size and shape of the elements to be generated on the surface are specified by metrics M3, which are converted into metrics M2 in the parametric domain using the first partial derivatives of sigma in order to control meshing in two dimensions. The initial metrics M3 on the surface are either isotropic or anisotropic, and can depend from geometric features to generate a "geometric mesh". To obtain an isotropic geometric mesh, the prescribed size is proportional to the minimum radius of curvature at any point of the surface. For an anisotropic geometric mesh, the prescribed size depends on the two principal curvatures and the two corresponding principal directions. All these geometric features are based on the first and second partial derivatives of sigma. To summarize, the indirect approach for surface meshing needs the evaluation of sigma, its first partial derivatives and possibly its second partial derivatives. In practice, these computations are realized by CAD queries, which in some cases may be time-consuming, non-parallelizable, or return degenerate values. To overcome these problems, a discrete model made up of triangles can be built from the continuous CAD model. Then, an approximation of sigma and it first partial derivatives can be computed efficiently; since the first derivatives are constant on each triangle, they can even be stored to increase speed. The second partial derivatives are null but curvatures can be approximated at each vertex and then interpolated in order to generate an isotropic geometric mesh [2]. We now propose a new approach where the first and second derivatives of sigma are computed at the vertices of the discrete model, based on Taylor expansions. This approach improves efficiency to generate isotropic geometric meshes, and makes it possible to generate anisotropic geometric meshes. [1] H. Borouchaki, P. Laug, P.L. George, Parametric Surface Meshing using a Combined Advancing- Front / Generalized-Delaunay Approach, Int. J. Numer. Meth. Eng. (49)1-2:233-259, 2000. [2] P. Laug, H. Borouchaki, Discrete CAD Model for Visualization and Meshing, 25th International Meshing Roundtable, Procedia Engineering (163)149–161, 2016.
Document type :
Conference papers
Complete list of metadata
Contributor : Patrick Laug Connect in order to contact the contributor
Submitted on : Friday, September 22, 2017 - 10:14:00 AM
Last modification on : Friday, January 21, 2022 - 3:19:28 AM


  • HAL Id : hal-01591550, version 1


Patrick Laug, Houman Borouchaki. Geometric Meshing of Discrete Surfaces. 14th U.S. National Congress on Computational Mechanics (USNCCM) - Symposium on Trends in Unstructured Mesh Generation (STUMG), Jul 2017, Montreal, QC, Canada. ⟨hal-01591550⟩



Record views


Files downloads