This article analyses stationary scalar mixing downstream an open flow Couette device operating in the creeping flow regime. The device consists of two coaxial cylinders of finite length Lz, and radii ¿ R and R (¿ <1), which can rotate independently. At relatively large values of the aspect ratio a = Lz/R 1, and of the Péclet number Pe, the stationary response of the system can be accurately described by enforcing the simplifying assumption of negligible axial diffusion. With this approximation, homogenization along the device axis can be described by a family of generalized one-dimensional eigenvalue problems with the radial coordinate as independent variable. A variety of mixing regimes can be observed by varying the geometric and operating parameters. These regimes are characterized by different localization properties of the eigenfunctions and by different scaling laws of the real part of the eigenvalues with the Péclet number. The analysis of this model flow reveals the occurrence of sharp transitions between mixing regimes, e.g. controlled by the geometric parameter ¿. The eigenvalue scalings can be theoretically predicted by enforcing eigenfunction localization and simple functional equalities relating the behaviour of the eigenvalues to the functional form of the associated eigenfunctions. Several flow protocols corresponding to different geometric and operating conditions are considered. Among these protocols, the case where the inner and the outer cylinders counter-rotate exhibits a peculiar intermediate scaling regime where the real part of the dominant eigenvalue is independent of Pe over more than two decades of Pe. This case is thoroughly analysed by means of scaling analysis. The practical relevance of the results deriving from spectral analysis for fluid mixing problems in finite-length Couette devices is addressed in detail.