While acoustic streaming patterns in microfluidic channels with rectangular cross-sections have been widely shown in the literature, boundary-driven streaming fields in non-rectangular channels have not been well studied. In this paper, a two-dimensional numerical model was developed to simulate the boundary-driven streaming fields on cross-sections of cylindrical fluid channels. Firstly, the linear acoustic pressure fields at the resonant frequencies were solved from the Helmholtz equation. Subsequently, the outer boundary-driven streaming fields in the bulk of fluid were modelled while using Nyborg’s limiting velocity method, of which the limiting velocity equations were extended to be applicable for cylindrical surfaces in this work. In particular, acoustic streaming fields in the primary (1, 0) mode were presented. The results are expected to be valuable to the study of basic physical aspects of microparticle acoustophoresis in microfluidic channels with circular cross-sections and the design of acoustofluidic devices for micromanipulation.