The simulation of the transition sequence of superheated Type I superconducting granules (SSG) in disordered suspensions when an external magnetic field is slowly increased from zero has been studied. Simulation takes into account diamagnetic interactions and the presence of surface defects. Results have been obtained for the transition sequence and surface fields distribution covering a wide range of densities. These results are compared with previous analytical perturbative theory, which provides qualitative information on transitions and surface magnetic fields during transitions, but with a range of validity apparently limited to extremely dilute samples. Simulations taking into account the complete diamagnetic interactions between spheres appear to be a promising tool in interpreting SSG experiments, in applications such as particle detectors, and in some fundamental calculations of Solid State Physics.