Next Article in Journal
Rotational Diode: Clockwise/Counterclockwise Asymmetry in Conducting and Mechanical Properties of Rotating (semi)Conductors
Previous Article in Journal
Temporal Behavior of Local Characteristics in Complex Networks with Preferential Attachment-Based Growth
Previous Article in Special Issue
Periodic Solutions of Nonlinear Relative Motion Satellites
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamics of Hyperbolically Symmetric Fluids

1
Instituto Universitario de Física Fundamental y Matemáticas, Universidad de Salamanca, 37007 Salamanca, Spain
2
Escuela de Física, Facultad de Ciencias, Universidad Central de Venezuela, Caracas 1050, Venezuela
3
Departamento de Matemática Aplicada and Instituto Universitario de Física Fundamental y Matemáticas, Universidad de Salamanca, 37007 Salamanca, Spain
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Symmetry 2021, 13(9), 1568; https://doi.org/10.3390/sym13091568
Submission received: 18 May 2021 / Revised: 9 July 2021 / Accepted: 2 August 2021 / Published: 26 August 2021
(This article belongs to the Special Issue Recent Advance in Astrophysics)

Abstract

:
We study the general properties of dissipative fluid distributions endowed with hyperbolical symmetry. Their physical properties are analyzed in detail. It is shown that the energy density is necessarily negative, and the central region cannot be attained by any fluid element. We describe this inner region by a vacuum cavity around the center. By assuming a causal transport equation some interesting thermodynamical properties of these fluids are found. Several exact analytical solutions, which evolve in the quasi–homologous regime and satisfy the vanishing complexity factor condition, are exhibited.

1. Introduction

In a recent paper [1], we have presented an approach to describe static, hyperbolically symmetric fluids. The main motivation (but not the only one) behind such an endeavor was the necessity to provide a rigorous description of fluid distributions sourcing the line element
d s 2 = 2 M R 1 d t 2 + d R 2 2 M R 1 + R 2 d Ω 2 , d Ω 2 = d θ 2 + sinh 2 θ d ϕ 2 ,
which in its turn is assumed to be the line element at the interior of the horizon, proposed in [2,3] as an alternative global description of the Schwarzschild black hole.
Such a proposal is motivated by the fact that it is impossible to remove the coordinate singularity in the line element, keeping at the same time the static form of the Schwarzschild metric (in the whole space–time) [4]. Thus, the regular extension of the Schwarzschild metric to the whole space–time may be achieved but at the price to admit a non-static space–time inside the horizon [5,6].
Then, from the belief that any dynamic regime should eventually lead to an equilibrium final state, a static solution has to be expected in the whole space–time.
Accordingly, the model proposed in [2] describes the space time as consisting of two four-dimensional manifolds, the outer one described by the usual Schwarzschild metric on the exterior side of the horizon and the inner one described by (1). A change in signature as well as a change in the symmetry at the horizon are required.
The metric (1) is a static solution admitting the four Killing vectors
K ( 0 ) = t ,
and
K ( 2 ) = cos ϕ θ + coth θ sin ϕ ϕ , K ( 1 ) = ϕ , K ( 3 ) = sin ϕ θ + coth θ cos ϕ ϕ .
Solutions to the Einstein equations endowed with the hyperbolic symmetry (3) have been the subject of research by several authors (see [7,8,9,10,11,12,13,14] and references therein).
Since the fluid that sources the line element (1) is considered as the final state ensuing from a dynamical regime, the obvious question is: What are the general properties of the fluid distribution during this evolving regime, before reaching the equilibrium?
Our purpose in this work is to answer to the above question by carrying on a comprehensive study on the physical properties of evolving fluid distributions in the region inner to the horizon, endowed with the hyperbolical symmetry (3) and that eventually may converge to the static fluid distributions described in [1].
We shall deploy all required equations for a full description of the fluid distribution, including a transport equation. Some specific analytical solutions to these equations will be exhibited. The solutions will be obtained assuming the quasi-homologous condition for their evolution, and the vanishing of the complexity factor.
As we shall see below within the region r < 2 m , where m ( t , r ) is a suitable definition of the mass function, the energy density is negative, and the central region cannot be filled with our fluid distribution. Thus, either the center is surrounded by an empty cavity or by a fluid distribution not endowed with hyperbolical symmetry. A discussion about the physical meaning of the obtained results is presented.

2. The General Setup of the Problem: Notation, Variables and Equations

We consider hyperbolically symmetric distributions of evolving fluids, which may be (or not) bounded from outside by a surface Σ e and, in the case when a cavity is present, are necessarily bounded from inside by a surface Σ i . The fluid is assumed to be locally anisotropic (principal stresses unequal) and undergoing dissipation in the form of heat flow (diffusion approximation).
Having chosen co-moving coordinates, the general interior metric can be written as
d s 2 = A 2 d t 2 + B 2 d r 2 + R 2 ( d θ 2 + sinh 2 θ d ϕ 2 ) ,
where A, B and R are assumed positive, and due to the symmetry (3) are functions of t and r. We number the coordinates x 0 = t , x 1 = r , x 2 = θ and x 3 = ϕ . A and B are dimensionless, whereas R has the same dimension as r.
The energy momentum tensor T α β of the fluid distribution may be written as
T α β = ( μ + P ) V α V β + P g α β + ( P r P ) χ α χ β + q α V β + V α q β ,
where μ , P r , P , q α , V α have the usual meaning, and χ α is unit four–vector along the radial direction. Besides the four–vectors V α , q α and χ α satisfy
V α V α = 1 , V α q α = 0 , χ α χ α = 1 , χ α V α = 0 .
Since the Lie derivative and the partial derivative commute, then
L K ( R α β 1 2 g α β R ) = 8 π L K T α β = 0 ,
implying because of (3) that all physical variables only depend on t and r.
We may also express the energy momentum tensor (5) in the equivalent (canonical) form
T α β = μ V α V β + P h α β + Π α β + q V α χ β + χ α V β
with
P = P r + 2 P 3 , h α β = g α β + V α V β ,
Π α β = Π χ α χ β 1 3 h α β , Π = P r P .
For our comoving observers, we have
V α = A 1 δ 0 α , q α = q B 1 δ 1 α , χ α = B 1 δ 1 α .
It is worth noticing that bulk or shear viscosity could be introduced by redefining the radial and tangential pressures. In addition, dissipation in the free streaming approximation can be absorbed in μ , P r and q.

2.1. Einstein Equations and Conservation Laws

The Einstein equations for (4) and (8) are
8 π μ = 1 R 2 1 B 2 2 B B R R + R R 2 + 2 R R + 1 A 2 2 B ˙ B R ˙ R + R ˙ 2 R 2 ,
4 π q = 1 A B R R B ˙ B + A A R ˙ R R ˙ R ,
8 π P r = 1 R 2 + 1 B 2 2 A A R R + R R 2 + 1 A 2 2 A ˙ A R ˙ R R ˙ 2 R 2 2 R ¨ R ,
8 π P = 1 B 2 A A B B + A A R R B B R R + A A + R R + 1 A 2 A ˙ A B ˙ B + A ˙ A R ˙ R B ˙ B R ˙ R B ¨ B R ¨ R ,
where dots and primes denote the derivative with respect to t and r, respectively. The difference between these equations and the corresponding to the spherically symmetric case are easily identified (see for example Equations (7)–(10) in [15]).
The conservation laws T ; μ μ ν = 0 , as in the spherically symmetric case, have only two independent components, which are displayed in Appendix A.

2.2. Kinematical Variables

The four–acceleration a α and the expansion Θ of the fluid are given by
a α = V α ; β V β , Θ = V α ; α .
From which we obtain for the four–acceleration and its scalar a,
a 1 = A A , a = a α a α = A A B a α = a χ α ,
and for the expansion
Θ = 1 A B ˙ B + 2 R ˙ R .
The shear tensor σ α β is defined by (the vorticity vanishes identically)
σ α β = V ( α ; β ) + a ( α V β ) 1 3 Θ h α β ,
its non zero components are
σ 11 = 2 3 B 2 σ , σ 22 = σ 33 sinh 2 θ = 1 3 R 2 σ ,
and its scalar
3 2 σ α β σ α β = σ 2 ,
reads
σ = 1 A B ˙ B R ˙ R .
All the expressions above are the same, in terms of the metric functions, as in the spherically symmetric case.

2.3. The Weyl Tensor

Using Maple we may easily obtain the Weyl tensor corresponding to our metric (4). Thus, the magnetic part of the Weyl tensor vanishes, whereas its electric part may be written as
E α β = E χ α χ β 1 3 h α β ,
with
E = 1 2 B 2 A A R R A A B B + R R B B + R R 2 + A A R R + 1 2 A 2 A ˙ A B ˙ B A ˙ A R ˙ R + R ˙ R B ˙ B R ˙ R 2 B ¨ B + R ¨ R + 1 2 R 2 .

2.4. The Mass Function

Following [16] we may define the mass function as
m ( r , t ) = R 2 R 232 3 = R 2 R B 2 R ˙ A 2 + 1 ,
where the Riemann tensor component R 232 3 is now calculated for (4).
Introducing the proper time derivative D T , and the proper radial derivative D R by
D T = 1 A t ,
D R = 1 R r ,
we can define the velocity U as
U = D T R ,
which must be smaller than 1 (in relativistic units).
Indeed, in Gaussian coordinates, the position of each fluid element may be given as
x α = x α ( y a , s ) ,
where s is the proper time along the world line of the particle, and y a determines the position of the fluid element on any three-dimensional hypersurface Σ (latin letters running from 1 to 3).
Next, for an infinitesimal variation of the world line we have
δ x α = x α y a δ y a ,
from which it follows
D T ( δ x α ) = V ; β α δ x β .
The position vector of the particle y a + δ y a relative to the particle y a on Σ is defined by
δ x α = h β α δ x β ,
implying that the relative velocity between these two particles is
u α = h β α D T ( δ x β ) ,
then (29) and (30) imply
u α = V ; β α δ x β .
Defining the infinitesimal distance between two neighboring points on Σ by
δ l 2 = g α β δ x β δ x α ,
then it can be shown (see [17] for details) that
δ l D T ( δ l ) = δ x β δ x α σ α β + 1 3 h α β Θ ,
or, introducing the spacelike unit vector
e α = δ x α δ l ,
D T ( δ l ) δ l = e α e β σ α β + Θ 3 .
The above expressions are completely general. Let us now consider our hyperbolically symmetric line element and apply (36) to two neighboring points on a closed curve (S) along the ϕ direction ( r = c o n s t a n t ; θ = c o n s t a n t ). In this case we have e α ( 0 , 0 , 0 , 1 R sinh θ ) and using (16), (18) and (20) in (36) we obtain
D T ( δ l ) δ l = R ˙ A R = U R .
Now, the D T ( δ l ) above is the relative velocity between two neighboring points on S. This quantity of course must be smaller than one (in relativistic units). On the other hand the rate of variation of the total length (L) of S per unit of proper time (say V L ) is also a velocity, and thereof must be smaller than 1, and because of the axial symmetry it is just the sum of (37) over all the curve S. Thus, we have
R L V L = U .
For any value of θ , R L < 1 and thereof U < 1 .
Then, since U < 1 , it follows at once from (23) that m is a positive defined quantity. Additionally, (23) can be rewritten as
E R B = 2 m R + U 2 1 1 / 2 .
Using (23) with (24) and (25) we obtain
D T m = 4 π P r U + q E R 2 ,
and
D R m = 4 π μ + q U E R 2 ,
which implies
m = 4 π 0 r μ + q U E R 2 R d r ,
satisfying the regular condition m ( t , 0 ) = 0 .
Integrating (42) we find
3 m R 3 = 4 π μ + 4 π R 3 0 r R 3 D R μ 3 q U R E R d r .
Since any causal transport equation is based on the assumption that the fluid is not very far from thermal equilibrium, then q < < | μ | . This implies from (42) that μ is necessarily negative, if we assume the condition R > 0 to avoid shell crossing, and remind that m > 0 and E is a regular function within the fluid distribution.
Furthermore, it follows from (42) that whenever the energy density is regular, then m r 3 as r tends to zero. However, in this same limit U 0 , and R r implying because of (39) that the central region cannot be filled with the fluid distribution under consideration. Among the many possible scenarios we shall assume here that the center is surrounded by a vacuum cavity. However, it should be clear that this is just one of the possible choices that, even if having implications on specific models, does not affect the general properties of the fluids endowed with hyperbolical symmetry.
The two above-mentioned features of the fluid appear also in the static case [1].
Before concluding this section it is worth discussing some detail of Equation (A7) and comparing it with the corresponding equation for the spherically symmetric case (see Equation (C6) in [15]).
First of all let us notice that it has the “Newtonian” form F o r c e = M a s s d e n s i t y × A c c e l e r a t i o n . Let us next analyze the different terms in (A7). The first term on the right represents the gravitational interaction, it is the product of the passive gravitational mass density (p.g.m.d) ( μ + P r ), which due to the fact that the energy density is negative, would be negative, and the active gravitational mass (a.g.m) ( 4 π P r R 3 m ) which would also be negative for most equations of state. Thus the gravitational term has the same sign as in the spherically symmetric case. However, its effect is the inverse of the latter case. Indeed, since the p.g.m.d is negative, so is the inertial mass density as consequence of the equivalence principle; thus, the gravitational term tends to increase D T U , i.e., it acts as a repulsive force instead of an attractive one, as in (C6) of [15]. In the same order of ideas we see that a negative pressure gradient would tend to push any fluid element inwardly, i.e., everything happens as if force terms switch their roles, as compared with the spherically symmetric case.

3. The Transport Equation

The treatment of dissipative processes requires the adoption of a heat transport equation. In order to ensure causality we shall resort to the transport equation obtained form the Müller–Israel–Stewart theory [18,19,20].
Then, the corresponding transport equation for the heat flux reads
τ h α β V γ q β ; γ + q α = κ h α β ( T , β + T a β ) 1 2 κ T 2 τ V β κ T 2 ; β q α ,
where κ and τ are two parameters denoting the thermal conductivity and relaxation time, respectively, and T denotes temperature.
There is only one non-vanishing independent component of Equation (44), which may be written as
τ D T q = q κ A B ( A T ) 1 2 τ Θ q 1 2 κ T 2 D T τ κ T 2 q .
In the case τ = 0 we recover the Eckart–Landau equation [21].
Under some circumstances it is possible to adopt the so-called “truncated” version, in which case the last term in (44) is neglected [22],
τ h α β V γ q β ; γ + q α = κ h α β ( T , β + T a β ) ,
and whose only non-vanishing independent component becomes
τ q ˙ + q A = κ B ( T A ) .
Let us now analyze in some detail the changes appearing in the condition for thermal equilibrium, as compared with the spherically symmetric case.
As it was pointed out by Tolman many years ago [23], since all forms of energy have inertia, so should be for thermal energy. This implies, because of the equivalence principle, that thermal energy would tend to displace to regions of lower gravitational potential, thereby modifying the condition of thermal equilibrium in the presence of a gravitational field. In other words, now a temperature gradient is necessary to ensure thermal equilibrium.
Thus, the Tolman condition reads as
T A = 0 T = T A A = T a B .
However, as it follows from (A3), if m > 4 π P r R 3 , in equilibrium a < 0 (the four–acceleration is now directed radially inwardly), implying the existence of a repulsive gravitational force, leading to a positive temperature gradient in order to assure thermal equilibrium. This situation is different from the one observed in the spherically symmetric case, where a negative temperature gradient is required to assure thermal equilibrium.
Before concluding this section it is worth discussing the physical implications of (A8). This equation comes out from the combination of the dynamical Equation (A7) and the transport equation. It brings out the thermal effect on the p.g.m.d., and by virtue of the equivalence principle, on the effective inertial mass density as well. A similar effect was pointed out for the first time for the spherically symmetric case in [24] (see also [25] for a discussion on this effect). In our case the term κ T τ increases the absolute value of the effective p.g.m.d (which is negative), thereby increasing the absolute value of the effective inertial mass density (the term in the bracket on the left of (A8)), as a result of which any hydrodynamic force directed outward tends to push the fluid element inward, weaker than it does in the non-dissipative case, due to the term κ T τ . On the other hand the gravitational term, which is negative, pushes any fluid element as it does in the non-dissipative case. Overall, the thermal effect enhances the tendency to expansion, as in the spherically symmetric case, but different terms in the equation playing different roles as compared with this latter case.
In order to obtain specific solutions to the Einstein equations we shall need to impose additional restrictions. In this work we shall assume that the fluid evolves in the quasi–homologous regime and satisfies the vanishing complexity factor condition. The next two sections are devoted to explain these conditions in some detail.

4. The Structure Scalars and the Complexity Factor

The complexity factor is a scalar function intended to measure the degree of complexity of a self-gravitating system (in some cases more than one scalar function may be required). For a static, hyperbolically symmetric fluid distribution it was assumed in [1] (following the arguments developed in [26]) that the simplest system corresponds to a homogeneous (in the energy density), locally isotropic fluid distribution (principal stresses equal). Thus, a zero value of the complexity factor was assumed for such a distribution. Furthermore, it was shown that a single scalar function (hereafter referred to as Y T F ) describes the modifications introduced by the energy density inhomogeneity and pressure anisotropy, to the Tolman mass, with respect to its value for the zero complexity case.
This scalar belongs to a set of variables named structure scalars and defined in [27], and which appear in the orthogonal splitting of the Riemann tensor [28,29,30,31]. For the sake of completeness we shall highlight the main steps leading to their acquisition (for the spherically symmetric case see [27,32] for details). For our purpose here, we shall need only one of the five structure scalars characterizing our fluid distribution.
The first step consists in defining the tensor Y α β by
Y α β = R α γ β δ V γ V δ ,
which may be split in terms of its trace and its trace-free part as
Y α β = 1 3 Y T h α β + Y T F χ α χ β 1 3 h α β .
Then, using the field equations and (22) the following expressions can be obtained
Y T = 4 π ( μ + 3 P r 2 Π ) , Y T F = E 4 π Π .
On the other hand, combining (12), (13), (22) and (23) we obtain
3 m R 3 = 4 π μ + 4 π Π + E ,
or using (43) and (51)
Y T F = 8 π Π + 4 π R 3 0 r R 3 D R μ 3 q U R E R d r .
Using (12), (13) and (22), we may express Y T F in terms of the metric functions and their derivatives
Y T F = 1 B 2 A A A A R R A A B B + 1 A 2 A ˙ A B ˙ B A ˙ A R ˙ R B ¨ B + R ¨ R .
Following the arguments presented in [1] we shall choose Y T F as the complexity factor. In the dynamic case, however, we still need to provide a criterion for the definition of complexity of the pattern of evolution.
We shall assume here that Y T F is identified with the complexity factor, and we shall consider the quasi-homologous evolution defined in [15] as the simplest mode of evolution.

5. The Quasi-Homologous Condition

In order to provide a rigorous definition of quasi-homologous evolution, let us write (11) as
U R = 4 π q B + σ R R ,
whose general solution is
U = a ˜ ( t ) R + R 0 r 4 π q E + σ R R d r ,
where a ˜ ( t ) is an integration function and (39) has been used.
If the fluid distribution is bounded by a surface Σ e defined by the equation r = r Σ e = c o n s t a n t , we may write
U = R U Σ e R Σ e R r r Σ e 4 π q E + σ R R d r .
The quasi-homologous condition reads
U = R U Σ e R Σ e ,
implying
4 π q E + σ R = 0 .
The above condition will be used to obtain specific models, and its assumption is supported, on the one hand, by the fact that it is the relativistic version of the well-known homologous condition widely used in classic astrophysics, and on the other hand by the fact that it qualifies as one of the simplest patterns of evolution (see [15,32] for a discussion on this point).

6. The Exterior Spacetime and Junction Conditions

In the case that the fluid is bounded, then junction conditions on the boundary have to be imposed [33] in order to avoid the presence of thin shells on the boundary. If any specific model violates Darmois conditions then we should relax the continuity of the second fundamental form, implying the presence of thin shells [34].
Thus, outside Σ e (but inside the horizon) we assume that we have the hyperbolic version of the Vaidya spacetime, described by
d s 2 = 2 M ( v ) r 1 d v 2 2 d r d v + r 2 ( d θ 2 + sinh 2 θ d ϕ 2 ) ,
where M ( v ) denotes the total mass, and v is the retarded time.
The continuity of the first fundamental form reads
( d s 2 ) Σ e = ( d s 2 ) Σ e + ,
where −, + means from the inner or the outer side of the boundary surface, respectively.
At the outer side of the boundary, the surface equation reads
Ψ r r Σ e ( v ) = 0 ,
whose unit normal vector is defined by
n μ + = μ Ψ | α Ψ β Ψ g α β | ,
with components
n μ + = β d r Σ e d v , β , 0 , 0 ,
where
β = 1 | 2 M ( v ) r Σ e 1 + 2 d r Σ e d v | .
At the inner side, the surface equation reads
Φ r r Σ i = 0 ,
whose normal unit vector is defined by
n μ = μ Φ | α Φ β Φ g α β | ,
with components
n μ = 0 , B Σ , 0 , 0 ,
observe that n μ = ( χ μ ) Σ .
From (61) it follows that
R ( t , r Σ e ) = r Σ e ( v ) .
Next, instead of calculating the second fundamental form at both sides of the boundary surface we shall impose the continuity of the flux of energy–momentum across Σ e , which of course implies the absence of thin shells on the boundary surface. For doing so we have to calculate
( T μ ν n ν n μ ) Σ e + ; ( T μ ν n ν n μ ) Σ e ; ( T μ ν n ν V μ ) Σ e + ; ( T μ ν n ν V μ ) Σ e .
where the vectors ( V μ ) + , ( V μ ) have components
( V μ ) + = β , β d r Σ e ( v ) d v , 0 , 0 ,
and
( V μ ) = 1 A , 0 , 0 , 0 .
Next, we have to calculate the energy–momentum tensor corresponding to the line element (60), we obtain
T μ ν ( + ) = 1 4 π r 2 d M d v δ μ 0 δ ν 0 .
From the above expression it follows at once that the energy density of the null fluid sourcing (60) would be negative for an outgoing flux, which is exactly the inverse of what happens for the usual Vaidya metric.
We can now evaluate (70) to obtain
( T μ ν n μ n ν ) Σ e = [ P r ] Σ e ,
( T μ ν n μ V ν ) Σ e = [ q ] Σ e ,
( T μ ν n μ n ν ) Σ e + = β 2 4 π r Σ e 2 d M d v ,
( T μ ν n μ V ν ) Σ e + = β 2 4 π r Σ e 2 d M d v .
Then, imposing the continuity of the flux of energy–momentum across Σ e , it follows that
q = Σ e P r .
where = Σ e means that both sides of the equation are evaluated on Σ e .
Finally, following the usual procedure used in the spherically symmetric case, it is a simple matter to check that the continuity of the second fundamental form implies
m ( t , r ) = Σ e M ( v ) .
In the cases where the central region is surrounded by an empty vacuole bounded by a surface Σ i , junction conditions should be considered also at the inner boundary of the fluid distribution. Then, following the same steps as before we find
P r = Σ i 0 .
and
m ( t , r ) = Σ i 0 .

7. Some Models

In the following subsections we shall exhibit several families of solutions to the Einstein equations for the hyperbolically symmetric fluids. These solutions will be obtained by assuming quasi-homologous evolution and the vanishing of the complexity factor. This choice, justified by the comments in previous sections, will allow us to compare the behavior of our fluid distributions with similar solutions found for the spherically symmetric case in [15]. Besides, different types of additional restrictions will be imposed in order to obtain specific models. It should clear that the purpose of the presentation of these models, besides the possible potential of some of them in the study of specific astrophysical scenario, is to illustrate the richness of fluid distributions endowed with hyperbolical symmetry.

7.1. Non-Dissipative Case

Excluding dissipative processes, and assuming the quasi-homologous condition (59), we may write
q = 0 σ = 0 B ˙ B = R ˙ R R = r B ,
and using (56)
U = R ˙ A = r B ˙ A = a ˜ ( t ) r B .
Imposing next the condition Y T F = 0 we have
A A A A B B A A R R = 0 .
In order to exhibit specific solutions, we shall further assume some additional restrictions.

7.1.1. E = 0 , Π = 0

We shall assume here that the fluid is conformally flat ( E = 0 ), and the pressure is isotropic ( Π = 0 ), which combined with Y T F = 0 produces μ = 0 (i.e., the energy density is homogeneous).
From the conditions E = 0 and Π = 0 we obtain
1 R 2 + 1 B 2 R R 2 + B B R R R R 1 A 2 R ˙ 2 R 2 B ˙ B R ˙ R = 0 .
Using (82) in (85) and (84) produces
1 + r 2 2 R R 2 1 r R R R R = 0 ,
and
A A A A 2 R R 1 r = 0 .
The solution to the system (86) and (87) is easily found to be
R = R ˜ ( t ) cos [ c 1 ( t ) + ln r ] ,
B = R ˜ ( t ) r cos [ c 1 ( t ) + ln r ] ,
A = γ ( t ) R ˜ 2 ( t ) tan [ c 1 ( t ) + ln r ] + b ( t ) ,
where R ˜ ( t ) , c 1 ( t ) , γ ( t ) , b ( t ) are arbitrary functions of their argument. The reader can easily check, using Maple or Mathematica, that the line element (4) with (88)–(90) produces E = 0 = Π = 0 .
To specify further the solution we shall choose the above functions as follows
c ˙ 1 = R ˜ ˙ R ˜ , b ( t ) = γ ( t ) R ˜ 2 ,
producing
R ˙ R = R ˜ ˙ R ˜ ( 1 + tan u ) ,
A = γ ( t ) R ˜ 2 ( 1 + tan u ) , A = a ˜ R ˙ R ,
with a ˜ = γ ( t ) R ˜ 3 R ˜ ˙ and u = c 1 ( t ) + ln r . From the above expressions we find for the physical variables and the mass function
8 π μ = 3 R ˜ 2 + 3 a ˜ 2 , 8 π P r = 8 π P = 3 a ˜ 2 + 3 tan u + 1 R ˜ 2 ( tan u + 1 )
+ 2 R ˜ a ˜ ˙ a ˜ 3 R ˜ ˙ ( tan u + 1 ) ,
m = R ˜ 2 cos 3 u 1 R ˜ 2 a ˜ 2 .
It a simple matter to check that this solution does not satisfy the Darmois conditions at either boundary surfaces, and therefore we must assume the presence of thin shells there.
If we choose R ˜ ( t ) , c 1 ( t ) , γ ( t ) such that they tend to a constant as t , then the above solution tends to the incompressible isotropic solution found in [14], which is a particular case of the hyperbolically symmetric Bowers–Liang solution found in [1].
The above solution might be considered as a version of the Friedman–Robertson–Walker space–time (FRW) for the hyperbolically symmetric case since they share some similar properties, e.g., E = Π = μ = σ = 0 . However it is not geodesic as in the spherically symmetric case. Therefore, we shall next find another version of the hyperbolically symmetric FRW space–time, but satisfying the geodesic condition A = 0 .

7.1.2. Geodesic Solutions

If we further impose the geodesic condition on the fluid, then we may put without loss of generality A = 1 , and the quasi-homologous condition also implies
R I R I I = c o n s t a n t ,
where R I and R I I denote the areal radii of two shells ( I , I I ) described by r = r I = c o n s t a n t , and r = r I I = c o n s t a n t , respectively.
From (97) it follows at once that R is a separable function. In the notation of [32], conditions (83) and (97) define the homologous evolution.
The conditions A = 1 and q = 0 imply
B ˙ B = R ˙ R ,
where (11) has been used. Since the fluid is shear-free we have R = B r , and since R is separable so is B. However, if B is separable, then by a simple reparametrization of r it becomes a function of t alone B = B ( t ) , i.e.,
R = r B ( t ) .
Then (98) is automatically satisfied, as well as Y T F = 0 as it follows from (84). In this case we may write the physical variables and the mass function as
8 π μ = 2 r 2 B 2 + 3 B ˙ 2 B 2 ,
8 π P r = 2 r 2 B 2 B ˙ 2 B 2 2 B ¨ B ,
8 π P = 2 B ˙ 2 B 2 2 B ¨ B ,
m = r B 2 ( 2 r 2 B ˙ 2 ) .
Thus, the fluid is conformally flat, shear-free, geodesic, evolves homologously and satisfies the vanishing complexity factor condition. In this sense it could be considered also as a version of the hyperbolically symmetric FRW space–time. However, unlike the spherically symmetric case, it is anisotropic and the energy density is inhomogeneous.
As in the previous solution, by simple inspection of (101), (103) it can be checked that Darmois conditions cannot be satisfied at either boundary surface.
It is worth analyzing with some detail the differences between this case and the situation in the spherically symmetric case (the usual one). In the latter case we have seen [32] that for a non-dissipative fluid satisfying the homologous condition, the complexity factor vanishes, and there is a single solution characterized by Π = μ = a = E = 0 (FRW).
However, in the present case, imposing homologous conditions on a geodesic non-dissipative fluid we get a conformally flat, shear-free geodesic fluid with Π , μ 0 . If we want to describe an isotropic, homogeneous, shear-free, non-dissipative fluid, then we have to relax the geodesic condition.
Finally, it is instructive to build up a toy model with the above solution, by choosing a particular form for the function B such that asymptotically it leads to a static regime. Thus, let us assume
B = β 1 + e α t ,
where α , β are two positive constants.
Then, it is a simple matter to check that as t we get
8 π μ = 2 r 2 β 2 ,
8 π P r = 2 r 2 β 2 ,
8 π P = 0 ,
and for the mass function we get asymptotically m = r β .
Thus, our toy model converges to the static solution corresponding to the stiff equation of state ( P r = | μ | ) found in [1] (Equations (138)–(139) in that reference).
We shall next consider dissipative solutions.

7.2. Dissipative Case with B = 1

Let us now consider dissipative solutions satisfying the condition B = 1 . As shown in [35], such a condition is particularly suitable for describing situations when the center is surrounded by an empty cavity, a scenario we expect for the kind of fluid distributions we are dealing with in this work.
Thus, the metric functions for this case read
B = 1 , A = R ˙ a ˜ ( t ) R ,
and the corresponding Einstein equations may be written as
8 π μ = 1 R 2 2 R R R R 2 + a ˜ 2 ,
4 π q = a ˜ ( t ) R R ,
8 π P r = 1 R 2 R R 2 + 2 R ˙ R R ˙ R 2 a ˜ ˙ a ˜ R R ˙ 3 a ˜ 2 ,
8 π P = R ˙ R ˙ R ˙ R ˙ R R + R R a ˜ ˙ a ˜ R R ˙ a ˜ 2 .
We may formally integrate (47) producing for the temperature
T ( t , r ) = a ˜ R R ˙ f ( t ) τ a ˜ ˙ 4 π κ ln R 1 4 π κ R ˙ R R R d r τ a ˜ 2 4 π κ ,
where f ( t ) is a function of integration.
On the other hand the condition Y T F = 0 and (59) now read
A A R R + A σ 2 = σ ˙ ,
R ˙ σ R = A .
Introducing the intermediate variables ( X , Y ) ,
A = X + σ ˙ σ 2 and R = X Y ,
(114) and (115) become
X X Y Y + σ 2 = 0 ,
X ˙ X + Y ˙ Y = σ X σ ˙ σ .
Thus, we have a large family of dissipative solutions, among which we shall select some specific ones, by imposing additional restrictions allowing us to integrate (117) and (118).

7.2.1. X Is a Separable Function

If we assume the function X to be separable, then we can integrate the system (117) and (118), obtaining
A = σ ˙ 2 β 2 σ 2 2 β 2 σ 2 ( β r + c 1 ) 2 ,
R = R ˜ 0 σ ( β r + c 1 ) e σ 2 4 β 2 ( β r + c 1 ) 2 ,
a ˜ = σ ,
where β , R 0 ˜ and c 1 are constants.
The above expressions allow us to write for the physical variables
8 π μ = σ 2 e σ 2 2 β 2 ( β r + c 1 ) 2 R ˜ 0 2 ( β r + c 1 ) 2 β 2 ( β r + c 1 ) 2 3 σ 4 4 β 2 ( β r + c 1 ) 2 3 σ 2 ,
4 π q = σ [ 2 β 2 + σ 2 ( β r + c 1 ) 2 ] 2 β ( β r + c 1 ) ,
8 π P r = σ 2 e σ 2 2 β 2 ( β r + c 1 ) 2 R ˜ 0 2 ( β r + c 1 ) 2 4 σ 2 β 2 2 β 2 σ 2 ( β r + c 1 ) 2 + β 2 ( β r + c 1 ) 2 + σ 4 4 β 2 ( β r + c 1 ) 2 ,
8 π P = σ 2 [ 2 β 2 + σ 2 ( β r + c 1 ) 2 ] 2 4 β 2 ( 2 β 2 σ 2 ( β r + c 1 ) 2 ) ,
m = R ˜ 0 ( β r + c 1 ) 2 σ e σ 2 4 β 2 ( β r + c 1 ) 2 1 + R ˜ 0 2 4 σ 2 β 2 4 β 4 + σ 4 ( β r + c 1 ) 4 e σ 2 2 β 2 ( β r + c 1 ) 2 ,
while the expression for the temperature reads in this case as
T ( t , r ) = 2 β 2 σ 2 σ ˙ 2 β 2 σ 2 ( β r + c 1 ) 2 f ( t ) + σ ˙ τ 4 π κ σ 2 4 β 2 ( β r + c 1 ) 2 + ln R ˜ 0 σ ( β r + c 1 ) + σ ˙ 4 π σ κ ln ( β r + c 1 ) σ ˙ σ 3 64 π β 4 κ ( β r + c 1 ) 4 τ σ 2 4 π κ .

7.2.2. A = A ( r )

Another sub-family of solutions may be obtained by assuming that A only depends on r, then the solution to the system (117) and (118) produces
A = 1 4 ( 2 σ 0 r + c 1 ) 2 , a ˜ = σ 0 t σ 1 ,
R = R ˜ ( r ) e 1 4 ( 2 σ 0 r + c 1 ) 2 ( σ 0 2 t 2 + σ 1 t ) ,
where R ˜ ( r ) is an arbitrary function of its argument, and σ 0 , σ 1 , c 1 are constants. To obtain a specific model, we shall further assume R ˜ = R ˜ 0 = constant , in which case we find for the physical variables
8 π μ = σ 1 2 3 σ 0 2 ( 2 σ 0 r + c 1 ) 2 σ 0 2 t 2 + σ 1 t 2 1 R ˜ 0 2 e 1 2 ( 2 σ 0 r + c 1 ) 2 ( σ 0 2 t 2 + σ 1 t ) ,
4 π q = 2 σ 0 2 ( 2 σ 0 r + c 1 ) σ 0 2 t 2 + σ 1 t ( σ 0 t + σ 1 ) ,
8 π P r = 1 R ˜ 0 2 e 1 2 ( 2 σ 0 r + c 1 ) 2 ( σ 0 2 t 2 + σ 1 t ) t 2 σ 0 2 + 2 t σ 0 σ 1 3 σ 1 2 8 σ 0 ( 2 σ 0 r + c 1 ) 2 + σ 0 2 ( 2 σ 0 r + c 1 ) 2 σ 0 2 t 2 + σ 1 t 2 ,
8 π P = 1 2 σ 0 2 t 2 t σ 0 σ 1 σ 1 2 + σ 0 2 ( 2 σ 0 r + c 1 ) 2 σ 0 2 t 2 + σ 1 t 2 ,
m = R ˜ 0 2 e 1 4 ( 2 σ 0 r + c 1 ) 2 ( σ 0 2 t 2 + σ 1 t ) × 1 + R ˜ 0 2 σ 0 2 ( 2 σ 0 r + c 1 ) 2 ( σ 0 2 t 2 + σ 1 t ) 2 ( σ 0 t σ 1 ) 2 e 1 2 ( 2 σ 0 r + c 1 ) 2 ( σ 0 2 t 2 + σ 1 t ) .
For the temperature the corresponding expression reads
T ( t , r ) = 4 ( 2 σ 0 r + c 1 ) 2 f ( t ) τ σ 0 4 π κ ln R ˜ 0 1 4 ( σ 0 2 t 2 + σ 1 t ) ( 2 σ 0 r + c 1 ) 2 ( σ 0 2 t 2 + σ 1 t ) ( σ 0 t + σ 1 ) ( 2 σ 0 r + c 1 ) 2 32 π κ τ ( σ 0 t + σ 1 ) 2 4 π κ .

7.2.3. σ ˙ = 0

Finally, we shall obtain a class of solutions by assuming that the shear scalar is constant, in which case the integration of the system (117) and (118) produces
A = β r β 2 σ t + β 0 , a ˜ = σ = c o n s t .
R = R ˜ 0 β e σ 2 2 r 2 σ β t r + σ 2 β 0 β r + β 2 2 t 2 σ β 0 t ,
where R ˜ 0 , β , β 0 are constants. The physical variables for this case read
8 π μ = σ 2 3 σ 2 r + β 0 β σ β t 2 e 2 σ 2 2 r 2 σ β t r + σ 2 β 0 β r + β 2 2 t 2 σ β 0 t R ˜ 0 2 β 2 ,
4 π q = σ 3 r β σ t + β 0 β ,
8 π P r = σ 2 + σ 4 r β σ t + β 0 β 2 + e 2 σ 2 2 r 2 σ β t r + σ 2 β 0 β r + β 2 2 t 2 σ β 0 t R ˜ 0 2 β 2 ,
8 π P = σ 2 + σ 2 r + β 0 β σ β t 2 ,
and
m = R ˜ 0 β 2 e σ 2 2 r 2 σ β t r + σ 2 β 0 β r + β 2 2 t 2 σ β 0 t × 1 + R ˜ 0 2 β 2 σ 2 σ 2 β 2 ( β r β 2 t σ + β 0 ) 2 1 ) e 2 σ 2 2 r 2 σ β t r + σ 2 β 0 β r + β 2 2 t 2 σ β 0 t ,
T ( t , r ) = f ( t ) ( β r β 2 σ t + β 0 ) + σ 3 12 β 2 π κ ( β r β 2 σ t + β 0 ) 2 τ σ 2 4 π κ .
It can be easily verified that none of the above solutions can be matched smoothly on either of the boundary surfaces.

8. Discussion and Conclusions

We have presented a general approach to describe the dynamics of hyperbolically symmetric fluids, including dissipative processes. Although our main motivation was (and still is) to provide a formalism allowing us to study the dynamic regime leading to a static source of the line element (1), the obtained results are sufficiently general as to be applied to any other scenario where we expect hyperbolical symmetry to play a relevant role.
The four more remarkable features of hyperbolically symmetric fluids are as follows:
  • The energy density is necessarily negative.
  • The fluid cannot fill the central region.
  • The Tolman condition for thermodynamic equilibrium implies in this case the presence of a positive temperature gradient.
  • The thermal modification of the inertial mass density reported for the spherically symmetric case in [24], produces an effect that is similar to the one obtained in this latter case (to enhance the tendency to expansion) but comes about through different terms in the equation.
It should be reminded that the first two properties are common to the static and the dynamic regimes.
With respect to the violation of the weak energy condition ( μ < 0 ) it should be stressed that while it is true that at the classic level we do not expect negative energy density in a realistic fluid, the situation is quite different in the quantum regime, where the appearance of negative energy density is possible (see [36,37,38,39,40] and references therein). This confirms our believe that the type of fluids considered in this manuscript might be useful for studying systems where quantum effects are expected to be relevant.
As mentioned in Section 3, this negative energy density implies the appearance of a repulsive gravitational force, which has two important thermodynamic consequences mentioned in point 3 above.
Next, the impossibility of the fluid distribution to fill the central region leaves several possible scenarios. We lean to assume the existence of an empty vacuole surrounding the center; however, many other scenarios may be regarded as well, such as filling the central region with a fluid endowed with a different type of symmetry. At any rate, this impossibility is consistent with the result obtained in [3], according to which test particles are not allowed to reach the center for the line element (1).
The final description of the central region, as well as the fulfillment or not of the Darmois conditions at both interfaces, would depend on the specific system under consideration.
After having deployed the set of equations for describing the dynamics of hyperbolically symmetric fluids, we have presented several exact solutions. These were found under the condition of the vanishing complexity factor defined in [26] ( Y T F = 0 ) and the quasi-homologous evolution defined by (59).
We first considered the non-dissipative case. Two exact solutions were found for this case. One of them (88)–(95) describes a fluid distribution satisfying conditions Y T F = E = σ = 0 = Π = μ = 0 , which is a reminiscence of the usual FRW space–time. However, unlike the latter it is not geodesic. If we impose the geodesic condition, then the quasi-homologous condition becomes homologous, and the solution is described by (99)–(103). This is a geodesic fluid, satisfying also the conditions Y T F = E = σ = 0 , and therefore is also a good candidate to be regarded as the hyperbolical version of the FRW space–time; however, unlike the latter, it is anisotropic in the pressure and inhomogeneous in the energy density.
In both cases, if the arbitrary functions appearing in the solutions are chosen such that the system tends to a static situation in the limit t , then these solutions tend to the static solutions studied in [1].
Thus, alternative cosmological models emerge from the study of hyperbolically symmetric fluids, which could be of interest when seeking for more sophisticated models of the universe (see for example [41] and references therein).
Finally, we have considered the dissipative case. In order to obtain specific models we have restricted ourselves to the case where the condition B = 1 is satisfied. Such a condition is suggested by the fact that it appears to be suitable for the description of fluids whose central region is surrounded by a vacuum cavity [35]. The purpose of these solutions, as well as the non-dissipative ones, is not the modeling of any specific astrophysical scenario, but just to illustrate a possible way of finding solutions, some of which might be used for the modeling of hyperbolically symmetric fluids required for describing specific physical situations. Neither of the exhibited models matches smoothly on the boundary surfaces. In order to obtain models satisfying Darmois conditions, one could try to extend the general methods developed for the spherically symmetric case in [42,43,44,45] to the hyperbolically symmetric case.
In the temperature profiles exhibited for each solution we may identify two types of contributions: the contribution in the stationary dissipative regime (not containing τ ) and the contribution from the transient regime (terms proportional to τ ).

Author Contributions

All authors contributed equally to this work. Conceptualization, L.H.; methodology, L.H., A.D.P. and J.O.; software, J.O.; formal analysis, L.H., A.D.P. and J.O.; writing—original draft preparation, L.H. writing—review and editing, L.H., A.D.P. and J.O.; funding acquisition, L.H. and J.O. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Ministerio de Ciencia, Innovación y Universidades. Grant number: PGC2018–096038–B–I00, and Junta de Castilla y León. Grant number: SA096P20.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Appendix A. Conservation Laws T ν ; μ μ = 0

In our case the conservation laws have only two independent components which read
μ ˙ + ( μ + P r ) B ˙ B + 2 ( μ + P ) R ˙ R + q A B + 2 q A B A A + R R = 0 ,
and
P r + ( μ + P r ) A A + 2 ( P r P ) R R + q ˙ B A + 2 q B A B ˙ B + R ˙ R = 0 .
Using (12) and (11) we may write
D T U = m R 2 4 π R P r + a E ,
D R U R = 4 π q E + σ R ,
which allows to rewrite (A1) and (A2) as
D T μ + 1 3 ( 3 μ + P r + 2 P ) Θ + 2 3 ( P r P ) σ + E D R q + 2 q a + E R = 0 ,
and
E D R P r + ( μ + P r ) a + 2 ( P r P ) E R + D T q + 2 3 q ( 2 Θ + σ ) = 0 .
Finally, combining (A3) with (A6) we find
( μ + P r ) D T U = ( μ + P r ) ( 4 π P r R 3 m ) 1 R 2 E 2 D R P r + 2 R ( P r P ) E D T q + 2 3 q ( 2 Θ + σ ) .
The above equation may be transformed further by replacing (45) in (A7), and using (A3)
μ + P r κ T τ D T U = μ + P r κ T τ 4 π R 3 P r m 1 R 2 E 2 D R P r + 2 R ( P r P ) κ τ D R T + E q 1 τ + 1 2 D T ln τ κ T 2 5 6 Θ 2 3 σ

References

  1. Herrera, L.; Prisco, A.D.; Ospino, J. Hyperbolically symmetric static fluids: A general study. Phys. Rev. D 2021, 103, 024037. [Google Scholar] [CrossRef]
  2. Herrera, L.; Witten, L. An alternative approach to the static spherically symmetric vacuum global solutions to the Einstein’s equations. Adv. High Energy Phys. 2018, 2018, 8839103. [Google Scholar] [CrossRef] [Green Version]
  3. Herrera, L.; Prisco, A.D.; Ospino, J.; Witten, L. Geodesics of the hyperbolically symmetric black hole. Phys. Rev. D 2020, 101, 064071. [Google Scholar] [CrossRef] [Green Version]
  4. Rosen, N. The nature of Schwarzschild singularity. In Relativity, Proceedings of the Relativity Conference in the Midwest, Cincinnati, OH, USA, 2–6 June 1969; Carmeli, M., Fickler, S.I., Witten, L., Eds.; Plenum Press: New York, NY, USA, 1970; pp. 229–258. [Google Scholar]
  5. Rindler, W. Relativity. Special, General and Cosmological; Oxford University Press: New York, NY, USA, 2001; pp. 260–261. [Google Scholar]
  6. Caroll, S. Spacetime and Geometry. An Introduction to General Relativity; Addison Wesley: San Francisco, CA, USA, 2004; pp. 218–246. [Google Scholar]
  7. Harrison, B.K. Exact Three-Variable Solutions of the Field Equations of General Relativity. Phys. Rev. 1959, 116, 1285–1296. [Google Scholar] [CrossRef]
  8. Stephani, H.; Kramer, D.; MacCallum, M.; Honselaers, C.; Herlt, E. Exact Solutions to Einsteins Field Equations, 2nd ed.; Cambridge University Press: Cambridge, UK, 2003. [Google Scholar]
  9. Gaudin, M.; Gorini, V.; Kamenshchik, A.; Moschella, U.; Pasquier, V. Gravity of a static massless scalar field and a limiting Schwarzschild-like geometry. Int. J. Mod. Phys. D 2006, 15, 1387–1399. [Google Scholar] [CrossRef] [Green Version]
  10. Rizzi, L.; Cacciatori, S.L.; Gorini, V.; Kamenshchik, A.; Piattella, O.F. Dark matter effects in vacuum spacetime. Phys. Rev. D 2010, 82, 027301. [Google Scholar] [CrossRef] [Green Version]
  11. Lobo, F.S.N.; Mimoso, J.P. Possibility of hyperbolic tunneling. Phys. Rev. D 2010, 82, 044034. [Google Scholar] [CrossRef] [Green Version]
  12. Kamenshchik, A.Y.; Pozdeeva, E.O.; Starobinsky, A.A.; Tronconi, A.; Vardanyan, T.; Venturi, G.; Yu, S. Verno. Duality between static spherically or hyperbolically symmetric solutions and cosmological solutions in scalar-tensor gravity. Phys. Rev. D 2018, 98, 124028. [Google Scholar] [CrossRef] [Green Version]
  13. Madler, T. On the affine-null metric formulation of General Relativity. Phys. Rev. D 2019, 99, 104048. [Google Scholar] [CrossRef] [Green Version]
  14. Maciel, A.; Delliou, M.L.; Mimoso, J.P. New perspectives on the TOV equilibrium from a dual null approach. Class. Quantum Gravity 2020, 37, 125005. [Google Scholar] [CrossRef]
  15. Herrera, L.; Prisco, A.D.; Ospino, J. Quasi–homologous evolution of self–gravitating systems with vanishing complexity factor. Eur. Phys. J. C 2020, 80, 631. [Google Scholar] [CrossRef]
  16. Misner, C.; Sharp, D. Relativistic Equations for Adiabatic, Spherically Symmetric Gravitational Collapse. Phys. Rev. 1964, 136, B571. [Google Scholar] [CrossRef]
  17. Herrera, L.; Santos, N.O.; Wang, A. Shearing expansion-free spherical anisotropic fluid evolution. Phys. Rev. D 2008, 78, 084026-10. [Google Scholar] [CrossRef] [Green Version]
  18. Israel, W. Nonstationary irreversible thermodynamics: A causal relativistic theory. Ann. Phys. 1976, 100, 310–331. [Google Scholar] [CrossRef]
  19. Israel, W.; Stewart, J.M. Thermodynamics of nonstationary and transient effects in a relativistic gas. Phys. Lett. A 1976, 58, 213–215. [Google Scholar] [CrossRef]
  20. Israel, W.; Stewart, J.M. Transient relativistic thermodynamics and kinetic theory. Ann. Phys. 1979, 118, 341–372. [Google Scholar] [CrossRef]
  21. Eckart, C. The Thermodynamics of Irreversible Processes. III. Relativistic Theory of the Simple Fluid. Phys. Rev. 1940, 58, 919. [Google Scholar] [CrossRef]
  22. Triginer, J.; Pavón, D. Heat transport in an inhomogeneous spherically symmetric universe. Class. Quantum Gravity 1995, 12, 689–698. [Google Scholar] [CrossRef]
  23. Tolman, R.C. On the weight of heat and thermal equilibrium in general relativity. Phys. Rev. 1930, 35, 904–924. [Google Scholar] [CrossRef] [Green Version]
  24. Herrera, L.; Prisco, A.D.; Hernández-Pastora, J.L.; Martín, J.; Martínez, J. Thermal conduction in systems out of hydrostatic equilibrium. Class. Quantum. Gravity 1997, 14, 2239–2247. [Google Scholar] [CrossRef] [Green Version]
  25. Herrera, L. The inertia of heat and its role in the dynamics of dissipative collapse. Int. J. Mod. Phys. D 2006, 15, 2197–2202. [Google Scholar] [CrossRef] [Green Version]
  26. Herrera, L. New definition of complexity for self–gravitating fluid distributions: The spherically symmetric static case. Phys. Rev. D 2018, 97, 044010. [Google Scholar] [CrossRef] [Green Version]
  27. Herrera, L.; Ospino, J.; Prisco, A.D.; Fuenmayor, E.; Troconis, O. Structure and evolution of self-gravitating objects and the orthogonal splitting of the Riemann tensor. Phys. Rev. D 2009, 79, 064025. [Google Scholar] [CrossRef] [Green Version]
  28. Bel, L. Sur la radiation gravitationelle. C. R. Acad. Sci. Paris 1958, 247, 1094–1096. [Google Scholar]
  29. Bel, L. Radiation states and the problem of energy in general relativity. Gen. Rel. Gravit. 2000, 32, 2047–2078. [Google Scholar] [CrossRef]
  30. Bel, L. Introduction d’un tenseur du quatrieme order. C. R. Acad. Sci. Paris 1959, 248, 1297–1300. [Google Scholar]
  31. Gomez-Lobo, A.G.-P. Dynamical laws of superenergy in general relativity. Class. Quantum Gravity 2008, 25, 015006. [Google Scholar] [CrossRef]
  32. Herrera, L.; Prisco, A.D.; Ospino, J. Definition of complexity for dynamical spherically symmetric dissipative self–gravitating fluid distributions. Phys. Rev. D 2018, 98, 104059. [Google Scholar] [CrossRef] [Green Version]
  33. Darmois, G. Mémorial des Sciences Mathématiques; Gauthier-Villars: Paris, France, 1927; p. 25. [Google Scholar]
  34. Israel, W. Singular hypersurfaces and thin shells in general relativity. Il Nuovo Cimento B 1966, 44, 1–14. [Google Scholar] [CrossRef]
  35. Herrera, L.; Denmat, G.L.; Santos, N.O. Cavity evolution in relativistic self–gravitating fluids. Class. Quantum Gravity 2010, 27, 135017. [Google Scholar] [CrossRef] [Green Version]
  36. Capozziello, S.; Lobo, F.S.N.; Mimoso, J.P. Energy conditions in modified gravity. Phys. Lett. B 2014, 730, 280–283. [Google Scholar] [CrossRef] [Green Version]
  37. Capozziello, S.; Lobo, F.S.N.; Mimoso, J.P. Generalized energy conditions in extended theories of gravity. Phys. Rev. D 2015, 91, 124019. [Google Scholar] [CrossRef] [Green Version]
  38. Barcelo, C.; Visser, M. Twilight for the Energy Conditions? Int. J. Mod. Phys. D 2002, 11, 1553–1560. [Google Scholar] [CrossRef] [Green Version]
  39. Kontou, E.A.; Sanders, K. Energy conditions in general relativity and quantum field theory. Class. Quantum Gravity 2020, 37, 193001. [Google Scholar] [CrossRef]
  40. Pavsic, M. On negative energies, strings, branes, and braneworlds: A review of novel approaches. Int. J. Mod. Phys. A 2020, 35, 2030020. [Google Scholar] [CrossRef]
  41. Kopteva, E.; Bormotova, I.; Churilova, M.; Stuchlik, Z. Accelerated Expansion of the Universe in the Model with Nonuniform Pressure. Astrophys. J. 2019, 887, 98. [Google Scholar] [CrossRef]
  42. Thirukkanesh, S.; Maharaj, S.D. Radiating relativistic matter in geodesic motion. J. Math. Phys. 2009, 50, 022502. [Google Scholar] [CrossRef] [Green Version]
  43. Thirukkanesh, S.; Maharaj, S.D. Mixed potentials in radiative stellar collapse. J. Math. Phys. 2010, 51, 072502. [Google Scholar] [CrossRef] [Green Version]
  44. Ivanov, B. All solutions for geodesic anisotropic spherical collapse with shear and heat radiation. Astrophys. Space Sci. 2016, 361, 18. [Google Scholar] [CrossRef] [Green Version]
  45. Ivanov, B. A different approach to anisotropic spherical collapse with shear and heat radiation. Int. J. Mod. Phys. D 2016, 25, 1650049. [Google Scholar] [CrossRef] [Green Version]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Herrera, L.; Di Prisco, A.; Ospino, J. Dynamics of Hyperbolically Symmetric Fluids. Symmetry 2021, 13, 1568. https://doi.org/10.3390/sym13091568

AMA Style

Herrera L, Di Prisco A, Ospino J. Dynamics of Hyperbolically Symmetric Fluids. Symmetry. 2021; 13(9):1568. https://doi.org/10.3390/sym13091568

Chicago/Turabian Style

Herrera, Luis, Alicia Di Prisco, and Justo Ospino. 2021. "Dynamics of Hyperbolically Symmetric Fluids" Symmetry 13, no. 9: 1568. https://doi.org/10.3390/sym13091568

APA Style

Herrera, L., Di Prisco, A., & Ospino, J. (2021). Dynamics of Hyperbolically Symmetric Fluids. Symmetry, 13(9), 1568. https://doi.org/10.3390/sym13091568

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop