  We describe a lattice Boltzmann algorithm to simulate liquid crystal
hydrodynamics in three dimensions. 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 viscoelastic effects
such as shear-thinning and shear-banding. We describe the implementation of
velocity boundary conditions and show that the algorithm can be used to
describe optical bounce in twisted nematic devices and secondary flow in
sheared nematics with an imposed twist.
