Abstract
In this paper, a new solving method for the velocity components of the three-dimensional κ-ε model equations based on generalized curvilinear coordinates are presented, which imposes the Dirichlet-type boundary condition on the momentum equations automatically through the procedure of iterative calculations using contravariant vector components. Three-dimensional turbulent flows are analyzed using the present method. Two examples of numerical simulations are presented. First, to confirm the applicability of the present method, numerical simulations of the room air flow are conducted. The room space is cubic. The simulations are then compared with results obtained by the existing numerical method, which is based on Cartesian coordinates using a staggered grid system, as well as with the experimental results. Good agreement is shown in these comparisons. Second, a simulation of air flow in a gymnasium, which has a complicated globular boundary configuration, is performed. The results clearly demonstrate the practicability of the present method based on the curvilinear coordinate system. It gives reasonable air flow distributions for a three-dimensional non-rectangular space.
