This paper presents a simulation method for evaluating the accuracy of assembly of two mating freeform surfaces based on the spatial CAD model of the gap between the surfaces. CAD models of actual surfaces are the basis for determining the gap model. Digital measurement data representing form profiles are determined in measurements on a coordinate measuring machine according to a regular grid of points, and data representing waviness and surface roughness on a profilometer. The SGP components separated in the measurement process are summed up in the CAD modeling process. In the next step, a closure simulation is carried out by docking the inverted complex model of one surface to the surface model of the other surface. In this way, a spatial complex CAD model of the gap between the surfaces is obtained, mapping the effect of all the components of the irregularity on the tightness of the assembly. Based on the obtained model, it is possible to determine its maximum dimension (estimate the local width of the gap between actual surfaces) and the average of local dimensions, but also to illustrate the detailed distribution of local dimensions of this model in the form of a color map. The effectiveness of the proposed method was verified experimentally for the closing surfaces of the elements of the injection mold (cavity and core).