A general timing condition for consecutive collision orbits in the limiting case µ = 0 of the elliptic restricted problem

Consecutive collision orbits in the limiting case µ = 0 of the elliptic restricted three-body problem are investigated. in particular those in which the infinitesimal mass collides twice with the smaller (massless) primary. A timing condition is presented that allows the extension of previous results to the case of arbitrary relative orientation of the orbits of the infinitesimal mass and the smaller primary. The timing condition is expressed in two general forms - in terms of orbit parameters and eccentric (or hyperbolic) anomalies at the times of collision - for the specific cases of elliptic. parabolic or hyperbolic orbits of the infinitesimal mass. Some families of solutions are presented.     

The theory of asynchronous relative motion II: universal and regular solutions

Two fully regular and universal solutions to the problem of spacecraft relative motion are derived from the Sperling–Burdet (SB) and the Kustaanheimo–Stiefel (KS) regularizations. There are no singularities in the resulting solutions, and their form is not affected by the type of reference orbit (circular, elliptic, parabolic, or hyperbolic). In addition, the solutions to the problem are given in compact tensorial expressions and directly referred to the initial state vector of the leader spacecraft. The SB and KS formulations introduce a fictitious time by means of the Sundman transformation. Because of using an alternative independent variable, the solutions are built based on the theory of asynchronous relative motion. This technique simplifies the required derivations. Closed-form expressions of the partial derivatives of orbital motion with respect to the initial state are provided explicitly. Numerical experiments show that the performance of a given representation of the dynamics depends strongly on the time transformation, whereas it is virtually independent from the choice of variables to parameterize orbital motion. In the circular and elliptic cases, the linear solutions coincide exactly with the results obtained with the Clohessy–Wiltshire and Yamanaka–Ankersen state-transition matrices. Examples of relative orbits about parabolic and hyperbolic reference orbits are also presented. Finally, the theory of asynchronous relative motion provides a simple mechanism to introduce nonlinearities in the solution, improving its accuracy.     

Modified Newton-Raphson methods for preliminary orbit determination

Preliminary orbit determination is a multipoint boundary value problem which may be solved by the generalized Newton-Raphson iteration. When applied formally the method suffers from extensive computer storage requirements, fairly long execution times and in some cases, insufficient accuracy. In this work we seek to remove these practical difficulties via modification of the computational algorithm in such a way that solution storage is eliminated for the most part and computational speed and tolerance to imprecise integration algorithms is improved. The modified methods are applied to nine typical preliminary orbit determination problems to demonstrate fast convergence and short computation times, even with very poor starting values for the iteration. Excellent precision of the resulting solution is also demonstrated as well as the algorithm's ability to handle circular, elliptic, parabolic and hyperbolic orbits.     

Visual binary star orbits in universal variables

In this paper, we make use of the Stumpff's functions to solve the problem of determining the orbit of a visual binary star in universal variables. The method is thus valid for all types of orbits: hyperbolic, parabolic and elliptic.     

月球卫星最优小推力变轨研究

对利用小推力发动机将月球探器从双曲线轨道转移到圆轨的燃料最省转移问题,进行了研究,首先,将问题分解为双曲线到椭圆的转移和从椭圆到目标圆轨道的转移两步,然后,分别利用遗传算法解决了冲量假设下的最估转移、小推力加速民政部下从双曲线到椭圆的转移轨道优化,以及转移时间有约束情况下的从椭圆到圆轨道的转移轨道优化问题。     

The stability of masses during three-body encounters

The dynamical interactio of a binary system and a third body not moving on a closed orbit arises in a large number of physical situations. The C²H condition for determining Hill stability of coplanar bound three-body systems is extended to cover situations where the outer body moves on a parabolic or hyperbolic orbit. Regions where such a body is stable against exchange or collision with other components of the system are determined for a number of important cases where closed solutions are possible.     

The transition from elliptic to hyperbolic orbits in the two-body problem by slow loss of mass

Transition from elliptic to hyperbolic orbits in the two-body problem with slowly decreasing mass is investigated by means of asymptotic approximations.Analytical results by Verhulst and Eckhaus are extended to construct approximate solutions for the true anomaly and the eccentricity of the osculating orbit if the initial conditions are nearly-parabolic. It becomes clear that the eccentricity will monotonously increase with time for all mass functions satisfying a Jeans-Eddington relation and even for a larger set of functions. To illustrate these results quantitatively we calculate the eccentricity as a function of time for Jeans-Eddington functionsn=0(1) 5 and 18 nearly-parabolic initial conditions to find that 93 out of 108 elliptic orbits become hyperbolic.     

Steady state obliquity of a rigid body in the spin–orbit resonant problem: application to Mercury

In this paper, we consider the elliptic collinear solutions of the classical n-body problem, where the n bodies always stay on a straight line, and each of them moves on its own elliptic orbit with the same eccentricity. Such a motion is called an elliptic Euler–Moulton collinear solution. Here we prove that the corresponding linearized Hamiltonian system at such an elliptic Euler–Moulton collinear solution of n-bodies splits into \((n-1)\) independent linear Hamiltonian systems, the first one is the linearized Hamiltonian system of the Kepler 2-body problem at Kepler elliptic orbit, and each of the other \((n-2)\) systems is the essential part of the linearized Hamiltonian system at an elliptic Euler collinear solution of a 3-body problem whose mass parameter is modified. Then the linear stability of such a solution in the n-body problem is reduced to those of the corresponding elliptic Euler collinear solutions of the 3-body problems, which for example then can be further understood using numerical results of Martínez et al. on 3-body Euler solutions in 2004–2006. As an example, we carry out the detailed derivation of the linear stability for an elliptic Euler–Moulton solution of the 4-body problem with two small masses in the middle.     

Bounds on the Solution to kEPLER'S EQUATION:II. UNIVERSAL AND OPTIMAL STARTING POINTS

In this paper we find bounds on the solution to Kepler's equation for hyperbolic and parabolic motions. Two general concepts introduced here may be proved useful in similar numerical problems. Moreover, we give optimal starting points for Kepler's equation in hyperbolic and elliptic motions with particular attention to nearly parabolic orbits. It allows to expand the accepted earlier interval |e - 1| ≤ 0.01 for nearly parabolic orbits to the interval |e - 1| ≤ 0.05. This revised version was published online in July 2006 with corrections to the Cover Date.     

Quasiintegrals of the photogravitational eccentric restricted three-body problem with Poynting–Robertson drag

A restricted three-body problem for a dust particle, in presence of a spherical cometary nucleus in an eccentric (elliptic, parabolic or hyperbolic) orbit about the Sun, is considered. The force of radiation pressure and the Poynting– Robertson effect are taken into account. The differential equations of the particle’s non-inertial spatial motion are investigated both analytically and numerically. With the help of a complex representation, a new single equation of the motion is obtained. Conversion of the equations of motion system into a single equation allows the derivation of simple expressions similar to the integral of energy and integrals of areas. The derived expressions are named quasiintegrals. Relative values of terms of the energy quasiintegral for a smallest, largest, and a mean comet are calculated. We have found that in a number of cases the quasiintegrals are related to the regular integrals of motion, and discuss how the quasiintegrals may be applied to find some significant constraints on the motion of a body of infinitesimal mass.     

Fast Procedure Solving Universal Kepler's Equation

We developed a procedure to solve a modification of the standard form of the universal Kepler’s equation, which is expressed as a nondimensional equation with respect to a nondimensional variable. After reducing the domain of the variable and the argument by using the symmetry and the periodicity of the equation, the method first separates the case where the solution is so small that it is given an inverted series. Second, it separates the cases where the elliptic, parabolic, or hyperbolic standard forms of Kepler’s equation are suitable. Here the separation is done by judging whether detouring these nonuniversal equations will cause a 1-bit loss of information to their nonuniversal solutions or not. Then the nonuniversal equations are solved by the author’s procedures to solve the elliptic Kepler’s equation (Fukushima, 1997a), Barker’s equation (Fukushima, 1998), and the hyperbolic Kepler’s equation (Fukushima, 1997b), respectively. And their nonuniversal solutions are transformed back to the solution of the universal equation. For the rest of the case, we obtain an approximate solution by solving roughly the approximated cubic equation as we did in solving Barker’s equation. Then the correction to the approximate solution is obtained by Halley’s method precisely. There the special function appeared in the universal equation is rewritten into a combination of similar special functions of small arguments, so that they are efficiently evaluated by their Taylor series. Numerical measurements showed that, in the case of Intel Pentium II processor, the new method is 10–25 times as fast as Shepperd’s method (Shepperd, 1985) and 7–13 times as fast as the standard Newton method. This revised version was published online in July 2006 with corrections to the Cover Date.     

Relative motion of satellites exploiting the super-integrability of Kepler’s problem

This paper builds upon the work of Palmer and Imre exploring the relative motion of satellites on neighbouring Keplerian orbits. We make use of a general geometrical setting from Hamiltonian systems theory to obtain analytical solutions of the variational Kepler equations in an Earth centred inertial coordinate frame in terms of the relevant conserved quantities: relative energy, relative angular momentum and the relative eccentricity vector. The paper extends the work on relative satellite motion by providing solutions about any elliptic, parabolic or hyperbolic reference trajectory, including the zero angular momentum case. The geometrical framework assists the design of complex formation flying trajectories. This is demonstrated by the construction of a tetrahedral formation, described through the relevant conserved quantities, for which the satellites are on highly eccentric orbits around the Sun to visit the Kuiper belt.     

Secular Motion in a 2nd Degree and Order-Gravity Field with no Rotation

The motion of a particle about a non-rotating 2nd degree and order-gravity field is investigated. Averaging conditions are applied to the particle motion and a qualitative analysis which reveals the general character of motion in this system is given. It is shown that the orbit plane will either be stationary or precess about the body's axis of minimum or maximum moment of inertia. It is also shown that the secular equations for this system can be integrated in terms of trigonometric, hyperbolic or elliptic functions. The explicit solutions are derived in all cases of interest.This revised version was published online in October 2005 with corrections to the Cover Date.     

Application of Hamiltonian structure-preserving control to formation flying on a J 2-perturbed mean circular orbit

The bounded quasi-periodic relative trajectories are investigated in this paper for on-orbit surveillance, inspection or repair, which requires rapid changes in formation configuration for full three-dimensional imaging and unpredictable evolutions of relative trajectories for non-allied spacecraft. A linearized differential equation for modeling J ₂ perturbed relative dynamics is derived without any simplified treatment of full short-period effects. The equation serves as a nominal reference model for stationkeeping controller to generate the quasi-periodic trajectories near the equilibrium, i.e., the location of the chief. The developed model exhibits good numerical accuracy and is applicable to an elliptic orbit with small eccentricity inheriting from the osculating conversion of orbital elements. A Hamiltonian structure-preserving controller is derived for the three-dimensional time-periodic system that models the J ₂-perturbed relative dynamics on a mean circular orbit. The equilibrium of the system has time-varying topological types and no fixed-dimensional unstable/stable/center manifolds, which are quite different from the two-dimensional time-independent system with a permanent pair of hyperbolic eigenvalues and fixed-dimensions of unstable/stable/ center manifolds. The unstable and stable manifolds are employed to change the hyperbolic equilibrium to elliptic one with the poles assigned on the imaginary axis. The detailed investigations are conducted on the critical controller gain for Floquet stability and the optimal gain for the fuel cost, respectively. Any initial relative position and velocity leads to a bounded trajectory around the controlled elliptic equilibrium. The numerical simulation indicates that the controller effectively stabilizes motions relative to the perturbed elliptic orbit with small eccentricity and unperturbed elliptic orbit with arbitrary eccentricity. The developed controller stabilizes the quasi-periodic relative trajectories involved in six foundational motions with different frequencies generated by the eigenvectors of the Floquet multipliers, rather than to track a reference relative configuration. Only the relative positions are employed for the feedback without the information from the direct measurement or the filter estimation of relative velocity. So the current controller has potential applications in formation flying for its less computation overload for on-board computer, less constraint on the measurements, and easily-achievable quasi-periodic relative trajectories.     

Analytical solutions for Euler parameters

Two new analytical solutions for Poinsot motion in terms of Euler parameters are derived. The first solution is a straightforward ‘universal’ (no branches) time series practical for short time motion calculations or as a basis for analytical continuation. The second, more involved solution is also universal but is not restricted to short times; it is in terms of circular, hyperbolic, and elliptic functions and elliptic integrals.     

An exact analytical solution of Kepler's equation

Complex-variable analysis is used to develop an exact solution to Kepler's equation, for both elliptic and hyperbolic orbits. The method is based on basic properties of canonical solutions to appropriately posed Riemann problems, and the final results are expressed in terms of elementary quadratures.     

N-Body Simulation of a Deeply Penetrating Collision Between Unequal Galaxies

We study the tidal effects of a deeply penetrating collision between two spherical galaxies, one twice massive but less dense than the other, by numerical simulations. We consider the relative motion of the galaxies to be initially in a hyperbolic orbit. The collision parameters are so chosen that the primary (bigger) galaxy is just below the limit of disruption and the relative velocity of the pair is slightly in excess of the escape limit and the primary suffer greater tidal damage than the secondary. The primary develops a core halo structure and shows over all expansion while the secondary while the secondary shows contraction in the inner region and less significant expansion in the outer parts. The initially hyperbolic orbit is transformed into a parabolic orbit as a result of the collision. The result also indicate that the tidal interaction does not induce appreciable rotation in hyperbolic collision. We calculate the angle of deflection of the orbit and compare it with that computed using analytical work. The numerical work shows larger angle of deflection which is attributed to the large tidal effects of the bigger galaxy in the interpenetrating collision. This revised version was published online in July 2006 with corrections to the Cover Date.     

Out-of-plane librations of spinning tethered satellite systems

We analyze the out-of-plane librations of a tethered satellite system that is nominally rotating in the orbit plane. To isolate the librational dynamics, the system is modeled as two point masses connected by a rigid rod with the system mass center constrained to an unperturbed circular orbit. For small out-of-plane librations, the in-plane motion is unaffected by the out-of-plane librations and a solution for the in-plane motion is determined in terms of Jacobi elliptic functions. This solution is used in the linearized equation for the out-of-plane librations, resulting in a Hill’s equation. Floquet theory is used to analyze the Hill’s equation, and we show that the out-of-plane librations are unstable for certain ranges of in-plane spin rate. For relatively high in-plane spin rates, the out-of-plane librations are stable, and the Hill’s equation can be approximated by a Mathieu’s equation. Approximate solutions to the Mathieu’s equation are determined, and we analyze the dominant characteristics of the out-of-plane librations for high in-plane spin rates. The results obtained from the analysis of the linearized equations of motion are compared to numerical simulations of the nonlinear equations of motion, as well as numerical simulations of a more realistic system model that accounts for tether flexibility. The instabilities discovered from the linear analysis are present in both the nonlinear system and the more realistic system model. The approximate solutions for the out-of-plane librations compare well to the nonlinear system for relatively high in-plane rotation rates, and also capture the significant qualitative behavior of the flexible system.     

Chaotic scattering in the restricted three-body problem I. The Copenhagen problem

We consider the scattering motion of the planar restricted three-body problem with two equal masses on a circular orbit. Using the methods of chaotic scattering we present results on the structure of scattering functions. Their connection with primitive periodic orbits and the underlying chaotic saddle are studied. Numerical evidence is presented which suggests that in some intervals of the Jacobi integral the system is hyperbolic. The Smale horseshoe found there is built from a countable infinite number of primitive periodic orbits, where the parabolic orbits play a fundamental role.