This paper presents the application of spectral element methods to simulate the time-dependent flow of viscoelastic fluids in non-trivial geometries using a closed-form differential constitutive equation. As an example, results relative to the flow of a FENE-CR fluid in a two-dimensional four-to-one contraction are given.