We investigate the migration of bubbles in several flow patterns occurring within the gap between a rotating inner cylinder and a concentric fixed outer cylinder. The time-dependent evolution of the two-phase flow is predicted through three-dimensional Euler-Lagrange simulations. Lagrangian tracking of spherical bubbles is coupled with direct numerical simulation of the Navier-Stokes equations. We assume that bubbles do not influence the background flow ͑one-way coupling simulations͒. The force balance on each bubble takes into account buoyancy, added-mass, viscous drag, and shear-induced lift forces. For increasing velocities of the rotating inner cylinder, the flow in the fluid gap evolves from the purely azimuthal steady Couette flow to Taylor toroidal vortices and eventually a wavy vortex flow. The migration of bubbles is highly dependent on the balance between buoyancy and centripetal forces ͑mostly due to the centripetal pressure gradient͒ directed toward the inner cylinder and the vortex cores. Depending on the rotation rate of the inner cylinder, bubbles tend to accumulate alternatively along the inner wall, inside the core of Taylor vortices or at particular locations within the wavy vortices. A stability analysis of the fixed points associated with bubble trajectories provides a clear understanding of their migration and preferential accumulation. The location of the accumulation points is parameterized by two dimensionless parameters expressing the balance of buoyancy, centripetal attraction toward the inner rotating cylinder, and entrapment in Taylor vortices. A complete phase diagram summarizing the various regimes of bubble migration is built. Several experimental conditions considered by Djéridi, Gabillet, and Billard ͓Phys.