Growth of the structure in the Universe manifest as accretion flows of galaxies onto groups and clusters. Thus, the present day properties of groups and their member galaxies are influenced by the characteristics of this continuous infall pattern. Several works both theoretical, in numerical simulations, and in observations, study this process and provide useful steps for a better understanding of galaxy systems and their evolution. We aim at exploring the streaming flow of galaxies onto groups using observational peculiar velocity data. The effects of distance uncertainties are also analyzed as well as the relation between the infall pattern and group and environment properties.This work deals with analysis of peculiar velocity data and their projection on the direction to group centers, to determine the mean galaxy infall flow. We applied this analysis to the galaxies and groups extracted from the Cosmicflows-3 catalog. We also use mock catalogs derived from numerical simulations to explore the effects of distance uncertainties on the derivation of the galaxy velocity flow onto groups. We determine the infalling velocity field onto galaxy groups with cz < 0.033 using peculiar velocity data. We measure the mean infall velocity onto group samples of different mass range, and also explore the impact of the environment where the group reside. Well beyond the group virial radius, the surrounding large-scale galaxy overdensity may impose additional infalling streaming amplitudes in the range 200 to 400 km s$^{-1}$. Also, we find that groups in samples with a well controlled galaxy density environment show an increasing infalling velocity amplitude with group mass, consistent with the predictions of the linear model. These results from observational data are in excellent agreement with those derived from the mock catalogs.

We perform statistical analyses to study the infall of galaxies onto groups and clusters in the nearby Universe. The study is based on the UZC and SSRS2 group catalogs and peculiar velocity samples. We find a clear signature of infall of galaxies onto groups over a wide range of scales 5 h^{-1} Mpc<r<30 h^{-1} Mpc, with an infall amplitude on the order of a few hundred kilometers per second. We obtain a significant increase in the infall amplitude with group virial mass (M_{V}) and luminosity of group member galaxies (L_{g}). Groups with M_{V}<10^{13} M_{\odot} show infall velocities V_{infall} \simeq 150 km s^{-1} whereas for M_{V}>10^{13} M_{\odot} a larger infall is observed, V_{infall} \simeq 200 km s^{-1}. Similarly, we find that galaxies surrounding groups with L_{g}<10^{15} L_{\odot} have V_{infall} \simeq 100 km s^{-1}, whereas for L_{g}>10^{15} L_{\odot} groups, the amplitude of the galaxy infall can be as large as V_{infall} \simeq 250 km s^{-1}. The observational results are compared with the results obtained from mock group and galaxy samples constructed from numerical simulations, which include galaxy formation through semianalytical models. We obtain a general agreement between the results from the mock catalogs and the observations. The infall of galaxies onto groups is suitably reproduced in the simulations and, as in the observations, larger virial mass and luminosity groups exhibit the largest galaxy infall amplitudes. We derive estimates of the integrated mass overdensities associated with groups by applying linear theory to the infall velocities after correcting for the effects of distance uncertainties obtained using the mock catalogs. The resulting overdensities are consistent with a power law with \delta \sim 1 at r \sim 10 h^{-1}Mpc.