Linking ghost penalty and aggregated unfitted methods

Abstract

In this work, we analyse the links between ghost penalty stabilisation and aggregation-based discrete extension operators for the numerical approximation of elliptic partial differential equations on unfitted meshes. We explore the behavior of ghost penalty methods in the limit as the penalty parameter goes to infinity, which returns a strong version of these methods. We observe that these methods suffer locking in that limit. On the contrary, aggregated finite element spaces are locking-free because they can be expressed as an extension operator from well-posed to ill-posed degrees of freedom. Next, we propose novel ghost penalty methods that penalise the distance between the solution and its aggregation-based discrete extension. These methods are locking-free and converge to aggregated finite element methods in the infinite penalty parameter limit. We include an exhaustive set of numerical experiments in which we compare weak (ghost penalty) and strong (aggregated finite elements) schemes in terms of error quantities, condition numbers and sensitivity with respect to penalty coefficients on different geometries, intersection locations and mesh topologies.

References in corpus (16)

• The deal.II finite element library: design, features, and insights
• Condition number analysis and preconditioning of the finite cell method
• High order unfitted finite element methods on level set domains using isoparametric mappings
• The aggregated unfitted finite element method for elliptic problems
• Shape Optimization Using the Cut Finite Element Method
• A Nitsche-based cut finite element method for a fluid--structure interaction problem
• A note on the penalty parameter in Nitsche's method for unfitted boundary value problems
• CutIGA with Basis Function Removal
• Skeleton-stabilized ImmersoGeometric Analysis for incompressible viscous flow problems
• Mixed aggregated finite element methods for the unfitted discretization of the Stokes problem
• Robust and scalable h-adaptive aggregated unfitted finite elements for interface elliptic problems
• Distributed-memory parallelization of the aggregated unfitted finite element method
• Robust and scalable domain decomposition solvers for unfitted finite element methods
• Numerical modelling of heat transfer and experimental validation in Powder-Bed Fusion with the Virtual Domain Approximation
• The aggregated unfitted finite element method on parallel tree-based adaptive meshes
• A robust and scalable unfitted adaptive finite element framework for nonlinear solid mechanics

Cited by in corpus (3)

• An Arbitrarily High Order Unfitted Finite Element Method for Elliptic Interface Problems with Automatic Mesh Generation
• Residual-based error estimation and adaptivity for stabilized immersed isogeometric analysis using truncated hierarchical B-splines
• Unfitted Trefftz discontinuous Galerkin methods for elliptic boundary value problems