In the present work, we propose a new approach, the so-called compressed adaptive integration scheme (C-AIS), for the computation of the stiffness and mass matrices in fictitious domain methods requiring the integration of discontinuous functions. The novel approach extends the conventional quadtree-decomposition-based adaptive integration scheme (AIS) by an additional step, in which established image-compression techniques are exploited to decrease the number of integration sub-cells. The benefits of the C-AIS are manifold: First, the compression of the sub-cells inevitably leads to significant savings in terms of computational time required by the numerical integration. Second, the compression procedure, which is executed directly after the quadtree-decomposition algorithm, can be easily included in existing codes. Third, if applied to polynomial integrands, the C-AIS yields exactly the same accuracy as the conventional AIS. Finally, the fourth advantage is seen in the fact that the C-AIS can readily be combined with other approaches seeking a reduction of the number of integration points such as the Boolean-FCM. The efficiency of the C-AIS approach is presented in the context of the FCM based on Cartesian meshes applied to problems of linear elastostatics and modal analysis, while it is also suitable for the quadrature in other fictitious domain approaches, e.g., CutFEM and cgFEM.

No abstract

Polygonal finite elements offer an increased freedom in terms of mesh generation at the price of more complex, often rational, shape functions. Thus, the numerical integration of rational interpolants over polygonal domains is one of the challenges that needs to be solved. If, additionally, strong discontinuities are present in the integrand, e.g., when employing fictitious domain methods, special integration procedures must be developed. Therefore, we propose to extend the conventional quadtree-decomposition-based integration approach by image compression techniques. In this context, our focus is on unfitted polygonal elements using Wachspress shape functions. In order to assess the performance of the novel integration scheme, we investigate the integration error and the compression rate being related to the reduction in integration points. To this end, the area and the stiffness matrix of a single element are computed using different formulations of the shape functions, i.e., global and local, and partitioning schemes. Finally, the performance of the proposed integration scheme is evaluated by investigating two problems of linear elasticity.

E. Enhanced numerical integration scheme based on image-compression techniques: application to fictitious domain methods. Adv Model Simul Eng Sci. 2020;7:21. Publisher's NoteSpringer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

In this contribution a novel integration scheme, extending the conventional quadtree-decomposition-based approach by image compression techniques, is investigated for unfitted polygonal meshes with a particular focus on the rational Wachspress shape functions. It is shown that a meaningful reduction of integration points can be achieved without a significant loss in accuracy. However, the full potential of the method in terms of time savings can only be leveraged when applied to higher order polynomial elements. For more information on this topic see Enhanced Numerical Integration Scheme Based on Image Compression Techniques: Application to Rational Polygonal Interpolants by Petö et al. [1].