An adaptive homotopy multi-grid method for molecule orientations of high dimensional liquid crystals

The liquid crystal molecule orientation is arranged by minimizing the so-called Oseen-Frank energy functional. For a better understanding of these complicated orientation singularities, simplified models resulting from specific choices of elastic constants are always of interest. In this paper a pseudo Newton method together with a multi-grid linear system solver or preconditioner is used to compute the orientation of liquid crystal molecules based on a simplified Oseen-Frank energy functional. The penalty method is used to deal with the unit-length constraint of liquid crystal molecules. The Newton and multi-grid methods do not converge when some parameters are small. A homotopy algorithm combined with mesh refinement strategies in order to deal with small parameter cases is studied and is found to be very robust in computing the solution of the model. The method is implemented to compute the orientation of liquid crystal molecules in domains of typical shapes and with various rotational boundary conditions in 2D and 3D. Interesting singularity patterns are observed.

[1]  P. Gennes,et al.  The physics of liquid crystals , 1974 .

[2]  T. Richter,et al.  SOLUTIONS OF 3D NAVIER-STOKES BENCHMARK PROBLEMS WITH ADAPTIVE FINITE ELEMENTS , 2006 .

[3]  R. Glowinski,et al.  Augmented Lagrangian and Operator-Splitting Methods in Nonlinear Mechanics , 1987 .

[4]  David Kinderlehrer,et al.  Mathematical Questions of Liquid Crystal Theory , 1987 .

[5]  A. Isihara,et al.  Theory of Liquid Crystals , 1972 .

[6]  F. Alouges.  A New Algorithm For Computing Liquid Crystal Stable Configurations: The Harmonic Mapping Case , 1997 .

[7]  Rolf Rannacher,et al.  A Feed-Back Approach to Error Control in Finite Element Methods: Basic Analysis and Examples , 1996 .

[8]  Robert V. Kohn,et al.  Singular Perturbation and the Energy of Folds , 2000, J. Nonlinear Sci..

[9]  Rolf Rannacher,et al.  An optimal control approach to a posteriori error estimation in finite element methods , 2001, Acta Numerica.

[10]  Chun Liu,et al.  Simulations of singularity dynamics in liquid crystal flows: A C0 finite element approach , 2006, J. Comput. Phys..

[11]  F. C. Frank,et al.  I. Liquid crystals. On the theory of liquid crystals , 1958 .

[12]  Chun Liu,et al.  Approximation of Liquid Crystal Flows , 2000, SIAM J. Numer. Anal..

[13]  Yoshikazu Giga,et al.  A mathematical problem related to the physical theory of liquid crystal configurations , 1987 .

[14]  Xing-Bin Pan,et al.  Asymptotics of minimizers of variational problems involving curl functional , 2000 .

[15]  Mitchell Luskin,et al.  Relaxation and gradient methods for molecular orientation in liquid crystals , 1989 .

[16]  Roland Glowinski,et al.  An operator-splitting method for a liquid crystal model , 2003 .

[17]  David Kinderlehrer,et al.  Existence and partial regularity of static liquid crystal configurations , 1986 .

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

[19]  F. Lin.  Nonlinear theory of defects in nematic liquid crystals; Phase transition and flow phenomena , 1989 .

[20]  Yoshikazu Giga,et al.  On lower semicontinuity of a defect energy obtained by a singular limit of the Ginzburg–Landau type energy for gradient fields , 1997, Proceedings of the Royal Society of Edinburgh: Section A Mathematics.

[21]  Jie Shen,et al.  Corrigendum: Fourier Spectral Approximation to a Dissipative System Modeling the Flow of Liquid Crystals , 2001, SIAM J. Numer. Anal..