Robert Nürnberg

We consider a sharp interface formulation for an anisotropic multi-phase Mullins–Sekerka problem with kinetic undercooling. The flow is characterized by a cluster of surfaces evolving such that the total surface energy plus a weighted sum of the volumes of the enclosed phases decreases in time. Upon deriving a suitable variational formulation, we introduce a fully discrete unfitted finite element method. In this approach, the approximations of the moving interfaces are independent of the triangulations used for the equations in the bulk. Our method can be shown to be unconditionally stable. Several numerical examples demonstrate the capabilities of the proposed method. In particular, it is demonstrated that the evolution of multiple ice crystals with junctions can be modeled using the proposed approach.

Mullins–Sekerka models with chemical reactions can lead to scenarios where droplets grow, become unstable, split, grow and undergo further division. These grow and division cycles have been proposed as a model for protocells and are believed to play a fundamental role in living systems by providing chemical compartments which are important in the organization of living systems. This paper analyzes chemically active Mullins–Sekerka models. Existence of radially symmetric solutions is shown and a detailed stability analysis in radial as well as planar situations is given. In particular, we also analyze multilayered solutions leading to shell-type situations. Finally, we introduce a numerical method based on a parametric finite element approach that explicitly accounts for topological changes, thereby allowing for droplet splitting and merging. Several numerical simulations verify the findings of the theoretical stability analysis and show complex dynamical behavior, including multiple instabilities, splittings of droplets and appearance of shell-type solutions.

In this article we investigate the question of finding a network configuration of minimal length connecting three given points in the Heisenberg group. After proving existence of (possibly degenerate) minimal horizontal triods, we investigate their characterization. We then formulate a horizontal curve shortening flow that deforms any given suitable initial triod into a critical point for the length functional. Numerical experiments based on a stable fully discrete finite element scheme provide useful insights into the rich landscape of this sub-Riemannian geometry.

We consider a sharp interface formulation for an anisotropic multi-phase Mullins-Sekerka problem with kinetic undercooling. The flow is characterized by a cluster of surfaces evolving such that the total surface energy plus a weighted sum of the volumes of the enclosed phases decreases in time. Upon deriving a suitable variational formulation, we introduce a fully discrete unfitted finite element method. In this approach, the approximations of the moving interfaces are independent of the triangulations used for the equations in the bulk. Our method can be shown to be unconditionally stable. Several numerical examples demonstrate the capabilities of the introduced method. In particular, it is demonstrated that the evolution of multiple ice crystals with junctions can be modeled using the proposed approach.

We prove optimal error bounds for a second order in time finite element approximation of curve shortening flow in possibly higher codimension. In addition, we introduce a second order in time method for curve diffusion. Both schemes are based on variational formulations of strictly parabolic systems of partial differential equations that feature a tangential velocity which under discretization is beneficial for the mesh quality. In each time step only two linear systems need to be solved. Numerical experiments demonstrate second order convergence as well as asymptotic equidistribution.

Mullins-Sekerka models with chemical reactions can lead to scenarios where droplets grow, become unstable, split, grow and undergo further division. These grow and division cycles have been proposed as a model for protocells and are believed to play a fundamental role in living systems by providing chemical compartments which are important in the organization of living systems. This paper analyses chemically active Mullins-Sekerka models. Existence of radially symmetric solutions is shown and a detailed stability analysis in radial as well as planar situations is given. In particular, we also analyze multilayered solutions leading to shell-type situations. Finally, we introduce a numerical method based on a parametric finite element approach that explicitly accounts for topological changes, thereby allowing for droplet splitting and merging. Several numerical simulations verify the findings of the theoretical stability analysis and show complex dynamical behavior, including multiple instabilities, splittings of droplets and appearance of shell-type solutions.

The Cahn–Hilliard model with reaction terms can lead to situations in which no coarsening is taking place and, in contrast, growth and division of droplets occur which all do not grow larger than a certain size. This phenomenon has been suggested as a model for protocells, and a model based on the modified Cahn–Hilliard equation has been formulated. We introduce this equation and show the existence and uniqueness of solutions. Then, formally matched asymptotic expansions are used to identify a sharp interface limit using a scaling of the reaction term, which becomes singular when the interfacial thickness tends to zero. We compute planar solutions and study their stability under non-planar perturbations. Numerical computations for the suggested model are used to validate the sharp interface asymptotics. In addition, the numerical simulations show that the reaction terms lead to diverse phenomena such as growth and division of droplets in the obtained solutions, as well as the formation of shell-like structures.