With a Yang-Mills flux-tube initial state and a high-resolution (3+1)D particle-in-cell relativistic (PICR) hydrodynamics simulation, we calculate the polarization for different energies. The origination of polarization in high energy collisions is discussed, and we find linear impact parameter dependence of the global polarization. Furthermore, the global polarization in our model decreases very quickly in the low energy domain, and the decline curve fits well the recent results of Beam Energy Scan (BES) program launched by the STAR Collaboration at the Relativistic Heavy Ion Collider (RHIC). The time evolution of polarization is also discussed.