No abstract

The hydrodynamic Q-tensor model has been used for studying flows of liquid crystals and liquid crystal polymers. It can be derived from a variational approach together with the generalized Onsager principle, in which the total energy decreases in time. In this paper, we develop a novel, linear, second order semi-discrete scheme in time to solve the governing system. The scheme is developed following the so called ' energy quadratization ' strategy so that it is linear and unconditionally energy stable at the semi-discrete level. This scheme is further discretized in space using a second order finite difference method and implemented on a GPU for high performance computing. The convergence rate in time is established using a mesh refinement test. Several numerical examples are presented to demonstrate the usefulness of the model and the effectiveness of the numerical scheme in simulating defect dynamics in flows of liquid crystals.

No abstract