A numerical study of cluster center formation in neutron-irradiated silicon