A second-order SO(3)-preserving and energy-stable scheme for orthonormal frame gradient flow model of biaxial nematic liquid crystals

Hanbin Wang,Jie Xu,Zhiguo Yang
2023-10-16
Abstract:In this paper, we present a novel second-order generalised rotational discrete gradient scheme for numerically approximating the orthonormal frame gradient flow of biaxial nematic liquid crystals. This scheme relies on reformulating the original gradient flow system into an equivalent generalised "rotational" form. A second-order discrete gradient approximation of the energy variation is then devised such that it satisfies an energy difference relation. The proposed numerical scheme has two remarkable properties: (i) it strictly obeys the orthonormal property of the tensor field and (ii) it satisfies the energy dissipation law at the discrete level, regardless of the time step sizes. We provide ample numerical results to validate the accuracy, efficiency, unconditional stability and SO(3)-preserving property of this scheme. In addition, comparisons of the simulation results between the biaxial orthonormal frame gradient flow model and uniaxial Oseen-Frank gradient flow are made to demonstrate the ability of the former to characterize non-axisymmetric local anisotropy.
Numerical Analysis,Mathematical Physics
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: **How to design an efficient and robust numerical method to solve the orthogonal frame gradient flow model of biaxial nematic liquid crystals while maintaining orthogonality and the energy dissipation law**. Specifically, the author focuses on the following problems: 1. **Challenges of nonlinear and coupled systems**: The orthogonal frame gradient flow system of biaxial nematic liquid crystals is a highly nonlinear and coupled system, involving three unknown variables $\mathbf{n}_1, \mathbf{n}_2, \mathbf{n}_3$ and six point - constraint conditions. These nonlinearities and couplings come not only from the nonlinear terms in the energy variation of the biaxial elastic energy functional, but also from the explicit orthogonality constraints. 2. **Anisotropy of elastic coefficients**: The elastic coefficients in the free - energy density derived from molecular parameters may become highly anisotropic with the change of molecular structure, resulting in strong nonlinearity and coupling problems in numerical calculations. 3. **Maintenance of physical properties**: This system obeys two basic physical properties: the orthogonality constraint at points and the energy dissipation law. Therefore, the numerical method needs to maintain these two properties simultaneously at the discrete level. If the orthogonality constraint is violated, it may lead to unpredictable calculation results; and simple projection onto SO(3) may introduce numerical instability and cause an unreasonable increase in the energy of the gradient flow system. To solve these problems, the author proposes a new numerical method - **the generalised rotational discrete gradient (gRdg) method**, which can strictly maintain orthogonality and the energy dissipation law at the discrete level and has unconditional energy stability. ### Core ideas of the method - **Reformulating the gradient flow system**: Reformulate the gradient flow system with orthogonality constraints into a generalized "rotational" form. - **Discrete gradient approximation**: Develop a discrete gradient approximation method so that it satisfies a certain energy difference relationship. - **Time - centered Crank - Nicolson scheme**: Combine the time - centered Crank - Nicolson scheme to form the complete gRdg method. Through this method, the author proves that the gRdg scheme can strictly maintain the orthogonality constraint and satisfy the discrete energy dissipation law. In addition, they also propose an adaptive time - step strategy to improve the efficiency and accuracy of solving nonlinear equations. ### Numerical experiment verification To verify the effectiveness of this method, the author has carried out a large number of numerical experiments, including: - Tests of spatial and temporal convergence rates - Changes in orthogonality error and energy curves over time - Comparison with the uniaxial Oseen - Frank gradient flow model - Research on dynamic behaviors under highly anisotropic elastic coefficients These experimental results show that the gRdg method has high accuracy and efficiency in dealing with the complex dynamic behaviors of biaxial nematic liquid crystals.