Stable approximations for axisymmetric Willmore flow for closed and open surfaces

For a hypersurface in ${\mathbb R}^3$, Willmore flow is defined as the $L^2$--gradient flow of the classical Willmore energy: the integral of the squared mean curvature. This geometric evolution law is of interest in differential geometry, image reconstruction and mathematical biology. In this paper, we propose novel numerical approximations for the Willmore flow of axisymmetric hypersurfaces. For the semidiscrete continuous-in-time variants we prove a stability result. We consider both closed surfaces, and surfaces with a boundary. In the latter case, we carefully derive weak formulations of suitable boundary conditions. Furthermore, we consider many generalizations of the classical Willmore energy, particularly those that play a role in the study of biomembranes. In the generalized models we include spontaneous curvature and area difference elasticity (ADE) effects, Gaussian curvature and line energy contributions. Several numerical experiments demonstrate the efficiency and robustness of our developed numerical methods.Comment: 52 pages, 19 figure

Numerical approximation of boundary value problems for curvature flow and elastic flow in Riemannian manifolds

We present variational approximations of boundary value problems for curvature flow (curve shortening flow) and elastic flow (curve straightening flow) in two-dimensional Riemannian manifolds that are conformally flat. For the evolving open curves we propose natural boundary conditions that respect the appropriate gradient flow structure. Based on suitable weak formulations we introduce finite element approximations using piecewise linear elements. For some of the schemes a stability result can be shown. The derived schemes can be employed in very different contexts. For example, we apply the schemes to the Angenent metric in order to numerically compute rotationally symmetric self-shrinkers for the mean curvature flow. Furthermore, we utilise the schemes to compute geodesics that are relevant for optimal interface profiles in multi-component phase field models.Comment: 42 pages, 21 figure

Finite element methods for fourth order axisymmetric geometric evolution equations

Fourth order curvature driven interface evolution equations frequently appear in the natural sciences. Often axisymmetric geometries are of interest, and in this situation numerical computations are much more efficient. We will introduce and analyze several new finite element schemes for fourth order geometric evolution equations in an axisymmetric setting, and for selected schemes we will show existence, uniqueness and stability results. The presented schemes have very good mesh and stability properties, as will be demonstrated by several numerical examples

Stable approximations for axisymmetric Willmore flow for closed and open surfaces

For a hypersurface in ℝ3, Willmore flow is defined as the L2-gradient flow of the classical Willmore energy: the integral of the squared mean curvature. This geometric evolution law is of interest in differential geometry, image reconstruction and mathematical biology. In this paper, we propose novel numerical approximations for the Willmore flow of axisymmetric hypersurfaces. For the semidiscrete continuous-in-time variants we prove a stability result. We consider both closed surfaces, and surfaces with a boundary. In the latter case, we carefully derive weak formulations of suitable boundary conditions. Furthermore, we consider many generalizations of the classical Willmore energy, particularly those that play a role in the study of biomembranes. In the generalized models we include spontaneous curvature and area difference elasticity (ADE) effects, Gaussian curvature and line energy contributions. Several numerical experiments demonstrate the efficiency and robustness of our developed numerical methods

Stable approximations for axisymmetric Willmore flow for closed and open surfaces

For a hypersurface in Double-struck capital R-3, Willmore flow is defined as the L-2-gradient flow of the classical Willmore energy: the integral of the squared mean curvature. This geometric evolution law is of interest in differential geometry, image reconstruction and mathematical biology. In this paper, we propose novel numerical approximations for the Willmore flow of axisymmetric hypersurfaces. For the semidiscrete continuous-in-time variants we prove a stability result. We consider both closed surfaces, and surfaces with a boundary. In the latter case, we carefully derive weak formulations of suitable boundary conditions. Furthermore, we consider many generalizations of the classical Willmore energy, particularly those that play a role in the study of biomembranes. In the generalized models we include spontaneous curvature and area difference elasticity (ADE) effects, Gaussian curvature and line energy contributions. Several numerical experiments demonstrate the efficiency and robustness of our developed numerical methods

Stable approximations for axisymmetric Willmore flow for closed and open surfaces

For a hypersurface in ${\mathbb R}^3$, Willmore flow is defined as the $L^2$--gradient flow of the classical Willmore energy: the integral of the squared mean curvature. This geometric evolution law is of interest in differential geometry, image reconstruction and mathematical biology. In this paper, we propose novel numerical approximations for the Willmore flow of axisymmetric hypersurfaces. For the semidiscrete continuous-in-time variants we prove a stability result. We consider both closed surfaces, and surfaces with a boundary. In the latter case, we carefully derive suitable boundary conditions. Furthermore, we consider many generalizations of the classical Willmore energy, particularly those that play a role in the study of biomembranes. In the generalized models we include spontaneous curvature and area difference elasticity (ADE) effects, Gaussian curvature and line energy contributions. Several numerical experiments demonstrate the efficiency and robustness of our developed numerical methods

Point particle interactions on surface biomembranes: second order splitting and surface finite elements

We study the well-posedness and approximation of mathematical models for small deformations of biological membranes where the deformations are due to point constraints. The differentiability of of the membrane energy with respect to the movement of the point constraints is studied. We begin by reviewing mathematical theory related to the shape of biomembranes and embedded proteins. We show that modifications of established theory hold and introduce notation which allows us to easily discuss the movement of many proteins embedded into the surface. We then discuss the well-posedness of an abstract second order splitting method with linear constraints, which we will apply to the energy minimising biomembrane with embedded proteins. We also consider a penalty method to weakly enforce the constraints. It is shown that the solution of this penalty method converges strongly to the solution of the constrained problem. We consider the abstract numerical analysis of these problems. Numerical experiments are given, demonstrating the convergence theory presented. After this, we consider the differentiability of the energy of the optimal membrane with point constraints with respect to a tangential movement of the points. We demonstrate that the energy is differentiable and give a convenient characterisation of the derivative which is efficient to evaluate. This numerically accessible derivative is employed in some numerical experiments. We conclude by discussing some directions to extend the theory presented, or ideas which are highly related to the studied theory. In particular, we discuss the extension to consider small deformations of a near-tube membrane