An unstructured mesh generation method for three dimensional region is described. The boundary surfaces are triangulated by a two dimensional anisotropic triangulation method. Then the three dimensional meshes are generated by means of a constrained Delaunay algorithm with automatic field point creation. Methods for mesh quality optimization are also discussed. Several examples are given to show the applications of the method.