Stokes flow in a cylindrical column of fluid, with a stress-free cylindrical sidewall, is considered. The motion is assumed to be generated by the linear, uniform motion of either or both of the flat endwalls. The field is obtained by a vector eigenfunction expansion procedure. If the field is assumed to have a theta- and z-dependence of the type exp(i theta + kz), k has to satisfy the equation (1 - k(2)) J(1)(3) + j(1) [2j(1)(2) - J(1) {2J(1) + (1 + k(2)) J(1)}] = 0, where J(1) = k J(0) - J(1) the argument of each Bessel function-is k. This equation admits, unlike in the plane case and with important consequences, not just a real sequence {lambda(n)} of eigenvalues but also a complex one [,). Using a least squares procedure to satisfy the boundary conditions on the top and bottom walls, the three-dimensional velocity field in the column is determined for various values of column height h and wall speed ratio S. Detailed computations show that there are strong effects of both the stress-free boundary and three-dimensionality. The principal effect of the former is to permit motion on that boundary leading to large azimuthal motions and of the latter, unlike in the plane flow, to multiple primary eddies when h is sufficiently large. A number of new eddy structures are also found, which demonstrate that three-dimensionality often leads to the elimination of compactness found in plane flows. It is finally shown that the flow fields exhibit interesting bifurcations as S and h are varied. (C) 2007 The Japan Society of Fluid Mechanics and Elsevier B.V. All rights reserved.