Time distributions in satellite constellation design

The aim of the time distribution methodology presented in this paper is to generate constellations whose satellites share a set of relative trajectories in a given time, and maintain that property over time without orbit corrections. The model takes into account a series of orbital perturbations such as the gravitational potential of the Earth, the atmospheric drag, the Sun and the Moon as disturbing third bodies and the solar radiation pressure. These perturbations are included in the design process of the constellation. Moreover, the whole methodology allows to design constellations with multiple relative trajectories that can be distributed in a minimum number of inertial orbits.     

基于iHCO通信卫星的CAPS星座优化研究

中国区域定位系统(Chinese Area Positioning System,CAPS)把寿命末期的地球静止轨道(Geostationary Earth Orbit,GEO)通信卫星推到比GEO轨道高约200 km的倾斜高圆轨道(inclined Highly Circular Orbit,iHCO),卫星相对地球向西漂移。利用该类卫星组建CAPS导航星座,可以实现全球范围内的导航通信覆盖。重点开展基于iHCO通信卫星的CAPS星座优化研究,结果表明:利用GEO通信卫星和iHCO通信卫星组成的星座可以实现较好的空间星座布局,可以满足一般导航用户的需要。     

CEMP星HE 1005-1439中子俘获元素天体物理来源的研究

HE1005-1439是一颗金属丰度极低([Fe/H]～-3.0)的碳增丰贫金属星(Carbon Enhanced Metal-Poor,CEMP),该星的s-过程元素显著超丰([Ba/Fe]=1.16±0.31,[Pb/Fe]=1.98±0.19),而r-过程元素温和超丰([Eu/Fe]=0.46±0.22),使用单一的s-过程模型和i-过程模型均不能拟合该星中子俘获丰度分布.采用丰度分解的方法探究该星化学元素的天体物理来源可有助于理解CEMP星的形成和化学演化.利用s-过程和r-过程的混合模型对其中子俘获元素的丰度分布进行拟合,发现该星的中子俘获元素主要来源于低质量低金属丰度AGB伴星的s-过程核合成,而r-过程核合成也有贡献.     

The 3-D lattice theory of Flower Constellations

Flower Constellations (FCs) have been extensively studied for use in optimal constellation design. The Harmonic FCs (HFCs) subset, representing the symmetric configurations, have recently been reformulated into 2-D Lattice Flower Constellations (2D-LFCs), encompassing the complete set of HFCs. Elliptic orbits are generally avoided due to the deleterious effects of Earth’s oblateness on the constellation, but here we present a novel concept for avoiding this problem and enabling more effective global coverage utilizing elliptic orbits. This new 3D Lattice Flower Constellations (3D-LFCs) framework generalizes the 2D-LFCs, Walker constellations, elliptical Walker constellations, and many of Draim’s global coverage constellations. Previous studies have shown FCs can provide improved performance in global navigation over existing Global Navigation Satellite Systems (GNSS). We found a 3D-LFC design that improved the average positioning accuracy by 3.5 % while reducing launch $\varDelta v$ Δ v requirements when compared to the existing Galileo GNSS constellation.     

The Characteristics and Related Problems of the Orbits Around the Earth-Moon Libration Points

Due to the specific dynamics, the probes located at the halo orbits or Lissajous orbits around the Earth-Moon collinear libration point L₁ or L₂ are always studied in the synodic system to understand their trajectories. In fact, they are also orbiting the Earth in a distant Keplerian ellipse. Because of their intrinsic orbital instability, in the orbit prediction the initial errors propagate more prominently than those of the normal orbiting satellites, this requires special attention in the orbit design, maneuver, and control. Despite of all this, they are similar to the normal orbiting satellites in orbit determination and hardly require other special attentions. In this paper, the quantitative results of error propagation under the unstable dynamics, together with the theoretical analysis are presented. The results of precise orbit determination and short-arc orbit predictions are also shown, and compared with the results from the Beijing Aerospace Control Center.     

A motivating exploration on lunar craters and low-energy dynamics in the Earth–Moon system

It is known that most of the craters on the surface of the Moon were created by the collision of minor bodies of the Solar System. Main Belt Asteroids, which can approach the terrestrial planets as a consequence of different types of resonance, are actually the main responsible for this phenomenon. Our aim is to investigate the impact distributions on the lunar surface that low-energy dynamics can provide. As a first approximation, we exploit the hyberbolic invariant manifolds associated with the central invariant manifold around the equilibrium point L ₂ of the Earth–Moon system within the framework of the Circular Restricted Three-Body Problem. Taking transit trajectories at several energy levels, we look for orbits intersecting the surface of the Moon and we attempt to define a relationship between longitude and latitude of arrival and lunar craters density. Then, we add the gravitational effect of the Sun by considering the Bicircular Restricted Four-Body Problem. In the former case, as main outcome, we observe a more relevant bombardment at the apex of the lunar surface, and a percentage of impact which is almost constant and whose value depends on the assumed Earth–Moon distance d_EM. In the latter, it seems that the Earth–Moon and Earth–Moon–Sun relative distances and the initial phase of the Sun θ ₀ play a crucial role on the impact distribution. The leading side focusing becomes more and more evident as d_EM decreases and there seems to exist values of θ ₀ more favorable to produce impacts with the Moon. Moreover, the presence of the Sun makes some trajectories to collide with the Earth. The corresponding quantity floats between 1 and 5 percent. As further exploration, we assume an uniform density of impact on the lunar surface, looking for the regions in the Earth–Moon neighbourhood these colliding trajectories have to come from. It turns out that low-energy ejecta originated from high-energy impacts are also responsible of the phenomenon we are considering.     

基于不同轨道高度iHCO通信卫星的CAPS星座优化研究

基于通信卫星的导航系统可以利用比地球静止轨道（Geostationary Earth Orbit,GEO）高约200 km的倾斜高圆轨道（inclined Highly Circular Orbit,iHCO）通信卫星组成导航星座.结合两种轨道高度的倾斜高圆轨道通信卫星,仿真分析了利用倾斜高圆轨道卫星组成的中国区域定位系统（Chinese Area Positioning System,CAPS）的导航性能,并讨论了利用倾斜高圆轨道卫星组成的中国区域定位系统实现中国区域覆盖的最佳星座布局.     

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)     

Origin of the moon by tidal capture and some geophysical consequences

Of the many proposed modes of origin of the Moon, some violate physical laws; many are in conflict with observations; all are improbable. Perhaps the least improbable - based on recent tidal theory calculations and on the interpretation of lunar rock data - is capture of the Moon as it passed near the Earth in adirect (prograde) orbit, shortly after the formation of Moon and Earth, about 4.5 billion years ago. (Capture of the Moon from an initiallyretrograde orbit which had been proposed some years ago, leads to physically unacceptable consequences.) The effects of capture on the Earth would have been cataclysmic, leading to intensive heating of its interior, to volcanism, and to the immediate formation of an atmosphere and hydrosphere. Thus capture of a Moon may have given rise to the unique properties of the Earth (in the Solar System) and to the early evolution of life, about 3.5 billion years ago.Presented at the NATO Advanced Study Institute on Lunar Studies in Patras, Greece, September, 1971.     

Weak stability boundary trajectories for the deployment of lunar spacecraft constellations

Suitable lunar constellation coverage can be obtained by separating the satellites in inclinations and node angles. It is shown in the paper that a relevant saving of velocity variation ΔV can be achieved using weak stability boundary trajectories. The weakly stable dynamics of such transfers allows the separation of the satellites from the nominal orbit to the required orbit planes with a small amount of ΔV. This paper also shows that only one different set of orbital parameters at Moon can be reached with the same ΔV manoeuvre starting from a nominal trajectory and ending at a fixed periselenium altitude. In fact, such a feature is proved to be common to other simpler dynamical systems, such as the two- and three-body problems.     

Recent results on the mass,gravitational field and moments of inertia of the moon

Doppler tracking data from the Lunar Orbiter series of spacecraft have been used in a more complete analysis of the spherical harmonic coefficients of the lunar gravitational field through thirteenth degree and order. The value obtained for the mass of the Moon,GM = 4902.84 km³ s^–2, is in good agreement with previous results and with results obtained by alternate procedures. Acceleration contour plots, derived from the gravitational coefficients, show correlations with surface features on the near side of the Moon, but are of questionable validity for the far side because of the lack of direct tracking data on the far side. Based on the most recent gravitational field data, the current estimate for the polar moment of inertia of the Moon isC/Ma ² = 0.4019_-0.002 ^+0.004. This value indicates that the interior of the Moon can be homogeneous, but some results presented strongly suggest that the Moon is differentiated, with an excess of mass in the direction toward the Earth.Paper presented at the NATO Advanced Study Institute on Lunar Studies, Patras, Greece, September, 1971.     

Regional positioning using a low Earth orbit satellite constellation

Global and regional satellite navigation systems are constellations orbiting the Earth and transmitting radio signals for determining position and velocity of users around the globe. The state-of-the-art navigation satellite systems are located in medium Earth orbits and geosynchronous Earth orbits and are characterized by high launching, building and maintenance costs. For applications that require only regional coverage, the continuous and global coverage that existing systems provide may be unnecessary. Thus, a nano-satellites-based regional navigation satellite system in Low Earth Orbit (LEO), with significantly reduced launching, building and maintenance costs, can be considered. Thus, this paper is aimed at developing a LEO constellation optimization and design method, using genetic algorithms and gradient-based optimization. The preliminary results of this study include 268 LEO constellations, aimed at regional navigation in an approximately 1000 km \(\times \) 1000 km area centered at the geographic coordinates [30, 30] degrees. The constellations performance is examined using simulations, and the figures of merit include total coverage time, revisit time, and geometric dilution of precision (GDOP) percentiles. The GDOP is a quantity that determines the positioning solution accuracy and solely depends on the spatial geometry of the satellites. Whereas the optimization method takes into account only the Earth’s second zonal harmonic coefficient, the simulations include the Earth’s gravitational field with zonal and tesseral harmonics up to degree 10 and order 10, Solar radiation pressure, drag, and the lunisolar gravitational perturbation.     

On quasi-periodic motions around the collinear libration points in the real Earth–Moon system

Due to various perturbations, the collinear libration points of the real Earth–Moon system are not equilibrium points anymore. Under the assumption that the Moon’s motion is quasi-periodic, special quasi-periodic orbits called dynamical substitutes exist. These dynamical substitutes replace the geometrical collinear libration points as time-varying equilibrium points. In the paper, the dynamical substitutes of the three collinear libration points in the real Earth–Moon system are computed. For the points L ₁ and L ₂, linearized motions around the dynamical substitutes are described, and the variational equations of the dynamical substitutes are reduced to a form with a near constant coefficient matrix. Then higher order analytical formulae of the central manifolds are constructed. Using these analytical solutions as initial seeds, Lissajous orbits and halo orbits are computed with numerical algorithms.     

A note on the dynamics around the Lagrange collinear points of the Earth–Moon system in a complete Solar System model

In this paper we study the dynamics of a massless particle around the L _1,2 libration points of the Earth–Moon system in a full Solar System gravitational model. The study is based on the analysis of the quasi-periodic solutions around the two collinear equilibrium points. For the analysis and computation of the quasi-periodic orbits, a new iterative algorithm is introduced which is a combination of a multiple shooting method with a refined Fourier analysis of the orbits computed with the multiple shooting. Using as initial seeds for the algorithm the libration point orbits of Circular Restricted Three Body Problem, determined by Lindstedt-Poincaré methods, the procedure is able to refine them in the Solar System force-field model for large time-spans, that cover most of the relevant Sun–Earth–Moon periods.     

On the construction of low-energy cislunar and translunar transfers based on the libration points

There exist cislunar and translunar libration points near the Moon, which are referred to as the LL ₁ and LL ₂ points, respectively. They can generate the different types of low-energy trajectories transferring from Earth to Moon. The time-dependent analytic model including the gravitational forces from the Sun, Earth, and Moon is employed to investigate the energy-minimal and practical transfer trajectories. However, different from the circular restricted three-body problem, the equivalent gravitational equilibria are defined according to the geometry of the instantaneous Hill boundary due to the gravitational perturbation from the Sun. The relationship between the altitudes of periapsis and eccentricities is achieved from the Poincaré mapping for all the captured lunar trajectories, which presents the statistical feature of the fuel cost and captured orbital elements rather than generating a specified Moon-captured segment. The minimum energy required by the captured trajectory on a lunar circular orbit is deduced in the spatial bi-circular model. The idea is presented that the asymptotical behaviors of invariant manifolds approaching to/traveling from the libration points or halo orbits are destroyed by the solar perturbation. In fact, the energy-minimal cislunar transfer trajectory is acquired by transiting the LL ₁ point, while the energy-minimal translunar transfer trajectory is obtained by transiting the LL ₂ point. Finally, the transfer opportunities for the practical trajectories that have escaped from the Earth and have been captured by the Moon are yielded by the transiting halo orbits near the LL ₁ and LL ₂ points, which can be used to generate the whole of the trajectories.     

The fate of primordial lunar Trojans

Multiple large impact basins on the lunar nearside formed in a relatively-short interval around 3.8-3.9 Gyr ago, in what is known as the Lunar Cataclysm (LC; also known as Late Heavy Bombardment). It is widely thought that this impact bombardment has affected the whole Solar System or at least all the inner planets. But with non-lunar evidence for the cataclysm being relatively weak, a geocentric cause of the Lunar Cataclysm cannot yet be completely ruled out [Ryder, G., 1990. Eos 71, 313, 322-323]. In principle, late destabilization of an additional Earth satellite could result in its tidal disruption during a close lunar encounter (cf. [Asphaug, E., Agnor, C.B., Williams, Q., 2006. Nature 439, 155-160]). If the lost satellite had D>500 km, the resulting debris can form multiple impact basins in a relatively short time, possibly explaining the LC. Canup et al. [Canup, R.M., Levison, H.F., Stewart, G.R., 1999. Astron. J. 117, 603-620] have shown that any additional satellites of Earth formed together with (and external to) the Moon would be unable to survive the rapid initial tidally-driven expansion of lunar orbit. Here we explore the fate of objects trapped in the lunar Trojan points, and find that small lunar Trojans can survive the Moon's orbital evolution until they and the Moon reach 38 Earth radii, at which point they are destabilized by a strong solar resonance. However, the dynamics of Trojans containing enough mass to cause the LC (diameters >150 km) is more complex; we find that such objects do not survive the passage through a weaker solar resonance at 27 Earth radii. This distance was very likely reached by the Moon long before the LC, which seems to rule out the disruption of lunar Trojans as a cause of the LC.     

Effect of the Earth’s compression on the physical libration of the Moon

The effect of the Earth??s compression on the physical libration of the Moon is studied using a new vector method. The moment of gravitational forces exerted on the Moon by the oblate Earth is derived considering second order harmonics. The terms in the expression for this moment are arranged according to their order of magnitude. The contribution due to a spherically symmetric Earth proves to be greater by a factor of 1.34 × 10⁶ than a typical term allowing for the oblateness. A linearized Euler system of equations to describe the Moon??s rotation with allowance for external gravitational forces is given. A full solution of the differential equation describing the Moon??s libration in longitude is derived. This solution includes both arbitrary and forced oscillation harmonics that we studied earlier (perturbations due to a spherically symmetric Earth and the Sun) and new harmonics due to the Earth??s compression. We posed and solved the problem of spinorbital motion considering the orientation of the Earth??s rotation axis with regard to the axes of inertia of the Moon when it is at a random point in its orbit. The rotation axes of the Earth and the Moon are shown to become coplanar with each other when the orbiting Moon has an ecliptic longitude of L _? = 90° or L _? = 270°. The famous Cassini??s laws describing the motion of the Moon are supplemented by the rule for coplanarity when proper rotations in the Earth-Moon system are taken into account. When we consider the effect of the Earth??s compression on the Moon??s libration in longitude, a harmonic with an amplitude of 0.03?? and period of T ₈ = 9.300 Julian years appears. This amplitude exceeds the most noticeable harmonic due to the Sun by a factor of nearly 2.7. The effect of the Earth??s compression on the variation in spin angular velocity of the Moon proves to be negligible.     

Strategies for plane change of Earth orbits using lunar gravity and derived trajectories of family G

The dynamics of the circular restricted three-body Earth-Moon-particle problem predicts the existence of the retrograde periodic orbits around the Lagrangian equilibrium point L1. Such orbits belong to the so-called family G (Broucke, Periodic orbits in the restricted three-body problem with Earth-Moon masses, JPL Technical Report 32–1168, 1968) and starting from them it is possible to define a set of trajectories that form round trip links between the Earth and the Moon. These links occur even with more complex dynamical systems as the complete Sun-Earth-Moon-particle problem. One of the most remarkable properties of these trajectories, observed for the four-body problem, is a meaningful inclination gain when they penetrate into the lunar sphere of influence and accomplish a swing-by with the Moon. This way, when one of these trajectories returns to the proximities of the Earth, it will be in a different orbital plane from its initial Earth orbit. In this work, we present studies that show the possibility of using this property mainly to accomplish transfer maneuvers between two Earth orbits with different altitudes and inclinations, with low cost, taking into account the dynamics of the four-body problem and of the swing-by as well. The results show that it is possible to design a set of nominal transfer trajectories that require ΔV _Total less than conventional methods like Hohmann, bi-elliptic and bi-parabolic transfer with plane change.     

Minimum fuel control of the planar circular restricted three-body problem

The circular restricted three-body problem is considered to model the dynamics of an artificial body submitted to the attraction of two planets. Minimization of the fuel consumption of the spacecraft during the transfer, e.g. from the Earth to the Moon, is considered. In the light of the controllability results of Caillau and Daoud (SIAM J Control Optim, 2012), existence for this optimal control problem is discussed under simplifying assumptions. Thanks to Pontryagin maximum principle, the properties of fuel minimizing controls is detailed, revealing a bang-bang structure which is typical of L¹-minimization problems. Because of the resulting non-smoothness of the Hamiltonian two-point boundary value problem, it is difficult to use shooting methods to compute numerical solutions (even with multiple shooting, as many switchings on the control occur when low thrusts are considered). To overcome these difficulties, two homotopies are introduced: One connects the investigated problem to the minimization of the L²-norm of the control, while the other introduces an interior penalization in the form of a logarithmic barrier. The combination of shooting with these continuation procedures allows to compute fuel optimal transfers for medium or low thrusts in the Earth–Moon system from a geostationary orbit, either towards the L ₁ Lagrange point or towards a circular orbit around the Moon. To ensure local optimality of the computed trajectories, second order conditions are evaluated using conjugate point tests.     

Utilization of GPS satellites for precise irradiation measurement and monitoring

Precise measurement of irradiance over the earth under various circumstances like solar flares, coronal mass ejections, over an 11-year solar cycle, etc. leads to better understanding of Sun-earth relationship. To continuously monitor the irradiance over earth-space regions several satellites at several positions are required. For that continuous and multiple satellite monitoring we can use GPS (Global Positioning System) satellites (like GLONASS, GALILEO, future satellites) installed with irradiance measuring and monitoring instruments. GPS satellite system consists of 24 constellations of satellites. Therefore usage of all the satellites leads to 24 measurements of irradiance at the top of the atmosphere (or 12 measurements of those satellites which are pointing towards the Sun) at an instant. Therefore in one day, numerous irradiance observations can be obtained for the whole globe, which will be very helpful for several applications like Albedo calculation, Earth Radiation Budget calculation, monitoring of near earth-space atmosphere, etc. Moreover, measuring irradiance both in ground (using ground instruments) and in space at the same instant of time over a same place, leads to numerous advantages. That is, for a single position we obtain irradiance at the top of the atmosphere, irradiance at ground and the difference in irradiance from over top of the atmosphere to the ground. Measurement of irradiance over the atmosphere and in ground at a precise location gives more fine details about the solar irradiance influence over the earth, path loss and interaction of irradiance with the atmosphere.