Computational fluid dynamics for nematic liquid crystals

Due to recent advances in fast iterative solvers in the field of computational fluid dynamics, more complex problems which were previously beyond the scope of standard techniques can be tackled. In this paper, we describe one such situation, namely, modelling the interaction of flow and molecular orientation in a complex fluid such as a liquid crystal. Specifically, we consider a nematic liquid crystal in a spatially inhomogeneous flow situation where the orientational order is described by a second rank alignment tensor. The evolution is determined by two coupled equations: a generalised Navier-Stokes equation for flow in which the divergence of the stress tensor also depends on the alignment tensor and its time derivative, and a convection-diffusion type equation with non-linear terms that stem from a Landau-Ginzburg-DeGennes potential for the alignment. In this paper, we use a specific model with three viscosity coefficients that allows the contribution of the orientation to the viscous stress to be cast in the form of an orientation-dependent force. This effectively decouples the flow and orientation, with each appearing only on the right-hand side of the other equation. In this way, difficulties associated with solving the fully coupled problem are circumvented and a stand-alone fast solver, such as the state-of-the-art preconditioned iterative solver implemented here, can be used for the flow equation. A time-discretised strategy for solving the flow-orientation problem is illustrated using the example of Stokes flow in a lid-driven cavity.