An Energy-Stable Parametric Finite Element Method for Simulating Solid-state Dewetting Problems in Three Dimensions

We propose an accurate and energy-stable parametric finite element method for solving the sharp-interface continuum model of solid-state dewetting in three-dimensional space. The model describes the motion of the film\slash vapor interface with contact line migration and is governed by the surface diffusion equation with proper boundary conditions at the contact line. We present a new weak formulation for the problem, in which the interface and its contact line are evolved simultaneously. By using piecewise linear elements in space and backward Euler in time, we then discretize the weak formulation to obtain a fully discretized parametric finite element approximation. The resulting numerical method is shown to be well-posed and unconditionally energy-stable. Furthermore, the numerical method is extended for solving the sharp interface model of solid-state dewetting with anisotropic surface energies in the Riemmanian metric form. Numerical results are reported to show the convergence and efficiency of the proposed numerical method as well as the anisotropic effects on the morphological evolution of thin films in solid-state dewetting.Comment: 25 pages, 11 figure

A perimeter-decreasing and area-conserving algorithm for surface diffusion flow of curves

A fully discrete finite element method, based on a new weak formulation and a new time-stepping scheme, is proposed for the surface diffusion flow of closed curves in the two-dimensional plane. It is proved that the proposed method can preserve two geometric structures simultaneously at the discrete level, i.e., the perimeter of the curve decreases in time while the area enclosed by the curve is conserved. Numerical examples are provided to demonstrate the convergence of the proposed method and the effectiveness of the method in preserving the two geometric structures