SUMMARYIncompressible turbulent ow past a long square cylinder is investigated using large eddy simulations (LES). Results are presented and compared with available experimental databases for a Reynolds number Re d = 22 000. The problem served as one of the validation cases for the development of a numerical code designed for e cient, parallel, three-dimensional N-S computations in complex geometrical conÿgurations. In contrast with previous studies, the geometrical deÿnition of the problem is established by the immersed boundary concept (IMB) while pressure solution is performed by a fast, fully parallel direct pressure solver. Calculations were performed with the widely applied Smagorinsky turbulence model and the ÿltered structure function model (FSF) which has not been previously applied to the ow case under consideration. In order to assess the potential of LES at its lowest (RANS), and highest (DNS) limit, di erent numerical resolutions were examined. Depending on the available resolution, either no-slip conditions or a modiÿed Werner and Wengle approximate wall boundary condition was used. The predicted mean velocity and uctuation proÿles, force statistics and Strouhal numbers were found to be in very good agreement with the experimental data sets. Analysis of the results indicates that for time varying blu -bodies ows that involve complex ow phenomena, successful large eddy simulations are not just possible, but can also achieve an excellent quality of results at a relatively low cost.