Accurate computation of highly eccentric satellite orbits

For computing highly eccentric (e0.9) Earth satellite orbits with special perturbation methods, a comparison is made between different schemes, namely the direct integration of the equations of motion in Cartesian coordinates, changes of the independent variable, use of a time element, stabilization and use of regular elements. A one-step and a multi-step integration are also compared.It is shown that stabilization and regularization procedures are very helpful for non or smoothly perturbed orbits. In practical cases for space research where all perturbations are considered, these procedures are no longer so efficient. The recommended method in these cases is a multi-step integration of the Cartesian coordinates with a change of the independent variable defining an analytical step size regulation. However, the use of a time element and a stabilization procedure for the equations of motion improves the accuracy, except when a small step size is chosen.     

Analytical linear perturbation theory for highly eccentric satellite orbits

Kaula's satellite linear perturbation theory has been extended for the case of highly eccentric orbits by using elliptic function expansions.On leave from Institute of Applied Astronomy, St.-Petersburg     

Length of arc as independent argument for highly eccentric orbits

For analytic step regulation in numerical integration of highly eccentric orbits it is proposed to use the orbital arc length of a moving particle as independent argument.     

Tidal interactions and disruptions of giant planets on highly eccentric orbits

We calculate the evolution of planets undergoing a strong tidal encounter using smoothed particle hydrodynamics (SPH), for a range of periastron separations. We find that outside the Roche limit, the evolution of the planet is well-described by the standard model of linear, non-radial, adiabatic oscillations. If the planet passes within the Roche limit at periastron, however, mass can be stripped from it, but in no case do we find enough energy transferred to the planet to lead to complete disruption. In light of the three new extrasolar planets discovered with periods shorter than two days, we argue that the shortest-period cases observed in the period-mass relation may be explained by a model whereby planets undergo strong tidal encounters with stars, after either being scattered by dynamical interactions into highly eccentric orbits, or tidally captured from nearly parabolic orbits. Although this scenario does provide a natural explanation for the edge found for planets at twice the Roche limit, it does not explain how such planets will survive the inevitable expansion that results from energy injection during tidal circularization.     

Charged-particle induced background expected in X-ray detecting CCDs in highly eccentric earth orbits

Charged particle fluxes on satellites in highly eccentric Earth orbits are studied, using data from ISEE1 and EXOSAT missions. The effects of minimum-ionizing particles, low-energy electrons and secondary radiation produced by these particles on CCDs used as X-ray detectors are identified. A technique to discriminate events caused by minimum-ionizing particles from X-ray events is proposed. The need for a deflecting system against low-energy electrons with focusing X-ray optics is demonstrated. An expected secondary radiation spectrum behind a typical detector shielding based on Monte-Carlo-simulations is presented. For the EPIC detectors on XMM, detection rates due to this radiation environment are estimated.     

Optimal periodic relative orbit and rectilinear relative orbits with eccentric reference orbits

The problem of two-body linearized periodic relative orbits with eccentric reference orbits is studied in this paper. The periodic relative orbit in the target-orbital coordinate system can be used in fly-around and formation-flying orbit design. Based on the closed-form solutions to the Tschauner–Hempel equations, the initial condition for periodic relative orbits is obtained. Then the minimum-fuel periodic-orbit condition with a single impulse is analytically derived for given initial position and velocity vectors. When considering the initial coasting time, the impulse position of the global minimum-fuel periodic orbit is proved to be near to the perigee of the target and can be obtained by numerical optimization algorithms. Moreover, the condition for a special periodic orbit, i.e., the rectilinear relative orbit in the target-orbital frame, is obtained. Numerical simulations are used to demonstrate the efficacy of the method, and show the geometry of the periodic relative orbit and the rectilinear relative orbit.     

Prediction of satellites in Earth's gravitational field with axial symmetry using Burdet's regularized theory

In this paper, economical and stable recurrence formulae for the Earth's zonal potential and its gradient for Burdet's regularized theory will be established for any number N of the zonal harmonic coefficients. A general recursive computational algorithm based on these formulae is also established for the initial value problem of Burdet oscillator for the prediction of artificial satellites in the Earth's gravitational field with axial symmetry. Applications of the algorithm for the problem of the final state prediction are illustrated by numerical examples of three test orbits each for two geopotential models corresponding to N = 2 and N = 36.A final state of any desired accuracy is obtained for each case study, a result which shows the flexibility of the algorithm.Dept. of Astronomy, KAU     

Main belt binary asteroidal systems with eccentric mutual orbits

Using 8-10-m class telescopes and their Adaptive Optics (AO) systems, we conducted a long-term adaptive optics campaign initiated in 2003 focusing on four binary asteroid systems: (130) Elektra, (283) Emma, (379) Huenna, and (3749) Balam. The analysis of these data confirms the presence of their asteroidal satellite. We did not detect any additional satellite around these systems even though we have the capability of detecting a loosely-bound fragment (located at 1/4×RHill) ∼40 times smaller in diameter than the primary. The orbits derived for their satellites display significant eccentricity, ranging from 0.1 to 0.9, suggesting a different origin. Based on AO size estimate, we show that (130) Elektra and (283) Emma, G-type and P-type asteroids, respectively, have a significant porosity (30-60% considering CI-CO meteorites as analogs) and their satellite's eccentricities (e∼0.1) are possibly due to excitation by tidal effects. (379) Huenna and (3749) Balam, two loosely bound binary systems, are most likely formed by mutual capture. (3749) Balam's possible high bulk density is similar to (433) Eros, another S-type asteroid, and should be poorly fractured as well. (379) Huenna seems to display both characteristics: the moonlet orbits far away from the primary in term of stability (20%×RHill), but the primary's porosity is significant (30-60%).     

Accurate integration of geostationary orbits with Burdet's focal elements

The theory of Burdet's focal elements is outlined. The differential equations are presented, and the initial value problem is described together with the transformation to rectangular coordinates and classical elements. The focal elements are well defined for zero eccentricity and inclination. They can be adopted for the computation of elliptic, parabolic and hyperbolic motion. For the numerical integration of near-geostationary orbits a comparison of the efficiency is made between focal elements, KS theory and rectangular coordinates. For this class of orbits, a higher accuracy has been obtained by integrating elements than integrating rectangular coordinates.     

Stability of axial orbits in galactic potentials

We investigate the dynamics in a galactic potential with two reflection symmetries. The phase-space structure of the real system is approximated with a resonant detuned normal form constructed with the method based on the Lie transform. Attention is focused on the stability properties of the axial periodic orbits that play an important role in galactic models. Using energy and ellipticity as parameters, we find analytical expressions of bifurcations and compare them with numerical results available in the literature.     

Prediction of trajectories in the Earth's gravitational field with axial symmetry using KS uniform regular elements

In this paper, a motion prediction algorithm based on the KS regular elements is developed for the motion in the Earth's gravitational field with axial symmetry. The algorithm is of recursive nature and general in the sense that it could be applied for any conic motion whatever the number of zonal harmonic coefficientsN 2 may be. Applications of the algorithm for the problem of the final state prediction are illustrated by numerical examples of eight typical ballistic missiles for geopotential model with zonal harmonic terms up toJ ₃₆. A final state of any desired accuracy is obtained for each case study, a result which shows the efficiency and the flexibility of the algorithm.     

On the rotation of co-orbital bodies in eccentric orbits

We investigate the resonant rotation of co-orbital bodies in eccentric and planar orbits. We develop a simple analytical model to study the impact of the eccentricity and orbital perturbations on the spin dynamics. This model is relevant in the entire domain of horseshoe and tadpole orbit, for moderate eccentricities. We show that there are three different families of spin–orbit resonances, one depending on the eccentricity, one depending on the orbital libration frequency, and another depending on the pericenter’s dynamics. We can estimate the width and the location of the different resonant islands in the phase space, predicting which are the more likely to capture the spin of the rotating body. In some regions of the phase space the resonant islands may overlap, giving rise to chaotic rotation.     

Prediction of satellite motions in the Earth's gravitational field with axial symmetry by the KS regularized theory

In this paper, economical and stable recurrence formulae for the Earth's zonal potential and its gradient for the KS regularized theory will be established for any numberN of the zonal harmonic coefficient. A general recursive computational algorithm based on these formulae is also established for the initial value problem of the KS theory for the prediction of artificial satellites in the Earth's gravitational field with axial symmetry. Applications of the algorithm for the problem of the final state prediction are illustrated by numerical examples of three test orbits each for two geopotential models corresponding toN=2 andN=36. A final state of any desired accuracy is obtained for each case study, a result which shows the flexibility of the algorithm.     

Stability of axial orbits in analytical galactic potentials

We study the stability of axial orbits in analytical galactic potentials as a function of the energy of the orbit and the ellipticity of the potential. The problem is solved by an analytical method, the validity of which is not limited to small amplitudes. The lines of neutral stability divide the parameter space in regions corresponding to different organizations of the main families of orbits in the symmetry planes.     

Coordinates for perturbed keplerian systems with axial symmetry

A coordinate system is defined on the phase space of a perturbed Keplerian system after the mean anomaly has been averaged out, for the purpose of explaining how eliminating the longitude of the ascending node reduces the orbital space to a two-dimensional sphere in case the system admits an axial symmetry. Concomitantly, on the submanifold of direct osculating ellipses, the CDM variables are replaced by functions which form the basis of a Poisson algebra isomorphic to the Lie algebra so(3) of the rotation group SO(3); furthermore, in these variables, the doubly reduced phase flow appears like a rotation of the reduced phase space.     

Exact solution to the relative orbital motion in eccentric orbits

This paper studies the relative orbital motion between arbitrary Keplerian trajectories. A closed-form vectorial solution to the nonlinear initial value problem that models this type of motion with respect to a noninertial reference frame is offered. Without imposing any particular conditions on the leader or the deputy satellites trajectories, exact expressions for the relative law of motion and relative velocity are obtained in a closed form. This solution allows the parameterization of the relative motion manifold and offers new methods to study its geometrical and topological properties. The result presented in this paper opens the way to obtain new classes of approximate solutions to the equations of relative motion with time, an eccentric or true anomaly as independent variables. Published in Russian in Solar System Research, 2009, Vol. 43, No. 1, pp. 44–55. The text was submitted by the autors in English.     

Broadband gravitational-wave pulses from binary neutron stars in eccentric orbits

The gravitational-wave radiation from binary stars in elliptical orbits peaks at times close to the periastron passage. For a stationary distribution of binary neutron stars in the Galaxy, there are several systems with large orbital eccentricities and periods in the range from several tens of minutes to several days from which gravitational-wave radiation at periastron will be observed as a broad pulse in the frequency range 1–100 mHz. The LISA space interferometer will be able to record pulsed signals from these systems at a signal-to-noise ratio $S/N > 5\sqrt 5$ in the frequency range ～10^?3–10^?1 Hz. Algorithms for detecting such signals are discussed.     

Decay of a highly eccentric satellite

The re-entry phase of a highly eccentric satellite is discussed. Numerical simulations allowing the prediction of the exact date of re-entry of a highly eccentric satellite are exposed.It is shown that under very particular circumstances the life of the satellite can be extended by a few days. The number of final revolutions of the rapidly contracting orbit depends critically on the air density between 70 km and 100 km.Re-entry of the European scientific satellite HEOS-1 predicted for 28 October, 1975 is near such a situation.