Enriched immersed finite element and isogeometric analysis

algorithms and data structures

Journal Article (2025)
Author(s)

Nils Wunsch (University of Colorado Boulder)

K. Doble (University of Colorado Boulder)

Mathias R. Schmidt (Lawrence Livermore National Laboratory)

Lise Noel (TU Delft - Computational Design and Mechanics)

John A. Evans (University of Colorado Boulder)

Kurt Maute (University of Colorado Boulder)

Research Group
Computational Design and Mechanics
DOI related publication
https://doi.org/10.1007/s00366-025-02163-7
More Info
expand_more
Publication Year
2025
Language
English
Research Group
Computational Design and Mechanics
Reuse Rights

Other than for strictly personal use, it is not permitted to download, forward or distribute the text or part of it, without the consent of the author(s) and/or copyright holder(s), unless the work is under an open content license such as Creative Commons.

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.