Abstract: In order to solve the problem of multi-domain adjacency in the modeling of complex orebodies, we implemented an adaptive meshing method for complex orebodies based on the Delaunay refinement method, which could generate multi-domain surface models and volume models for the finite element simulation in the later stage. The method relies on the constraint Delaunay triangulation to approximate the domain and the surface, and the Delaunay refinement to ensure the approximate accuracy of the domain and the grid quality of the cell. On the basis of constraint Delaunay triangulation, the multi-domain polyhedron method of prediction was constructed by separating the intermediate domain or fitting the intermediate domain by distance field. The method of adaptive meshing of multi-domain was implemented. Experimental tests of Beiminghe Iron Mine show that the method can generate a seamless and continuous multi-domain model with high-quality tetrahedral and triangular meshes by mesh optimization. And it can avoid many degraded and singular triangles by contour-matching. It is significant to apply this method in the reverse engineering, finite element analysis, scientific computing visualization and other fields.