Two‐body system with a rotating central body (e.g. Earth‐Moon system) within the Projective Unified Field Theory

In this treatise the well‐known 2‐body problem with a rotating central body is systematically reinvestigated on the basis of the Projective Unified Field Theory (PUFT) under the following aspects (including the special case of the Newton mechanics): First, equation of motion with abstract additional terms being appropriate for the interpretation of the various effects under discussion: tidal friction effect as well as non‐tidal effects (e.g. rebound effect as temporal variation of the moment of inertia of the rotating body, general‐relativistic Lense‐Thirring effect, new scalaric effects of cosmological origin, being an outcome of the scalarity phenomenon of matter (PUFT). Second, numerical evaluation of the theory. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Origin and evolution of the earth-moon system

The origin and evolution of the Earth-Moon system is studied by comparing it to the satellite systems of other planets. The normal structure of a system of secondary bodies orbiting around a central body depends essentially on the mass of the central body. The Earth with a mass intermediate between Uranus and Mars should have a normal satellite system that consists of about half a dozen satellites each with a mass of a fraction of a percent of the lunar mass. Hence, the Moon is not likely to have been generated in the environment of the Earth by a normal accretion process as is claimed by some authors.Capture of satellites is quite a common process as shown by the fact that there are six satellites in the solar system which, because they are retrograde, must have been captured. There is little doubt that the Moon is also a captured satellite, but its capture orbit and tidal evolution are still incompletely understood.The Earth and the Moon are likely to have been formed from planetesimals accreting in particle swarms in Kepler orbits (jet streams). This process leads to the formation of a cool lunar interior with an outer layer accreted at increasingly higher temperatures. The primeval Earth should similarly have formed with a cool inner core surrounded in this case by a very strongly heated outer core and with a mantle accreted slowly and with a low average temperature but with intense transient heating at each individual impact site.     

Control of chaos in the vicinity of the Earth‐Moon L5 Lagrangian point to keep a spacecraft in orbit

The concept of Space Manifold Dynamics is a new method of space research. We have applied it along with the basic idea of the method of Ott, Grebogi, and York (OGY method) to stabilize the motion of a spacecraft around the triangular Lagrange point L₅ of the Earth‐Moon system. We have determined the escape rate of the trajectories in the general three‐ and four‐body problem and estimated the average lifetime of the particles. Integrating the two models we mapped in detail the phase space around the L₅ point of the Earth‐Moon system. Using the phase space portrait our next goal was to apply a modified OGY method to keep a spacecraft close to the vicinity of L₅. We modified the equation of motions with the addition of a time dependent force to the motion of the spacecraft. In our orbit‐keeping procedure there are three free parameters: (i) the magnitude of the thrust, (ii) the start time, and (iii) the length of the control. Based on our numerical experiments we were able to determine possible values for these parameters and successfully apply a control phase to a spacecraft to keep it on orbit around L₅. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Current bombardment of the Earth-Moon system: Emphasis on cratering asymmetries

We calculate the current spatial distribution of projectile delivery to the Earth and Moon using numerical orbital dynamics simulations of candidate impactors drawn from a debiased Near-Earth Object (NEO) model. We examine the latitude distribution of impactor sites and find that for both the Earth and Moon there is a small deficiency of time-averaged impact rates at the poles. The ratio between deliveries within 30° of the pole to that of a 30° band centered on the equator is small for Earth (<5%) (0.958±0.001) and somewhat greater for the Moon (∼10%) (0.903±0.005). The terrestrial arrival results are examined to determine the degree of AM/PM asymmetry to compare with the PM excess shown in meteorite fall times. We find that the average lunar impact velocity is 20 km/s, which has ramifications in converting observed crater densities to impactor size distributions. We determine that current crater production on the leading hemisphere of the Moon is 1.28±0.01 that of the trailing when considering the ratio of craters within 30° of the apex to those within 30° of the antapex and that there is virtually no nearside-farside asymmetry, in agreement with observations of rayed craters. As expected, the degree of leading-trailing asymmetry increases when the Moon's orbital distance is decreased.     

Differential Rotation Driven by Precession

The tidal force in the Earth–Moon system exerted on the Earth's equatorial bulge results in the Earth's precession. It was proposed a long time ago that the strong shear flow driven by the precession of the Earth may power the Earth's dynamo in its liquid core. We present a nonlinear analytical study investigating how the Poincaré force in a rotating, precessing spherical system drives a large-amplitude differential rotation which plays a major role in the modern theory of the geodynamo. The analysis is based on a perturbation approach in terms of the small Poincaré force parameter. It is found that the amplitude of the precession-driven differential rotation is consistent with that estimated from the geomagnetic secular variation.     

月球探测器过渡轨道的短弧定轨方法

论述的短弧定轨,是指在无先验信息情况下又避开多变元迭代的初轨计算方法,它需要相应的动力学问题有一能反映短弧内达到一定精度的近似分析解．探测器进入月球引力作用范围后接近月球时可以处理成相对月球的受摄二体问题,而在地球附近,则可处理成相对地球的受摄二体问题,但在整个过渡段的力模型只能处理成一个受摄的限制性三体问题．而限制性三体问题无分析解,即使在月球引力作用范围外,对于大推力脉冲式的过渡方式,相对地球的变化椭圆轨道的偏心率很大(超过Laplace极限),在考虑月球引力摄动时亦无法构造摄动分析解．就此问题,考虑在地球非球形引力(只包含J2项)和月球引力共同作用下,构造了探测器飞抵月球过渡轨道段的时间幂级数解,在此基础上给出一种受摄二体问题意义下的初轨计算方法,经数值验证,定轨方法有效,可供地面测控系统参考．     

Did the moon form as a result of a thermonuclear explosion?

The proposed model explains the Moon formation as a result of a thermo-nuclear explosion due to which a big land mass was torn off from the Earth. Within the model framework, on the one hand, the data on the Moon’s physical and chemical parameters are in good agreement. On the other hand, this model corresponds to modern ideas about the dynamism of the Earth’s geological structure which presupposes the presence of a powerful energy source in the Earth’s core, which might have thermonuclear origin.     

Co-accretion of the Earth-Moon system after the giant impact: reduction of the total angular momentum by lunar impact ejecta

Though the Moon is considered to have been formed by the so-called giant impact, the mass of the Earth immediately after the impact is still controversial. If the Moon was formed during the Earth's accretion, a subsequent accretion of residual heliocentric planetesimals onto the protoearth and the protomoon must have occurred. In this co-accretion stage, a significant amount of lunar-impact-ejecta would be ejected to circumterrestrial orbits, since the mean impact velocity of the planetesimals with the protomoon is much larger than the escape velocity of the protomoon. Orbital calculations of test particles ejected from the protomoon, whose semimajor axis is smaller than that of the present Moon, reveal that most of the particles escaping from the protomoon also escape from the Hill sphere of the protoearth and reduce the planetocentric angular momentum of the primordial Earth-Moon system. Using the results of the ejecta simulations, we investigate the evolution of the mass ratio and the total angular momentum (Earth's spin angular momentum + Moon's orbital angular momentum) of the Earth-Moon system during the co-accretion. We find that the mass of the protomoon is almost constant or rather decreases and the total angular momentum decreases significantly, if the random velocity of planetesimals is as large as the escape velocity of the protoearth. On the other hand, if the random velocity is the half of the escape velocity of the protoearth, the mass ratio is kept to be almost as large as the present value and the decrease of the total angular momentum is not so significant. Comparing with the results of giant impact simulations, we find that the mass of the protoearth immediately after the Moon-forming impact was 0.7-0.8 times the present value if the impactor-to-target mass ratio was 3:7, whereas the giant impact occurred almost in the end of the Earth's accretion if the impactor-to-target mass ratio was 1:9.     

月球卫星轨道力学综述

月球探测器的运动通常可分为3个阶段，这3个阶段分别对应3种不同类型的轨道：近地停泊轨道、向月飞行的过渡轨道与环月飞行的月球卫星轨道。近地停泊轨道实为一种地球卫星轨道；过渡轨道则涉及不同的过渡方式(大推力或小推力等)；环月飞行的月球卫星轨道则与地球卫星轨道有很多不同之处，它决不是地球卫星轨道的简单克隆。针对这一点，全面阐述月球卫星的轨道力学问题，特别是环月飞行中的一些热点问题，如轨道摄动解的构造、近月点高度的下降及其涉及的卫星轨道寿命、各种特殊卫星(如太阳同步卫星和冻结轨道卫星等)的轨道特征、月球卫星定轨等。     

Analytical theory of a lunar artificial satellite with third body perturbations

We present here the first numerical results of our analytical theory of an artificial satellite of the Moon. The perturbation method used is the Lie Transform for averaging the Hamiltonian of the problem, in canonical variables: short-period terms (linked to l, the mean anomaly) are eliminated first. We achieved a quite complete averaged model with the main four perturbations, which are: the synchronous rotation of the Moon (rate ), the oblateness J ₂ of the Moon, the triaxiality C ₂₂ of the Moon ( ) and the major third body effect of the Earth (ELP2000). The solution is developed in powers of small factors linked to these perturbations up to second-order; the initial perturbations being sorted ( is first-order while the others are second-order). The results are obtained in a closed form, without any series developments in eccentricity nor inclination, so the solution apply for a wide range of values. Numerical integrations are performed in order to validate our analytical theory. The effect of each perturbation is presented progressively and separately as far as possible, in order to achieve a better understanding of the underlying mechanisms. We also highlight the important fact that it is necessary to adapt the initial conditions from averaged to osculating values in order to validate our averaged model dedicated to mission analysis purposes.     

Near-Earth object velocity distributions and consequences for the Chicxulub impactor

An Öpik-based geometric algorithm is used to compute impact probabilities and velocity distributions for various near-Earth object (NEO) populations. The resulting crater size distributions for the Earth and Moon are calculated by combining these distributions with assumed NEO size distributions and a selection of crater scaling laws. This crater probability distribution indicates that the largest craters on both the Earth and the Moon are dominated by comets. However, from a calculation of the fractional probabilities of iridium deposition, and the velocity distributions at impact of each NEO population, the only realistic possibilities for the Chicxulub impactor are a short-period comet (possibly inactive) or a near-Earth asteroid. For these classes of object, sufficiently large impacts have mean intervals of 100 and 300 Myr respectively, slightly favouring the cometary hypothesis.     

Oceanic and solid Earth angular momentum

We reconsider two hypotheses used in calculating the transfer of angular momentum between the oceans and the solid Earth: (1) The locked-cean-ypothesis was already given up some time ago; here we provide a simple manner of understanding the relative importance of the motion and matter term. (2) The isolation hypothesis implied the isolation of the whole Earth in short timescales with regard to angular momentum exchange, and consequently, the neglection of the exchange with the tide-enerating body. It is shown that for present accuracy requirements this exchange has to be taken into account.     

RDAN97: An Analytical Development of Rigid Earth Nutation Series Using the Torque Approach

New series of rigid Earth nutations for the angular momemtum axis, the rotation axis and the figure axis, named RDAN97, are computed using the torque approach. Besides the classical J2 terms coming from the Moon and the Sun, we also consider several additional effects: terms coming from J3 and J4 in the case of the Moon, direct and indirect planetary effects, lunar inequality, J2 tilt, planetary‐tilt, effects of the precession and nutations on the nutations, secular variations of the amplitudes, effects due to the triaxiality of the Earth, new additional out‐of‐phase terms coming from second order effect and relativistic effects. Finally, we obtain rigid Earth nutation series of 1529 terms in longitude and 984 terms in obliquity with a truncation level of 0.1 μ (microarcsecond) and 8 significant digits. The value of the dynamical flattening used in this theory is HD=(C-A)/C=0.0032737674 computed from the initial value pa=50′.2877/yr for the precession rate. These new rigid Earth nutation series are then compared with the most recent models (Hartmann et al., 1998; Souchay and Kinoshita, 1996, 1997; Bretagnon et al., 1997, 1998. We also compute a benchmark series (RDNN97) from the numerical ephemerides DE403/LE403 (Standish et al., 1995) in order to test our model. The comparison between our model (RDAN97) and the benchmark series (RDNN97) shows a maximum difference, in the time domain, of 69 μas in longitude and 29 μas in obliquity. In the frequency domain, the maximum differences are 6 μas in longitude and 4 μ as in obliquity which is below the level of precision of the most recent observations (0.2 mas in time domain (temporal resolution of 1 day) and 0.02 mas in frequency domain). This revised version was published online in July 2006 with corrections to the Cover Date.     

The evolution of the moon

The thermal evolution of the Moon as it can be defined by the available data and theoretical calculations is discussed. A wide assortment of geological, geochemical and geophysical data constrain both the present-day temperatures and the thermal history of the lunar interior. On the basis of these data, the Moon is characterized as a differentiated body with a crust, a 1000-km-thick solid mantle (lithosphere) and an interior region (core) which may be partially molten. The presence of a crust indicates extensive melting and differentiation early in the lunar history. The ages of lunar samples define the chronology of igneous activity on the lunar surface. This covers a time span of about 1.5 billion yr, from the origin to about 3.16 billion yr ago. Most theoretical models require extensive melting early in the lunar history, and the outward differentiation of radioactive heat sources.Thermal history calculations, whether based on conductive or convective computation codes define relatively narrow bounds for the present day temperatures in the lunar mantle. In the inner region of the 700 km radius, the temperature limits are wider and are between about 100 and 1600°C at the center of the Moon. This central region could have a partially or totally molten core.The lunar heat flow values (about 30 ergs/cm²s) restrict the present day average uranium abundance to 60 ± 15 ppb (averaged for the whole Moon) with typical ratios of K/U = 2000 and Th/U = 3.5. This is consistent with an achondritic bulk composition for the Moon.The Moon, because of its smaller size, evolved rapidly as compared to the Earth and Mars. The lunar interior is cooling everywhere at the present and the Moon is tectonically inactive while Mars could be and the Earth is definitely active.     

On the fissions of a solid body under influence of tidal force; with application to the problem of twin craters on the moon

The internal strain due to the tidal force in the proximity of a tide-generating body (in the present case, the Moon) is calculated according to the Lord Kelvin theory of Earth tides. The conditions for which uniform elastic sphere possessing a definite tensile strength is crushed near the surface of the Moon is investigated. The state of internal stress is almost independent of the value of elastic constants. Many lunar features, such as twin craters, craterous walled plains of irregular forms, compound craters, may be explained by fission of the meteoritic material before impact.     

Direct and indirect capture of near-Earth asteroids in the Earth–Moon system

Near-Earth asteroids have attracted attention for both scientific and commercial mission applications. Due to the fact that the Earth–Moon \(\hbox {L}_{1}\) and \(\hbox {L}_{2}\) points are candidates for gateway stations for lunar exploration, and an ideal location for space science, capturing asteroids and inserting them into periodic orbits around these points is of significant interest for the future. In this paper, we define a new type of lunar asteroid capture, termed direct capture. In this capture strategy, the candidate asteroid leaves its heliocentric orbit after an initial impulse, with its dynamics modeled using the Sun–Earth–Moon restricted four-body problem until its insertion, with a second impulse, onto the \(\hbox {L}_{2}\) stable manifold in the Earth–Moon circular restricted three-body problem. A Lambert arc in the Sun-asteroid two-body problem is used as an initial guess and a differential corrector used to generate the transfer trajectory from the asteroid’s initial obit to the stable manifold associated with Earth–Moon \(\hbox {L}_{2}\) point. Results show that the direct asteroid capture strategy needs a shorter flight time compared to an indirect asteroid capture, which couples capture in the Sun–Earth circular restricted three-body problem and subsequent transfer to the Earth–Moon circular restricted three-body problem. Finally, the direct and indirect asteroid capture strategies are also applied to consider capture of asteroids at the triangular libration points in the Earth–Moon system.     

Effect of Moon perturbation on the energy curves and equilibrium points in the Sun–Earth–Moon system

In this paper, we have considered that the Moon motion around the Earth is a source of a perturbation for the infinitesimal body motion in the Sun–Earth system. The perturbation effect is analyzed by using the Sun–Earth–Moon bi–circular model (BCM). We have determined the effect of this perturbation on the Lagrangian points and zero velocity curves. We have obtained the motion of infinitesimal body in the neighborhood of the equivalent equilibria of the triangular equilibrium points. Moreover, to know the nature of the trajectory, we have estimated the first order Lyapunov characteristic exponents of the trajectory emanating from the vicinity of the triangular equilibrium point in the proposed system. It is noticed that due to the generated perturbation by the Moon motion, the results are affected significantly, and the Jacobian constant is fluctuated periodically as the Moon is moving around the Earth. Finally, we emphasize that this model could be applicable to send either satellite or telescope for deep space exploration.     

Study of the transfer between libration point orbits and lunar orbits in Earth–Moon system

This paper is devoted to the study of the transfer problem from a libration point orbit of the Earth–Moon system to an orbit around the Moon. The transfer procedure analysed has two legs: the first one is an orbit of the unstable manifold of the libration orbit and the second one is a transfer orbit between a certain point on the manifold and the final lunar orbit. There are only two manoeuvres involved in the method and they are applied at the beginning and at the end of the second leg. Although the numerical results given in this paper correspond to transfers between halo orbits around the \(L_1\) point (of several amplitudes) and lunar polar orbits with altitudes varying between 100 and 500 km, the procedure we develop can be applied to any kind of lunar orbits, libration orbits around the \(L_1\) or \(L_2\) points of the Earth–Moon system, or to other similar cases with different values of the mass ratio.     

Analytical methods for the orbits of artificial satellites of the Moon

The motion of a close artificial satellite of the Moon is considered. The principal perturbations taken into account are caused by the nonsphericity of the Moon and the attraction of the Earth and the Sun. To begin with, the expansions of the disturbing functions due to the nonsphericity of the primary body and the action of the disturbing mass-point body have been derived. The second expansion is produced in terms of the Keplerian elements of a satellite and the spherical coordinates of the disturbing body. Both expansions are valid for an arbitrary reference plane. The motion of a satellite of the Moon is studied in the selenocentric coordinate system referred to the Lunar equator and rotating with respect to the fixed ecliptic system. However, the coordinate exes in the equatorial plane are chosen so that the angular speed of rotation of the system is small. The motion of the satellite is described by means of the contact elements which enable one to utilize the conventional Lagrange's planetary equations and may be regarded as the generalization of the notion of the osculating elements to the case of the disturbing function depending not only o the coordinates and the time but on the velocities as well. Two methods are proposed to represent the motion of Lunar satellites over long intervals of time: the von Zeipel method and the Euler method of analytical integration with application of the variation-of-elements technique at every step of integration. The second method is exposed in great detail.Presented at the Meeting of Commission 7 of the IAU on Analytical Methods for the Orbits of Artificial Celestial Objects 14-th General Assembly of the IAU, Brighton, 1970.