A fast and robust method of grid generation to multiple functions has been developed for flow analysis in three dimensional space. It is based on the Non-Uniform Rational B-Spline(NURBS) of an approximation method. Many of NURBS intrinsic properties are introduced and much more easily understood. The grid generation method, details of numerical implementation. examples of application, and potential extensions of the current method are illustrated in this paper. The object of this study is to develop the surface grid generation and the grid cluster techniques capable of resolving complex flows with shock waves, expansion waves, shear layers. The knot insert method of Non-Uniform Rational B-Spline seems well worked. In addition, NURBS has been widely utilized to generate grids in the computational fluid dynamics community. Computational examples associated with practical configurations have shown the utilization of the algorithm.