A finite volume method for numerical grid generation

A novel method to generate body‐fitted grids based on the direct solution for three scalar functions is derived. The solution for scalar variables ξ, η and ζ is obtained with a conventional finite volume method based on a physical space formulation. The grid is adapted or re‐zoned to eliminate the residual error between the current solution and the desired solution, by means of an implicit grid‐correction procedure. The scalar variables are re‐mapped and the process is reiterated until convergence is obtained. Calculations are performed for a variety of problems by assuming combined Dirichlet–Neumann and pure Dirichlet boundary conditions involving the use of transcendental control functions, as well as functions designed to effect grid control automatically on the basis of boundary values. The use of dimensional analysis to build stable exponential functions and other control functions is demonstrated. Automatic procedures are implemented: one based on a finite difference approximation to the Cristoffel terms assuming local‐boundary orthogonality, and another designed to procure boundary orthogonality. The performance of the new scheme is shown to be comparable with that of conventional inverse methods when calculations are performed on benchmark problems through the application of point‐by‐point and whole‐field solution schemes. Advantages and disadvantages of the present method are critically appraised. Copyright © 1999 National Research Council of Canada.