Optimization of transfers under constraints on the thrust direction: II

This paper completes the study of optimal transfers with constraints imposed on the thrust vector direction that was opened by paper [1]. The linear inhomogeneous and homogeneous constraints on the thrust direction are considered (specified either by equalities or inequalities), as well as mixed constraints. Some examples of the constraints are presented. A modified method of the transporting trajectory is applied in order to find the optimal transfer under the linear constraints on the thrust direction. This method also gives the necessary condition for a transfer possibility at a given constraint on the thrust direction. A numerical example is considered, in which the propellant consumption is analyzed for the cases of transfers with and without constraints.     

Deflecting APOPHIS with a flotilla of solar shields

The possibility to use the photonic pressure from the Sun for acting upon the orbit of a man-made object is well known.     

Optimization of transfers under constraints on the thrust direction: I

Transfers with a low thrust are considered under constraints imposed on the thrust vector direction. These constraints can be caused by peculiarities of the attitude control system and the mode of stabilization of a spacecraft, and, in the general case, they are functions of the time and state vector. The constraints specified by equalities and inequalities are investigated. It is shown that the optimal thrust is directed along the projection of the Lawden’s primer vector onto the restricting set.     

Analysis of the orbit lifetime of CubeSats in low Earth orbits including periodic variation in drag due to attitude motion

It is estimated that more than 22,300 human-made objects are in orbit around the Earth, with a total mass above 8,400,000 kg. Around 89% of these objects are non-operational and without control, which makes them to be considered orbital debris. These numbers consider only objects with dimensions larger than 10 cm. Besides those numbers, there are also about 2000 operational satellites in orbit nowadays. The space debris represents a hazard to operational satellites and to the space operations. A major concern is that this number is growing, due to new launches and particles generated by collisions. Another important point is that the development of CubeSats has increased exponentially in the last years, increasing the number of objects in space, mainly in the Low Earth Orbits (LEO). Due to the short operational time, CubeSats boost the debris population. One of the requirements for space debris mitigation in LEO is the limitation of the orbital lifetime of the satellites, which needs to be lower than 25 years. However, there are space debris with longer estimated decay time. In LEÓs, the influence of the atmospheric drag is the main orbital perturbation, and is used in maneuvers to increment the losses in the satellite orbital energy, to locate satellites in constellations and to accelerate the decay.The goal of the present research is to study the influence of aerodynamic rotational maneuver in the CubeSat?s orbital lifetime. The rotational axis is orthogonal to the orbital plane of the CubeSat, which generates variations in the ballistic coefficient along the trajectory. The maneuver is proposed to accelerate the decay and to mitigate orbital debris generated by non-operational CubeSats. The panel method is selected to determine the drag coefficient as a function of the flow incident angle and the spinning rate. The pressure distribution is integrated from the satellite faces at hypersonic rarefied flow to calculate the drag coefficient. The mathematical model considers the gravitational potential of the Earth and the deceleration due to drag. To analyze the effects of the rotation during the decay, multiple trajectories were propagated, comparing the results obtained assuming a constant drag coefficient with trajectories where the drag coefficient changes periodically. The initial perigees selected were lower than 400 km of altitude with eccentricities ranging from 0.00 to 0.02. Six values for the angular velocity were applied in the maneuver. The technique of rotating the spacecraft is an interesting solution to increase the orbit decay of a CubeSat without implementing additional de-orbit devices. Significant changes in the decay time are presented due to the increase of the mean drag coefficient calculated by the panel method, when the maneuver is applied, reducing the orbital lifetime, however the results are independent of the angular velocity of the satellite.     

A Modification of the Method of Transporting Trajectory

Further development of an approximate method for optimizing a flight with an ideally controlled small thrust is proposed. The method is based on the employment of the transporting trajectory and considered in [1–3]. A detailed analysis of the means of improving the accuracy of this method suggested in [2, 3] is carried out, and the solution is presented in finite form. The proposed approach is applied to the flights making flybys of many celestial bodies. In the case of small bodies the solution is also obtained in finite form. A numerical example is considered confirming the high efficiency of this method.     

Mapping orbits around the asteroid 2001SN263

The present paper has the goal of mapping orbits, with respect to the perturbations, for a spacecraft traveling around the asteroid 2001SN₂₆₃. This asteroid is a triple system, which center of mass is in an elliptic orbit around the Sun. The perturbations considered in the present model are the ones due to the oblateness of the central body, the gravity field of the two satellite bodies (Beta and Gamma), the Sun, the Moon, the asteroids Vesta, Pallas and Ceres and all the planets of the Solar System. This mapping is important, because it shows the relative importance of each force for a given orbit for the spacecraft, helping to make a decision about which forces need to be included in the model for a given accuracy and nominal orbit. Another important application of this type of mapping is to find orbits that are less perturbed, since it is expected that those orbits have good potential to require a smaller number of station-keeping maneuvers. Simulations under different conditions are made to find those orbits. The main reason to study those trajectories is that, currently, there are several institutions in Brazil studying the possibility to make a mission to send a spacecraft to this asteroid (the so-called ASTER mission), because there are many important scientific studies that can be performed in that system. The results showed that Gamma is the main perturbing body, followed by Beta (10 times smaller) and the group Sun–Mars-oblateness of Alpha, with perturbations 1000 times weaker than the effects of Gamma. The other bodies have perturbations 10⁷ times smaller. The results also showed that circular and polar orbits are less perturbed, when compared to elliptical and equatorial orbits. Regarding the semi-major axis, an internal orbit is the best choice, followed by a larger external orbit. The inclination of the orbit plays an important role, and there are values for the inclination where the perturbations show minimum and maximum values, so it is important to make a good decision on those values.     

Optimization of low-thrust transfers in the three body problem

We consider transfers with low thrust in an arbitrary field of forces. The modified method of transporting trajectory [1–4] is used for optimization of the transfers. The complexity of finding the transporting trajectory of a preset type can be the main obstacle to application of this method. This challenge is solved for the three-body problem in the Hill motion model. Numerical analysis of the method is performed using an example of the transfers to halo-orbits around the solar-terrestrial libration points.     

SMESE: A SMall Explorer for Solar Eruptions

The SMall Explorer for Solar Eruptions (SMESE) mission is a microsatellite proposed by France and China. The payload of SMESE consists of three packages: LYOT (a Lyman imager and a Lyman coronagraph), DESIR (an Infra-Red Telescope working at 35–80 and 100–250 μm), and HEBS (a High-Energy Burst Spectrometer working in X- and γ-rays).

The status of research on flares and coronal mass ejections is briefly reviewed in the context of on-going missions such as SOHO, TRACE and RHESSI. The scientific objectives and the profile of the mission are described. With a launch around 2012–2013, SMESE will provide a unique tool for detecting and understanding eruptions (flares and coronal mass ejections) close to the maximum phase of activity.     

On one approach to the optimization of low-thrust station keeping manoeuvres

A simple and effective mathematical method to calculate optimal station-keeping manoeuvres by means of electric propulsion is suggested. The method is based on a linearization of the satellite motion near a reference orbit. Two versions of the method allow station keeping both in an assigned position and in any position of the orbit. The method is fully analytical for the two-body problem and takes a simpler form for the circular assigned orbit. The suggested method may also be used in the case when constraints are imposed on the thrust direction due to specific features of the satellite stabilization mode. An application of the method to any force field is shown. Illustrative examples of satellite station keeping in a circular orbit are given. Both cases of the station keeping, i.e., in an assigned position and in an assigned orbit, are considered without and with a constraint on the thrust direction.     

Mapping orbits with low station keeping costs for constellations of satellites based on the integral over the time of the perturbing forces

The present paper is concerned with the search for orbits that have potential to require low fuel consumption for station-keeping maneuvers for constellations of satellites. The method used to study this problem is based on the integral over the time of the undesired perturbing forces. This integral measures the change of velocity caused by the perturbation forces acting on the satellite, so mapping orbits that are less perturbed, which generates good candidates for orbits that requires low fuel consumption for station-keeping maneuvers. The integral over the time depends only on the orbit of the spacecraft and the dynamical system considered. The type of engine and the control technique applied to the spacecraft are not considered to search for those orbits. It can be a good strategy to be applied for a first mapping of orbits. For this search, it is analyzed the integral of orbits with different values of the Keplerian elements in order to find the best ones with respect to this criterion. The perturbations considered are the ones caused by the third body, which includes the Sun and the Moon, and the J₂ term of the geopotential. The results presented here show numerical simulations to obtain the integral of those perturbing forces for different orbits. The GPS and the Molniya constellations are used as examples for those calculations.