Max Wardetzky

Julian Eckhardt, Ralf Hiptmair, Thorsten Hohage et al.

By introducing a shape manifold as a solution set to solve inverse obstacle scattering problems we allow the reconstruction of general, not necessarily star-shaped curves. The bending energy is used as a stabilizing term in Tikhonov regularization to gain independence of the parametrization. Moreover, we discuss how self-intersections can be avoided by penalization with the Möbius energy and prove the regularizing property of our approach as well as convergence rates under variational source conditions. In the second part of the paper the discrete setting is introduced, and we describe a numerical method for finding the minimizer of the Tikhonov functional on a shape-manifold. Numerical examples demonstrate the feasibility of reconstructing non-star-shaped obstacles.

Henrik Schumacher, Max Wardetzky

Building on and extending tools from variational analysis, we prove Kuratowski convergence of sets of simplicial area minimizers to minimizers of the smooth Douglas-Plateau problem under simplicial refinement. This convergence is with respect to a topology that is stronger than uniform convergence of both positions and surface normals.

Stefan W. von Deylen, David Glickenstein, Max Wardetzky

We define barycentric coordinates on a Riemannian manifold using Karcher's center of mass technique applied to point masses for n+1 sufficiently close points, determining an n-dimensional Riemannian simplex defined as a "Karcher simplex." Specifically, a set of weights is mapped to the Riemannian center of mass for the corresponding point measures on the manifold with the given weights. If the points lie sufficiently close and in general position, this map is smooth and injective, giving a coordinate chart. We are then able to compute first and second derivative estimates of the coordinate chart. These estimates allow us to compare the Riemannian metric with the Euclidean metric induced on a simplex with edge lengths determined by the distances between the points. We show that these metrics differ by an error that shrinks quadratically with the maximum edge length. With such estimates, one can deduce convergence results for finite element approximations of problems on Riemannian manifolds.

Tim Brüers, Christoph Lehrenfeld, Tim van Beeck et al.

We present a discrete Helmholtz--Hodge decomposition for H(div)-conforming Brezzi--Douglas--Marini (BDM) finite elements on triangulated surfaces of arbitrary topology. The divergence-free BDM subspace is split L2-orthogonally into rotated gradients of a continuous streamfunction space and a finite-dimensional space of discrete harmonic fields whose dimension equals the first Betti number of the surface. Consequently, any incompressible flow discretized on this subspace can be reformulated with a scalar streamfunction and finitely many harmonic coefficients as the only unknowns. This eliminates the pressure and the saddle-point structure while ensuring exact tangentiality, pointwise divergence-freeness, and pressure-robustness. We present a randomized algorithm for constructing the harmonic basis and discuss implementation aspects including hybridization, efficient treatment of the harmonic unknowns, and pressure reconstruction. Numerical experiments for unsteady surface Navier--Stokes equations on a trefoil knot and a multiply-connected sculpture surface demonstrate the method and illustrate the physical role of the harmonic velocity component.