Numerical Analysis of Convex Splitting Schemes for Cahn-Hilliard and Coupled Cahn-Hilliard-Fluid-Flow Equations

This dissertation investigates numerical schemes for the Cahn-Hilliard equation and the Cahn-Hilliard equation coupled with a Darcy-Stokes flow. Considered independently, the Cahn-Hilliard equation is a model for spinodal decomposition and domain coarsening. When coupled with a Darcy-Stokes flow, the resulting system describes the flow of a very viscous block copolymer fluid. Challenges in creating numerical schemes for these equations arise due to the nonlinear nature and high derivative order of the Cahn-Hilliard equation. Further challenges arise during the coupling process as the coupling terms tend to be nonlinear as well. The numerical schemes presented herein preserve the energy dissipative structure of the Cahn- Hilliard equation while maintaining unique solvability and optimal error bounds. Specifically, we devise and analyze two mixed finite element schemes: a first order in time numerical scheme for a modified Cahn-Hilliard equation coupled with a non- steady Darcy-Stokes flow and a second order in time numerical scheme for the Cahn- Hilliard equation in two and three dimensions. The time discretizations are based on a convex splitting of the energy of the systems. We prove that our schemes are unconditionally energy stable with respect to a spatially discrete analogue of the continuous free energies and unconditionally uniquely solvable. For each system, we prove that the discrete phase variable is essentially bounded in both time and space with respect to the Lebesque integral and the discrete chemical potential is Lesbegue square integrable in space and essentially bounded in time. We show these bounds are completely independent of the time and space step sizes in two and three dimensions. We subsequently prove that these variables converge with optimal rates in the appropriate energy norms. The analyses included in this dissertation will provide a bridge to the development of stable, efficient, and optimally convergent numerical schemes for more robust and descriptive coupled Cahn-Hilliard-Fluid-Flow systems