A structure-preserving finite element approximation of surface diffusion for curve networks and surface clusters

Abstract

We consider the evolution of curve networks in two dimensions (2d) and surface clusters in three dimensions (3d). The motion of the interfaces is described by surface diffusion, with boundary conditions at the triple junction points/lines, where three interfaces meet, and at the boundary points/lines, where an interface meets a fixed planar boundary. We propose a parametric finite element method based on a suitable variational formulation. The constructed method is semi-implicit and can be shown to satisfy the volume conservation of each enclosed bubble and the unconditional energy-stability, thus preserving the two fundamental geometric structures of the flow. Besides, the method has very good properties with respect to the distribution of mesh points, thus no mesh smoothing or regularization technique is required. A generalization of the introduced scheme to the case of anisotropic surface energies and non-neutral external boundaries is also considered. Numerical results are presented for the evolution of two-dimensional curve networks and three-dimensional surface clusters in the cases of both isotropic and anisotropic surface energies.

arXiv comment29 figures

References in corpus (11)

• Parametric finite element approximations of curvature driven interface evolutions
• A parametric finite element method for solid-state dewetting problems with anisotropic surface energies
• A structure-preserving parametric finite element method for surface diffusion
• A parametric finite element method for solid-state dewetting problems in three dimensions
• A perimeter-decreasing and area-conserving algorithm for surface diffusion flow of curves
• An energy-stable parametric finite element method for anisotropic surface diffusion
• Regular three-dimensional bubble clusters: shape, packing and growth-rate
• A structure preserving front tracking finite element method for the Mullins--Sekerka problem
• Are large two-dimensional clusters of perimeter-minimizing bubbles of equal-area hexagonal or circular?
• An Energy-Stable Parametric Finite Element Method for Simulating Solid-state Dewetting Problems in Three Dimensions
• Standard Planar Double Bubbles are Stable under Surface Diffusion Flow

Cited by in corpus (4)

• Structure-preserving discretizations of two-phase Navier-Stokes flow using fitted and unfitted approaches
• A diffuse-interface approach for solid-state dewetting with anisotropic surface energies
• A parametric finite element method for a degenerate multi-phase Stefan problem with triple junctions
• A structure-preserving finite element method for the multi-phase Mullins-Sekerka problem with triple junctions