Numerical Integration on Hyperrectangles in Isoparametric Unfitted Finite Elements

摘要

We consider the recently introduced idea of isoparametric unfitted finite element methods and extend it from simplicial meshes to quadrilateral and hexahedral meshes. The concept of the isoparametric unfitted finite element method is the construction of a mapping from a reference configuration to a higher order accurate configuration where the reference configuration is much more accessible for higher order quadrature. The mapping is based on a level set description of the geometry and the reference configuration is a lowest order level set approximation. On simplices this results in a piecewise planar and continuous approximation of the interface. With a simple geometry decomposition quadrature rules can easily be applied based on a tesselation. On hyperrectangles the reference configuration corresponds to the zero level of a multilinear level set function which is not piecewise planar. In this work we explain how to achieve higher order accurate quadrature with only positive quadrature weights also in this case.