We study the assembly of globular clusters (GCs) in 9 galaxy clusters using the cosmological simulation Illustris. GCs are tagged to individual galaxies at their infall time. The tidal removal of GCs from their galaxies and the distribution of the GCs within the cluster is later followed self-consistently by the simulation. The method relies on the simple assumption of a single power-law relation between halo mass (M_vir) and mass in GCs (M_GC) as found in observations. We find that the GCs specific frequency $S_N$ as a function of V-band magnitude naturally reproduces the observed "U"-shape, due to the combination of a power law M_GC-M_vir relation and the non-linear M_*-M_vir relation from the simulation. Additional scatter in the $S_N$ values are traced back to galaxies with early infall times due to the evolution in the M_*-M_vir relation with redshift. GCs that have been tidally removed from their galaxies form today the intra-cluster component from which about ~60% were brought in by galaxies that orbit today within the cluster potential. The remaining "orphan" GCs are contributed by satellite galaxies with a wide range of stellar masses that are fully tidally disrupted at z=0. This intra-cluster component is a good dynamical tracer of the dark matter potential. As a consequence of the accreted nature of most intra-cluster GCs, their orbits are fairly radial with a predicted orbital anisotropy \beta >= 0.5. However, local tangential motions may appear as a consequence of localized substructure, providing a possible interpretation to the \beta<0 values suggested in observations of M87.

Using a cosmological dark matter simulation of a galaxy-cluster halo, we follow the temporal evolution of its globular cluster population. To mimic the red and blue globular cluster populations, we select at high redshift $(z\sim 1)$ two sets of particles from individual galactic halos constrained by the fact that, at redshift $z=0$, they have density profiles similar to observed ones. At redshift $z=0$, approximately 60\% of our selected globular clusters were removed from their original halos building up the intra-cluster globular cluster population, while the remaining 40\% are still gravitationally bound to their original galactic halos. Since the blue population is more extended than the red one, the intra-cluster globular cluster population is dominated by blue globular clusters, with a relative fraction that grows from 60\% at redshift $z=0$ up to 83\% for redshift $z\sim 2$. In agreement with observational results for the Virgo galaxy cluster, the blue intra-cluster globular cluster population is more spatially extended than the red one, pointing to a tidally disrupted origin.