Exact subdomain and embedded interface polynomial integration in finite elements with planar cuts

The implementation of discontinuous functions occurs in many of today’s state-of-the-art partial differential equation solvers. However, in finite element methods, this poses an inherent difficulty: efficient quadrature rules available when integrating functions whose discontinuity falls in the elem...

Full description

Saved in:
Bibliographic Details
Published in:Numerical algorithms 2023-09, Vol.94 (1), p.315-350
Main Authors: Aulisa, Eugenio, Loftin, Jonathon
Format: Article
Language:English
Subjects:
Algebra
Algorithms
Automation
Complexity
Computer Science
Discontinuity
Equivalence
Finite element method
Hypercubes
Mathematical analysis
Methods
Numeric Computing
Numerical Analysis
Original Paper
Partial differential equations
Polynomials
Quadratures
Tetrahedra
Theory of Computation
Triangles
Citations: Items that this one cites
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:The implementation of discontinuous functions occurs in many of today’s state-of-the-art partial differential equation solvers. However, in finite element methods, this poses an inherent difficulty: efficient quadrature rules available when integrating functions whose discontinuity falls in the element’s interior are for low order degree polynomials, not easily extended to higher order degree polynomials, and cover a restricted set of geometries. Many approaches to this issue have been developed in recent years. Among them, one of the most elegant and versatile is the equivalent polynomial technique. This method replaces the discontinuous function with a polynomial, allowing integration to occur over the entire domain rather than integrating over complex subdomains. Although eliminating the issues involved with discontinuous function integration, the equivalent polynomial tactic introduces its problems. The exact subdomain integration requires a machinery that quickly grows in complexity when increasing the polynomial degree and the geometry dimension, restricting its applicability to lower order degree finite element families. The current work eliminates this issue. We provide algebraic expressions to exactly evaluate the subdomain integral of any degree polynomial on parent finite element shapes cut by a planar interface. These formulas also apply to the exact evaluation of the embedded interface integral. We provide recursive algorithms that avoid overflow in computer arithmetic for standard finite element geometries: triangle, square, cube, tetrahedron, and prism, along with a hypercube of arbitrary dimensions.
ISSN:1017-1398
1572-9265
DOI:10.1007/s11075-023-01502-3