An EigenValue Stabilization Technique for Immersed Boundary Finite Element Methods in Explicit Dynamics

Abstract

The application of immersed boundary methods in static analyses is often impeded by poorly cut elements (small cut elements problem), leading to ill-conditioned linear systems of equations and stability problems. While these concerns may not be paramount in explicit dynamics, a substantial reduction in the critical time step size based on the smallest volume fraction $\chi$ of a cut element is observed. This reduction can be so drastic that it renders explicit time integration schemes impractical. To tackle this challenge, we propose the use of a dedicated eigenvalue stabilization (EVS) technique. The EVS-technique serves a dual purpose. Beyond merely improving the condition number of system matrices, it plays a pivotal role in extending the critical time increment, effectively broadening the stability region in explicit dynamics. As a result, our approach enables robust and efficient analyses of high-frequency transient problems using immersed boundary methods. A key advantage of the stabilization method lies in the fact that only element-level operations are required. This is accomplished by computing all eigenvalues of the element matrices and subsequently introducing a stabilization term that mitigates the adverse effects of cutting. Notably, the stabilization of the mass matrix ${\mathbf{M}}_{\mathrm{c}}$ of cut elements -- especially for high polynomial orders $p$ of the shape functions -- leads to a significant raise in the critical time step size $\Delta t_{\mathrm{cr}}$. To demonstrate the efficacy of our technique, we present two specifically selected dynamic benchmark examples related to wave propagation analysis, where an explicit time integration scheme must be employed to leverage the increase in the critical time step size.

arXiv comment45 pages, 25 figures

References in corpus (13)

• Stable Generalized Finite Element Method (SGFEM)
• Condition number analysis and preconditioning of the finite cell method
• Robust and parallel scalable iterative solutions for large-scale finite cell analyses
• Scalable multigrid methods for immersed finite element methods and immersed isogeometric analysis
• A mathematical theory for mass lumping and its generalization with applications to isogeometric analysis
• High-order implicit time integration scheme based on Padé expansions
• Towards higher-order accurate mass lumping in explicit isogeometric analysis for structural dynamics
• Moment fitted cut spectral elements for explicit analysis of guided wave propagation
• The aggregated unfitted finite element method on parallel tree-based adaptive meshes
• Critical time-step size analysis and mass scaling by ghost-penalty for immersogeometric explicit dynamics
• Hierarchical multigrid approaches for the finite cell method on uniform and multi-level hp-refined grids
• High-order implicit time integration scheme with controllable numerical dissipation based on mixed-order Padé expansions
• Implicit-Explicit Time Integration for the Immersed Wave Equation

Cited by in corpus (4)

• On the efficiency of explicit and semi-explicit immersed boundary finite element methods for wave propagation problems
• A CFL condition for the finite cell method
• A Non-intrusive Approach for the Imposition of Strong Dirichlet Boundary Conditions in Unfitted Boundary Meshes
• Generalized eigenvalue stabilization for immersed explicit dynamics