This paper presents the development of spectral element methods to simulate unsteady flows of viscoelastic fluids using a differential constitutive equation. Preliminary results for simple test problems are discussed as well as possible techniques to improve the numerical stability of the computational algorithms.