A V-cycle multigrid method is developed for a time-dependent viscoelastic fluid flow satisfying an Oldroyd-B-type constitutive equation in two-dimensional domains. Also existence, uniqueness, and error estimates of an approximate solution are discussed. The approximate stress, velocity, and pressure are, respectively, σ
k
-discontinuous, u
k
-continuous, and p
k
-continuous.
