This paper describes the fundamental aspects of turbulence modelling for incompressible fluid flow and corresponding numerical models based on the boundarydomain integral equations. The velocity-vorticity formulation of the Reynolds averaged Navier-Stokes (RANS) equations is used, while the averaged pressure field is determined by solving the Poisson velocity equation. Several lowReynolds-number turbulent models, e.g. Launder-Sharma, are applied.