English

Full Gradient Stabilized Cut Finite Element Methods for Surface Partial Differential Equations

Numerical Analysis 2016-08-24 v1

Abstract

We propose and analyze a new stabilized cut finite element method for the Laplace-Beltrami operator on a closed surface. The new stabilization term provides control of the full R3\mathbb{R}^3 gradient on the active mesh consisting of the elements that intersect the surface. Compared to face stabilization, based on controlling the jumps in the normal gradient across faces between elements in the active mesh, the full gradient stabilization is easier to implement and does not significantly increase the number of nonzero elements in the mass and stiffness matrices. The full gradient stabilization term may be combined with a variational formulation of the Laplace-Beltrami operator based on tangential or full gradients and we present a simple and unified analysis that covers both cases. The full gradient stabilization term gives rise to a consistency error which, however, is of optimal order for piecewise linear elements, and we obtain optimal order a priori error estimates in the energy and L2L^2 norms as well as an optimal bound of the condition number. Finally, we present detailed numerical examples where we in particular study the sensitivity of the condition number and error on the stabilization parameter.

Keywords

Cite

@article{arxiv.1602.01512,
  title  = {Full Gradient Stabilized Cut Finite Element Methods for Surface Partial Differential Equations},
  author = {Erik Burman and Peter Hansbo and Mats G. Larson and André Massing and Sara Zahedi},
  journal= {arXiv preprint arXiv:1602.01512},
  year   = {2016}
}

Comments

20 pages, 4 figures, 5 tables. arXiv admin note: text overlap with arXiv:1507.05835

R2 v1 2026-06-22T12:43:13.467Z