A method for nearly orthogonal grid generation is presented in this study. The generating system is based on solution of a system of partial differential equations with finite difference discretization. The influence of the number of grid points, type of boundary, and intensity of the grid quality control function and grid properties are investigated. Specification of both boundary point distribution on all sides and moving boundaries is used. The proposed method is applied to various test problems£¬which shows this method provides a good balance between controlling grid orthogonality and cell aspect ratio.