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. .