After a brief review of the more popular turbulence models, the author presents and discusses accurate and efficient numerical methods for solving the boundary-layer equations with turbulence models based on algebraic formulas (mixing length, eddy viscosity) or partial-differential transport equations. A computer program employing the Cebeci-Smith model and the k-e model for obtaining the solution of two-dimensional incompressible turbulent flows without separation is discussed in detail and is presented in the accompanying CD.