The backreaction of inhomogeneities describes the effect of inhomogeneous structure on average properties of the Universe. We investigate this approach by testing the consistency of cosmological $N$-body simulations as non-linear structure evolves. Using the Delaunay Tessellation Field Estimator (DTFE), we calculate the kinematical backreaction $Q$ from simulations on different scales in order to measure how much $N$-body simulations should be corrected for this effect. This is the first step towards creating fully relativistic and inhomogeneous $N$-body simulations. In this paper we compare the interpolation techniques available in DTFE and illustrate the statistical dependence of $Q$ as a function of length scale.

The 3-torus (T^3) FLRW model better fits the nearly zero large-scale auto-correlation of the WMAP CMB sky maps than the infinite flat model. The T^3 model's WMAP parameters imply approximately equal-redshift topological lensing at z \sim 6. We investigate observational strategies for rejecting the T^3 solution or providing candidate topologically lensed galaxy pairs. T^3 holonomies are applied to (i) existing z \sim 6 observations and (ii) simulated observations, creating multiply connected catalogues. Corresponding simply connected catalogues are generated. Each catalogue is analysed using a successive filter method and collecting matched quadruples. Quadruple statistics between the multiply and simply connected catalogues are compared. The expected rejection of the hypothesis, or detection of candidate topologically lensed galaxies, is possible at a significance of 5% for a pair of T^3 axis-centred northern and southern surveys if photometric redshift accuracy is \sigma(\zphot) < 0.01 for a pair of nearly complete 100 deg^2 surveys with a total of > 500 galaxies over 4.3 < z < 6.6, or for a pair of 196 deg^2 surveys with > 400 galaxies and \sigma(\zphot) < 0.02 over 4<z<7. Dropping the maximum time interval in a pair from \Delta t =1 Gyr/h to \Delta t =0.1 Gyr/h yields a requirement of \sigma(\zphot) < 0.005 or \sigma(\zphot) < 0.01, respectively. Millions of z \sim 6 galaxies will be observed over fields of these sizes during the coming decades, implying much stronger constraints. The question is not if the hypothesis will be rejected or confirmed, it is when.

One of the FLRW models that best fits the WMAP sky maps of the CMB is the Poincare dodecahedral space. The optimal fit of this model to WMAP data was recently found using an optimal cross-correlation method, but the systematic error in the estimate of the matched-circle angular radius \alpha, or equivalently, the (comoving) size of the Universe 2\rinj (twice the injectivity radius), might be much higher than the random error. In order to increase the falsifiability of the model, it would be useful to reduce the uncertainty in this estimate and to estimate the fraction of the sky where multiply imaged gravitationally bound objects should potentially be detectable. "Matched discs" are defined in order to describe a useful subset of multiply imaged objects. The cross-correlation method at \ltapprox 1 \hGpc is applied to WMAP 7-year data in order to improve the estimate of \alpha. The improved matched-circle radius estimate is \alpha = 23 \pm 1.4 deg, where the uncertainty represents systematic error dependent on the choices of galactic mask and all-sky map. This is equivalent to 2\rinj = 18.2\pm 0.5 \hGpc for matter density parameter Omega_m=0.28\pm 0.02. The lowest redshift of multiply imaged objects is z=106\pm18. Multiply imaged high overdensity (rare) peaks visible during 200>z>106 should be present in matched discs of radius 14.8\pm2.3 deg. The accuracy in the matched circle radius estimate is considerably improved by using the higher resolution signal. The predicted matched discs (over 200>z>106) project to about 20% of the full sky. Since any object located exactly in the discs would be multiply imaged at equal redshifts, evolutionary effects would be small for objects that are nearly located in the discs.