Optimal two-impulse rendezvous using constrained multiple-revolution Lambert solutions

A solution to the fixed-time minimum-fuel two-impulse rendezvous problem for the general non-coplanar elliptical orbits is provided. The optimal transfer orbit is obtained using the constrained multiple-revolution Lambert solution. Constraints consist of lower bound for perigee altitude and upper bound for apogee altitude. The optimal time-free two-impulse transfer problem between two fixed endpoints implies finding the roots of an eighth order polynomial, which is done using a numerical iterative technique. The set of feasible solutions is determined by using the constraints conditions to solve for the short-path and long-path orbits semimajor axis ranges. Then, by comparing the optimal time-free solution with the feasible solutions, the optimal semimajor axis for the two fixed-endpoints transfer is identified. Based on the proposed solution procedure for the optimal two fixed-endpoints transfer, a contour of the minimum cost for different initial and final coasting parameters is obtained. Finally, a numerical optimization algorithm (e.g., evolutionary algorithm) can be used to solve this global minimization problem. A numerical example is provided to show how to apply the proposed technique.     

Optimal transfers between unstable periodic orbits using invariant manifolds

This paper presents a method to construct optimal transfers between unstable periodic orbits of differing energies using invariant manifolds. The transfers constructed in this method asymptotically depart the initial orbit on a trajectory contained within the unstable manifold of the initial orbit and later, asymptotically arrive at the final orbit on a trajectory contained within the stable manifold of the final orbit. Primer vector theory is applied to a transfer to determine the optimal maneuvers required to create the bridging trajectory that connects the unstable and stable manifold trajectories. Transfers are constructed between unstable periodic orbits in the Sun–Earth, Earth–Moon, and Jupiter-Europa three-body systems. Multiple solutions are found between the same initial and final orbits, where certain solutions retrace interior portions of the trajectory. All transfers created satisfy the conditions for optimality. The costs of transfers constructed using manifolds are compared to the costs of transfers constructed without the use of manifolds. In all cases, the total cost of the transfer is significantly lower when invariant manifolds are used in the transfer construction. In many cases, the transfers that employ invariant manifolds are three times more efficient, in terms of fuel expenditure, than the transfer that do not. The decrease in transfer cost is accompanied by an increase in transfer time of flight.     

Mars interplanetary trajectory design via Lagrangian points

With the increase in complexities of interplanetary missions, the main focus has shifted to reducing the total delta-V for the entire mission and hence increasing the payload capacity of the spacecraft. This paper develops a trajectory to Mars using the Lagrangian points of the Sun-Earth system and the Sun-Mars system. The whole trajectory can be broadly divided into three stages: (1) Trajectory from a near-Earth circular parking orbit to a halo orbit around Sun-Earth Lagrangian point L₂. (2) Trajectory from Sun-Earth L₂ halo orbit to Sun-Mars L₁ halo orbit. (3) Sun-Mars L₁ halo orbit to a circular orbit around Mars. The stable and unstable manifolds of the halo orbits are used for halo orbit insertion. The intermediate transfer arcs are designed using two-body Lambert’s problem. The total delta-V for the whole trajectory is computed and found to be lesser than that for the conventional trajectories. For a 480 km Earth parking orbit, the total delta-V is found to be 4.6203 km/s. Another advantage in the present approach is that delta-V does not depend upon the synodic period of Earth with respect to Mars.     

A Hamiltonian approach to the planar optimization of mid-course corrections

Lawden’s primer vector theory gives a set of necessary conditions that characterize the optimality of a transfer orbit, defined accordingly to the possibility of adding mid-course corrections. In this paper a novel approach is proposed where, through a polar coordinates transformation, the primer vector components decouple. Furthermore, the case when transfer, departure and arrival orbits are coplanar is analyzed using a Hamiltonian approach. This procedure leads to approximate analytic solutions for the in-plane components of the primer vector. Moreover, the solution for the circular transfer case is proven to be the Hill’s solution. The novel procedure reduces the mathematical and computational complexity of the original case study. It is shown that the primer vector is independent of the semi-major axis of the transfer orbit. The case with a fixed transfer trajectory and variable initial and final thrust impulses is studied. The acquired related optimality maps are presented and analyzed and they express the likelihood of a set of trajectories to be optimal. Furthermore, it is presented which kind of requirements have to be fulfilled by a set of departure and arrival orbits to have the same profile of primer vector.     

Halo orbit transfer trajectory design using invariant manifold in the Sun-Earth system accounting radiation pressure and oblateness

In this paper, we study the invariant manifold and its application in transfer trajectory problem from a low Earth parking orbit to the Sun-Earth \(L_{1}\) and \(L_{2}\)-halo orbits with the inclusion of radiation pressure and oblateness. Invariant manifold of the halo orbit provides a natural entrance to travel the spacecraft in the solar system along some specific paths due to its strong hyperbolic character. In this regard, the halo orbits near both collinear Lagrangian points are computed first. The manifold’s approximation near the nominal halo orbit is computed using the eigenvectors of the monodromy matrix. The obtained local approximation provides globalization of the manifold by applying backward time propagation to the governing equations of motion. The desired transfer trajectory well suited for the transfer is explored by looking at a possible intersection between the Earth’s parking orbit of the spacecraft and the manifold.     

A multi-objective approach to the design of low thrust space trajectories using optimal control

In this article, we introduce a novel three-step approach for solving optimal control problems in space mission design. We demonstrate its potential by the example task of sending a group of spacecraft to a specific Earth L ₂ halo orbit. In each of the three steps we make use of recently developed optimization methods and the result of one step serves as input data for the subsequent one. Firstly, we perform a global and multi-objective optimization on a restricted class of control functions. The solutions of this problem are (Pareto-)optimal with respect to ΔV and flight time. Based on the solution set, a compromise trajectory can be chosen suited to the mission goals. In the second step, this selected trajectory serves as initial guess for a direct local optimization. We construct a trajectory using a more flexible control law and, hence, the obtained solutions are improved with respect to control effort. Finally, we consider the improved result as a reference trajectory for a formation flight task and compute trajectories for several spacecraft such that these arrive at the halo orbit in a prescribed relative configuration. The strong points of our three-step approach are that the challenging design of good initial guesses is handled numerically by the global optimization tool and afterwards, the last two steps only have to be performed for one reference trajectory.     

The use of invariant manifolds for transfers between unstable periodic orbits of different energies

Techniques from dynamical systems theory have been applied to the construction of transfers between unstable periodic orbits that have different energies. Invariant manifolds, trajectories that asymptotically depart or approach unstable periodic orbits, are used to connect the initial and final orbits. The transfer asymptotically departs the initial orbit on a trajectory contained within the initial orbit’s unstable manifold and later asymptotically approaches the final orbit on a trajectory contained within the stable manifold of the final orbit. The manifold trajectories are connected by the execution of impulsive maneuvers. Two-body parameters dictate the selection of the individual manifold trajectories used to construct efficient transfers. A bounding sphere centered on the secondary, with a radius less than the sphere of influence of the secondary, is used to study the manifold trajectories. A two-body parameter, κ, is computed within the bounding sphere, where the gravitational effects of the secondary dominate. The parameter κ is defined as the sum of two quantities: the difference in the normalized angular momentum vectors and eccentricity vectors between a point on the unstable manifold and a point on the stable manifold. It is numerically demonstrated that as the κ parameter decreases, the total cost to complete the transfer decreases. Preliminary results indicate that this method of constructing transfers produces a significant cost savings over methods that do not employ the use of invariant manifolds.     

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.     

Lambert problem solution in the hill model of motion

The goal of this paper is obtaining a solution of the Lambert problem in the restricted three-body problem described by the Hill equations. This solution is based on the use of pre determinate reference orbits of different types giving the first guess and defining the sought-for transfer type. A mathematical procedure giving the Lambert problem solution is described. This procedure provides step-by-step transformation of the reference orbit to the sought-for transfer orbit. Numerical examples of the procedure application to the transfers in the Sun–Earth system are considered. These examples include transfer between two specified positions in a given time, a periodic orbit design, a halo orbit design, halo-to-halo transfers, LEO-to-halo transfer, analysis of a family of the halo-to-halo transfer orbits. The proposed method of the Lambert problem solution can be used for the two-point boundary value problem solution in any model of motion if a set of typical reference orbits can be found.     

Three-dimensional time optimal double angular momentum reversal trajectory using solar sails

A new concept of three dimensional non-Keplerian trajectories with double angular momentum reversal is investigated with high performance solar sails. The main discussion of this paper is about such 3D solar inverse orbits with inner constraints. The problem is addressed in a time optimal control framework solved by an indirect method. Two typical solar inverse orbits have been achieved and presented in a 3D non-dimensional dynamic model in the Heliocentric Inertial Frame. Starting from the Earth orbit ecliptic plane, a sailcraft in the inverse orbit exhibits a butterfly shape trajectory. As such, the new orbits are symmetrical with respect to a plane which contains the Sun-perihelion line. The relation of the sail attitude angles between the two symmetrical parts of the orbits are used to reduce the simulation effort. The quasi-heliostationary property at its aphelia is demonstrated with variation of the orbital radius. Evolutions of the orbital velocity and optimal sail orientations are also outlined and discussed to benefit future design work. As is suited for space observation guaranteed by its butterfly shape, the inverse orbits are thoroughly studied in terms of the concerned parameters. The discussion of the parametric influence is ranked in order as perihelion distance r _E , required maximum position z _max, perihelion position z _f and the sail lightness number β. Suitable ranges of each parameter are adopted to illustrate the orbital variation trend. Through numerical simulations the features of such inverse orbits are further emphasized to provide an initial reference for future researchers.     

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.     

Transfers between libration-point orbits in the elliptic restricted problem

A strategy is formulated to design optimal time-fixed impulsive transfers between three-dimensional libration-point orbits in the vicinity of the interiorL ₁ libration point of the Sun-Earth/Moon barycenter system. The adjoint equation in terms of rotating coordinates in the elliptic restricted three-body problem is shown to be of a distinctly different form from that obtained in the analysis of trajectories in the two-body problem. Also, the necessary conditions for a time-fixed two-impulse transfer to be optimal are stated in terms of the primer vector. Primer vector theory is then extended to non-optimal impulsive trajectories in order to establish a criterion whereby the addition of an interior impulse reduces total fuel expenditure. The necessary conditions for the local optimality of a transfer containing additional impulses are satisfied by requiring continuity of the Hamiltonian and the derivative of the primer vector at all interior impulses. Determination of the location, orientation, and magnitude of each additional impulse is accomplished by the unconstrained minimization of the cost function using a multivariable search method. Results indicate that substantial savings in fuel can be achieved by the addition of interior impulsive maneuvers on transfers between libration-point orbits.An earlier version was presented as Paper AAS 92–126 at the AAS/AIAA Spaceflight Mechanics Meeting, Colorado Springs, Colorado, February 24–26, 1992.     

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.     

Periodic orbits of solar sail equipped with reflectance control device in Earth–Moon system

In this paper, families of Lyapunov and halo orbits are presented with a solar sail equipped with a reflectance control device in the Earth–Moon system. System dynamical model is established considering solar sail acceleration, and four solar sail steering laws and two initial Sun-sail configurations are introduced. The initial natural periodic orbits with suitable periods are firstly identified. Subsequently, families of solar sail Lyapunov and halo orbits around the \(L_{1}\) and \(L_{2}\) points are designed with fixed solar sail characteristic acceleration and varying reflectivity rate and pitching angle by the combination of the modified differential correction method and continuation approach. The linear stabilities of solar sail periodic orbits are investigated, and a nonlinear sliding model controller is designed for station keeping. In addition, orbit transfer between the same family of solar sail orbits is investigated preliminarily to showcase reflectance control device solar sail maneuver capability.     

The haloes of merger remnants

We perform collisionless N -body simulations of 1:1 galaxy mergers, using models which include a galaxy halo, disc and bulge, focusing on the behaviour of the halo component. The galaxy models are constructed without recourse to a Maxwellian approximation. We investigate the effect of varying the galaxies' orientation, their mutual orbit and the initial velocity anisotropy or cusp strength of the haloes upon the remnant halo density profiles and shape, as well as on the kinematics. We observe that the halo density profile (determined as a spherical average, an approximation we find appropriate) is exceptionally robust in mergers, and that the velocity anisotropy of our remnant haloes is nearly independent of the orbits or initial anisotropy of the haloes. The remnants follow the halo anisotropy – local density slope (β–γ) relation suggested by Hansen & Moore in the inner parts of the halo, but β is systematically lower than this relation predicts in the outer parts. Remnant halo axis ratios are strongly dependent on the initial parameters of the haloes and on their orbits. We also find that the remnant haloes are significantly less spherical than those described in studies of simulations which include gas cooling.     

Study of the transfer from the Earth to a halo orbit around the equilibrium pointL 1

The purpose of this paper is to study a transfer strategy from the vicinity of the Earth to a halo orbit around the equilibrium pointL ₁ of the Earth-Sun system. The study is done in the real solar system (we use the DE-118 JPL ephemeris in the simulations of motion) although some simplified models, such as the restricted three body problem (RTBP) and the bicircular problem, have been also used in order to clarify the geometrical aspects of the problem. The approach used in the paper makes use of the hyperbolic character of the halo orbits under consideration. The invariant stable manifold of these orbits enables the transfer to be achieved with, theoretically, only one manoeuvre: the one of insertion into the stable manifold. For the total v required, the figures obtained are similar to the ones given by the standard procedures of optimization.     

A closed-form solution to the minimum $${\Delta V_{\rm tot}^2}$$ Lambert’s problem

A closed form solution to the minimum DV_tot²{\Delta V_{\rm tot}^2} Lambert problem between two assigned positions in two distinct orbits is presented. Motivation comes from the need of computing optimal orbit transfer matrices to solve re-configuration problems of satellite constellations and the complexity associated in facing these problems with the minimization of DV_tot{\Delta V_{\rm tot}}. Extensive numerical tests show that the difference in fuel consumption between the solutions obtained by minimizing DV_tot²{\Delta V_{\rm tot}^2} and DV_tot{\Delta V_{\rm tot}} does not exceed 17%. The DV_tot²{\Delta V_{\rm tot}^2} solution can be adopted as starting point to find the minimum DV_tot{\Delta V_{\rm tot}}. The solving equation for minimum DV_tot²{\Delta V_{\rm tot}^2} Lambert problem is a quartic polynomial in term of the angular momentum modulus of the optimal transfer orbit. The root selection is discussed and the singular case, occurring when the initial and final radii are parallel, is analytically solved. A numerical example for the general case (orbit transfer “pork-chop” between two non-coplanar elliptical orbits) and two examples for the singular case (Hohmann and GTO transfers) are provided.     

Solar Sail Halo Orbits at the Sun–Earth Artificial L 1 Point

Halo orbits for solar sails at artificial Sun–Earth L₁ points are investigated by a third order approximate solution. Two families of halo orbits are explored as defined by the sail attitude. Case I: the sail normal is directed along the Sun-sail line. Case II: the sail normal is directed along the Sun–Earth line. In both cases the minimum amplitude of a halo orbit increases as the lightness number of the solar sail increases. The effect of the z-direction amplitude on x- or y-direction amplitude is also investigated and the results show that the effect is relatively small. In case I, the orbit period increases as the sail lightness number increases, while in case II, as the lightness number increases, the orbit period increases first and then decreases after the lightness number exceeds ~0.01.     

Superosculating intermediate orbits: Fourth-and fifth-order tangencies to trajectories of perturbed motion

The theory of superosculating intermediate orbits previously suggested by the author is developed. A new class of orbits with a fourth-order tangency to the actual trajectory of a celestial body at the initial time is constructed. Orbits with a fifth-order tangency have been constructed for the first time. The motion in the constructed orbits is represented as a combination of two motions: the motion of a fictitious attracting center with a variable mass and the motion relative to this center. The first motion is generally parabolic, while the second motion is described by the equations of the Gylden—Mestschersky problem. The variation in the mass of the fictitious center obeys Mestschersky’s first and combined laws. The new orbits represent more accurately the actual motion in the initial segment of the trajectory than an osculating Keplerian orbit and other existing analogues. Encke’s generalized methods of special perturbations in which the constructed intermediate orbits are used as reference orbits are presented. Numerical simulations using the approximations of the motions of Asteroid Toutatis and Comet P/Honda—Mrkos—Pajdu?áková as examples confirm that the constructed orbits are highly efficient. Their application is particularly beneficial in investigating strongly perturbed motion.     

Optimum impulsive transfers between elliptic and non-coplanar circular orbits

A study has been made of optimum transfers between elliptic and non-coplanar circular orbits having a common centre of attraction and whose planes intersect along the major axis of the ellipse. Elliptic transfer paths with up to three apsidal impulses are considered, with the whole plane change taking place at the (coincident) apocentres of these paths.

For three-impulse transfers, the optimum mode is always to transfer from the pericentre of the elliptic orbit to the circular orbit, or vice-versa. For the two-impulse “tilted-Hohmann” type of transfer, the optimum mode for the ellipse-in-circle arrangement of initial and final orbits is also to transfer from pericentre to circle; but the optimum mode is from apocentre to circle for both the circle-in-ellipse and overlap arrangements.