We model here the behavior of spin transition compounds, considering molecules arranged in a 2D hexagonal lattice and interacting via springs. The role of impurities in the clustering and nucleation phenomena is analyzed, as well as the manner in which the impurities affect the relaxation curves. The switching of the individual molecules is checked using a Monte Carlo procedure. When a molecule changes its state, it also modifies its volume, and the new equilibrium positions of all the molecules are calculated. As in previously reported experiments, it is found here that bigger impurities slow down the relaxation from the metastable high-spin state to the low-spin state, while smaller impurities act in an opposite way. It is shown that if the concentration of the impurities is higher than a certain threshold, then they act as a barrier, trammeling the fast evolution of domains developing from the edges.