Computing confined elasticae

Abstract

AbstractWe devise a numerical scheme for computing arc-length parameterized curves of low bending energy that are confined to convex domains. We address the convergence of the discrete formulations to a continuous model and the unconditional stability of an iterative scheme. Numerical simulations confirm the theoretical results and lead to a classification of observed optimal curves within spheres.