We investigate Kantowski-Sachs models in Einstein-{\ae}ther theory with a perfect fluid source using the singularity analysis to prove the integrability of the field equations and dynamical system tools to study the evolution. We find an inflationary source at early times, and an inflationary sink at late times, for a wide region in the parameter space. The results by A. A. Coley, G. Leon, P. Sandin and J. Latta (JCAP 12, 010, 2015), are then re-obtained as particular cases. Additionally, we select other values for the non-GR parameters which are consistent with current constraints, getting a very rich phenomenology. In particular, we find solutions with infinite shear, zero curvature, and infinite matter energy density in comparison with the Hubble scalar. We also have stiff-like future attractors, anisotropic late-time attractors, or both, in some special cases. Such results are developed analytically, and then verified by numerics. Finally, the physical interpretation of the new critical points is discussed.

We investigate spherically symmetric cosmological models in Einstein-aether theory with a tilted (non-comoving) perfect fluid source. We use a 1+3 frame formalism and adopt the comoving aether gauge to derive the evolution equations, which form a well-posed system of first order partial differential equations in two variables. We then introduce normalized variables. The formalism is particularly well-suited for numerical computations and the study of the qualitative properties of the models, which are also solutions of Horava gravity. We study the local stability of the equilibrium points of the resulting dynamical system corresponding to physically realistic inhomogeneous cosmological models and astrophysical objects with values for the parameters which are consistent with current constraints. In particular, we consider dust models in ($\beta-$) normalized variables and derive a reduced (closed) evolution system and we obtain the general evolution equations for the spatially homogeneous Kantowski-Sachs models using appropriate bounded normalized variables. We then analyse these models, with special emphasis on the future asymptotic behaviour for different values of the parameters. Finally, we investigate static models for a mixture of a (necessarily non-tilted) perfect fluid with a barotropic equations of state and a scalar field.

We use cosmological perturbation theory to study the backreaction effects of a self-consistent and well-defined cosmological averaging on the dynamics and the evolution of the Universe. Working with a perturbed Friedman-Lemaitre-Robertson-Walker Einstein-de Sitter cosmological solution in a comoving volume-preserving gauge, we compute the expressions for the expansion scalar and deceleration parameter to second order, which we use to characterize the backreaction. We find that the fractional shift in the Hubble parameter with respect to the input background cosmological model is Delta~10^{-5}, which leads to an effective energy density of the order of a few times 10^{-5}. In addition, we find that an appropriate measure of the fractional shift in the deceleration parameter Q is very large.

The cosmological backreaction from perturbations is clearly gauge-dependent, and obviously depends on the choice of averaged Hubble rate. We consider two common choices of Hubble rate and advocate the use of comoving volume-preserving gauges. We highlight two examples valid to an appropriate order in perturbation theory, uniform curvature gauge, which is as close to volume-preserving as possible, and a spatially-traceless uniform cold dark matter gauge which preserves the volume to linear order. We demonstrate the strong gauge- and frame-dependences in averaging. In traceless uniform CDM gauge the backreaction exhibits a strong ultra-violet divergence and can be tuned to an arbitrary magnitude with an appropriate choice of smoothing scale. In uniform curvature gauge we find that for a choice of Hubble rate locked to the spatial surface the backreaction vanishes identically, while for a Hubble rate defined from a fluid's expansion scalar the effective energy density at the current epoch in an Einstein-de Sitter universe is Omega_eff~5e-4, slightly bigger than but in broad agreement with previous results in conformal Newtonian gauge.

The Observed Perlmutter-Riess Acceleration (OPRA) implies that the expansion of the Universe is currently increasing and is motivation for the so-called phantom energy models. We consider the dynamics of phantom scalar field models. An important physical time constraint, which can be used to rule out many cosmological models, is obtained from the condition that all forms of energy density, including the field causing OPRA (e.g., the phantom field), must be non-negligible for an extended period, which is conservatively estimated to be of the order of a few $Gyr$. We find that this physical time constraint cannot be satisfied in conventional phantom cosmological models.