The Orbital Evolution of Two Sounding Satellites and Analysis of the Accuracy of Orbit Determination

The satellites TC-1 and TC-2 are the two Chinese satellites with great elliptical orbits which are still in orbit around the earth at present. Since the launch the orbits of the two satellites have continuously evolved, which has a certain effect on the orbit determination and prediction precision. The regularities of the orbital evolution of the two sounding satellites are qualitatively and quantitatively analyzed. Under the current tracking mode the corresponding prediction precision of orbit determination is analyzed based on the different stages of the orbital evolution, thereby providing the basis for the adjustment of planning mode by the satellite application departments and the guarantee of normal satellite payload. Finally, the orbital lifetimes of the two satellites are predicted through the trend of the orbital evolution.     

Mass of Mercury from observations of asteroids

The Sun-to-Mercury mass ratio adopted by the International Astronomical Union (6023600 ± 250) was obtained in 1987 by analyzing Mariner 10 observations (Anderson et al. 1987) and since then has not been improved. The large number of asteroids in Mercury-approaching orbits and the ever-increasing accuracy of their observations allow the mass of Mercury to be estimated by a different method. We have improved the orbital parameters of 43 asteroids and obtained 6023440 ± 530 for the Sun-to-Mercury mass ratio through a simultaneous solution based on their optical and radar observations. A further improvement in this estimate is possible in the immediate future owing to the rapid increase in the number of known asteroids whose observations can be used to solve this problem.     

Two-integral distribution functions for axisymmetric stellar systems with separable densities

We show different expressions of distribution functions (DFs) which depend only on the two classical integrals of the energy and the magnitude of the angular momentum with respect to the axis of symmetry for stellar systems with known axisymmetric densities. The density of the system is required to be a product of functions separable in the potential and the radial coordinates, where the functions of the radial coordinate are powers of a sum of a square of the radial coordinate and its unit scale. The even part of the two-integral DF corresponding to this type of density is in turn a sum or an infinite series of products of functions of the energy and of the magnitude of the angular momentum about the axis of symmetry. A similar expression of its odd part can be also obtained under the assumption of the rotation laws. It can be further shown that these expressions are in fact equivalent to those obtained by using Hunter & Qian's contour integral formulae for the system. This method is generally computationally preferable to the contour integral method. Two examples are given to obtain the even and odd parts of their two-integral DFs. One is for the prolate Jaffe model and the other for the prolate Plummer model.
It can be also found that the Hunter–Qian contour integral formulae of the two-integral even DF for axisymmetric systems can be recovered by use of the Laplace–Mellin integral transformation originally developed by Dejonghe.     

Three-dimensional numerical simulation of a nonstationary gravitating <Emphasis Type="Italic">N</Emphasis>-body system with gas

We developed a three-dimensional numerical model to investigate nonstationary processes in gravitating N-body systems with gas. We used efficient algorithms for solving the Vlasov and Poisson equations that included the evolutionary processes under consideration, which ensures rapid convergence at high accuracy. We give examples of the numerical solution of the problem on the growth of physical instability in the model of a flat rotating disk with a gaseous component and its three-dimensional dynamics under various initial conditions including a nonzero velocity dispersion along the rotation axis.     

Galactic Sun's motion in the cold dark matter,MOdified Newtonian Dynamics and modified gravity scenarios

We numerically integrate the equations of motion of the Sun in Galactocentric Cartesian rectangular coordinates for –4.5 Gyr ≤ t ≤ 0 in Newtonian mechanics with two different models for the Cold Dark Matter (CDM) halo, in MOdified Newtonian Dynamics (MOND) and in MOdified Gravity (MOG) without resorting to CDM. The initial conditions used come from the latest kinematical determination of the 3D Sun's motion in the Milky Way (MW) by assuming for the rotation speed of the Local Standard of Rest (LSR) the recent value Θ₀ = 268 km s^–1 and the IAU recommended value Θ₀ = 220 km s^–1; the Sun is assumed located at 8.5 kpc from the Galactic Center (GC). For Θ₀ = 268 km s^–1 the birth of the Sun, 4.5 Gyr ago, would have occurred at large Galactocentric distances (12–27 kpc depending on the model used), while for Θ₀ = 220 km s^–1 it would have occurred at about 8.8–9.3 kpc for almost all the models used. The integrated trajectories are far from being circular, especially for Θ₀ = 268 km s^–1, and differ each other with the CDM models yielding the widest spatial extensions for the Sun's orbital path (© 2009 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

On the habitability of the OGLE-2006-BLG-109L planetary system

We investigate the dynamics of putative Earth-mass planets in the habitable zone (HZ) of the extrasolar planetary system OGLE-2006-BLG-109L, a close analogue of the Solar system. Our work is inspired by the work of Malhotra & Minton. Using the linear Laplace–Lagrange theory, they identified a strong secular resonance that may excite large eccentricity of orbits in the HZ. However, due to uncertain or unconstrained orbital parameters, the subsystem of Jupiters may be found in a dynamically active region of the phase space spanned by low-order mean-motion resonances. To generalize this secular model, we construct a semi-analytical averaging method in terms of the restricted problem. The secular orbits of large planets are approximated by numerically averaged osculating elements. They are used to calculate the mean orbits of terrestrial planets by means of a high-order analytic secular theory developed in our previous works. We found regions in the parameter space of the problem in which stable, quasi-circular orbits in the HZ are permitted. The excitation of eccentricity in the HZ strongly depends on the apsidal angle of jovian orbits. For some combinations of that angle, eccentricities and semimajor axes consistent with the observations, a terrestrial planet may survive in low eccentric orbits. We also study the effect of post-Newtonian gravity correction on the innermost secular resonance.     

The Data Depth-weight-kernel Estimation of the Model Error of Satellite Orbit Determination

It is an objective fact that there exists error in the satellite dynamic model and it will be transferred to satellite orbit determination algorithm, forming a part of the connotative model error. Mixed with the systematic error and random error of the measurements, they form the unitive model error and badly restrict the precision of the orbit determination. We deduce in detail the equations of orbit improvement for a system with dynamic model error, construct the parametric model for the explicit part of the model and nonparametric model for the error that can not be explicitly described. We also construct the partially linear orbit determination model, estimate and fit the model error using a two-stage estimation and a kernel function estimation, and finally make the corresponding compensation in the orbit determination. Beginning from the data depth theory, a data depth weight kernel estimator for model error is proposed for the sake of promoting the steadiness of model error estimation. Simulation experiments of SBSS are performed. The results show clearly that the model error is one of the most important effects that will influence the precision of the orbit determination. The kernel function method can effectively estimate the model error, with the window width as a major restrict parameter. A data depth-weight-kernel estimation, however, can improve largely the robustness of the kernel function and therefore improve the precision of orbit determination.