Ausgewählte Vorträge

S. Bartels

Babuska's paradox in linear and nonlinear bending problems

ICERM Programme on Numerical PDEs
ICERM, Brown University, Providence, RI
March 19, 2024

Abstract. The plate bending or Babuska paradox refers to the failure of convergence when a linear bending problem with simple support boundary conditions is approximated using polygonal domain approximations. We provide an explanation based on a variational viewpoint and identify sufficient conditions that avoid the paradox and which show that boundary conditions have to be suitably modified. We show that the paradox also matters in nonlinear thin-sheet folding problems and devise approximations that correctly converge to the original problem.

S. Bartels

Paper Folding: Modeling and Simulation

SFB 1481 Sparsity and Singular Structures
Kickoff Workshop, Aachen
September 22, 2022

Abstract. The folding of a thin elastic sheet along a curved arc has various applications including the construction of bistable devices. We discuss the derivation of a plate model from three-dimensional hyperelasticity and rigidity properties of admissible deformations and minimizers. The numerical solution is based on an isoparametric discontinuous Galerkin finite element method that provides a suitable geometric approximation of the folding arc. Error estimates are presented for a linearized version of the model problem. The talk is based on joint work with Andrea Bonito (Texas A&M), Peter Hornung (U Dresden), and Philipp Tscherner (U Freiburg).

S. Bartels

Simulation of Free Boundary Problems in the Nonlinear Bending of Elastic Rods and Plates

15th International Conference on Free Boundary Problems
Online Conference, Berlin
September 15, 2021

Abstract. The mathematical description of large bending deformations of thin elastic rods and plates leads to fourth order problems with nonlinear pointwise constraints that give rise to various free boundary phenomena. The free boundary may separate regions of trivial and large deformations, describe the contact zone in the presence of an obstacle, or be related to the occurrence of self-contact. We devise and analyse numerical methods that are capable of reliably capturing these effects under minimal regularity assumptions and which allow us to experimentally study topological transitions of free boundaries.

S. Bartels et al.

Error estimates for nonconforming and discontinuous discretizations of nonsmooth problems via convex duality

Sayas Numerics Seminar
Joint online seminar by UMD College Park & Baltimore, George Mason U, U Delaware
September 7, 2021

Abstract. Nonconforming and discontinuous finite elements are attractive for discretising variational problems with limited regularity properties, in particular, when discontinuitiesmay occur. The possible lack of differentiability of related functionals prohibits the use of classical arguments to derive error estimates. As an alternative we make use of discrete andcontinuous convex duality relations and of quasi-interpolation operators with suitable projection properties. For total variation regularised minimisation problems a quasi-optimal error estimate is derived which is not available for standard finite element methods. Our results use and extend recent ideas by Chambolle and Pock.

S. Bartels and P. Reiter

Approximating gradient flow evolutions of self-avoiding inextensible curves and elastic knots

Thematic Programme on Numerical analysis of complex PDE models in the sciences
Erwin-Schrödinger-Institut, Vienna, Austria, July 18, 2018

Abstract. We discuss a semi-implict scheme that allows for minimizing the bending energy of curves within certain isotopy classes. To this end we consider a weighted sum of a bending energy and a tangent-point functional. The latter is defined via the iterated integral of a negative power of the tangent-point radius which is for a given pair of points the radius of the circle that is tangent to the curve in one point and intersects the curve in the other point. Evolutions are defined via gradient flows for the total energy within classes of arclength-parametrized curves. The numerical scheme is obtained by a semi-implicit discretization which treats the nonlinear and nonlocal tangent-point functional explicitly. This avoids the occurrence of fully populated matrices and allows for a straightforward parallelization of its computation. Based on estimates for the second derivative of the tangent-point functional and a uniform bi-Lipschitz radius, we prove a stability result implying energy decay during the evolution as well as maintenance of arclength parametrization. We present some numerical experiments exploring the energy landscape, targeted to the question how to obtain global minimizers of the bending energy in knot classes, so-called elastic knots.

S. Bartels et al.

Numerical solution of nonsmooth problems and application to damage evolution and optimal insulation

Foundations of Computational Mathematics
Barcelona, Spain, July 17, 2017

Abstract. Nonsmooth minimization problems and singular partial differential equations arise in the description of inelastic material behavior, image processing, and modeling of non-Newtonian fluids. The numerical discretization and iterative solution is often based on regularizing or stabilizing terms. In this talk we address the influence of such modifications on error estimates and the robustness of iterative solution methods. In particular, we present an unconditional stability result for semi-implicit discretizations of a class of singular flows and devise a variant of the alternating direction of multipliers method with variable step sizes. The results and methods are illustrated by numerical experiments for a damage evolution model and a problem of optimal insulation leading to a break of symmetry.

S. Bartels, P. Reiter, J. Riege

Approximation of self-avoiding inextensible curves

8th Workshop on Numerical Methods for Evolution Equations
Heraklion, Crete, September 23, 2016

Abstract. The prediction of the behavior of elastic curves arises in knot theory and the modeling of filament networks in biological cells. When large deformations are considered or when one is interested in determining homotopy classes, injectivity of deformations has to be imposed. A practical realization consists in including an appropriate potential in the energy functional. Its numerical treatment is difficult due to the singular and nonlocal nature. We present a discretization of a tangent-point potential that leads to a simple and stable scheme and which allows for the use of large step sizes.

S. Bartels, P. Schön

Adaptive approximation of the Monge-Kantorovich problem

Mafelap 2016
Brunel, UK, June 15, 2016

Abstract. Optimal transportation problems define high-dimensional linear programs. An efficient approach to their numerical solution is based on reformulations as nonlinear partial differential equations. If transportation cost is proportional to distance this leads to the Monge--Kantorovich problem which is a constrained minimization problem on Lipschitz continuous functions. We discuss the iterative solution via splitting methods and devise an adaptive mesh refinement strategy based on an a~posteriori error estimate for the primal-dual gap.