Sub-Optimal Orbital Maneuvers for Artificial Satellites

The determination of a specific orbit and the procedure to calculate orbital maneuvers of artificial satellites are problems of extreme importance in the study of orbital mechanics. Therefore, the transferring problem of a spaceship from one orbit to another, and the attention due to this subject has in increased during the last years. Many applications can be found in several space activities, for example, to put a satellite in a geostationary orbit, to change the position of a spaceship, to maintain a specific satellite's orbit, in the design of an interplanetary mission, and others. The Brazilian Satellite SCD-1 (Data Collecting Satellite) will be used as example in this paper. It is the first satellite developed entirely in Brazil, and it remains in operation to this date. SCD-1 was designed, developed, built, and tested by Brazilian scientists, engineers, and technicians working at INPE (National Institute for Space Research, and in Brazilian Industries. During the lifetime, it might be necessary do some complementary maneuvers, being this one either an orbital transferring, or just to make periodical corrections. The purpose of transferring problem is to change the position, velocity and the satellite's mass to a new pre determined state. This transfer can be totally linked (in the case of "Rendezvous") or partially free (free time, free final velocity, etc). In the global case, the direction, the orientation and the magnitude of the thrust to be applied must be chosen, respecting the equipment's limit. In order to make this transferring, either sub-optimal or optimal maneuvers may be used. In the present study, only the sub-optimal will be shown. Hence, this method will simplify the direction of thrust application, to allow a fast calculation that may be used in real time, with a very fast processing. The thrust application direction to be applied will be assumed small and constant, and the purpose of this paper is to find the time interval that the thrust is applied. This paper is basically divided into three parts: during the first one the sub-optimal maneuver is explained and detailed, the second presents the Satellite SCD-1, and finally the last part shows the results using the sub-optimal maneuver applied to the Brazilian Satellite.

Near-Space Wide-Swath Radar Imaging With Multiaperture Antenna

Near-space, defined as the altitude region between 20 and 100 km, offers many capabilities that are not accessible for low Earth-orbit (LEO) satellites or airplanes because it is above storm and not constrained by orbital mechanics and high fuel consumption. Hence, a high flying speed can be obtained for the maneuvering vehicles operating in near-space. This offers a promising solution to simultaneous high-resolution and wide-swath synthetic aperture radar (SAR) imaging. As such, one near-space wide-swath SAR imaging technique is presented in this letter. The system configuration, signal model, and imaging scheme are described. An example near-space SAR system is designed, and its imaging performance is analyzed. Simulation results show that near-space maneuvering vehicle SAR indeed seems to be a promising solution to wide-swath SAR imaging.

Families of periodic horseshoe orbits in the restricted three-body problem

We compute families of symmetric periodic horseshoe orbits in the restricted three-body problem. Both the planar and three-dimensional cases are considered and several families are found.We describe how these families are organized as well as the behavior along and among the families of parameters such as the Jacobi constant or the eccentricity. We also determine the stability properties of individual orbits along the families. Interestingly, we find stable horseshoe-shaped orbit up to the quite high inclination of 17◦

On the effects of each term of the geopotential perturbation along the time I: Quasi-circular orbits

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Revisiting the ISN Flow Parameters, using a variable IBEX pointing strategy

The Interstellar Boundary Explorer (IBEX) has observed the interstellar neutral (ISN) gas flow over the past 6 yr during winter/spring when the Earth's motion opposes the ISN flow. Since IBEX observes the interstellar atom trajectories near their perihelion, we can use an analytical model based upon orbital mechanics to determine the interstellar parameters. Interstellar flow latitude, velocity, and temperature are coupled to the flow longitude and are restricted by the IBEX observations to a narrow tube in this parameter space. In our original analysis we found that pointing the spacecraft spin axis slightly out of the ecliptic plane significantly influences the ISN flow vector determination. Introducing the spacecraft spin axis tilt into the analytical model has shown that IBEX observations with various spin axis tilt orientations can substantially reduce the range of acceptable solutions to the ISN flow parameters as a function of flow longitude. The IBEX operations team pointed the IBEX spin axis almost exactly within the ecliptic plane during the 2012-2014 seasons, and about 5° below the ecliptic for half of the 2014 season. In its current implementation the analytical model describes the ISN flow most precisely for the spin axis orientation exactly in the ecliptic. This analysis refines the derived ISN flow parameters with a possible reconciliation between velocity vectors found with IBEX and Ulysses, resulting in a flow longitude lambda∞ = 74.°5 ± 1.°7 and latitude beta∞ = -5.°2 ± 0.°3, but at a substantially higher ISN temperature than previously reported.

Reliability of orbital fits for resonant extrasolar planetary systems: the case of HD82943

We perform a statistical study of the process of orbital determination of the HD82943 extrasolar planetary system, using the current observational data set of N = 165 radial velocity (RV) measurements. Our aim is to analyse the dispersion of possible orbital fits leading to residuals compatible with the best solution, and to discuss the sensitivity of the results with respect to both the data set and the error distribution around the best fit. Although some orbital parameters (e.g. semimajor axis) appear well constrained, we show that the best fits for the HD82943 system are not robust, and at present it is not possible to estimate reliable solutions for these bodies. Finally, we discuss the possibility of a third planet, with a mass of 0.35M(Jup) and an orbital period of 900 d. Stability analysis and simulations of planetary migration indicate that such a hypothetical three-planet system could be locked in a double 2/1 mean-motion resonance, similar to the so-called Laplace resonance of the three inner Galilean satellites of Jupiter.

Dynamics of two planets in co-orbital motion

We study the stability regions and families of periodic orbits of two planets locked in a co-orbital configuration. We consider different ratios of planetary masses and orbital eccentricities; we also assume that both planets share the same orbital plane. Initially, we perform numerical simulations over a grid of osculating initial conditions to map the regions of stable/chaotic motion and identify equilibrium solutions. These results are later analysed in more detail using a semi-analytical model. Apart from the well-known quasi-satellite orbits and the classical equilibrium Lagrangian points L(4) and L(5), we also find a new regime of asymmetric periodic solutions. For low eccentricities these are located at (delta lambda, delta pi) = (+/- 60 degrees, -/+ 120 degrees), where delta lambda is the difference in mean longitudes and delta pi is the difference in longitudes of pericentre. The position of these anti-Lagrangian solutions changes with the mass ratio and the orbital eccentricities and are found for eccentricities as high as similar to 0.7. Finally, we also applied a slow mass variation to one of the planets and analysed its effect on an initially asymmetric periodic orbit. We found that the resonant solution is preserved as long as the mass variation is adiabatic, with practically no change in the equilibrium values of the angles.

Some orbital characteristics of lunar artificial satellites

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Nebular gas drag and co-orbital system dynamics

Aims. We study trajectories of planetesimals whose orbits decay due to gas drag in a primordial solar nebula and are perturbed by the gravity of the secondary body on an eccentric orbit whose mass ratio takes values from mu(2) = 10(-7) to mu(2) = 10(-3) increasing ten times at each step. Each planetesimal ultimately suffers one of the three possible fates: (1) trapping in a mean motion resonance with the secondary body; (2) collision with the secondary body and consequent increase of its mass; or (3) diffusion after crossing the orbit of the secondary body.Methods. We take the Burlirsh-Stoer numerical algorithm in order to integrate the Newtonian equations of the planar, elliptical restricted three-body problem with the secondary body and the planetesimal orbiting the primary. It is assumed that there is no interaction among planetesimals, and also that the gas does not affect the orbit of the secondary body.Results. The results show that the optimal value of the gas drag constant k for the 1: 1 resonance is between 0.9 and 1.25, representing a meter size planetesimal for each AU of orbital radius. In this study, the conditions of the gas drag are such that in theory, L4 no longer exists in the circular case for a critical value of k that defines a limit size of the planetesimal, but for a secondary body with an eccentricity larger than 0.05 when mu(2) = 10(-6), it reappears. The decrease of the cutoff collision radius increase the difusions but does not affect the distribution of trapping. The contribution to the mass accretion of the secondary body is over 40% with a collision radius 0.05R(Hill) and less than 15% with 0.005R(Hill) for mu(2) = 10(-7). The trappings no longer occur when the drag constant k reachs 30. That means that the size limit of planetesimal trapping is 0.2 m per AU of orbital radius. In most cases, this accretion occurs for a weak gas drag and small secondary eccentricity. The diffusions represent most of the simulations showing that gas drag is an efficient process in scattering planetesimals and that the trapping of planetesimals in the 1: 1 resonance is a less probable fate. These results depend on the specific drag force chosen.

Determination of onboard coplanar orbital maneuvers with orbits determined using GPS

In the present paper a study is made in order to find an algorithm that can calculate coplanar orbital maneuvers for an artificial satellite. The idea is to find a method that is fast enough to be combined with onboard orbit determination using GPS data collected from a receiver that is located in the satellite. After a search in the literature, three algorithms are selected to be tested. Preliminary studies show that one of them (the so called Minimum Delta-V Lambert Problem) has several advantages over the two others, both in terms of accuracy and time required for processing. So, this algorithm is implemented and tested numerically combined with the orbit determination procedure. Some adjustments are performed in this algorithm in the present paper to allow its use in real-time onboard applications. Considering the whole maneuver, first of all a simplified and compact algorithm is used to estimate in real-time and onboard the artificial satellite orbit using the GPS measurements. By using the estimated orbit as the initial one and the information of the final desired orbit (from the specification of the mission) as the final one, a coplanar bi-impulsive maneuver is calculated. This maneuver searches for the minimum fuel consumption. Two kinds of maneuvers are performed, one varying only the semi major axis and the other varying the semi major axis and the eccentricity of the orbit, simultaneously. The possibilities of restrictions in the locations to apply the impulses are included, as well as the possibility to control the relation between the processing time and the solution accuracy. Those are the two main reasons to recommend this method for use in the proposed application.

Analytical direct planetary perturbations on the orbital motion of satellites

An analytical expansion of the disturbing function arising from direct planetary perturbations on the motion of satellites is derived. As a Fourier series, it allows the investigation of the secular effects of these direct perturbations, as well as of every argument present in the perturbation. In particular, we construct an analytical model describing the evection resonance between the longitude of pericenter of the satellite orbit and the longitude of a planet, and study briefly its dynamic. The expansion developed in this paper is valid in the case of planar and circular planetary orbits, but not limited in eccentricity or inclination of the satellite orbit. © 2012 Springer Science+Business Media Dordrecht.

Exterior orientation of CBERS-2B imagery using multi-feature control and orbital data

The major contribution of this paper relates to the practical advantages of combining Ground Control Points (GCPs), Ground Control Lines (GCLs) and orbital data to estimate the exterior orientation parameters of images collected by CBERS-2B (China-Brazil Earth Resources Satellite) HRC (High-resolution Camera) and CCD (High-resolution CCD Camera) sensors. Although the CBERS-2B is no longer operational, its images are still being used in Brazil, and the next generations of the CBERS satellite will have sensors with similar technical features, which motivates the study presented in this paper. The mathematical models that relate the object and image spaces are based on collinearity (for points) and coplanarity (for lines) conditions. These models were created in an in-house developed software package called TMS (Triangulation with Multiple Sensors) with multi-feature control (GCPs and GCLs). Experiments on a block of four CBERS-2B HRC images and on one CBERS-2B CCD image were performed using both models. It was observed that the combination of GCPs and GCLs provided better bundle block adjustment results than conventional bundle adjustment using only GCPs. The results also demonstrate the advantages of using primarily orbital data when the number of control entities is reduced. © 2013 International Society for Photogrammetry and Remote Sensing, Inc. (ISPRS).

An analysis of the Hygiea asteroid family orbital region

(10) Hygiea is the fourth largest asteroid of the main belt, by volume and mass, and it is the largest member of its family, that is made mostly by low-albedo, C-type asteroids, typical of the outer main belt. Like many other large families, it is associated with a 'halo' of objects, that extends far beyond the boundary of the core family, as detected by traditional hierarchical clustering methods (HCM) in proper element domains. Numerical simulations of the orbital evolution of family members may help in estimating the family and halo family age, and the original ejection velocity field. But, in order to minimize the errors associated with including too many interlopers, it is important to have good estimates of family membership that include available data on local asteroid taxonomy, geometrical albedo and local dynamics. For this purpose, we obtained synthetic proper elements and frequencies of asteroids in the Hygiea orbital region, with their errors. We revised the current knowledge on asteroid taxonomy, including Sloan Digital Sky Survey-Moving Object Catalog 4th release (SDSS-MOC 4) data, and geometric albedo data from Wide-field Infrared Survey Explorer (WISE) and Near-Earth Object WISE (NEOWISE). We identified asteroid family members using HCM in the domain of proper elements (a, e, sin (i)) and in the domains of proper frequencies most appropriate to study diffusion in the local web of secular resonances, and eliminated possible interlopers based on taxonomic and geometrical albedo considerations. To identify the family halo, we devised a new hierarchical clustering method in an extended domain that includes proper elements, principal components PC1, PC2 obtained based on SDSS photometric data and, for the first time, WISE and NEOWISE geometric albedo. Data on asteroid size distribution, light curves and rotations were also revised for the Hygiea family. The Hygiea family is the largest group in its region, with two smaller families in proper element domain and 18 families in various frequencies domains identified in this work for the first time. Frequency groups tend to extend vertically in the (a, sin (i)) plane and cross not only the Hygiea family but also the near C-type families of Themis and Veritas, causing a mixture of objects all of relatively low albedo in the Hygiea family area. A few high-albedo asteroids, most likely associated with the Eos family, are also present in the region. Finally, the new multidomains hierarchical clustering method allowed us to obtain a good and robust estimate of the membership of the Hygiea family halo, quite separated from other asteroids families halo in the region, and with a very limited (about 3 per cent) presence of likely interlopers. © 2013 The Author Published by Oxford University Press on behalf of the Royal Astronomical Society.

Planet formation in a triple stellar system: implications of the third star's orbital inclination

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Celestial mechanics: from the errant stars to guidance of spacecrafts

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)