A procedure for tetrahedral boundary layer mesh generation

We develop a methodology for introducing regions of high anisotropy in existing isotropic unstructured grids in complex, curved, three-dimensional domains. The new procedures are here applied to the capturing of solution features in the proximity of model boundaries (e.g. boundary layers). Suitable voids are created in an existing grid in the regions of localization using a mesh motion algorithm that solves a fictitious elasticity problem. The voids are then filled with stacks of prisms that are subsequently tetrahedronized to yield a simplicial mesh. The mesh motion algorithm allows us to deal in a simple and effective manner with the problem of self-intersection of elements in concave regions of the model boundaries, and in the case of closely spaced model faces, avoiding the need for cross-over checks and complex grid correction procedures. The capabilities and performance of the proposed methodology are illustrated with the help of practical examples.