Some orbital characteristics of lunar artificial satellites

In this paper we present an analytical theory with numerical simulations to study the orbital motion of lunar artificial satellites. We consider the problem of an artificial satellite perturbed by the non-uniform distribution of mass of the Moon and by a third-body in elliptical orbit (Earth is considered). Legendre polynomials are expanded in powers of the eccentricity up to the degree four and are used for the disturbing potential due to the third-body. We show a new approximated equation to compute the critical semi-major axis for the orbit of the satellite. Lie-Hori perturbation method up to the second-order is applied to eliminate the terms of short-period of the disturbing potential. Coupling terms are analyzed. Emphasis is given to the case of frozen orbits and critical inclination. Numerical simulations for hypothetical lunar artificial satellites are performed, considering that the perturbations are acting together or one at a time.     

Luni-solar perturbations of an Earth satellite

Luni-solar perturbations of the orbit of an artificial Earth satellite are given by modifying the analytical theory of an artificial lunar satellite derived by the author in recent papers. Expressions for the first-order changes, both secular and periodic, in the elements of the geocentric Keplerian orbit of the earth satellite are given, the moon's geocentric orbit, including solar perturbations in it, being found by using Brown's lunar theory.The effects of Sun and Moon on the satellite orbit are described to a high order of accuracy so that the theory may be used for distant earth satellites.     

Gravitational perturbations of equatorial orbits

Asymptotic solutions are developed for the motion of a geocentric satellite in the equatorial plane due to gravitational perturbations such as nonsphericity (especially oblateness) of the primary body. Axisymmetric potentials are considered. A class of transformations is developed and the equations of motion are solved by the method of generalized multiple scales. Further it is shown that the equations of motion can be transformed into the required form to within any specified degree of accuracy. The transformations form an Abelian group of infinite order which leaves the differential equations of motion invariant. Solutions are developed in terms of elementary functions instead of elliptic or other higher transcendental functions and are shown to agree with known results.This investigation was carried out under NASA Grant NGR-31-001-152 with the author as a consultant to Princeton University.     

Ballooning modes in the inner magnetosphere of the Earth

In this paper the low-frequency ideal MHD (magnetohydrodynamical) perturbations in the inner magnetosphere of the Earth are studied. The set of partial differential equations obtained from the MHD equations in the ballooning approximation and the dipole model of the geomagnetic field is used for this purpose. These equations describe both small-scale and large-scale perturbations in the magnetospheric plasmas. In the “cold” plasma approximation the obtained equations describe poloidal and toroidal standing Alfvén modes. The account of plasma pressure leads to the appearance of an additional type of oscillations—the slow magnetosonic modes. The stability of the magnetospheric plasma with respect to the ballooning perturbations was analyzed. We describe the ballooning perturbations taking into account a coupling between the poloidal Alfvén modes and the slow magnetosonic modes.     

Nonlinear theory of secular perturbations of satellites of an oblate planet

Anonlinear analytical theory of secular perturbations in the problem of the motion of a systemof small bodies around a major attractive center has been developed. Themutual perturbations of the satellites and the influence of the oblateness of the central body are taken into account in the model. In contrast to the classical Laplace-Lagrange theory based on linear equations for Lagrange elements, the third-degree terms in orbital eccentricities and inclinations are taken into account in the equations. The corresponding improvement of the solution turns out to be essential in studying the evolution of orbits over long time intervals. A program inC has been written to calculate the corrections to the fundamental frequencies of the solution and the third-degree secular perturbations in orbital eccentricities and inclinations. The proposed method has been applied to investigate the motion of the major Uranian satellites. Over time intervals longer than 100 years, allowance for the nonlinear terms in the equations is shown to give corrections to the coordinates of Miranda on the order of the orbital eccentricity, which is several thousand kilometers in linear measure. For other satellites, the effect of allowance for the nonlinear terms turns out to be smaller. Obviously, when a general analytical theory of motion for the major Uranian satellites is constructed, the nonlinear terms in the equations for the secular perturbations should be taken into account.     

Approximate analytic method for high-apogee twelve-hour orbits of artificial Earth’s satellites

We propose an approach to the study of the evolution of high-apogee twelve-hour orbits of artificial Earth’s satellites. We describe parameters of the motion model used for the artificial Earth’s satellite such that the principal gravitational perturbations of the Moon and Sun, nonsphericity of the Earth, and perturbations from the light pressure force are approximately taken into account. To solve the system of averaged equations describing the evolution of the orbit parameters of an artificial satellite, we use both numeric and analytic methods. To select initial parameters of the twelve-hour orbit, we assume that the path of the satellite along the surface of the Earth is stable. Results obtained by the analytic method and by the numerical integration of the evolving system are compared. For intervals of several years, we obtain estimates of oscillation periods and amplitudes for orbital elements. To verify the results and estimate the precision of the method, we use the numerical integration of rigorous (not averaged) equations of motion of the artificial satellite: they take into account forces acting on the satellite substantially more completely and precisely. The described method can be applied not only to the investigation of orbit evolutions of artificial satellites of the Earth; it can be applied to the investigation of the orbit evolution for other planets of the Solar system provided that the corresponding research problem will arise in the future and the considered special class of resonance orbits of satellites will be used for that purpose.     

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.     

Relativistic effects for near-earth satellite orbit determination

The relativistic formulations for the equations which describe the motion of a near-Earth satellite are compared for two commonly used coordinate reference systems (RS). The discussion describes the transformation between the solar system barycentric RS and both the non-inertial and inertial geocentric RSs. A relativistic correction for the Earth's geopotential expressed in the solar system barycentric RS and the effect of geodesic precession on the satellite orbit in the geocentric RS are derived in detail. The effect of the definition of coordinate time on scale is also examined. A long-arc solution using 3 years of laser range measurements of the motion of the Lageos satellite is used to demonstrate that the effects of relativity formulated in the geocentric RS and in the solar system barycentric RS are equivalent to a high degree of accuracy.     

Numerical integration of dynamical systems with Lie series

The integration of the equations of motion in gravitational dynamical systems—either in our Solar System or for extra-solar planetary systems—being non integrable in the global case, is usually performed by means of numerical integration. Among the different numerical techniques available for solving ordinary differential equations, the numerical integration using Lie series has shown some advantages. In its original form (Hanslmeier and Dvorak, Astron Astrophys 132, 203 1984), it was limited to the N-body problem where only gravitational interactions are taken into account. We present in this paper a generalisation of the method by deriving an expression of the Lie terms when other major forces are considered. As a matter of fact, previous studies have been done but only for objects moving under gravitational attraction. If other perturbations are added, the Lie integrator has to be re-built. In the present work we consider two cases involving position and position-velocity dependent perturbations: relativistic acceleration in the framework of General Relativity and a simplified force for the Yarkovsky effect. A general iteration procedure is applied to derive the Lie series to any order and precision. We then give an application to the integration of the equation of motions for typical Near-Earth objects and planet Mercury.     

The theory of asynchronous relative motion I: time transformations and nonlinear corrections

Using alternative independent variables in lieu of time has important advantages when propagating the partial derivatives of the trajectory. This paper focuses on spacecraft relative motion, but the concepts presented here can be extended to any problem involving the variational equations of orbital motion. A usual approach for modeling the relative dynamics is to evaluate how the reference orbit changes when modifying the initial conditions slightly. But when the time is a mere dependent variable, changes in the initial conditions will result in changes in time as well: a time delay between the reference and the neighbor solution will appear. The theory of asynchronous relative motion shows how the time delay can be corrected to recover the physical sense of the solution and, more importantly, how this correction can be used to improve significantly the accuracy of the linear solutions to relative motion found in the literature. As an example, an improved version of the Clohessy-Wiltshire (CW) solution is presented explicitly. The correcting terms are extremely compact, and the solution proves more accurate than the second and even third order CW equations for long propagations. The application to the elliptic case is also discussed. The theory is not restricted to Keplerian orbits, as it holds under any perturbation. To prove this statement, two examples of realistic trajectories are presented: a pair of spacecraft orbiting the Earth and perturbed by a realistic force model; and two probes describing a quasi-periodic orbit in the Jupiter-Europa system subject to third-body perturbations. The numerical examples show that the new theory yields reductions in the propagation error of several orders of magnitude, both in position and velocity, when compared to the linear approach.     

Kinematical modeling of the Earth rotation,focusing on the Oppolzer terms in a rigid Earth and the Oppolzer-like terms in an elastic Earth

Under perturbations from outer bodies, the Earth experiences changes of its angular momentum axis, figure axis and rotational axis. In the theory of the rigid Earth, in addition to the precession and nutation of the angular momentum axis given by the Poisson terms, both the figure axis and the rotational axis suffer forced deviation from the angular momentum axis. This deviation is expressed by the so-called Oppolzer terms describing separation of the averaged figure axis, called CIP (Celestial Intermediate Pole) or CEP (Celestial Ephemeris Pole), and the mathematically defined rotational axis, from the angular momentum axis. The CIP is the rotational axis in a frame subject to both precession and nutation, while the mathematical rotational axis is that in the inertial (non-rotating) frame. We investigate, kinematically, the origin of the separation between these two axes—both for the rigid Earth and an elastic Earth. In the case of an elastic Earth perturbed by the same outer bodies, there appear further deviations of the figure and rotational axes from the angular momentum axis. These deviations, though similar to the Oppolzer terms in the rigid Earth, are produced by quite a different physical mechanism. Analysing this mechanism, we derive an expression for the Oppolzer-like terms in an elastic Earth. From this expression we demonstrate that, under a certain approximation (in neglect of the motion of the perturbing outer bodies), the sum of the direct and convective perturbations of the spin axis coincides with the direct perturbation of the figure axis. This equality, which is approximate, gets violated when the motion of the outer bodies is taken into account.     

General relativity and satellite orbits: The motion of a test particle in the Schwarzschild metric

The motion of a satellite with negligible mass in the Schwarzschild metric is treated as a problem in Newtonian physics. The relativistic equations of motion are formally identical with those of the Newtonian case of a particle moving in the ordinary inverse-square law field acted upon by a disturbing function which varies asr ^?3. Accordingly, the relativistic motion is treated with the methods of celestial mechanics. The disturbing function is expressed in terms of the Keplerian elements of the orbit and substituted into Lagrange's planetary equations. Integration of the equations shows that a typical Earth satellite with small orbital eccentricity is displaced by about 17 cm from its unperturbed position after a single orbit, while the periodic displacement over the orbit reaches a maximum of about 3 cm. Application of the equations to the planet Mercury gives the advance of the perihelion and a total displacement of about 85 km after one orbit, with a maximum periodic displacement of about 13 km.     

Why Newtonian gravity is reliable in large-scale cosmological simulations

Until now, it has been common to use Newtonian gravity to study the non-linear clustering properties of large-scale structures. Without confirmation from Einstein's theory, however, it has been unclear whether we can rely on the analysis (e.g. near the horizon scale). In this work we will provide confirmation of the use of Newtonian gravity in cosmology, based on the relativistic analysis of weakly non-linear situations to third order in perturbations. We will show that, except for the gravitational-wave contribution, the relativistic zero-pressure fluid equations perturbed to second order in a flat Friedmann background coincide exactly with the Newtonian results. We will also present the pure relativistic correction terms appearing in the third order. The third-order correction terms show that these terms are the linear-order curvature perturbation times the second-order relativistic/Newtonian terms. Thus, the pure general relativistic corrections in the third order are independent of the horizon scale and are small when considering the large-scale structure of the Universe because of the low-level temperature anisotropy of the cosmic microwave background radiation. Since we include the cosmological constant, our results are relevant to currently favoured cosmology. As we prove that the Newtonian hydrodynamic equations are valid in all cosmological scales to second order, and that the third-order correction terms are small, our result has the important practical implication that one can now use the large-scale Newtonian numerical simulation more reliably as the simulation scale approaches and even goes beyond the horizon. In a complementary situation, where the system is weakly relativistic (i.e. far inside the horizon) but fully non-linear, we can employ the post-Newtonian approximation. We also show that in large-scale structures, the post-Newtonian effects are quite small.     

The Artificial Earth Satellite Motion in an Oblate,Rotating Atmosphere with Symmetrical Diurnal Effect

The equations describing the disturbed motion of an artificial Earth satellite in the atmosphere are integrated by using a new variable instead of the true anomaly. The atmospheric flattening and rotation, the linear variation of the density scale height with the altitude, the symmetrical diurnal effect and the variation of its amplitude with the heigth are taken into account. Approximate formulae for the perturbations in all the orbital elements over a revolution of the satellite are given.     

Lunar frozen orbits revisited

Lunar frozen orbits, characterized by constant orbital elements on average, have been previously found using various dynamical models, incorporating the gravitational field of the Moon and the third-body perturbation exerted by the Earth. The resulting mean orbital elements must be converted to osculating elements to initialize the orbiter position and velocity in the lunar frame. Thus far, however, there has not been an explicit transformation from mean to osculating elements, which includes the zonal harmonic \(J_2\), the sectorial harmonic \(C_{22}\), and the Earth third-body effect. In the current paper, we derive the dynamics of a lunar orbiter under the mentioned perturbations, which are shown to be dominant for the evolution of circumlunar orbits, and use von Zeipel’s method to obtain a transformation between mean and osculating elements. Whereas the dynamics of the mean elements do not include \(C_{22}\), and hence does not affect the equilibria leading to frozen orbits, \(C_{22}\) is present in the mean-to-osculating transformation, hence affecting the initialization of the physical circumlunar orbit. Simulations show that by using the newly-derived transformation, frozen orbits exhibit better behavior in terms of long-term stability about the mean values of eccentricity and argument of periapsis, especially for high orbits.     

Special Solutions in a Simplified Restricted Three-Body Problem with Gravitational Radiation Taken into Account

The distinctive feature of the relativistic restricted three-body problem within the c ^–5 order of accuracy (2 post-Newtonian approximation) is the presence of the gravitational radiation. To simplify the problem the motion of the massive binary components is assumed to be quasi-circular. In terms of time these orbits have linearly changing radii and quadratically changing phase angles. By substituting this motion into the Newtonian-like equations of motion one gets the quasi-Newtonian restricted quasi-circular three-body problem sufficient to take into account the main indirect perturbations caused by the binary radiation terms. Such problem admits the Lagrange-like quasi-libration solutions and rather simple quasi-circular orbits lying at large distance from the binary.     

Innovative methods of correlation and orbit determination for space debris

We propose two algorithms to provide a full preliminary orbit of an Earth-orbiting object with a number of observations lower than the classical methods, such as those by Laplace and Gauss. The first one is the Virtual debris algorithm, based upon the admissible region, that is the set of the unknown quantities corresponding to possible orbits for a given observation for objects in Earth orbit (as opposed to both interplanetary orbits and ballistic ones). A similar method has already been successfully used in recent years for the asteroidal case. The second algorithm uses the integrals of the geocentric 2-body motion, which must have the same values at the times of the different observations for a common orbit to exist. We also discuss how to account for the perturbations of the 2-body motion, e.g., the J ₂ effect.     

天文地球动力学中的相对论效应

天文地球动力学利用空间与地面观测手段 ,监测和研究地球整体与各圈层的物质运动以及它们间的相互作用 ,这都离不开广义相对论涉及的时间与空间。随着空间对地观测精度的提高 ,为了充分利用高精度观测提供的信息 ,在天文地球动力学的研究中必须考虑相对论效应。所涉及的相对论效应包括 :( 1 )相对论参考系的建立 ,( 2 )在恰当的参考系中对观测者和被观测对象的相对论运动方程 (平动和自转 )的描述 ,( 3 )观测者和被观测对象间的电磁信号传播 ,( 4 )依赖于坐标选择的结果与具有物理意义的可观测量间的转换 ,( 5)某些基本概念与定义在广义相对论框架下的重新确认。本文对天文地球动力学中的这些相对论效应作了简要的评述。     

Toward the construction of a new precession-nutation theory of nonrigid Earth

The accuracy of the rigid Earth solution SMART97 is 2?μ as over the time interval (1968, 2023), accuracy showed by the comparison with a numerical integration using the positions of the Moon, the Sun, and the planets given by DE403. To obtain a nonrigid Earth solution, we use the transfer function of Mathews et al. (2000) and , to keep the precision of our rigid Earth solution in the computation of the geophysical effects, we apply this transfer function to the Earth's angular velocity vector in order to avoid the inherent approximations of the classical methods. Moreover the perturbations of the third component of the angular velocity vector are taken into account. Lastly, we take into account, in an iterative process, the second order perturbations due to the geophysical effects. The results are compared with the Herring solution (1996) published in the IERS Conventions.