A hybrid differential‐algebraic method for three‐dimensional grid generation

An effective method is presented for generating three-dimensional boundary-fitted grids. The technique is based on a two-step procedure involving the solution of differential equations and algebraic interpolations. A coarse grid with an arbitrarily given point distribution at all the boundaries is first generated by solving the Poisson equations, with the forcing functions being automatically adjusted during the grid generation to achieve orthogonality at the boundaries. A fine grid with a desired concentration of grid lines is then obtained by fitting tricubic spline functions on the coarse grid. A primary feature of the method is that the mesh spacing can be quickly and explicitly changed while retaining the smoothness over the entire field and the boundary orthogonality. Applications to a variety of geometries demonstrate the performance of the method.