Numerical analysis for constrained and unconstrained Q-tensor energies for liquid crystals

TL;DR

This paper develops a finite element framework for accurately approximating 3D Landau-de Gennes Q-tensor energies in nematic liquid crystals, accounting for anisotropic elasticity and physical eigenvalue constraints.

Contribution

It introduces a novel finite element approach that handles eigenvalue constraints and anisotropic energies, with rigorous analysis of well-posedness and convergence.

Findings

Proves well-posedness of the discrete problem

Establishes optimal convergence rates in energy norm

Analyzes impact of eigenvalue constraints on error estimates

Abstract

This paper introduces a comprehensive finite element approximation framework for three-dimensional Landau-de Gennes $Q$-tensor energies for nematic liquid crystals, with a particular focus on the anisotropy of the elastic energy and the Ball-Majumdar singular potential. This potential imposes essential physical constraints on the eigenvalues of the $Q$-tensor, ensuring realistic modeling. We address the approximation of regular solutions to nonlinear elliptic partial differential equations with non-homogeneous boundary conditions associated with Landau-de Gennes energies. The well-posedness of the discrete linearized problem is rigorously demonstrated. The existence and local uniqueness of the discrete solution is derived using the Newton-Kantorovich theorem. Furthermore, we demonstrate an optimal order convergence rate in the energy norm and discuss the impact of eigenvalue constraints on the a priori error analysis.