Black hole acoustics in the minimal geometric deformation of a de Laval nozzle

The correspondence between sound waves, in a de Laval propelling nozzle, and quasinormal modes emitted by brane-world black holes deformed by a 5D bulk Weyl fluid are here explored and scrutinized. The analysis of sound waves patterns in a de Laval nozzle in the laboratory, reciprocally, is here shown to provide relevant data about the 5D bulk Weyl fluid and its on-brane projection, comprised by the minimal geometrically deformed compact stellar distribution on the brane. Acoustic perturbations of the gas fluid flow in the de Laval nozzle are proved to coincide with the quasinormal modes of black holes solutions deformed by the 5D Weyl fluid, in the geometric deformation procedure. Hence, in a phenomenological Eötvös–Friedmann fluid brane-world model, the realistic shape of a de Laval nozzle is derived and its consequences studied.


Introduction
General relativity (GR) is a successful theory, widely tested by experiments and observations, however a limited one as regards some recent questions, like the nature of dark energy/dark matter. GR can be recovered from models involving higher dimensions as a very restricted case [1,2]. In brane-world models, the brane self-energy density is manifest as the brane tension (σ ), which is assumed to be infinite in the GR limit. Nevertheless, a finite value for the brane tension, in codimension one models, is an ubiquitous feature of the brane that also comprises the warped five-dimensional (5D) geometry, besides the brane self-gravity.
Although an enormous brane tension value recovers GR at low energy regimes, phenomenological evidence indicates a variable brane tension [3][4][5][6]. At an immensely hot universe, the tension of the brane had attained a nugatory value. Afterwards, the brane tension increased and the universe cooled down [3][4][5]. Fluid membranes play a central role in modela e-mail: roldao.rocha@ufabc.edu.br ing this scenario, wherein a temperature-dependent tension is ruled by the Eötvös principle [6], which governs a Friedmann brane with a scale factor that drives the expansion of the universe [7].
4D gravity can be effectively formulated on a brane-world by two complementary methods. The first one comprises the Shiromizu-Maeda-Sasaki implementation of the Gauss-Codazzi on-brane projection routine [8]. Nevertheless, this method does not represent a consistent system of equations, since the 5D bulk Weyl tensor cannot be determined from data on the brane. In fact, there is no action whose projected Euler-Lagrange equations onto the brane can be derived [2]. A complementary technique does involve an action that, at low energy regimes, derives the effective 4D theory described by the respective Euler-Lagrange equations [9]. Hence, these two complementary procedures can be employed in the construction of an effective 4D theory [5].
Among successful efforts to formulate theories beyond GR, the procedure consisting in accomplishing a minimal geometric deformation of the Schwarzschild solution in a brane-world was derived [10][11][12]. It comprises exact solutions of the Einstein equations for the 4D effective theory on the brane [10][11][12]. The minimal geometric deformation procedure incorporates high-energy improvements to GR, when the (brane) vacuum state is percolated by a 5D Weyl fluid in the bulk [10,13,14]. These analytical solutions -of the 4D brane effective Einstein's field equations -encode compact distributions supporting stellar systems that can even exhibit solid crusts [13,14], driven by the 5D bulk Weyl fluid, and also peculiar generalizations [15]. The deformation itself comprises the brane tension as the managing parameter of high-energy regimes. This setup has GR corresponding to an ideally rigid brane (σ → ∞), at low energies. A more refined setup can be implemented by considering a variable tension fluid brane [3,16]. This approach has been comprehensively and successfully constrained by experimental and observational bounds, provided by the perihelion precession of Mer-cury, the deflection of light by the Sun, the gravitational redshift and the radar echo delay, recently obtained in Ref. [17]. Besides, observational lensing effects by minimal geometrically deformed black holes have a typical signature that may soon be probed by the European Space Agency satellite mission [18]. Despite this comprehensive list of experimentally and observationally possible signatures, regarding the brane-world black hole that is geometrically deformed by a 5D Weyl fluid, any study of the quasinormal modes produced by this kind of black holes is still lacking.
On the other hand, acoustic perturbations of a gas flow, in the so-called de Laval nozzle, have been shown to correspond to the general form of perturbations of Schwarzschild black holes [19][20][21]. It introduces the feasibility to produce and observe their quasinormal resonances in the laboratory. The de Laval nozzle is an example among propelling nozzles, which are widely studied devices that turn a fluid (gas) turbine into a jet engine. A de Laval nozzle is constituted by an hour-glass-shaped tube, strained in the middle, utilized to accelerate a hot pressurized gas to a high supersonic speed into the thrust direction. The fluid thermal energy is commuted into kinetic energy, and the fluid velocity increases. The energy to accelerate the gas stream induces the gas to adiabatically expand with high efficiency to a final -transonic, supersonic, or even hypersonic speed -propelling jet. The de Laval nozzle is constructed upon the theory of quasi-1D flows, where a fluid moves at the magnitude of the speed of sound. In this regime, the changes in the fluid density turn significant and compress the flow. The de Laval nozzle is based upon the Venturi effect.
Fluid flows have been studied, in this context, in a de Laval nozzle, aiming to observe acoustic black holes in the Schwarzschild setup [19][20][21]. The acoustic black hole surface gravity was experimentally derived in the laboratory, in Ref. [22]. As argued in Refs. [20,21], sonic regions in a fluid can induce a surface for the sound waves, known as the acoustic horizon, that emulates a black hole event horizon. Perturbations of sound waves have been shown to be analog to the quasinormal modes, corresponding to black holes gravitational excitations [23,24], being, moreover, successfully explored in different contexts [25][26][27].
Our point here is to derive and analyze the correction to a de Laval nozzle trend, using also its analogy to a braneworld black hole in the minimal geometric deformation setup, regarding a variable brane tension. Besides, another goal here is to study the analogy between waves in a de Laval nozzle in the laboratory and quasinormal modes of minimal geometrically deformed brane-world black holes. Hence, to scrutinize sonic waves in a de Laval nozzle can circumvent the indeterminacy of the 5D Weyl tensor on the bulk [2,28], which encrypts the bulk geometry. Since the minimal geometric deformation is generated by a 5D Weyl fluid on the bulk, experiments regarding a de Laval nozzle in the labo-ratory may reciprocally provide relevant data about the 5D bulk Weyl fluid. This paper is organized as follows: Sect. 2 is devoted to a brief review, regarding minimal geometrically deformed compact systems, further refined by phenomenological Eötvös-Friedmann fluid branes. In Sect. 3, a gas flow is perturbed in a de Laval nozzle, whose wave equation is analog to the wave equation regarding spin-s perturbations of minimal geometrically deformed brane-world black holes. Hence, the current bound for the variable brane tension provides corrections to the expression for the trend of de Laval nozzles. Moreover, we here propose to study quasinormal modes from these black holes in the laboratory, when the wave equation in a de Laval nozzle equals the wave equation of spin-s perturbations of brane-world black holes undergoing a minimal geometric deformation. Section 4 is dedicated to a discussion of our results, to summarizing the conclusions and to providing relevant perspectives.

The minimal geometric deformation setup and fluid branes
Acoustic analogs of brane-world black holes have been reported [29], in the context of setups due to both Randall and Sundrum [30] and to Dvali et al. [31]. However, no approach has regarded realistic data on the brane yet, comprising the variable brane tension paradigm on fluid branes. Employing the minimal geometric deformation technique incorporates high-energy refinements to general relativity, by permeating the brane vacuum outer of a compact distribution with a Weyl fluid in the 5D bulk [10,13,14]. Brane-world models encompassing a variable brane tension are best implemented by Eötvös-Friedmann fluid branes, where the brane temperature drives the brane tension, according to the Eötvös' rule across the universe expansion [3,6]. The most stringent brane tension bound σ 3.2 × 10 −6 GeV 4 has been provided in the context of the minimal geometric deformation of black holes, formed as Bose-Einstein condensates of gravitons that weakly interact among themselves [32]. The associated entropic information content can also predict the Chandrasekhar critical density of compact objects in this paradigm [32,33].
The effective 4D Einstein equations can be derived by the Gauss-Codazzi on-brane projection method, from the 5D bulk equations [2,8]. In natural units, the 4D Einstein effective equations were obtained in Ref. [8] (Greek indices run in the set of Minkowski space-time indices): where G μν denotes the Einstein tensor and the cosmological parameter on the brane is denoted by . The effective stress-energy tensor T μν = T μν + σ −1 S μν + E μν + P μν + L μν encrypts the matter stress-energy tensor on the brane (T μν ), the 5D bulk Weyl tensor electric projection on the brane (E μν ) -that comprise data (constituting a Weyl fluid) about the gravitational field out of the brane -and S μν is the traceless irreducible component, proportional to the brane extrinsic curvature [8] -regards 5D effects onto the brane; the tensor components L μν encrypt the asymmetric embedding of the brane, into the bulk, and P μν stands for the pullback onto the brane of the stress-energy tensor, designating eventual 5D non-standard model fields, comprising radiation of quantum origin, dilatonic, and even moduli fields [2,7,8]. Deviations from the usual Einstein standard equations in GR can also be generated by excitations of 5D gravitons, whose effects are encompassed in the P μν tensor. Eventually, some of the terms constituting T μν can equal zero. Exact solutions of the 4D effective Einstein equations are rare, due to the intricacy of the system of equations and to initial data out of the brane as well. Compact distributions, modeling stellar structures, are spherically symmetric, static solutions of Eq. (1), of type for d 2 representing the 2-sphere surface element. The deformation on the radial component in Eq. (2) is caused by the bulk constituents, encrypting not only anisotropic effects originating from the bulk gravity but also the effects of a 5D Weyl fluid in the bulk, whose brand permeates the brane vacuum.
The minimal geometric deformation is implemented by fixing the temporal component in (2) and deforming the outer radial component [12,14], where for [12] The prime denotes the derivative with respect to the radial coordinate, and R denotes the compact distribution effective radius [10]. The ς parameter in (3) regards the Weyl fluid in the bulk and its induced deformation of the brane 4D vacuum [17]. The region inner to the stellar distribution is regular at the origin. The inner and outer region to a star have a shared boundary constituted of a solid incrustation. The (variable) brane tension and the stellar effective radius are parameters that determine the star crust width [13,14]. The outer region r > R [10] then promotes the deformed metric [12] References [10,16] show that the metric radial component (6a) can be split as follows -in this section the subindex "0" refers to the GR limit σ → ∞ (or, equivalently, ς = 0): where where B ς (r ) evinces a high-energy correction to the Schwarzschild solution to order The parameter ς is proportional to the stellar distribution compactness and drives the geometric deformation of the Schwarzschild solution, having an explicit expression in terms of the brane tension [10,12]: for α ≡ d 0 R 2 and Typically d 0 0.27492 R 2 . The GR σ → ∞ limit yields ς = 0 in Eqs. (6a) and (6b), leading to the standard Schwarzschild metric solution of Einstein's equations.
The current experimental and observational data was shown to enforce the strongest bound |ς | 6.1 × 10 −11 (obtained by the perihelion precession classical test of GR) and the weakest bound |ς | 8.2 × 10 −5 (derived from the radar echo delay classical test of GR) on the adimensional deformation parameter, in Ref. [17]. Besides, the most recent brane tension bound σ 3.2×10 −6 GeV 4 has been obtained by the informational entropy of the minimal geometrically deformed Bose-Einstein condensate of gravitons [32]. Equation (9) implies a negative value for the parameter ς . Hence, the gravitational field strength is mitigated by the finite value of the brane tension and by the 5D encompassing scenario [17], which attains a maximum at the stellar surface r = R. Denoting the stellar distribution density by ρ star , the bound σ 3.2 × 10 −6 GeV 4 complies with the condition ρ star σ 1 [2,3].
Eötvös-Friedmann fluid branes are known to have a variable brane tension proportional to the temperature of the universe, T . The Eötvös law asserts that σ ≈ T − τ , [3,34], for τ a crucial constant parameter that drives σ into positive values, subsequently to the Big Bang [3,6]. The scale factor constant value a 0 fixes the beginning of the universe at a τ temperature [3,34]. Reference [3] computed the temperature dependence upon the scale factor as T (t) ≈ 1 a(t) [3], yielding a variable brane tension that is time dependent [3,4], where κ 2 denotes the 4D coupling "constant" and κ 2 0 = (8π G) −1 is the late-time coupling constant, G denoting the Newton's constant.
At early times, until the radiation density equated the brane matter density, the brane tension could be taken as having a low value. However, the brane tension and the 4D coupling parameter likewise magnified as the universe expanded. The time-dependent brane tension expression yields 4D [3,4]. Black string and black brane solutions in a variable brane-world context were studied in Refs. [16,35,36].

The minimal geometric deformation of a de Laval nozzle
A de Laval nozzle is intrinsically based on the Venturi effect. When the flow of a (gas) fluid passes through a constricted part of a tube with variable cross section A(x), it generates a reduction in the fluid pressure, whereas the fluid velocity increases. Modeling de Laval nozzles considers quasi-1D flows, which are isentropic, adiabatic, and frictionless ones; these shall be considered here. The fluid studied can consist of an ideal gas and expressed by the equation of state p = ρ RT , where p denotes the fluid pressure, T is the absolute temperature, and R is the universal gas constant. An ideal gas is well known to have a constant heat capacity, at constant pressure and constant volume -respectively denoted by C p and C V . Hence R = C p −C V and the specific heat ratio reads γ = C p /C V . For instance, the heat capacity ratio for helium is γ = 1.66, whereas nitrogen has γ = 1.4. Here diatomic gas molecules shall be regarded. An isentropic gas flow, from an initial to a final state, has the property [37] where these quantities shall be normalized by the initial state. Prominent properties of isentropic flows comprise the uniform expansion of the gas, composing then a shock-free, continuous, flow. A relevant parameter of a compressible flow is the Mach number, , where c 2 s = d p dρ = γ RT is the local speed of sound and x denotes the transversal nozzle coordinate, namely, the coordinate along the de Laval nozzle, and v is the local flow velocity. The Mach number is employed to categorize the distinct regimes of flow. 1 Besides, the mass flow rate dm dt is the flux per unit throat area ρ Av, meaning the mass of the gas that passes through a cross section of the tube per unit time, also known as the fluid discharge [37]. A quasi-1D fluid flow is ruled by the Euler-Lagrange equations and the continuity relation in fluid mechanics, given by [37] ∂ ∂t Instead of Eq. (13b), one can use the Euler equation associated with the Bernoulli equation where the last term in Eq. (15) represents the heat function of a barotropic fluid, identified with the enthalpy, and = v dx denotes the velocity potential. From Eq. (15), the linearized equation for sound waves can then be obtained, considering perturbations φ and δρ, respectively, around and ρ [19,20].
A quasi-1D fluid flow in a de Laval nozzle has a stagnation state, which is a state attained by the fluid if it is conveyed to rest into an isentropic state and without work. The stagnation speed of sound is denoted by c s0 . Now, the acoustic analog of the tortoise coordinate, x , is defined by Perturbing the system of equations (13a)-(13c) in a nozzle yields [19] where the associated potential, representing the sound waves curvature scattering on the acoustic black hole, reads for [19,20] where A de Laval nozzle trend is constructed upon the nozzle throat cross-sectional area. Dimensionless quantities for ρ(x) and A(x) are obtained by measuring them in units of the stagnation gas density ρ 0 and of the throat nozzle crosssectional area, respectively. Moreover [20,21], which by Eq. (19) yields g = ρ 1−γ 2(ρ 1−γ −1) 1/2 , following which we have yielding The Mach number equals the unit at the horizon, wherein thus g has to be finite, Replacing Eq. (22) into Eq. (19) implies the cross-sectional nozzle area expressed in terms of g [19], On the other hand, the analogy between fluid flows in a de Laval nozzle and the brane-world black hole -undergoing a minimal geometric deformation -can be implemented. In fact, scalar field perturbations in the minimally geometric deformed brane-world black hole background are known to yield the wave-like equation [25] where dr * = dr B(r ) and the effective potential for the quasinormal ringing of the brane-world black hole under the minimal geometric deformation reads Equation (27) is analog to Eq. (17). In fact, to find a scalar function g that produces the same potential, the tortoise coordinate of the black hole solution is identified with the de Laval nozzle, dr * = dx , yielding The differential equation for g(r ) then reads One can elect an unit event horizon radius, so that the nozzle coordinate to be written with respect to the event horizon.
In the limit ς → 0, Eq. (7) is reduced to Eq. (30) in Ref. [20], which has the solution in Eq. (31) in that reference. Substituting Eq. (7) into (30) yields an intricate equation that cannot be analytically solved. However, by splitting the solution of (30) into a sum of a purely GR component (g 0 (r ) ≡ lim σ →∞ g(r )) and a component that is induced by the 5D Weyl fluid, we can substitute the solution of (30) for ς = 0, obtained in Ref. [20], to find the g ς (r ) function, iteratively solving Eq. (30). The solution of Eq. (30) has two integration constants, determined by Eq. (24). Equations (29) and (6b) then provide the nozzle coordinate x with respect to r , The integral lower limit can be made consistent with the fact that the coordinate x is null at the sonic point, by imposing r 0 = 1. Hence, the de Laval nozzle cross section A(x) can finally be derived, modeling the nozzle shape. In fact, Eq. (32) can be put into Eq. (30), whose numerical solutions for g(r ) yield the corrections (due to the 5D Weyl fluid) g ς (r ), in Eq. (31). Subsequently, we rewrite the numerical solution g(r ) in Eq. (31) with respect to the transversal nozzle coordinate x in Eq. (32), substituting it into the expression for the nozzle cross section A(x) in Eq. (26).
In what follows, the solid gray areas in Figs. 1, 2 and 3 indicate the cross section A(x) of the de Laval nozzle and its shape, in the σ → ∞ GR limit (gray area limited by the continuous gray line) and its minimal geometric deformation due to a 5D bulk Weyl fluid (gray area limited by the dotted gray line). The black strips, respectively, represent the effective potential for those two cases. The brane tension bound adopted, σ ≈ 3.2 × 10 −6 GeV 4 , was derived in Ref. [32] through the numerical computation of the minimal geometric deformation case, as well as the weakest experimental bound |ς | 8.2 × 10 −5 on the minimal geometric deformation parameter [17].
The minimal geometric deformation constricts the nozzle throat cross-sectional area. In addition, by comparing the Schwarzschild versus minimal geometric deformation, the nozzle corrections due to the influence of a 5D bulk Weyl fluid that permeates the vacuum on the brane are notorious. Hence, the quasinormal modes of black holes solutions deformed by the 5D bulk Weyl fluid can be probed by their analogy with acoustic perturbations of a diatomic gas fluid flow in a de Laval nozzle, using the geometric deformation technique. Reciprocally, the signature and the lacking data about the Weyl tensor in the bulk can be probed in the laboratory, by analyzing sonic waves throughout a minimal geometrically deformed de Laval nozzle. The next section is devoted to further exploring and analyzing the consequences of our results.

Concluding remarks and outlook
The minimal geometric deformation of a de Laval nozzle can have double-handed applications. The de Laval nozzle associated with black hole analogs produced in the laboratory can present their trend slightly modified by 5D bulk Weyl fluid effects. On the other hand, 5D effects can also be probed in the laboratory, due to the analogy presented here.
Using a phenomenological Eötvös-Friedmann fluid brane setup, describing an inflationary brane-world universe, the perturbation of a fluid flow in a de Laval nozzle was considered, providing a wave equation that is similar to the wave equation regarding perturbations of minimal geometrically deformed brane-world black holes. The precise bounds for the variable brane tension value provided corrections to the shape of de Laval nozzles in this context. Figures 1, 2 and 3 plot the de Laval nozzle profile and also the analysis of the nozzle deformation with respect to the Schwarzschild solution. Such a deformation is generated by a 5D Weyl fluid permeating a compact distribution described by the Schwarzschild metric solution of the 4D brane effective Einstein equations.
Moreover, quasinormal modes of brane-world black holes undergoing a minimal geometric deformation can then be produced and observed in the laboratory, by analyzing the sonic waves throughout the associated deformed Laval nozzle. Hence, the solution for the inverse technique, consisting of the correspondence between the shape of the de Laval nozzle and the general trend of perturbations in brane-world black holes deformed by a minimal geometric deformation of Schwarzschild black holes, has been implemented here. The corrections to the Schwarzschild solution on the brane, permeated by a 5D bulk Weyl fluid, affect how the pressure is dispersed across the deformed de Laval nozzle. The finite brane tension then specifies a protocol for the analysis telling us whether the de Laval nozzle highest thrust can be achieved and where we search for the best flow properties that are being attained, for the derived de Laval nozzle shape.
Using the sonic analog to black holes [38], the thermal spectrum of sound waves was given from the sonic horizon in transonic fluid flows, also in the context of analog gravity [39][40][41]. These approaches can be further explored, using the methods here introduced, together with more fluid analogy phenomena regarding black holes in the laboratory [42,43]. Still, further types of black holes can be studied [27,44]. Finally, the extended MGD approach [45] can also be used to derive further corrections to the de Laval nozzle profile.