A finite element method for nematic liquid crystals with variable degree of orientation

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

DOI: https://doi.org/10.1137/15M103844X

2017-08-03

Abstract: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 $\Gamma$-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)," <a class="link-external link-http" href="http://www.youtube.com/watch?v=pWWw7" rel="external noopener nofollow">this http URL</a>\_6cQ-U

Numerical Analysis,Analysis of PDEs

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: for nematic liquid crystal materials with a variable degree of orientation, how to design a finite element method (FEM) without regularization to simulate their equilibrium state configurations and deal with the resulting line and surface defects. Specifically, the paper focuses on the simplest one - constant model proposed by Ericksen, in which the equilibrium state of the liquid crystal is jointly determined by the director field \( \mathbf{n} \) that describes the molecular orientation and the scalar function \( s \) that represents the degree of orientation. The pair \((s, \mathbf{n})\) minimizes a total energy that includes Frank - like energy and double - well potential energy.

### Key Problem Analysis

1. **Energy Minimization Problem**:
- The equilibrium state of the liquid crystal is obtained by minimizing the following energy functional:
\[
E[s, \mathbf{n}] := \int_{\Omega} \left( \kappa |\nabla s|^2 + s^2 |\nabla \mathbf{n}|^2 \right) \, dx + \int_{\Omega} \psi(s) \, dx
\]
where \( \kappa > 0 \), and \( \psi(s) \) is a double - well potential function that satisfies specific conditions (such as having a minimum at \( s = 0 \)).

2. **Degenerate Elliptic Equation**:
- The \( \mathbf{n} \) equation in the Euler - Lagrange equation is a degenerate elliptic equation, which allows the existence of line and surface defects. These defects have finite energy when \( s = 0 \).

3. **Finite Element Discretization**:
- The paper proposes a structure - preserving discretization method. It uses piecewise linear finite elements to handle the degenerate elliptic part without regularization. This method ensures the consistency and stability of the discrete energy and proves that the discrete global minimum solution converges to the continuous solution as the mesh size approaches zero.

4. **Gradient Flow Scheme**:
- A quasi - gradient flow scheme is developed to calculate the discrete equilibrium solution, and it is proved that this scheme has a strictly monotonically decreasing energy property.

5. **Numerical Experiments**:
- Two - dimensional and three - dimensional numerical simulations are carried out, showing that this method can handle non - trivial defects and verifying its effectiveness and accuracy.

### Summary

The main contribution of the paper is to provide an effective finite element method that can handle complex defect structures in nematic liquid crystal materials with a variable degree of orientation without introducing additional regularization. This method not only theoretically guarantees the correctness of the energy - minimizing solution but also shows good numerical performance in practical applications.