A Moving Mesh Method for Modelling Defects in Nematic Liquid Crystals

The properties of liquid crystals can be modelled using an order parameter which describes the variability of the local orientation of rod-like molecules. Defects in the director field can arise due to external factors such as applied electric or magnetic fields, or the constraining geometry of the cell containing the liquid crystal material. Understanding the formation and dynamics of defects is important in the design and control of liquid crystal devices, and poses significant challenges for numerical modelling. In this paper we consider the numerical solution of a $\bf{Q}$-tensor model of a nematic liquid crystal, where defects arise through rapid changes in the $\bf{Q}$-tensor over a very small physical region in relation to the dimensions of the liquid crystal device. The efficient solution of the resulting six coupled partial differential equations is achieved using a finite element based adaptive moving mesh approach, where an unstructured triangular mesh is adapted towards high activity regions, including those around defects. Spatial convergence studies are presented using a stationary defect as a model test case, and the adaptive method is shown to be optimally convergent using quadratic triangular finite elements. The full effectiveness of the method is then demonstrated using a challenging two-dimensional dynamic Pi-cell problem involving the creation, movement, and annihilation of defects.

[1]  A. Fernández-Nieves,et al.  Drops and shells of liquid crystal , 2011 .

[2]  M. Schadt,et al.  Voltage-Dependent Optical Activity of a Twisted Nematic Liquid Crystal , 1971 .

[3]  A. Sonnet,et al.  Dissipative Ordered Fluids: Theories for Liquid Crystals , 2012 .

[4]  Weizhang Huang,et al.  On the mesh nonsingularity of the moving mesh PDE method , 2018, Math. Comput..

[5]  W. Marsden I and J , 2012 .

[6]  A. Sonnet,et al.  Elastic forces on nematic point defects , 2002 .

[7]  Alternative dynamics for the symmetric splay to bend transition in a nematic liquid crystal layer , 2007 .

[8]  S. Hess,et al.  Alignment tensor versus director: Description of defects in nematic liquid crystals. , 1995, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[9]  Riccardo Barberi,et al.  Electric field-induced fast nematic order dynamics , 2011 .

[10]  T. Sluckin,et al.  Hedgehog structure in nematic and magnetic systems , 1988 .

[11]  Robert D. Russell,et al.  Adaptive Moving Mesh Methods , 2010 .

[12]  P. Bos,et al.  Fast Q-tensor method for modeling liquid crystal director configurations with defects , 2002 .

[13]  S. Griffis EDITOR , 1997, Journal of Navigation.

[14]  T. Sluckin,et al.  Defect core structure in nematic liquid crystals. , 1987, Physical review letters.

[15]  Alison Ramage,et al.  Adaptive solution of a one‐dimensional order reconstruction problem in Q‐tensor theory of liquid crystals , 2007 .

[16]  Robert D. Russell,et al.  Moving Mesh Strategy Based on a Gradient Flow Equation for Two-Dimensional Problems , 1998, SIAM J. Sci. Comput..

[17]  Giuseppe Lombardo,et al.  Moving mesh partial differential equations to describe nematic order dynamics , 2010, Comput. Math. Appl..

[18]  Timothy A. Davis,et al.  Finite Element Analysis of the Landau--de Gennes Minimization Problem for Liquid Crystals , 1998 .

[19]  R. Varga,et al.  Numerical Minimization of the Landau-de Gennes Free Energy: Defects in Cylindrical Capillaries , 1991 .

[20]  P. G. de Gennes,et al.  Short Range Order Effects in the Isotropic Phase of Nematics and Cholesterics , 1971 .

[21]  J. Mackenzie,et al.  On a uniformly accurate finite difference approximation of a singulary peturbed reaction-diffusion problem using grid equidistribution , 2001 .

[22]  H. Trebin,et al.  Fine structure of point defects and soliton decay in nematic liquid crystals , 1989 .

[23]  Jun-ichi Fukuda,et al.  Director configuration and dynamics of a nematic liquid crystal around a two-dimensional spherical particle: Numerical analysis using adaptive grids , 2001 .

[24]  A. Ramage,et al.  Computational solution of two-dimensional unsteady PDEs using moving mesh methods , 2002 .

[25]  John A. Mackenzie,et al.  Robust adaptive computation of a one-dimensional Q-tensor model of nematic liquid crystals , 2012, Comput. Math. Appl..

[26]  F. Ciuchi,et al.  Electric field induced order reconstruction in a nematic cell , 2004, The European physical journal. E, Soft matter.

[27]  E. Virga,et al.  Biaxial torus around nematic point defects. , 1999, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[28]  Alison Ramage,et al.  Adaptive Grid Methods for Q-Tensor Theory of Liquid Crystals: A One-Dimensional Feasibility Study , 2008 .

[29]  Robert D. Russell,et al.  Adaptivity with moving grids , 2009, Acta Numerica.

[30]  Pingwen Zhang,et al.  Moving mesh methods in multiple dimensions based on harmonic maps , 2001 .

[31]  Nigel J. Mottram,et al.  Introduction to Q-tensor theory , 2014, 1409.3542.

[32]  P. D. Gennes,et al.  An analogy between superconductors and smectics A , 1972 .

[33]  M. Yoneya,et al.  Physics of Liquid Crystals , 2014 .

[34]  S. Hess,et al.  Relaxation of the Alignment Tensor in the Isotropic and Nematic Phases of Liquid Crystals in the Presence of an External Field , 1991 .

[35]  G. Boyd,et al.  Liquid‐crystal orientational bistability and nematic storage effects , 1980 .

[36]  Henk A. van der Vorst,et al.  Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems , 1992, SIAM J. Sci. Comput..

[37]  A. Ramage,et al.  On the numerical solution of one-dimensional PDEs using adaptive methods based on equidistribution , 2001 .

[38]  P. Bos,et al.  The pi-Cell: A Fast Liquid-Crystal Optical-Switching Device , 1984 .

[39]  Robert D. Russell,et al.  A comparative numerical study of meshing functionals for variational mesh adaptation , 2015 .

[40]  John A. Mackenzie,et al.  Efficient Moving Mesh Methods for Q-Tensor Models of Nematic Liquid Crystals , 2015, SIAM J. Sci. Comput..

[41]  A. Dvinsky Adaptive grid generation from harmonic maps on Reimannian manifolds , 1991 .

[42]  E. Virga,et al.  Perspectives in active liquid crystals , 2014, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[43]  E. Willman,et al.  Finite-element modeling of liquid-crystal hydrodynamics with a variable degree of order , 2006, IEEE Transactions on Electron Devices.

[44]  J. Meijerink,et al.  An iterative solution method for linear systems of which the coefficient matrix is a symmetric -matrix , 1977 .

[45]  A. M. Winslow Numerical Solution of the Quasilinear Poisson Equation in a Nonuniform Triangle Mesh , 1997 .

[46]  Jun-ichi Fukuda,et al.  Defect structure of a nematic liquid crystal around a spherical particle: adaptive mesh refinement approach. , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

[47]  George Beckett,et al.  Convergence analysis of finite difference approximations on equidistributed grids to a singularly perturbed boundary value problem , 2000 .