The expanding complex pattern of filaments, walls and voids build the evolving cosmic web with material flowing from underdense onto high density regions. Here we explore the dynamical behaviour of voids and galaxies in void shells relative to neighboring overdense superstructures, using the Millenium Simulation and the main galaxy catalogue in Sloan Digital Sky Survey data. We define a correlation measure to estimate the tendency of voids to be located at a given distance from a superstructure. We find voids-in-clouds (S-types) preferentially located closer to superstructures than voids-in-voids (R-types) although we obtain that voids within $\sim40~\mathrm{Mpc}\,\mathrm{h}^{-1}$ of superstructures are infalling in a similar fashion independently of void type. Galaxies residing in void shells show infall towards the closest superstructure, along with the void global motion, with a differential velocity component depending on their relative position in the shell with respect to the direction to the superstructure. This effect is produced by void expansion and therefore is stronger for R-types. We also find that galaxies in void shells facing the superstrucure flow towards the overdensities faster than galaxies elsewere at the same relative distance to the superstructure. The results obtained for the simulation are also reproduced for the SDSS data with a linearized velocity field implementation.

We perform a statistical study of the global motion of cosmic voids using both a numerical simulation and observational data. We analyse their relation to large-scale mass flows and the physical effects that drive those motions. We analyse the bulk motions of voids, defined by the mean velocity of haloes in the surrounding shells in the numerical simulation, and by galaxies in the Sloan Digital Sky Survey Data Release 7. We find void mean bulk velocities close to 400 km/s, comparable to those of haloes (~ 500-600 km/s), depending on void size and the large-scale environment. Statistically, small voids move faster than large ones, and voids in relatively higher density environments have higher bulk velocities than those placed in large underdense regions. Also, we analyze the mean mass density around voids finding, as expected, large-scale overdensities (underdensities) along (opposite to) the void motion direction, suggesting that void motions respond to a pull-push mechanism. This contrasts with massive cluster motions who are mainly governed by the pull of the large-scale overdense regions. Our analysis of void pairwise velocities shows how their relative motions are generated by large-scale density fluctuations. In agreement with linear theory, voids embedded in low (high) density regions mutually recede (attract) each other, providing the general mechanism to understand the bimodal behavior of void motions. In order to compare the theoretical results and the observations we have inferred void motions in the SDSS using linear theory, finding that the estimated observational void motions are in statisticalagreement with the results of the simulation. Regarding large-scale flows, our results suggest a scenario of galaxies and galaxy systems flowing away from void centers with the additional, and morerelevant, contribution of the void bulk motion to the total velocity.

We compute the bulk motions of cosmic voids, using a $\Lambda$CDM numerical simulation considering the mean velocities of the dark matter inside the void itself and that of the haloes in the surrounding shell. We find coincident values of these two measures in the range $\sim$ 300-400 km/s, not far from the expected mean peculiar velocities of groups and galaxy clusters. When analysing the distribution of the pairwise relative velocities of voids, we find a remarkable bimodal behaviour consistent with an excess of both systematically approaching and receding voids. We determine that the origin of this bimodality resides in the void large scale environment, since once voids are classified into void-in-void (R-type) or void-in-cloud (S-type), R-types are found mutually receding away, while S-types approach each other. The magnitude of these systematic relative velocities account for more than 100 km/s, reaching large coherence lengths of up to 200 h$^{-1}$ Mpc . We have used samples of voids from the Sloan Digital Sky Survey Data Release 7 (SDSS-DR7) and the peculiar velocity field inferred from linear theory, finding fully consistent results with the simulation predictions. Thus, their relative motion suggests a scenario of a sparkling Universe, with approaching and receding voids according to their local environment.