Critical time-step size analysis and mass scaling by ghost-penalty for immersogeometric explicit dynamics

Abstract

In this article, we study the effect of small-cut elements on the critical time-step size in an immersogeometric context. We analyze different formulations for second-order (membrane) and fourth-order (shell-type) equations, and derive scaling relations between the critical time-step size and the cut-element size for various types of cuts. In particular, we focus on different approaches for the weak imposition of Dirichlet conditions: by penalty enforcement and with Nitsche's method. The stability requirement for Nitsche's method necessitates either a cut-size dependent penalty parameter, or an additional ghost-penalty stabilization term is necessary. Our findings show that both techniques suffer from cut-size dependent critical time-step sizes, but the addition of a ghost-penalty term to the mass matrix serves to mitigate this issue. We confirm that this form of `mass-scaling' does not adversely affect error and convergence characteristics for a transient membrane example, and has the potential to increase the critical time-step size by orders of magnitude. Finally, for a prototypical simulation of a Kirchhoff-Love shell, our stabilized Nitsche formulation reduces the solution error by well over an order of magnitude compared to a penalty formulation at equal time-step size.

References in corpus (12)

• Condition number analysis and preconditioning of the finite cell method
• The aggregated unfitted finite element method for elliptic problems
• Weak Boundary Condition Enforcement for Linear Kirchhoff-Love Shells: Formulation, Error Analysis, and Verification
• A note on the penalty parameter in Nitsche's method for unfitted boundary value problems
• A boundary penalization technique to remove outliers from isogeometric analysis on tensor-product meshes
• Spectral approximation properties of isogeometric analysis with variable continuity
• Cut Finite Element Methods for Elliptic Problems on Multipatch Parametric Surfaces
• Coupling of non-conforming trimmed isogeometric Kirchhoff-Love shells via a projected super-penalty approach
• Residual-based error estimation and adaptivity for stabilized immersed isogeometric analysis using truncated hierarchical B-splines
• Variationally consistent mass scaling for explicit time-integration schemes of lower- and higher-order finite element methods
• A variational approach based on perturbed eigenvalue analysis for improving spectral properties of isogeometric multipatch discretizations
• Conservative Discontinuous Cut Finite Element Methods

Cited by in corpus (5)

• An EigenValue Stabilization Technique for Immersed Boundary Finite Element Methods in Explicit Dynamics
• Isogeometric Multi-Resolution Full Waveform Inversion based on the Finite Cell Method
• Implicit-Explicit Time Integration for the Immersed Wave Equation
• A theoretical study on the effect of mass lumping on the discrete frequencies in immersogeometric analysis
• A CFL condition for the finite cell method