This paper presents a finite volume numerical method to calculate turb
ulent flows using second-moment turbulence closures. The numerical met
hod is formulated for a collocated grid arrangement and on body fitted
coordinate grids. A special interpolation scheme is derived to resolv
e the potential decoupling between momentum and Reynolds stress transp
ort equations. This interpolation scheme naturally leads to stability-
promoting terms in the momentum equations in the form of anisotropic e
ddy viscosities. Issues such as the implementation of boundary conditi
ons with a body-fitted coordinate grid are also discussed. Application
s are given to demonstrate the robust ness of the numerical method and
the performance of a second moment turbulence closure. Computational
results are compared with available experimental data and to those of
the standard two equation turbulence model. Good results have been obt
ained with the second moment turbulence closure.