Turbulent duct flows are investigated using large eddy simulation at bulk Reynolds numbers, from 4410 to 250000. Mean secondary flow are found to reveal the existence of two steamwise counterrotating vortices in each corner of the duct. Turbulence-driven secondary motions that arise in duct flows act to transfer fluid momentum from the centre of the duct to its corners, thereby causing a bulging of the streamwise velocity contours towards the corners. As Reynolds number increases, the ratio of centerline streamwise velocity to the bulk velocity decreases and all turbulent components increase. In addition, the core of the secondary vortex in the lower corner-bisector tends to approach the wall and the corner with increasing Reynolds number. The turbulence intensity profiles for the low Reynolds number flows are quite different from those for the high Reynolds number flows.Typical turbulence structures in duct flows are found to be responsible for the interactions between ejections from wall and this interaction results in the bending of the ejection stems, which indicates that the existence of streaky wall structures is much like in a channel flow.