Analytical method for perturbed frozen orbit around an Asteroid in highly inhomogeneous gravitational fields : A first approach

This article provides a method for nding initial conditions for perturbed frozen orbits around inhomogeneous fast rotating asteroids. These orbits can be used as reference trajectories in missions that require close inspection of any rigid body. The generalized perturbative procedure followed exploits the analytical methods of relegation of the argument of node and Delaunay normalisation to arbitrary order. These analytical methods are extremely powerful but highly computational. The gravitational potential of the heterogeneous body is rstly stated, in polar-nodal coordinates, which takes into account the coecients of the spherical harmonics up to an arbitrary order. Through the relegation of the argument of node and the Delaunay normalization, a series of canonical transformations of coordinates is found, which reduces the Hamiltonian describing the system to a integrable, two degrees of freedom Hamiltonian plus a truncated reminder of higher order. Setting eccentricity, argument of pericenter and inclination of the orbit of the truncated system to be constant, initial conditions are found, which evolve into frozen orbits for the truncated system. Using the same initial conditions yields perturbed frozen orbits for the full system, whose perturbation decreases with the consideration of arbitrary homologic equations in the relegation and normalization procedures. Such procedure can be automated for the first homologic equation up to the consideration of any arbitrary number of spherical harmonics coefficients. The project has been developed in collaboration with the European Space Agency (ESA)

Estabilidad orbital de satélites estacionarios

Se establecen las condiciones, para un satélite alrededor de un planeta en rotación alrededor de su eje de inercia de mayor mo inento, bajo las que puede aparecer, en un sistema de referencia fijo en el planeta, dos posiciones de equilibrio con exponentes carac terísticos imaginarios puros. En este caso, despu& de la apropiada normalización mediante una transformación de Lic, efectuada de forma mecánica con un procesador algebraico simbólico, se aplica el teorema de Arnold sobre formas cuadráticas no definidas. Se concluye que los equilibrios son estables en el sentido de Liapunov. Las condiciones de estabilidad se verifican en el caso de la Tierra.. For a satellite about an oblate planet in rotation about its axis of greatest inertia, conditions are given under which there may appear, in a frame fixed in the planet, two positions of equilib ría with characteristic exponents that are purely imaginary. In which case, after appropriate nonnalization by Lic transformation executed mechanically through a symbolic algebraic processor, the theorem of Amold about non definite quadratie forms is applied. It is concluded that the equilibria are stable in the sense of Lia punov. The conditions for stability are verified in the case of tbe earth. 1991 Mathematics Subject Classiflcation 68Q40, 70-04, 70F15, 701115, 70325 1985 A. C. M. Classiñcation: 1.1.4 Servicio Publicaciones Univ. Complutense. Madrid, 1996. *flnwciado parcialmente por el Ministerio Español de Educación y Ciencia, Proyecto DGICYT # PB93-1236-C02-02. 312 A. Depñt y It López Moratail

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 normalised the Hamiltonian using Lie transform as in Coppola and Rand (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_6$ is stable. We plotted graphs for $(\omega_1, D_2).$ They are rectangular hyperbola.Comment: Accepted for publication in Astrophysics & Space Scienc

Production of trans-Neptunian binaries through chaos-assisted capture

The recent discovery of binary objects in the Kuiper-belt opens an invaluable window into past and present conditions in the trans-Neptunian part of the Solar System. For example, knowledge of how these objects formed can be used to impose constraints on planetary formation theories. We have recently proposed a binary-object formation model based on the notion of chaos-assisted capture. Here we present a more detailed analysis with calculations performed in the spatial (three-dimensional) three- and four-body Hill approximations. It is assumed that the potential binary partners are initially following heliocentric Keplerian orbits and that their relative motion becomes perturbed as these objects undergo close encounters. First, the mass, velocity, and orbital element distribu- tions which favour binary formation are identified in the circular and elliptical Hill limits. We then consider intruder scattering in the circular Hill four-body problem and find that the chaos-assisted capture mechanism is consistent with observed, apparently randomly distributed, binary mutual orbit inclinations. It also predicts asymmetric distributions of retrograde versus prograde orbits. The time-delay induced by chaos on particle transport through the Hill sphere is analogous to the formation of a resonance in a chemical reaction. Implications for binary formation rates are considered and the 'fine-tuning' problem recently identified by Noll et al. (2007) is also addressed.Comment: submitted to MNRA

Accurate free and forced rotational motions of rigid Venus

% context :The precise and accurate modelling of a terrestrial planet like Venus is an exciting and challenging topic, all the more interesting since it can be compared with that of the Earth for which such a modelling has already been achieved at the milliarcsecond level % aims: We want to complete a previous study (Cottereau and Souchay, 2009), by determining at the milliarcsecond level the polhody, i.e. the torque-free motion of the axis of angular momentum of a rigid Venus in a body-fixed frame, as well as the nutation of its third axis of figure in space, which is fundamental from an observational point of view. results :In a first part we have computed the polhody, i.e. the respective free rotational motion of the axis of angular momentum of Venus with respect to a body-fixed frame. We have shown that this motion is highly elliptical, with a very long period of 525 cy to be compared with 430 d for the Earth. This is due to the very small dynamical flattening of Venus in comparison with our planet. In a second part we have computed precisely the Oppolzer terms which allow to represent the motion in space of the third Venus figure axis with respect to Venus angular momentum axis, under the influence of the solar gravitational torque. We have determined the corresponding tables of coefficients of nutation of the third figure axis both in longitude and in obliquity due to the Sun, which are of the same order of amplitude as for the Earth. We have shown that the coefficients of nutation for the third figure axis are significantly different from those of the angular momentum axis on the contrary of the Earth. Our analytical results have been validated by a numerical integration which revealed the indirect planetary effects.Comment: 14 pages, 11 figures, accepted for publication in section 11. Celestial mechanics and astrometry of Astronomy and Astrophysics (27/02/2010

Families of Canonical Transformations by Hamilton-Jacobi-Poincar\'e equation. Application to Rotational and Orbital Motion

The Hamilton-Jacobi equation in the sense of Poincar\'e, i.e. formulated in the extended phase space and including regularization, is revisited building canonical transformations with the purpose of Hamiltonian reduction. We illustrate our approach dealing with orbital and attitude dynamics. Based on the use of Whittaker and Andoyer symplectic charts, for which all but one coordinates are cyclic in the Hamilton-Jacobi equation, we provide whole families of canonical transformations, among which one recognizes the familiar ones used in orbital and attitude dynamics. In addition, new canonical transformations are demonstrated.Comment: 21 page

Stress field and spin axis relaxation for inelastic triaxial ellipsoids

A compact formula for the stress tensor inside a self-gravitating, triaxial ellipsoid in an arbitrary rotation state is given. It contains no singularity in the incompressible medium limit. The stress tensor and the quality factor model are used to derive a solution for the energy dissipation resulting in the damping (short axis mode) or excitation (long axis) of wobbling. In the limit of an ellipsoid of revolution, we compare our solution with earlier ones and show that, with appropriate corrections, the differences in damping times estimates are much smaller than it has been claimed. This version implements corrections of misprints found in the MNRAS published text.Comment: 14 pages, 6 figures, published in Monthly Notices RAS (containing misprints

Resonances in a spring-pendulum: algorithms for equivariant singularity theory

A spring-pendulum in resonance is a time-independent Hamiltonian model system for formal reduction to one degree of freedom, where some symmetry (reversibility) is maintained. The reduction is handled by equivariant singularity theory with a distinguished parameter, yielding an integrable approximation of the Poincaré map. This makes a concise description of certain bifurcations possible. The computation of reparametrizations from normal form to the actual system is performed by Gröbner basis techniques.

The long-term stability of extrasolar system HD 37124. Numerical study of resonance effects

We describe numerical tools for the stability analysis of extrasolar planetary systems. In particular, we consider the relative Poincare variables and symplectic integration of the equations of motion. We apply the tangent map to derive a numerically efficient algorithm of the fast indicator MEGNO (a measure of the maximal Lyapunov exponent) that helps to distinguish chaotic and regular configurations. The results concerning the three-planet extrasolar system HD 37124 are presented and discussed. The best fit solutions found in earlier works are studied more closely. The system involves Jovian planets with similar masses. The orbits have moderate eccentricities, nevertheless the best fit solutions are found in dynamically active region of the phase space. The long term stability of the system is determined by a net of low-order two-body and three-body mean motion resonances. In particular, the three-body resonances may induce strong chaos that leads to self-destruction of the system after Myrs of apparently stable and bounded evolution. In such a case, numerically efficient dynamical maps are useful to resolve the fine structure of the phase space and to identify the sources of unstable behavior.Comment: 11 pages (total), 8 figures. Accepted for publication in MNRAS. The definitive version will be/is available at http://www.blackwellpublishing.com. The astro-ph version is prepared with low resolution figures. To obtain the manuscript with full-resolution figures, please visit http://www.astri.uni.torun.pl/~chris/mnrasIII.ps.g

Secondary resonances of co-orbital motions

The size distribution of the stability region around the Lagrangian point L4 is investigated in the elliptic restricted three-body problem as the function of the mass parameter and the orbital eccentricity of the primaries. It is shown that there are minimum zones in the size distribution of the stability regions, and these zones are connected with secondary resonances between the frequencies of librational motions around L4. The results can be applied to hypothetical Trojan planets for predicting values of the mass parameter and the eccentricity for which such objects can be expected or their existence is less probable.Comment: 9 pages, 7 figures, accepted for publication in MNRA