A finite-element model of liquid-crystal hydrodynamics based on the Qian and Sheng formulation has been developed. This formulation is a generalization of the Ericksen-Leslie theory to include variations in the order parameter, allowing for a proper description of disclinations. The present implementation is well suited to treat properly the various length scales necessary to model large regions yet resolve the rapid variations in the order parameter in proximity to disclinations.Index Terms-Finite elements, liquid-crystal (LC) modeling, variable order parameter.