Consistent discretization of higher-order interface models for thin layers and elastic material surfaces, enabled by isogeometric cut-cell methods

Abstract

Many interface formulations, e.g. based on asymptotic thin interphase models or material surface theories, involve higher-order differential operators and discontinuous solution fields. In this article, we are taking first steps towards a variationally consistent discretization framework that naturally accommodates these two challenges by synergistically combining recent developments in isogeometric analysis and cut-cell finite element methods. Its basis is the mixed variational formulation of the elastic interface problem that provides access to jumps in displacements and stresses for incorporating general interface conditions. Upon discretization with smooth splines, derivatives of arbitrary order can be consistently evaluated, while cut-cell meshes enable discontinuous solutions at potentially complex interfaces. We demonstrate via numerical tests for three specific nontrivial interfaces (two regimes of the Benveniste-Miloh classification of thin layers and the Gurtin-Murdoch material surface model) that our framework is geometrically flexible and provides optimal higher-order accuracy in the bulk and at the interface.

A computational approach based on extended finite element method for thin porous layers in acoustic problems

“…These methods allow meshing to be independent of the interface geometry, leading to flexible approaches. Burman et al 13 discussed the implementation of cut‐FEM for first‐order discontinuity in various problems, and Han et al 14 employed isogeometric cut‐cell methods for high‐order interface models. Yvonnet et al, 15 Zhu et al, 16 and Benvenuti et al 17 used the extended finite element method (X‐FEM) to implement cohesive interface and spring‐layer models in composites.…”

An ALE method for simulations of axisymmetric elastic surfaces in flow

“…Main problems of such methods are the loss in accuracy associated to this interpolation between grid and the handling of high viscosity ratios [26,27]. These problems can be overcome by coupling immersed boundaries with technically far more complicated cut cell methods [28].…”