High-Order Quadrature on Multi-Component Domains Implicitly Defined by Multivariate Polynomials

by   Robert I. Saye, et al.

A high-order quadrature algorithm is presented for computing integrals over curved surfaces and volumes whose geometry is implicitly defined by the level sets of (one or more) multivariate polynomials. The algorithm recasts the implicitly defined geometry as the graph of an implicitly defined, multi-valued height function, and applies a dimension reduction approach needing only one-dimensional quadrature. In particular, we explore the use of Gauss-Legendre and tanh-sinh methods and demonstrate that the quadrature algorithm inherits their high-order convergence rates. Under the action of h-refinement with q fixed, the quadrature schemes yield an order of accuracy of 2q, where q is the one-dimensional node count; numerical experiments demonstrate up to 22nd order. Under the action of q-refinement with the geometry fixed, the convergence is approximately exponential, i.e., doubling q approximately doubles the number of accurate digits of the computed integral. Complex geometry is automatically handled by the algorithm, including, e.g., multi-component domains, tunnels, and junctions arising from multiple polynomial level sets, as well as self-intersections, cusps, and other kinds of singularities. A variety of accompanying numerical experiments demonstrates the quadrature algorithm on two- and three-dimensional problems, including randomly generated geometry involving multiple high curvature pieces; challenging examples involving high degree singularities such as cusps; adaptation to simplex constraint cells in addition to hyper-rectangular constraint cells; and boolean operations to compute integrals on overlapping domains.


page 1

page 2

page 3

page 4


A Connected Component Labeling Algorithm for Implicitly-Defined Domains

A connected component labeling algorithm is developed for implicitly-def...

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

We present a high-order method that provides numerical integration on vo...

A coupled discontinuous Galerkin-Finite Volume framework for solving gas dynamics over embedded geometries

We present a computational framework for solving the equations of invisc...

A priori error analysis of high-order LL* (FOSLL*) finite element methods

A number of non-standard finite element methods have been proposed in re...

Convergent adaptive hybrid higher-order schemes for convex minimization

This paper proposes two convergent adaptive mesh-refining algorithms for...

Tailored meshing for parallel 3D electromagnetic modeling using high-order edge elements

We present numerical experiments for geophysics electromagnetic (EM) mod...

Please sign up or login with your details

Forgot password? Click here to reset