Unfitted finite element methods using bulk meshes for surface partial differential equations
type: Item
instance of: scholarly article
zbMATH Open document ID: 1307.65164
DOI: 10.1137/130948641
author: Charles M. Elliott
author: Thomas Ranner
author: Klaus Deckelnick
published in: SIAM Journal on Numerical Analysis
publication date: 2014-11-04T00:00:00Z
full work available at URL: https://arxiv.org/abs/1312.2905
Mathematics Subject Classification ID: 65N30
Mathematics Subject Classification ID: 65M60
Mathematics Subject Classification ID: 35R01
Mathematics Subject Classification ID: 65N15
Mathematics Subject Classification ID: 65M15
Mathematics Subject Classification ID: 58J05
Mathematics Subject Classification ID: 58J35
zbMATH DE Number: 6365796
zbMATH Keywords: elliptic and parabolic surface equations
zbMATH Keywords: unfitted finite elements
zbMATH Keywords: cut cells
zbMATH Keywords: error analysis
zbMATH Keywords: narrow band
zbMATH Keywords: sharp interface
zbMATH Keywords: bulk finite element
zbMATH Keywords: advection-diffusion conservation law
zbMATH Keywords: numerical result
Wikidata QID: Q59898944
describes a project that uses: TetGen
MaRDI profile type: Publication
OpenAlex ID: W1970875985
arXiv ID: 1312.2905
title: Unfitted finite element methods using bulk meshes for surface partial differential equations
review text: The authors consider the finite element approximation of surface partial differential equations. They define new unfitted finite element methods for numerically approximating the exact solution using bulk finite elements. The key idea is that the \(n\)-dimensional hypersurface, \(\Gamma \subset \mathbb{R}^{n+1}\), is embedded in a polyhedral domain in \(\mathbb R^{n+1}\), consisting of a union, \(\mathcal{T}_h\), of \((n+1)\)-simplices. The unifying feature of the methodological approach is that the finite element approximating space is based on continuous piecewise linear finite element functions on the bulk triangulation \(\mathcal{T}_h\) which is independent of \(\Gamma\). The first method is a sharp interface method (\textit{SIF}) which uses the bulk finite element space in an approximating weak formulation obtained from integration on a polygonal approximation, \(\Gamma_{h}\), of \(\Gamma\). The full gradient is used rather than the projected tangential gradient and it is this which distinguishes \textit{SIF} from some previous methods. The second method is a narrow band method (\textit{NBM}) in which the region of integration is a narrow band of width \(O(h)\). \textit{NBM} is similar to the method of \textit{K. Deckelnick} et al. [IMA J. Numer. Anal. 30, No. 2, 351--376 (2010; Zbl 1191.65151)], but again, the full gradient is used in the discrete weak formulation. An a priori error analysis is presented and it shows that the methods are of optimal order in the surface \(L^{2}\) and \(H^{1}\) norms. The third method combines bulk finite elements, discrete sharp interfaces, and narrow bands in order to give an unfitted finite element method for parabolic equations on evolving surfaces. It is shown that the method is conservative so that it preserves mass in the case of an advection-diffusion conservation law. Some numerical results are presented to illustrate the rates of convergence.
reviewed by: Abdallah Bradji
description: scientific article; zbMATH DE number 6365796