Compact feature-aware Hermite-style high-order surface reconstruction

High-order surface reconstruction is now a commonly used technique for mesh generation and adaptation, geometric and physical modeling, and high-order numerical methods for solving partial differential equations (PDEs). However, surface reconstruction from a relatively coarse mesh remained a challenging problem, especially for surfaces with sharp features. In this paper, we introduce a new method to address this challenge by improving the previous state of the art, including continuous moving frames (CMF) and weighted averaging of local fittings (WALF) (Eng Comput 28 (2012)), in two aspects. First, we significantly improve the robustness of reconstruction from coarse meshes using a Hermite-style least squares approximation to incorporate normals of the surface and tangents of the feature curves. Second, we ensure both $$G^{0}$$ G 0 continuity and high-order accuracy of the reconstruction near sharp features by generating parametric surface elements with an iterative feature-aware parameterization. We present the theoretical framework of our method and compare it against point-based methods in terms of accuracy and stability. We demonstrate the use of the proposed technique in generating high-order meshes for finite element methods, and show that it enables nearly identical solutions as using the meshes generated from exact geometry, while allowing additional flexibility.

[1]  Holger Wendland,et al.  Piecewise polynomial, positive definite and compactly supported radial functions of minimal degree , 1995, Adv. Comput. Math..

[2]  Gene H. Golub,et al.  Matrix computations , 1983 .

[3]  Dereck S. Meek,et al.  A triangular G1 patch from boundary curves , 1996, Comput. Aided Des..

[4]  David Levin,et al.  The approximation power of moving least-squares , 1998, Math. Comput..

[5]  Xiangmin Jiao,et al.  Identification of C1 and C2 discontinuities for surface meshes in CAD , 2008, Comput. Aided Des..

[6]  P. G. Ciarlet,et al.  Interpolation theory over curved elements, with applications to finite element methods , 1972 .

[7]  Christophe Geuzaine,et al.  Gmsh: A 3‐D finite element mesh generator with built‐in pre‐ and post‐processing facilities , 2009 .

[8]  A. Huerta,et al.  Arbitrary Lagrangian–Eulerian Methods , 2004 .

[9]  J. Remacle,et al.  Gmsh: A 3‐D finite element mesh generator with built‐in pre‐ and post‐processing facilities , 2009 .

[10]  Ahmad H. Nasri,et al.  T-splines and T-NURCCs , 2003, ACM Trans. Graph..

[11]  John C. Young Higher-Order Mesh Generation Using Linear Meshes [EM Programmer's Notebook] , 2019, IEEE Antennas and Propagation Magazine.

[12]  Yi Liu,et al.  Accuracy-preserving source term quadrature for third-order edge-based discretization , 2017, J. Comput. Phys..

[13]  M. Lenoir Optimal isoparametric finite elements and error estimates for domains involving curved boundaries , 1986 .

[14]  Piecewise Polynomial , 2014, Computer Vision, A Reference Guide.

[15]  Antonio Huerta,et al.  3D NURBS‐enhanced finite element method (NEFEM) , 2008 .

[16]  Victoria Interrante,et al.  A novel cubic-order algorithm for approximating principal direction vectors , 2004, TOGS.

[17]  L. R. Scott,et al.  The Mathematical Theory of Finite Element Methods , 1994 .

[18]  Jeremy Ims,et al.  Automated low-order to high-order mesh conversion , 2018, Engineering with Computers.

[19]  Francesca Rapetti,et al.  On the generation of symmetric Lebesgue-like points in the triangle , 2012, J. Comput. Appl. Math..

[20]  D. Cohen-Or,et al.  Robust moving least-squares fitting with sharp features , 2005, ACM Trans. Graph..

[21]  Philippe G. Ciarlet,et al.  The finite element method for elliptic problems , 2002, Classics in applied mathematics.

[22]  Gerald Farin,et al.  Curves and surfaces for computer aided geometric design , 1990 .

[23]  Xiangmin Jiao,et al.  Anisotropic mesh adaptation for evolving triangulated surfaces , 2006, Engineering with Computers.

[24]  Abdul J. Jerri The Gibbs Phenomenon in Fourier Analysis, Splines and Wavelet Approximations , 1998 .

[25]  P. George,et al.  Mesh Generation: Application to Finite Elements , 2007 .

[26]  Hongyuan Zha,et al.  Consistent computation of first- and second-order differential quantities for surface meshes , 2008, SPM '08.

[27]  T. Hughes,et al.  Isogeometric analysis : CAD, finite elements, NURBS, exact geometry and mesh refinement , 2005 .

[28]  Jörg Peters,et al.  Curved PN triangles , 2001, I3D '01.

[29]  Hongyuan Zha,et al.  Simple and effective variational optimization of surface and volume triangulations , 2010, Engineering with Computers.

[30]  A. Huerta,et al.  NURBS-Enhanced Finite Element Method (NEFEM) , 2011 .

[31]  Xiangmin Jiao,et al.  Reconstructing high-order surfaces for meshing , 2011, Engineering with Computers.

[32]  A. Sluis Condition numbers and equilibration of matrices , 1969 .

[33]  Sophia Blau,et al.  Analysis Of The Finite Element Method , 2016 .

[34]  Timothy J. Tautges,et al.  AHF: array-based half-facet data structure for mixed-dimensional and non-manifold meshes , 2015, Engineering with Computers.

[35]  G. Karniadakis,et al.  Spectral/hp Element Methods for Computational Fluid Dynamics , 2005 .