AbstractThis paper describes a fully automatic 3D anisotropic mesh generation method for domains of arbitrary shape. The spacing of the boundary mesh is computed by the analysis of the principal curvatures and directions of the boundary surfaces. The spacing in the domain is obtained by interpolation of the spacing at the boundaries on a suitably constructed background mesh. Examples which illustrate the performance of the proposed methodology are presented. Copyright © 2005 John Wiley & Sons, Ltd.