Variational Convergence of Discrete Elasticae

Sebastian Scholtes,Henrik Schumacher,Max Wardetzky
DOI: https://doi.org/10.48550/arXiv.1901.02228
2019-01-08
Abstract:We discuss a discretization by polygonal lines of the Euler-Bernoulli bending energy and of Euler elasticae under clamped boundary conditions. We show Hausdorff convergence of the set of almost minimizers of the discrete bending energy to the set of smooth Euler elasticae under mesh refinement in (i) the $W^{1,\infty}$-topology for piecewise-linear interpolation and in (ii) the $W^{2,p}$-topology, $p \in{[2,\infty[}$, using a suitable smoothing operator to create $W^{2,p}$-curves from polygons.
Numerical Analysis
What problem does this paper attempt to address?
The problem that this paper attempts to solve is about how discrete elasticae approximate their corresponding smooth Euler elasticae under specific boundary conditions. Specifically, the authors focus on the polygonal line discretization of the Euler - Bernoulli bending energy and explore how the almost - minimizing sets of the discrete bending energy converge to the sets of smooth Euler elasticae during the process of mesh refinement. They prove that in the \(W^{1,\infty}\) topology and the \(W^{2,p}\) topology (where \(p\in[2,\infty[\)), by generating \(W^{2,p}\) curves from polygons through an appropriate smoothing operator, the almost - minimizing sets of the discrete bending energy converge to the sets of smooth Euler elasticae in the Hausdorff distance. ### Specific Problem Description 1. **Definition and Background**: - The Euler - Bernoulli bending energy is used to model the elastic stored energy in the bent part of a thin, flexible but inextensible material. - For a compact interval \(\Sigma\subset\mathbb{R}\) and a curve \(\gamma\in W^{2,2}(\Sigma;\mathbb{R}^m)\), the Euler - Bernoulli bending energy is defined as the integral of the square of the curvature: \[ E(\gamma)=\frac{1}{2}\int_{\Sigma}|\kappa_{\gamma}|^2\omega_{\gamma} \] - The classical Euler elasticae problem is to find the energy - minimizing curve under a given fixed curve length \(L\) and fixed first - order boundary conditions (i.e., the positions and tangent vectors at both ends). 2. **Main Challenges**: - The fixed curve length and fixed boundary conditions are the main difficulties of the problem. - To deal with these constraints, the paper proposes a polygonal model to approximate the smooth curve and studies the convergence of these discrete models in different topological spaces. 3. **Research Objectives**: - Prove how the discrete elasticae converge to the smooth elasticae during the mesh refinement process. - Ensure that during the convergence process, the discrete model can maintain the accuracy and stability of the energy. ### Solution The authors solve this problem through the following steps: - **Discretization Method**: Introduce the polygonal line discretization method and define the discrete Euler - Bernoulli energy: \[ E_T(P)=\frac{1}{2}\sum_{i\in V_{\text{int}}(T)}\left(\frac{\alpha_P(i)}{\bar{\ell}_P(i)}\right)^2\bar{\ell}_P(i) \] where \(\alpha_P(i)\) is the turning angle at the internal vertex \(i\) and \(\bar{\ell}_P(i)\) is the double - edge length. - **Convergence Analysis**: Prove that in the \(W^{1,\infty}\) and \(W^{2,p}\) topologies, the almost - minimizing sets of the discrete bending energy converge to the sets of smooth Euler elasticae in the Hausdorff distance. - **Key Results**: - The Hausdorff convergence theorem (Theorem 1.1) shows that, under a given curve length \(L\) and compatible boundary conditions, for each \(p\in[2,\infty[\), there exists a constant \(C\geq0\) such that: \[ \lim_{h(T)\to0}\text{dist}_{W^{2,p}}(M,R_T(A_T\cap M_{Ch(T)})) = 0 \] and \[ \lim_{h(T)\to0}\text{dist}_{W^{1,\infty}}