Backreaction in Numerical Relativity: Averaging on Newtonian gauge-like hypersurfaces in Einstein Toolkit cosmological simulations

We introduce a spatial averaging scheme and use it to study the evolution of spatial averages in large-scale simulations of cosmological structure formation performed with the Einstein Toolkit. The averages are performed on the spatial hypersurfaces of the simulation setup which, at least initially, represent the hypersurfaces of statistical homogeneity and isotropy. We find only negligible cosmic backreaction on these hypersurfaces even on very small scales, but find significant curvature fluctuations of up to $10\%$ in $\Omega_R$ for sub-volumes with radius $\sim200$ Mpc and even larger fluctuations in smaller sub-volumes. In addition, we quantify fluid flow in and out of these sub-volumes. We find this to be significant, up to a $5\%$ change in the density between redshift $z=1$ and $z=0$ of a single sphere of radius $\sim200$ Mpc (and larger for smaller spheres). We suggest this may be important for studies basing averages on volumes co-moving with the simulation hypersurfaces.