English

Finite element approximation of scalar curvature in arbitrary dimension

Numerical Analysis 2023-01-06 v1 Numerical Analysis Differential Geometry

Abstract

We analyze finite element discretizations of scalar curvature in dimension N2N \ge 2. Our analysis focuses on piecewise polynomial interpolants of a smooth Riemannian metric gg on a simplicial triangulation of a polyhedral domain ΩRN\Omega \subset \mathbb{R}^N having maximum element diameter hh. We show that if such an interpolant ghg_h has polynomial degree r0r \ge 0 and possesses single-valued tangential-tangential components on codimension-1 simplices, then it admits a natural notion of (densitized) scalar curvature that converges in the H2(Ω)H^{-2}(\Omega)-norm to the (densitized) scalar curvature of gg at a rate of O(hr+1)O(h^{r+1}) as h0h \to 0, provided that either N=2N = 2 or r1r \ge 1. As a special case, our result implies the convergence in H2(Ω)H^{-2}(\Omega) of the widely used "angle defect" approximation of Gaussian curvature on two-dimensional triangulations, without stringent assumptions on the interpolated metric ghg_h. We present numerical experiments that indicate that our analytical estimates are sharp.

Keywords

Cite

@article{arxiv.2301.02159,
  title  = {Finite element approximation of scalar curvature in arbitrary dimension},
  author = {Evan S. Gawlik and Michael Neunteufel},
  journal= {arXiv preprint arXiv:2301.02159},
  year   = {2023}
}