Virial oscillations of celestial bodies

A general approach to the solution of the perturbed oscillation problem for celestial bodies is considered. The solution sought describes unperturbed virial oscillations (zero approximation) affected by external perturbing effects. In the general case, these perturbations can be expressed by an arbitrary given function of time, Jacobi's function and its first derivative. Standard methods and modes of perturbation theory are used for solution of the problem.It is shown that while studying the evolution of a celestial body as a dissipative system in the framework of perturbed virial oscillations, the analytical expression for perturbing function can be derived, assuming the celestial body to be an oscillating electrical dipole emitting electromagnetic energy.The general covariant form of Jacobi's equation is derived and its spur is examined. It is shown that the scalar form of Jacobi's equation appears to be more universal than Newton's laws of motion from which it is derived.     

Determination of Large Asteroid Masses by the Dynamical Method

Masses of 19 asteroids have been determined from the analysis of their gravitational effect on the motion of perturbed bodies. The following asteroids were selected as perturbed bodies: (1) those which had single close encounters with the perturbing asteroid; (2) those whose mean motion was in a 1 : 1 commensurability with that of the perturber and which had close or moderate recurrent encounters with the perturber. The perturber mass was determined from observations of several tens of perturbed asteroids that were selected from these two groups. The selection criterion was the error of the mass determined from observations of only one asteroid. Positional observations of the asteroids on the interval 1900–2002 were used. The masses were determined with errors by an order-half an order of magnitude smaller than the masses found. The results are compared with those of other authors.     

Determining asteroid masses from perturbations on Mars

The orbit of Mars is perturbed more than 5 m, a value compatible with the accuracy of the Viking lander ranging data, by about three dozen asteroids. In addition to larger asteroids throughout the belt, significant perturbations of long period are generated by smaller objects near commensurabilities with Mars. The largest periodic terms induced by 1 Ceres and 2 Pallas have amplitudes of 0.8 and 0.2 km, respectively, both with 10-year periods. Due to a near commensurability, 4 Vesta causes a 5-km, 52-year term. While the Viking ranges will yield significant mass determinations for the largest three asteroids, and some of the smaller bodies should be detectable, it will be difficult to seperate the smaller bodies with useful accuracies. Accurate discrimination must await range data from future missions to Mars or other bodies in the neighborhood of the asteroid belt. The Viking ranges can also yield improved masses for the outer planets (except Pluto), an application which is being exploited by groups analyzing these data. Uncertainties in the asteroid masses limit the ultimate accuracy of the Viking determinations of both the long time scale motion of the system the inner four planets with respect to an inertial frame and the rate of change of the gravitational constant.     

Astrometric studies of the results of a new reduction of old photographic observations of the Saturnian System based on the comparison with the modern theories of satellite motion

The paper shows the possibility of increasing the accuracy of the results of photographic observations of Saturn and its moons made in the 1970s and reduced using the old reference star catalogues and semiautomatic measurements. New celestial coordinates of the moons (from the third to the eighth), “satellite minus satellite” relative moon coordinates, and Saturn coordinates by positions of satellites are obtained without measuring its images. The results are stored in the Pulkovo Observatory database on the Solar System bodies and are available online at www.puldb.ru. The efficiency of the reduction method based on digitizing of astronegatives using 21 Mpx Canon digital camera and IZMCCD software is shown. The comparison of new results of old observations with the latest theories of moon motion has revealed a significant increase in satellite positioning accuracy. The investigation of the differences (O–C) of celestial coordinates from satellite positions in their apparent Saturn-centric orbits has revealed a noticeable motion of the differences (O–C) in right ascension depending on their distances from Saturn for all moons.     

第三体摄动分析解的一种表达式

在太阳系中,大行星、小行星和卫星（包括自然卫星和人造卫星）等对应的运动问题,都可以处理成受摄二体问题,而摄动源又多为第三体,作为第三体的摄动天体,有的比运动天体离中心天体近,有的则相反,前者称为内摄内体,全者则称为外摄天体,对一个具体的运动天体,可以同时出现这两个摄动天体,但是,只要运动天体与摄动天体的轨道都建立在以中心天体（质心）为坐标原点的同一坐标系内,那么在一定条件下（即除运动天体与摄动天体     

Dynamical parameters of spacecraft motion near small celestial body

The paper considers how a spacecraft can be put into orbit around a small asteroid to function as its artificial moon. We study the general behavior of perturbations that affect the current coordinates of an orbiting spacecraft and estimate the perturbations caused by the main perturbing factors, i.e., (1) the irregular shape of an asteroid and (2) celestial bodies of the Solar System. With specific orbital parameters, a long-term targeted operation of a spacecraft can be actualized in a mission to the asteroid Apophis where the spacecraft will carry a radio beacon transponder.     

Computer simulation of the formation of asteroid belt structure

The system of two gravitational centers with variable separation between components one of which (the primary) loses its mass onto another (the secondary) is investigated under condition of total mass and angular momentum conservation. When the primary/secondary mass ratio becomes about that of Jupiter/Sun the small bodies ejected with the gaseous matter through the inner Lagrange point from the Roche lobe of the primary form a ring similar to the asteroid belt of the solar system. The formation of ring structure is calculated by numerical integration of Newtonian equations of N-body problem in orbital plane of the gravitational centers. The results are compared with the planar subsystem of the asteroid belt. The presence of the main gaps in the distribution of their mean motions at 2/1, 3/1, 5/2 and some other commensurabilities with the primary mean motion is found. More fine details of the belt structure are obtained, e.g. the gap asymmetry and a qualitative agreement with the eccentricity distribution. Within the scope of the same model the external part of the ring is investigated all the pairwise interactions being included. The clustering of bodies near 3/2 commensurability isolated from the main belt by the wide gap centered at 5/3 commensurability is obtained. It is supposed that the ring structure and the interplanetary spacing law for the terrestrial planets are due to the same mechanism.     

Planets in the asteroid belt

Abstract— The main asteroid belt has lost >99.9% of its solid mass since the time at which the planets were forming, according to models for the protoplanetary nebula. Here we show that the primordial asteroid belt could have been cleared efficiently if much of the original mass accreted to form planetsized bodies, which were capable of perturbing one another into unstable orbits. We provide results from 25 N‐body integrations of up to 200 planets in the asteroid belt, with individual masses in the range 0.017–0.33 Earth masses. In the simulations, these bodies undergo repeated close encounters which scatter one another into unstable resonances with the giant planets, leading to collision with the Sun or ejection from the solar system. In response, the giant planets' orbits migrate radially and become more circular. This reduces the size of the main‐belt resonances and the clearing rate, although clearing continues. If ～3 Earth masses of material was removed from the belt this way, Jupiter and Saturn would initially have had orbital eccentricities almost twice their current values. Such orbits would have made Jupiter and Saturn 10–100x more effective at clearing material from the belt than they are on their current orbits. The time required to remove 90% of the initial mass from the belt depends sensitively on the giant planets' orbits, and weakly on the masses of the asteroidal planets. 18 of the 25 simulations end with no planets left in the belt, and the clearing takes up to several hundred million years. Typically, the last one or two asteroidal planets are removed by interactions with planets in the terrestrial region     

A Method for the Determination of an Intermediate Orbit from Three Positions of the Small Body on the Celestial Sphere

A new method is suggested for finding the preliminary orbit from three complete measurements of the angular coordinates of a celestial body developed by analogy with the classic Lagrange–Gauss method. The proposed method uses the intermediate orbit that we had constructed in an earlier paper based on two position vectors and the corresponding time interval. This intermediate orbit allows for most of the perturbations in the motion of the body. Using the orbital motion of asteroid 1566 Icarus as an example, we compare the results obtained by applying the classic and the new method. The comparison shows the new method to be highly efficient for studying perturbed motion. It is especially efficient if applied to high-precision observational data covering short orbital arcs.     

Approximate solutions of non-linear circular orbit relative motion in curvilinear coordinates

A compact, time-explicit, approximate solution of the highly non-linear relative motion in curvilinear coordinates is provided under the assumption of circular orbit for the chief spacecraft. The rather compact, three-dimensional solution is obtained by algebraic manipulation of the individual Keplerian motions in curvilinear, rather than Cartesian coordinates, and provides analytical expressions for the secular, constant and periodic terms of each coordinate as a function of the initial relative motion conditions or relative orbital elements. Numerical test cases are conducted to show that the approximate solution can be effectively employed to extend the classical linear Clohessy–Wiltshire solution to include non-linear relative motion without significant loss of accuracy up to a limit of 0.4–0.45 in eccentricity and 40–45\(^\circ \) in relative inclination for the follower. A very simple, quadratic extension of the classical Clohessy–Wiltshire solution in curvilinear coordinates is also presented.     

Equations of motion of the restricted problem of three bodies with variable mass

In this paper we have deduced the differential equations of motion of the restricted problem of three bodies with decreasing mass, under the assumption that the mass of the satellite varies with respect to time. We have applied Jeans law and the space time transformation contrast to the transformation of Meshcherskii. The space time transformation is applicable only in the special casen=1,k=0,q=1/2. The equations of motion of our problem differ from the equations of motion of the restricted three body problem with constant mass only by small perturbing forces.     

On the motion of a mass point inside a rotating inhomogeneous ellipsoidal body

The plane motion of a mass point inside an inhomogeneous rotating ellipsoidal body with a homothetic density distribution is considered. The force function of the problem is expanded in terms of the ellipsoid's second eccentricities up to the fourth order, which are taken as small parameters. We derive an expression for the perturbing function and solve the equations of perturbed motion in orbital elements.     

Equations of motion of elliptic restricted problem of three bodies with variable mass

The differential equations of motion of the elliptic restricted problem of three bodies with decreasing mass are derived. The mass of the infinitesimal body varies with time. We have applied Jeans' law and the space-time transformation of Meshcherskii. In this problem the space-time transformation is applicable only in the special case whenn=1,k=0,q=1/2. We have applied Nechvile's transformation for the elliptic problem. We find that the equations of motion of our problem differ from that of constant mass only by a small perturbing force.     

The orbit of asteroid (99942) Apophis as determined from optical and radar observations

The results of improving the orbit accuracy for the asteroid Apophis and the circumstances of its approach to Earth in 2029 are described. Gravitational perturbations from all of the major planets and Pluto, Ceres, Pallas, and Vesta are taken into account in the equations of motion of the asteroid. Relativistic perturbations from the Sun and perturbations due to the oblateness of the Sun and Earth and due to the light pressure are also included in the model. Perturbations from the Earth and Moon are considered separately. The coordinates of the perturbing bodies are calculated using DE405. The phase correction and the gravitational deflection of light are taken into account. The numerical integration of the equations of motion and equations in variations is performed by the 15th-order Everhart method. The error of the numerical integration over the 2005–2029 interval, estimated using forward and backward computations, is not more than 3 × 10^?11 AU. Improved coordinates and velocities at epoch JD2454200.5 (April 10, 2007) were obtained applying the weighted leastsquares fit. For the period from March 15, 2004, to August 16, 2006, 989 optical and 7 radar observations were used. The resulting system represents the optical observations with an error of 0.37 (66 conditional equations were rejected). The residuals of the radar observations are an order, or more, smaller than their errors. The system of Apophis’ elements and the estimates of their precision obtained in this study are in perfect agreement with the results published by other authors. The minimum Apophis-Earth distance is about 38 200 km on April 13, 2029. This estimate agrees to within 20 km with those calculated based on other published systems of elements. The effect of some model components on the minimum distance is estimated.     

Accurate computation of highly eccentric satellite orbits

For computing highly eccentric (e0.9) Earth satellite orbits with special perturbation methods, a comparison is made between different schemes, namely the direct integration of the equations of motion in Cartesian coordinates, changes of the independent variable, use of a time element, stabilization and use of regular elements. A one-step and a multi-step integration are also compared.It is shown that stabilization and regularization procedures are very helpful for non or smoothly perturbed orbits. In practical cases for space research where all perturbations are considered, these procedures are no longer so efficient. The recommended method in these cases is a multi-step integration of the Cartesian coordinates with a change of the independent variable defining an analytical step size regulation. However, the use of a time element and a stabilization procedure for the equations of motion improves the accuracy, except when a small step size is chosen.     

Time to Perform the Physical Observations of Asteroids

The astrometric observations of asteroids are often performednear their opposition, where they have faster apparent motions andbrighter visual magnitudes. However, the physical observations ofasteroids (photometric and spectral) often require longer exposuretimes as well as brighter magnitudes for better signal-to-noiseratio, which are combined effects of both apparent moving speedsand visual magnitudes. We derive the equations of the apparentmotion of asteroid in geocentric ecliptic coordinates. Comparisonof the apparent magnitudes of asteroid at opposition and atstationary shows that the magnitude differences in the two casesare around 1 magnitude in average for most main belt asteroids,but are much larger for Near Earth Asteroids. Combining withcomparison of asteroid apparent motion, the proper time forasteroid physical observation in different cases are discussed.     

Mutual gravitational potential and torque of solid bodies via inertia integrals

The mutual gravitational potential and the mutual gravitational torque of two bodies of arbitrary shape are expanded to the fourth order. The derivations are based on Cartesian coordinates, inertia integrals with relation to the principal reference frames of each body, and the relative rotation matrix. The current formulation is convenient to utilize in high precision problems in rotational dynamics.     

Numerical stabilization of the differential equations of Keplerian motion

A stabilization of the classical equations of two-body motion is offered. It is characterized by the use of the regularizing independent variable (eccentric anomaly) and by the addition of a control-term to the differential equations. This method is related to the KS-theory (Stiefel, 1970) which performed for the first time a stabilization of the Kepler motion. But in contrast to the KS-theory our method does not transform the coordinates of the particle. As far as the theory of stability and the numerical experiments are concerned we restrict ourselves to thepure Kepler motion. But, of course, the stabilizing devices will also improve the accuracy of the computation of perturbed orbits. We list, therefore, also the equations of the perturbed motion.     

Computing secular motion under slowly rotating quadratic perturbation

We consider secular perturbations of nearly Keplerian two-body motion under a perturbing potential that can be approximated to sufficient accuracy by expanding it to second order in the coordinates. After averaging over time to obtain the secular Hamiltonian, we use angular momentum and eccentricity vectors as elements. The method of variation of constants then leads to a set of equations of motion that are simple and regular, thus allowing efficient numerical integration. Some possible applications are briefly described.     

The general motion of three mutually-interacting,axisymmetric spherical rigid bodies

In this paper, the translational-rotational motions of an axisymmetric rigid body and two spherical rigid bodies under the influence of their mutual gravitational attraction are considered. The equations of motion in the canonical elements of Delaunay-Andoyer are obtained. The elements of motion in the zero and first approximations can be determined.