Next Article in Journal
Avoidance of Singularity during the Gravitational Collapse with String T-Duality Effects
Previous Article in Journal
Recent Cross-Section Measurements of Top-Quark Pair Production in Association with Gauge Bosons
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hot Spots in Sgr A* Accretion Disk: Hydrodynamic Insights

by
Elizabeth P. Tito
1,*,
Victor P. Goncharov
1,2 and
Vadim I. Pavlov
1,3
1
Scientific Advisory Group, Pasadena, CA 91125, USA
2
A. M. Obukhov Institute of Atmospheric Physics RAS, Moscow 109017, Russia
3
Faculté des Sciences et Technologies, Université de Lille, F-59000 Lille, France
*
Author to whom correspondence should be addressed.
Universe 2023, 9(1), 40; https://doi.org/10.3390/universe9010040
Submission received: 27 November 2022 / Revised: 2 January 2023 / Accepted: 4 January 2023 / Published: 8 January 2023
(This article belongs to the Section Compact Objects)

Abstract

:
The recent image of our galaxy’s supermassive black hole Sgr A* derived from the 7 April 2017 data of the Event Horizon Telescope Collaboration shows multiple hot spots in its accretion disk. Using the analytical framework, we demonstrate that the observed hot spots may not be disjoint elements but causally linked components (“petals”) of one rotating quasi-stationary macro-structure formed in the thermo-vorticial field within the accretion disk.

1. Introduction

The image of our galaxy’s central supermassive black hole Sagittarius A* (Sgr A*), derived recently by the Event Horizon Telescope (EHT) Collaboration (see Figure 1A and Refs. [1,2,3,4,5,6]) shows a multi-spot structure of its accretion disk. The disk structure is a product of complex state-of-the-art data analysis rather than a direct observation.
The EHT—a collection of radio-telescopes scattered around the Earth—operates in the digital interferometer mode: the signal from each antenna is recorded, and then the image of the object is restored using correlation analysis. Sophisticated data-processing algorithms have permitted the EHT to achieve angular resolution on the order of 20 microarcseconds. At the level of sensations, this is equivalent to the ability to read newspaper headlines on the Moon. However, as Figure 1A indicates, this resolution scale is comparable to the size of Sgr A* itself; the accretion disk is slightly greater (∼ 50 μ as ). Furthermore, the EHT telescopes could only record data from a small study area for a short period of time (see colored zones in Figure 2). Many (white) parts have remained unexplored. To restore the full mosaic, the algorithms had to fill the gaps.
The shape of the observed structure (Figure 1A)—even if the structure is short-lived—appears to indicate that it is likely to be not an artifact of image-reconstruction algorithms, but a real phenomenon. Using the analytical framework, we demonstrate that the observed hot spots may be not disjoint but causally linked components (“petals”) of one rotating quasi-stationary macro-structure formed in the thermo-vorticial field within the accretion disk.
Indeed, when a black hole’s accretion disk—whose rotation axis is perpendicular to the disk plane—is heated non-homogeneously (so temperatures are higher near the outer edge of the disk), then, in the field of the centrifugal force, spontaneously self-formed hot “bubbles” (composed of locally clustered plasma with temperatures in excess of the “average”, hence with lower densities) should move towards the axis of the disk rotation. However, when the hot “bubbles” are also vortices, then each such vortex (via the induced velocity field) “forces” other vortices to rotate around itself, hence diverting their motion “sideways”, curtailing the movement towards the central axis of accretion-disk rotation. All these vortices are subject to the influence of the cumulative velocity field induced by all other vortices. Thus, the radial motion of the vortices towards the axis becomes suppressed. As a result, if stabilized, the vortices take positions equidistantly from the axis and self-organize into a symmetric thermo-vorticial macro-structure that rotates as a whole around the mutual center (like in Figure 1B). The dynamics and longevity of this structure are linked to the thermal and vortical properties of the system and its elements. Visually, if observed, the petals of this structure look like bright “hot spots”.
In this paper, we consider the EHT image from the perspective of theoretical hydrodynamics. In particular, we describe a model of large-scale stationary rotating heated vortices. However, we consider not the usual hydrodynamic field of vorticity but a complex thermo-hydrodynamic field system that under certain circumstances may self-organize into regular structures. The physical and mathematical underpinnings of this analytical approach are elaborated in the references provided in the relevant places. The explanation of their details is beyond the scope of this paper. To avoid any confusion, let us also emphasize upfront that we work with the field, not with individual particles (their trajectories or orbits). Perhaps what may help the reader grasp this nuance better is the reminder that the velocity of displacement of electrons in a usual house wire is not the same thing as the speed of propagation of the electro-magnetic field perturbation along that same wire. As the result of our analysis, we show that multi-hot-spot thermo-vorticial macro-structures may indeed self-organize in the accretion disk. The model makes it possible to determine basic characteristics of such structures, for example, the horizontal space-scale, the period of proper rotation, and the peak temperature magnitude in the vortex.
The paper is organized as follows: Section 2 presents the model, Section 3 presents the results, and Section 4 summarizes the conclusions.

2. Model

The model setup is straightforward: a black hole pulls in and crushes the matter from the surrounding space; the particles are then accelerated to near-light velocities and twisted around the black hole, forming a flattened accretion plasma disk in the equatorial plane.
We will use the spacetime metric entirely characterized by the black hole mass parameter and its “spin” (described in our Refs. [7,8]; for more details, see also [9,10,11,12,13,14], and bibliographies therein). We will assume that the mass of the accretion disk is negligible compared to the black hole “mass” M, probably ( 10 5 ÷ 10 4 ) M ; the radiative cooling does not strongly affect the dynamics of fluid motion; and the electrons and ions are very weakly coupled by Coulomb interaction and hence ions and electrons plasmas components have different temperatures, T e T i (see Ref. [15]), and thus it is the electron component that contributes the most to the equation of state of the accretion disk matter. Due to the large difference in the masses of electrons and protons, electrons are highly mobile and provide quasi-neutrality of the plasma. Due to the high conductivity of the plasma, its own magnetic field can be considered as a field “frozen” into the plasma.
Generally speaking, equations of fluid motion in the vicinity of a black hole must be written using the concept of relativistic dynamics. The key points are as follows: We suppose that the space-time near the (non-charged) black hole Sgr A* is described by the Kerr metric—an exact, singular, stationary, and axially symmetric solution of the Einstein–Hilbert equations of the gravitational “field” in vacuum. Next, we introduce the Boyer–Lindquist 4-coordinates, q α = ( t , r , θ , ϕ ) (it is well known that besides the Boyer–Lindquist coordinate representation, other representations of space-time locations exist). In terms of the Boyer–Lindquist coordinates, the square of interval is written as d s 2 = g α β ( r / r g , θ ) d q α d q β with α , β = 0 , 1 , 2 , 3 , i.e., the components of g α β depend only on the dimensionless combination r g / r and θ . Here, r g = 2 G M / c 2 is the Schwarzschild radius, c is the speed of light, G is the gravitational constant, and M is the “mass” of the black hole. The off-diagonal term g 03 in the metric tensor is proportional to the rate of the black hole’s own rotation and to 1 / r . For the Minkowski tensor, we use the metric signature d i a g ( + ) (see, for example, Ref. [8], and Refs. therein). To satisfy the principle of causality for moving material objects, obviously, d s 2 > 0 . The four time-space coordinates q α = ( t , r , θ , ϕ ) give the location of a world-event from the viewpoint of a remote observer. The meaning of space coordinates r , θ , ϕ is clear once transitioned to the limit r r g , r ω r g 2 / c . When the square of the interval becomes d s 2 c 2 d t 2 d r 2 r 2 ( d θ 2 + s i n 2 θ d ϕ 2 ) , i.e., at infinity, parameters r , θ , ϕ may be interpreted as the standard spherical coordinates in flat space-time. As for the parameter r, strictly speaking, note that it is not the “distance” in the usual meaning from the center of black hole. This is because, for any material object, in the space-time defined by equation d s 2 = g α β d q α d q β , no central point r = 0 exists in the sense of a world-event on a valid world-line.
Next, we consider the motion of the medium far away from the event horizon, i.e., when parameter r is meaningfully greater than r g . For the flow at r > 3 r g , in the expansion of the metric tensor, we may neglect the terms of order ( r g / r ) 2 and greater. They contribute less than ( 1 / 3 ) 2 10 % to the components of the metric tensor. Such omission of smaller terms makes our approximation Newtonian or post-Newtonian. The non-diagonal metric-tensor term (describing involvement of the medium in the rotation of space-time in the vicinity of the black hole) gives rise to the “force” analogous to the traditional Coriolis force in the equations for medium flow.
Hence, we write the system of equations of relativistic fluid dynamics in the curved space-time and expand the metric tensor and the fluid energy–momentum tensor into a series with respect to small parameters r g / r < 1 and ∼ c 1 . We keep only the leading terms in the equations of fluid motion.
Next, the fluid is presumed to be localized near surface θ = π / 2 (in a pancake-like accretion disk)—the flows of the disk medium are considered only near this surface. Then, we can transition to cylindrical coordinates q i = ( r , ϕ , z ) and presume that the gas particles orbit near the z = 0 plane and the vertical component of their velocity v z m a x ( v r , v ϕ ) .
In this model, when r exceeds the Schwarzschild radius r g (specifically, r > 3 r g ), then the vertical component of gravity acceleration g z = G M z / ( r 2 + z 2 ) 3 / 2 K 2 z ( 1 + z 2 / r 2 ) 3 / 2 is balanced by the pressure gradient. Here, G is the gravitational constant, c is the speed of light, K 2 = G M / r 3 = c 2 r g / ( 2 r 3 ) , and K r 3 / 2 is the radius depending so-called Kepler parameter: the angular velocity of a test particle in a circular orbit at distance r from a point mass M in Newtonian approximation of gravity.
Gas pressure along the direction perpendicular to the disk plane ( x , y ) is determined by the hydrostatic equilibrium, d P = ρ g z d z . Generally speaking, pressure may have important contributions from electromagnetic radiation and induced magnetic field. The simplest case, however, is when the pressure is dominated by gas pressure, and the vertical temperature distribution is isothermal, which is roughly appropriate when the disk is optically thick and externally heated. The equation of state of gas/plasma is then P = s 2 ρ , where s is the “isothermal sound speed” (which is not a function of the transversal coordinate z).
If we further assume that z r , the equation of hydrostatic equilibrium becomes s 2 d ρ = s 2 h 2 ρ z d z , the solution of which is ρ ( z ) = ρ 0 exp ( z 2 / 2 h 2 ) . Here, the transversal space scale h is introduced, h = s r 3 / 2 / ( G M ) 1 / 2 . Parameter ρ 0 has the meaning of density at the disk mid-plane (at z = 0 ), and parameter h is the characteristic local thickness of the accretion disk.
In order for the accretion disk to be considered thin, it is necessary that h / r s r 1 / 2 / ( G M ) 1 / 2 ( s / c ) ( r / r g ) 1 / 2 1 . On the other hand, for great distances away from the black hole ( r > r g ), specific relativistic effects may be neglected or parametrized. Thus, we obtain natural bounds: r g < r < r g ( c / s ) 2 .
Equations of Motion: The equations of motion for an inviscid incompressible fluid express the laws of conservation for mass, momentum, and energy. The flow of gas may be considered incompressible (see, for example, [16]) when its velocity v satisfies condition v 2 s 2 and the characteristic time scale of the flow change T L / s where L is the characteristic spatial scale where the flow characteristics change substantially. In this case, the mass conservation law is expressed as d i v v = 0 , i.e., the law is not ρ = C o n s t , but t ρ + v j j ρ = 0 .
We assume that the domain of the disk where the vortex structure of interest is formed, is characterized by an approximately constant angular velocity Ω . In a coordinate system rotating with the angular velocity Ω = ( 0 , 0 , Ω ) , where x x 1 and y x 2 are the axes in the horizontal plane and z is vertically upwards, the basic equations of motion become
j v j = 0 ,
ρ ( t v i + v j j v i + 2 ϵ i j k Ω j v k ) = i P ρ i Φ ,
t θ + v j j θ = 0 .
Here, v i are the components of the velocity field, ρ is density, P is pressure, and θ is temperature. The potential Φ of the force field is Φ K 2 r 2 + ( 1 / 2 ) K 2 z 2 ( 1 / 2 ) | [ Ω , r ] | 2 . Furthermore, here, [ a , b ] is the cross product of vectors a and b ; and ϵ i j k is the alternating tensor ( ϵ 123 = 1 , zero for any two indices being equal, + 1 for any even number of permutations from ϵ 123 , and 1 for any odd number of permutations).
Due to the assumption of an incompressible medium (implying s 2 0 ), the equation of state ρ = ρ ( θ , P ) turns into the density expression, which only depends on the temperature (and not on the pressure). In fact, ρ = ( ρ / θ ) θ + ( ρ / P ) P ρ β θ . The coefficient of thermal expansion, β = ρ 1 ( ρ / θ ) , may be assumed constant and positive. (For gases, β = 1 / θ 0 .) Thus, we can set ρ = ρ 0 ( 1 β ( θ θ 0 ) ) , where subscript zero denotes the reference values. Due to the fact that β ( θ θ 0 ) is generally significantly less then one, one may neglect the density variations in all principal terms and hence replace ρ with the constant value ρ 0 , except in the “buoyancy” term, which is proportional to j Φ (see, for example, [16,17]).
Next, we apply the curl operator × to the linear momentum conservation Equation (2). This gives
[ β θ , d v d t + 2 [ Ω , v ] ] = β [ ( τ s + τ ) , Φ ] ,
where brackets symbolize cross-product. Since β | θ | 1 , term β [ θ , A ] is small with respect to [ , A ] at the horizontal scales typical for any hydrodynamical vector A and may be not taken into account.
Consider now temperature θ as θ = θ 0 + τ s + τ , where θ 0 is the (constant) baseline temperature of the accretion disk, quantity τ s is an axially symmetrical part of the temperature distribution that is not time-dependent, and τ = τ ( t , x j ) is the dynamical quantity related to the vortex structures in the fluid. When r = ( x , z ) , v j = ( v , w ) , then ω i = ϵ i j k j v k + 2 Ω i , and Equation (4) can be rewritten in the tensorial form as
( D t + w 3 ) ω i = ω j j v i + β ϵ i j k ( j ( τ s + τ ) ) k Φ ,
where D t = t + ( v · ) is the substantial derivative, v = ( v 1 , v 2 ) is the flow velocity, and = ( 1 , 2 ) is the gradient operator with components i in the ( x , y ) plane.
In the simplest model—in which the vortex structures are realized in a ”thin” flat sheet of an incompressible inviscid fluid (i.e., when h L , V 2 s 2 c 2 , R e 1 )—the z-component of velocity, w, vanishes and may be dropped in all formulas (see, for example, [18]). Note also that ϵ 3 j k j τ s k Φ 0 because of an axial symmetry of both τ s and Φ . Thus, the set of equations for the z component of vorticity ω 3 and for variation of temperature τ becomes
j v j = 0 ,
D t ω 3 = β ϵ 3 j k j τ k Φ ,
D t τ = v i i τ s ,
where indices j , k = 1 , 2 . Equation (6) permits the introduction of the stream function ψ :
v i = ϵ 3 i j j ψ ϵ i j j ψ ,
where tensor ϵ i j is the antisymmetric unit tensor of the second order, ϵ 12 = ϵ 21 , and diagonal components are zero. In this case, vorticity ω 3 = Δ ψ + 2 Ω 3 with Δ = 1 2 + 2 2 .
Temperature Stratification: To describe the effect of temperature stratification—i.e., to show how stationary temperature increases as the distance from the axis of rotation increases—we express (using r = | x | , x = ( x 1 , x 2 ) and Taylor series expansion) the background distribution of the temperature τ s in the form
τ s = α 2 r 2 , α > 0 ,
where α is the parameter characterizing the “rapidity” of increase in temperature with distance from the disk rotation axis. (Generally speaking, the question of what shape the temperature profile takes within a black hole’s accretion disk is an open one. Experimental measurements remain challenging, despite significant progress. See, for example, [19,20].) When the leading contribution to potential Φ comes from the centrifugal effects, we can express Φ = K 2 | x | 2 ( 1 / 2 ) [ Ω , x ] 2 ( 1 / 2 ) [ Ω , x ] 2 . Thus, when Ω is presumed constant in a band of r where the vortex hotspots are forming, the set of coupled nonlinear evolution equations becomes
t Δ ψ + ϵ i k i ψ k Δ ψ = β Ω 2 ϵ i k x i k τ ,
t τ + ϵ i k i ψ k τ = α ϵ i k x i k ψ .
For a general case, Δ ψ should be replaced: Δ ψ Δ ψ 2 Ω .
Equations (11) and (12) have a transparent physical meaning. Their left sides describe transport of dynamic quantities: vortex Δ ψ and temperature perturbation τ . Their right sides describe “sources” that generate the vortices and temperature perturbations.
In other words, the vortices are generated by the source (the right part of Equation (11) which is effective (non-zero) only when there exists an inhomogeneous gravity-like force field Φ Ω 2 with which temperature perturbation τ interacts via the equation of state (when β 0 ). The quantity ψ —which characterizes the vortex field in the fluid—is transported by the self-induced flow (the left part of Equation (11)). On the other hand, in Equation (12), this temperature perturbation τ is transported by the self-generated flow ( ψ 0 ); the temperature perturbation τ ( t , x ) is generated by the source (the right part of Equation (12)), which is non-zero only when there exists a spatial and time-independent temperature gradient (i.e., when τ s 0 ). The processes are interlinked because the “source” of one dynamic quantity depends on the complex combination involving another dynamic quantity. The set of Equations (11) and (12) shows that when temperature stratification is absent ( α = 0 ) and there is no disk rotation (disk Ω = 0 , i.e., centrifugal force is zero), then Equations (11) and (12) degenerate into the traditional equations for vortex evolution and transport in a two-dimensional ideal fluid.
Linear Approximation: Assuming that an excess of temperature τ above the basic level of temperature is not too large, in view of the link between fields τ and ψ via Equation (12), we consider a simple linear dependence between the excess of temperature τ and the vorticity ψ that generates this excess:
τ = C ψ .
Obviously, it follows from Equation (13) and from the meaning of quantities τ and ψ that the dimension of the coefficient of proportionality is [ C ] = θ L 2 T [ t e m p e r a t u r e ] × [ l e n g h ] 2 × [ t i m e ] . By substituting Equation (13) into Equation (12), we obtain
t Δ ψ + ϵ i k i ψ k Δ ψ = C β Ω 2 ϵ i k x i k ψ ,
C β α Ω 2 × C ( t ψ + ϵ i k i ψ k ψ ) = C β α Ω 2 × α ϵ i k x i k ψ .
Hence, we conclude that both the first and second equations in Equations (14) and (15) describe the evolution of the same physical quantity. The equations will be consistent when the following condition is imposed on the current function ψ :
( t + ϵ i k i ψ k ) ( Δ ψ R 2 ψ ) = 0 .
Here, parameter R 2 = C 2 ( β / α ) Ω 2 . The dimension of this quantity is ( θ 1 L 2 T 1 ) 2 × θ 1 L 2 × θ 1 × T 2 = L 2 ; i.e., R is a space scale factor. Quantity q = Δ ψ R 2 ψ describes the distribution of generalized vorticity. Its change in time and in space has to satisfy the evolution equation Equation (16).
The stream function ψ is found via the Green function approach. In the symbolic integral form, in the boundless space, it is
ψ ( x , t ) = ( Δ R 2 ) 1 ( x , x ) q ( x , t ) .
The subsequent calculation procedure is as follows: (i) The initial vorticity distribution is set; the stream function ψ is found from Equation (17); together with it, the non-linear evolution Equation (16) is numerically solved. (ii) The distribution of vorticity is set in the form of a macrostructure with petals (with constant vorticity inside) whose moving boundaries evolve according to Equation (16); i.e., we are considering a region bounded by some closed contour (with possibly a rather complex shape) such that quantity q ( x , t ) takes a constant value inside and zero outside. For stationary dynamical regimes, this can be accomplished analytically (details of the contour dynamics method and of the operator techniques can be found, for example, in Refs. [21,22]), as well as Refs. [18,23]); (iii) For a strongly localized vortex, quantity q ( x , t ) can be parameterized by the function in form F ( x x 0 ( t ) ) , i.e., the one with the center at coordinate x = x 0 ( t ) , which satisfies the equation of transport x ˙ 0 i ( t ) + ϵ i k i ψ = 0 , i.e., when the center of vortex moves according to x ˙ 0 ( t ) v s . [ x 0 ( t ) ] = 0 . (Here, the symbol “dot” signifies the derivative with respect to time).
Stationary Vortex Structures: Stationary vortex structures—the ones rotating with constant angular velocity ω —are simpler to consider in a rotating coordinate system where the structures appear immovable; i.e., when rotation direction is co-aligned with z axis, the derivative with respect to time becomes t = ω ϵ i k x i k . (The procedure is laid out, for example, in Ref. [24].) Indeed, when f = f ( ρ , ϕ ω t ), then the calculation relying on the properties of Jacobeans produces t f = ω ϕ f = ω ( f , ρ ) / ( ϕ , ρ ) = ω ( ( f , ρ ) / ( x 1 , x 2 ) ) ( ( x 1 , x 2 ) / ( ϕ , ρ ) ) = ω ( ϵ i k ρ 1 x i k f ) ( ρ ) . Then, Equation (12) may be rewritten as:
ϵ i k i ω 2 | x | 2 + ψ k τ + α 2 | x | 2 = 0 ,
which is satisfied by the ansatz
τ + α 2 | x | 2 = F ( ω 2 | x | 2 + ψ ) .
The explicit expression of function F can be found from the obvious fact that temperature-driven flow perturbations must vanish when temperature perturbations vanish themselves. The suitable expressions is function F ( u ) = α u / ω . Then, temperature fluctuations are expressed via the stream function
τ = α ω ψ .
Combining this with Equation (11) where t = ω ϵ i k x i k is taken into account, we obtain the second evolution equation: the one that describes rotation of a stationary vortex structure with angular velocity ω caused by the self-induced field of hydrodynamical velocity:
ϵ i k ( ω x i + i ψ ) k Δ ψ R 2 ψ = 0 .
Once Equation (20) is written out, parameter R 2 becomes R 2 = ( α β ) 1 ( ω / Ω ) 2 . Here, β is the coefficient of thermal expansion. For almost all physically realizable situations, β > 0 (the well-known exception is water in the temperature range between 0 and + 4 ° C). Parameter α —characterizing the “rapidity” of increase in temperature with distance from the disk rotation axis—can be either positive or negative. When α > 0 , the periphery of the accretion disk is heated more than the central zone; when α < 0 , the central zone of the disk is hotter than the periphery. Thermal length scale λ θ = ( | α | β ) 1 / 2 is determined by the background state of the disk in the framework of the model. Obviously, space scale parameter R characterizes the “rapidity” of the decrease in temperature (and stream function) with distance from the disk rotation axis. Below, we write R 2 R 2 ( Θ ( α ) Θ ( α ) ) , where Θ ( α ) is the Heaviside step-function (equal to unit for positive argument and zero for negative argument), and, leaving the old notations, R 2 = ( | α | β ) 1 ( ω / Ω ) 2 .

3. Results

We use the model described above to gain insights into the three-spot structure in the EHT image of Sgr A* accretion disk (Figure 1A). The bright spots are the zones with higher (relative to the base level) temperature τ and (since τ ψ ) with higher vorticity Δ ψ . In this framework, thermo-hydrodynamical vorticity-field structures may be conceptually analyzed with two limit approaches: the method of localized vortices and the method of thermo-vorticial spots (blobs).
Localized Vortices: The simplest consideration for a three-spot stationary thermo-vortex structure (symmetric with respect to the rotation axis of a pancake-like thin accretion disk) is to analytically treat the hot zones as narrowly localized formations modeled as 2D delta-functions at the limit case. Doing so would allow us to find the required conditions for the existence of the observed phenomenon and the resulting relationships between key characteristics of the structure. (This approach is not limited to only structures with three spots, but can be generalized to structures with any vortex-number N.) Hence, we write:
Δ ψ R 2 ( Θ ( α ) Θ ( α ) ) ψ = j = 1 3 κ R 2 δ ( 2 ) ( x x ( j ) ) ,
where κ characterizes the intensity of one localized vortex, and δ ( 2 ) ( x x ( j ) ) δ ( 1 ) ( x x ( j ) ) × δ ( 1 ) ( y y ( j ) ) is the two-dimensional Dirac function in the x y -disk plane. The appearance of the factor R 2 in the right side of Equation (22) is due to the following reasoning: the delta-function is a dimensional function. Its dimension is [ δ ( 2 ) ] = [ l e n g t h ] 2 for 2D space. Obviously, the argument of the delta-function must be dimensionless. (Indeed, there is no such thing as s i n of “one inch” or t a n of “ten gallons”.) Therefore, δ ( 2 ) δ ( 1 ) ( R 1 x R 1 x ( j ) ) × δ ( 1 ) ( R 1 y R 1 y ( j ) ) = R 2 δ ( 2 ) ( x x ( j ) ) . In view of the meaning of the delta function in the right part of Equation (22) and of the structure of the left part of Equation (22), the scale-factor R is the same R that was introduced above—the characteristic size of the vortex kernel. Although alternatives may exist, following the Occam’s Razor principle—the problem-solving principle of parsimony that “entities should not be multiplied beyond necessity”—we choose the simplest option. Indeed, among the options to use R or R multiplied by dimensionless ( κ / ω ) n , we choose the simple R.
We further note that x ( j ) and y ( j ) may be expressed as components of the complex quantity l × exp ( i 2 π ( j 1 ) / 3 ) in the plane x + i y (where l is the radius of a circumscribed circle).
When a vortex is modeled by a delta-function, a weak (logarithmic) singularity appears in the distribution of stream function ψ (which is the solution of Equation (22)) when α > 0 . If α < 0 , or a smoother distribution is adopted for α > 0 , any singularities in ψ function then disappear.
The solution to Equation (22) for unbounded space is
ψ = κ R 2 2 π j = 1 3 Z 0 ( | x x ( j ) | R ) .
Here, Z 0 ( ξ ) = J 0 ( ξ ) Θ ( α ) + K 0 ( ξ ) Θ ( α ) , J 0 ( ξ ) is the Bessel function of order n = 0 , K 0 ( ξ ) is the modified Bessel (McDonald) function of order n = 0 .
Recall that for small values of the argument, function J 0 ( ξ ) 1 , function K 0 ( ξ ) behaves logarithmically: K 0 ( ξ ) log ξ + 0.1156 . For large values, function J 0 ( ξ ) , as known, behaves as J 0 ( ξ ) 2 / π ξ cos [ ξ + π / 4 ] , function K 0 ( ξ ) behaves as K 0 ( ξ ) π / ( 2 ξ ) exp ( ξ ) . This means that for α > 0 and r R ; i.e., when ξ r / R 1 , function K 0 ( ξ ) exponentially quickly tends to zero, and therefore, the stream function ψ , and consequently the local vortex magnitude and temperature excess τ , all will also tend to zero. This is the consequence of the initial choice to model the vortices via delta-functions. Numerically, for some characteristic points, K 0 ( 1 ) 0.421 and K 0 ( 0.4569 ) 1 . The derivative of Z 0 with respect to argument is Z 0 ( ξ ) = Z 1 ( ξ ) .
By substituting Equation (22) into Equation (21), and then taking into account Equation (23), and by setting the terms with delta function and their derivatives as equal to zero, we find the set of conditions for the existence of the modeled vortex structure:
ϵ i k n = 1 3 ω x i ( n ) + m = 1 3 κ R 2 2 π R ( x i ( n ) x i ( m ) ) | z ( n m ) | Z 1 ( | z ( n m ) | R κ R 2 k δ ( z ( n ) ) = 0 .
In this expression, symbol “prime” in the second summation indicates that the effect of vortex self-action is excluded. When notation z ( j ) = x x ( j ) is used, then z ( i j ) = x ( i ) x ( j ) and | z ( 12 ) | = | z ( 23 ) | = | z ( 31 ) | = 3 l (where l is the radius of a circumscribed circle).
To satisfy Equation (24), the coefficients before every singularity must equal zero. Thus, the following condition for the parameters of the vortex configuration arises:
ω l 2 + 2 κ R 2 2 π R l 2 ( 1 cos 2 π / 3 ) l 3 Z 1 ( 3 l R ) = 0 .
Recall that R = ( | α | β ) 1 / 2 | ω / Ω ) | is a dimensional function of four dimensional parameters.
Equation (25) can be rewritten in a compact dimensionless form, where ξ = 3 ( l / R ) :
κ ω = 2 π 3 ξ Z 1 ( ξ )
which says that for any specific ξ , ratio κ / ω is uniquely defined. Figure 3 plots Equation (26).
For the presented model, the temperature excess τ follows from Equation (20) and Equation (23):
τ = 1 2 π β ( κ ω ) ( ω Ω ) 2 j = 1 3 Z 0 ( | x x ( j ) | R ) .
The characteristic temperature scale is defined thus as τ 0 = ( 2 π β ) 1 ( κ / ω ) ( ω / Ω ) 2 . When α < 0 (disk periphery is cooler), both “anticyclones” (temperatures excesses) and “cyclones” (temperature depressions) are possible.
Figure 1B (on the front page of the article) illustrates the temperature distribution for a system with three localized vortices, calculated using Equation (27), for a case when α > 0 (disk periphery is hotter). In Figure 1B, the unit for the color scale is the characteristic magnitude of temperature excess τ 0 = ( 2 π β ) 1 ( κ / ω ) ( ω / Ω ) 2 . For the x and y axes, the unit is R = λ θ ( ω / Ω ) ; ω is the angular velocity of rotation of the vortex structure as a whole; Ω is the characteristic angular velocity of rotation of the accretion disk; and λ θ = ( α β ) 1 / 2 . In the expression for λ θ , the parameter α (characterizing the “rapidity” of the increase in temperature with distance from the disk rotation axis) may be estimated as α ( Δ θ / R * 2 ) . Here, R * is the characteristic external radius of the active part of the accretion disk (obviously not equal to R) is presumed to be greater than the size of the macrostructure, and Δ θ may be estimated as the temperature difference across the disk (from the periphery to the center). In the expression for λ θ , the parameter β = ρ 1 ( ρ / θ ) is the coefficient of thermal expansion, which (for the model whose equation of state is approximated by the equation for an ideal gas) may be written as β θ 1 , where θ is the averaged temperature of the accretion disk.
Thermo-Vorticial Spots (“Blobs”): Another fruitful approach is to use the concept of thermo-vortical spots, for which the vorticity is constant inside domains bounded by movable boundaries (contours) and is zero outside. Then, the problem of vortex evolution becomes the problem of examination of the movement of the contours and determination of their final macro-configuration. Examples of such macro-configurations are depicted in Figure 4 and Figure 5 (discussed below).
Within this framework—called the contour dynamics method (CDM) (see Refs. [18,22,23,24,25,26,27])—the distribution of vorticity q for a vortex structure may be written as q = q 0 Θ ( x , y ) . Function Θ ( x , y ) is a two-dimensional Heaviside step-function equal to one inside the domain in consideration and zero outside. Parameter q 0 describes the vorticity inside the spot, which is presumed constant. Obviously, the dimension of this parameter is [ q 0 ] = [ t i m e ] 1 .
Using CDM (see the physical foundation and the details on how to perform calculations using CDM in Ref. [22]), we find the general expression for temperature distribution:
τ ( x , y ) = ( q 0 α R L 2 π ω ) D d s K 1 ( | z ^ z | μ ) μ | z ^ z | z ^ , s ( z ^ z * ) | z ^ z | .
Here, the integral is taken along the contour (whose shape is previously found by solving the problem of the spot configuration); z = ( x + i y ) / L is the complex variable of the x y -position, normalized by the characteristic space-scale L = ( R / 2 ) ( q 0 / ω ) 1 / 3 , which is linked to the size of the macro-configuration; parameter μ = R / L is dimensionless; and domain D in the complex plane z is a region filled with the vorticity q 0 . The domain is bounded by a closed contour D ; this boundary is expressed in the parametric form z = z ^ ( s ) ; parameter s is the contour arc length beginning from some initial point; and derivative z ^ , s z ^ / s is a unit vector that is obviously tangential to the contour D .
Based on the contour-integral Equation (28), the temperature distribution τ ( x ) / τ 0 is solved (and plotted in Figure 4B) for the contour boundary (black line) in Figure 4A. Shaded gray in Figure 4A is the constant vorticity distribution q ( x ) = q 0 .
The key parameters here are τ 0 = q 0 α R L / 2 π ω , R = λ θ ( Ω / ω ) , λ θ = ( α β ) 1 / 2 , β = θ 1 , ( q 0 / ω ) = 8 ( L / R ) 3 and α = ( 8 / π 2 ) ( Δ θ / R * 2 ) , where θ is the averaged temperature of the accretion disk, R * ( R ) is the characteristic radius of the accretion disk, and Δ θ is the temperature difference across the disk (from the periphery to the center). Thus, the characteristic level of the temperature excess in this limit case (i.e., in the model of thermo-vorticial spots, not in the model of localized thermo-vortices) can be estimated by the expression
τ 0 = 1 4 π β ( q 0 ω ) 4 / 3 ( ω Ω ) 2 .
The above-mentioned parameters in Equation (29) are the very parameters that must be measured experimentally to be able to comprehend the phenomena in images such as Figure 1A.
Equation (29) is functionally similar to Equation (27). The difference is in the index of power dependence of the ratio of the vortex intensity and its angular velocity. This is an insignificant difference considering the fact that the two models describe radically different limit cases.

4. Conclusions

Large-scale long-lived vortices are found in many types of hydrodynamic flows. Large vortices in turbulent flows are called coherent structures. They are observed, for example, in the planetary atmospheres and in the oceans. The horizontal scales of the vortices are much greater than the atmospheric or oceanic thicknesses. In the simplest geophysical context, examples of such structures are the Gulf Stream rings, the vortices shed from coastal currents, the cyclones and anti-cyclones, the Antarctic Polar Vortex, etc. A very well-known example is Jupiter’s Great Red Spot: a huge vortex plunged in the equatorial flow, which has persisted for more than three centuries; the presence of intense small-scale turbulence around it does not destroy it. This prompts the question: under what conditions do large-scale vortex structures form?
In traditional 3D space, the turbulent motion is usually considered to be homogeneous and isotropic. Common sense and the laws of thermodynamics show that it is very difficult to extract energy from a fully chaotic system, and only with some additional specific properties of such systems is this possible to realize. Homogeneous, isotropic turbulence, which does not possess any preferred directions or preferred scales, is extremely symmetric, giving birth to large-scale vortices; self-organization seems to be quasi-improbable in this case. It is evident that the breaking of some structural symmetry is one of the necessary conditions for the possibility of self-organization. It becomes clear that turbulent motion of fluid/gas with broken spherical symmetry can be a candidate for a system where self-organization of 2D turbulence into large-scale 2D vortex structures can take place.
Notably, quasi-2D flows in “thin” pancake-like accretion disks (see, e.g., [28,29] on the physics of accretion), where components of hydrodynamical flows that are perpendicular to the disk plane are strongly suppressed, can be a place where large-scale vortices can be self-organized. However, another condition is an insignificant influence of a dissipative process on large-scale motions in the system. As is known, the role of dissipation in a typical hydrodynamical process is characterized by the so-called number of Reynolds, which is determined as the ratio of magnitude of the inertial term to the dissipative term in the equation of fluid motion. For “smooth” flows described by the models of classical hydrodynamics, the introduction of such a number is not a problem; however, for flows of plasma, the determination of the predominant mechanism of dissipation is not apparent. However, the three-spot structure in the EHT image of Sgr A* accretion disk (seen in Figure 1A) is a clear example of self-organization in plasma.
To examine the observed phenomenon in the Sgr A* disk from the perspective of theoretical hydrodynamics, we first considered a simplified thermo-hydrodynamic model that permits analytical consideration. In the model, the vorticity clusters are approximated via delta-functions. As the result, we established the condition of existence of a regular thermo-vorticial structure (Equation (26); see also Figure 2). We also spelled out the relationships that should take place between, on the one hand, the parameters that determine the vortex-structure dynamics—each vortex size (l), its period of proper revolution ( 2 π ω 1 ), and temperature excess τ in the vortex—and, on the other hand, the accretion disk characteristics ( α and β , i.e., λ θ ) and the angular velocity of the entire accretion disk Ω .
The necessary conditions for the formation of large quasi-stationary symmetric thermo-vorticial structures in the plasma disk are as follows: (1) the accretion disk has to be pancake-like thin (i.e., h R * , where h is the thickness and R * is the characteristic radius) and rotate with non-zero angular velocity ( Ω 0 ); and (2) the disk temperature has to decrease towards the center (i.e., parameter α > 0 ). A multi-spot thermo-vortex structure forms only when key system parameters fall within the ranges captured by the dimensionless relationship: for a three-spot structure, κ / ω = ( 2 π / 3 ) ( ξ / K 1 ( ξ ) ) , where ξ = 3 ( l / λ θ ) ( Ω / ω ) . The temperature of the vortices τ (i.e., the excess over the base level) is also linked to another parameter of the system, vortex intensity κ of the hot spot: τ κ . Because the bright spots are hot, i.e., τ > 0 , this result also means that the rotation directions of the vortices ω and the entire accretion disk Ω must co-align.
In the framework of localized vortices, the estimate for temperature excess of the bright hot spots is given by Equation (27), and in the framework of thermo-vorticial spots (with sharp contour boundaries), by Equation (29).
For quasi-2D flows in thin layers of ideal fluid, another fruitful approach to analysis also exists (for deeper insight, see, Refs. [18,22,23,24,25,26,27], and applications and references therein). This approach is called the contour dynamics method (CDM). The gist of the method is that the continuous hydrodynamical velocity distribution may be treated as a set of patches with movable boundaries (contours) and constant vorticity inside. As the underpinnings of the CDM show, such vortex-patch approximation correctly grasps the general tendency in the dynamics and evolution of large-scale flows when the large-scale motions are weakly sensitive to a fine structure in the hydrodynamical velocity field. The contours of the vortex structures within the CDM framework are determined by spatially one-dimensional integro-differential nonlinear equations. Equations of the contour dynamics—which describe the self-induced motion of the vorticity-discontinuity boundaries, or “contours”, in an inviscid, incompressible, two-dimensional fluid with piecewise constant vorticity distribution—may be effectively resolved by either numerical or analytical approaches. Figure 5 illustrates the variability of shapes that a symmetric stationary macro-structure may take, depending on one of its guiding parameters. The shape may range from a weakly deformed circle to a sharply pronounced three-petal “flower”. (Other parameters define how many petals appear: two, three, etc.)
In view of the visual similarity between the hot-spot arrangement in the Sgr A* accretion disk revealed in Figure 1A and the thermo-vortex-structure plotted in Figure 5, we conclude that the observed hot spots in the Sgr A* accretion disk are highly likely to be large-scale quasi-2D quasi-stationary vortices in their nature.
In conclusion, let us emphasize that further improvement in understanding undoubtedly depends on the progress in numerical simulations, for which solid experimental knowledge of parameters defining the processes is paramount. Specifically, as discussed above, these parameters are the background temperature, its contrast, the size of the macro-structure, the number of petals, the angular velocity of rotation of the structure as a whole, the angular velocity of rotation of the accretion disk, gradients of velocity of the collective flows of plasma in the disk, and the characteristic time of existence of this quasi-stationary structure.

Author Contributions

Conceptualization and Writing, E.P.T., V.P.G. and V.I.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare that there is no conflict of interests regarding the publication of this article.

References

  1. Akiyama, K. et al. [Event Horizon Telescope Collaboration]. First Sagittarius A* Event Horizon Telescope Results. I. The Shadow of the Supermassive Black Hole in the Center of the Milky Way. Astrophys. J. Lett. 2022, 930, L12. [Google Scholar] [CrossRef]
  2. Akiyama, K. et al. [Event Horizon Telescope Collaboration]. First Sagittarius A* Event Horizon Telescope Results. II. EHT and Multiwavelength Observations, Data Processing, and Calibration. Astrophys. J. Lett. 2022, 930, L13. [Google Scholar] [CrossRef]
  3. Akiyama, K. et al. [Event Horizon Telescope Collaboration]. First Sagittarius A* Event Horizon Telescope Results. III. Imaging of the Galactic Center Supermassive Black Hole. Astrophys. J. Lett. 2022, 930, L14. [Google Scholar] [CrossRef]
  4. Akiyama, K. et al. [Event Horizon Telescope Collaboration]. First Sagittarius A* Event Horizon Telescope Results. IV. Variability, Morphology, and Black Hole Mass. Astrophys. J. Lett. 2022, 930, L15. [Google Scholar] [CrossRef]
  5. Akiyama, K. et al. [Event Horizon Telescope Collaboration]. First Sagittarius A* Event Horizon Telescope Results. V. Testing Astrophysical Models of the Galactic Center Black Hole. Astrophys. J. Lett. 2022, 930, L16. [Google Scholar] [CrossRef]
  6. Akiyama, K. et al. [Event Horizon Telescope Collaboration]. First Sagittarius A* Event Horizon Telescope Results. VI. Testing the Black Hole Metric. Astrophys. J. Lett. 2022, 930, L17. [Google Scholar] [CrossRef]
  7. Tito, E.P.; Pavlov, V.I. Relativistic Motion of Stars near Rotating Black Holes. Galaxies 2018, 6, 61. [Google Scholar] [CrossRef] [Green Version]
  8. Tito, E.P.; Pavlov, V.I. Black Hole Spin and Stellar Flyby Periastron Shift. Universe 2021, 7, 364. [Google Scholar] [CrossRef]
  9. Landau, L.D.; Lifshitz, E.M. The Classical Theory of Fields; Elsevier: Amsterdam, The Netherlands, 2013. [Google Scholar]
  10. Misner, C.W.; Thorne, K.S.; Wheeler, J.A. Gravitation; W. H. Freeman and Company, Princeton University Press: Princeton, NJ, USA, 1973. [Google Scholar]
  11. Shapiro, S.L.; Teukolsky, S.A. Black Holes, White Dwarfs, and Neutron Stars; Wiley: Hoboken, NJ, USA, 1983. [Google Scholar] [CrossRef]
  12. Visser, M. The Kerr Spacetime: A Brief Introduction. arXiv 2008, arXiv:0706.0622v3. [Google Scholar]
  13. Frolov, V.P.; Zelnikov, A. Introduction to Black Hole Physics; Oxford University Press: Oxford, UK, 2011. [Google Scholar]
  14. Weinberg, S. Gravitation and Cosmology: Principles and Applications of the General Theory of Relativity; Wiley: New York, NY, USA, 1972. [Google Scholar]
  15. Kadomtsev, B.B. Collective Phenomena in Plasma, 2nd ed.; USSR: Moscow, Russia, 1988. [Google Scholar]
  16. Landau, L.D.; Lifshitz, E.M. Fluid Mechanics, 2nd ed.; Pergamon Press: Oxford, NY, USA, 1987. [Google Scholar]
  17. Turner, J.S. Buoyancy Effects in Fluids; Cambridge University Press: Cambridge, UK, 1979. [Google Scholar]
  18. Goncharov, V.P.; Pavlov, V.I. Two-dimensional vortex motions of fluid at large Reynolds numbers. Phys. Fluids 1998, 10, 2384–2395. [Google Scholar] [CrossRef]
  19. Runge, J.; Walker, S.A. Probing within the Bondi radius of the ultramassive black hole in NGC 1600. arXiv 2021, arXiv:2102.06216. [Google Scholar] [CrossRef]
  20. Wong, K.-W.; Irwin, J.A.; Yukita, M.; Million, E.T.; Mathews, W.G.; Bregman, J.N. Resolving the Bondi accretion flow toward the supermassive black hole of NGC 3115 with Chandra. Astrophys. J. Lett. 2011, 736, L23. [Google Scholar] [CrossRef] [Green Version]
  21. Goncharov, V.P.; Pavlov, V.I. Problemy Gidrodinamiki v Gamil’tonovom Opisanii; Izd-vo Moskovskogo Universiteta: Moskva, Russia, 1993. [Google Scholar]
  22. Goncharov, V.P.; Pavlov, V.I. Hamiltonian Vortex and Wave Dynamics; Geos: Moscow, Russia, 2008. [Google Scholar]
  23. Goncharov, V.P.; Pavlov, V.I. Large-scale vortex structures in shear flow. Eur. J. Mech. B/Fluids 2000, 19, 831. [Google Scholar] [CrossRef]
  24. Goncharov, V.P.; Gryanik, V.M.; Pavlov, V.I. Venusian “hot spots”: Physical phenomenon and its quantification. Phys. Rev. E 2002, 66, 066304. [Google Scholar] [CrossRef]
  25. Goncharov, V.P.; Pavlov, V.I. Multipetal Vortex Structures in Two-Dimensional Models of Geophysical Fluid Dynamics and Plasma. Soviet J. Exp. Theor. Phys. 2001, 92, 594–607. [Google Scholar] [CrossRef]
  26. Goncharov, V.P.; Pavlov, V.I. Cyclostrophic vortices in polar regions of rotating planets. Nonlin. Process. Geophys. 2001, 8, 301–311. [Google Scholar] [CrossRef] [Green Version]
  27. Goncharov, V.P.; Pavlov, V.I. Hamitonian contour dynamics. In Fundamental and Applied Problems of the Vortex Theory; Borisov, A.V., Mamaev, I.S., Sokolovskiy, M.A., Eds.; Institute of Computer Science: Moscow/Izhevsk, Russia, 2003; pp. 179–237. (In Russian) [Google Scholar]
  28. Liu, B.F.; Qiao, E. Accretion around black holes: The geometry and spectra. iScience 2022, 25, 103544. [Google Scholar] [CrossRef] [PubMed]
  29. Armitage, P.J. Lecture notes on accretion disk physics. arXiv 2022, arXiv:2201.07262. [Google Scholar]
Figure 1. Left panel (A): Image of Sgr A* from Ref. [1]. Representative EHT image of Sgr A* from observations on 7 April 2017. This image is an average over different reconstruction methodologies (CLEAN, RML, and Bayesian) and reconstructed morphologies. Color denotes the specific intensity, shown in units of brightness temperature. The inset circle shows the restoring beam used for CLEAN image reconstructions ( 20 μ as FWHM). The bottom panels show average images within subsets with similar morphologies, with their prevalence indicated by the inset bars. Right panel (B): Normalized distribution of temperature-excess in an accretion disk for the model of localized vortices (Section 3).
Figure 1. Left panel (A): Image of Sgr A* from Ref. [1]. Representative EHT image of Sgr A* from observations on 7 April 2017. This image is an average over different reconstruction methodologies (CLEAN, RML, and Bayesian) and reconstructed morphologies. Color denotes the specific intensity, shown in units of brightness temperature. The inset circle shows the restoring beam used for CLEAN image reconstructions ( 20 μ as FWHM). The bottom panels show average images within subsets with similar morphologies, with their prevalence indicated by the inset bars. Right panel (B): Normalized distribution of temperature-excess in an accretion disk for the model of localized vortices (Section 3).
Universe 09 00040 g001
Figure 2. From Ref. [1] (one panel from original Figure 2). EHT baseline coverage, where dimensionless coordinates u = ( u , v ) give the projected baseline vector for each antenna pair in units of the observing wavelength.
Figure 2. From Ref. [1] (one panel from original Figure 2). EHT baseline coverage, where dimensionless coordinates u = ( u , v ) give the projected baseline vector for each antenna pair in units of the observing wavelength.
Universe 09 00040 g002
Figure 3. Model of Localized Vortices. Left panel (A): Both axes are in linear scale. Right panel (B): Vertical axis is in log-scale. Both panels (A,B): A structure with three vortices forms only when κ / ω = ( 2 π / 3 ) ( ξ / Z 1 ( ξ ) ) , where ξ = 3 ( l / R ) , l is the radius of a circumscribed circle, R = ( | α | β ) 1 / 2 | ω / Ω ) | , α is the parameter characterizing the “rapidity” of increase in temperature with distance from the disk rotation axis, β = ρ 1 ( ρ / θ ) is the coefficient of thermal expansion, ω is the angular velocity of rotation of the vortex structure, Ω is the angular velocity of rotation of the accretion disk, and κ is the characteristic intensity of one localized vortex. The blue curve is for α > 0 (the disk periphery is hotter). The red curve is for α < 0 (the disk periphery is cooler): intensity κ as ξ 3 ( l / R ) ξ N (red dot) at which the Bessel function J 1 ( ξ N ) = 0 . (Here ξ 1 3.8317 ).
Figure 3. Model of Localized Vortices. Left panel (A): Both axes are in linear scale. Right panel (B): Vertical axis is in log-scale. Both panels (A,B): A structure with three vortices forms only when κ / ω = ( 2 π / 3 ) ( ξ / Z 1 ( ξ ) ) , where ξ = 3 ( l / R ) , l is the radius of a circumscribed circle, R = ( | α | β ) 1 / 2 | ω / Ω ) | , α is the parameter characterizing the “rapidity” of increase in temperature with distance from the disk rotation axis, β = ρ 1 ( ρ / θ ) is the coefficient of thermal expansion, ω is the angular velocity of rotation of the vortex structure, Ω is the angular velocity of rotation of the accretion disk, and κ is the characteristic intensity of one localized vortex. The blue curve is for α > 0 (the disk periphery is hotter). The red curve is for α < 0 (the disk periphery is cooler): intensity κ as ξ 3 ( l / R ) ξ N (red dot) at which the Bessel function J 1 ( ξ N ) = 0 . (Here ξ 1 3.8317 ).
Universe 09 00040 g003
Figure 4. Illustration of analysis via method of thermo-vorticial spots (“blobs”). Left panel (A): A three-petal vorticity distribution is specified as q ( x ) = q 0 Θ ( x , y ) , where vorticity q 0 is constant and function Θ ( x , y ) is two-dimensional Heaviside step-function equal to one inside the shaded domain and zero outside. (Ref. [22] explains in detail the entire theory and methodology. This depicted three-petal thermo-vorticial structure takes this particular shape when one of the guiding parameters of the macro-configuration reaches its limit case.) Right Panel (B): For the vorticity distribution specified in Panel (A), the resulting temperature distribution τ ( x ) is obtained per Equation (28). The unit along the vertical axis is the characteristic magnitude of temperature excess τ 0 = q 0 α R L / 2 π ω . (The central τ peak is a consequence of the vorticity parametrization for the geophysical application in Ref. [22], which for a black hole accretion disk should be obviously nil.) Both Panels (A,B): For x and y axes, the unit is R = λ θ ( ω / Ω ) ; ω is the angular velocity of rotation of the vortex structure as a whole; Ω is the characteristic angular velocity of rotation of the accretion disk; and λ θ = ( α β ) 1 / 2 . In the expression for λ θ , the parameter α (characterizing the “rapidity” of increase in temperature with distance from the disk rotation axis) may be estimated as α = ( 8 / π 2 ) ( Δ θ / R * 2 ) . Here R * is the characteristic radius of the accretion disk (obviously not equal to R) is presumed to be greater than size of the macro-structure, and Δ θ may be estimated as the temperature difference across the disk (from the periphery to the center). In the expression for λ θ , the parameter β = ρ 1 ( ρ / θ ) is the coefficient of thermal expansion, which (for the model whose equation of state is approximated by the equation for an ideal gas) may be written as β θ 1 , where θ is the averaged temperature of the accretion disk.
Figure 4. Illustration of analysis via method of thermo-vorticial spots (“blobs”). Left panel (A): A three-petal vorticity distribution is specified as q ( x ) = q 0 Θ ( x , y ) , where vorticity q 0 is constant and function Θ ( x , y ) is two-dimensional Heaviside step-function equal to one inside the shaded domain and zero outside. (Ref. [22] explains in detail the entire theory and methodology. This depicted three-petal thermo-vorticial structure takes this particular shape when one of the guiding parameters of the macro-configuration reaches its limit case.) Right Panel (B): For the vorticity distribution specified in Panel (A), the resulting temperature distribution τ ( x ) is obtained per Equation (28). The unit along the vertical axis is the characteristic magnitude of temperature excess τ 0 = q 0 α R L / 2 π ω . (The central τ peak is a consequence of the vorticity parametrization for the geophysical application in Ref. [22], which for a black hole accretion disk should be obviously nil.) Both Panels (A,B): For x and y axes, the unit is R = λ θ ( ω / Ω ) ; ω is the angular velocity of rotation of the vortex structure as a whole; Ω is the characteristic angular velocity of rotation of the accretion disk; and λ θ = ( α β ) 1 / 2 . In the expression for λ θ , the parameter α (characterizing the “rapidity” of increase in temperature with distance from the disk rotation axis) may be estimated as α = ( 8 / π 2 ) ( Δ θ / R * 2 ) . Here R * is the characteristic radius of the accretion disk (obviously not equal to R) is presumed to be greater than size of the macro-structure, and Δ θ may be estimated as the temperature difference across the disk (from the periphery to the center). In the expression for λ θ , the parameter β = ρ 1 ( ρ / θ ) is the coefficient of thermal expansion, which (for the model whose equation of state is approximated by the equation for an ideal gas) may be written as β θ 1 , where θ is the averaged temperature of the accretion disk.
Universe 09 00040 g004
Figure 5. The four panels (obtained in Ref. [22] via the method of thermo-vorticial spots) illustrate that uniformly rotating macro-structures may take various shapes depending on guiding parameters: the shape may range from a weakly deformed circle to a sharply pronounced three-petal “flower”. Other guiding parameters define how many petals appear: two, three, etc. The theory and models are explained in detail in Ref. [22]. The units for axes are the same as in Figure 4.
Figure 5. The four panels (obtained in Ref. [22] via the method of thermo-vorticial spots) illustrate that uniformly rotating macro-structures may take various shapes depending on guiding parameters: the shape may range from a weakly deformed circle to a sharply pronounced three-petal “flower”. Other guiding parameters define how many petals appear: two, three, etc. The theory and models are explained in detail in Ref. [22]. The units for axes are the same as in Figure 4.
Universe 09 00040 g005
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Tito, E.P.; Goncharov, V.P.; Pavlov, V.I. Hot Spots in Sgr A* Accretion Disk: Hydrodynamic Insights. Universe 2023, 9, 40. https://doi.org/10.3390/universe9010040

AMA Style

Tito EP, Goncharov VP, Pavlov VI. Hot Spots in Sgr A* Accretion Disk: Hydrodynamic Insights. Universe. 2023; 9(1):40. https://doi.org/10.3390/universe9010040

Chicago/Turabian Style

Tito, Elizabeth P., Victor P. Goncharov, and Vadim I. Pavlov. 2023. "Hot Spots in Sgr A* Accretion Disk: Hydrodynamic Insights" Universe 9, no. 1: 40. https://doi.org/10.3390/universe9010040

APA Style

Tito, E. P., Goncharov, V. P., & Pavlov, V. I. (2023). Hot Spots in Sgr A* Accretion Disk: Hydrodynamic Insights. Universe, 9(1), 40. https://doi.org/10.3390/universe9010040

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