Notes on Hori Method for Non-Canonical Systems

In this paper the new approach for the integration theory of the canonical version of Hori method recently proposed is extended to the non-canonical one. It will be shown that the non-homogeneous ordinary differential equation with an auxiliary parameter t* associated with the mth order equation of the algorithm can also be replaced by a non-homogeneous partial differential equation in the time t. Using a generalized canonical approach, the general algorithm proposed by Sessin is then revised; as well as the Lagrange variational equations for the non-canonical version of Hori method. A simplified algorithm derived from Sessin's algorithm is presented for non-linear oscillations problem.     

On stationary solutions of some canonical systems of differential equations

In our article (Zhuravlev, 1979) a formal method of constructing conditionally periodic solutions of canonical systems of differential equations with a quick-rotating phase in the case of sharp commensurability was presented. The existence of stationary (or periodic) solutions of an averaged system of differential equations corresponding to the initial system of differential equations is necessary for an effective application of the method for different problems.Evidently, the stationary solutions do not always exist but in numerous papers on stationary solutions (oscillations or motions), the conditions of existence of such solutions are very often not considered at all. Usually a simple assumption is used that the stationary solutions do exist.Otherwise it is well known that Poincaré's theory of periodic solutions (Poincaré, 1892) let one set up conditions of existence of periodic solutions in different systems of differential equations. Particularly, in papers,Mah (1949, 1956), see alsoexmah (1971), the necessary and sufficient conditions of the existence of periodic solutions of (non-canonical) systems of differential equations which are close to arbitrary non-linear systems are given. For canonical autonomous systems of differential equations the conditions of existence of periodic solutions and a method of calculation are presented in the paperMepmah (1952).In our paper another approach is given and the conditions of existence of stationary solutions of canonical systems of differential equations with a quick-rotating phase are proved. For this purpose Delaunay-Zeipel's transformation and Poincaré's small parameter method are used.     

An analytic method to account for drag in the Vinti satellite theory

In order to retain separability in the Vinti theory of Earth satellite motion when a nonconservative force such as air drag is considered, a set of variational equations for the orbital elements are introduced, and expressed as functions of the transverse, radial, and normal components of the nonconservative forces acting on the system. In this approach, the Hamiltonian is preserved in form, and remains the total energy, but the initial or boundary conditions and hence the Jacobi constants of the motion advance with time through the variational equations. In particular, the atmospheric density profile is written as a fitted exponential function of the eccentric anomaly, which adheres to tabular data at all, altitudes and simultaneously reduces the variational equations to definite integrals with closed form evaluations whose limits are in terms of the eccentric anomaly. The values of the limits for any arbitrary time interval are obtained from the Vinti program.Results of this technique for the case of the intense air drag satellites San Marco-2 and Air Force Cannonball are given. These results indicate that the satellite ephemerides produced by this theory in conjunction with the Vinti program are of very high accuracy. In addition, since the program is entirely analytic, several months of ephemerides can be obtained within a few seconds of computer time.     

A note on the solution of the variational equations of a class of dynamical systems

Some properties are derived for the solutions of the variational equations of a class of dynamical systems. It is shown that in rather general conditions the matrix of the linearized Lagrangian equations of motion have an important property for which the word skew-symplectic has been introduced. It is also shown that the fundamental matrix of solutions is symplectic, the word symplectic being used here in a more general sense than in the classical literature. Two consequences of the symplectic property are that the fundamental matrix is easily invertible and that the eigenvalues appear in reciprocal pairs. The effect of coordinate transformations is also analyzed; in particular the change from Lagrangian to canonical systems.     

Similarity considerations and conservation laws for magneto-static atmospheres

The equations of magnetohydrostatic equilibria for a plasma in a gravitational field are investigated analytically. For equilibria with one ignorable spatial coordinate, the equations reduce to a single nonlinear elliptic equation for the magnetic potential A. Similarity solutions of the elliptic equation are obtained for the case of an isothermal atmosphere in a uniform gravitational field. The solutions are obtained from a consideration of the invariance group of the elliptic equation. The importance of symmetries of the elliptic equation also appears in the determination of conservation laws. It turns out that the elliptic equation can be written as a variational principle, and the symmetries of the variational functional lead (via Noether's theorem) to conservation laws for the equation. As an example of the application of the similarity solutions, we construct a model magnetostatic atmosphere in which the current density J is proportional to the cube of the magnetic potential, and falls off exponentially with distance vertical to the base, with an e-folding distance equal to the gravitational scale height. The solutions show the interplay between the gravitational force, the J × B force (B, magnetic field induction) and the gas pressure gradient.     

Notes on Hori Method for Canonical Systems

In this paper a slightly different approach is proposed for the process of determining the functions S _m and H _m ^* of the algorithm of the canonical version of Hori method. This process will be referred to as integration theory of the mth order equation of the method. It will be shown that the ordinary differential equation with an auxiliary parameter t ^* as independent variable, introduced through Hori auxiliary system, can be replaced by a partial differential equation in the time t. In this way, the mth order equation of the algorithm assumes a form very similar to the one of other perturbation methods. In virtue of this new approach of the integration theory for Hori method, Lagrange's variational equations introduced by Sessin are revised. As an example, the Duffing equation is solved through this new approach.     

The elimination of the critical terms of a first order Uranus-Neptune theory by Hori's method

In this part we determine the value ofS ₁, and in terms of the canonical variables of H. Poincaré. A complete solution of the auxiliary system of equations generated by the Hamiltonian is presented.     

Canonical solution of the two critical argument problem

The solution by Sessin and Ferraz-Mello (Celes. Mech. 32, 307–332) of the Hori auxiliary system for the motion of two planets with periods nearly commensurate in the ratio 21 is considerably simplified by the introduction of canonical variables. An analogous canonical transformation simplifies the elliptic restricted problem.     

Radiative Transfer in Inhomogeneous Atmospheres. I

This series of papers is devoted to multiple scattering of light in plane parallel, inhomogeneous atmospheres. The approach proposed here is based on Ambartsumyan's method of adding layers. The main purpose is to show that one can avoid difficulties with solving various boundary value problems in the theory of radiative transfer, including some standard problems, by reducing them to initial value problems. In this paper the simplest one dimensional problem of diffuse reflection and transmission of radiation in inhomogeneous atmospheres with finite optical thicknesses is considered as an example. This approach essentially involves first determining the reflection and transmission coefficients of the atmosphere, which, as is known, are a solution of the Cauchy problem for a system of nonlinear differential equations. In particular, it is shown that this system can be replaced with a system of linear equations by introducing auxiliary functions P and S. After the reflectivity and transmissivity of the atmosphere are determined, the radiation field in it is found directly without solving any new equations. We note that this approach can be used to obtain the required intensities simultaneously for a family of atmospheres with different optical thicknesses. Two special cases of the functional dependence of the scattering coefficient on the optical thickness, for which the solutions of the corresponding equations can be expressed in terms of elementary functions, are examined in detail. Some numerical calculations are presented and interpreted physically to illustrate specific features of radiative transport in inhomogeneous atmospheres.     

Forbidden Zones for Circular Regular Orbits of the Moons in Solar System,R3BP

Previously, we have considered the equations of motion of the three-body problem in a Lagrange form (which means a consideration of relative motions of 3-bodies in regard to each other). Analysing such a system of equations, we considered the case of small-body motion of negligible mass m ₃ around the second of two giant-bodies m ₁, m ₂ (which are rotating around their common centre of masses on Kepler’s trajectories), the mass of which is assumed to be less than the mass of central body. In the current development, we have derived a key parameter η that determines the character of quasi-circular motion of the small third body m ₃ relative to the second body m ₂ (planet). Namely, by making several approximations in the equations of motion of the three-body problem, such the system could be reduced to the key governing Riccati-type ordinary differential equations. Under assumptions of R3BP (restricted three-body problem), we additionally note that Riccati-type ODEs above should have the invariant form if the key governing (dimensionless) parameter η remains in the range 10^?2 Open image in new window 10^?3. Such an amazing fact let us evaluate the forbidden zones for Moon’s orbits in the inner solar system or the zones of distances (between Moon and Planet) for which the motion of small body could be predicted to be unstable according to basic features of the solutions of Riccati-type.     

Exact anisotropic polytropic cylindrical solutions

In this paper, we study anisotropic compact stars with static cylindrically symmetric anisotropic matter distribution satisfying polytropic equation of state. We formulate the field equations as well as the corresponding mass function for the particular form of gravitational potential \(z(x)=(1+bx)^{\eta }~(\eta =1,~2,~3)\) and explore exact solutions of the field equations for different values of the polytropic index. The values of arbitrary constants are determined by taking mass and radius of compact star (Her X-1). We find that resulting solutions show viable behavior of physical parameters (density, radial as well as tangential pressure, anisotropy) and satisfy the stability condition. It is concluded that physically acceptable solutions exist only for \(\eta =1,~2\).     

Introduction to the restricted Jupiter orbiter problem

We consider a restricted six-body problem, consisting of Jupiter, the four Galilean satellites, and an orbiter. The Galilean satellites' orbits are circular and coplanar; Io, Europa, and Ganymede are in exact resonance; their mean longitudes obey the Laplace relation. We seek periodic orbits which avoid close approaches to any satellite; such orbits are of interest for mission planning. They are approximated as equilibrium points of sets of variational equations associated with time-averaged disturbing functions. Stability of the solutions is also determined. The orbits of greatest interest are:Planar: twice Callisto's period, eccentricity0.6Planar: four times Callisto's period, eccentricity0.75Slightly inclined: twice Callisto's period, eccentricity arbitraryPlanar: 4/5 or 5/4 Europa's period.     

Slow rotational perturbations of Friedmann universes

In this work the rotational perturbations of the Friedmann universes are investigated. In the general case where none of the terms including (r, t) are neglected, for perfect fluid, the field equations belonging to the perturbed metric give =(t). In this case, since the condition =0 can be accomplished by a coordinate transformation, the solutions of the field equations reduce to those of the classical Friedmann equations. For this reason, the approximate solutions obtained by other authors become formal solutions without physical interest.     

Solution of canonical equations which result from elimination of short-periodic terms in second-order planetary theory

We solve the first order non-linear differential equation and we calculate the two quadratures to which are reduced the canonical differential equations resulting from the elimination of the short period terms in a second order planetary theory carried out through Hori's method and slow Delaunay canonical variables when powers of eccentricities and the sines of semi-inclinations which are >3 are neglected and the eccentricity of the disturbing planet is identically equal to zero. The procedure can be extended to the case when the eccentricity of the disturbing planet is not identically equal to zero. In this latter general case, we calculatedthe two quadratures expressing angular slow Delaunay canonical variable ₁ of the disturbed planet and angular slow Delaunay canonical variable ₂ of the disturbing planet in terms of timet.     

Tidal evolution in close binary systems,II

The aim of the present investigation will be to determine the explicit forms of differential equations which govern secular perturbations of the orbital elements of close binary systems in the plane of the orbit (i.e., of the semi-major axisA, eccentricitye, and longitude of the periastron ), arising from the lag of dynamical tides due to viscosity of stellar material. The results obtained are exact for any value of orbital eccentricity comprised between 0e<1; and include the effects produced by the second, third and fourth-harmonic dynamical tides, as well as by axial rotation with arbitrary inclination of the equator to the orbital plane.In Section 2 following brief introductory remarks the variational equations of the problem of plane motion will be set up in terms of the rectangular componentsR, S, W of disturbing accelerations with respect to a revolving system of coordinates. The explicit form of these coefficients will be established in Section 3 to the degree of accuracy to which squares and higher powers of quantities of the order of superficial distortion can be ignored. Section 4 will be devoted to a derivation of the explicit form of the variational equations for the case of a perturbing function arising from axial rotation; and in Section 5 we shall derive variational equations which govern the perturbation of orbital elements caused by lagging dynamical tides.Numerical integrations of these equations, which govern the tidal evolution of close binary systems prompted by viscous friction at constant mass, are being postponed for subsequent investigations.Prepared at the Lunar Science Institute, Houston, Texas, under the joint support of the Universities Space Research Association, Charlottesville, Virginia, and the National Aeronautics and Space Administration Manned Spacecraft Center, Houston, Texas, under Contract No. NSR 09-051-001. This paper constitutes Lunar Science Institute Contribution no. 100.Normally at the Department of Astronomy, University of Manchester, England.     

On the ejection solutions in a central force field

In this paper we consider the problem concerning the reduction of the two-body motion to that of a single particle in a central field. As a force function we takeU(r)=r ^–, where is some positive real number. Making use of the variational equations we study the ejection solutions of the differential equations of motion.

---

Resumé Nous considérons dans cet article le problème concernant la réduction du mouvement de deux corps à celui d'une particule dans un champ de forces central. Comme fonction de forces nous prenonsU(r)=r ^–; où est un réel positif. Nous étudions à l'aide des équations aux variations les solutions d'éjection des équations du mouvement.

     

Analytic contributions to the Störmer problem

The purpose of the paper is to furnish two analytic contributions to the Störmer problem (motion of a charged particle in a magnetic dipole field). The fundamental role in these considerations is played by a quantity which depends on the constant velocity of the particle: by treating this quantity or its reciprocal as a small parameter, power series solutions can be obtained for the orbit projection upon a meridian plane valid for both low-energy and high-energy particles.For high-energy particles, the zero-order approximation is in general an ellipse, but will reduce to a straightline for particular values of the integration constants.For low-energy particles, the zero-order approximation is given as an infinite series in the variablew=sin², being the latitude, and it is shown that this basic orbit is not a magnetic line of force.Leading terms of the expansions have been given for the first-order approximations in both cases. Higher-order approximations can be obtained recurrently by solving linear, first-order differential equations having the same integrating factor, which depends on the zero-order approximation alone.The method is suitable for extensive numerical work in conjunction with a computer program of the formac class.     

The stability parameters of three-dimensional periodic three-body orbits

Formulae containing the elements of the variational matrix are obtained which determine the linear isoenergetic stability parameters of three-dimensional periodic orbits of the general three-boy problem. This requires the numerical integration of the variational equations but produces the stability parameters with the effective accuracy of the numerical integration. The conditions for stability, criticality, and bifurcations are briefly examined and the stability determination procedure is tested in the determination of some three-dimensional periodic orbits of low inclination bifurcating from vertical-critical coplanar orbits.     

Computer implementation of a new approach to the ideal resonance problem

A new approach to the librational solution of the Ideal Resonance Problem has been devised--one in which a non-canonical transformation is applied to the classical Hamiltonian to bring it to the form of the simple harmonic oscillator. Although the traditional form of the canonical equations of motion no longer holds, a quasi-canonical form is retained in this single-degree-of-freedom system, with the customary equations being multiplied by a non-constant factor. While this makes the resulting system amenable to traditional transformation techniques, it must then be integrated directly. Singularities of the transformation in the circulation region limit application of the method to the librational region of motion.Computer-assisted algebra has been used in all three stages of the solution to fourth order of this problem: using a general-purpose FORTRAN program for the quadratic analytical solution of Hamiltonians in action-angle variables, the initial transformation is carried out by direct substitution and the resulting Hamiltonian transformed to eliminate angular variables. The resulting system of differential equations, requiring the expected elliptic functions as part of their solution, is currently in the process of being integrated using the LISP-based REDUCE software, by programming the required recursive rules for elliptic integration.Basic theory of this approach and the computer implementation of all these techniques is described. Extension to higher order of the solution is also discussed.     

The precession and nutation of deformable bodies,III

In preceding papers of this series (Kopal, 1968; 1969) the Eulerian equations have been set up which govern the precession and nutation of self-gravitating fluid globes of arbitrary structures in inertial coordinates (space-axes) as well as with respect to the rotating body axes; with due account being taken of the effects arising from equilibrium as well as dynamical tides.In Section 1 of the present paper, the explicit form of these equations is recapitulated for subsequent solations. Section 2 contains then a detailed discussion of the coplanar case (in which the equation of the rotating configuration and the plane of its orbit coincide with the invariable plane of the system); and small fluctuations in the angular velocity of axial rotation arising from the tidal breathing in eccentric binary systems are investigated.In Section 3, we consider the angular velocity of rotation about theZ-axis to be constant, but allow for finite inclination of the equator to the orbital plane. The differential equations governing such a problem are set up exactly in terms of the time-dependent Eulerian angles and , and their coefficients averaged over a cycle. In Section 4, these equations are linearized by the assumption that the inclinations of the equator and the orbit to the invariable plane of the system are small enough for their squares to be negligible; and the equations of motion reduced to their canonical form.The solution of these equations — giving the periods of precession and nutation of rotating components of close binary systems, as well as the rate of nodal regression which is synchronised with precession — are expressed in terms of the physical properties of the respective system and of its constituent components; while the concluding Section 6 contains a discussion of the results, in which the differences between the precession and nutation of rigid and fluid bodies are pointed out.