The theory of asynchronous relative motion II: universal and regular solutions

Two fully regular and universal solutions to the problem of spacecraft relative motion are derived from the Sperling–Burdet (SB) and the Kustaanheimo–Stiefel (KS) regularizations. There are no singularities in the resulting solutions, and their form is not affected by the type of reference orbit (circular, elliptic, parabolic, or hyperbolic). In addition, the solutions to the problem are given in compact tensorial expressions and directly referred to the initial state vector of the leader spacecraft. The SB and KS formulations introduce a fictitious time by means of the Sundman transformation. Because of using an alternative independent variable, the solutions are built based on the theory of asynchronous relative motion. This technique simplifies the required derivations. Closed-form expressions of the partial derivatives of orbital motion with respect to the initial state are provided explicitly. Numerical experiments show that the performance of a given representation of the dynamics depends strongly on the time transformation, whereas it is virtually independent from the choice of variables to parameterize orbital motion. In the circular and elliptic cases, the linear solutions coincide exactly with the results obtained with the Clohessy–Wiltshire and Yamanaka–Ankersen state-transition matrices. Examples of relative orbits about parabolic and hyperbolic reference orbits are also presented. Finally, the theory of asynchronous relative motion provides a simple mechanism to introduce nonlinearities in the solution, improving its accuracy.     

Two-body motion in the post-Newtonian approximation

The motion of two massive particles is considered within the framework of the first post-Newtonian approximation. The system Hamiltonian is constructed and normalized through first order using a canonical transformation method of implicit variables. Closed-form solutions for the Delaunay elements in the phase space are obtained. The bridge between the phase space and the state space of the Lagrangian of the motion is provided by a velocity-dependent Legendre transformation. By explicit inversion of this transformation, expressions for the Keplerian elements in the state space are obtained from the Delaunay element solutions.     

Exact solution to the relative orbital motion in eccentric orbits

This paper studies the relative orbital motion between arbitrary Keplerian trajectories. A closed-form vectorial solution to the nonlinear initial value problem that models this type of motion with respect to a noninertial reference frame is offered. Without imposing any particular conditions on the leader or the deputy satellites trajectories, exact expressions for the relative law of motion and relative velocity are obtained in a closed form. This solution allows the parameterization of the relative motion manifold and offers new methods to study its geometrical and topological properties. The result presented in this paper opens the way to obtain new classes of approximate solutions to the equations of relative motion with time, an eccentric or true anomaly as independent variables. Published in Russian in Solar System Research, 2009, Vol. 43, No. 1, pp. 44–55. The text was submitted by the autors in English.     

Analytical solutions around planar equilibria of three dipole problem

Proceeding with our investigation into the motion of a particle influenced by the electromagnetic field of three celestial bodies of a magnetic-dipole nature we give here for the first time the analytical expressions of periodic solutions around a planar equilibrium point. These relations are expansions of the planar equations of motion in series of second order power of a parameter in the vicinity of equilibria. The above analytical expressions of periodic solutions give the first members of the family of periodic orbits which emanate from a stable equilibrium point. The whole family can then be calculated using a predictor-corrector algorithm.     

Quasiintegrals of the photogravitational eccentric restricted three-body problem with Poynting–Robertson drag

A restricted three-body problem for a dust particle, in presence of a spherical cometary nucleus in an eccentric (elliptic, parabolic or hyperbolic) orbit about the Sun, is considered. The force of radiation pressure and the Poynting– Robertson effect are taken into account. The differential equations of the particle’s non-inertial spatial motion are investigated both analytically and numerically. With the help of a complex representation, a new single equation of the motion is obtained. Conversion of the equations of motion system into a single equation allows the derivation of simple expressions similar to the integral of energy and integrals of areas. The derived expressions are named quasiintegrals. Relative values of terms of the energy quasiintegral for a smallest, largest, and a mean comet are calculated. We have found that in a number of cases the quasiintegrals are related to the regular integrals of motion, and discuss how the quasiintegrals may be applied to find some significant constraints on the motion of a body of infinitesimal mass.     

Analytical solutions for J 2-perturbed unbounded equatorial orbits

While solutions for bounded orbits about oblate spheroidal planets have been presented before, similar solutions for unbounded motion are scarce. This paper develops solutions for unbounded motion in the equatorial plane of an oblate spheroidal planet, while taking into account only the J ₂ harmonic in the gravitational potential. Two cases are distinguished: A pseudo-parabolic motion, obtained for zero total specific energy, and a pseudo-hyperbolic motion, characterized by positive total specific energy. The solutions to the equations of motion are expressed using elliptic integrals. The pseudo-parabolic motion unveils a new orbit, termed herein the fish orbit, which has not been observed thus far in the perturbed two-body problem. The pseudo-hyperbolic solutions show that significant differences exist between the Keplerian flyby and the flyby performed under the the J ₂ zonal harmonic. Numerical simulations are used to quantify these differences.     

Periodic solutions about the collinear lagrangian solution in the general problem of three bodies

The article describes the solutions near Lagrange's circular collinear configuration in the planar problem of three bodies with three finite masses. The article begins with a detailed review of the properties of Lagrange's collinear solution. Lagrange's quintic equation is derived and several expressions are given for the angular velocity of the rotating frame.The equations of motion are then linearized near the circular collinear solution, and the characteristic equation is also derived in detail. The different types of roots and their corresponding solutions are discussed. The special case of two equal outer masses receives special attention, as well as the special case of two small outer masses.Finally, the fundamental family of periodic solutions is extended by numerical integration all the wap up to and past a binary collision orbit. The stability and the bifurcations of this family are briefly enumerated.     

Lagrange’s planetary equations for the motion of electrostatically charged spacecraft

A spacecraft that generates an electrostatic charge on its surface in a planetary magnetic field will be subject to a perturbative Lorentz force. Active modulation of the surface charge can take advantage of this electromagnetic perturbation to modify or to do work on the spacecraft’s orbit. Lagrange’s planetary equations are derived using the Lorentz force as the perturbation on a Keplerian orbit, incorporating orbital inclination and true anomaly for the first time for an electrostatically charged vehicle. The planetary equations reveal that orbital inclination is a second-order effect on the perturbation, explaining results found in earlier studies through numerical integration. All of the orbital elements are coupled, but the coupling notably does not depend on the magnitude of the electrostatic charge or on the strength of the magnetic field. Analytical expressions that characterize this coupling are tested with a propellantless escape example at Jupiter. A closed-form solution exists that constrains the set of equatorial orbits for which planetary escape is possible, and a sufficient condition is identified for escape from inclined orbits. The analytical solutions agree with results from the numerically integrated equations of motion to within a fraction of a percent.     

Analytical solutions for the equations of motion of a space vehicle during the atmospheric re-entry phase on a 2-D trajectory

A practical and important problem encountered during the atmospheric re-entry phase is to determine analytical solutions for the space vehicle dynamical equations of motion. The author proposes new solutions for the equations of trajectory and flight-path angle of the space vehicle during the re-entry phase in Earth’s atmosphere. Explicit analytical solutions for the aerodynamic equations of motion can be effectively applied to investigate and control the rocket flight characteristics. Setting the initial conditions for the speed, re-entering flight-path angle, altitude, atmosphere density, lift and drag coefficients, the nonlinear differential equations of motion are linearized by a proper choice of the re-entry range angles. After integration, the solutions are expressed with the Exponential Integral, and Generalized Exponential Integral functions. Theoretical frameworks for proposed solutions as well as, several numerical examples, are presented.     

Orbits in an anisotropic radiation field

The elliptic-type motion around a source whose luminosity is anisotropic is being perturbatively treated. Using a Fourier expansion for the perturbing force, exact expressions describing the evolution of each orbital element are determined from Newton-Euler equations. For near-unit frequency the resonant solutions are pointed out. An approximate expression for the variation of the nodal period is determined, too. The behaviour of the radius vector of the unstable orbit (resonance case) is pointed out, showing periodic variations of constantly increasing amplitude. The (both angular and physical) time scale for escape is estimated. Concrete astronomical situations modellable in this way are mentioned.     

第三体摄动分析解的一种表达式

在太阳系中,大行星、小行星和卫星（包括自然卫星和人造卫星）等对应的运动问题,都可以处理成受摄二体问题,而摄动源又多为第三体,作为第三体的摄动天体,有的比运动天体离中心天体近,有的则相反,前者称为内摄内体,全者则称为外摄天体,对一个具体的运动天体,可以同时出现这两个摄动天体,但是,只要运动天体与摄动天体的轨道都建立在以中心天体（质心）为坐标原点的同一坐标系内,那么在一定条件下（即除运动天体与摄动天体     

A note on a Lagrangian formulation for motion about the collinear points

A lagrangian formulation for the three-dimensional motion of a satellite in the vicinity of the collinear points of the circular-restricted problem is reconsidered. It is shown that the influence of the primaries can be expressed in the form of two third-body disturbing functions. By use of this approach, the equations for the Lagrangian and for the motion itself are readily developed into highly compact expressions. All orders of the non-linear developments are shown to be easily obtainable using well-known recursive relationships. The resulting forms for these equations are well suited for use in the initial phase of canonical or non-canonical investigations.     

Tidal interactions in binary stellar systems

The tidal interactions in binary stellar systems are studied under the assumption that the orbital motion of the binary is negligible in comparison with the stellar motion. By integrating over time the tidal forces acting on the stars, the energy changes are derived. These are used to obtain simple analytical expressions for the rates of disruption and merging. This method gives appropriate value for the Roche density ρ_r and it is found that the disruption rate of a satellite of density ρ changes drastically at ρ ρ_R A comparison is made with earlier results obtained under the simplifying assumption that stellar motion is negligible in comparison with the orbital motion of the binary and its implications are discussed.     

Periodic solution of the nonlinear Sitnikov restricted three-body problem

The objective of this paper is to find periodic solutions of the circular Sitnikov problem by the multiple scales method which is used to remove the secular terms and find the periodic approximated solutions in closed forms. Comparisons among a numerical solution (NS), the first approximated solution (FA) and the second approximated solution (SA) via multiple scales method are investigated graphically under different initial conditions. We observe that the initial conditions play a vital role in the numerical and approximated solutions behaviour. The obtained motion is periodic, but the difference of its amplitude is directly proportional with the initial conditions. We prove that the obtained motion by the numerical or the second approximated solutions is a regular and periodic, when the infinitesimal body starts its motion from a nearer position to the common center of primaries. Otherwise when the start point distance of motion is far from this center, the numerical solution may not be represent a periodic motion for along time, while the second approximated solution may present a chaotic motion, however it is always periodic all time. But the obtained motion by the first approximated solution is periodic and has regularity in its periodicity all time. Finally we remark that the provided solutions by multiple scales methods reflect the true motion of the Sitnikov restricted three–body problem, and the second approximation has more accuracy than the first approximation. Moreover the solutions of multiple scales technique are more realistic than the numerical solution because there is always a warranty that the motion is periodic all time.     

On the transverse vibrations of a revolving tether

We analyse the transverse vibrations of a tether, modelled as an inextensible cable, and revolving at an average rate equal to the orbital rate. The reference motion is a revolving rigid tether. During this motion the force in the tether (time and location dependent) remains, in a first approximation, aligned with the tether axis. Separation of variables for the vibrations about this motion gives a Legendre equation for the spatial dependency of the deformations and Hill's equations for time dependency of the in- and out-of-plane deformations. The boundary conditions on the Legendre equation generate a series of admissible values of the separation constant that become equidistant. The two Hill's equations generate a series of intervals, contracting to equidistant critical values, where the solutions are unbounded. The admissible values of the separation constant must avoid these intervals. Asymptotic expressions for the separation constant and the critical values are given. The first and second in-plane deformation mode arc unstable for zero end masses. By increasing the ratio of the concentrated over the distributed mass the deformation modes can be stabilised and the values of the separation constant can be made a multiple of the distribution of the critical points. Introducing unequal tip masses does not affect this result.     

Exact analytic solutions for the rotation of an axially symmetric rigid body subjected to a constant torque

New exact analytic solutions are introduced for the rotational motion of a rigid body having two equal principal moments of inertia and subjected to an external torque which is constant in magnitude. In particular, the solutions are obtained for the following cases: (1) Torque parallel to the symmetry axis and arbitrary initial angular velocity; (2) Torque perpendicular to the symmetry axis and such that the torque is rotating at a constant rate about the symmetry axis, and arbitrary initial angular velocity; (3) Torque and initial angular velocity perpendicular to the symmetry axis, with the torque being fixed with the body. In addition to the solutions for these three forced cases, an original solution is introduced for the case of torque-free motion, which is simpler than the classical solution as regards its derivation and uses the rotation matrix in order to describe the body orientation. This paper builds upon the recently discovered exact solution for the motion of a rigid body with a spherical ellipsoid of inertia. In particular, by following Hestenes’ theory, the rotational motion of an axially symmetric rigid body is seen at any instant in time as the combination of the motion of a “virtual” spherical body with respect to the inertial frame and the motion of the axially symmetric body with respect to this “virtual” body. The kinematic solutions are presented in terms of the rotation matrix. The newly found exact analytic solutions are valid for any motion time length and rotation amplitude. The present paper adds further elements to the small set of special cases for which an exact solution of the rotational motion of a rigid body exists.     

Approximate solutions of non-linear circular orbit relative motion in curvilinear coordinates

A compact, time-explicit, approximate solution of the highly non-linear relative motion in curvilinear coordinates is provided under the assumption of circular orbit for the chief spacecraft. The rather compact, three-dimensional solution is obtained by algebraic manipulation of the individual Keplerian motions in curvilinear, rather than Cartesian coordinates, and provides analytical expressions for the secular, constant and periodic terms of each coordinate as a function of the initial relative motion conditions or relative orbital elements. Numerical test cases are conducted to show that the approximate solution can be effectively employed to extend the classical linear Clohessy–Wiltshire solution to include non-linear relative motion without significant loss of accuracy up to a limit of 0.4–0.45 in eccentricity and 40–45\(^\circ \) in relative inclination for the follower. A very simple, quadratic extension of the classical Clohessy–Wiltshire solution in curvilinear coordinates is also presented.     

Anharmonicity effects in bounce motion of particles

Equations for parallel motion for a particle trapped in a magnetic field have been considered and improved solutions of differential equations have been derived. The expressions for the change in energy of the particle (Δw) and diffusion coefficient (D_ww) have been presented in a simple form using the improved solution.     

A core-mantle interaction in the rotation of the Earth

Effects of an interaction between the mantle and the core of the Earth on its rotational motion are investigated. Assuming that the Earth consists of a rigid mantle and a rigid core with a frictional coupling and a kind of inertial coupling between them, the equations of motion are derived, and they are solved in a close approximation. The solution gives the expressions for the precession, the nutation, the secular changes in the obliquity and the rotational speed, the polar motion and so on as functions of the magnitudes of these forces. A numerical estimation shows that the effect of the friction on the amplitude and phase of the nutation is small for a reasonable intensity of the friction while inertial coupling force has a decisive influence on the amplitude, and an appropriately chosen value of the latter force gives a nutation which closely agrees with observations. It is also indicated that this torque remarkably lessens the rates of the secular changes in the obliquity and the rotational speed. The possibility of a periodical change in the amplitude of the polar motion is suggested as a result of the interaction between the two consituents.