Using elliptic equation, the best control function is Thomas and midlecoff, and modified Hilgenstock, some author used the anisotropic transfinite interpolation for forcing function for increaing the decay rate of the function to improve convergence and stability, but this forcing function does not fitted for coner mesh distribution because it is very difficult to distribute the grid line angle. up to now, I do not find a method to distribute the angle, So I need some guy's help. I am very appreciate that some guy give me a good suggestion

(1). The best place to start is to take a look at the book free on Internet by Prof. Joe Thompson, "Numerical Grid Generation", if you have not done so. (2). You can also try to e-mail the researchers there to discuss your problems. Numerical grid generation is hard to control,so, if possible, try some simple algebraic methods.

To find a "method" to make coner mesh distribution is very difficult. It differs of the case on which you work but like says Jhon you must contact some one who is a specialist.

If you can email me the boundary points for your corner grid (perhaps in PLOT3D format) I would be happy to try creating your grid in Gridgen.

Gridgen uses a combined method that includes Thomas & Middlecoff with either Hilgenstock-White or Steger-Sorenson to obtain orthogonality at the boundary.