Discontinuous Galerkin finite element methods for the Landau–de Gennes minimization problem of liquid crystals

Abstract

We consider a system of second-order nonlinear elliptic partial differential equations that models the equilibrium configurations of a two-dimensional planar bistable nematic liquid crystal device. Discontinuous Galerkin (dG) finite element methods are used to approximate the solutions of this nonlinear problem with nonhomogeneous Dirichlet boundary conditions. A discrete inf–sup condition demonstrates the stability of the dG discretization of a well-posed linear problem. We then establish the existence and local uniqueness of the discrete solution of the nonlinear problem. A priori error estimates in the energy and L2 norms are derived and a best approximation property is demonstrated. Further, we prove the quadratic convergence of the Newton iterates along with complementary numerical experiments.