A method for orthogonal grid generation is presented. 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 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.