Solution (source code)

= Solution

Steadiness, axial and rotational symmetry, and the absence of any imposed radial or azimuthal motion make the velocity independent of $z$ and $\theta$ and directed along the cylinders. Incompressibility is then automatic, so
$$
\boxed{\mathbf u=u(r)\mathbf e_z}.
$$