Details of the turbulent flow in a 1:8 aspect ratio rectangular duct at a Reynolds number of approximately 5800 were investigated both numerically and experimentally. The three‐dimensional mean velocity field and the normal stresses were measured at a position 50 hydraulic diameters downstream from the inlet using laser doppler velocimetry (LDV). Numerical simulations were carried out for the same flow case assuming fully developed conditions by imposing cyclic boundary conditions in the main flow direction. The numerical approach was based on the finite volume technique with a non‐staggered grid arrangement and the SIMPLEC algorithm. Results have been obtained with a linear and a non‐linear (Speziale) k–ε model, combined with the Lam–Bremhorst damping functions for low Reynolds numbers. The secondary flow patterns, as well as the magnitude of the main flow and overall parameters predicted by the non‐linear k–ε model, show good agreement with the experimental results. However, the simulations provide less anisotropy in the normal stresses than the measurements. Also, the magnitudes of the secondary velocities close to the duct corners are underestimated. © 1998 John Wiley & Sons, Ltd.