Allen-Cahn and Cahn-Hilliard variational inequalities solved with Optimization Techniques

Parabolic variational inequalities of Allen-Cahn and Cahn- Hilliard type are solved using methods involving constrained optimization. Time discrete variants are formulated with the help of Lagrange multipliers for local and non-local equality and inequality constraints. Fully discrete problems resulting from finite element discretizations in space are solved with the help of a primal-dual active set approach. We show several numerical computations also involving systems of parabolic variational inequalities

Fast solution of Cahn-Hilliard variational inequalities using implicit time discretization and finite elements

We consider the e�cient solution of the Cahn-Hilliard variational inequality using an implicit time discretization, which is formulated as an optimal control problem with pointwise constraints on the control. By applying a semi-smooth Newton method combined with a Moreau-Yosida regularization technique for handling the control constraints we show superlinear convergence in function space. At the heart of this method lies the solution of large and sparse linear systems for which we propose the use of preconditioned Krylov subspace solvers using an e�ective Schur complement approximation. Numerical results illustrate the competitiveness of this approach

Phase-field approaches to structural topology optimization

The mean compliance minimization in structural topology optimization is solved with the help of a phase field approach. Two steepest descent approaches based on L2- and H-1 gradient flow dynamics are discussed. The resulting flows are given by Allen-Cahn and Cahn-Hilliard type dynamics coupled to a linear elasticity system. We finally compare numerical results obtained from the two different approaches

Optimal control of Allen-Cahn systems

Optimization problems governed by Allen-Cahn systems including elastic effects are formulated and first-order necessary optimality conditions are presented. Smooth as well as obstacle potentials are considered, where the latter leads to an MPEC. Numerically, for smooth potential the problem is solved efficiently by the Trust-Region-Newton-Steihaug-cg method. In case of an obstacle potential first numerical results are presented

Optimal distributed control of a nonlocal convective Cahn-Hilliard equation by the velocity in 3D

In this paper we study a distributed optimal control problem for a nonlocal convective Cahn--Hilliard equation with degenerate mobility and singular potential in three dimensions of space. While the cost functional is of standard tracking type, the control problem under investigation cannot easily be treated via standard techniques for two reasons: the state system is a highly nonlinear system of PDEs containing singular and degenerating terms, and the control variable, which is given by the velocity of the motion occurring in the convective term, is nonlinearly coupled to the state variable. The latter fact makes it necessary to state rather special regularity assumptions for the admissible controls, which, while looking a bit nonstandard, are however quite natural in the corresponding analytical framework. In fact, they are indispensable prerequisites to guarantee the well-posedness of the associated state system. In this contribution, we employ recently proved existence, uniqueness and regularity results for the solution to the associated state system in order to establish the existence of optimal controls and appropriate first-order necessary optimality conditions for the optimal control problem

Preconditioning for Allen-Cahn variational inequalities with non-local constraints

The solution of Allen-Cahn variational inequalities with mass constraints is of interest in many applications. This problem can be solved both in its scalar and vector-valued form as a PDE-constrained optimization problem by means of a primal-dual active set method. At the heart of this method lies the solution of linear systems in saddle point form. In this paper we propose the use of Krylov-subspace solvers and suitable preconditioners for the saddle point systems. Numerical results illustrate the competitiveness of this approach

A boundary control problem for the viscous Cahn-Hilliard equation with dynamic boundary conditions

A boundary control problem for the viscous Cahn-Hilliard equations with possibly singular potentials and dynamic boundary conditions is studied and first order necessary conditions for optimality are proved. Key words: Cahn-Hilliard equation, dynamic boundary conditions, phase separation, singular potentials, optimal control, optimality conditions, adjoint state syste

Primal-dual active set methods for Allen-Cahn variational inequalities

This thesis aims to introduce and analyse a primal-dual active set strategy for solving Allen-Cahn variational inequalities. We consider the standard Allen-Cahn equation with non-local constraints and a vector-valued Allen-Cahn equation with and without non-local constraints. Existence and uniqueness results are derived in a formulation involving Lagrange multipliers for local and non-local constraints. Local Convergence is shown by interpreting the primal-dual active set approach as a semi-smooth Newton method. Properties of the method are discussed and several numerical simulations in two and three space dimensions demonstrate its efficiency. In the second part of the thesis various applications of the Allen-Cahn equation are discussed. The non-local Allen-Cahn equation can be coupled with an elasticity equation to solve problems in structural topology optimisation. The model can be extended to handle multiple structures by using the vector-valued Allen-Cahn variational inequality with non-local constraints. Since many applications of the Allen-Cahn equation involve evolution of interfaces in materials an important extension of the standard Allen-Cahn model is to allow materials to exhibit anisotropic behaviour. We introduce an anisotropic version of the Allen-Cahn variational inequality and we show that it is possible to apply the primal-dual active set strategy efficiently to this model. Finally, the Allen-Cahn model is applied to problems in image processing, such as segmentation, denoising and inpainting. The primal-dual active set method proves exible and reliable for all the applications considered in this thesis

Emerging Developments in Interfaces and Free Boundaries

The field of the mathematical and numerical analysis of systems of nonlinear partial differential equations involving interfaces and free boundaries is a well established and flourishing area of research. This workshop focused on recent developments and emerging new themes. By bringing together experts in these fields we achieved progress in open questions and developed novel research directions in mathematics related to interfaces and free boundaries. This interdisciplinary workshop brought together researchers from distinct mathematical fields such as analysis, computation, optimisation and modelling to discuss emerging challenges

New Directions in Simulation, Control and Analysis for Interfaces and Free Boundaries

The field of mathematical and numerical analysis of systems of nonlinear partial differential equations involving interfaces and free boundaries is a flourishing area of research. Many such systems arise from mathematical models in material science, fluid dynamics and biology, for example phase separation in alloys, epitaxial growth, dynamics of multiphase fluids, evolution of cell membranes and in industrial processes such as crystal growth. The governing equations for the dynamics of the interfaces in many of these applications involve surface tension expressed in terms of the mean curvature and a driving force. Here the forcing terms depend on variables that are solutions of additional partial differential equations which hold either on the interface itself or in the surrounding bulk regions. Often in applications of these mathematical models, suitable performance indices and appropriate control actions have to be specified. Mathematically this leads to optimization problems with partial differential equation constraints including free boundaries. Because of the maturity of the field of computational free boundary problems it is now timely to consider such control problems