TY - JOUR
T1 - A sharp interface method using enriched finite elements for elliptic interface problems
AU - Hollbacher, Susanne
AU - Wittum, Gabriel
N1 - KAUST Repository Item: Exported on 2021-03-26
Acknowledgements: Open Access funding enabled and organized by Projekt DEAL.
PY - 2021/3/2
Y1 - 2021/3/2
N2 - We present an immersed boundary method for the solution of elliptic interface problems with discontinuous coefficients which provides a second-order approximation of the solution. The proposed method can be categorised as an extended or enriched finite element method. In contrast to other extended FEM approaches, the new shape functions get projected in order to satisfy the Kronecker-delta property with respect to the interface. The resulting combination of projection and restriction was already derived in Höllbacher and Wittum (TBA, 2019a) for application to particulate flows. The crucial benefits are the preservation of the symmetry and positive definiteness of the continuous bilinear operator. Besides, no additional stabilisation terms are necessary. Furthermore, since our enrichment can be interpreted as adaptive mesh refinement, the standard integration schemes can be applied on the cut elements. Finally, small cut elements do not impair the condition of the scheme and we propose a simple procedure to ensure good conditioning independent of the location of the interface. The stability and convergence of the solution will be proven and the numerical tests demonstrate optimal order of convergence.
AB - We present an immersed boundary method for the solution of elliptic interface problems with discontinuous coefficients which provides a second-order approximation of the solution. The proposed method can be categorised as an extended or enriched finite element method. In contrast to other extended FEM approaches, the new shape functions get projected in order to satisfy the Kronecker-delta property with respect to the interface. The resulting combination of projection and restriction was already derived in Höllbacher and Wittum (TBA, 2019a) for application to particulate flows. The crucial benefits are the preservation of the symmetry and positive definiteness of the continuous bilinear operator. Besides, no additional stabilisation terms are necessary. Furthermore, since our enrichment can be interpreted as adaptive mesh refinement, the standard integration schemes can be applied on the cut elements. Finally, small cut elements do not impair the condition of the scheme and we propose a simple procedure to ensure good conditioning independent of the location of the interface. The stability and convergence of the solution will be proven and the numerical tests demonstrate optimal order of convergence.
UR - http://hdl.handle.net/10754/668267
UR - http://link.springer.com/10.1007/s00211-021-01180-0
UR - http://www.scopus.com/inward/record.url?scp=85102068962&partnerID=8YFLogxK
U2 - 10.1007/s00211-021-01180-0
DO - 10.1007/s00211-021-01180-0
M3 - Article
SN - 0945-3245
JO - Numerische Mathematik
JF - Numerische Mathematik
ER -