A number of codes for general-relativistic simulations of cosmological structure formation have been developed in recent years. Here we demonstrate that a sample of these codes produce consistent results beyond the Newtonian regime. We simulate solutions to Einstein's equations dominated by gravitomagnetism -- a vector-type gravitational field that doesn't exist in Newtonian gravity and produces frame-dragging, the leading-order post-Newtonian effect. We calculate the coordinate-invariant effect on intersecting null geodesics by performing ray tracing in each independent code. With this observable quantity, we assess and compare each code's ability to compute relativistic effects.

With the era of precision cosmology upon us, and upcoming surveys expected to further improve the precision of our observations below the percent level, ensuring the accuracy of our theoretical cosmological model is of the utmost importance. Current tensions between our observations and predictions from the standard cosmological model have sparked curiosity in extending the model to include new physics. Although, some suggestions include simply accounting for aspects of our Universe that are ignored in the standard model. One example acknowledges the fact that our Universe contains significant density contrasts on small scales; in the form of galaxies, galaxy clusters, filaments, and voids. This small-scale structure is smoothed out in the standard model, by assuming large-scale homogeneity of the matter distribution, which could have a measurable effect due to the nonlinearity of Einstein's equations. This backreaction of small-scale structures on the large-scale dynamics has been suggested to explain the measured accelerating expansion rate of the Universe. Current standard cosmological simulations ignore the effects of General Relativity by assuming purely Newtonian dynamics. In this thesis, we take the first steps towards quantifying the backreaction of small-scale structures by performing cosmological simulations that solve Einstein's equations directly. Simulations like these will allow us to quantify potentially important effects on our observations that could become measurable as the precision of these observations increases into the future.

We perform large-scale cosmological simulations that solve Einstein's equations directly via numerical relativity. Starting with initial conditions sampled from the cosmic microwave background, we track the emergence of a cosmic web without the need for a background cosmology. We measure the backreaction of large-scale structure on the evolution of averaged quantities in a matter-dominated universe. Although our results are preliminary, we find the global backreaction energy density is of order $10^{-8}$ compared to the energy density of matter in our simulations, and is thus unlikely to explain accelerating expansion under our assumptions. Sampling scales above the homogeneity scale of the Universe ($100-180\,h^{-1}$Mpc), in our chosen gauge, we find $2-3\%$ variations in local spatial curvature.

In a fully inhomogeneous, anisotropic cosmological simulation performed by solving Einstein's equations with numerical relativity, we find a local measurement of the effective Hubble parameter differs by less than 1\% compared to the global value. This variance is consistent with predictions from Newtonian gravity. We analyse the averaged local expansion rate on scales comparable to Type 1a supernova surveys, and find that local variance cannot resolve the tension between the \citet{riess2018b} and \citet{planck2018a} measurements.

We perform three-dimensional numerical relativity simulations of homogeneous and inhomogeneous expanding spacetimes, with a view towards quantifying non-linear effects from cosmological inhomogeneities. We demonstrate fourth-order convergence with errors less than one part in 10^6 in evolving a flat, dust Friedmann-Lemaitre-Roberston-Walker (FLRW) spacetime using the Einstein Toolkit within the Cactus framework. We also demonstrate agreement to within one part in 10^3 between the numerical relativity solution and the linear solution for density, velocity and metric perturbations in the Hubble flow over a factor of ~350 change in scale factor (redshift). We simulate the growth of linear perturbations into the non-linear regime, where effects such as gravitational slip and tensor perturbations appear. We therefore show that numerical relativity is a viable tool for investigating nonlinear effects in cosmology.

Using photometry and high resolution spectroscopy we investigate for the first time the physical connection between the open clusters NGC 5617 and Trumpler 22. Based on new CCD photometry we report their spatial proximity and common age of ~70 Myr. Based on high resolution spectra collected using the HERMES and UCLES spectrographs on the Anglo-Australian telescope, we present radial velocities and abundances for Fe, Na, Mg, Al, Si, Ca and Ni. The measured radial velocities are -38.63 +/-2.25 km/s for NGC 5617 and -38.46 +/-2.08 km/s for Trumpler 22. The mean metallicity of NGC 5617 was found to be [Fe/H] =-0.18 +/-0.02 and for Trumpler 22 was found to be [Fe/H] = -0.17 +/-0.04. The two clusters share similar abundances across the other elements, indicative of a common chemical enrichment history of these clusters. Together with common motions and ages we confirm that NGC 5617 and Trumpler 22 are a primordial binary cluster pair in the Milky Way.