Spectral direct numerical simulations (DNS) are carried out for a source–sink flow in an annular cavity between two co-rotating disks. When the Reynolds number based on the forced inflow is increased, a self-sustained crossflow instability of finite amplitude is observed. We show that this nonlinear global mode is made up of a front located at the upstream boundary of the absolutely unstable domain, followed by a saturated spiral mode, and that its properties are in good agreement with results of the local stability theory. This structure is characteristic of the so-called elephant mode of Pier & Huerre (J. Fluid Mech. vol. 435, 2001, p. 145). The global bifurcation is subcritical since only large-amplitude initial perturbations are fou...