Tetrahedron formation of nanosatellites with single-input control

The present paper studies the formation flight of four nanosatellites forming a tetrahedron. The main goal of this research is to find the relative orbits of these satellites that, at least in the linear Hill–Clohessy–Wiltshire model, ensure finite relative motion and keep the volume and shape of the tetrahedron configuration. Since real motions of these satellites will differ from the linear ones, especially under the influence of the \(J_{2}\) perturbation, active control is necessary. In addition, the limited size of the satellites does not allow us to use a complex 3-axis attitude control system. In the present paper we consider the passive magnetic attitude control system and suppose that the thrust direction is always aligned with the local geomagnetic field. In order to increase mission lifetime the control algorithm that minimizes the propellant consumption and keeps the tetrahedron volume and shape is investigated.     

Bounded relative motions for formation flying in displaced orbits by low-thrust propulsion

This paper discusses a numerical searching approach for the relative motion of formation flying in displaced orbits by spacecraft with low-thrust propulsion. The nonlinear dynamical model of spacecraft is established in a two-body rotating reference frame with arbitrary polar component of momentum and thrust-induced acceleration. Motions near the stable equilibria are distinguished from each other by means of five-dimensional variables, which can then be compressed uniquely into two-dimensional mapping images characterized by the crossing interval and the angle drifts. The surjective but not injective mapping makes the generation of three configurations of the relative motions possible. The corresponding relative orbits for three kinds of two-spacecraft formation flying are searched and exemplified based on the formation conditions formulized as functions of the crossing interval and the angle drifts. Furthermore, based on the assignment of displaced relative orbits to five-dimensional vector, the working orbit of the deputy for a specific chief can also be searched via the optimization algorithm to generate the bounded relative motion with the minimum thrust acceleration magnitude, which is of certain significance in reducing fuel consumption of formations.     

The tidal stripping of satellites

We present an improved analytic calculation for the tidal radius of satellites and test our results against N -body simulations.
The tidal radius in general depends upon four factors: the potential of the host galaxy, the potential of the satellite, the orbit of the satellite and the orbit of the star within the satellite . We demonstrate that this last point is critical and suggest using three tidal radii to cover the range of orbits of stars within the satellite. In this way we show explicitly that prograde star orbits will be more easily stripped than radial orbits; while radial orbits are more easily stripped than retrograde ones. This result has previously been established by several authors numerically, but can now be understood analytically. For point mass, power-law (which includes the isothermal sphere), and a restricted class of split power-law potentials our solution is fully analytic. For more general potentials, we provide an equation which may be rapidly solved numerically.
Over short times (≲1–2 Gyr ∼1 satellite orbit), we find excellent agreement between our analytic and numerical models. Over longer times, star orbits within the satellite are transformed by the tidal field of the host galaxy. In a Hubble time, this causes a convergence of the three limiting tidal radii towards the prograde stripping radius. Beyond the prograde stripping radius, the velocity dispersion will be tangentially anisotropic.     

Origin of orbits of secondaries in the discovered trans-Neptunian binaries

The dependences of inclinations of orbits of secondaries in the discovered trans-Neptunian binaries on the distance between the primary and the secondary, on the eccentricity of orbits of the secondary around the primary, on the ratio of diameters of the secondary and the primary, and on the elements of heliocentric orbits of these binaries are studied. These dependences are interpreted using the model of formation of a satellite system in a collision of two rarefied condensations composed of dust and/or objects less than 1 m in diameter. It is assumed in this model that a satellite system forms in the process of compression of a condensation produced in such a collision. The model of formation of a satellite system in a collision of two condensations agrees with the results of observations: according to observational data, approximately 40% of trans-Neptunian binaries have a negative angular momentum relative to their centers of mass.     

Resonant orbits of grains and the formation of satellites

Grains, an abundant constituent of the former solar system, will have had a high probability of being driven into orbits resonant with major bodies already formed. This arises because of the presence of gas drag and Poynting-Robertson drag on small grains, providing the dissipation necessary to concentrate matter into special orbits. Since the mean density in resonant orbits can be built up by such a process without limit, these may become the favored orbits for gravitational contraction to gather material into major bodies. Satellite formation processes may therefore depend upon the buildup of resonant lanes of dust grains around the parent body. Saturn's rings are possibly one example of such lanes, though an unsuitable one for the final step of satellite formation on account of their being too close to Saturn.     

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.     

Exploiting Unstable Periodic Orbits of a Chaotic Invariant Set for Spacecraft Control

The purpose of this work is to show that chaos control techniques (OGY, in special) can be used to efficiently keep a spacecraft around another body performing elaborate orbits. We consider a satellite and a spacecraft moving initially in coplanar and circular orbits, with slightly different radii, around a heavy central planet. The spacecraft, which is the inner body, has a slightly larger angular velocity than the satellite so that, after some time, they eventually go to a situation in which the distance between them becomes sufficiently small, so that they start to interact with one another. This situation is called as an encounter. In previous work we have shown that this scenario is a typical situation of a chaotic scattering for some well-defined range of parameters. Considering this scenario, we first show how it is possible to find the unstable periodic orbits that are located in the chaotic invariant set. From the set of unstable periodic orbits, we select the ones that can be combined to provide the desired elaborate orbit. Then, chaos control technique based on the OGY method is used to keep the spacecraft in the desired orbit. Finally, we analyze the results and make considerations regarding a realistic scenario of space exploration.     

Orbital evolution of long-period comets I. Trivariate Monte Carlo simulation

The mean orbital evolution of long-period comets for 16 representative initial orbits to short-period comets is calculated by a Monte Carlo method. First, trivariate perturbation distributions of barycentric Kepler energy, total angular momentum, and its z component in single encounters of comets with Jupiter are obtained numerically. Their characteristics are examined in detail and the distributions are found to be simple, symmetric, and easy to handle. Second, utilizing these distributions, we have done trivariate Monte Carlo simulations of the orbital evolution of long-period comets, with special emphasis on high-inclination orbits. About half of the 16 initial orbits are traced up to 5000 returns. For each of these orbits, the mean values of semimajor axis, perihelion distance, and inclination; their standard deviations, survival, and capture rates; as well as time scales of orbital evolution are calculated as functions of return number. Survival rates of the initial orbits with high inclination (～90°) and small perihelion distance (～1–2 AU) have been found to be only two or three times smaller than those of the main-source orbits of short-period comets established quantitatively by Everhart. The time scales of orbitsl evolution of the former, however, are nearly 10 times longer than the latter. There is a general trend that, for smaller perihelion distance, the survival efficiency becomes higher. The results of this paper should be considered a basis for a succeeding paper (Paper II) in which the physical lifetime of comets will be determined, and a comparison with the orbital data will be done.     

定点在日-地(月)系L1点附近的探测器的发射及维持

在限制性三体问题中共线平动点附近的运动虽然是不稳定的,但可以是有条件稳定的,该动力学特征使得一些有特殊目的的探测器只需消耗较少的能量即可定点在这些点附近(如ISEE-3、SOHO).以日-地(月)系的L1点为例,根据其附近的运动特征,探讨定点探测器的发射与轨道控制问题,给出了相应的数值模拟结果,为工程上的实现提供理论依据.     

The future of artificial satellite theories Hybrid ephemeris compression model

Since the time of Newton, astrodynamics has focused on the analytical solution of orbital problems. This was driven by the desire to obtain a theoretical understanding of the motion and the practical desire to be able to produce a computational result, Only with the advent of the computer did numerical integration become a practical consideration for solving dynamical problems. Although computer technology is not yet to the point of being able to provide numerical integration support for all satellite orbits, we are in a transition period which is being driven by the unprecedented increase in computational power, This transition will affect the future of analytical, semi-analytical and numerical artificial satellite theories in a dramatic way, In fact, the role for semi-analytical theories may disappear. During the time of transition, a central site may have the capacity to maintain the orbits using numerical integration, but the user may not have such a capacity or may need results in a more timely manner, One way to provide for this transition need is through the use of some type of satellite ephemeris compression. Through the combined use of a power series and a Fourier series, good quality ephemeris compression has been achieved for 7 day periods, The ephemeris compression requires less than 40 terms and is valid for all eccentricities and inclinations.     

Celestial-mechanical peculiarities of Uranus’s satellite system

We consider the structural peculiarities of Uranus’s satellite system associated with its separation into two groups: inner equatorial satellites moving in nearly circular orbits and distant irregular satellites with retrograde motion in highly elliptical orbits. The intermediate region is free from satellites in a wide range of semimajor axes. By analyzing the evolution of satellite orbits under the combined effect of solar attraction and Uranus’s oblateness, we offer a celestial-mechanical explanation for the absence of equatorial satellites in this region. M.L. Lidov’s studies during 1961–1963 have served as a basis for our analysis.     

月球卫星轨道力学综述

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

Possible satellites of Mercury and Venus

On the basis of works of King and Innanen, the limiting direct and retrograde orbits around the planets Mercury and Venus have been calculated. Synthesizing this concept with the concept of synchronous orbits around the planets and tidal drags acting within them it is shown that Venus may not have retained any satellite direct or retrograde but Mercury may have retained a retrograde satellite at a distance between 225000 and 252700 km from its center. It is urged that this satellite may be investigated observationally.     

Numerical investigation of the orbit of Interball‐1

We investigated the motion of the Earth's artificial satellite Interball‐1 by using a method suitable for the computation of large eccentricity orbits. Though the measured and the computed orbital elements differ from each other within the measured error bound, we found a slight tendency for secular decreasing in the semi‐major axis, caused probably by electromagnetic drag. We analysed the dominant role of the Moon in the variations of the orbital eccentricity, leading to zero perigee height and the end of the lifetime of the satellite. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Satellite Relative Motion Using Differential Equinoctial Elements

For precise control, to minimize the fuel consumption, and to maximize the lifetime of satellite formations a precise analytic solution is needed for the relative motion of satellites. Based on the relationship between the relative states and the differential orbital elements, the state transition matrix for the linearized relative motion that includes the effects due to the reference orbit eccentricity and the gravitational perturbations is derived. This method is called the Geometric Method. To avoid any singularities at zero eccentricity and zero inclination, equinoctial variables are used to derive the relative motion state transition matrices for both mean and osculating elements. This approach can be extended easily to include other perturbing forces.     

Constructive analytical solution of the evolution hill problem

A new analytical solution of the system of differential equations describing secular perturbations and long-period solar perturbations of mean orbits of outer satellites of giant planets was obtained. As distinct from other solutions, the solution constructed using von Zeipel’s method approximately takes into account, in the secular part of the perturbing function, the totality of fourth order with respect to the small parameter m of the ratio of the mean motions of the primary planet and the satellite. This enables us to describe more accurately the evolution of satellite orbits with large apocentric distances, which in the course of evolution may exceed the halved radius of the Hill sphere of the planet with respect to the Sun. Among these are the orbits of the two outermost Neptunian satellites N10 (Psamathe) and N13 (Neso). For these satellites, the parameter m amounts to 0.152 and 0.165, respectively. Different from a purely analytical solution, the proposed solution requires preliminary calculations for each satellite. More precisely, in doing so, we need to construct some simple functions to approximate more complex ones. This is why we use the phrase “constructive analytical.” To illustrate the solution, we compare it with the results of the numerical integration of the strict motion equations of the satellites N10 and N13 over time intervals 5–15 thousand years.     

Determination of the optimal conditions for inclination maneuvers using a Swing-by

The search for methods to reduce the fuel consumption in orbital transfers is something relevant and always current in astrodynamics. Therefore, the maneuvers assisted by the gravity, also called Swing-by maneuvers, can be an advantageous option to save fuel. The proposal of the present research is to explore the influence of some parameters in a Swing-by of an artificial satellite orbiting a planet with one of the moons of this mother planet, with the goal of changing the inclination of the artificial satellite around the main body of the system. The fuel consumption of this maneuver is compared with the required consumption to perform the same change of inclination using the classical approach of impulsive maneuvers.     

A numerical-analytical method for studying the orbital evolution of distant planetary satellites

We describe an approximate numerical-analytical method for calculating the perturbations of the elements of distant satellite orbits. The model for the motion of a distant satellite includes the solar attraction and the eccentricity and ecliptic inclination of the orbit of the central planet. In addition, we take into account the variations in planetary orbital elements with time due to secular perturbations. Our work is based on Zeipel’s method for constructing the canonical transformations that relate osculating satellite orbital elements to the mean ones. The corresponding transformation of the Hamiltonian is used to construct an evolution system of equations for mean elements. The numerical solution of this system free from rapidly oscillating functions and the inverse transformation from the mean to osculating elements allows the evolution of distant satellite orbits to be studied on long time scales on the order of several hundred or thousand satellite orbital periods.