Hayden Liu Weng did his honours’ project at the Chair for Computation in Engineering (TUM) and
was supervised by László Kudela M.Sc. and Prof. Dr. rer. nat. Ernst Rank.
Efficient and accurate methods for numerical integration are an essential element in many computational algorithms. This is especially important for the solution of partial differential equations on domains with internal interfaces where discontinuous integrands may have to be considered, such as in the Finite Cell Method, an immersed boundary method relying on high-order finite elements. One technique which successfully deals with the integration requirements is the smart octree introduced by Kudela in . This space decomposition method creates boundary-conforming integration sub-cells for a composed numerical quadrature. The technique, however, is limited to boundary representation (B-Rep) geometries such as the one that is available from a CAD model.
Given that this is not always the case, the aim of this project was the implementation of a version of the method compatible with implicit geometries – such as those originating from Constructive Solid Geometry (CSG). While the general scheme translates smoothly from one representation to the other, it is not trivial to deal with (sharp) feature recognition or to directly recover the non-linear character of the surface. Since the boundary is no longer explicitly stored as part of the representation, recovery techniques to identify vertices and edges become necessary. The present work uses the ideas of Kobbelt et al. from . To verify the resulting setup, a set of benchmark examples representing these difficulties were computed using both the regular- and the smart octree partitioner. Figures 1 and 2 show the integration meshes and relative error in volume respectively for a non-boundary conforming cuboid. In Figures 3 and 4, the process is repeated with a torus. Finally, Figure 5 presensts a full FCM example of a linear thermal problem making use of the new setup.