Enriched immersed finite element and isogeometric analysis: algorithms and data structures

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,...

Full description

Saved in:
Bibliographic Details
Published in:Engineering with computers
Main Authors: Wunsch, Nils, Doble, Keenan, Schmidt, Mathias R., Noël, Lise, Evans, John A., Maute, Kurt
Format: Journal Article
Language:English
Published: United States Springer Science and Business Media LLC 18.08.2025
Subjects:
Computer implementation
Ghost stabilization
Heaviside enrichment
Immersed finite element method
Multi-material problems
XIGA
ISSN:0177-0667, 1435-5663, 1435-5663
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary: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.
Bibliography:AC52-07NA27344
LLNL--JRNL-871539
National Science Foundation (NSF)
Defense Advanced Research Projects Agency (DARPA)
USDOE National Nuclear Security Administration (NNSA)
US Air Force Office of Scientific Research (AFOSR)
ISSN:0177-0667
1435-5663
1435-5663
DOI:10.1007/s00366-025-02163-7