Necessary conditions for the existence of particular algebraic first integrals

Existence of algebraic particular first integrals for a class of dynamical systems is discussed in connection with the nature of the singularities of solutions. It is shown that under some conditions, the existence of algebraic particular first integrals controls a quantity characterising a singularity (Kowaleveski's exponent) which can be calculated in a finite procedure. Four examples are given which illustrate how main theorems work.     

The eleventh motion constant of the two-body problem

The two-body problem is a twelfth-order time-invariant dynamic system, and therefore has eleven mutually-independent time-independent integrals, here referred to as motion constants. Some of these motion constants are related to the ten mutually-independent algebraic integrals of the n-body problem, whereas some are particular to the two-body problem. The problem can be decomposed into mass-center and relative-motion subsystems, each being sixth-order and each having five mutually-independent motion constants. This paper presents solutions for the eleventh motion constant, which relates the behavior of the two subsystems. The complete set of mutually-independent motion constants describes the shape of the state-space trajectories. The use of the eleventh motion constant is demonstrated in computing a solution to a two-point boundary-value problem.     

Bruns' Theorem: The Proof and Some Generalizations

We give here a proof of Bruns’ Theorem which is both complete and as general as possible: Generalized Bruns’ Theorem.In the Newtonian (n+1)-body problem in ℝ ^p with n≥2 and 1≤p≤n+1, every first integral which is algebraic with respect to positions, linear momenta and time, is an algebraic function of the classical first integrals: the energy, the p(p−1)/2 components of angular momentum and the 2p integrals that come from the uniform linear motion of the center of mass. Bruns’ Theorem only dealt with the Newtonian three-body problem in ℝ³; we have generalized the proof to n+1 bodies in ℝ^p with p≤n+1. The whole proof is much more rigorous than the previous versions (Bruns, Painlevé, Forsyth, Whittaker and Hagiara). Poincaré had picked out a mistake in the proof; we have understood and developed Poincaré’s instructions in order to correct this point (see Subsection 3.1). We have added a new paragraph on time dependence which fills in an up to now unnoticed mistake (see Section 6). We also wrote a complete proof of a relation which was wrongly considered as obvious (see Section 3.3). Lastly, the generalization, obvious in some parts, sometimes needed significant modifications, especially for the case p=1 (see Section 4). This revised version was published online in July 2006 with corrections to the Cover Date.     

About the first integrals of the generalized problem of translatory-rotary motion of rigid bodies

The existence of ten first integrals for the classical problem of the motion of a system of material points, mutually attracting according to Newtonian law, is well known.The existence of the analogous ten first integrals for the more complicated problem of the motion of a system of absolutely rigid bodies, whose elementary particles mutually attract according to the Newtonian law, was established by the author (Duboshin, 1958, 1963, 1968).In his later papers (Duboshin, 1969, 1970), the problem of the motion of a system of material points, attracting each other according to a more general law, was considered and, in particular, it was shown under what conditions the ten first integrals, analogous to the classical integrals, may exist for this problem.In the present paper, the generalized problem of translatory-rotatory motion of rigid bodies, whose elementary particles acting upon each other according to arbitrary laws of forces along the straight line joining them, is discussed.The author has shown that the first integrals for this general problem, analogous to the integrals of the problem of the translatory-rotatory motion of rigid bodies, whose elementary particles acting according to the Newtonian law, exist under certain well known conditions.That is, it has been established that if the third axiom of dynamics (action = reaction) is satisfied, then the integrals of the motion of centre of inertia and the integrals of the moment of momentum exist for this generalized problem.If the third axiom is not satisfied, then the above mentioned integrals do not exist.The third axiom is a necessary but not a sufficient condition for the existence of the tenth integral-the energy integral. The tenth integral always exists if the elementary particles of the bodies acting with a force, depend only on the mutual distances between them. In this case the force function exists for the problem and the energy integral can be expressed in a well known form.The tenth integral may exist for some more general case, without expressing the principle of conservation of energy, but permitting calculation of the kinetic energy, if the configuration of a system is given.The problem, in which the elementary particles acting according to the generalized Veber's law (Tisserand, 1896) has been cited as an example of this more general case.     

An extension of Newton's apsidal precession theorem

Newton's apsidal precession theorem in Proposition 45 of Book I of the 'Principia' has great mathematical, physical, astronomical and historical interest. The lunar theory and the precession of the perihelion of the planet Mercury are but two examples of the applications of this theorem. We have examined the precession of orbits under varying force laws as measured by the apsidal angle θ( N , e ), where N is the index for the centripetal force law, for varying eccentricity e . The paper derives a general function for the apsidal angle, dependent only on e and N as the potential is spherically symmetric. Further, we explore approximate ways of the solution of this equation, in the neighbourhood of   N = 2  which happens to be the case of greatest historical interest. Exact solutions are derived where they are possible. The first derivatives  ∂θ/∂ N   and  ∂θ/∂ h   [where h ( N , e ) is the angular momentum] are analytically expressed in the neighbourhood of   N = 2  (case of the inverse square law). The value of  ∂θ/∂ N   is computed numerically as well for  1 ≤ N < 3  . The resulting integrals are interesting improper integrals with singularities at both limits. Some of the integrals, especially for   N = 2  , can be given in closed form in terms of generalized hypergeometric functions which are reducible in terms of algebraic and logarithmic functions. No evidence was found for isolated cases of zero precession as e was increased. The   N = 1  case of the logarithmic potential is also briefly discussed in view of its interest for the dynamics of eccentric orbits and its relevance to realistic galaxy models. The possibility of apsidal precession was also examined for a few cases of high-eccentricity asteroids and extrasolar planets. We find that these systems may provide interesting new laboratories for studies of gravity.     

Solution of non-coherent spectral line transfer problem for finite atmosphere by theF n -method

By use of the orthogonality and normalization integrals developed by McCormick and Siewert (1970) a set of singular integral equations suitable forF _n -method is derived for non-coherent spectral line formation problem in finite media.F _n -equations for exit distributions are used to develop some algebraic equations with suitable recursion relations.     

The effect of orbital eccentricity and nodal regression on spin-stabilized satellites

The gyroscopic motion of a spin-stabilized satellite due to gravity gradient torques in a circular orbit has been analyzed to varying degrees in numerous publications. This paper shows that the restriction to a circular orbit is, in fact, not essential and with a slight increase in complexity, noncircular orbits may be treated. More importantly, a uniform regression of the orbital node can also be accounted for. The general results are expressed in closed form using Jacobian elliptic functions. Finally, and this is perhaps most important, certain algebraic integrals of the precession are given which can be used to place limits on the excursions of the spin axis without actually solving for the motion. This allows one to design orientations such that the maximum angle between the orbit normal and spin axis never exceeds a specific amount even though the orbit normal is in precession.     

On an analytical solution of the optimum trajectory problem in a gravitational field

Mayer's variational problem for a point with a limited mass flow rate is described by differential equations of the fourteenth order, allowing for a few first integrals. By reducing the equations to closed canonical form, these integrals are analyzed from the viewpoint of finding a possible solution to the problem via quadratures on zero, intermediate, and maximum thrust sections. In addition to confirming well-known cases of total integrability, this approach enabled us to establish that the essential difficulty of the solution of the space problem with intermediate thrust is reduced to finding one integral, and the solution of the problem with maximum thrust requires two integrals in involution. It is shown that these integrals can be applied to find particular solutions.     

Solar phase curves and phase integrals for the leading and trailing hemispheres of Iapetus from the Cassini Visual Infrared Mapping Spectrometer

We performed photometry of Cassini Visual Infrared Mapping Spectrometer observations of Iapetus to produce the first phase integrals calculated directly from solar phase curves of Iapetus for the leading hemisphere and to estimate the phase integrals for the trailing hemisphere. We also explored the phase integral dependence on wavelength and geometric albedo. The extreme dichotomy of the brightness of the leading and trailing sides of Iapetus is reflected in their phase integrals. Our phase integrals, which are lower than the results of Morrison et al. (Morrison, D., Jones, T.J., Cruikshank, D.P., Murphy, R.E. [1975]. Icarus 24, 157-171) and Squyres et al. (Squyres, S.W., Buratti, B.J., Veverka, J., Sagan, C. [1984]. Icarus 59, 426-435), have profound implications on the energy balance and volatile transport on this icy satellite.     

The formal integrals of the regularized restricted planar problem of three bodies

The algorithm for constructing the first integrals of motion of the regularized restricted planar problem of three bodies is proposed. The integrals are constructed as the formal power series in one from variables. It is shown that coefficients of these series are trigonometric polynomials of the other variable. The proposed algorithm can be realized on a computer both analytically and numerically.     

Elliptic integrals for cosmological constant cosmologies

Following the discussion of some general properties and analytical formulae for cosmological models with non-zero cosmological constant, we show how the elliptic integrals for comoving distance and light travel times as function of redshift can be expressed through Legendre integrals of the first and third kind, for which standard implementations are readily available. Observational properties are then illustrated for selected but typical models using the previously derived formulae.     

Equivalence of some integrals of the radiation theory

A definite integral which occurs in radiation theory is shown to be equal in value to another definite integral by evaluating the flux from a spherically symmetrical radiating sphere in two ways. As a corollary, an alternate proof of the invariance of the specific intensity of a ray in empty space along its path is presented.Furthermore, the equality of these two indefinite integrals leads to the conversion of members of a class of indefinite and definite integrals involving arbitrary functions of angle into other integrals. These transformations facilitate the calculation of some of these integrals which arise not only in the theory of radiation, but in other physical situations with spherical or axial symmetry — especially those in which inverse-square laws are involved.     

The formal integrals of the restricted planar problem of the three bodies

The first integrals of motion of the restricted planar circular problem of three bodies are constructed as the formal power series in r^1/2, r being the distance of a moving particle from the primary. It is shown that the coefficients of these series are trigonometric polynomials of an angular variable. Some particular solutions have been found in a closed form. The proposed method for constructing the formal integrals can be generalized to a spatial problem of three bodies.     

Volume integrals of the products of spherical harmonics and their application to viscous dissipation phenomena in fluids

The equations governing the conversion of kinetic energy into heat in moving viscous media are formulated as volume integrals of products of spherical harmonics. Although the formulation of the fundamental equations is classical, difficulties in the integration of certain products of generalized spherical harmonics over a sphere have permitted heretofore the treatment of only two cases. The closed, form evaluation of eight fundamental types of definite integrals of the product of spherical harmonics, some of them new, or at least missing in the literature, makes possible for the first time the evaluation of these volume integrals in closed form for arbitrary order and index. Explicit details are given for the rates of energy dissipation produced by viscous motions characterized by spheroidal as well as toroidal symmetry.     

Association of spherical and ellipsoidal gravity coefficients of the Earth's potential

A comparison is drawn between the expansion of the potential in spherical harmonics on the one hand and in ellipsoidal harmonics on the other, with the objective of associating the spherical and ellipsoidal gravity coefficients of the Earth's potential.For this purpose the properties of orthogonality of the Lamé functions of the first kind have been tailored to this subject of investigation and become instrumental in establishing the mathematical expressions which relate the two classes of gravity coefficients to each other. In deriving the elements of the transition matrices elliptic integrals have been encountered whose reduction to the three kinds of canonical elliptic integrals is discussed.Presented at the Conference on Celestial Mechanics, Oberwolfach, Germany, August 17–23, 1969.     

Influence of the resonance in a gravity-gradient stabilized satellite

Using Hamiltonian formalism the translational-rotational motion of a satellite is studied near a resonance considering the orbital and rotational motions. A first order perturbation theory is derived by Hori's transformation in order to eliminate short and long periodic terms, preserving in the new Hamiltonian secular and resonant terms. This theory is again applied to study the resonant system whose analysis lead us to a system of equations equivalent to the equations of a simple pendulum which is integrable in terms of elliptical integrals.     

Orbit determination with the two-body integrals. II

The first integrals of the Kepler problem are used to compute preliminary orbits starting from two short observed arcs of a celestial body, which may be obtained either by optical or by radar observations. We write polynomial equations for this problem, which can be solved using the powerful tools of computational Algebra. An algorithm to decide if the linkage of two short arcs is successful, i.e. if they belong to the same observed body, is proposed and tested numerically. This paper continues the research started in Gronchi et al. (Celest. Mech. Dyn. Astron. 107(3):299–318, 2010), where the angular momentum and the energy integrals were used. The use of a suitable component of the Laplace–Lenz vector in place of the energy turns out to be convenient, in fact the degree of the resulting system is reduced to less than half.     

Some remarks on the generalized many-body problem

This paper deals with the generalized problem of motion of a system of a finite number of bodies (material points).We suppose here that every point of the system acts on another one with a force (attractive or repulsive) directed along the straight line connecting these two points, and proportional to the product of their masses and a certain function of time, mutual distance and its derivatives of the first and second order (Duboshin, 1970).The laws of forces are different for different pairs of points and, generally speaking, the validity of the third axiom of dynamics (law of action and reaction) is not assumed in advance.With these general assumptions we find the conditions for the laws of the forces under which the problem admits the first integrals, analogous to the classic integrals of the many-body problem with the Newton's law of attraction.It is shown furthermore, that in this generalized problem it is possible to obtain an equation, analogous to the classic equation of Lagrange-Jacobi and deduce the conditions of stability or instability of the system in Lagrange's sense.The results obtained may be applied for the investigation of motion in some isolated stellar systems, where the laws of mechanics may be different from the laws in our solar system.     

The use of integrals in numerical integrations of theN-body problem

The numerical integration of systems of differential equations that possess integrals is often approached by using the integrals to reduce the number of degrees of freedom or by using the integrals as a partial check on the resulting solution, retaining the original number of degrees of freedom.Another use of the integrals is presented here. If the integrals have not been used to reduce the system, the solution of a numerical integration may be constrained to remain on the integral surfaces by a method that applies corrections to the solution at each integration step. The corrections are determined by using linearized forms of the integrals in a least-squares procedure.The results of an application of the method to numerical integrations of a gravitational system of 25-bodies are given. It is shown that by using the method to satisfy exactly the integrals of energy, angular momentum, and center of mass, a solution is obtained that is more accurate while using less time of calculation than if the integrals are not satisfied exactly. The relative accuracy is ascertained by forward and backward integrations of both the corrected and uncorrected solutions and by comparison with more accurate integrations using reduced step-sizes.     

Bilinear integrals of the radiative transfer equation

It is shown that the group of problems in the theory of radiative transfer that are reducible to the sourcefree problem admits a class of integrals involving quadratic moments of the intensity of arbitrarily high orders. Based on a variational principle, it is found that these integrals, which include the R-integral, follow from the corresponding conservation laws. Some of the results are generalized to the case of anisotropic scattering.