Two-dimensional elliptic grid solver using boundary grid control and curvature correction

The Thompson method of explicitly evaluating control function of the ellipt ic grid solver is employed to improve local grid quality around convex and concave boundaries. The curvature correction is principally added to the gr id lines around a boundary and is quickly attenuated as the grid points mov e inward. The algebraic advancing front method is employed to estimate two levels of orthogonal grids for evaluating the control functions at the boun dary. If the grid equation based on the Cauchy-Riemann relation with xi and eta as independent variables is employed, numerical examinations show that a lengthy trial-and-error procedure is required to get a satisfactory grid distribution. In contrast, the method based on the Cauchy-Riemann relation with x and y as independent variables, which benefits from the maximum pri nciple, not only removes the undesired grid clustering and diluting easily but also effectively improves overall grid smoothness and slightly enhances boundary grid orthogonality.