Level set methods for higher order evolution laws

A numerical treatment of non-linear higher-order geometric evolution equations with the level set and the finite element method is presented. The isotropic, weak anisotropic and strong anisotropic situation is discussed. Most of the equations considered in this work arise from the field of thin film growth. A short introduction to the subject is given. Four different models are discussed: mean curvature flow, surface diffusion, a kinetic model, which combines the effects of mean curvature flow and surface diffusion and includes a further kinetic component, and an adatom model, which incorporates in addition free adatoms. As an introduction to the numerical schemes, first the isotropic and weak anisotropic situation is considered. Then strong anisotropies (non-convex anisotropies) are used to simulate the phenomena of faceting and coarsening. The experimentally observed effect of corner and edge roundings is reached in the simulation through the regularization of the strong anisotropy with a higher-order curvature term. The curvature regularization leads to an increase by two in the order of the equations, which results in highly non-linear equations of up to 6th order. For the numerical solution, the equations are transformed into systems of second order equations, which are solved with a Schur complement approach. The adatom model constitutes a diffusion equation on a moving surface. An operator splitting approach is used for the numerical solution. In difference to other works, which restrict to the isotropic situation, also the anisotropic situation is discussed and solved numerically. Furthermore, a treatment of geometric evolution equations on implicitly given curved surfaces with the level set method is given. In particular, the numerical solution of surface diffusion on curved surfaces is presented. The equations are discretized in space by standard linear finite elements. For the time discretization a semi-implicit discretization scheme is employed. The derivation of the numerical schemes is presented in detail, and numerous computational results are given for the 2D and 3D situation. To keep computational costs low, the finite element grid is adaptively refined near the moving curves and surfaces resp. A redistancing algorithm based on a local Hopf-Lax formula is used. The algorithm has been extended by the authors to the 3D case. A detailed description of the algorithm in 3D is presented in this work.

Level thesis:degree_level

thesis.doctoral

Grantor dc:publisher

Forschungszentrum caesar

Year

2008

Author dc:creator

• Stöcker, Christina

Contributors dc:contributor

• Voigt, Axel
• Burger, Martin
• Smereka, Peter

• applied mathematics
• geometric evolution equations
• level set method
• finite element method
• free boundary problem
• crystal growth
• thin film growth
• strong anisotropy
• non-convex anisotropy
• curvature regularization
• faceting
• coarsening
• PDE on moving s
• Numerische Mathematik
• geometrische Evolutionsgleichungen
• Levelset-Methode
• Finite-Elemente-Methode
• freies Randwertproblem
• Kristallwachstum
• Dünnschichtwachstum
• starke Anisotropie
• nicht-konvexe Anisotropie
• Krümmungsregularisierung
• Facettierung
• Ver