Smart octrees: Accurately integrating discontinuous functions in 3D

Abstract This paper presents an efficient and accurate method for the integration of discontinuous functions on a background mesh in three dimensions. This task is important in computational mechanics applications where internal interfaces are present in the computational domain. The proposed method creates boundary-conforming integration subcells for composed numerical quadrature, even in the presence of sharp geometric features (e.g. edges or vertices). Similar to the octree procedure, the algorithm subdivides the cut elements into eight octants. However, the octant nodes are moved onto the interface, which allows for a robust resolution of the intersection topology in the element while maintaining algorithmic simplicity. Numerical examples demonstrate that the method is able to deliver highly accurate domain integrals with a minimal number of quadrature points. Further examples show that the proposed method provides a viable alternative to standard octree-based approaches when combined with the Finite Cell Method.

[1]  D. W. Arthur,et al.  Methods of numerical Integration (2nd edition), by Philip J. Davis and Philip Rabinowitz. Pp 612. £36·50. 1984. ISBN 0-12-206360-0 (Academic Press) , 1986, The Mathematical Gazette.

[2]  Wolfgang A. Wall,et al.  Quadrature schemes for arbitrary convex/concave volumes and integration of weak form in enriched partition of unity methods , 2013 .

[3]  Ernst Rank,et al.  Efficient and accurate numerical quadrature for immersed boundary methods , 2015, Advanced Modeling and Simulation in Engineering Sciences.

[4]  Cv Clemens Verhoosel,et al.  Image-based goal-oriented adaptive isogeometric analysis with application to the micro-mechanical modeling of trabecular bone , 2015 .

[5]  Ernst Rank,et al.  The finite cell method for three-dimensional problems of solid mechanics , 2008 .

[6]  R. Cárdenas.  NURBS-Enhanced Finite Element Method (NEFEM) , 2009 .

[7]  I. Babuska,et al.  The partition of unity finite element method: Basic theory and applications , 1996 .

[8]  Giulio Ventura,et al.  Equivalent polynomials for quadrature in Heaviside function enriched elements , 2015 .

[9]  William E. Lorensen,et al.  Marching cubes: A high resolution 3D surface construction algorithm , 1987, SIGGRAPH.

[10]  T. Fries,et al.  Higher‐order accurate integration of implicit geometries , 2016 .

[11]  Barna A. Szabó,et al.  Quasi-regional mapping for the p-version of the finite element method , 1997 .

[12]  Ulrich Gabbert,et al.  Efficient integration method for fictitious domain approaches , 2015 .

[13]  Alexander Düster,et al.  Finite and spectral cell method for wave propagation in heterogeneous materials , 2014, Computational Mechanics.

[14]  Gernot Beer,et al.  A simple approach to the numerical simulation with trimmed CAD surfaces , 2015, ArXiv.

[15]  Z. Pammer,et al.  The p–version of the finite–element method , 2014 .

[16]  Thomas-Peter Fries,et al.  Higher‐order XFEM for curved strong and weak discontinuities , 2009 .

[17]  C. Peskin.  The immersed boundary method , 2002, Acta Numerica.

[18]  Michael E. Mortenson,et al.  Geometric Modeling , 2008, Encyclopedia of GIS.

[19]  Ming-Chen Hsu,et al.  The tetrahedral finite cell method: Higher‐order immersogeometric analysis on adaptive non‐boundary‐fitted meshes , 2016 .

[20]  Juan José Ródenas,et al.  Efficient finite element methodology based on cartesian grids: application to structural shape optimization , 2013 .

[21]  Gianluca Iaccarino,et al.  IMMERSED BOUNDARY METHODS , 2005 .

[22]  Ernst Rank,et al.  Geometric modeling, isogeometric analysis and the finite cell method , 2012 .

[23]  Hong Yi,et al.  A survey of the marching cubes algorithm , 2006, Comput. Graph..

[24]  Ivo Babuška,et al.  Approximate optimal points for polynomial interpolation of real functions in an interval and in a triangle , 1995 .

[25]  Ernst Rank,et al.  Finite Cell Method: High-Order Structural Dynamics for Complex Geometries , 2015 .

[26]  W. J. Gordon,et al.  Transfinite element methods: Blending-function interpolation over arbitrary curved element domains , 1973 .

[27]  L. Kudela,et al.  Highly Accurate Subcell Integration in the Context of The Finite Cell Method , 2013 .

[28]  Dominik Schillinger,et al.  The Finite Cell Method: A Review in the Context of Higher-Order Structural Analysis of CAD and Image-Based Geometric Models , 2015 .

[29]  Grégory Legrain,et al.  A NURBS enhanced extended finite element approach for unfitted CAD analysis , 2013 .

[30]  Y. Bazilevs,et al.  Small and large deformation analysis with the p- and B-spline versions of the Finite Cell Method , 2012 .

[31]  Ernst Rank,et al.  The finite cell method for bone simulations: verification and validation , 2012, Biomechanics and modeling in mechanobiology.

[32]  T. Belytschko,et al.  A review of extended/generalized finite element methods for material modeling , 2009 .

[33]  T. Belytschko,et al.  Element‐free Galerkin methods , 1994 .

[34]  Juan José Ródenas,et al.  Exact 3D boundary representation in finite element analysis based on Cartesian grids independent of the geometry , 2015 .

[35]  N. Sukumar,et al.  Numerical integration of polynomials and discontinuous functions on irregular convex polygons and polyhedrons , 2011 .

[36]  Hyun-Jung Kim,et al.  Isogeometric analysis for trimmed CAD surfaces , 2009 .

[37]  Ernst Rank,et al.  Finite cell method , 2007 .

[38]  G. Ventura.  On the elimination of quadrature subcells for discontinuous functions in the eXtended Finite‐Element Method , 2006 .

[39]  Ernst Rank,et al.  The Finite Cell Method for linear thermoelasticity , 2012, Comput. Math. Appl..

[40]  Ming-Chen Hsu,et al.  The tetrahedral finite cell method for fluids: Immersogeometric analysis of turbulent flow around complex geometries , 2016 .