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

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 boundary control of a simplified Ericksen--Leslie system for nematic liquid crystal flows in 2D2D

In this paper, we investigate an optimal boundary control problem for a two dimensional simplified Ericksen--Leslie system modelling the incompressible nematic liquid crystal flows. The hydrodynamic system consists of the Navier--Stokes equations for the fluid velocity coupled with a convective Ginzburg--Landau type equation for the averaged molecular orientation. The fluid velocity is assumed to satisfy a no-slip boundary condition, while the molecular orientation is subject to a time-dependent Dirichlet boundary condition that corresponds to the strong anchoring condition for liquid crystals. We first establish the existence of optimal boundary controls. Then we show that the control-to-state operator is Fr\'echet differentiable between appropriate Banach spaces and derive first-order necessary optimality conditions in terms of a variational inequality involving the adjoint state variables

Relating phase field and sharp interface approaches to structural topology optimization

A phase field approach for structural topology optimization which allows for topology changes and multiple materials is analyzed. First order optimality conditions are rigorously derived and it is shown via formally matched asymptotic expansions that these conditions converge to classical first order conditions obtained in the context of shape calculus. We also discuss how to deal with triple junctions where e.g. two materials and the void meet. Finally, we present several numerical results for mean compliance problems and a cost involving the least square error to a target displacement

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

Surface, Bulk, and Geometric Partial Differential Equations: Interfacial, stochastic, non-local and discrete structures

Partial differential equations in complex domains with free\linebreak boundaries and interfaces continue to be flourishing research areas at the interfaces between PDE theory, differential geometry, numerical analysis and applications. Main themes of the workshop have been PDEs on evolving domains, phase field approaches, interactions of bulk and surface PDEs, curvature driven evolution equations. Applications particular from biology, such as cell and cancer modelling and fluid as well solid mechanics have been subjects of the conference

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