Spectral collocation/domain decomposition method for viscoelastic flow simulations in model porous geometries