Convex minimization and phase field models

Abstract

Phase field models are a widely used approach to describe physical processes that are characterized by thin interfacial regions between large almost homogeneous domains. Important application areas of phase field models are transition processes of the state of matter and the separation of alloys. A fundamental property of these models is, that the transition and separation of phases is driven by a double-well potential with distinct minima for the different phases. Already the pioneering work of Cahn and Hilliard used a temperature dependent logarithmic potential that is differentiable with singular derivatives. If the temperature tends to zero it degenerates to the non-differentiable obstacle potential. The goal of this thesis is to develop methods for the efficient numerical solution of such equations that are also robust for nonsmooth potentials and anisotropic surface energies. These methods are derived for the Cahn-Hilliard equation that are prototypic for a multitude of such models. The main result of the thesis is the development of a fast iterative solver for nonlinear saddle point problems like the ones that arise from finite element discretization of Cahn-Hilliard equations. The solver relies on a reformulation of the problem as dual minimization problem whose energy functional is differentiable. The gradient of this functional turns out to be the nonlinear Schur complement of the saddle point problem. Generalized linearizations for the Schur complement are derived and used for a nonsmooth Newton method. Global convergence for this 'Schur Nonsmooth Newton' method and inexact versions is proved using the fact the equivalence to a descent method for the dual minimization problem. Each step of this method requires the solution of a nonlinear convex minimization problem. To tackle this problem the 'Truncated Nonsmooth Newton Multigrid' (TNNMG) method is developed. In contrast to other nonlinear multigrid methods the TNNMG method is significantly easier to implement and can also be applied to anisotropic problems while its convergence speed is in general comparable or sometimes even faster. Numerical examples show that the derived methods exhibit mesh independent convergence. Furthermore they turn out to be robust with respect to the temperature including the limiting case zero. The reason for this robustness is, that all methods do not rely on smoothness but on the inherent convex structure of the problems.

Phasenfeldmodelle sind ein weit verbreiteter Ansatz zur Beschreibung von Prozessen, die sich wesentlich durch dünne Interfaceregionen zwischen weitestgehend homogenen Bereichen auszeichnen. Ein wichtiges Anwendungsfeld von Phasenfeldmodellen ist die Modellierung von physikalischen Phasenübergangs- und Phasenseparationvorgängen. Eine wesentliche Eigenschaften ist dabei, dass die Trennung der Phasen durch ein Doppelmuldenpotential getrieben ist, welches voneinander getrennte Minima für jede Phase besitzt. Bereits 1958 haben Cahn und Hilliard ein logarithmisches Potential vorgeschlagen, das zwar differenzierbar ist, aber singuläre Ableitungen besitzt. Geht die Temperatur gegen Null, so degeneriert das temperaturabhängige logarithmische Potential gegen das nichtdifferenzierbare Hindernispotential. Ziel der vorliegenden Arbeit ist die Entwicklung von Methoden zur effizienten numerischen Lösung solcher Gleichungen, die auch im Fall nichtglatter Potentiale und anisotroper Oberflächenenergien robust sind. Diese Methoden werden für die Cahn-Hilliard-Gleichung entwickelt, die prototypisch für eine Vielzahl solcher Modelle ist. Das Hauptresultat der Arbeit ist die Entwicklung eines schnellen iterativen Verfahrens zur Lösung nichtlinearer Sattelpunktprobleme, wie sie bei der Diskretisierung anisotroper Cahn-Hilliard-Gleichungen mit Finite-Elemente-Methoden entstehen. Die Grundlage dieses Verfahrens ist eine Umformulierung des Sattelpunktproblems als äquivalentes duales Minimierungsproblem. Das Energiefunktional dieses Minimierungsproblems ist differenzierbar und seine Ableitung ist das nichtlineare Schur-Komplement des Sattelpunktproblems. Für dieses Schur- Komplement wird eine verallgemeinerte Linearisierung hergeleitet, die im Rahmen eines nichtglatten Newton-Verfahrens Verwendung findet. Für dieses so genannte "Schur Nonsmooth Newton" Verfahren wird die globale Konvergenz mittels der Äquivalenz zu Abstiegsverfahren für das duale Minimierungsproblem gezeigt. Ferner wird bewiesen, dass die globale Konvergenz auch bei inexaktem Lösen der linearen Newton Probleme erhalten bleibt. In jedem Schritt dieses Verfahrens ist ein nichtlineares konvexes Minimierungsproblem zu lösen. Für die effiziente Behandlung dieser Teilprobleme wird das so genannte "Truncated Nonsmooth Newton Multigrid" (TNNMG) Verfahren entwickelt. Dieses nichtlineare Mehrgitterverfahren zeichnet sich im Gegensatz zu verwandten Verfahren durch seine Einfachheit und die Anwendbarkeit auf anisotrope Probleme aus. Hinsichtlich der Konvergenzgeschwindigkeit ist es bereits bekannten Verfahren ebenbürtig. Numerische Beispiele zeigen, dass die entwickelten Verfahren gitterunabhängig konvergieren. Auch erweisen sich die Verfahren als robust für verschiedene Temperaturen, einschließlich dem Grenzfall Null. Der Grund für diese Robustheit ist, dass die Verfahren nicht auf Differenzierbarkeit sondern der konvexen Struktur der Probleme basieren.