Generalized Scalar Auxiliary Variable Exponential Integrator for A Modified Landau-de Gennes Theory for Smectic Liquid Crystals
This work provides a robust numerical method for simulating complex liquid crystal dynamics, addressing stability and accuracy issues in a specific domain (smectic liquid crystals).
The authors propose a novel, efficient, and unconditionally energy-stable numerical scheme for a coupled system modeling Smectic-A liquid crystals, achieving optimal error estimates and demonstrating accuracy and efficiency through numerical experiments.
The Smectic-A (SmA) phase is modeled by a modified Landau-de Gennes (mLdG) model proposed by Xia et al. [Phys. Rev. Lett., 126 (2021), 177801], in which a tensor order parameter $\mathbf{Q}$ for the orientational order is coupled with a real scalar $u$ characterizing the positional order. In this paper, we propose and analyze a novel, highly efficient, and unconditionally energy-stable numerical scheme for this coupled system by combining the generalized scalar auxiliary variable-exponential integrator (GSAV-EI) approach with a relaxed correction strategy. In particular, we reformulate the exponential time differencing time discretization into an equivalent quasi-implicit backward Euler-type structure, a pivotal step that eliminates the restrictive CFL mesh-ratio conditions of the original GSAV-EI method and enables a rigorous fully discrete error analysis. Theoretically, we rigorously establish the unconditional energy stability with respect to a modified discrete energy and the uniform boundedness of the numerical solutions $\mathbf{Q}$, along with optimal error estimates in both time and space. Comprehensive numerical experiments are presented to demonstrate the accuracy, efficiency, and structural preservation of the algorithm, as well as its capability in capturing complex topological defect dynamics.