Artificial halo orbits for low-thrust propulsion spacecraft

We consider periodic halo orbits about artificial equilibrium points (AEP) near to the Lagrange points L ₁ and L ₂ in the circular restricted three body problem, where the third body is a low-thrust propulsion spacecraft in the Sun–Earth system. Although such halo orbits about artificial equilibrium points can be generated using a solar sail, there are points inside L ₁ and beyond L ₂ where a solar sail cannot be placed, so low-thrust, such as solar electric propulsion, is the only option to generate artificial halo orbits around points inaccessible to a solar sail. Analytical and numerical halo orbits for such low-thrust propulsion systems are obtained by using the Lindstedt Poincaré and differential corrector method respectively. Both the period and minimum amplitude of halo orbits about artificial equilibrium points inside L ₁ decreases with an increase in low-thrust acceleration. The halo orbits about artificial equilibrium points beyond L ₂ in contrast show an increase in period with an increase in low-thrust acceleration. However, the minimum amplitude first increases and then decreases after the thrust acceleration exceeds 0.415 mm/s². Using a continuation method, we also find stable artificial halo orbits which can be sustained for long integration times and require a reasonably small low-thrust acceleration 0.0593 mm/s².     

Selection and trajectory design to mission secondary targets

Recently, with new trajectory design techniques and use of low-thrust propulsion systems, missions have become more efficient and cheaper with respect to propellant. As a way to increase the mission’s value and scientific return, secondary targets close to the main trajectory are often added with a small change in the transfer trajectory. As a result of their large number, importance and facility to perform a flyby, asteroids are commonly used as such targets. This work uses the Primer Vector theory to define the direction and magnitude of the thrust for a minimum fuel consumption problem. The design of a low-thrust trajectory with a midcourse asteroid flyby is not only challenging for the low-thrust problem solution, but also with respect to the selection of a target and its flyby point. Currently more than 700,000 minor bodies have been identified, which generates a very large number of possible flyby points. This work uses a combination of reachability, reference orbit, and linear theory to select appropriate candidates, drastically reducing the simulation time, to be later included in the main trajectory and optimized. Two test cases are presented using the aforementioned selection process and optimization to add and design a secondary flyby to a mission with the primary objective of 3200 Phaethon flyby and 25143 Itokawa rendezvous.     

Periodic orbits in the photogravitational restricted problem with the smaller primary an oblate body

In this paper, we have studied periodic orbits generated by Lagrangian solutions of the restricted three body problem when more massive body is a source of radiation and the smaller primary is an oblate body. We have determined periodic orbits for fixed values of μ, σ and different values of p and h (μ mass ratio of the two primaries, σ oblate parameter, p radiation parameter and h energy constant). These orbits have been determined by giving displacements along the tangent and normal to the mobile co-ordinates as defined by Karimov and Sokolsky (in Celest. Mech. 46:335, 1989). These orbits have been drawn by using the predictor-corrector method. We have also studied the effect of radiation pressure on the periodic orbits by taking some fixed values of μ and σ.     

Periodic orbits generated by Lagrangian solutions of the restricted three body problem when one of the primaries is an oblate body

We have studied periodic orbits generated by Lagrangian solutions of the restricted three body problem when one of the primaries is an oblate body. We have determined the periodic orbits for different values of μ, h and A (h is energy constant, μ is mass ratio of the two primaries and A is an oblateness factor). These orbits have been determined by giving displacements along the tangent and normal to the mobile coordinates as defined by Karimov and Sokolsky (Celest. Mech. 46:335, 1989). These orbits have been drawn by using the predictor-corrector method. We have also studied the effect of oblateness by taking some fixed values of μ, A and h. As starters for our method, we use some known periodic orbits in the classical restricted three body problem.     

Stable Manifolds and Homoclinic Points Near Resonances in the Restricted Three-Body Problem

The restricted three-body problem describes the motion of a massless particle under the influence of two primaries of masses 1− μ and μ that circle each other with period equal to 2π. For small μ, a resonant periodic motion of the massless particle in the rotating frame can be described by relatively prime integers p and q, if its period around the heavier primary is approximately 2π p/q, and by its approximate eccentricity e. We give a method for the formal development of the stable and unstable manifolds associated with these resonant motions. We prove the validity of this formal development and the existence of homoclinic points in the resonant region. In the study of the Kirkwood gaps in the asteroid belt, the separatrices of the averaged equations of the restricted three-body problem are commonly used to derive analytical approximations to the boundaries of the resonances. We use the unaveraged equations to find values of asteroid eccentricity below which these approximations will not hold for the Kirkwood gaps with q/p equal to 2/1, 7/3, 5/2, 3/1, and 4/1. Another application is to the existence of asymmetric librations in the exterior resonances. We give values of asteroid eccentricity below which asymmetric librations will not exist for the 1/7, 1/6, 1/5, 1/4, 1/3, and 1/2 resonances for any μ however small. But if the eccentricity exceeds these thresholds, asymmetric librations will exist for μ small enough in the unaveraged restricted three-body problem.     

Combined low-thrust propulsion and invariant manifold trajectories to capture NEOs in the Sun–Earth circular restricted three-body problem

In this paper, a method to capture near-Earth objects (NEOs) incorporating low-thrust propulsion into the invariant manifolds technique is investigated. Assuming that a tugboat-spacecraft is in a rendez-vous condition with the candidate asteroid, the aim is to take the joint spacecraft-asteroid system to a selected periodic orbit of the Sun–Earth restricted three-body system: the orbit can be either a libration point periodic orbit (LPO) or a distant prograde periodic orbit (DPO) around the Earth. In detail, low-thrust propulsion is used to bring the joint spacecraft-asteroid system from the initial condition to a point belonging to the stable manifold associated to the final periodic orbit: from here onward, thanks to the intrinsic dynamics of the physical model adopted, the flight is purely ballistic. Dedicated guided and capture sets are introduced to exploit the combined use of low-thrust propulsion with stable manifolds trajectories, aiming at defining feasible first guess solutions. Then, an optimal control problem is formulated to refine and improve them. This approach enables a new class of missions, whose solutions are not obtainable neither through the patched-conics method nor through the classic invariant manifolds technique.     

Cosmic‐ray exposure age and heliocentric distance of the parent body of the Rumuruti chondrite PRE 95410

We measured concentrations and isotopic ratios of noble gases in the Rumuruti (R) chondrite Mount Prestrud (PRE) 95410, a regolith breccia exhibiting dark/light structures. The meteorite contains solar and cosmogenic noble gases. Based on the solar and cosmogenic noble gas compositions, we calculated a heliocentric distance of its parent body, a cosmic‐ray exposure age on the parent body regolith (parent body exposure age), and a cosmic‐ray exposure age in interplanetary space (space exposure age) of the meteorite. Assuming a constant solar wind flux, the estimated heliocentric distance was smaller than 1.4 ± 0.3 au, suggesting inward migration from the asteroid belt regions where the parent body formed. The largest known Mars Trojan 5261 Eureka is a potential parent body of PRE 95410. Alternatively, it is possible that the solar wind flux at the time of the parent body exposure was higher by a factor of 2–3 compared to the lunar regolith exposure. In this case, the estimated heliocentric distance is within the asteroid belt region. The parent body exposure age is longer than 19.1 Ma. This result indicates frequent impact events on the parent body like that recorded for other solar‐gas‐rich meteorites. Assuming single‐stage exposure after an ejection event from the parent body, the space exposure age is 11.0 ± 1.1 Ma, which is close to the peak of ~10 Ma in the exposure age distribution for the solar‐gas‐free R chondrites.     

Hill stability of the satellite in the elliptic restricted four-body problem

We consider an elliptic restricted four-body system including three primaries and a massless particle. The orbits of the primaries are elliptic, and the massless particle moves under the mutual gravitational attraction. From the dynamic equations, a quasi-integral is obtained, which is similar to the Jacobi integral in the circular restricted three-body problem (CRTBP). The energy constant \(C\) determines the topology of zero velocity surfaces, which bifurcate at the equilibrium point. We define the concept of Hill stability in this problem, and a criterion for stability is deduced. If the actual energy constant \(C_{\mathrm{ac}}\ ( {>} 0 ) \) is bigger than or equal to the critical energy constant \(C_{\mathrm{cr}}\), the particle will be Hill stable. The critical energy constant is determined by the mass and orbits of the primaries. The criterion provides a way to capture an asteroid into the Earth–Moon system.     

Analysis of impulsive maneuvers to keep orbits around the asteroid 2001SN<Subscript>263</Subscript>

The strongly perturbed environment of a small body, such as an asteroid, can complicate the prediction of orbits used for close proximity operations. Inaccurate predictions may make the spacecraft collide with the asteroid or escape to the deep space. The main forces acting in the dynamics come from the solar radiation pressure and from the body’s weak gravity field. This paper investigates the feasibility of using bi-impulsive maneuvers to avoid the aforementioned non-desired phenomena (collisions and escapes) by connecting orbits around the triple system asteroid 2001SN₂₆₃, which is the target of a proposed Brazilian space mission. In terms of a mathematical formulation, a recently presented rotating dipole model is considered with oblateness in both primaries. In addition, a “two-point boundary value problem” is solved to find a proper transfer trajectory. The results presented here give support to identifying the best strategy to find orbits for close proximity operations, in terms of long orbital lifetimes and low delta-\(V\) consumptions. Numerical results have also demonstrated the significant influence of the spacecraft orbital elements (semi-major axis and eccentricity), angular position of the Sun and spacecraft area-to-mass ratio, in the performance of the bi-impulsive maneuver.     

Accurate analytical approximation of asteroid deflection with constant tangential thrust

We present analytical formulas to estimate the variation of achieved deflection for an Earth-impacting asteroid following a continuous tangential low-thrust deflection strategy. Relatively simple analytical expressions are obtained with the aid of asymptotic theory and the use of Peláez orbital elements set, an approach that is particularly suitable to the asteroid deflection problem and is not limited to small eccentricities. The accuracy of the proposed formulas is evaluated numerically showing negligible error for both early and late deflection campaigns. The results will be of aid in planning future low-thrust asteroid deflection missions.     

THE PHASE SPACE STRUCTURE OF THE EXTENDED SITNIKOV-PROBLEM

In this article we treat the 'Extended Sitnikov Problem' where three bodies of equal masses stay always in the Sitnikov configuration. One of the bodies is confined to a motion perpendicular to the instantaneous plane of motion of the two other bodies (called the primaries), which are always equally far away from the barycenter of the system (and from the third body). In contrary to the Sitnikov Problem with one mass less body the primaries are not moving on Keplerian orbits. After a qualitative analysis of possible motions in the 'Extended Sitnikov Problem' we explore the structure of phase space with the aid of properly chosen surfaces of section. It turns out that for very small energies H the motion is possible only in small region of phase space and only thin layers of chaos appear in this region of mostly regular motion. We have chosen the plane ( ) as surface of section, where r is the distance between the primaries; we plot the respective points when the three bodies are 'aligned'. The fixed point which corresponds to the 1 : 2 resonant orbit between the primaries' period and the period of motion of the third mass is in the middle of the region of motion. For low energies this fixed point is stable, then for an increased value of the energy splits into an unstable and two stable fixed points. The unstable fixed point splits again for larger energies into a stable and two unstable ones. For energies close toH = 0 the stable center splits one more time into an unstable and two stable ones. With increasing energy more and more of the phase space is filled with chaotic orbits with very long intermediate time intervals in between two crossings of the surface of section. We also checked the rotation numbers for some specific orbits. This revised version was published online in July 2006 with corrections to the Cover Date.     

Non-linear stability in photogravitational non-planar restricted three body problem with oblate smaller primary

We have discussed non-linear stability in photogravitational non-planar restricted three body problem with oblate smaller primary. By photogravitational we mean that both primaries are radiating. We normalized the Hamiltonian using Lie transform as in Coppola and Rand (Celest. Mech. 45:103, 1989). We transformed the system into Birkhoff’s normal form. Lie transforms reduce the system to an equivalent simpler system which is immediately solvable. Applying Arnold’s theorem, we have found non-linear stability criteria. We conclude that L ₆ is stable. We plotted graphs for (ω ₁,D ₂). They are rectangular hyperbola.     

Periodic orbits in the generalized perturbed restricted three-body problem

This paper studies the existence and stability of equilibrium points under the influence of small perturbations in the Coriolis and the centrifugal forces, together with the non-sphericity of the primaries. The problem is generalized in the sense that the bigger and smaller primaries are respectively triaxial and oblate spheroidal bodies. It is found that the locations of equilibrium points are affected by the non-sphericity of the bodies and the change in the centrifugal force. It is also seen that the triangular points are stable for 0<μ<μ _c and unstable for m_c ￡ m < \frac12\mu_{c}\le\mu <\frac{1}{2}, where μ _c is the critical mass parameter depending on the above perturbations, triaxiality and oblateness. It is further observed that collinear points remain unstable.     

Binary near-Earth asteroid formation: Rubble pile model of tidal disruptions

We present numerical simulations of near-Earth asteroid (NEA) tidal disruption resulting in bound, mutually orbiting systems. Using a rubble pile model we have constrained the relative likelihoods for possible physical and dynamical properties of the binaries created. Overall 110,500 simulations were run, with each body consisting of ∼1000 particles. The encounter parameters of close approach distance and velocity were varied, as were the bodies' spin, elongation and spin axis direction. The binary production rate increases for closer encounters, at lower speeds, for more elongated bodies, and for bodies with greater spin. The semimajor axes for resultant binaries are peaked between 5 to 20 primary radii, and there is an overall trend for high eccentricity, with 97% of binaries having e > 0.1. The secondary-to-primary size ratios of the simulated binaries are peaked between 0.1 and 0.2, similar to trends among observed asteroid binaries. The spin rates of the primary bodies are narrowly distributed between 3.5- and 6-h periods, whereas the secondaries' periods are more evenly distributed and can exceed 15-h periods. The spin axes of the primary bodies are very closely aligned with the angular momenta of the binary orbits, whereas the secondary spin axes are nearly random. The shapes of the primaries show a large distribution of axis ratios, where those with low elongation (ratio of long and short axis) are both oblate and prolate, and nearly all with large elongation are prolate. This work presents results that suggest tidal disruption of gravitational aggregates can make binaries physically similar to those currently observed in the NEA population. As well, tidal disruption may create an equal number of binaries with qualities different from those observed, mostly binaries with large separation and with elongated primaries.     

Detection by MEGNO of the gravitational resonances between a rotating ellipsoid and a point mass satellite

Nowadays the scientific community considers that more than a third of the asteroids are double. The study of the stability of these systems is quite complex, because of their irregular shapes and tumbling rotations, and requires a full body–full body approach. A particular case is analysed here, when the secondary body is sufficiently small and distant from the primary to be considered as a point mass satellite. Gravitational resonances (between the revolution of the satellite and the rotation of the asteroid) of a small body in fast or slow rotation around a rigid ellipsoid are studied. The same model can be used for the motion of a probe around an irregular asteroid. The gravitational potential induced by the primary body is modelled by the MacMillan potential. The stability of the satellite is measured thanks to the MEGNO indicator (Mean Exponential Growth Factor of Nearby Orbits). We present stability maps in the plane (\fracbd, \fraccd){\left(\frac{b}{d}, \frac{c}{d}\right)} where d, b, and c are the three semi-axes of the ellipsoid shaping the asteroid. Special stable conic-like curves are detected on these maps and explained by an analytical model, based on a simplification of the MacMillan potential for some specific resonances (1 : 1 and 2 : 1). The efficiency of the MEGNO to detect stability is confirmed.     

Periodic orbits of the elliptic restricted problem for the Sun-Jupiter-Saturn system

A systematic approach to generate periodic orbits in the elliptic restricted problem of three bodies in introduced. The approach is based on (numerical) continuation from periodic orbits of the first and second kind in the circular restricted problem to periodic orbits in the elliptic restricted problem. Two families of periodic orbits of the elliptic restricted problem are found by this approach. The mass ratio of the primaries of these orbits is equal to that of the Sun-Jupiter system. The sidereal mean motions between the infinitesimal body and the smaller primary are in a 2:5 resonance, so as to approximate the Sun-Jupiter-Saturn system. The linear stability of these periodic orbits are studied as functions of the eccentricities of the primaries and of the infinitesimal body. The results show that both stable and unstable periodic orbits exist in the elliptic restricted problem that are close to the actual Sun-Jupiter-Saturn system. However, the periodic orbit closest to the actual Sun-Jupiter-Saturn system is (linearly) stable.     

Optimal low-thrust transfers between libration point orbits

Over the past three decades, ballistic and impulsive trajectories between libration point orbits (LPOs) in the Sun–Earth–Moon system have been investigated to a large extent. It is known that coupling invariant manifolds of LPOs of two different circular restricted three-body problems (i.e., the Sun–Earth and the Earth–Moon systems) can lead to significant mass savings in specific transfers, such as from a low Earth orbit to the Moon’s vicinity. Previous investigations on this issue mainly considered the use of impulsive maneuvers along the trajectory. Here we investigate the dynamical effects of replacing impulsive ΔV’s with low-thrust trajectory arcs to connect LPOs using invariant manifold dynamics. Our investigation shows that the use of low-thrust propulsion in a particular phase of the transfer and the adoption of a more realistic Sun–Earth–Moon four-body model can provide better and more propellant-efficient solution. For this purpose, methods have been developed to compute the invariant tori and their manifolds in this dynamical model.     

Parametric dependence of the stationary solutions in the restricted 2 + 2 body problem

The restricted gravitational 2 + 2 body problem, is a particular case of the N body problem and it may be used to approximate the dynamical behaviour of binary asteroids or dual sattelites moving in the gravitational field of two primaries Pi, i = 1,2. By considering oblate primaries, five parameters are needed to describe the model, namely the reduced mass μ of the primary P2, the reduced masses μ1 and μ2 of the minor bodies and the oblatenesses Ii, i = 1,2 of the primaries. This work deals with the effect of those parameters on the location of the stationary solutions. This revised version was published online in July 2006 with corrections to the Cover Date.     

Analysis on the stability of triangular points in the perturbed photogravitational restricted three-body problem with variable masses

This paper examines the effect of a constant κ of a particular integral of the Gylden-Meshcherskii problem on the stability of the triangular points in the restricted three-body problem under the influence of small perturbations in the Coriolis and centrifugal forces, together with the effects of radiation pressure of the bigger primary, when the masses of the primaries vary in accordance with the unified Meshcherskii law. The triangular points of the autonomized system are found to be conditionally stable due to κ. We observed further that the stabilizing or destabilizing tendency of the Coriolis and centrifugal forces is controlled by κ, while the destabilizing effects of the radiation pressure remain unchanged but can be made strong or weak due to κ. The condition that the region of stability is increasing, decreasing or does not exist depend on this constant. The motion around the triangular points L _4,5 varying with time is studied using the Lyapunov Characteristic Numbers, and are found to be generally unstable.