DeepAI AI Chat
Log In Sign Up

Exact CutFEM Polynomial Integration

by   Eugenio Aulisa, et al.

The implementation of discontinuous functions occurs in many of today's state of the art partial differential equation solvers. In finite element methods this poses an inherent difficulty: there are no quadrature rules readily available, when integrating functions whose discontinuity falls in the interior of the element. Many approaches to this issue have been developed in recent years, among them is the equivalent polynomial technique. This method replaces the discontinuous function with a polynomial, potentially allowing for the 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 own set of problems. In particular, either adaptivity is required to capture the discontinuity or error is introduced when regularization of the discontinuous function is implemented. In the current work we eliminate both of these issues. The results of this work provide exact algebraic expressions for subdomain and interface polynomial integration, where the interface represents the boundary of the cut domain. We also provide algorithms for the implementation of these expressions for standard finite element shapes in one, two, and three dimensions, along with a hypercube of arbitrary dimension .


page 1

page 2

page 3

page 4


Conservative Discontinuous Cut Finite Element Methods

We develop a conservative cut finite element method for an elliptic coup...

Implementation of hp-adaptive discontinuous finite element methods in Dune-Fem

In this paper we describe generic algorithms and data structures for the...

Higher-Order SGFEM for One-Dimensional Interface Elliptic Problems with Discontinuous Solutions

We study a class of enriched unfitted finite element or generalized fini...

High order discontinuous cut finite element methods for linear hyperbolic conservation laws with an interface

We develop a family of cut finite element methods of different orders ba...

Efficient quadrature rules for finite element discretizations of nonlocal equations

In this paper we design efficient quadrature rules for finite element di...

Convex Optimization-Based Structure-Preserving Filter For Multidimensional Finite Element Simulations

In simulation sciences, it is desirable to capture the real-world proble...