We discuss the construction of cosmological models within the framework of Macroscopic Gravity (MG), which is a theory that models the effects of averaging the geometry of space-time on large scales. We find new exact spatially homogeneous and isotropic FLRW solutions to the MG field equations, and investigate large-scale perturbations around them. We find that any inhomogeneous perturbations to the averaged geometry are severely restricted, but that possible anisotropies in the correlation tensor can have dramatic consequences for the measurement of distances. These calculations are a first step within the MG approach toward developing averaged cosmological models to a point where they can be used to interpret real cosmological data, and hence to provide a working alternative to the "concordance" LCDM model.

Schwarzschild's solution to the Einstein Field Equations was one of the first and most important solutions that lead to the understanding and important experimental tests of Einstein's theory of General Relativity. However, Schwarzschild's solution is essentially based on an ideal theory of gravitation, where all inhomogeneities are ignored. Therefore, any generalization of the Schwarzschild solution should take into account the effects of small perturbations that may be present in the gravitational field. The theory of Macroscopic Gravity characterizes the effects of the inhomogeneities through a non-perturbative and covariant averaging procedure. With similar assumptions on the geometry and matter content, a solution to the averaged field equations as dictated by Macroscopic Gravity are derived. The resulting solution provides a possible explanation for the flattening of galactic rotation curves, illustrating that Dark Matter is not real but may only be the result of averaging inhomogeneities in a spherically symmetric background.

We study tilted perfect fluid cosmological models with a constant equation of state parameter in spatially homogeneous models of Bianchi type VI$_{-1/9}$ using dynamical systems methods and numerical simulations. We study models with and without vorticity, with an emphasis on their future asymptotic evolution. We show that for models with vorticity there exists, in a small region of parameter space, a closed curve acting as the attractor.

We study tilted perfect fluid cosmological models with a constant equation of state parameter in spatially homogeneous models of Bianchi type VI_h using dynamical systems methods and numerical experimentation, with an emphasis on their future asymptotic evolution. We determine all of the equilibrium points of the type VI_h state space (which correspond to exact self-similar solutions of the Einstein equations, some of which are new), and their stability is investigated. We find that there are vacuum plane-wave solutions that act as future attractors. In the parameter space, a `loophole' is shown to exist in which there are no stable equilibrium points. We then show that a Hopf-bifurcation can occur resulting in a stable closed orbit (which we refer to as the Mussel attractor) corresponding to points both inside the loophole and points just outside the loophole; in the former case the closed curves act as late-time attractors while in the latter case these attracting curves will co-exist with attracting equilibrium points. In the special Bianchi type III case, centre manifold theory is required to determine the future attractors. Comprehensive numerical experiments are carried out to complement and confirm the analytical results presented. We note that the Bianchi type VI_h case is of particular interest in that it contains many different subcases which exhibit many of the different possible future asymptotic behaviours of Bianchi cosmological models.

We use expansion-normalised variables to investigate the Bianchi type VII$_0$ model with a tilted $\gamma$-law perfect fluid. We emphasize the late-time asymptotic dynamical behaviour of the models and determine their asymptotic states. Unlike the other Bianchi models of solvable type, the type VII$_0$ state space is unbounded. Consequently we show that, for a general non-inflationary perfect fluid, one of the curvature variables diverges at late times, which implies that the type VII$_0$ model is not asymptotically self-similar to the future. Regarding the tilt velocity, we show that for fluids with $\gamma<4/3$ (which includes the important case of dust, $\gamma=1$) the tilt velocity tends to zero at late times, while for a radiation fluid, $\gamma=4/3$, the fluid is tilted and its vorticity is dynamically significant at late times. For fluids stiffer than radiation ($\gamma>4/3$), the future asymptotic state is an extremely tilted spacetime with vorticity.

Using dynamical systems theory and a detailed numerical analysis, the late-time behaviour of tilting perfect fluid Bianchi models of types IV and VII$_h$ are investigated. In particular, vacuum plane-wave spacetimes are studied and the important result that the only future attracting equilibrium points for non-inflationary fluids are the plane-wave solutions in Bianchi type VII$_h$ models is discussed. A tiny region of parameter space (the loophole) in the Bianchi type IV model is shown to contain a closed orbit which is found to act as an attractor (the Mussel attractor). From an extensive numerical analysis it is found that at late times the normalised energy-density tends to zero and the normalised variables 'freeze' into their asymptotic values. A detailed numerical analysis of the type VII$_h$ models then shows that there is an open set of parameter space in which solution curves approach a compact surface that is topologically a torus.

We investigate the dynamical properties of a class of spatially homogeneous and isotropic cosmological models containing a barotropic perfect fluid and multiple scalar fields with independent exponential potentials. We show that the assisted inflationary scaling solution is the global late-time attractor for the parameter values for which the model is inflationary, even when curvature and barotropic matter are included. For all other parameter values the multi-field curvature scaling solution is the global late-time attractor (in these solutions asymptotically the curvature is not dynamically negligible). Consequently, we find that in general all of the scalar fields in multi-field models with exponential potentials are non-negligible in late-time behaviour, contrary to what is commonly believed. The early-time and intermediate behaviour of the models is also studied. In particular, n-scalar field models are investigated and the structure of the saddle equilibrium points corresponding to inflationary m-field scaling solutions and non-inflationary m-field matter scaling solutions are also studied (where m<n), leading to interesting transient dynamical behaviour with new physical scenarios of potential importance.