Periodic Orbits in Rotating Second Degree and Order Gravity Fields

Periodic orbits in an arbitrary 2nd degree and order uniformly rotating gravity field are studied. We investigate the four equilibrium points in this gravity field. We see that close relation exists between the stability of these equilibria and the existence and stability of their nearby periodic orbits. We check the periodic orbits with non-zero periods. In our searching procedure for these periodic orbits, we remove the two unity eigenvalues from the state transition matrix to find a robust, non-singular linear map to solve for the periodic orbits. The algorithm converges well, especially for stable periodic orbits. Using the searching procedure, which is relatively automatic, we find five basic families of periodic orbits in the rotating second degree and order gravity field for planar motion, and discuss their existence and stability at different central body rotation rates.     

Periodic Orbits in Rotating Second Degree and Order Gravity Fields

Periodic orbits in an arbitrary 2nd degree and order uniformly rotating gravity field are studied. We investigate the four equilibrium points in this gravity field. We see that close relation exists between the stability of these equilibria and the existence and stability of their nearby periodic orbits. We check the periodic orbits with non-zero periods. In our searching procedure for these periodic orbits, we remove the two unity eigenvalues from the state transition matrix to find a robust, non-singular linear map to solve for the periodic orbits. The algorithm converges well, especially for stable periodic orbits. Using the searching procedure, which is relatively automatic, we find five basic families of periodic orbits in the rotating second degree and order gravity field for planar motion, and discuss their existence and stability at different central body rotation rates.     

Families of orbits in planar anisotropic fields

The aim of the planar inverse problem of dynamics is: given a monoparametric family of curves f(x, y) = c, find the potential V (x, y) under whose action a material point of unit mass can describe the curves of the family. In this study we look for V in the class of the anisotropic potentials V(x, y) = v(a²x² + y²), (a=constant). These potentials have been used lately in the search of connections between classical, quantum, and relativistic mechanics. We establish a general condition which must be satisfied by all the families produced by an anisotropic potential. We treat special cases regarding the families (e. g. families traced isoenergetically) and we present certain pertinent examples of compatible pairs of families of curves and anisotropic potentials. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

A basic system of equations in the field of a rotationally symmetric potential

We consider system of equations of motion in the field of a rotationally symmetric potential.     

Groups of symmetries in the two‐body problem associated to Einstein's PN field

The two‐body problem associated to a spherical post‐Newtonian (PN) field with Einsteinian parameterization is revisited from the single standpoint of symmetries. The corresponding vector fields, in Hamiltonian and standard polar coordinates, or in collision‐blow‐up and infinity‐blow‐upMcGehee‐type coordinates, present symmetries that form diffeomorphic commutative groups endowed with a Boolean structure. The existence of such symmetries is of much help in understanding characteristics of the global flow, or in finding symmetric periodic orbits in more complex problems depending on a small parameter.     

On the periodically evolving orbits in the singly averaged Hill problem

We continue to analyze the periodic solutions of the singly averaged Hill problem. We have numerically constructed the families of solutions that correspond to periodically evolving satellite orbits for arbitrary initial values of their eccentricities and inclinations to the plane of motion of the perturbing body. The solutions obtained are compared with the numerical solutions of the rigorous (nonaveraged) equations of the restricted circular three-body problem. In particular, we have constructed a periodically evolving orbit for which the well-known Lidov-Kozai mechanism manifests itself, just as in the doubly averaged problem.     

Eccentricity Generation in Hierarchical Triple Systems with Non-coplanar and Initially Circular Orbits

In a previous paper, we developed a technique for estimating the inner eccentricity in coplanar hierarchical triple systems on initially circular orbits, with comparable masses and with well-separated components, based on an expansion of the rate of change of the Runge-Lenz vector. Now, the same technique is extended to non-coplanar orbits. However, it can only be applied to systems with I ₀ < 39.23° or I ₀ > 140.77°, where I is the inclination of the two orbits, because of complications arising from the so-called ‘Kozai effect’. The theoretical model is tested against results from numerical integrations of the full equations of motion. This revised version was published online in July 2006 with corrections to the Cover Date.     

Gravothermal catastrophe in anisotropic spherical systems

In this paper we investigate the gravothermal instability of spherical stellar systems endowed with a radially anisotropic velocity distribution. We focus our attention on the effects of anisotropy on the conditions for the onset of instability and in particular we study the dependence of the spatial structure of critical models on the amount of anisotropy present in a system. The investigation has been carried out by the method of linear series which has already been used in the past to study the gravothermal instability of isotropic systems._   We consider models described by King, Wilson and Woolley–Dickens distribution functions. In the case of King and Woolley–Dickens models, our results show that, for quite a wide range of the amount of anisotropy in the system, the critical value of the concentration of the system (defined as the ratio of the tidal to the King core radius of the system) is approximately constant and equal to the corresponding value for isotropic systems. Only for very anisotropic systems does the critical value of the concentration start to change and it decreases significantly as the anisotropy increases and penetrates the inner parts of the system. For Wilson models the decrease of the concentration of critical models is preceded by an intermediate regime in which critical concentration increases, reaches a maximum and then starts to decrease. The critical value of the central potential always decreases as the anisotropy increases.     

Secular Effects on Orbits of Binary Stars Induced by Temporal Variation of Gravitational Constant (Case for Elliptical Orbit)

In this paper we investigate the influence of a varying gravitation constant on the orbits of celestial bodies. Regarding the eccentric anomaly as an independent variable, we find the solutions to the perturbed equations of motion. In the first order solutions, we find the secular and periodic variations in semi-major axis. For the other orbital elements only periodic variations exhibit. However in the second order solutions, the longitude of periastron and the mean longitude have secular terms. Applying the calculations to six selected binaries, we give the numerical estimations of the variations of orbits. These results are then carefully compared and discussed.     

Using wavelet approximation to study the influence of Jupiter on the orbital evolution of a distant satellite of Saturn

We suggest a nonstandard methodology for studying the influence of Jupiter on the secular orbital evolution of a distant satellite of Saturn. This influence is tangible only in short time spans near the times of the smallest separation between Jupiter and Saturn, i.e., when the heliocentric longitudes of the two planets coincide. These times are spaced about 20 years apart. To describe the jumplike behavior of perturbations, we suggest approximating the principal part of the perturbing function averaged over the satellite’s motion by a two-parameter exponential wavelet-type (burst) function. The subsequent averaging (smoothing) of the perturbing function allows us to eliminate the 20-year-period terms and obtain an approximate analytical solution in a special case of the problem. The results are illustrated by plots of the variations in the averaged perturbing function and the orbital eccentricity of Saturn’s outer satellite S/2000 S1, which is most strongly perturbed by Jupiter.     

Periodic solutions near a central libration point in the problem of the motion of a star inside an elliptical galaxy

We consider the three-dimensional problem of the motion of a star inside an inhomogeneous rotating elliptical galaxy with a homothetic density distribution. We construct and analyze the periodic solutions near a central libration point by using Lyapunov’s method.     

Particular solutions of the singly averaged Hill problem

We consider the case of averaging the perturbing function of the Hill problem over the fastest variable, the mean anomaly of the satellite. In integrable special cases, we found solutions to the evolutionary system of equations in elements.     

Gravitational radiation damping and the three-body problem

A model of three-body motion is developed which includes the effects of gravitational radiation reaction. The radiation reaction due to the emission of gravitational waves is the only post-Newtonian effect that is included here. For simplicity, all of the motion is taken to be planar. Two of the masses are viewed as a binary system, and the third mass, whose motion will be a fixed orbit around the centre-of-mass of the binary system, is viewed as a perturbation. This model aims to describe the motion of a relativistic binary pulsar that is perturbed by a third mass. Numerical integration of this simplified model reveals that, given the right initial conditions and parameters, one can see resonances. These ( m , n ) resonances are defined by the resonance condition,   mω =2 n Ω  , where m and n are relatively prime integers, and ω and Ω are the angular frequencies of the binary orbit and third mass orbit (around the centre-of-mass of the binary), respectively. The resonance condition consequently fixes a value for the semimajor axis of the binary orbit for the duration of the resonance; therefore the binary energy remains constant on average, while its angular momentum changes during the resonance.     

On the stability of particular solutions of the singly averaged Hill problem

We consider the particular solutions of the evolutionary system of equations in elements that correspond to planar and spatial circular orbits of the singly averaged Hill problem. We analyze the stability of planar and spatial circular orbits to inclination and eccentricity, respectively. We construct the instability regions of both particular solutions in the plane of parameters of the problem.     

On a two-body problem with periodically changing equivalent gravitational parameter

Assigning to the equivalent gravitational parameter of a two-body dynamic system, a periodic change of a small amplitude B and arbitrary frequency and phase, the behaviour of an elliptic-type orbit is studied. The first order (in B) perturbations of the orbital elements are determined by using Delaunay's canonical variables. According to the value of the ratio between oscillation frequency and dynamic frequency, three cases (non-resonant (NR), quasi-resonant (QR), and resonant (R) ones) are pointed out. The solution of motion equations shows that only in the QR and R cases there are elements (argument of pericentre and mean anomaly) affected by secular perturbations. The solutions are valid over prediction times of order of pericentre and mean anomaly) affected by secular perturbations. The solutions are valid over prediction times of order B⁻¹ in the NR case and B^−1/2 in the QR and R cases.     

天体轨道长期数值积分的误差估计方法

数值积分方法是进行天体力学研究的重要工具, 尤其对于行星历表的研究工作而言. 由于在使用数值方法计算天体轨道时, 最终误差通常是难以预知的, 所以在面对精度要求较高或者积分时间较长的工作时具体积分方案的设计---尤其是当使用定步长方法时的步长选择---需要十分谨慎, 因为这将意味着是否能在时间成本可以被接受的范围内使解的精度达到要求. 因此, 在使用数值方法解决实际问题时如何快速寻找效率与精度之间的最佳平衡点是每一个数值积分方法的设计者与使用者都会面临的难题. 为解决这一问题, 在定步长条件下对数值积分方法的舍入误差概率分布函数以及截断误差积累量对步长的依赖关系和随时间的增长关系进行了深入研究. 基于所得结论, 提出了一种仅需较少的数值实验资料即可对选择任意时间步长积分至任意积分时刻时的舍入误差概率分布函数与截断误差积累量进行准确估计的方法, 并使用Adams-Cowell方法对该误差估计方法在圆周期轨道条件下进行了验证. 该误差估计方法在未来有望用于不同数值算法的性能对比研究, 同时也可以对数值积分方法求解实际轨道问题时的决策工作带来重要帮助.     

引力常数随时间变化对双星轨道演变的长期效应(椭圆轨道情形)

给出了以偏近点角为自变量的变引力常数的摄动方程组的解.解包括轨道半长轴的长期和周期变化项,其他轨道根数在一阶解中无长期项,只有周期项.近星点经度和平经度在二阶解中显示长期项变化.给出了由于引力常数变化对双星轨道演变情况的数值估计,对结果做了讨论并给出结论.     

Electromagnetic Radiation and Motion of a Particle

We consider the motion of uncharged dust grains of arbitrary shape including the effects of electromagnetic radiation and thermal emission. The resulting relativistically covariant equation of motion is expressed in terms of standard optical parameters. Explicit expressions for secular changes of osculating orbital elements are derived in detail for the special case of the Poynting-Robertson effect. Two subcases are considered: (i) central acceleration due to gravity and the radial component of radiation pressure independent of the particle velocity, (ii) central acceleration given by gravity and the radiation force as the disturbing force. The latter case yields results which may be compared with secular orbital evolution in terms of orbital elements for an arbitrarily shaped dust particle. The effects of solar wind are also presented. This revised version was published online in July 2006 with corrections to the Cover Date.     

Elliptic-type motion in Fock's gravitational field

The elliptic-type motion in the gravitational field found by Fock as exact solution of Einstein's vacuum equations in the case of spherical symmetry (Solution called here Fock's gravitational field) is studied by means of a classic method based on the perturbation theory. Regarding the deviations of the orbit from a Kepleian orbit as perturbations, the first and second order variations of the Keplerian orbital elements over one nodal period as well as those of the nodal period itself are determined.