We examine the fully relativistic evolution of cosmic voids constituted by baryons and cold dark matter (CDM), represented by two non-comoving dust sources in a $\Lambda$CDM background. For this purpose, we consider numerical solutions of Einstein's field equations in a fluid-flow representation adapted to spherical symmetry and multiple components. We present a simple example that explores the frame-dependence of the local expansion and the Hubble flow for this mixture of two dusts, revealing that the relative velocity between the sources yields a significantly different evolution in comparison with that of the two sources in a common 4-velocity (which reduces to a Lemaitre-Tolman-Bondi model). In particular, significant modifications arise for the density contrast depth and void size, as well as in the amplitude of the surrounding over-densities. We show that an adequate model of a frame-dependent evolution that incorporates initial conditions from peculiar velocities and large-scale density contrast observations may contribute to understand the discrepancy between the local value of $H_0$ and that inferred from the CMB.

We examine the gravitational collapse and black hole formation of multiple non--spherical configurations constructed from Szekeres dust models with positive spatial curvature that smoothly match to a Schwarzschild exterior. These configurations are made of an almost spherical central core region surrounded by a network of "pancake-like" overdensities and voids with spatial positions prescribed through standard initial conditions. We show that a full collapse into a focusing singularity, without shell crossings appearing before the formation of an apparent horizon, is not possible unless the full configuration becomes exactly or almost spherical. Seeking for black hole formation, we demand that shell crossings are covered by the apparent horizon. This requires very special fine-tuned initial conditions that impose very strong and unrealistic constraints on the total black hole mass and full collapse time. As a consequence, non-spherical non-rotating dust sources cannot furnish even minimally realistic toy models of black hole formation at astrophysical scales: demanding realistic collapse time scales yields huge unrealistic black hole masses, while simulations of typical astrophysical black hole masses collapse in unrealistically small times. We note, however, that the resulting time--mass constraint is compatible with early Universe models of primordial black hole formation, suitable in early dust-like environments. Finally, we argue that the shell crossings appearing when non-spherical dust structures collapse are an indicator that such structures do not form galactic mass black holes but virialise into stable stationary objects.

We study the differences and equivalences between the non-perturbative description of the evolution of cosmic structure furnished by the Szekeres dust models (a non-spherical exact solution of Einstein's equations) and the dynamics of Cosmological Perturbation Theory (CPT) for dust sources in a $\Lambda$CDM background. We show how the dynamics of Szekeres models can be described by evolution equations given in terms of "exact fluctuations" that identically reduce (at all orders) to evolution equations of CPT in the comoving isochronous gauge. We explicitly show how Szekeres linearised exact fluctuations are specific (deterministic) realisations of standard linear perturbations of CPT given as random fields but, as opposed to the latter perturbations, they can be evolved exactly into the full non-linear regime. We prove two important results: (i) the conservation of the curvature perturbation (at all scales) also holds for the appropriate approximation of the exact Szekeres fluctuations in a $\Lambda$CDM background, and (ii) the different collapse morphologies of Szekeres models yields, at nonlinear order, different functional forms for the growth factor that follows from the study of redshift space distortions. The metric based potentials used in linear CPT are computed in terms of the parameters of the linearised Szekeres models, thus allowing us to relate our results to linear CPT results in other gauges. We believe that these results provide a solid starting stage to examine the role of non-perturbative General Relativity in current cosmological research.

We show that the full dynamical freedom of the well known Szekeres models allows for the description of elaborated 3--dimensional networks of cold dark matter structures (over--densities and/or density voids) undergoing "pancake" collapse. By reducing Einstein's field equations to a set of evolution equations, which themselves reduce in the linear limit to evolution equations for linear perturbations, we determine the dynamics of such structures, with the spatial comoving location of each structure uniquely specified by standard early Universe initial conditions. By means of a representative example we examine in detail the density contrast, the Hubble flow and peculiar velocities of structures that evolved, from linear initial data at the last scattering surface, to fully non--linear 10--20 Mpc. scale configurations today. To motivate further research, we provide a qualitative discussion on the connection of Szekeres models with linear perturbations and the pancake collapse of the Zeldovich approximation. This type of structure modelling provides a coarse grained -- but fully relativistic non--linear and non--perturbative -- description of evolving large scale cosmic structures before their virialisation, and as such it has an enormous potential for applications in cosmological research.

We undertake a hydrodynamical study of a mixture of tightly coupled primordial radiation, neutrinos, baryons, electrons and positrons, together with a gas of already decoupled dark matter WIMPS and an already existing "frozen" magnetic field in the infinite conductivity regime. Considering this cosmic fluid as the source of a homogeneous but anisotropic Bianchi I model, we describe its interaction with the magnetic field by means of suitable equations of state that are appropriate for the particle species of the mixture between the end of the leptonic era and the beginning of the radiation-dominated epoch. Fulfilment of observational bounds on the magnetic field intensity yields a "near FLRW" (but strictly non-perturbative) evolution of the geometric, kinematic and thermodynamical variables. This evolution is roughly comparable to the weak field approximation in linear perturbations on a spatially flat FLRW background of sources in which the frozen magnetic fields are coherent over very large supra-horizon scales. Our approach and results may provide interesting guidelines in potential situations in which non-perturbative methods are required to study the interaction between magnetic fields and the cosmic fluid.

We examine the spatial extrema (local maxima, minima and saddle points) of the covariant scalars (density, Hubble expansion, spatial curvature and eigenvalues of the shear and electric Weyl tensors) of the quasi-spherical Szekeres dust models. Sufficient conditions are obtained for the existence of distributions of multiple extrema in spatial comoving locations that can be prescribed through initial conditions. These distributions evolve without shell crossing singularities at least for ever expanding models (with or without cosmological constant) in the full evolution range where the models are valid. By considering the local maxima and minima of the density, our results allow for setting up elaborated networks of "pancake" shaped evolving cold dark matter over-densities and density voids whose spatial distribution and amplitudes can be controlled from initial data compatible with standard early Universe initial conditions. We believe that these results have an enormous range of potential application by providing a fully relativistic non-perturbative coarse grained modelling of cosmic structure at all scales.