Patrick Dondl

Patrick Dondl, Behrend Heeren, Martin Rumpf

Branching can be observed at the austenite-martensite interface of martensitic phase transformations. For a model problem, Kohn and Müller studied a branching pattern with optimal scaling of the energy with respect to its parameters. Here, we present finite element simulations that suggest a topologically different class of branching patterns and derive a novel, low dimensional family of patterns. After a geometric optimization within this family, the resulting pattern bears a striking resemblance to our simulation. The novel microstructure admits the same scaling exponents but results in a significantly lower upper energy bound.

Patrick Dondl, Stephan Wojtowytsch

We describe the implementation of a topological constraint in finite element simulations of phase field models which ensures path-connectedness of preimages of intervals in the phase field variable. Two main applications of our method are presented. First, a discrete steepest decent of a phase field version of a bending energy with spontaneous curvature and additional surface area penalty is shown, which leads to disconnected surfaces without our topological constraint but connected surfaces with the constraint. The second application is the segmentation of an image into a connected component and its exterior. Numerically, our constraint is treated using a suitable geodesic distance function which is computed using Dijkstra's algorithm.

Patrick Dondl, Patrina S. P. Poh, Martin Rumpf et al.

This paper deals with the simulateneous optimization of a subset $\mathcal{O}_0$ of some domain $Ω$ and its complement $\mathcal{O}_1 = Ω\setminus \overline{\mathcal{O}}_0$ both considered as separate elastic objects subject to a set of loading scenarios. If one asks for a configuration which minimizes the maximal elastic cost functional both phases compete for space since elastic shapes usually get mechanically more stable when being enlarged. Such a problem arises in biomechanics where a bioresorbable polymer scaffold is implanted in place of lost bone tissue and in a regeneration phase new bone tissue grows in the scaffold complement via osteogenesis. In fact, the polymer scaffold should be mechanically stable to bear loading in the early stage regeneration phase and at the same time the new bone tissue grown in the complement of this scaffold should as well bear the loading. Here, this optimal subdomain splitting problem with appropriate elastic cost functionals is introduced and existence of optimal two phase configurations is established for a regularized formulation. Furthermore, based on a phase field approximation a finite element discretization is derived. Numerical experiments are presented for the design of optimal periodic scaffold microstructure.