To simulate the incompressible flow in complex three-dimensional geometry efficiently and accurately, a solver based on solution of the Navier-Stokes equations in the generalized curvilinear coordinate system was developed. The system of equations in three-dimension are solved simultaneously by the artificial compressibility method. The convective terms are differenced using a flux difference splitting approach. The viscous terms are differenced using second-order accurate central differences. An implicit line relaxation scheme is employed to solve the numerical system of equations. The solver was tested for two cases including flow past a circular cylinder and flow around a hemispherical head of a cylindrical object.
