A numerical method for 2D orthogonal grid generation with control of the boundary point distribution is presented. The Ryskin and Leal classified their method in two categomethod is based on the solution of a system of partial differential ries: the ''weak constraint'' and ''strong constraint'' methequations. The grid cell aspect ratio, the so-called distortion ods, according to the type of boundary conditions adopted function, is calculated from its definition equation in the entire [4]. The ''strong constraint'' method applies to domains domain. The method allows the specification of the boundary whose shape is not known in advance and the ''weak conpoint distribution in all the boundaries and also the distance of the first grid node to the boundary. The method is successfully straint'' method, to domains with a known shape. The applied to several geometries, including geometries with grid latter one has been the subject of several papers [5, 7, singularities and grids around airfoils and cross-sections of ship 12, 13] which investigated the conditions under which the sterns.