Regular n-gon as a model of discrete gravitational system Rosaev A.E. OAO NPC NEDRA, Jaroslavl Russia, E-mail: hegem@mail.ru
Introduction A system of N points, each having mass m, forming a planar regular polygon (N-gon), and a central mass M, are considered. The motion equations for a testing particle are given in different coordinates. To provide stability, N-gon must always rotate, including case zero central mass. The velocity of rotation is calculated. N stationary points (libration points) appear in a system. Studying of libration point behavior is a main result of this paper. It is shown, that the stationary solution (libration points) in considered system may be determined from algebraic equation of 5-th degree. In case small m / M coordinates of these points may be calculated as a generalization of classical gravitation -body problem. In present research, the dependence of libration points coordinates on mass and number of particles is studied. Obtained results are discussed and compared with another authors. Area, where potential allows a linearization is described in small vicinity of N-gon. However, oscillations even in this small area become strongly nonlinear due to differential rotation of system. In fact, approximation ring potential by N-gon is an example of a pointmass modeling. But specific of potential expansion gives an ability to study more complex (non-homogenous) systems on base of considered model.
Definitions Definition [1]: A configuration of the n particles is a central configuration at time t if there exist some scalar such that r r for i=1,...,n. U U mir i mri 0 r r i i i i In particular every body can move on a circle around the common center of masses. In these solutions the motion behaves as if the bodies form a rigid body. These solutions are also denoted as relative equilibrium solutions, being a fixed point of (1) if we use a rotating frame. U mim r r i j j
Main target to describe motion infinitesimal particle perturbed by n-gon Determine stationary points in the system Show difference between three body and continue case Consider some generalizations
Libration points
A system of N points, each having mass m, placed in vertex of regular polygon, and a central mass M, is considered. Such system, forming a planar central configuration, called a relative equilibrium system [1]. There are N-1 stationary points (libration points) appear in a system (Fig.1). Fig. 1. Libration points positions for -body (left) and N-body (right) relative equilibrium configurations.
[]: The libration point coordinates can be determined from equations R GM R N j1 Gm x R sin j x R sin 1 x R j R GM R Gm x N 1 j 1 Gm x R sin j x R sin j 1 x R j N j 0 Here G gravity constant, R is the central configuration s radius, j is the angle between particles, and x is a distance between the test particle and the central configuration. In the first case 0 = /N, x=0 and we have non-collinear libration point. In the second case 0 = 0, x 0 and the x-coordinate of the collinear points must be determined.
Regular n-gon configuration begets a periodic solution in which the bodies rotate uniformly about the central mass with rotation speed (Roberts[6], Rosaev[]): N 1 1 0 1 sin j 4M j1 m 1/ 0 GM R The equations for the collinear libration points are the most interesting. They can be reduced to a fifth degree polynomial: ( +A)x 5 + ( R + AR + B)x 4 + ( R + R A+ BR)x + (BR Gm)x GmR =0 The sign + is for libration point inside central configuration, and for outside ones. Equations for non-collinear libration points have a similar form.
Coefficients are: Gm A R N 1 j1 1 sin 8sin j j B Gm 4R 1 N 1 1 j sin j When A=0 and B=0 we have well known case -body problem (Szebehely1967). Considered equation for determination libration points coordinates always have only one real root.
Summation Coefficients A of expansion depending on number of particles. For collinear case: N N 1/(sin( i / N)) N /( ) / i.404101... N /( ) i1 i1 N 100 1/(sin( i / N)).408618 / i.40401 500.40458.40411 1000.40418.40411 1500.404146.40411 000.4041.404114 500.40416.404114 000.4041.404114 For non-collinear case:: Gm A R N 1 j1 1 sin 1.654117... N j N 8sin j N /( ) Gm R N1 N Gm 1 Gm B N /( ) 1/ i 4R sin R j1 j N i1
By analogy with three body problem, solution of equations for collinear libration points can be calculated as a series:... 9 1 1 / 1/ 1 B A M m B A M m B A M m x... 9 1 1 / 1/ B A M m B A M m B A M m x for inner point L 1 and for outer point L respectively. These equations are valid for small m/m ratio and when N is not large than few hundreds. But in all cases, solutions of fifth degree equation can be obtained for example, with Maple.
F := 0.80881970885 10 14 x 5 0.7661775810708 10 14 x 4 809057078071 x 0.9861775186 10 7 x.9559 m/m=0.1000000, ln N=5 x=-0.00004699 m/m=0.1000000, ln N=5 x=0.000046875
Libration point coordinates on mass ratio dependence
Fig. To a non-collinear points coordinates determination Equation for non-collinear libration points has trivial solution x=0 and non-trivial solution of cubic equation: ( ( A) x R ( R R AR B) x A BR) x BR which have at least one real root. For case B=0 we have three real roots if A sufficiently large (small central mass): 0 J := 0, A 4 4 A, ( A) A 4 4 A ( A)
Motion in field of N-gon
In restricted three body problem (RTBP), each k-th of it s period, test particle has strong perturbation close encounter with planet. For case N+ body problem strong perturbed encounters take place in N-times more frequently. There are more significant variations in orbital elements expected in this case. At limit large N there are no any resonant phenomena. At small N we have chaotic motion due to close encounters. If N sufficiently large, effect of initial phase not significant. Conjunction. At fixed total mass of N-gon, m=m i N, interaction test particle with regular N- gon strongly depends on number of N. At conditions: ( r / R) / p / P n : m 1 outer ( r / R) / p / P n : m 1 inner where r, R and p, P central distances and periods of test particle and particle of ring respectively, resonance effects in motion test particle take place.
The longitude (solid line) and eccentricity (dashed line) on number of ring particles N dependence.
Fig.. Difference between perturbed and unperturbed positions for resonance case N=19 and for case N=100, close to continuous limit. Orbit with a=1.07
Generalization to Low-Elliptic case
Put ring of N particles m with small eccentricity in gravity field of central mass M. It is naturally to find stationary distribution of particles by true anomaly and about influence of mutual perturbations of particles on ring s eccentricity. There are two possibility. For the low-massive particles m<<m, their orbital distribution determined by keplerian angular velocity (L - angular momentum): d const dt L (1 0 ecos( )) mr Fig.. Orbital distribution of particles along eccentric ring (small mass particle case)
Fig. Declination of the shape of ring (from best fitted rotated ellipse) as function of true anomaly. Fig.. Eccentric ring shape change due to gravitational interactios of particles.
Other applications
U Thompson Heptagon i j j i log x x y y i j i j 1/ Model of N point-vortex on sphere by 6.Boatto S, Simo C.: Thompson heptagon: a case of bifurcation ad infinity., Physica D: Nonlinear Phenomena., Volume 7, Issue 14-17, p. 051-055. 008. Real stable vortex on Saturn pole
Results and conclusions In result, the dependence of libration points coordinates on mass and number of particles is studied. Coefficients A and B have a limit at large N, depends on N and m/m ratio. Accordingly, libration point coordinates close to a fixed value, slowly depends on mass ratio. On the other hand, at large m coordinates have a maximal value, which respect a case infinitesimal central mass. In case small m and small N we have a -body problem. Positions of non-collinear libration points are slightly differ from n-gon described circle. In conclusions, some generalizations of problem and possible way of future work are considered. Area, where potential allows a linearization exists in small vicinity of N-gon. However, oscillations even in this small area become strongly nonlinear due to differential rotation of system. In fact, approximation ring potential by N-gon is an example of a point-mass modeling. But specific of potential expansion gives an ability to study more complex (non-homogenous) systems on base of considered model.
Few References 1.Wintner, A. The Analytical Foundation of Celestial Mechanics, Princeton University Press (1941)..Rosaev A.E. The investigation of stationary points in central configuration dynamics. Proceedings of Libration Point Orbits and Applications Conference, Parador d'aiguablava, Girona, Spain 10-14 June, 00 Eds. by G. Gomez, M. W. Lo and.j.j. Masdemont, World Scientific Pub, New Jersey - London Singapore - Hon Kong, (00). p. 6-67..Rosaev A.E. The application of Computer Algebra to Central configuration dynamics. Abstr. of 6-th IMACS Int Conf on Applications of Computer Algebra (IMACS ACA 000), June 5-8, 000 St.Petersburg, Russia, p.51. 4.Szebehely, V. Theory of Orbits. The Restricted Problem of Three Bodies. Academic Press (1967). 5.Ollongren A. (1981). On a restricted (n+) body problem. Celestial Mechanics, 45,16-168. 6.Roberts G. Linear Stability of the 1+n-gon Relative Equilibrium Hamiltonian Systems and Celestial Mechanics (HAMSYS-98), World Scientific Monograph Series in Mathematics 6, 0-0, 000. 7.Boatto S, Simo C.: Thompson heptagon: a case of bifurcation ad infinity., Physica D: Nonlinear Phenomena., Volume 7, Issue 14-17, p. 051-055. 008.
Future works Test particle motion in N-gon potential Relative equilibrium in case non-equal masses Approximation of oblate planet potential by system of N-gons N-gon as a Toda chain with non-zero curvature
Thank You very much