Abstract:We present a second-order strictly length-preserving and unconditionally energy-stable rotational discrete gradient (Rdg) scheme for the numerical approximation of the Oseen-Frank gradient flows with anisotropic elastic energy functional. Two essential ingredients of the Rdg method are reformulation of the length constrained gradient flow into an unconstrained rotational form and discrete gradient discretization for the energy variation. Besides the well-known mean-value and Gonzalez discrete gradients, we propose a novel Oseen-Frank discrete gradient, specifically designed for the solution of Oseen-Frank gradient flow. We prove that the proposed Oseen-Frank discrete gradient satisfies the energy difference relation, thus the resultant Rdg scheme is energy stable. Numerical experiments demonstrate the efficiency and accuracy of the proposed Rdg method and its capability for providing reliable simulation results with highly disparate elastic coefficients.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to design a numerical method that can strictly preserve the vector length and is unconditionally energy - stable for the Oseen - Frank gradient flow with anisotropic elastic energy functionals. Specifically, the author proposes a new second - order Rotational Discrete Gradient (Rdg) method for numerically approximating the Oseen - Frank gradient flow. ### Background and Problem Description - **Research Background**: Liquid crystal dynamics involves the local anisotropic evolution generated by the non - uniform orientation distribution. For uniaxial nematic liquid crystals, the equilibrium orientation distribution is symmetric and can freely rotate as a whole. This evolution can be described by the unit vector field \( \mathbf{n}(x) \). - **Deficiencies of Existing Methods**: Most of the existing numerical methods mainly study isotropic elasticity (i.e., \( k_1 = k_2 = k_3 = k \) and \( k_4 = 0 \)), and there are fewer studies on the case of anisotropic elasticity. Even under the assumption of isotropic elasticity, it is very difficult to achieve a numerical method that simultaneously satisfies the vector length constraint and energy dissipation. - **Specific Problem**: How to design a numerical method that can both preserve the vector length and ensure energy stability to deal with the Oseen - Frank gradient flow with anisotropic elasticity. ### Solution - **Method Overview**: The author proposes a new second - order Rotational Discrete Gradient (Rdg) method. By reformulating the gradient flow with unit vector constraints into an unconstrained rotational form and using the discrete gradient to handle the energy variation. - **Key Technologies**: - **Rotational Form**: Reformulate the gradient flow with unit vector constraints into an unconstrained rotational form, making the length constraint intrinsic. - **Discrete Gradient**: Propose three different discrete gradient methods, including the commonly used mean discrete gradient and Gonzalez discrete gradient, as well as the discrete gradient specifically designed for the Oseen - Frank energy. - **Energy Stability and Length Preservation**: Prove that the proposed Rdg method is strictly length - preserving and unconditionally energy - stable. ### Numerical Experiments - **Convergence Test**: Verify the convergence rate of the Rdg method in space and time by the method of manufactured solutions. - **Property Preservation Test**: Verify the energy stability and length - preservation of the Rdg method through numerical experiments. - **Adaptive Time - step Strategy**: Propose an adaptive time - step strategy to improve computational efficiency and accuracy. ### Conclusion - **Advantages of the Method**: The proposed Rdg method can not only strictly preserve the vector length but also unconditionally maintain energy stability, and is suitable for dealing with the Oseen - Frank gradient flow with anisotropic elasticity. - **Future Work**: This method is expected to become the basis of an efficient and robust numerical method for the Ericksen - Leslie equations and be further applied to more complex liquid crystal dynamics problems. Through these techniques, the paper provides a reliable and efficient numerical method, providing a powerful tool for studying the dynamics of nematic liquid crystals formed by different materials.

A second-order length-preserving and unconditionally energy stable rotational discrete gradient method for Oseen-Frank gradient flows

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

Energy diminishing implicit-explicit Runge–Kutta methods for gradient flows

Energy diminishing implicit-explicit Runge--Kutta methods for gradient flows

High-order energy stable discrete variational derivative schemes for gradient flows

Two novel numerical methods for gradient flows: generalizations of the Invariant Energy Quadratization method

Energy-dissipative spectral renormalization exponential integrator method for gradient flow problems

Higher-Order Energy-Decreasing Exponential Time Differencing Runge-Kutta methods for Gradient Flows

A new flow dynamic approach for Wasserstein gradient flows

Projection-Free Method for the Full Frank-Oseen Model of Liquid Crystals

Variational Extrapolation of Implicit Schemes for General Gradient Flows

A weighted scalar auxiliary variable method for solving gradient flows: bridging the nonlinear energy-based and Lagrange multiplier approaches

A General Class of Linear Unconditionally Energy Stable Schemes for the Gradient Flows

Unconditionally energy stable DG schemes for the Swift-Hohenberg equation

Solution of the Oseen Flow With the Multigrid Method Based on Alternating-Line Relaxation

Global-in-time energy stability: a powerful analysis tool for the gradient flow problem without maximum principle or Lipschitz assumption

A novel high-order linearly implicit and energy-stable additive Runge-Kutta methods for gradient flow models

The scalar auxiliary variable (SAV) approach for gradient flows

Uniformly high-order bound-preserving OEDG schemes for two-phase flows

A New Scalar Auxiliary Variable Approach for Gradient Flows