The non-linear k-.epsilon. model developed by Speziale was employed for the prediction of developing turbulent flow in a square duct. The numerical procedure incorporated a finite volume method using a strong conservation form of the partially-parabolized Navier-Stokes equation. Results of the calculation were compared with available experimental data on the mean velocity field and turbulent kinetic energy, and was found to be in favorable agreement.