Three numerical schemes and one approximate model are developed to compute the overall rate constants for diffusion and incorporation of small entities in clusters of spheres. These include the Brownian dynamic simulation, multipole expansion, boundary co