English

Cut Finite Element Methods for Partial Differential Equations on Embedded Manifolds of Arbitrary Codimensions

Numerical Analysis 2016-10-07 v1

Abstract

We develop a theoretical framework for the analysis of stabilized cut finite element methods for the Laplace-Beltrami operator on a manifold embedded in Rd\mathbb{R}^d of arbitrary codimension. The method is based on using continuous piecewise polynomials on a background mesh in the embedding space for approximation together with a stabilizing form that ensures that the resulting problem is stable. The discrete manifold is represented using a triangulation which does not match the background mesh and does not need to be shape-regular, which includes level set descriptions of codimension one manifolds and the non-matching embedding of independently triangulated manifolds as special cases. We identify abstract key assumptions on the stabilizing form which allow us to prove a bound on the condition number of the stiffness matrix and optimal order a priori estimates. The key assumptions are verified for three different realizations of the stabilizing form including a novel stabilization approach based on penalizing the surface normal gradient on the background mesh. Finally, we present numerical results illustrating our results for a curve and a surface embedded in R3\mathbb{R}^3.

Keywords

Cite

@article{arxiv.1610.01660,
  title  = {Cut Finite Element Methods for Partial Differential Equations on Embedded Manifolds of Arbitrary Codimensions},
  author = {Erik Burman and Peter Hansbo and Mats G. Larson and Andre Massing},
  journal= {arXiv preprint arXiv:1610.01660},
  year   = {2016}
}

Comments

30 pages, 4 figures, 3 tables

R2 v1 2026-06-22T16:12:29.834Z