Presented in the paper is an efficient and accurate numerical method for simulation of ventilation duct flow. The mathematical method is based on the three-dimensional incompressible RANS equations with isotropic k-w near-wall turbulence closures, written in generalized curvilinear coordinates in strong conservation form. The numerical method presented here is used to calculate the turbulent flow through a bend of rectangular ventilation duct featuring pressure induced secondary motions and rotation effects on turbulence.