Solid-state dewetting of axisymmetric thin film on axisymmetric curved-surface substrates: Modeling and simulation

Abstract

In this work, we consider the solid-state dewetting of an axisymmetric thin film on a curved-surface substrate, with the assumption that the substrate morphology is also axisymmetric. Under the assumptions of axisymmetry, the surface evolution problem on a curved-surface substrate can be reduced to a curve evolution problem on a static curved substrate. Based on the thermodynamic variation of the anisotropic surface energy, we thoroughly derive a sharp-interface model that is governed by anisotropic surface diffusion, along with appropriate boundary conditions. The continuum system satisfies the laws of energy decay and volume conservation, which motivates the design of a structure-preserving numerical algorithm for simulating the mathematical model. We introduce an arclength parameterization of the generated curve on the axisymmetric curved substrate surface, which plays a crucial role in the subsequent construction of the structure-preserving approximation. By introducing a symmetrized surface energy matrix, we derive a novel symmetrized variational formulation. Then, by carefully discretizing the boundary terms of the variational formulation, we establish an unconditionally energy-stable parametric finite element approximation of the axisymmetric system. By applying an ingenious correction method, we further develop another structure-preserving method that can preserve both the energy stability and volume conservation properties. Finally, we present extensive numerical examples to demonstrate the convergence and structure-preserving properties of our proposed numerical scheme. Additionally, several interesting phenomena are explored, including the migration of ‘small’ particles on a curved-surface substrate generated by curves with positive or negative curvature, pinch-off events, and edge retraction.