Design of transfer trajectories between resonant orbits in the Earth–Moon restricted problem

The application of dynamical systems techniques to mission design has demonstrated that employing invariant manifolds and resonant flybys enables previously unknown trajectory options and potentially reduces the ΔV$\mathrm{\Delta }V$ requirements. In this investigation, planar and three-dimensional resonant orbits are analyzed and cataloged in the Earth–Moon system and the associated invariant manifold structures are computed and visualized with the aid of higher-dimensional Poincaré maps. The relationship between the manifold trajectories associated with multiple resonant orbits is explored through the maps with the objective of constructing resonant transfer arcs. As a result, planar and three-dimensional homoclinic- and heteroclinic-type trajectories between unstable periodic resonant orbits are identified in the Earth–Moon system. To further illustrate the applicability of 2D and 3D resonant orbits in preliminary trajectory design, planar transfers to the vicinity of L₅ and an out-of-plane transfer to a 3D periodic orbit, one that tours the entire Earth–Moon system, are constructed. The design process exploits the invariant manifolds associated with orbits in resonance with the Moon as transfer mechanisms.     

限制性四体问题中的时间相关不变流形

借助有限时间Lyapunov指数（FTLE）定义了拉格朗日拟序结构（LCS）,并将LCS作为不变流形的替代物。针对日-地-月双圆模型（BCM）,利用LCS研究了限制性四体问题（R4BP）中的时间相关不变流形（TDIM）的性质。采用数值方法验证了TDIM是运动分界面和轨道不变集。继而,利用二分法对给定Poincare截面上的LCS进行了精确提取,通过一系列等能量面上的LCS描绘出TDIM在给定截面上的构形。最后,借助TDIM,初步研究了低能奔月轨道在非自治系统BCM中的直接构建。     

Manifold dynamics in the Earth–Moon system via isomorphic mapping with application to spacecraft end-of-life strategies

Recently, manifold dynamics has assumed an increasing relevance for analysis and design of low-energy missions, both in the Earth–Moon system and in alternative multibody environments. With regard to lunar missions, exterior and interior transfers, based on the transit through the regions where the collinear libration points L₁ and L₂ are located, have been studied for a long time and some space missions have already taken advantage of the results of these studies. This paper is focused on the definition and use of a special isomorphic mapping for low-energy mission analysis. A convenient set of cylindrical coordinates is employed to describe the spacecraft dynamics (i.e. position and velocity), in the context of the circular restricted three-body problem, used to model the spacecraft motion in the Earth–Moon system. This isomorphic mapping of trajectories allows the identification and intuitive representation of periodic orbits and of the related invariant manifolds, which correspond to tubes that emanate from the curve associated with the periodic orbit. Heteroclinic connections, i.e. the trajectories that belong to both the stable and the unstable manifolds of two distinct periodic orbits, can be easily detected by means of this representation. This paper illustrates the use of isomorphic mapping for finding (a) periodic orbits, (b) heteroclinic connections between trajectories emanating from two Lyapunov orbits, the first at L₁, and the second at L₂, and (c) heteroclinic connections between trajectories emanating from the Lyapunov orbit at L₁ and from a particular unstable lunar orbit. Heteroclinic trajectories are asymptotic trajectories that travels at zero-propellant cost. In practical situations, a modest delta-v budget is required to perform transfers along the manifolds. This circumstance implies the possibility of performing complex missions, by combining different types of trajectory arcs belonging to the manifolds. This work studies also the possible application of manifold dynamics to defining suitable, convenient end-of-life strategies for spacecraft orbiting the Earth. Seven distinct options are identified, and lead to placing the spacecraft into the final disposal orbit, which is either (a) a lunar capture orbit, (b) a lunar impact trajectory, (c) a stable lunar periodic orbit, or (d) an outer orbit, never approaching the Earth or the Moon. Two remarkable properties that relate the velocity variations with the spacecraft energy are employed for the purpose of identifying the optimal locations, magnitudes, and directions of the velocity impulses needed to perform the seven transfer trajectories. The overall performance of each end-of-life strategy is evaluated in terms of time of flight and propellant budget.     

双三体系统不变流形拼接成的低成本探月轨道

传统的探月轨道设计原理为二体模型框架下的Hohmann变轨理论,但1991年日本的Hiten探月器利用太阳的摄动,用比传统的方法更少的燃料完成了探月任务。利用三体问题非线性系统的不变流形设计了节省燃料的探月轨道。沿用JPL研究组的思路,将太阳-地球-月亮-航天器四体问题分解成太阳-地球-航天器和地球-月亮-航天器两个共面的圆形限制性三体问题,对Hiten类的探月轨道给出了更深刻的数学、力学解释;给出了流形的结构以及更合理的拼接方式;找到了发射位置、发射速度和拼接点;设计出了类似Hiten探月器的探月轨道,可比传统方法节省速度增量12%左右。结果证明了三体系统不变流形在登月轨道设计研究中的可行性和优越性。     

基于拉格朗日方法的通气云状空泡三维非定常脱落特性研究

采用均相流模型并结合FBM湍流模型,对绕轴对称回转体通气云状空泡流动特性进行了三维数值模拟,基于实验结果对数值方法进行验证,同时利用基于拉格朗日体系的有限时间李雅普诺夫指数(FTLE)、拉格朗日拟序结构(LCS)和粒子追踪方法分析了其三维非定常脱落特性。研究结果表明:纵截面上空泡覆盖区域的拉格朗日拟序结构整体呈椭球状分布,内部为不规律的复杂拟序结构;不同横截面上拟序结构分布存在很大差异。空泡内部的非对称流动结构和周向流动导致空泡呈现很强的三维运动特性。反向射流在周向上推进的不同步性,是造成空泡呈现不规则断裂和大尺度U型空泡团脱落的主要原因。     

Earth–Moon libration point orbit stationkeeping: Theory,modeling, and operations

Collinear Earth–Moon libration points have emerged as locations with immediate applications. These libration point orbits are inherently unstable and must be maintained regularly which constrains operations and maneuver locations. Stationkeeping is challenging due to relatively short time scales for divergence, effects of large orbital eccentricity of the secondary body, and third-body perturbations. Using the Acceleration Reconnection and Turbulence and Electrodynamics of the Moon's Interaction with the Sun (ARTEMIS) mission orbit as a platform, the fundamental behavior of the trajectories is explored using Poincaré maps in the circular restricted three-body problem. Operational stationkeeping results obtained using the Optimal Continuation Strategy are presented and compared to orbit stability information generated from mode analysis based in dynamical systems theory.     

Low-energy Earth–Moon transfers involving manifolds through isomorphic mapping

Analysis and design of low-energy transfers to the Moon has been a subject of great interest for decades. Exterior and interior transfers, based on the transit through the regions where the collinear libration points are located, have been studied for a long time and some space missions have already taken advantage of the results of these studies. This paper is concerned with a geometrical approach for low-energy Earth-to-Moon mission analysis, based on isomorphic mapping. The isomorphic mapping of trajectories allows a visual, intuitive representation of periodic orbits and of the related invariant manifolds, which correspond to tubes that emanate from the curve associated with the periodic orbit. Two types of Earth-to-Moon missions are considered. The first mission is composed of the following arcs: (i) transfer trajectory from a circular low Earth orbit to the stable invariant manifold associated with the Lyapunov orbit at L₁ (corresponding to a specified energy level) and (ii) transfer trajectory along the unstable manifold associated with the Lyapunov orbit at L₁, with final injection in a periodic orbit around the Moon. The second mission is composed of the following arcs: (i) transfer trajectory from a circular low Earth orbit to the stable invariant manifold associated with the Lyapunov orbit at L₁ (corresponding to a specified energy level) and (ii) transfer trajectory along the unstable manifold associated with the Lyapunov orbit at L₁, with final injection in a capture (non-periodic) orbit around the Moon. In both cases three velocity impulses are needed to perform the transfer: the first at an unknown initial point along the low Earth orbit, the second at injection on the stable manifold, the third at injection in the final (periodic or capture) orbit. The final goal is in finding the optimization parameters, which are represented by the locations, directions, and magnitudes of the velocity impulses such that the overall delta-v of the transfer is minimized. This work proves how isomorphic mapping (in two distinct forms) can be profitably employed to optimize such transfers, by determining in a geometrical fashion the desired optimization parameters that minimize the delta-v budget required to perform the transfer.     

基于地月信息的奔月探测器自主光学导航

基于地月信息提出了一种奔月转移轨道的自主光学导航算法。该算法首先利用导航相机得到的地球和月球边缘图像与地、月几何形状模型，运用扩展源的空间获取算法导出地心和月心对应的像素，然后利用地心和月心像素、探测器的惯性姿态和该观测历元的地、月星历信息，通过基于UD分解的递推加权最小二乘算法确定了探测器轨道。还引入导航系统的可观度定义来分析了导航系统的可观性，数学仿真结果表明导航的位置误差小于30km，速度误差小于0．3m／s，证明该自主光学导航算法是有效的。     

Interplanetary CubeSats system for space weather evaluations and technology demonstration

The paper deals with the mission analysis and conceptual design of an interplanetary 6U CubeSats system to be implemented in the L₁ Earth–Sun Lagrangian Point mission for solar observation and in-situ space weather measurements.     

地月平动点拟周期轨道探测器自主导航方法研究

地月L2点的拟周期轨道可以用于实现与月球背面的持续通信,具有重要的科学研究价值和广阔的应用前景。针对地月L2点探测器所处的弱稳定拟周期轨道,论证了基于日地月信息的自主导航方法的可行性,并进行了深入分析。首先,推导了会合坐标系下带有星历的精确导航动力学方程;其次,针对弱稳定轨道不同于近地强稳定轨道的特性,在基于日地月方位信息导航的基础上,提出了三种敏感器组合方案。使用迭代最小二乘方法给出导航仿真结果,并结合非线性可辨识性理论对这三种情况下历元状态的可辨识性及可辨识度进行分析。最后,仿真结果表明,使用日地月敏感器以及对地多普勒雷达可以满足历元状态的可辨识性、导航的收敛性以及系统经济性的要求。     

Potential evolution of an international moon base programme

The Moon is a major target in expanding human activity in Space. President Bush has called for a Space Exploration Initiative. European participation may depend on achieving an affordable programme and identifying distinct elements for non-U.S. participation. Affordability requires that all participants can influence the “cost to user” of Base operations. If lunar activity is to evolve towards resource exploitation, there will need to be a progressive reduction in operating costs. European interest would prefer participation that allowed longer-term independent interests. The paper addresses how non-U.S. agencies could contribute valuable elements to an International Moon Base while meeting three criteria:

• — Keep a core infrastructure under U.S. control.

• — Avoid a total reliance by the partner on U.S. services.

• — Allow the partner to evolve towards an eventual, semi-autonomous or autonomous capability.

The paper illustrates possible implications of meeting these constraints through “mini infrastructures” combining several elements to form a working architecture. It is concluded that any European participation in an International Moon Base Programme should contain both Space transport and surface elements.     

Motion of a Planet with a Complex Structure

The motion of a planet consisting of a mantle and a core (solid bodies) connected by a viscoelastic layer and interacting with each other and an external point mass by the law of gravitation is considered. The mutual motions of the core and mantle are investigated assuming that the centers of mass of the planet and external point mass moves along undisturbed Keplerian orbits around the common center of mass of the system. The planetary core and mantle have an axial symmetry and different principal moments of inertia, which leads to a displacement of the center of mantle relative to the center of core and to their mutual rotations. The results obtained on the basis of averaged equations are illustrated by the example of the Earth–Moon system.     

Evolution of almost circular orbits of satellites under the action of noncentral gravitational field of the Earth and lunisolar perturbations

Based on the results of paper [1] by G.V. Mozhaev, joint perturbations produced by nonsphericity of the Earth and by attraction of the Moon and the Sun are investigated using the method of averaging. Arbitrary number of spherical harmonics was taken into account in the force function of the Earth’s gravitational filed, and only the principal term was retained in the perturbing function of the Sun. In the perturbing function of the Moon two parallactic terms were considered in addition to the dominant term. The flight altitude was chosen in such a way that perturbations produced by the Sun and Moon would have the second order of smallness relative to the polar oblateness of the Earth. As a result, the formulas for calculation of satellite coordinates are derived that give a high precision on long time intervals.     

A gateway for human exploration of space? The weak stability boundary

NASA's plans for future human exploration of the Solar System describe only missions to Mars. Before such missions can be initiated, much study remains to be done in technology development, mission operations and human performance. While, for example, technology validation and operational experience could be gained in the context of lunar exploration missions, a NASA lunar program is seen as a competitor to a Mars mission rather than a step towards it. The recently characterized weak stability boundary in the Earth–Moon gravitational field may provide an operational approach to all types of planetary exploration, and infrastructure developed for a gateway to the Solar System may be a programmatic solution for exploration that avoids the fractious bickering between Mars and Moon advocates. This viewpoint proposes utilizing the concept of Greater Earth to educate policy makers, opinion makers and the public about these subtle attributes of our space neighborhood.     

采用Gauss伪谱法的小推力日-火Halo轨道转移优化设计

针对日-地Halo轨道到日-火Halo轨道的小推力轨道转移问题,给出一种基于不变流形理论和Gauss伪谱法的优化设计方法。首先,在日心惯性坐标系中建立小推力轨道优化模型,并基于不变流形理论给出轨道转移中流形出口和入口的选择原则,应用该原则在日-地系统中选择流形出口,在日-火系统中选择流形入口,并将其作为轨道转移的初末状态;然后基于Gauss伪谱法将最优控制问题离散化为非线性规划(NLP)问题,并采用基于逆多项式的形状算法给出了NLP初值的计算方法;最后对该轨道转移问题进行了数学仿真。仿真结果表明：Gauss伪谱法可有效用于小推力日-火Halo轨道转移的优化,且采用逆多项式形状算法得到的初值具有初始误差小,使得NLP收敛速度快的特点。     

Autonomous optical navigation for orbits around Earth–Moon collinear libration points

The analysis of optical navigation in an Earth–Moon libration point orbit is examined. Missions to libration points have been winning momentum during the last decades. Its unique characteristics make it suitable for a number of operational and scientific goals. Literature aimed to study dynamics, guidance and control of unstable orbits around collinear libration points is vast. In particular, several papers deal with the optimisation of the Δv budget associated to the station-keeping of these orbits. One of the results obtained in literature establishes the critical character of the Moon–Earth system in this aspect. The reason for this behaviour is twofold: high Δv cost and short optimal manoeuvre spacing. Optical autonomous navigation can address the issue of allowing a more flexible manoeuvre design. This technology has been selected to overcome similar difficulties in other critical scenarios. This paper analyses in detail this solution. A whole GNC system is defined to meet the requirements imposed by the unstable dynamic environment. Finally, a real simulation of a spacecraft following a halo orbit of the L2 Moon–Earth system is carried out to assess the actual capabilities of the optical navigation in this scenario.     

Spatial approaches to moons from resonance relative to invariant manifolds

In this study, the final approach to a moon or other body from resonance is explored and compared to the invariant manifolds of unstable periodic orbits. It is shown that the stable manifolds of planar Lyapunov orbits can act as a guide for the periods or resonances that are required for the final approach in both the planar and spatial problems. Previously developed techniques for the planar problem are expanded for use with resonances and used for comparison with trajectories approaching a moon from these resonances. A similar technique is then used for exploring the relationship of invariant manifolds to approach trajectories in the spatial problem. It is shown that the invariant manifolds of unstable periodic orbits provide insight into the trajectory design, and they can be used as a guide to the more direct approach trajectories.     

Statistical method for ground and sky coverage calculation

A novel statistical method has been devised for evaluating the ground and the sky coverage of an observation experiment on board a satellite. Owing to its unrivalled rapidity compared with other conventional calculation techniques, the method can be applied to evaluate the coverage percentages for the whole globe or any area on it, to calculate the visibility percentages for one or more ground stations and to determine the percentages of observation time of any given celestial direction including Sun, Moon, Earth and Ground Stations constraints. The orbits considered can be elliptical and account is taken of the drift due to the Earth's oblateness.     

Space technology transfer: Spin-off cases from Japan

Government organizations have to justify high expenditure during periods of financial crisis such as the one we are experiencing today. Space agencies have attempted to increase the returns on their investments in space missions by encouraging the commercial use of advanced technologies. This paper describes two technology transfer (TT) cases promoted by JAXA, in order to identify the organizational models and determinants of TT. The development of a TT process from space to Earth not only benefits the aerospace industry but also the network of national companies. The aim of the paper is to investigate who the actors are and the nature of their role, as well as the determinants of the TT process in the Japanese space sector. The case studies confirm the typical path of transfer as ‘Earth–space–Earth’.     

Planetary Defense From Space: Part 2 (Simple) Asteroid Deflection Law

A system of two space bases housing missiles for an efficient Planetary Defense of the Earth from asteroids and comets was firstly proposed by this author in 2002. It was then shown that the five Lagrangian points of the Earth–Moon system lead naturally to only two unmistakable locations of these two space bases within the sphere of influence of the Earth. These locations are the two Lagrangian points L1 (in between the Earth and the Moon) and L3 (in the direction opposite to the Moon from the Earth). In fact, placing missiles based at L1 and L3 would enable the missiles to deflect the trajectory of incoming asteroids by hitting them orthogonally to their impact trajectory toward the Earth, thus maximizing the deflection at best. It was also shown that confocal conics are the only class of missile trajectories fulfilling this “best orthogonal deflection” requirement.The mathematical theory developed by the author in the years 2002–2004 was just the beginning of a more expanded research program about the Planetary Defense. In fact, while those papers developed the formal Keplerian theory of the Optimal Planetary Defense achievable from the Earth–Moon Lagrangian points L1 and L3, this paper is devoted to the proof of a simple “(small) asteroid deflection law” relating directly the following variables to each other:

(1) the speed of the arriving asteroid with respect to the Earth (known from the astrometric observations);

(2) the asteroid's size and density (also supposed to be known from astronomical observations of various types);

(3) the “security radius” of the Earth, that is, the minimal sphere around the Earth outside which we must force the asteroid to fly if we want to be safe on Earth. Typically, we assume the security radius to equal about 10,000 km from the Earth center, but this number might be changed by more refined analyses, especially in the case of “rubble pile” asteroids;

(4) the distance from the Earth of the two Lagrangian points L1 and L3 where the defense missiles are to be housed;

(5) the deflecting missile's data, namely its mass and especially its “extra-boost”, that is, the extra-energy by which the missile must hit the asteroid to achieve the requested minimal deflection outside the security radius around the Earth.

This discovery of the simple “asteroid deflection law” presented in this paper was possible because:

(1) In the vicinity of the Earth, the hyperbola of the arriving asteroid is nearly the same as its own asymptote, namely, the asteroid's hyperbola is very much like a straight line. We call this approximation the line/circle approximation. Although “rough” compared to the ordinary Keplerian theory, this approximation simplifies the mathematical problem to such an extent that two simple, final equations can be derived.

(2) The confocal missile trajectory, orthogonal to this straight line, ceases then to be an ellipse to become just a circle centered at the Earth. This fact also simplifies things greatly. Our results are thus to be regarded as a good engineering approximation, valid for a preliminary astronautical design of the missiles and bases at L1 and L3.

Still, many more sophisticated refinements would be needed for a complete Planetary Defense System:

(1) taking into account many perturbation forces of all kinds acting on both the asteroids and missiles shot from L1 and L3;

(2) adding more (non-optimal) trajectories of missiles shot from either the Lagrangian points L4 and L5 of the Earth–Moon system or from the surface of the Moon itself;

(3) encompassing the full range of missiles currently available to the USA (and possibly other countries) so as to really see “which missiles could divert which asteroids”, even just within the very simplified scheme proposed in this paper.

In summary: outlined for the first time in February 2002, our Confocal Planetary Defense concept is a simplified Keplerian Theory that already proved simple enough to catch the attention of scholars, popular writers, and representatives of the US Military. These developments would hopefully mark the beginning of a general mathematical vision for building an efficient Planetary Defense System in space and in the vicinity of the Earth, although not on the surface of the Earth itself!We must make a real progress beyond academic papers, Hollywood movies and secret military plans, before asteroids like 99942 Apophis get close enough to destroy us in 2029 or a little later.