$Q$-tensor gradient flow with quasi-entropy and discretizations preserving physical constraints

Abstract

We propose and analyze numerical schemes for the gradient flow of $Q$-tensor with the quasi-entropy. The quasi-entropy is a strictly convex, rotationally invariant elementary function, giving a singular potential constraining the eigenvalues of $Q$ within the physical range $(-1/3,2/3)$. Compared with the potential derived from the Bingham distribution, the quasi-entropy has the same asymptotic behavior and underlying physics. Meanwhile, it is very easy to evaluate because of its simple expression. For the elastic energy, we include all the rotationally invariant terms. The numerical schemes for the gradient flow are built on the nice properties of the quasi-entropy. The first-order time discretization is uniquely solvable, keeping the physical constraints and energy dissipation, which are all independent of the time step. The second-order time discretization keeps the first two properties unconditionally, and the third with an $O(1)$ restriction on the time step. These results also hold when we further incorporate a second-order discretization in space. Error estimates are also established for time discretization and full discretization. Numerical examples about defect patterns are presented to validate the theoretical results.

arXiv comment25 pages, 7 figures

References in corpus (9)

• Construction of a pathway map on a complicated energy landscape
• Third order Maximum-Principle-Satisfying Direct discontinuous Galerkin methods for time dependent convection diffusion equations on unstructured triangle mesh
• The cutoff method for the numerical computation of nonnegative solutions of parabolic PDEs with application to anisotropic diffusion and lubrication-type equations
• Unconditionally positivity preserving and energy dissipative schemes for Poisson--Nernst--Planck equations
• On Minimizers of the Landau-de Gennes Energy Functional on Planar Domains
• A Fast Algorithm for the Moments of Bingham Distribution
• Quasi-entropy by log-determinant covariance matrix and application to liquid crystals
• Efficient Spectral Methods for Quasi-Equilibrium Closure Approximations of Symmetric Problems on Unit Circle and Sphere
• Regularity of a gradient flow generated by the anisotropic Landau-de Gennes energy with a singular potential