摘要

We present a preconditioned nullspace method for the numerical solution of large sparse linear systems that arise from discretizations of continuum models for the orientational properties of liquid crystals. The approach effectively deals with pointwise unit-vector constraints, which are prevalent in such models. The indefinite, saddle-point nature of such problems, which can arise from either or both of two sources (pointwise unit-vector constraints, coupled electric fields), is illustrated. Both analytical and numerical results are given for a model problem.

  • 出版日期2013