A finite difference numerical method for the solution of the two-dimensional incompressible, steady Navier-Stokes equation in general curvilinear coordinates is presented, and pressure has been selected as a primary dependent variable.The procedure is developed on the basis of the non-orthogonal, body-fitted coordinate. Because of using non-staggered grids, that all flow variables are stored at the same set of nodes. Further more the Cartesian velocity component has been used to make the difference approximation equations simple and to enhance the conservative property of the equation.To suppress the pressure oscillations while using non-staggered grid arrangement, the pressure-weighted interpolation method, proposed by Rhie and Chow, was conducted. And several common difference approximations of the convective-diffusive terms are compared in order to verify the method.The k - s model is utilized to describe the turbulent flow process.In addition to above, several methods of grid generation is also discussed. The governing equations used in elliptic partial differential equation system are given out in detail, and the geometry method of grid generation is used here.The results of calculations are compared with the available experimental data and other numerical results available in the literatures. As a conclusion ,all results show that the method presented is available and efficient in the Computational Fluid Dynamics. |