The Library
Unfitted finite element methods using bulk meshes for surface partial differential equations
Tools
Deckelnick, Klaus, Elliott, Charles M. and Ranner, Thomas (2014) Unfitted finite element methods using bulk meshes for surface partial differential equations. SIAM Journal on Numerical Analysis, Volume 52 (Number 4). pp. 2137-2162. doi:10.1137/130948641 ISSN 0036-1429.
Research output not available from this repository.
Request-a-Copy directly from author or use local Library Get it For Me service.
Official URL: http://dx.doi.org/10.1137/130948641
Abstract
In this paper, we define new unfitted finite element methods for numerically approximating the solution of surface partial differential equations 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, $\mathscr{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 $\mathscr{T}_h$ which is independent of $\Gamma$. Our first method is a sharp interface method (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 SIF from the method of [M. A. Olshanskii, A. Reusken, and J. Grande, SIAM J. Numer. Anal., 47 (2009), pp. 3339--3358]. The second method is a narrow band method (NBM) in which the region of integration is a narrow band of width $O(h)$. NBM is similar to the method of [K. Deckelnick et al., IMA J. Numer. Anal., 30 (2010), pp. 351--376] but again the full gradient is used in the discrete weak formulation. The a priori error analysis in this paper shows that the methods are of optimal order in the surface $L^{2}$ and $H^{1}$ norms and have the advantage that the normal derivative of the discrete solution is small and converges to zero. Our 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. We show that our method is conservative so that it preserves mass in the case of an advection-diffusion conservation law. Numerical results are given which illustrate the rates of convergence
Item Type: | Journal Article | ||||||||
---|---|---|---|---|---|---|---|---|---|
Divisions: | Faculty of Science, Engineering and Medicine > Science > Mathematics | ||||||||
Journal or Publication Title: | SIAM Journal on Numerical Analysis | ||||||||
Publisher: | Society for Industrial and Applied Mathematics | ||||||||
ISSN: | 0036-1429 | ||||||||
Official Date: | 19 August 2014 | ||||||||
Dates: |
|
||||||||
Volume: | Volume 52 | ||||||||
Number: | Number 4 | ||||||||
Page Range: | pp. 2137-2162 | ||||||||
DOI: | 10.1137/130948641 | ||||||||
Status: | Peer Reviewed | ||||||||
Publication Status: | Published | ||||||||
Access rights to Published version: | Restricted or Subscription Access | ||||||||
Adapted As: |
Request changes or add full text files to a record
Repository staff actions (login required)
View Item |