PIMS- SFU CSC Seminar: Jessica Bosch

The Cahn-Hilliard equation models the motion of interfaces between several phases. The underlying energy functional includes a potential for which different types were proposed in the literature. We consider smooth and nonsmooth potentials with a focus on the latter. In the nonsmooth case, we apply a function space-based algorithm, which combines a Moreau-Yosida regularization technique with a semismooth Newton method. We apply classical finite element methods to discretize the problems in space. At the heart of our method lies the solution of large and sparse fully discrete systems of linear equations. Block preconditioners using effective Schur complement approximations are presented. For the smooth systems, we derive optimal preconditioners, which are proven to be robust with respect to crucial model parameters. Further, we prove that the use of the same preconditioners give poor approximations for the nonsmooth formulations. The preconditioners we present for the nonsmooth problems incorporate the regularization terms.

This approach provides fairly robust results in practice. Our strategy applies to different Cahn-Hilliard problems including phase separation and coarsening processes and image inpainting. Additionally, our preconditioner can be embedded into a fully coupled Cahn-Hilliard Navier-Stokes-system for the simulation of two-phase flows.

This is joint work with Martin Stoll.