摘要

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 of equations in the model. The scheme is developed following the novel '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.

  • 出版日期2017-5-1
  • 单位北京计算科学研究中心; 南开大学