The TOV equation appears as the relativistic counterpart of the classical condition for hydrostatic equilibrium. In the present work we aim at showing that a generalised TOV equation also characterises the equilibrium of models endowed with other symmetries besides spherical. We apply the dual null formalism to spacetimes with two dimensional spherical, planar and hyperbolic symmetries with a perfect fluid as the source. We also assume a Killing vector field orthogonal to the surfaces of symmetry, which gives us static solutions, in the timelike Killing field case, and homogeneous dynamical solutions in the case the Killing field is spacelike. In order to treat equally all the aforementioned cases, we discuss the definition of a quasi-local energy for the spacetimes with planar and hyperbolic foliations, since the Hawking-Hayward definition only applies to compact foliations. After this procedure, we are able to translate our geometrical formalism to the fluid dynamics language in a unified way, to find the generalized TOV equation, for the three cases when the solution is static, and to obtain the evolution equation, for the homogeneous spacetime cases. Remarkably, we show that the static solutions which are not spherically symmetric violate the weak energy condition (WEC). We have also shown that the counterpart of the TOV equation for the spatially homogeneous models is just the familiar equation \r{ho} + P = 0, defining a cosmological constant-type behaviour, both in the hyperbolic and spherical cases. This implies a violation of the strong energy condition in both cases, added to the above mentioned violation of the weak energy condition in the hyperbolic case. We illustrate our unified treatment obtaining analogs of Schwarzschild interior solution, for an incompressible fluid $\rho = \rho_0$ constant.

The Birkhoff theorem is a well-known result in general relativity and it is used in many applications. However, its most general version, due to Bona, is almost unknown and presented in a form less accessible to the relativist and cosmologist community. Moreover, many wield it mistakenly as a simple transposition of Newton's iron sphere theorem. In the present work, we propose a modern, dual null, presentation --- useful in many explorations, including black holes --- of the theorem that renders accessible most of the results of Bona's version. In addition, we discuss the fluid contents admissible for the application of the theorem, beyond a vacuum, and we demonstrate how the formalism greatly simplifies solving the dynamical equations and allows one to express the solution as a power expansion in $r$. We present a family of solutions that share the properties predicted by the Birkhoff theorem and discuss the existence of trapped and antitrapped regions. The formalism manifestly shows how the type of region --- trapped or untrapped --- determines the character of the Killing vector.

A class of nonstationary spacetimes is obtained by means of a conformal transformation of the Schwarzschild metric, where the conformal factor $a(t)$ is an arbitrary function of the time coordinate only. We investigate several situations including some where the final state is a central object with constant mass. The metric is such that there is an initial big-bang type singularity and the final state depends on the chosen conformal factor. The Misner-Sharp mass is computed and a localized central object may be identified. The trapping horizons, geodesic and causal structure of the resulting spacetimes are investigated in detail. When $a(t)$ asymptotes to a constant in a short enough time scale, the spacetime presents an event horizon and its analytical extension reveals black-hole or white-hole regions. On the other hand, when $a(t)$ is unbounded from above as in cosmological models, the spacetime presents no event horizons and may present null singularities in the future. The energy-momentum content and other properties of the respective spacetimes are also investigated.

Theories of gravity with a preferred foliation usually display arbitrarily fast signal propagation, changing the black hole definition. A new inescapable barrier, the universal horizon, has been defined and many static and spherically symmetric examples have been studied in the literature. Here, we translate the usual definition of the universal horizon in terms of an optical scalar built with the preferred flow defined by the preferred spacetime foliation. The new expression has the advantages of being of quasilocal nature and independent of specific spacetime symmetries in order to be well defined. Therefore, we propose it as a definition for general quasilocal universal horizons. Using the new formalism we show that there are no universal analog of cosmological horizons for FLRW models for any scale factor function, and we also state that quasilocal universal horizons are restricted to trapped regions of the spacetime. Using the evolution equation, we analyze the formation of universal horizons under a truncated Horava-Lifshitz theory, in spherical symmetry, showing the existence of regions in parameter space where the universal horizon formation cannot be smooth from the center, under some physically reasonable assumptions. We conclude with our view on the next steps for the understanding of black holes in nonrelativistic gravity theories.

Granular materials of different sizes are present on the surface of several atmosphere-less Solar System bodies. The phenomena related to granular materials have been studied in the framework of the discipline called Granular Physics; that has been studied experimentally in the laboratory and, in the last decades, by performing numerical simulations. The Discrete Element Method simulates the mechanical behavior of a media formed by a set of particles which interact through their contact points. The difficulty in reproducing vacuum and low-gravity environments makes numerical simulations the most promising technique in the study of granular media under these conditions. In this work, relevant processes in minor bodies of the Solar System are studied using the Discrete Element Method. Results of simulations of size segregation in low-gravity environments in the cases of the asteroids Eros and Itokawa are presented. The segregation of particles with different densities was analysed, in particular, the case of comet P/Hartley 2. The surface shaking in these different gravity environments could produce the ejection of particles from the surface at very low relative velocities. The shaking causing the above processes is due to: impacts, explosions like the release of energy by the liberation of internal stresses or the re accommodation of material. Simulations of the passage of impact-induced seismic waves through a granular medium were also performed. We present several applications of the Discrete Element Methods for the study of the physical evolution of agglomerates of rocks under low-gravity environments.