General method for producing a boundary-fitted orthogonal curvilinear mesh

A general method is described for computing an orthogonal mesh fitted to a two-dimensional physical domain with arbitrary closed boundary. The method allows optimum control of mesh spacing through the introduction of arbitrary (with weak constraints) ‘packing’ functions into the elliptic governing equations. Two particular aspects are addressed: first, the condition on a scaling factor which normalizes the mesh aspect ratio; second, the condition for avoiding run-out of the mesh beyond the boundaries of the physical domain.Conversion of the equations to finite difference form and appropriate iterative techniques are discussed. Finally applications of the method in the context of flow across a bundle of rods are presented.