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.
翻译:在这项工作中,我们分析幽灵惩罚稳定化和基于聚合的离散扩展操作器之间的联系,以了解对不合适的 meshes 的椭圆部分偏差方程式的数字近似值。我们探讨当惩罚参数到达无限度时,幽灵惩罚方法在限值内的行为,这返回了这些方法的强大版本。我们观察到,这些方法被锁定在这个限度内。相反,合计的有限要素空间是无锁的,因为它们可以表现为从充分允许到错误自由度的扩展操作器。接下来,我们提出新的幽灵惩罚方法,以惩罚溶液与其基于聚合的离散扩展之间的距离。这些方法是无锁的,与无限惩罚参数限制中的总定值要素方法趋同。我们包括一系列详尽的数字实验,在差数数量、条件和敏感度方面将弱(宿主惩罚)和强(集中的有限要素)方法与不同地理特征、交叉地点和类比。