Out-of-plane librations of spinning tethered satellite systems

We analyze the out-of-plane librations of a tethered satellite system that is nominally rotating in the orbit plane. To isolate the librational dynamics, the system is modeled as two point masses connected by a rigid rod with the system mass center constrained to an unperturbed circular orbit. For small out-of-plane librations, the in-plane motion is unaffected by the out-of-plane librations and a solution for the in-plane motion is determined in terms of Jacobi elliptic functions. This solution is used in the linearized equation for the out-of-plane librations, resulting in a Hill’s equation. Floquet theory is used to analyze the Hill’s equation, and we show that the out-of-plane librations are unstable for certain ranges of in-plane spin rate. For relatively high in-plane spin rates, the out-of-plane librations are stable, and the Hill’s equation can be approximated by a Mathieu’s equation. Approximate solutions to the Mathieu’s equation are determined, and we analyze the dominant characteristics of the out-of-plane librations for high in-plane spin rates. The results obtained from the analysis of the linearized equations of motion are compared to numerical simulations of the nonlinear equations of motion, as well as numerical simulations of a more realistic system model that accounts for tether flexibility. The instabilities discovered from the linear analysis are present in both the nonlinear system and the more realistic system model. The approximate solutions for the out-of-plane librations compare well to the nonlinear system for relatively high in-plane rotation rates, and also capture the significant qualitative behavior of the flexible system.     

Integration theory for the restricted problem of three axisymmetric bodies

In this paper the circular planar restricted problem of three axisymmetric ellipsoids S _i (i = 1, 2, 3), such that their equatorial planes coincide with the orbital plane of the three centres of masses, be considered. The equations of motion of infinitesimal body S ₃ be obtained in the polar coordinates. Using iteration approach we have given an approximation for another integral, which independent of the Jacobian integral, in the case of P-type orbits (near circular orbits surrounding both primaries).     

Mass effects in the problem of three bodies

This is a study of the dynamical behavior of three point masses moving under their mutual gravitational attraction in a plane. The initial positions and velocities are identical for all cases studied and only the masses of the participating bodies change in the series of numerical experiments. In this way the effect of the coupling terms in the differential equations of motion are investigated. The motion in all 125 cases begins with an interplay between the three bodies, followed by temporary ejections or by an eventual escape. The total mass of the system is kept constant while the massratios change from 1 to 5. The initial velocities being zero, the total energy is negative in all cases.Approximately 74% of the cases disintegrated (i.e. two bodies formed a binary and the third body escaped) in less than 140 time units, 47% in less than 50 time units and 10% ended in escape in less than 10 time units. Considering three stars with total mass 12M , initially placed at 3, 4 and 5 parsec distances (or three galaxies with mass 2.4×10¹² M , initially placed 30, 40 and 50 kpc apart), the unit of time (approximately the crossing time) becomes 1.5×10⁷ y (3.2×10⁷ y). The average time of disintegration was found to be of the order of 10⁹ y. The average semi-major axis of the binaries left behind after disintegration was 0.7 parsec and the average value of the eccentricity was 0.76. The effect of the masses on the escapes was established and it was found that the bodynot with the smallest mass escaped in 13% of the disintegrated cases. The cases which did not disintegrate in 150 time units were analyzed in detail and the time of their eventual escape was estimated.The numerical results are tabulated regarding escape time, ejection period, total energy, escape energy, terminal velocity, semi-major axis, and eccentricity.The evolution of triple systems is followed from interplays through ejections to escapes and the orbital parameters for the separation of these classes are estimated.     

The onset of chaotic motion in the restricted problem of three bodies

A full characterization of a nonintegrable dynamical system requires an investigation into the chaotic properties of that system. One such system, the restricted problem of three bodies, has been studied for over two centuries, yet few studies have examined the chaotic nature of some ot its trajectories. This paper examines and classifies the onset of chaotic motion in the restricted three-body problem through the use of Poincaré surfaces of section, Liapunov characteristic numbers, power spectral density analysis and a newly developed technique called numerical irreversibility. The chaotic motion is found to be intermittent and becomes first evident when the Jacobian constant is slightly higher thanC ₂.     

Transfers between libration-point orbits in the elliptic restricted problem

A strategy is formulated to design optimal time-fixed impulsive transfers between three-dimensional libration-point orbits in the vicinity of the interiorL ₁ libration point of the Sun-Earth/Moon barycenter system. The adjoint equation in terms of rotating coordinates in the elliptic restricted three-body problem is shown to be of a distinctly different form from that obtained in the analysis of trajectories in the two-body problem. Also, the necessary conditions for a time-fixed two-impulse transfer to be optimal are stated in terms of the primer vector. Primer vector theory is then extended to non-optimal impulsive trajectories in order to establish a criterion whereby the addition of an interior impulse reduces total fuel expenditure. The necessary conditions for the local optimality of a transfer containing additional impulses are satisfied by requiring continuity of the Hamiltonian and the derivative of the primer vector at all interior impulses. Determination of the location, orientation, and magnitude of each additional impulse is accomplished by the unconstrained minimization of the cost function using a multivariable search method. Results indicate that substantial savings in fuel can be achieved by the addition of interior impulsive maneuvers on transfers between libration-point orbits.An earlier version was presented as Paper AAS 92–126 at the AAS/AIAA Spaceflight Mechanics Meeting, Colorado Springs, Colorado, February 24–26, 1992.     

Motion around the out-of-plane equilibrium points of the perturbed restricted three-body problem

This paper studies the motion of an infinitesimal body near the out-of-plane equilibrium points, L _6,7, in the perturbed restricted three-body problem. The problem is perturbed in the sense that the primaries of the system are oblate spheroids as well as sources of radiation and small perturbations are give to the Coriolis and centrifugal forces. It locates the positions and examines the stability of L _6,7 with a particular application to the binary system Struve 2398. It is observed that their positions are affected by the radiation, oblateness and a small perturbation in the centrifugal force, but is unaffected by that of the Coriolis force. They are also found to be unstable.     

The variational equation of the three-body problem

Letx ₀ (t),x ₀ ⁴ be a homothetic solution of the planar three-body problem with total energyh, described in relative coordinates with respect to one body. It is shown that the variational equation of the problem atx ₀ (t) can be solved explicitly in terms of hypergeometric functions. This is done by using the scaled true anomaly of the one-dimensional Kepler motion as the independent variable.The classical theorems about hypergeometric functions allow a simple calculation of all the values needed in applications. By means of this theory the past of a homothetic triple close encounter may be described in a linearized approximation.Proceedings of the Sixth Conference on Mathematical Methods in Celestial Mechanics held at Oberwolfach (West Germany) from 14 to 19 August, 1978.     

A new analytic approach to the Sitnikov problem

A new analytic approach to the solution of the Sitnikov Problem is introduced. It is valid for bounded small amplitude solutions (z _max = 0.20) (in dimensionless variables) and eccentricities of the primary bodies in the interval (–0.4 < e < 0.4). First solutions are searched for the limiting case of very small amplitudes for which it is possible to linearize the problem. The solution for this linear equation with a time dependent periodic coefficient is written up to the third order in the primaries eccentricity. After that the lowest order nonlinear amplitude contribution (being of order z ³) is dealt with as perturbation to the linear solution. We first introduce a transformation which reduces the linear part to a harmonic oscillator type equation. Then two near integrals for the nonlinear problem are derived in action angle notation and an analytic expression for the solution z(t) is derived from them. The so found analytic solution is compared to results obtained from numeric integration of the exact equation of motion and is found to be in very good agreement. CERN SL/AP     

On the solution of the problem of three fixed centres

For a special choice of the initial conditions a solution of the plane restricted problem of four bodies i.e. the problem of the motion of a passively gravitating material pointP attracted according to Newton's law by three fixed point massesP ₁,P ₂ andP ₃ has been obtained.     

The global flow of the parabolic restricted three-body problem

We have two mass points of equal masses m ₁=m ₂ > 0 moving under Newton’s law of attraction in a non-collision parabolic orbit while their center of mass is at rest. We consider a third mass point, of mass m ₃=0, moving on the straight line L perpendicular to the plane of motion of the first two mass points and passing through their center of mass. Since m ₃=0, the motion of m ₁ and m ₂ is not affected by the third and from the symmetry of the motion it is clear that m ₃ will remain on the line L. The parabolic restricted three-body problem describes the motion of m ₃. Our main result is the characterization of the global flow of this problem.     

Periodic orbits under combined effects of oblateness and radiation in the restricted problem of three bodies

In this paper, we study the existence of libration points and their linear stability when the three participating bodies are axisymmetric and the primaries are radiating, we found that the collinear points remain unstable, it is further seen that the triangular points are stable for 0<μ<μ _c , and unstable for where , it is also observed that for these points the range of stability will decrease. In addition to this we have studied periodic orbits around these points in the range 0<μ<μ _c , we found that these orbits are elliptical; the frequencies of long and short orbits of the periodic motion are affected by the terms which involve parameters that characterize the oblateness and radiation repulsive forces. The implication is that the period of long periodic orbits adjusts with the change in its frequency while the period of short periodic orbit will decrease.     

Generalized Hill's problem: Some cases of complete integrability

In this paper, we are investigating cases of integrability in the planar Hill's problem. The external potential U _extis supposed to be time independent in a given uniformly rotating frame. Cases of integrability of the relative motion of two interacting particles in the vicinity of an equilibrium solution of U _extare found. In all these cases, the form of the second integral is explicitly given, the first being the Jacobian one. Cases in which the interacting potential U between the two particles is of newtonian type are particularized.     

Effects of zonal harmonics on the out-of-plane equilibrium points in the generalized Robe's circular restricted three-body problem

This article examines the effects of the zonal harmonics on the out-of-plane equilibrium points of Robe's circular restricted three-body problem when the hydrostatic equilibrium shape of the first primary is an oblate spheroid, the shape of the second primary is an oblate spheroid with oblateness coefficients up to the second zonal harmonic, and the full buoyancy of the fluid is considered. It is observed that the size of the oblateness and the zonal harmonics affect the positions of the out-of-plane equilibrium points L₆ and L₇. It is also observed that these points within the possible region of motion are unstable.     

Out-of-plane equilibrium points and stability in the generalised photogravitational restricted three body problem

The restricted three body problem is generalised to include the effects of an inverse square distance radiation pressure force on the infinitesimal mass due to the primaries, which are both radiating. In this paper we investigate the stability of out-of-plane equilibrium points, based on equations in variations. We have found the characteristic equation for the complex normal frequencies which is a sixth order polynomial. Thus we conclude that out-of-plane equilibrium points are unstable due to positive real part in complex roots.     

Numerical explorations of the rectilinear problem of three bodies

We present some results of a numerical exploration of the rectilinear problem of three bodies, with the two outer masses equal. The equations of motion are first given in relative coordinates and in regularized variables, removing both binary collision singularities in a single coordinate transformation. Among our most important results are seven periodic solutions and three symmetric triple collision solutions. Two of these periodic solutions have been continued into families, the outer massm ₃ being the family parameter. One of these families exists for all masses while the second family is a branch of the first at a second-kind critical orbit. This last family ends in a triple collision orbit.Proceedings of the Sixth Conference on Mathematical Methods in Celestial Mechanics held at Oberwolfach (West Germany) from 14 to 19 August, 1978.     

The restricted planar isosceles three-body problem with non-negative energy

We consider a restricted three-body problem consisting of two positive equal masses m ₁ = m ₂ moving, under the mutual gravitational attraction, in a collision orbit and a third infinitesimal mass m ₃ moving in the plane P perpendicular to the line joining m ₁ and m ₂. The plane P is assumed to pass through the center of mass of m ₁ and m ₂. Since the motion of m ₁ and m ₂ is not affected by m ₃, from the symmetry of the configuration it is clear that m ₃ remains in the plane P and the three masses are at the vertices of an isosceles triangle for all time. The restricted planar isosceles three-body problem describes the motion of m ₃ when its angular momentum is different from zero and the motion of m ₁ and m ₂ is not periodic. Our main result is the characterization of the global flow of this problem.     

The uniform,regular differential equations of the KS transformed perturbed two-body problem

The Newtonian differential equations of motion for the two-body problem can be transformed into four, linear, harmonic oscillator equations by simultaneously applying the regularizing time transformation dt/ds=r and the Kustaanheimo-Stiefel (KS) coordinate transformation. The time transformation changes the independent variable from time to a new variables, and the KS transformation transforms the position and velocity vectors from Cartesian space into a four-dimensional space. This paper presents the derivation of uniform, regular equations for the perturbed twobody problem in the four-dimensional space. The variation of parameters technique is used to develop expressions for the derivatives of ten elements (which are constants in the unperturbed motion) for the general case that includes both perturbations which can arise from a potential and perturbations which cannot be derived from a potential. These element differential equations are slightly modified by introducing two additional elements for the time to further improve long term stability of numerical integration.Originally presented at the AAS/AIAA Astrodynamics Specialists Conference, Vail, Colorado, July 1973     

The problem of motions up to the second post-newtonian approximation in general relativity theory

In this paper we give the Hamiltonian function for aN-body system up to the 2-P.N.A. Then as an example, from the LagrangianL _m of a test particle we derive the equations of its motion up to the 2-P.N.A. in the field of a heavy bodym ₂at rest.     

The extended phase space formulation of the Vinti problem

The Vinti problem, motion about an oblate spheroid, is formulated using the extended phase space method. The new independent variable, similar to the true anomaly, decouples the radius and latitude equations into two perturbed harmonic oscillators whose solutions toO(J ₂ ⁴ ) are obtained using Lindstedt's method. From these solutions and the solution to the Hamilton-Jacobi equation suitable angle variables, their canonical conjugates and the new Hamiltonian are obtained. The new Hamiltonian, accurate toO(J ₂ ⁴ ) is function of only the momenta.     

The motion of vortices within a rotating,fluid shell

We consider a spherical, solid planet surrounded by a thin layer of an incompressible, inviscid fluid. The planet rotates with constant angular velocity.Within the constraints of the geostrophic approximation of hydrodynamics, we determine the equation that governs the motion of a vortex tube within this rotating ocean. This vorticity equation turns out to be a nonlinear partial differential equation of the third order for the stream function of the motion.We next examine the existence of particular solutions to the vorticity equation that represent travelling waves of permanent form but decaying at infinity. A particular solution is obtained in terms of I ₁ and k ₁, the modified Bessel functions of order one.The question whether these localized vortices that move like solitary waves could even be solitons depends on their behavior during and after collision with each other and has not yet been resolved.Retired, U.S. Naval Research Laboratory, Washington, D.C., U.S.A.