Resonances and bifurcations in axisymmetric scale-free potentials

We investigate an analytical treatment of bifurcations of families of resonant 'thin' tubes in axisymmetric galactic potentials. We verify that the most relevant bifurcations are due to the (1:1) resonance producing the 'inclined' orbits through two different mechanisms: from the disc orbit and from the 'thin' tube associated with the vertical oscillation. The closest resonances occurring after these are the (4:3) resonance in the oblate case and the (2:1) resonance in the prolate case. The (1:1) resonances are treated in a straightforward way using a second-order truncated normal form. The higher order resonances are instead cumbersome to investigate, because the normal form has to be truncated to a high degree and the number of terms grows very rapidly. We therefore adopt a further simplification giving analytical formulae for the values of the parameters at which bifurcations ensue and compare them with selected numerical results. Thanks to the asymptotic nature of the series involved, the predictions are reliable well beyond the convergence radius of the original series.     

Specialized celestial mechanics systems for symbolic manipulation

Construction and application of the current high accuracy analytical theories of motion of celestial bodies necessitates the development of specialized software for the implementation of analytical algorithms of celestial mechanics. This paper describes a typical software package of this kind. This package includes a universal Poisson processor for the rational functions of many variables, a tensorial processor for purposes of relativistic celestial mechanics, a Keplerian processor valid for the solutions of the two body problem in the form of a Poisson series, Taylor expansions in powers of time and closed expressions, and an analytical generator of celestial mechanics functions, facilitating the immediate implementation of the present analytical methods of celestial mechanics. The package is completed with a numerical-analytical interface designed, in particular, for the fast evaluation of the long Poisson series.     

Integration of the elliptic restricted three-body problem with Lie series

The method of Lie series is used to construct a solution for the elliptic restricted three body problem. In a synodic pulsating coordinate system, the Lie operator for the motion of the third infinitesimal body is derived as function of coordinates, velocities and true anomaly of the primaries. The terms of the Lie series for the solution are then calculated with recurrence formulae which enable a rapid successive calculation of any desired number of terms. This procedure gives a very useful analytical form for the series and allows a quick calculation of the orbit.The project is supported by the Austrian Fonds zur Förderung der wissénschaftlichen Forschung under Project No. 4471.     

Analytical theories of satellite motion constructed with a new package of formula manipulation and compared with numerical integration

Specialized to the Lie series based perturbation method of Kirchgraber and Stiefel (1978) a new computer algebra package called ANALYTOS has been developed for constructing analytical orbital theories either in noncanonical or canonical form. We present results on the (extended) Main Problem of orbital theory of artificial earth satellites and related issues. The order of the solutions achieved is generally one order higher than those known from literature. Moreover, the analytical orbits have been checked succesfully against precise numerical ephemerides. This revised version was published online in July 2006 with corrections to the Cover Date.     

Radius of the Roche equipotential surfaces

In literature, there is no exact analytical solution available for determining the radius of Roche equipotential surfaces of distorted close binary systems in synchronous rotation. However, Kopal (Roche Model and Its Application to Close Binary Systems, Advances in Astronomy and Astrophysics, Academic Press, New York 1972) and Morris (Publ. Astron. Soc. Pac. 106:154, 1994) have provided the approximate analytical solutions in the form of infinite mathematical series. These series expressions have been commonly used by various authors to determine the radius of the Roche equipotential surfaces, and hence the equilibrium structures of rotating stars and stars in the binary systems. However, numerical results obtained from these approximating series expressions are not very accurate. In the present paper, we have expanded these series expressions to higher orders so as to improve their accuracy. The objective of this paper is to check, whether, there is any effect on the accuracy of these series expressions when the terms of higher orders are considered. Our results show that in most of the cases these expanded series give better results than the earlier series. We have further used these expanded series to find numerically the volume radius of the Roche equipotential surfaces. The obtained results are in good agreement with the results available in literature. We have also presented simple and accurate approximating formulas to calculate the radius of the primary component in a close binary system. These formulas give very accurate results in a specified range of mass ratio.     

On the artificial satellite theory

Some of the basic ideas of an analytical orbiter theory which is being developed by Hubert Claes in Namur are presented.The theory is based on the Lie transform technique and will be expressed in a closed form up to second order. The inclusion of additional terms of the third order (expanded in power series of the eccentricity) will be considered.Special attention is being given to the choice of the elements and to the final form of the theory. Three main criteria are used. The removal of the virtual singularities of small inclination and eccentricity. The simplicity of the final form of the theory once the elements have been given their numerical values. The numerical stability of the evaluation of the theory.     

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.     

Normal forms for the epicyclic approximations of the perturbed Kepler problem

We compute the normal forms for the Hamiltonian leading to the epicyclic approximations of the (perturbed) Kepler problem in the plane. The Hamiltonian setting corresponds to the dynamics in the Hill synodic system where, by means of the tidal expansion of the potential, the equations of motion take the form of perturbed harmonic oscillators in a rotating frame. In the unperturbed, purely Keplerian case, the post-epicyclic solutions produced with the normal form coincide with those obtained from the expansion of the solution of the Kepler equation. In all cases where the perturbed problem can be cast in autonomous form, the solution is easily obtained as a perturbation series. The generalization to the spatial problem and/or the non-autonomous case is straightforward.     

On the expansion of the secular part of the perturbing function of mutual attraction in the satellite systems of planets

We propose a special representation for the secular part of the perturbing function describing the mutual attraction of satellites. In contrast to the known representations, it has a single analytical form for any ratio between the semimajor axes of the perturbed and perturbing satellites. The resulting expression is a partial sum of a power series with respect to the small eccentricities and planet-equatorial inclinations of the satellites’ orbits. This sum includes terms up to and including the fourth degree with respect to these small parameters. The proposed expansion is compared with one of the known expansions for the secular part of the perturbing function.     

YORP torque as the function of shape harmonics

The second-order analytical approximation of the mean Yarkovsky–O'Keefe–Radzievskii–Paddack (YORP) torque components is given as an explicit function of the shape spherical harmonics coefficients for a sufficiently regular minor body. The results are based upon a new expression for the insolation function, significantly simpler than in previous works. Linearized plane-parallel model of the temperature distribution derived from the insolation function allows us to take into account a non-zero conductivity. Final expressions for the three average components of the YORP torque related with rotation period, obliquity and precession are given in a form of the Legendre series of the cosine of obliquity. The series have good numerical properties and can be easily truncated according to the degree of the Legendre polynomials or associated functions, with first two terms playing the principal role.     

A New Echeloned Poisson Series Processor (EPSP)

A specialized Echeloned Poisson Series Processor (EPSP) is proposed. It is a typical software for the implementation of analytical algorithms of Celestial Mechanics. EPSP is designed for manipulating long polynomial-trigonometric series with literal divisors. The coefficients of these echeloned series are the rational or floating-point numbers. The Keplerian processor and analytical generator of special celestial mechanics functions based on the EPSP are also developed.     

The long period behavior of the orbits of the Galilean satellites of Jupiter

Some of the results of an investigation into the long period behavior of the orbits of the Galilean satellites of Jupiter are presented. Special purpose computer programs were used to perform all the algebraic manipulations and series expansions that are necessary to describe the mutual interactions among the satellites.The disturbing function was expanded as a Poisson series in the modified Keplerian elements referred to a Jovicentric coordinate system. The differential equations for the modified Keplerian elements were then formed, and all short period perturbations were removed using Kamel's perturbation method. Approximate analytical solutions for these differential equations are derived, and the general form of the solutions are given.     

Dynamical problems of thermoelastic solids

The dynamical problem for thermal stresses in an infinite isotropic elastic cylinder of radius a with its axis along the z-axis, subject to fixed boundary conditions is studied. The Fourier heat conduction equation has been solved applying the Fourier transform and the theory of complex variable. The thermoelastic equation of motion has been separated into two wave equations which can be solved separately. The temperature, the displacement and the stress components have been obtained in analytical form as series involving Bessel function of first kind and of order zero.     

Generation of a secular system in the analytical theory for the motion of the moon

In order to generate an analytical theory of the motion of the Moon by considering planetary perturbations, a procedure of general planetary theory (GPT) is used. In this case, the Moon is considered as an addition planet to the eight principal planets. Therefore, according to the GPT procedure, the theory of the Moon’s orbital motion can be presented in the form of series with respect to the evolution of eccentric and oblique variables with quasi-periodic coefficients, which are the functions of mean longitudes for principal planets and the Moon. The relationship between evolution variables and the time is determined by a trigonometric solution for the independent secular system that describes the secular motion of a perigee and the Moon node by considering secular planetary inequalities. Principal planetary coordinates required for generating the theory of the motion of the Moon includes only Keplerian terms, the intermediate orbit, and the linear theory with respect to eccentricities and inclinations in the first order relative to the masses. All analytical calculations are performed by means of the specialized echeloned Poisson Series Processor EPSP.     

Bias in velocity determinations from full-disk solar observations

Estimates for the global solar surface velocity field can be obtained from time series of full-disk solar images. However, bias is introduced by mistakes in the assumptions about the geometry of the observations, and by imperfections in the optical system. I investigate many sources of bias and determine their first-order influence on the measured velocity field (both the transverse and longitudinal components). Results are presented in analytical and pictorial form. By comparing bias velocities of unknown origin with the results, one may obtain clues to the cause of the bias.     

Analytical extension of lunar libration tables

Tables of lunar physical libration defining the analytical dependence upon the parameters of the lunar gravitational field are presented. The tables are obtained on the framework of the main problem in lunar libration by integration of the Hamilton equations reduced to the harmonic oscillator equations.The variables of physical libration have been obtained in the form of Poisson series. The distinguishing feature of the tables is that these series are the analytical extension of semianalytical solution computed for a number of dynamical parameters LURE2.A comparison with the Eckhardt's solution is briefly presented. The previously revealed disagreement of the mean inclination of lunar equator to ecliptic with that in Eckhardt's solution 500 has been maintained.     

Dynamo mechanism for turbulent wave in celestial bodies

It is well known that under cosmic conditions the various modes of plasma turbulence waves (including MHD waves) are easily excited. In this paper we are trying to show that the turbulent wave also generates a source-term for the magnetic induced equations as does the turbulent fluid with nonzero helicity. By expanding the turbulent field in Fourier series, we have obtained dynamo equation for turbulent wave and a reasonable solution which indicates that the poloidal field may be built-up in the turbulent source region. Perhaps, we may think that the poloidal field of Equation (9) is the analytical form of the magnetic field in a turbulent source region of celestial bodies.     

Accurate analytical representation of Pluto modern ephemeris

An accurate development of the latest JPL’s numerical ephemeris of Pluto, DE421, to compact analytical series is done. Rectangular barycentric ICRF coordinates of Pluto from DE421 are approximated by compact Fourier series with a maximum error of 1.3 km over 1900–2050 (the entire time interval covered by the ephemeris). To calculate Pluto positions relative to the Sun, a development of rectangular heliocentric ICRF coordinates of the Solar System barycenter to Poisson series is additionally made. As a result, DE421 Pluto heliocentric positions by the new analytical series are represented to an accuracy of better than 5 km over 1900–2050.     

Accurate analytical nutation series

In this letter we present the first accurate analytical nutation series, deduced from the Hamiltonian theory by the authors. They provide the highest accuracy ever obtained by any analytical nutation series, since the deviation in CEP (celestial ephemeris pole) offsets with respect to that of the IERS Conventions 1996 is kept below 1 mas in the time domain, in spite of still lacking ocanic corrections.     

Symbolic analytical developments of the zero pressure cosmological model of the universe

In the present paper, literal analytical solutions in power series forms are developed for the radius of curvature and the expansion velocity of the zero pressure cosmological models of the universe at any time t. Also, we develop literal analytical solutions in power series forms for the inverse problem of the zero pressure cosmological model, that is to find the time $t=\tilde{t}$ (say) at which the radius of curvature of the model $R=\tilde{R}$ (say) is known. The importance of these analytical power series representations is that, they are invariant under many operations because, addition, multiplication, exponent ion, integration, different ion, etc of a power series is also a power series. A fact which provides excellent flexibility in dealing with analytical as well as computational developments of the problems related to zero pressure cosmological models.For computational developments of these solutions, an efficient method using continued fraction theory is provided. By means of the present methods we able to analyze some known zero-pressure cosmological models, of these are Einstein and De Sitter models. In addition we also analyzed some other models by which one can know if the universe keep expanding forever, or will it reach a maximal size and then turn into contraction stage.