Shawn W. Walker

Ricardo H. Nochetto, Abner J. Salgado, Shawn W. Walker

We introduce a diffuse interface model for the phenomenon of electrowetting on dielectric and present an analysis of the arising system of equations. Moreover, we study discretization techniques for the problem. The model takes into account different material parameters on each phase and incorporates the most important physical processes, such as incompressibility, electrostatics and dynamic contact lines; necessary to properly reflect the relevant phenomena. The arising nonlinear system couples the variable density incompressible Navier-Stokes equations for velocity and pressure with a Cahn-Hilliard type equation for the phase variable and chemical potential, a convection diffusion equation for the electric charges and a Poisson equation for the electric potential. Numerical experiments are presented, which illustrate the wide range of effects the model is able to capture, such as splitting and coalescence of droplets.

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.

Ricardo H. Nochetto, Shawn W. Walker, Wujun Zhang

We consider the simplest one-constant model, put forward by J. Ericksen, for nematic liquid crystals with variable degree of orientation. The equilibrium state is described by a director field $\mathbf{n}$ and its degree of orientation $s$, where the pair $(s, \mathbf{n})$ minimizes a sum of Frank-like energies and a double well potential. In particular, the Euler-Lagrange equations for the minimizer contain a degenerate elliptic equation for $\mathbf{n}$, which allows for line and plane defects to have finite energy. We present a structure preserving discretization of the liquid crystal energy with piecewise linear finite elements that can handle the degenerate elliptic part without regularization, and show that it is consistent and stable. We prove $Γ$-convergence of discrete global minimizers to continuous ones as the mesh size goes to zero. We develop a quasi-gradient flow scheme for computing discrete equilibrium solutions and prove it has a strictly monotone energy decreasing property. We present simulations in two and three dimensions to illustrate the method's ability to handle non-trivial defects. A music video summary of the paper is available on YouTube: "Mathematical Modeling and Simulation of Nematic Liquid Crystals (A Montage)," http://www.youtube.com/watch?v=pWWw7\_6cQ-U

Ricardo H. Nochetto, Shawn W. Walker, Wujun Zhang

We present a robust discretization of the Ericksen model of liquid crystals with variable degree of orientation coupled with colloidal effects and electric fields. The total energy consists of the Ericksen energy, a weak anchoring (or penalized Dirichlet) energy to model colloids, and an electrical energy for a given electric field. We describe our special discretization of the total energy along with a method to compute minimizers via a discrete quasi-gradient flow algorithm which has a strictly monotone energy decreasing property. Numerical experiments are given in two and three dimensions to illustrate that the method is able to capture non-trivial defect patterns, such as the Saturn ring defect. We conclude with a rigorous proof of the Gamma-convergence of our discrete energy to the continuous energy.