We describe a lattice Boltzmann algorithm to simulate liquid crystal hydrodynamics. The equations of motion are written in terms of a tensor order parameter. This allows both the isotropic and the nematic phases to be considered. Backflow effects and the hydrodynamics of topological defects are naturally included in the simulations, as are non-Newtonian flow properties such as shear thinning and shear banding.

• Received 16 August 2000

DOI:https://doi.org/10.1103/PhysRevE.63.056702