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.     

Determination of the potential in a synodic system

Determination of the potential field in a fixed (inertial) system may be accomplished by the solution of a homogeneous linear partial differential equation when a family of orbits of a body moving in the field is given. This partial differential equation was presented and thoroughly analyzed earlier. The present paper discusses the same problem in a rotating system where the centrifugal and Coriolis effects render the pertinent partial differential equation in general non-homogeneous and non-linear. A linear, though non-homogeneous, partial differential equation for the determination of the synodic potential is obtained only in the special case of iso-energetic families of orbits.     

An approximate integral and solution for eccentric planetary type orbits in the restricted problem of three bodies

The circular restricted problem of three bodies is investigated analytically with respect to the problem of deriving a second integral of motion besides the well known Jacobian Integral. The second integral is searched for as a correction the angular momentum integral valid in the two body case. A partial differential equation equivalent to the problem is derived and solved approximately by an asymptotic Fourier method assuming either sufficiently small values for the dimensionless mass parameter or sufficiently large distances from the barycentre. The solution of the partial equation then leads to a function of the coordinates, velocities and time being nearly constant, which means that its variation with time is about 40–300 times less than that of the pure angular momentum. By averaging over the remaining fluctuating part of the quasi-integral we are able to integrate the first order equations using a renormalization transformation. This leads to an explicit expression for the approximate solution of the circular problem which describes the motion of the third body orbiting both primaries with nonvanishing initial eccentricity (eccentric planetary type orbits). One of the main results is an explicit formula for the frequency of the perihelion motion of the third body which depends on the mass parameter, the initial distance of the third body from the barycentre and the initial eccentricity. Finally we study orbits of the P-Type, being defined as solutions of the restricted problem with circular initial conditions (vanishing initial eccentricity).     

Effects of partial frequency redistribution on the level population densities in a resonance line

We have obtained a simultaneous solution of the statistical equilibrium equation for a non-LTE two-level atom and the radiative transfer equation in the comoving frames by employing the angle-averaged partial frequency redistribution.R _i with isotropic scattering. In the first iteration we have set the population density of the upper level equal to zero and allow it to be populated in the subsequent iterations. The solution converges within two to four iterations. The process of iteration is terminated when the ratios of population densities in two successive iterations at each radial point, attain an accuracy of 1%. The effects of partial frequency redistribution is to increase the population density of the upper level. Radial gas motions do not seem to have significant effects, although in highly extend geometries, velocity gradients change the population densities considerably.     

A theory of integration for the extended Delaunay method

In this paper a method for the integration of the equations of the extended Delaunay method is proposed. It is based on the equations of the characteristic curves associated with the partial differential equation of Delaunay-Poincaré. The use of the method of characteristics changes the partial differential equation for higher order approximations into a system of ordinary differential equations. The independent variable of the equations of the characteristics is used instead of the angular variables of the Jacobian methods and the averaging principle of Hori is applied to solve the equations for higher orders. It is well known that Jacobian methods applied to resonant problems generally lead to the singularity of Poincaré. In the ideal resonance problem, this singularity appears when higher order approximations of the librational motion are considered. The singularity of Poincaré is non-essential and is caused by the choice of the critical arguments as integration variables. The use of the independent variable of the equation of the characteristics in the place of the critical angles eliminates the singularity of Poincaré.     

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     

Inverse problem with two-parametric families of planar orbits

As a possible extension of recent work we study the following version of the inverse problem in dynamics: Given a two-parametric familyf(x, y, b)=c of plane curves, find an autonomous dynamical system for which these curves are orbits.We derive a new linear partial differential equation of the first order for the force componentsX(x, y) andY(x, y) corresponding to the given family. With the aid of this equation we find that, depending on the given functionf, the problem may or may not have a solution. Based on given criteria, we present a full classification of the various cases which may arise.     

The Fragmentation of Asteroids – A Synopsis of Analytical Methods (Mitteilungen der Universitäts-Sternwarte Jena Nr. 119)

Basing on the author's work a review of the possibilities as well as the limits of treating the problem of the collisional history of the asteroids by analytical methods is given. Using empirical data on rock fragmentation and general principles like symmetry and mass conservation the distribution function of the fragments arising from a single collision is analytically formulated. The size distribution of asteroids adjusting when crushing collisions have taken place a sufficiently long time can be obtained as the solution of an integrodifferential equation with partial derivatives (equation of fragmentation). Quasi-stationary solutions of the equation of fragmentation are discussed for particular cases. The problem of the steady state is reduced to the solution of a transcendental equation. The results obtained show that analytical methods already offer a good theoretical understanding of the observed size distribution of the asteroids. They should be, therefore, a useful basis of carrying out numerical calculations.     

Solution of the Sitnikov Problem

A new analytic expression for the position of the infinitesimal body in the elliptic Sitnikov problem is presented. This solution is valid for small bounded oscillations in cases of moderate primary eccentricities. We first linearize the problem and obtain solution to this Hill's type equation. After that the lowest order nonlinear force is added to the problem. The final solution to the equation with nonlinear force included is obtained through first the use of a Courant and Snyder transformation followed by the Lindstedt–Poincaré perturbation method and again an application of Courant and Snyder transformation. The solution thus obtained is compared with existing solutions, and satisfactory agreement is found.     

Vibrational stability of stars in thermal imbalance: A solution in terms of asymptotic expansions

The solution of the partial differential equation describing the ‘non-isentropic’ oscillations of a star in thermal imbalance has been obtained in terms of asymptotic expansions up to the first order in the parameterII/t _s, whereII is the adiabatic pulsation period for the fundamental mode andt _s , a secular time scale of the order of the Kelvin-Helmholtz time. Use has been made of the zeroth order ‘isentopic’ solution derived in I. The solution obtained allows one to derive unambiguously a general integral expression for the coefficient of vibrational stability for arbitrary stellar models in thermal imbalance. The physical interpretation of this stability coefficient is discussed and its generality and its simplicity are stressed. Application to some simple analytic stellar models in homologous and nonhomologous contraction enables one to recover, in a more straightforward manner, results obtained by Coxet al. (1973). Aizenman and Cox (1974) and Davey (1974). Finally, we emphasize that the inclusion of the effects of thermal imbalance in the stability calculations of realistic evolutionary sequences of stellar models, not considered up to now by the other authors, is quite easy and straightforward with the simple formula derived here.     

Solvable cases of Szebehely's equation

Szebehely’s equation is a first order partial differential equation relating a given family of orbits f (x, y) = q traced by a unit mass material point, the total energy E=E(f), and the unknown potential V=V (x, y) which produces the family. Although linear in V, this equation cannot generally be solved. In this paper we develop the reasoning for finding several cases for which Szebehely’s equation can be solved by quadratures. This revised version was published online in July 2006 with corrections to the Cover Date.     

Generalization of Szebehely's equation

Szebehely's partial differential equation for the force functionU=U(x,y) which gives rise to a given family of planar orbitsf(x,y)=Constant is generalized to account for velocity-dependent potentials V^*=V^*(x,y, ). The new partial differential equation is quasi-linear and of the first order. An example is given and a comparison is made of the two equations.     

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.     

Solitary waves of vortices

We study the propagation of solitary waves of vortices within a spherical shell which constitutes the uppermost layer of a solid planet. This solid-liquid configuration rotates with constant angular velocity about an axis which is fixed with respect to the solid surface. The fluid within the shell is inviscid, incompressible, and of constant density. The motion imparted by the planetary rotation upon this fluid mass is governed by the Laplace tidal equation from which the potential of the extraplanetary forces has been deleted. Consistent with this ocean model, we establish that the stream function of a solitary wave of vortices must satisfy a third-order partial differential equation. We obtain solutions to this wave equation by imposing the condition that the vertical component of vorticity be functionally related to the stream function. We find that this dependence must necessarily be of the exponential type and that the solution to the wave equation then reduces to a quadrature depending on some arbitrary parameters. We prove that we can always choose the values of these parameters in order to approximate the integral in question by means of an analytic function: we reach a representation of the stream function of a solitary wave of vortices in terms of hyperbolic functions of time and position.This paper is dedicated to the memory of Professor Zdenek Kopal.     

Decoupling treatment of linear non-adiabatic non-radial pulsations of stars

A decoupling method is developed in this paper to deal with linear non-adiabatic non-radial pulsations of stars. The sixth order differential equation of linear pulsation is decomposed into a fourth order and a second order differential equations. The decoupling overcomes such difficulties encountered in the numerical solution as small domain and slowness of convergence and provides a natural guess solution needed in Henyey's method.     

Sur un modele explicitant les differents types morphologiques des galaxies

In the now classical Lindblad-Lin density-wave theory, the linearization of the collisionless Boltzmann equation is made by assuming the potential functionU expressed in the formU=U ₀ + \(\tilde U\) +... WhereU ₀ is the background axisymmetric potential and \(\tilde U<< U_0 \) . Then the corresponding density distribution is \(\rho = \rho _0 + \tilde \rho (\tilde \rho<< \rho _0 )\) and the linearized equation connecting \(\tilde U\) and the component \(\tilde f\) of the distribution function is given by $$\frac{{\partial \tilde f}}{{\partial t}} + \upsilon \frac{{\partial \tilde f}}{{\partial x}} - \frac{{\partial U_0 }}{{\partial x}} \cdot \frac{{\partial \tilde f}}{{\partial \upsilon }} = \frac{{\partial \tilde U}}{{\partial x}}\frac{{\partial f_0 }}{{\partial \upsilon }}.$$ One looks for spiral self-consistent solutions which also satisfy Poisson's equation $$\nabla ^2 \tilde U = 4\pi G\tilde \rho = 4\pi G\int {\tilde f d\upsilon .} $$ Lin and Shu (1964) have shown that such solutions exist in special cases. In the present work, we adopt anopposite proceeding. Poisson's equation contains two unknown quantities \(\tilde U\) and \(\tilde \rho \) . It could be completelysolved if a second independent equation connecting \(\tilde U\) and \(\tilde \rho \) was known. Such an equation is hopelesslyobtained by direct observational means; the only way is to postulate it in a mathematical form. In a previouswork, Louise (1981) has shown that Poisson's equation accounted for distances of planets in the solar system(following to the Titius-Bode's law revised by Balsano and Hughes (1979)) if the following relation wasassumed $$\rho ^2 = k\frac{{\tilde U}}{{r^2 }} (k = cte).$$ We now postulate again this relation in order to solve Poisson's equation. Then, $$\nabla ^2 \tilde U - \frac{{\alpha ^2 }}{{r^2 }}\tilde U = 0, (\alpha ^2 = 4\pi Gk).$$ The solution is found in a classical way to be of the form $$\tilde U = cte J_v (pr)e^{ - pz} e^{jn\theta } $$ wheren = integer,p =cte andJ _v (pr) = Bessel function with indexv (v ² =n ² + α²). By use of the Hankel function instead ofJ _v (pr) for large values ofr, the spiral structure is found to be given by $$\tilde U = cte e^{ - pz} e^{j[\Phi _v (r) + n\theta ]} , \Phi _v (r) = pr - \pi /2(v + \tfrac{1}{2}).$$ For small values ofr, \(\tilde U\) = 0: the center of a galaxy is not affected by the density wave which is onlyresponsible of the spiral structure. For various values ofp,n andv, other forms of galaxies can be taken into account: Ring, barred and spiral-barred shapes etc. In order to generalize previous calculations, we further postulateρ ₀ =kU ₀/r ², leading to Poisson'sequation which accounts for the disc population $$\nabla ^2 U_0 - \frac{{\alpha ^2 }}{{r^2 }}U_0 = 0.$$ AsU ₀ is assumed axisymmetrical, the obvious solution is of the form $$U_0 = \frac{{cte}}{{r^v }}e^{ - pz} , \rho _0 = \frac{{cte}}{{r^{2 + v} }}e^{ - pz} .$$ Finally, Poisson's equation is completely solvable under the assumptionρ =k(U/r ². The general solution,valid for both disc and spiral arm populations, becomes $$U = cte e^{ - pz} \left\{ {r^{ - v} + } \right.\left. {cte e^{j[\Phi _v (r) + n\theta ]} } \right\},$$ The density distribution along the O _z axis is supported by Burstein's (1979) observations.     

A linear solution of the equations of motion of an earth-orbiting satellite based on a Lie-series

Using Hill's variables, an analytical solution of a canonical system of six differential equations describing the motion of a satellite in the gravitational field of the earth is derived. The gravity field, expanded into spherical harmonics, has to be expressed as a function of the Hill variables. The intermediary is chosen to include the main secular terms. The first order solution retains the highly practical formal structure of Kaula's linear solution, but is valid for circular orbits and provides of course a spectral decomposition of radius vector and radial velocity. The resulting eccentricity functions are much simpler than the Hansen functions, since a series evaluation of the Kepler equation is avoided. The present solution may be extended to higher order solutions by Hori's perturbation method.     

Polytropic equilibrium figures of equal mass

A numerical method for the solution of the (astrophysical) potential problem is presented. The problem is formulated as a free boundary problem for a mildly nonlinear elliptic partial differential equation and the method is obtained by combining Newton-Raphson's procedure and two different types of discretization. The performance of the method is discussed.     

Gravitational clustering of galaxies in an expanding universe

We inquire the phenomena of clustering of galaxies in an expanding universe from a theoretical point of view on the basis of thermodynamics and correlation functions. The partial differential equation is developed both for the point mass and extended mass structures of a two-point correlation function by using thermodynamic equations in combination with the equation of state taking gravitational interaction between particles into consideration. The unique solution physically satisfies a set of boundary conditions for correlated systems and provides a new insight into the gravitational clustering problem.     

恒星线性非绝热非径向脉动问题的退耦化处理方法

本文发展了一种解恒星线性非绝热非径向脉动问题的退耦化方法。这个方法把非绝热非径向脉动问题的六阶线性微分方程,分解为由一个代数方程联系起来的一个四阶线性微分方程和另一个二阶线性微分方程进行数值求解。这样的一个退耦处理,有利于克服以前在数值解这类问题时常常遇到的收敛域小和收敛速度慢等困难,并且为数值解方程时所采用的Henyey方法提供了一个自然和方便的初始猜测解。