Simulations of solutions of flexible polymer molecules during flow in simple or complex confined geometries are performed. Concentrations from ultradilute up to near the overlap concentration are considered. As concentration increases, the hydrodynamic migration effects observed in dilute solution unidirectional flows (Couette flow, Poiseuille flow) become less prominent, virtually vanishing as the overlap concentration is approached. In a grooved channel geometry, the groove is almost completely depleted of polymer chains at high Weissenberg number in the dilute limit, but at finite concentration this depletion effect is dramatically reduced. Only upon inclusion of hydrodynamic interactions can these phenomena be properly captured.