Numerical simulation of complex viscoelastic flows using discontinuous galerkin spectral/hp element methods