The axial dispersion of approximately monosized particles in rolling mode in rotating cylinders with bulk flow is examined using a Monte Carlo model and discrete element method (DEM) simulations. The Monte Carlo model predicts that the mean square displacement relative to the mean axial displacement of the bed undergoes oscillations in time. The nature of these oscillations depends on the fill level of the cylinder and the extent of particle mixing during avalanches. When the cylinder is half full the Monte Carlo model predicts undamped oscillations, whereas a filling fraction of 0.26 produces oscillations whose amplitude decreases with time. If mixing during avalanches is assumed to be perfect then the oscillations occur about a linear increase with time. In contrast, if it is assumed that the particles do not mix during avalanching, the oscillations occur about an increase with time which has a gradient which increases with time. There is good qualitative agreement between the Monte Carlo model with perfect mixing and the DEM when the filling fraction is 0.26. For a filling fraction of 0.5 the DEM data show oscillations about a faster than linear increase with time.