Enriched Immersed Finite Element and Isogeometric Analysis -- Algorithms and Data Structures

TL;DR

This paper presents a robust, efficient preprocessing framework for immersed finite element and isogeometric analysis, enabling easier handling of complex geometries through custom algorithms and data structures.

Contribution

It introduces a novel preprocessing framework that simplifies the implementation of immersed finite element methods with complex geometries, including algorithms for enrichment and stabilization.

Findings

Framework is robust in geometric edge cases

Preprocessor demonstrates good parallel scalability

Enrichment and stabilization strategies are effective

Abstract

Immersed finite element methods provide a convenient analysis framework for problems involving geometrically complex domains, such as those found in topology optimization and microstructures for engineered materials. However, their implementation remains a major challenge due to, among other things, the need to apply nontrivial stabilization schemes and generate custom quadrature rules. This article introduces the robust and computationally efficient algorithms and data structures comprising an immersed finite element preprocessing framework. The input to the preprocessor consists of a background mesh and one or more geometries defined on its domain. The output is structured into groups of elements with custom quadrature rules formatted such that common finite element assembly routines may be used without or with only minimal modifications. The key to the preprocessing framework is the construction of material topology information, concurrently with the generation of a quadrature rule, which is then used to perform enrichment and generate stabilization rules. While the algorithmic framework applies to a wide range of immersed finite element methods using different types of meshes, integration, and stabilization schemes, the preprocessor is presented within the context of the extended isogeometric analysis. This method utilizes a structured B-spline mesh, a generalized Heaviside enrichment strategy considering the material layout within individual basis functions$'$ supports, and face-oriented ghost stabilization. Using a set of examples, the effectiveness of the enrichment and stabilization strategies is demonstrated alongside the preprocessor$'$s robustness in geometric edge cases. Additionally, the performance and parallel scalability of the implementation are evaluated.