Abstract:We present a novel method to perform numerical integration over curved polyhedra enclosed by high-order parametric surfaces. Such a polyhedron is first decomposed into a set of triangular and/or rectangular pyramids, whose certain faces correspond to the given parametric surfaces. Each pyramid serves as an integration cell with a geometric mapping from a standard parent domain (e.g., a unit cube), where the tensor-product Gauss quadrature is adopted. As no constraint is imposed on the decomposition, certain resulting pyramids may intersect with themselves, and thus their geometric mappings may present negative Jacobian values. We call such cells the folded cells and refer to the corresponding decomposition as a folded decomposition. We show that folded cells do not cause any issues in practice as they are only used to numerically compute certain integrals of interest. The same idea can be applied to planar curved polygons as well. We demonstrate both theoretically and numerically that folded cells can retain the same accuracy as the cells with strictly positive Jacobians. On the other hand, folded cells allow for a much easier and much more flexible decomposition for general curved polyhedra, on which one can robustly compute integrals. In the end, we show that folded cells can flexibly and robustly accommodate real-world complex geometries by presenting several examples in the context of immersed isogeometric analysis, where involved sharp features can be well respected in generating integration cells.

High-Order Numerical Integration on Domains Bounded by Intersecting Level Sets

Computational-Complexity on Numerical Integrals

Geometric Integration Over Irregular Domains with topologic Guarantees

High-order numerical integration on regular embedded surfaces

High-Order Integration on regular triangulated manifolds reaches Super-Algebraic Approximation Rates through Cubical Re-parameterizations

Numerical Integration over Implicitly Defined Domains with Topological Guarantee

Numerical integration of homogeneous functions on convex and nonconvex polygons and polyhedra

Robust Numerical Integration on Curved Polyhedra Based on Folded Decompositions

Numerical integration of discontinuous functions: moment fitting and smart octree

High-order numerical integration on self-affine sets

Gauss-Type Quadrature for Highly Oscillatory Integrals with Algebraic Singularities and Applications

Numerical Integration Approach Based on Radial Integration Method for General 3D Polyhedral Finite Elements

A New Method for Numerical Evaluation of Nearly Singular Integrals over High-Order Geometry Elements in 3D BEM.

Integration Strategy for High-Order Discontinuous Galerkin(DG)Method for Solving Euler Equations

Numerical Integration Based on Bivariate Quartic Quasi-Interpolation Operators

Numerical calculation method for function integration on submanifolds of $\mathbb{R}^n$ or compact Riemannian manifolds

Interpolation by integrals on balls

An Effective Method for Numerical Evaluation of General 2D and 3D High Order Singular Boundary Integrals

Quadrature-free polytopic discontinuous Galerkin methods for transport problems

A note on the rate of convergence of integration schemes for closed surfaces

Numerical Integration Based On A Multilevel Quartic Quasi-Interpolation Operator