Exact CutFEM Polynomial Integration

10/01/2021
by   Eugenio Aulisa, et al.
0

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 .

READ FULL TEXT

page 1

page 2

page 3

page 4

research
05/05/2021

Conservative Discontinuous Cut Finite Element Methods

We develop a conservative cut finite element method for an elliptic coup...
research
04/25/2016

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

In this paper we describe generic algorithms and data structures for the...
research
04/15/2022

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

We study a class of enriched unfitted finite element or generalized fini...
research
01/18/2022

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...
research
08/05/2017

FEMPAR: An object-oriented parallel finite element framework

FEMPAR is an open source object oriented Fortran200X scientific software...
research
08/08/2022

Construction of Discontinuous Enrichment Functions for Enriched FEM's for Interface Elliptic Problems in 1D

We introduce an enriched unfitted finite element method to solve 1D elli...
research
03/12/2022

Arbitrary high-order methods for one-sided direct event location in discontinuous differential problems with nonlinear event function

In this paper we are concerned with numerical methods for the one-sided ...

Please sign up or login with your details

Forgot password? Click here to reset