Numerical investigation of collision orbits of lunar satellites

In the present study an investigation of the collision orbits of natural satellites of the Moon (considered to be of finite dimensions) is developed, and the tendency of natural satellites of the Moon to collide on the visible or the far side of the Moon is studied. The collision course of the satellite is studied up to its impact on the lunar surface for perturbations of its initial orbit arbitrarily induced, for example, by the explosion of a meteorite. Several initial conditions regarding the position of the satellite to collide with the Moon on its near (visible) or far (invisible) side is examined in connection to the initial conditions and the direction of the motion of the satellite. The distribution of the lunar craters-originating impact of lunar satellites or celestial bodies which followed a course around the Moon and lost their stability - is examined. First, we consider the planar motion of the natural satellite and its collision on the Moon's surface without the presence of the Earth and Sun. The initial velocities of the satellite are determined in such a way so its impact on the lunar surface takes place on the visible side of the Moon. Then, we continue imparting these velocities to the satellite, but now in the presence of the Earth and Sun; and study the forementioned impacts of the satellites but now in the Earth-Moon-Satellite system influenced also by the Sun. The initial distances of the satellite are taken as the distances which have been used to compute periodic orbits in the planar restricted three-body problem (cf. Gousidou-Koutita, 1980) and its direction takes different angles with the x-axis (Earth-Moon axis). Finally, we summarise the tendency of the satellite's impact on the visible or invisible side of the Moon.     

An averaging method to study the motion of lunar artificial satellites I: Disturbing Function

We describe a semi-analytical averaging method aimed at the computation of the motion of an artificial satellite of the Moon. In this paper, the first of the two part study, we expand the disturbing function with respect to the small parameters. In particular, a semi-analytic theory of the motion of the Moon around the Earth and the libration of the lunar equatorial plane using different reference frames are introduced. The second part of this article shows that the choice of the canonical Poincaré variables lead to equations in closed form without singularities in e = 0 or I = 0. We introduce new expressions that are sufficiently compact to be used for the study of any artificial satellite. This revised version was published online in July 2006 with corrections to the Cover Date.     

An averaging method to study the motion of lunar artificial satellites II: Averaging and Applications

On this, the second part of a two part study (Steichen, 1998) we further develop a semi-analytical theory for a lunar artificial satellite. This theory is obtained by averaging analytically the Hamiltonian function over period up to a month. The averaged equations are then numerically integrated. The solution is free from singularities at e = 0 and I = 0 and is not expanded in powers of these variables. In the last section, the analytic work is applied to characteristic examples to validate the method used. This revised version was published online in July 2006 with corrections to the Cover Date.     

Quasi-critical orbits for artificial lunar satellites

We study the problem of critical inclination orbits for artificial lunar satellites, when in the lunar potential we include, besides the Keplerian term, the J ₂ and C ₂₂ terms and lunar rotation. We show that, at the fixed points of the 1-D averaged Hamiltonian, the inclination and the argument of pericenter do not remain both constant at the same time, as is the case when only the J ₂ term is taken into account. Instead, there exist quasi-critical solutions, for which the argument of pericenter librates around a constant value. These solutions are represented by smooth curves in phase space, which determine the dependence of the quasi-critical inclination on the initial nodal phase. The amplitude of libration of both argument of pericenter and inclination would be quite large for a non-rotating Moon, but is reduced to <0°.1 for both quantities, when a uniform rotation of the Moon is taken into account. The values of J ₂, C ₂₂ and the rotation rate strongly affect the quasi-critical inclination and the libration amplitude of the argument of pericenter. Examples for other celestial bodies are given, showing the dependence of the results on J ₂, C ₂₂ and rotation rate.     

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.     

The motion of the Martian satellites

An extensive analysis of the motion of Phobos and Deimos from 1877 to 1973 has been fulfilled. The new values of the parameters of the orbital model first developed by Struve have been determined for both satellites. The new sets of the orbital parameters compete with the solutions of similar accuracy found by Wilkins and Sinclair. A secular acceleration in longitude of Phobos is found to be equal to +(0.107±0.011)×10^?7 deg day^?2. The value of the acceleration is little affected when one or another group of oppositions is omitted. The acceleration of Deimos is determined with great uncertainty: +(0.06±0.34)×10^?9 deg day^?2. Values found for the orbital parameters seem to be in good agreement since the mass, oblateness and coordinates of the pole of Mars inferred from the motion of each satellite have similar values in both cases.     

A semi-analytic theory for the motion of a lunar satellite

A semi-analytical solution to the problem of the motion of a satellite of the moon is presented. Perturbative effects which are considered include those due to the attraction of the moon, earth, and sun, the non-sphericity of the moon's gravitational field, coupling of lower-order terms, solar radiation pressure, and physical libration. Short-period terms and intermediate-period terms, terms with the period of the moon's longitude, are produced by means of von Zeipel's method; it is proposed to obtain the secular perturbations, and those depending only on the argument of perilune, by numerical integration of the equations of motions. The short-period terms and intermediate-period terms are developed up to second order, where first order is 10^–2. The secular perturbations and perturbations dependent on the argument of perilune are obtained to third order.     

Relative motion of near orbiting satellites

The relative motion of two particles on adjacent orbits about the same primary has been investigated under the condition that both motions have the same period. The geometrical properties of the relative displacement and velocity traces, on representative planes, are studied. A complete state of the motion is given; and, the range and range-rate variations, over one or more orbits, are described.It has been found that cusps appear on some of the traces provided that a proper relationship exists between the eccentricity and inclination. (Here, one particle moves on a circular path while the second moves on an ellipse.) The conditions for which cusps appear are given, and typical traces are shown.     

Two problems about the motion of low-moon-orbit satellites

A detailed theoretical analysis on the orbital lifetime and frozen orbit of low-moon-orbit satellites (LMOS) is carried out, and their relationships with the orbital inclination, as well as some mutual relationships are presented. Taking account of the main perturbing sources of low-orbit satellites, we carried out numerical simulations under a comprehensive force model, and the results not only confirm the correctness of the theoretical analysis, but also provide some valuable insights on the orbital design of LMOS.     

Application of two special orbits in the orbit determination of lunar satellites

Using inter-satellite range data,the combined autonomous orbit determination problem of a lunar satellite and a probe on some special orbits is studied in this paper.The problem is firstly studied in the circular restricted three-body problem,and then generalized to the real force model of the Earth-Moon system.Two kinds of special orbits are discussed:collinear libration point orbits and distant retrograde orbits.Studies show that the orbit determination accuracy in both cases can reach that of the observations.Some important properties of the system are carefully studied.These findings should be useful in the future engineering implementation of this conceptual study.     

On analytic modeling of lunar perturbations of artificial satellites of the earth

Two different procedures for analytically modeling the effects of the Moon's direct gravitational force on artificial Earth satellites are discussed from theoretical and numerical viewpoints. One is developed using classical series expansions of inclination and eccentricity for both the satellite and the Moon, and the other employs a method of averaging. Both solutions are seen to have advantages, but it is shown that while the former can be more accurate in special situations, the latter is quicker and more practical for the general orbit determination problem where observed data is used to correct the orbit in near real time.This work was sponsored with the support of the Department of the Air Force under contract F19628-85-C-0002. The views expressed are those of the author and do not reflect the official policy or position of the US Government.     

The motion of a lunar satellite

Presented in this theory is a semianalytical solution for the problem of the motion of a satellite in orbit around the moon. The principal perturbations on such a body are due to the nonspherical gravity field of the moon, the attraction of the earth, and, to a lesser degree, the attraction of the sun. The major part of the problem is solved by means of the celebrated von Zeipel Method, first successfully applied to the motion of an artificial earth satellite by Brouwer in 1959. After eliminating from the Hamiltonian all terms with the period of the satellite and those with the period of the moon, it is suggested to solve the remaining problem with the aid of numerical integration of the modified equations of motion.This theory was written in 1964 and presented as a dissertation to Yale University in 1965. Since then a great deal has been learned about the gravity field of the moon. It seems that quite a number of recently determined gravity coefficients would qualify as small quantities of order two. Hence, according to the truncation criteria employed, they should be considered in the present theory. However, the author has not endeavored to update the work accordingly. The final results, therefore, are incomplete in the lunar gravitational perturbations. Nevertheless, the theory does give the largest such variations and it does present the methods by which perturbations may be derived for any gravity terms not actually developed.     

Analytic short period lunar and solar perturbations of artificial satellites

The short period luni-solar theory of Kozai is generalized for arbitrary obliquity of the ecliptic and inclination of the moon's orbit to the ecliptic. Analytic first order lunar perturbations to the elements are derived. The theory is illustrated by an application to the communication satellite Intelsat 3F3.Presently at the Department of Environmental Sciences, University of Tel Aviv, Ramat Aviv, Israel.     

On the tidal effects in the motion of artificial satellites

The general perturbations in the elliptic and vectorial elements of a satellite as caused by the tidal deformations of the non-spherical Earth are developed into trigonometric series in the standard ecliptical arguments of Hill-Brown lunar theory and in the equatorinal elements ω and Ω of the satellite. The integration of the differential equations for variation of elements of the satellite in this theory is easy because all arguments are linear or nearly linear in time. The trigonometrical expansion permits a judgment about the relative significance of the amplitudes and periods of different tidal ‘waves’ over a long period of time. Graphs are presented of the tidal perturbations in the elliptic elements of the BE-C satellite which illustrate long term periodic behavior. The tidal effects are clearly noticeable in the observations and their comparison with the theory permits improvement of the ‘global’ Love numbers for the Earth.     

Insight-building models for lunar range and range rate

The analysis of range or Doppler data between sites on the Earth and Moon requires an accurate computation of the lunar orbit and detailed models of the orientation of the Earth and Moon. Models constructed to understand range and range rate can lack detail, but if they include the largest lunar orbit variations, tracking stations on a rotating Earth, and lunar sites on a synchronously rotating Moon, then they will display the largest effects for orbit elements, Earth orientation, tracking station locations, and lunar site coordinates. The range and range rate are expanded into periodic series. To understand accurate solutions, the largest periodic terms that are sensitive to various solution parameters indicate the sensitivity of data to solution parameters and the time spans needed for their determination. Conclusions include: cylindrical coordinates work well for sites on the rapidly rotating Earth, but Cartesian coordinates are more natural for the synchronously rotating Moon since the series for the three coordinate projections are distinct. For range and range rate data, daily, semimonthly, monthly, and longer periods are present. For Doppler data, the daily periods may be stronger and more useful than the long periods, particularly for terms associated with the terrestrial tracking station. Doppler data do not determine the lander coordinate toward the Earth well. Observational strategies for range and Doppler data are not identical. For all data types, one wishes a variety of hour angles, lunar declinations, times of month, and longer periods. A long span of high-quality range data can improve the lunar orbit, orientation of the Earth’s equator, and physical librations. The locations of new lunar sites or new tracking stations can be determined from shorter spans of data.     

Influence of geomagnetic field on the motion of low satellites

The first aim of the present work is to compute a more accurate and recent model for the Earth’s magnetic field. The second aim is to determine the effects of the Earth’s magnetic field on the motion of a charged artificial satellite to evaluate the variations of the orbital elements of the satellite due to these effects. The magnetic field and its variation with time have been studied at different heights, longitudes and latitudes. The geomagnetic field is considered as a multiple potential field and the electrical charge of the satellite is assumed to be constant. A new computer code has been constructed to follow the components of the magnetic field in spherical harmonic models. The Gauss equations are solved numerically. The results concentrate on the computation of the numerical values of orbital perturbation for the case of a low Earth satellite. RS-1 satellite and space craft gravity probe B (GPB) are chosen as cases of studies for a detailed numerical analysis.     

Morphology of lunar craters: A test of lunar erosional models

This paper is a test of published theoretical and experimental studies of crater erosion by micrometeorite bombardment which predict systematic variations in the morphology of lunar craters as a function of crater diameter and crater age. Numerical, ranking-type degradation classifications indicate that the craters on Mare Imbrium and Mare Tranquillitatus confirm these predictions by showing a systematic increase in degradation with decreasing diameter for craters smaller than a few kilometers in diameter but larger than the equilibrium diameter, and by showing fixed proportions of fresh, moderately degraded and very degraded craters under equilibrium conditions. Furthermore, the relative ages of the two mare surfaces may be determined using a diameter/mean-degradation-number curve. These determinations of relative age and process of crater erosion are both essentially independent of the traditionally studied crater diameter/frequency relationships. Morphologies of terra craters near Mare Humorum suggest a young, non-equilibrium crater population superposed on a perimordial population with about equilibrium proportions of fresh, moderately degraded and very degraded craters. The primordial population has been modified by pre-Imbrian or early Imbrian deposition of blanketing deposits. A comparative study of several crater degradation classifications indicates that all are essentially interchangeable.     

On a semi-inverse problem in the motion of gyrostat satellites

In this paper we study the equilibrium orientation of a gyrostat satellite in the gravity field of a point mass. Direct problem is to find all possible equilibrium orientation when the relative angular momentum vector is given. Inverse problem is to find this relative angular momentum in order to obtain equilibrium in a given orientation. Semi-inverse problem is solved here when some parameters (but not all) giving orientation of the satellite are chosen arbitrarily, giving for what choices real solutions occur.