Modeling the Evolution of a gravitating bodies cluster based on absolutely inelastic collisions Kiryan D.G., Kiryan G.V. Institute of Problems of Mechanical Engineering of RAS 61 Bolshoy Prospect V.O., 199178, Saint Petersburg, Russia e-mail: diki.ipme@gmail.com arxiv:1811.06372v1 [physics.class-ph] 14 Nov 2018 Numerical simulation of evolution of a cluster of a finite number of gravitating bodies has been accomplished in the scope of classical mechanics taking into account accretion. The goal of the study was to reveal the basic characteristic phases of the intra-cluster distribution of material bodies. In solving the problem, the possibility of interbody collisions was taken into account. The collisions were assumed to be absolutely inelastic. Non-gravitational forces external with respect to the body cluster in question were ignored. Among all the internal force factors acting within the cluster, only the gravitational interaction was taken into account. To check the process of solution, the so-called rotation curve was used which presents a current radial distribution of orbital velocities of the cluster bodies. The Cauchy problem was considered. The issues of defining natural initial characteristics of the cluster bodies were touched upon. Conditions for commencement of rotation of gravitating bodies comprising the cluster about their common instantaneous center of mass were investigated. The numerical analysis showed that the characteristic shape of the rotation curves of stars of some galaxies depends only on the current configuration of the material body orbits. The rotation curve plateau characterizes the current redistribution phase of the intra-cluster matter. This means that invariance of radial distribution of star linear velocities in some of the observed clusters can be explained without considering the hypothesis of the non-material gravitating dark matter or modifying the classical Newton s Law on gravitational interaction between two material bodies. 1 Objective To show via a numerical experiment in the scope of classical mechanics that the characteristic shape of the rotation curve 1 of a cluster of gravitating bodies can be explained without applying the dark matter hypothesis. 1 This is the curve representing orbital velocities of the galaxy stars versus the radius of rotation about the galaxy center of gravity. Just this specific feature of the curve, namely, invariance of the star orbital velocities with distance from the conditional center of the galaxy, gave rise to such a non-material essence as dark matter and also initiated attempts to modify the classical Newton s Law on the gravitational interaction between two material bodies. 1
2 Problem definition The classical law on gravitational interaction between material bodies states that the force of gravitational attraction of two homogeneous spheres (material points) is directly proportional to their gravitating masses and inversely proportional to the squared distance between their centers of symmetry (centers of mass). Let us consider dynamics of a closed system consisting of n bodies (homogeneous spheres) with masses m i,i = 1,...,n taking into consideration only gravitational accretion. No kinematic restrictions are imposed on the cluster components. External force factors are excluded. Internal force interactions between bodies are limited to only the gravitational interaction. Non-gravitational processes are ignored. The body-to-body collisions are assumed to be absolutely inelastic. Evolution of the cluster of gravitational bodies is a process of continuous gravitational interaction between the bodies leading to spatial rearrangement of the cluster structure and also to reduction of the total number of bodies due to absolute inelasticity of collisions. Thus, the problem will be defined as follows: computer simulation of evolution of a cluster of gravitating bodies based on only gravitational accretion. Let us construct a fixed Cartesian frame of reference Oxyz (Fig. 1) m j z m j+1 r j r j+1 f j f j+1 m i O r i y x r j 1 m j 1 f j 1 Figure 1: Schematic diagram of a gravitational impact upon a cluster body with mass m i from bodies with masses m j 1, m j and m j+1. Spatial location and velocity of each i-th body of the cluster are defined by vectors r i andv i, respectively. Each body interacts with all others. Due to additivity, we can sum up gravitational forces acting upon body m i from all 2
the bodies of the system. The gravitational force acting upon body m i from body m j is defined as a product of the body m j gravitational field intensity at the point where body m i is located and mass m i. Designate as f the sum of forces of other than gravitational nature acting on moving body m i from its material environment. This has been done only for generalization. In the scope of our task those forces are excluded from consideration. Let us construct a set of second-order differential equations modeling dynamics of gravitational interaction between n bodies comprising the cluster: m i d 2 r i dt 2 = n j=1, j i f j {}}{ 0 E(m j,r j r i ) m i + f (t,r i,ṙ i ), i = 1,...,n. (1) Here E is the gravitation field intensity of the j-th body at distance r j r i. Each gravitating body of the cluster in question possesses its own gravitation field characterized by field intensity vector E. E(m j,r j r i ) = G m j r j r i (r j r 2 i ), (2) where G is the coefficient matching scales and dimensions (scale-dimension factor) 2. In solving equation set (1), we will control the bodies approach to each other to the critical distance equal to the semi-sum of body diameters which determines the moment of the absolutely inelastic collision. Let us assume that the new body formed after the two-body collision has a diameter equal to that of the largest of the two bodies and mass equal to the sum of their masses. To our opinion, the absolutely inelastic collision is that the new body formed as a result of the contact of two bodies continues moving in the cluster gravitational field with the velocity dictated by the law of momentum conservation. For instance, in the case of collision of two bodies with masses m 1, m 2 and velocities v 1, v 2, respectively, velocity of the stuck bodies may be represented as follows: v = m 1 v 1 +m 2 v 2 (m 1 +m 2 ). (4) 2 At present, international Committee on Data for Science and Technology (CODATA) recommends the following value of the gravitational constant [1] G = 6.67384(80) 10 11 m 3 /(kg s 2 ). (3) Note that the accuracy of the given value gives rise to some doubts in its reliability. In reality, it is possible to speak about only two decimals in the SI system. This is described in more details in paper [2]. 3
Now, to solve the set of second-order differential equations (1), we need to define natural initial conditions. The next section is devoted to this problem. 3 Initial conditions Well, the study object is defined: this is a cluster of bodies interacting with each other only by gravity. Now consider the initial conditions. What spatial distribution of bodies should be chosen? What should be the preset values of velocities? Since we are going to simulate really observable evolution phases of the gravitating bodies (e.g., galaxies), the initial conditions should comply with the real system state at the chosen time moment. This is just what we call natural initial conditions. Let us superpose the cluster center of mass with point O that is the origin of reference frame Oxyz (Fig. 1). Assume that all the bodies of the cluster under consideration move counterclockwise about point O in the xoy plane. At some moment of the cluster evolution, when gravity forces get prevailing in it, there begins ordered rotation of bodies about their common center of mass. Thus, self-rotation of the gravitating bodies cluster about the instantaneous center of mass in any direction results from the combination of the bodies mutual attraction and curvilinear motion of the cluster as a whole in the external initially inhomogeneous gravitation field. Consider in more details the mechanism for arising of self-rotation of a two-body cluster (Fig. 2). Let two bodies m 1 and m 2 forming a cluster move along the same circular trajectory. This condition is not mandatory, but let us accept it for clarity. The trajectory circularity is caused by an external gravitating mass located at point O. Under the condition of joint curvilinear motion of bodies, their gravitational interaction initiates rotation about their common center of mass. Mutual attraction of bodies moving jointly counterclockwise along a curvilinear (circular) trajectory decelerates body m 1 and accelerates body m 2. As a result, the balance of gravitational f g and centrifugal f ω forces gets violated for each body. Increase in the body m 2 velocity results in lifting of its orbit, deceleration of body m 1 velocity causes its orbit lowering. Hence, the two-body system will begin rotating counterclockwise about the common center of mass that, in its turn, moves along a circular trajectory. Body trajectories presented in Fig. 2 are designated as ξ 1 and ξ 2. The clockwise rotation of the cluster bodies arises similarly. Therefore, it is not important what direction of rotation of bodies about the cluster s instantaneous center of mass has been chosen, the main point 4
ξ 1 v 1 f ω 1 f g 1 f 1 m 1 ξ 2 f 21 v 2 y r 1(t) f 12 f 2 f ω 2 f g 2 λ 1(t) λ 2(t) r 2(t) m 2 O x Figure 2: Conditions for self-rotation of a two-body system (cluster). is that rotation must exist. Let it be counterclockwise. In this study we consider that stage of the total matter evolution [3] in the Space which is caused only by the gravitational interaction between material bodies. One of specific features of the gravitating systems (galaxies) is a plateau in the rotation curve 3, that takes place when linear velocities of the bodies become invariant with body distance from the cluster center of mass. Fig. 3 presents in one and the same scale rotation curves of different galaxies containing the characteristic plateau. Earlier paper [5] showed that the rotation curve plateau (Fig. 3) merely reflects current evolution moments of some galaxies. Not going into details, we can assume that evolution of gravitating body clusters proceeds in three main stages: 1. Initial stage. 3D distribution of the gravitating bodies resulting from accretions of various natures. 2. Borderline stage. Well-pronounced ordered rotation of the cluster bodies about its center of mass. The presence of a plateau in the cluster rotation curve. 3 Based on systematic observations of the 21-th spiral galaxy (i.e., measurements of Doppler shifts of star spectral lines) [4], V.C. Rubin has obtained a characteristic radial distribution of the orbital velocity with the plateau-like section. 5
v, km/s 250 200 150 100 4 6 3 1 5 2 1. NGC6503 2. NGC3198 3. NGC2403 4. NGC598 (M33) 5. NGC7331 6. Milky Way 50 0 10000 30000 50000 70000 100000 r, year Figure 3: Examples of rotation curves of different galaxies taken from [5]. v is the radial velocity; r is the distance in time units. 3. Final stage. Minimal manifestation of the gravitational accretion. Based on the above, we should take as initial conditions for the problem defined (1) the cluster body coordinates and velocities at the moment preceding formation of the rotation curve plateau. Assume that initially all the cluster bodies move along circular orbits. But, what is mutual location of the body orbits, or, more exactly, what it should be? What should be the preset values of velocities? Our task is not only to calculate dynamics of the cluster bodies taking into account gravitational accretion but also to obtain in the process of the numerical experiment the rotation curve with the plateau and compare it with Doppler measurements of star radial velocities for some galaxies. From general considerations, let us use the following empiric formula to preset initial distribution of orbits radii: r i = r ( 3 i ) α, i = 1,...n, (5) where r is the minimal radius of the cluster body circular orbit, α is the parameter defining the character of the circular orbits distribution, i is the orbit number. Based on the orbit radii sequence (5) whose character is defined by parameter α, let us construct an ordered sequence of nested Spheres [5]. The mass of each Sphere is determined by the total mass of all the bodies whose orbit radii are shorter than the Sphere radius. Define the mass and density of the i-th Sphere as follows: i 1 m s i = m j, ρ s i = ms i j=1 / 4 3 πr3 i, i = 2,...,n, (6) 6
where m s i and ρ s i are the mass and density of the i-th Sphere, respectively. Superscript s means that we consider an ordered sequence of nested Spheres. Now consider the velocities. Each Sphere of the sequence has a unique characteristic that is velocity of a test body moving along a circular orbit whose radius is equal to the Sphere radius. Now let us write down the expression for the i-th body circular velocity: vi s = G ms i, i = 2,...,n, (7) r i where vi s is the circular velocity of the i-th body moving along the orbit with radius r i, m s i is the i-th Sphere mass. Fig. 4 demonstrates the variants of distribution of the cluster bodies orbits calculated via formula (5) for different values of α. To each ordered distribution of orbits m s mass α 1.3 v s velocity ρ s density β 0.52 α 1.9 β 0.60 α 2.7 β 0.71 α 3.5 β 0.76 α 7.0 β 0.89 radius, r radius, r radius, r radius, r Figure 4: Variants of the natural initial values of body orbit radii at different stages of the cluster evolution. β is the cluster evolution number. orbit sequence, an individual sequence of nested Spheres is assigned, as well as sequences of their masses m s, circular velocities v s and densities ρ s. 7
Then let us choose based on Fig. 4 data such an orbit distribution defined by parameter α that precedes formation of the rotation curve plateau. In our case α<2.7. Thus, we have obtained a procedure for defining initial conditions for the task (1). 4 Numerical simulation The computational model of the gravitating bodies cluster is n homogeneous spheres of one and the same radius r and density ρ. m i = m = 100kg, ρ i = ρ = 2500kg/m 3, i = 1,...,n = 100. (8) Here we solve (1) for the two-dimensional case, which means that all the body trajectories lie in the xoy plane. The i-th body coordinates are defined by radius r i and angle λ i. The angular coordinate is measured counterclockwise from axis Ox. Angle λ i will be chosen for each body randomly within the range of 0 to 2π. Initial distribution of the orbit radii will be found via formula (5) with parameters α = 0.97 and r = 500m. Using formula (7), calculate initial circular velocities v i of the cluster bodies. Now, as initial conditions for the set of second-order differential equations (1) are defined, let us solve the Cauchy problem continuously checking the inter-body distance in order to find out the moment when the collision conditions are fulfilled. Two bodies (spheres) will be regarded as collided if the distance between their centers of mass is shorter than or equal to the semi-sum of their diameters. Since we consider absolutely inelastic collisions, the two bodies continue moving after contacting with the same velocity and in the same direction. Fig. 5 presents design trajectories of the bodies before and after the absolutely inelastic collision. The result of simulating evolution of a cluster of gravitating bodies (Fig. 6) looks like a hodgepodge of trajectories but only at first sight. How can we quantitatively estimate the obtained result from the evolution point of view? For this purpose, let us calculate the cluster evolution number β introduced in paper [5]. Using β, it is possible to estimate the current evolution stage of the galaxy (a cluster of gravitating bodies). Evolution number β anges from 0 (the initial evolution phase) to 1 (the final evolution phase) and is invariant with respect to the cluster size and masses of bodies comprising it. Value of β will be calculated as follows. First construct an ordered sequence of nested Spheres using already known body coordinates and masses, 8
Figure 5: Calculated trajectories of the bodies before and after an absolutely inelastic collision. The bodies move jointly counterclockwise about the instantaneous cluster s center of mass. Figure 6: Examples of calculated trajectories of the cluster gravitating bodies. i.e., form a sequence of m s according to formula (6). Approximate the calculated distribution of the Spheres masses with a general-type power function. This allows maximal reduction of the evolution duration and initial number of the cluster bodies. Knowing the Spheres masses and radii, construct the sequence of densities ρ s. Using the least-square method, approximate the obtained sequence ρ s with the power function (9) thus determining evolution number β. ρ(r) = ar 3β +ρ 0, r > 0, ρ 0 0, 0 < β < 1. (9) 9
Here r is the Sphere radius, ρ is the Sphere density; β is the dimensionless coefficient, namely, evolution number, ρ 0 is the density of medium containing the cluster of gravitating bodies, a is the coefficient matching scales and dimensions (scale-dimension factor). The result of modeling dynamics of the cluster bodies for different evolution durations T = {1200; 1500; 2500; 3000; 5000} days is given in Fig. 7 in the form of a sequence of 5 evolution phases. Red color in the figure corresponds to the initial state of the system and its characteristics while blue color indicates the system state at moment T. The first column (top to bottom) reflects variations in the distribution of body masses and orbits over the duration of cluster evolution. Digits in the top-right corners of the first column plots present the ratios between the initial number of bodies and current one. The second column demonstrates dynamics of the nested Spheres masses. The number of Spheres is equal to that of gravitating bodies, while their radii depend on the current distance from the reference frame origin Oxyz. The sequence of nested Spheres gets formed after sorting the Spheres radii in the ascending order. Thus, we construct an ordered sequence of Spheres for each time moment. The mass of each Sphere is a sum of masses of all the bodies included into this Sphere. All the Spheres are defined in the Cartesian frame of references Oxyz and their centers coincide with point O. Blue circles represent the design distribution of the nested Spheres masses, while the solid blue line is its approximation with a power function. Quality of the sequence m s approximation is characterized by correlation coefficient R whose value allows us to use later the obtained power function with certain parameters. The third column demonstrates variations in the orbital velocity v s distribution (7) in the process of evolution of the gravitating bodies cluster and in evolution number β. 5 Conclusions 1. Numerical simulation of evolution of a gravitating bodies cluster taking into account accretions (absolutely inelastic collisions) has shown that mysterious plateaus in rotation curves of the observed galaxies characterize instantaneous distributions of the classical gravitating matter. 2. Evolution number β enables quantitative estimation of the current evolution phase of the gravitating bodies cluster (galaxy). 10
distribution of mass T =1200; 100/48 m s mass v s velocity R = 0.961 β = 0.27 T =1500; 100/43 R = 0.964 β = 0.44 T =2500; 100/36 R = 0.959 β = 0.67 T =3000; 100/31 R = 0.938 β = 0.80 T =5000; 100/24 r, radius of orbit R = 0.952 r, radius of orbit β = 0.92 r, radius of orbit Figure 7: Calculated evolution phases of the gravitating bodies cluster. Red color indicates the initial state of the cluster bodies, blue color is for the current state. T is the evolution duration in days, R is the correlation coefficient, β is the evolution number, 100/48 is the ratio of the initial number of bodies to the final one. 3. The paper has shown the efficiency of analyzing the current state of the gravitating bodies cluster by using the method of nested Spheres [5]. 11
6 Afterword Here we will briefly comment the rotation curve phenomenon based on the above. Consider the Fig. 8 curves. Here we see the galaxy rotation curve Figure 8: Characteristic distribution of the observed and calculated linear star velocities over their distances from the galaxy NGC3198 center [6]. and its theoretically expected Kepler s rotation curve. Evidently, they do not coincide. What can we conclude from this? The rotation curve represents Doppler measurements of radial velocities of the galaxy stars. The rotation curve plateau is an instrumentally observed fact. Each star is held at an almost circular orbit by the total mass of the matter enclosed in a Sphere whose radius is equal to that of the star orbit. The conduced numerical experiment showed that it is not vitally necessary to use the dark matter hypothesis or try to modify the Newton s Universal Gravitation Law in order to explain existence of the galaxy rotation curve plateaus (Figs. 3, 8). All can be explained in the scope of classical mechanics. To this we can add that as early as in 19-th century paper [7] clearly showed invalidity of the dark matter hypothesis. References [1] P. J. Mohr, B. N. Taylor, and D. B. Newell, CODATA recommended values of the fundamental physical constants: 2010, Reviews of Modern Physics, vol. 84, pp. 1527 1605, Nov. 2012. 12
[2] D. G. Kiryan and G. V. Kiryan, Decimal-by-Decimal Analysis of the Gravitational Constant Value as Exemplified by Torsion Balance, International Journal of Applied and Fundamental Research, vol. 4, p. 12, Nov. 2018. [3] H. Alfvén and G. Arrhenius, Evolution of the Solar System. NASA SP, Scientific and Technical Information Office, National Aeronautics and Space Administration : for sale by the Supt. of Docs., U.S. Govt. Print. Off., 1976. [4] V. C. Rubin, W. K. J. Ford, and N.. Thonnard, Rotational properties of 21 SC galaxies with a large range of luminosities and radii, from NGC 4605 /R=4kpc/ to UGC 2885 /R=122 kpc/, The Astrophysical Journal, vol. 238, pp. 471 487, Jun. 1980. [5] D. G. Kiryan and G. V. Kiryan, The Evolution of the System of Gravitating Bodies, International Journal of Applied and Fundamental Research, vol. 4, p. 19, Dec. 2017. [6] Y. Yang and W. B. Yeung, A new approach to the yang-mills gauge theory of gravity and its applications, ArXiv e-prints, 2013. [7] H. C. Andersen, Kejsarens nya kläder. 1837. 13