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

Fuente: arXiv
Enregistré dans:
Détails bibliographiques
Auteurs principaux: Wunsch, Nils, Doble, Keenan, Schmidt, Mathias R., Noël, Lise, Evans, John A., Maute, Kurt
Format: Preprint
Publié: 2025
Sujets:
Accès en ligne:
Tags: Ajouter un tag
Pas de tags, Soyez le premier à ajouter un tag!
_version_ 1866929691035697152
author Wunsch, Nils
Doble, Keenan
Schmidt, Mathias R.
Noël, Lise
Evans, John A.
Maute, Kurt
author_facet Wunsch, Nils
Doble, Keenan
Schmidt, Mathias R.
Noël, Lise
Evans, John A.
Maute, Kurt
contents 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.
format Preprint
id arxiv_https___arxiv_org_abs_2501_17853
institution arXiv
publishDate 2025
record_format arxiv
spellingShingle Enriched Immersed Finite Element and Isogeometric Analysis -- Algorithms and Data Structures
Wunsch, Nils
Doble, Keenan
Schmidt, Mathias R.
Noël, Lise
Evans, John A.
Maute, Kurt
Numerical Analysis
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.
title Enriched Immersed Finite Element and Isogeometric Analysis -- Algorithms and Data Structures
topic Numerical Analysis
url https://arxiv.org/abs/2501.17853