Numerical calculation of Reynolds stresses in a square duct with secondary flow

A new model for the Reynolds stress equations is presented. This model is used to obtain a theoretical solution for the problem of fully developed turbulent flow in a square duct. Nine governing equations for the axial velocity, lateral vorticity, lateral stream function and six components of the Reynolds stresses are simultaneously solved, by a finite-difference technique. To ensure numerical stability of the solution a special linearised implicit representation of the source terms is proposed, and simultaneous solution of the equations at each.mesh point is obtained. Near the wall a special procedure is used, by which the Reynolds stress equations are assumed to be in local equilibrium, and the velocity profile is assumed to be logarithmic. However, due to the secondary motion the logarithmic velocity profile is inclined to the axial direction. The results bear reasonable agreement with experimental data. Computer time requirements are moderate.