Turbulent flow past a square cylinder confined in a channel is numerically investigated by Large Eddy Simulation(LES). The main objectives of this study are to verify the experimental results of Nakagawa et al.[Exp. in Fluids, Vol. 27, 3, pp. 284∼294, 1999] by LES and to obtain related flow information in detail. The LES results obtained are in excellent agreement with the experiment both qualitatively and quantitatively. The passive paticles numerically released into the flow field clearly show the barman vortex street. However, the vortices shed from the cylinder are significantly affected by the presence of the plates. Futhermore, periodic and alternating vortex-rollups are observed in the vicinity of the plates. The rolled-up vortex is convected downstream together with the corresponding Karman vortex forming a counter-rotating vortex pair. It is also revealed that the cylinder greatly enhances mixing process of the flow.