Development of a computer program for three-dimensional viscous flow

This paper describes a numerical calculation procedure which has been embodied in a computer code for the computation of three-dimensional viscous flows governed by partial differential equations of the elliptic type. The solution technique employs velocities and pressure as the main hydrodynamic variables as well as an implicit, hybrid finite-difference/finite-volume method. It can handle steady incompressible flow with heat/mass transfer which is laminar or turbulent. In the case of turbulent flow problems, the 'k-epsilon' two-equation turbulence model is used. The applicability of this computational procedure is demonstrated by the solution of three flow problems; (1) laminar flow in a square duct, (2) a turbulent jet in a cross flow and (3) the calculation of the flow around a building with a short adjacent stack. The comparisons between predictions and experimental data are fairly good.