^{1}

^{*}

^{2}

^{*}

In this paper, the equations of motion for spatial restricted circular three body problem will be established using the cylindrical coordinates. Initial value procedure that can be used to compute both the cylindrical and Cartesian coordinates and velocities is also developed.

Depending on the application, a certain coordinate system may be simpler to use than the Cartesian coordinate system. As an example, a physical problem with spherical symmetry defined in R^{3} (e.g., motion in the field of a point mass), is usually easier to solve in spherical polar coordinates than in Cartesian coordinates. Also, for instance, in the galactic rotation, cylindrical coordinates are usually adopted, while the spherical coordinates are suitable for the dynamics of globular clusters. In fact, the change of the dependent and/or independent variables for the differential equations of motion becomes of the focal point of researches in space dynamics. Some authors proposed successful methods to change of the dependent and/or independent variables so as to regularize the differential equations of motion. Of these, the method established by Stiefel and Scheifele, in 1971 [

The great success of the these devices in regularizing the equations of motion for the perturbed two body problem ,and on the other hand, the importance of the three body problem in space dynamics (e.g [

The present paper represents the first phase of our studies towards establishing new differential equations for the different forms of the three body problem using some important coordinate systems.

The objective of the present paper, is to establish the equations of motion for spatial restricted circular three body problem in cylindrical coordinates system together with initial value procedure that can be used to compute both the cylindrical and Cartesian coordinates and velocities.

If two of the bodies, say and in the three-body problem move in circular, coplanar orbits about their common center of mass and the mass of the third body is too small to affect the motion of the other bodies, the problem of the motion of the third body is called the circular ,restricted ,three body problem. The two revolving bodies are called the primaries, their masses are arbitrary but have such internal mass distributions that they may be considered point masses. The equations of motion of the third body in a dimensionless synodic rotating coordinate system with the mean motion, are [

where is given as

denotes the mass of the smaller primary when the total mass of the primaries has been normalized to unity.

In what follows we shall establish ,the differential equations for the spatial circular restricted three body problem in the cylindrical coordinate system together with a list of the direct and the, inverse formulae of the transformation.

;; (7)

where

, ,.

From Equation (7) we have

;;. (9)

From Equation (8) we get

;;. (10)

where is given in terms of and y by the first of Equation (9).

The kinetic energy of a particle of unit mass in the cylindrical coordinate system is

By using the transformation equations (Equation (7)), the gravitational potential V could be expressed in term of.

Using Lagrange’s dynamical equations, we have

,

.

Consequently, we deduce for the equations of motion in the cylindrical coordinate the forms

where are given as

can be computed from Equation (7), while and can be computed from Equation (4), and we get

.

where

In what follows, we shall establish a procedure that can be used to compute (say) both:

1) the cylindrical coordinates and velocities

, and

2) the Cartesian coordinates and velocities

.

So, such procedure is a double usefulness computational algorithm, for which a differential solver can be used for the cylindrical six order system to obtain. While the Cartesian coordinates and velocities are obtained by the substitution in the direct transformation formulae (Equations (7) and (8)), rather than solving the six order system of Equations (1), (2) and (3). By this way, great time can be saved.

This initial value procedure using sidereal cylindrical coordinate system will be described through its basic points, input, output and computational steps.

Input: 1) at

2) the final time

3), ,

.

Output: 1) from to.

2) from to.

1) Using the inverse transformation Equations (9) and (10) to find the initial values using the given values at.

2) Using the partial derivatives (functions of) to construct the analytical forms of Equation (12) as first order system in the form:

3) Using the initial conditions from step 1 to solve numerically the above differential system for from to, (note that).

4) Using from step 3 and the direct transformations Equations (7) and (8) to compute numerically and to.

5) End.

In this paper, the equations of motion for spatial restricted circular three body problem was established in cylindrical coordinates system. Initial value procedure that can be used to compute both the cylindrical and Cartesian coordinates and velocities is also developed.