Abstract. The Cahn-Hilliard variational inequality is a non-standard parabolic variational inequality of fourth order for which straightforward numerical approaches cannot be applied. We propose a primal-dual active set method which can be interpreted as a semi-smooth Newton method as solution technique for the discretized Cahn-Hilliard variational inequality. A (semi-)implicit Euler discretization is used in time and a piecewise linear finite element discretization of splitting type is used in space leading to a discrete variational inequality of saddle point type in each time step. In each iteration of the primal-dual active set method a linearized system resulting from the discretization of two coupled elliptic equations which are defined on different sets has to be solved. We show local convergence of the primal-dual active set method and demonstrate its efficiency with several numerical simulations.Mathematics Subject Classification. 35K55, 35K85, 90C33, 49N90, 80A22, 82C26, 65M60.

Abstract. Parabolic variational inequalities of Allen-Cahn and CahnHilliard 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.Mathematics Subject Classification (2010). 35K55, 35S85, 65K10, 90C33, 90C53, 49N90, 65M60.