Amanda E. Diegel

Amanda E. Diegel, Cheng Wang, Xiaoming Wang et al.

In this paper, we present a novel second order in time mixed finite element scheme for the Cahn-Hilliard-Navier-Stokes equations with matched densities. The scheme combines a standard second order Crank-Nicholson method for the Navier-Stokes equations and a modification to the Crank-Nicholson method for the Cahn-Hilliard equation. In particular, a second order Adams-Bashforth extrapolation and a trapezoidal rule are included to help preserve the energy stability natural to the Cahn-Hilliard equation. We show that our scheme is unconditionally energy stable with respect to a modification of the continuous free energy of the PDE system. Specifically, the discrete phase variable is shown to be bounded in $\ell^\infty \left(0,T;L^\infty\right)$ and the discrete chemical potential bounded in $\ell^\infty \left(0,T;L^2\right)$, for any time and space step sizes, in two and three dimensions, and for any finite final time $T$. We subsequently prove that these variables along with the fluid velocity converge with optimal rates in the appropriate energy norms in both two and three dimensions.

Susanne C. Brenner, Amanda E. Diegel, Li-Yeng Sung

We develop a robust solver for a second order mixed finite element splitting scheme for the Cahn-Hilliard equation. This work is an extension of our previous work in which we developed a robust solver for a first order mixed finite element splitting scheme for the Cahn-Hilliard equaion. The key ingredient of the solver is a preconditioned minimal residual algorithm (with a multigrid preconditioner) whose performance is independent of the spacial mesh size and the time step size for a given interfacial width parameter. The dependence on the interfacial width parameter is also mild.

Amanda E. Diegel, Shawn W. Walker

We develop a novel finite element method for a phase field model of nematic liquid crystal droplets. The continuous model considers a free energy comprised of three components: the Ericksen's energy for liquid crystals, the Cahn-Hilliard energy representing the interfacial energy of the droplet, and a weak anchoring energy representing the interaction of the liquid crystal molecules with the surface tension on the interface (i.e. anisotropic surface tension). Applications of the model are for finding minimizers of the free energy and exploring gradient flow dynamics. We present a finite element method that utilizes a special discretization of the liquid crystal elastic energy, as well as mass-lumping to discretize the coupling terms for the anisotropic surface tension part. Next, we present a discrete gradient flow method and show that it is monotone energy decreasing. Furthermore, we show that global discrete energy minimizers $Γ$-converge to global minimizers of the continuous energy. We conclude with numerical experiments illustrating different gradient flow dynamics, including droplet coalescence and break-up.

Susanne C. Brenner, Amanda E. Diegel, Li-Yeng Sung

We develop a robust solver for a mixed finite element convex splitting scheme for the Cahn-Hilliard equation. The key ingredient of the solver is a preconditioned minimal residual algorithm (with a multigrid preconditioner) whose performance is independent of the spacial mesh size and the time step size for a given interfacial width parameter. The dependence on the interfacial width parameter is also mild.

Amanda E. Diegel, Cheng Wang, Steven M. Wise

In this paper we devise and analyze an unconditionally stable, second-order-in-time numerical scheme for the Cahn-Hilliard equation in two and three space dimensions. We prove that our two-step scheme is unconditionally energy stable and unconditionally uniquely solvable. Furthermore, we show that the discrete phase variable is bounded in $L^\infty (0,T;L^\infty)$ and the discrete chemical potential is bounded in $L^\infty (0,T;L^2)$, for any time and space step sizes, in two and three dimensions, and for any finite final time $T$. We subsequently prove that these variables converge with optimal rates in the appropriate energy norms in both two and three dimensions. We include in this work a detailed analysis of the initialization of the two-step scheme.