On an unharmonic oscillator potential for elliptical galaxies

Using a direct numerical procedure we approximate a ’realistic’ potential for an elliptical galaxy by that of a perturbed harmonic oscillator. The quality of fitting is checked using two criteria. First by computing the value of the fitting parameterf, and second by comparing the behaviour of orbits for the two potentials. Both criteria suggest that the fitting is good when the total energyh is smaller than the energy of escape in the unharmonic potential.     

Comparison of Numerical Methods for the Integration of Natural Satellite Systems

We present a new implementation of the recurrent power series (RPS) method which we have developed for the integration of the system of N satellites orbiting a point-mass planet. This implementation is proved to be more efficient than previously developed implementations of the same method. Furthermore, its comparison with two of the most popular numerical integration methods: the 10th-order Gauss–Jackson backward difference method and the Runge–Kutta–NystrRKN12(10)17M shows that the RPS method is more than one order of magnitude better in accuracy than the other two. Various test problems with one up to four satellites are used, with initial conditions obtained from ephemerides of the saturnian satellite system. For each of the three methods we find the values of the user-specified parameters (such as the method's step-size (h or tolerance (TOL)) that minimize the global error in the satellites' coordinates while keeping the computer time within reasonable limits. While the optimal values of the step-sizes for the methods GJ and RKN are all very small (less than T/100, the ones that are suitable for the RPS method are within the range: T/13<h<T/6 (T being the period of the innermost satellite of the problem). Comparing the results obtained by the three methods for these step-sizes and for the various test problems we observe the superiority of the RPS method over GJ in terms of accuracy and over RKN both in accuracy and in speed. This revised version was published online in July 2006 with corrections to the Cover Date.     

Numerical investigation of collision orbits of lunar satellites

In the present study an investigation of the collision orbits of natural satellites of the Moon (considered to be of finite dimensions) is developed, and the tendency of natural satellites of the Moon to collide on the visible or the far side of the Moon is studied. The collision course of the satellite is studied up to its impact on the lunar surface for perturbations of its initial orbit arbitrarily induced, for example, by the explosion of a meteorite. Several initial conditions regarding the position of the satellite to collide with the Moon on its near (visible) or far (invisible) side is examined in connection to the initial conditions and the direction of the motion of the satellite. The distribution of the lunar craters-originating impact of lunar satellites or celestial bodies which followed a course around the Moon and lost their stability - is examined. First, we consider the planar motion of the natural satellite and its collision on the Moon's surface without the presence of the Earth and Sun. The initial velocities of the satellite are determined in such a way so its impact on the lunar surface takes place on the visible side of the Moon. Then, we continue imparting these velocities to the satellite, but now in the presence of the Earth and Sun; and study the forementioned impacts of the satellites but now in the Earth-Moon-Satellite system influenced also by the Sun. The initial distances of the satellite are taken as the distances which have been used to compute periodic orbits in the planar restricted three-body problem (cf. Gousidou-Koutita, 1980) and its direction takes different angles with the x-axis (Earth-Moon axis). Finally, we summarise the tendency of the satellite's impact on the visible or invisible side of the Moon.     

Numerical models for the study of motion of lunar satellites

The present study deals with numerical modeling of the elliptic restricted three-body problem as well as of the perturbed elliptic restricted three-body (Earth-Moon-Satellite) problem by a fourth body (Sun). Two numerical algorithms are established and investigated. The first is based on the method of the series solution of the differential equations and the second is based on a 5th-order Runge-Kutta method. The applications concern the solution of the equations and integrals of motion of the circular and elliptical restricted three-body problem as well as the search for periodic orbits of the natural satellites of the Moon in the Earth-Moon system in both cases in which the Moon describes circular or elliptical orbit around the Earth before the perturbations induced by the Sun. After the introduction of the perturbations in the Earth-Moon-Satellite system the motions of the Moon and the Satellite are studied with the same initial conditions which give periodic orbits for the unperturbed elliptic problem.     

A study of collision orbits with the Moon by the method of surface of section

Non-periodic orbits of a natural satellite of the Moon are studied, for the case of the circular three-body problem with the method of surface of section. According to this method, each orbit is represented by a point, in the plane x0\.x, which corresponds to y = 0 and \.y > 0 and a fixed energy. Conclusions are deduced from the shape of this curve for probable collisions of the satellite on the lunar surface. This method of surface of section can be used for the study of orbits which collide with the Moon's surface after a large number of revolutions around the Moon and their study would be difficult to explore with other methods.