Perturbation theory in rectangular coordinates

For treating the perturbed two-body problem in rectangular coordinates a new method is developed. The method is based on the reduction of the variational equations of the two-body problem with arbitrary elements to the Jordan system. The equations of perturbed motion rewritten in the quasi-Jordan form are subjected to a transformation excluding fast variables and leading to a system governing the long term evolution of motion. The method may be easily extended to the problem of the heliocentric motion of the major planets. For performing this method on computer it is suitable to use facilities of Poissonian and Keplerian processors.     

Perturbations in rectangular coordinates by iteration

An investigation has been made on computing orbits with Picard's method of successive approximations. The perturbations are integrated in the form of a general displacement from a fixed Keplerian reference orbit. Several variation-of-parameters methods are obtained for the integration of the displacement equation. These variation-of-parameters methods could be used as special perturbation or general perturbation methods. The present paper investigates the applications as iterative numerical perturbation techniques. Four different formulations are proposed. They have been implemented on a computer with Chebychev series and their respective advantages and disadvantages are analyzed. Connections with other known perturbation methods are also described.This paper presents the results of one phase of research carried out at the Jet Propulsion Laboratory, California Institute of Technology, under Contract No. NAS 7-100, sponsored by the National Aeronautics and Space Administration.     

Time elements in Keplerian orbital elements

Time elements are introduced in terms of Keplerian (classical) orbital elements for use with time transformations of the Sundman type. Three different time elements are introduced. One time element is associated with the eccentric anomaly, a second time element is associated with the true anomaly, and a third time element is associated with theintermediate anomaly.Numerical results are presented that show accuracy improvements of from one to two orders of magnitude when time elements are employed along with Sundman time transformations, compared with using time transformations alone.     

A set of elements based on polar coordinates in the KS-space and its applications

A new set of canonical elements is introduced into the field of KS-theory. The close relationship of these elements with a set of elements proposed by Scheifele (1970) is analysed. Some applications are outlined.Presented at the Conference on Celestial Mechanics, Oberwolfach, Germany, August 27–September 2, 1972.     

Report of the IAU/IAG Working Group on cartographic coordinates and rotational elements: 2006

Every three years the IAU/IAG Working Group on Cartographic Coordinates and Rotational Elements revises tables giving the directions of the poles of rotation and the prime meridians of the planets, satellites, minor planets, and comets. This report introduces improved values for the pole and rotation rate of Pluto, Charon, and Phoebe, the pole of Jupiter, the sizes and shapes of Saturn satellites and Charon, and the poles, rotation rates, and sizes of some minor planets and comets. A high precision realization for the pole and rotation rate of the Moon is provided. The expression for the Sun’s rotation has been changed to be consistent with the planets and to account for light travel time     

The Roche coordinates in non-synchronous binaries

The Cayley-Darboux problem for the Roche model of binaries is reinvestigated. Generalised Roche coordinates are then defined and calculated in the form of power series of potential for the general case of non-synchronous binaries with eccentric orbits.     

Isovortical orbits in uniformly rotating coordinates

For the conservative, two degree-of-freedom system with autonomous potential functionV(x,y) in rotating coordinates; $$\dot u - 2n\upsilon = V_x , \dot \upsilon + 2nu = V_y $$ , vorticity (v _x -u _y ) is constant along the orbit when the relative velocity field is divergence-free such that: $$u(x,y,t) = \psi _y , \upsilon (x,y,t) = - \psi _x $$ . Unlike isoenergetic reduction using the Jacobi, integral and eliminating the time,non-singular reduction from fourth to second-order occurs when (u,v) are determined explicitly as functions of their arguments by solving for ψ (x, y, t). The orbit function ψ satisfies a second-order, non-linear partial differential equation of the Monge Ampere type: $$2(\psi _{xx} \psi _{yy} - \psi _{xy}^2 ) - 2(\psi _{xx} + \psi _{yy} ) + V_{xx} + V_{yy} = 0$$ . Isovortical orbits in the rotating frame arenot level curves of ψ because it contains time explicitly due to coriolis effects. Rather, (x, y) coordinates along the orbit are obtained, from (u, v) either by numerical integration of the kinematic equations, or by partial differentiation of the Legendre transform ? of ψ. In the latter case, ? is shown to satisfy a non-linear, second-order partial differential equation in three independent variables, derived from the Monge-Ampere Equation. Complete reduction to quadrature is possible when space-time symmetries exist, as in the case of central force motion.     

The celestial pole coordinates

The coordinates of the Celestial Ephemeris Pole in the Celestial Reference System (CRS) can advantageously replace the classical precession and nutation parameters in the matrix transformation of vector components from the CRS to the Terrestrial Reference System (TRS). This paper shows that the new matrix transformation using these coordinates in place of the preceding parameters would be conceptually more simple, especially when associated with the use of the non-rotating origin on the instantaneous equator (Guinot 1979, Capitaine et al. 1986) and of a celestial reference frame as realized by positions of extragalactic sources. In such a representation, the artificial separation between precession and nutation is avoided and the practical computation of the matrix transformation only requires the knowledge of the two celestial direction cosines of the pole, instead of the large number of the quantities generally considered. The development of these coordinates is given as function of time so that their use is equivalent (when using the CRS defined by the mean pole and mean equinox of epoch J2000.0, the 1976 IAU System of Astronomical Constants and the 1980 IAU theory of nutation) to the one of the conventional series for the precession (Lieske et al. 1977) and nutation (Seidelmann 1982) parameters. Such a theoretical development should also be used in order to derive more directly the numerical coefficients of the celestial motion of the instantaneous equator from very precise observations such as VLBI.

---

Résumé Les coordonnées du Pôle Céleste des Ephémerides dans le Systeme de Référence Céleste (CRS) pourraient remplacer avantageusement les paramètres classiques de precession et de nutation dans la matrice de transformation entre le CRS et le Système de Référence Terrestre (TRS). Cet article montre que la nouvelle matrice de transformation utilisant ces coordonnées à la place des paramètres classiques serait ainsi conceptuellement plus simple, en particulier lorsque l'on utilise l'origine non-tournante sur l'équa-teur instantané (Guinot 1979, Capitaine et al. 1986), ainsi que le repère de référence céleste réalisé par les positions des radiosources extragalactiques. Une telle representation évite la séparation artificielle entre précession et nutation et le calcul de la matrice de transformation correspondante ne nécessite que la connaissance des deux cosinus directeurs du pole dans le repère céleste, au lieu du grand nombre de paramètres considérés généralement. Le dèveloppement de ces coordonnées en fonction du temps est donné de façon à ce que leur usage soit équivalent (lorsque l'on se rapporte au CRS défine par le pôle et l'équinoxe moyens de l'époque J2000.0, au Système de Constantes Astronomiques IAU-1976, ainsi qu'au modèle UAI-1980 de la nutation) à celui des séries conventionnelles de la precession (Lieske et al. 1977) et de la nutation (Seidelmann 1982). Un tel développement théorique devrait également être utilise pour déterminer plus directement les coefficients numériques du déplacement céleste de l'équateur instantané, à partir des observations très précises, comme par exemple, les observations VLBI.

     

Ignorable coordinates in the ideal resonance problem

If a dynamical system ofN degrees of freedom is reduced to the Ideal Resonance Problem, the Hamiltonian takes the form $$F = B(y) + 2\mu ^2 A(y)\sin ^2 x_1 , \mu<< 1.$$ Herey is the momentum-vectory _k withk=1, 2,...,N, andx ₁ is thecritical argument. A first-orderglobal solution,x ₁(t) andy ₁(t), for theactive variables of the problem, has been given in Garfinkelet al. (1971). Sincex _k fork>1 are ignorable coordinates, it follows that $$y_\kappa = const., k > 1.$$ The solution is completed here by the construction of the functionsx _k(t) fork>1, derivable from the new HamiltonianF′(y′) and the generatorS(x, y′) of the von Zeipel canonical transformation used in the cited paper. The solution is subject to thenormality condition, derived in a previous paper fork=1, and extended here to 2≤k≤N. It is shown that the condition is satisfied in the problem of the critical inclination provided it is satisfied fork=1.     

Motion on two-dimensional manifolds in rotating coordinates

The two degree-of-freedom system in rotating coordinates: \.u – 2nv = V _x, \.v + 2nu = V _y, \.x = u, \.y = v and its Jacobi integral define a time-dependent velocity field on a differentiable, two-dimensional manifold of integral curves. Explicit time dependence is determined by the dynamical system, coordinate frame, and initial conditions. In the autonomous cases, orbits are level curves of an autonomous function satisfying a second-order, quasi-linear, partial differential equation of parabolic type. Important aspects of the theory are illustrated for the two-body problem in rotating coordinates.     

Micropulsations and orthogonal curvilinear coordinates

The theory of torsional hydromagnetic oscillations of the magnetosphere is usually cast in terms of orthogonal curvilinear coordinates. For a general magnetic field B with potential Ω it is shown that no coordinates exist in which a suitable solution may be found unless the Alfvén velocity V_A, together with B and Ω satisfy certain functional relationships. In the case V_A = constant, for example we must have

$\text{(}\text{B}\text{· ?)B =}\text{function of}\text{B}\text{and}\text{Ω}\text{only}$

. The relationships presented are in fact satisfied by all the magnetic fields considered to date.     

Accurate coordinates of planetary nebulae

Accurate optical coordinates of 734 PNe, measured on the charts of the Digitized Palomar Sky Survey, are presented. As a result of the discussion about the external accuracy the constants –0.8″ in RA and +0.8″ in DEC should be added to the coordinates measured by us. They were used but rounded off already in CGPN(2000). The list and measurements of new 31 candidates of central stars are given which might be interesting for stellar evolution.     

Which electromagnetic equations apply in rotating coordinates?

It was discovered some years ago by Schiff that the equations divE = 4πQ and curlB - (1/c) ∂E/∂t = (4π/c)J for fields in vacuum do not carry over without change from an inertial frame to a frame with rotating axes of space coordinates, even for a region with all velocities of orderv≪c. However, the belief that all four of the field equations are invariant under such conditions is still prevalent and causes misconceptions in physical applications, including astrophysical and geophysical ones. The purpose of the present paper is therefore to call attention to Schiff's discovery, discussing its basis and its extension to fields in material media, and to interpret the additional terms that must be added to the equations in order to obtain valid transformations to rotating axes of coordinates.     

Families of periodic orbits continued in regularizing coordinates

A predictor-corrector algorithm is proposed for continuing analytically families of periodic orbits beyond collision trajectories in the restricted problem of three bodies. It is based on Hill's equation for normal variations in Thiele's regularizing coordinates.     

The circular restricted three-body problem in curvilinear coordinates

We derive the exact equations of motion for the circular restricted three-body problem in cylindrical curvilinear coordinates together with a number of useful analytical relations linking curvilinear coordinates and classical orbital elements. The equations of motion can be seen as a generalization of Hill’s problem after including all neglected nonlinear terms. As an application of the method, we obtain a new expression for the averaged third-body disturbing function including eccentricity and inclination terms. We employ the latter to study the dynamics of the guiding center for the problem of circular coorbital motion providing an extension of some results in the literature.     

Report of the IAU/IAG/COSPAR working group on cartographic coordinates and rotational elements of the planets and satellites: 1991

Every three years the IAU/IAG/COSPAR Working Group on Cartographic Coordinates and Rotational Elements of the Planets and Satellites revises tables giving the directions of the north poles of rotation and the prime meridians of the planets and satellites. Also presented are revised tables giving their sizes and shapes.     

Conversion from geocentric to geodetic coordinates

The geodetic latitude and the height of a satellite are obtained as power series in the ellipsoid's eccentricity, the terms being Fourier sums in the geocentric latitude with polynomials in 1/r as coefficients, with a view of determining the height to the fullest accuracy required by altimetry and geodesy from satellites.     

Matrix perturbation methods using regularized coordinates

Matrix methods for computing perturbations of non-linear perturbed systems, as formulated by Alexeev, involve an expression for the full solution of the first variational equations of the system evaluated about a reference orbit. These cannot be immediately applied to a regularized system of equations where perturbations about Keplerian motion are considered since the solution of the variational equations of regularized Keplerian motion does not in general correspond to the solution of the variational equations of the unregularized equations. But, as Kustaanheimo and Stiefel have pointed out, the regularized equations of Keplerian motion should be excellent for the initiation of a perturbation theory since they are linear in form. This paper describes a method for applying Alexeev's theorem to a regularized system where full advantage is taken of the basic linear form of the unperturbed equations.Presented at the Conference on Celestial Mechanics, Oberwolfach, Germany, August 17–23, 1969.     

Non-singular coordinates of some Kiselev space-times

In this article, non-singular Kruskal-like coordinates of some Kiselev space-times are presented. Also, non-singular Carter-like coordinates are constructed for the extreme case of Kiselev space-time.