The two-equation k-e model, widely employed in the analysis of turbulent flows, is not capable of adequately modelling problems involving secondary and swirling flows in ducts, boudary-layer detachment and other situations in which the inherent anisotropy of turbulent flows must be taken into account. More complex models, that take this anisotropy into account - the so-called Reynolds-stress models - are employed with the purpose of producing numerical results closer to those obtained experimentally. A geometrically simple problem, the turbulent flow in a duct with a square cross-section, in which the presence of secondary flows was observed experimentally and documentd by several authors, was modelletd and resolved using the Finite Volume Method. Initially, the k-e model was implemtend, being proven not capable of predicting, due to its isotropic nature, the secondary flows in a duct cross-section. The Reynolds-stress model was then implemented. The validation of this model is obtained through comparison of the numerical resuls to experimental and numerical results found in the bibliography.