Beam arrays have great application value in free-space optical communication. In this paper, the light intensity evolution and the on-axis scintillation index of radial Gaussian vortex beam arrays propagating through atmospheric turbulence are analyzed using multi-phase screen simulation. The effect of initial beam parameters on the on-axis scintillation index of radial Gaussian vortex beam arrays is studied, and the on-axis scintillation index of radial Gaussian vortex beam arrays is compared with that of a single Gaussian vortex beam. The results indicate that in the weak fluctuation regime and when the rytov index is less than 0.5, the on-axis scintillation index of Gaussian vortex beams remains within a numerical range of less than 1, while the on-axis scintillation index of radial Gaussian vortex beam arrays is around 1. In the medium fluctuation regime, the on-axis scintillation index of the radial Gaussian vortex beam arrays is smaller than that of a single Gaussian vortex beam. The on-axis scintillation index of radial Gaussian vortex beam arrays decreases with the decrease of orbital angular momentum and the increase of radial array radius. The research results hold theoretical significance and application value for vortex optical communication in turbulent atmospheric environments.