Isovortical orbits of autonomous,conservative, two degree-of-freedom dynamical systems

For an autonomous, conservative, two degree-of-freedom dynamical system, vorticity (the curl of velocity) is constant along the orbit if the velocity field is divergence-free such that: $$u\left( {x, v} \right) - \psi _y , v\left( {x, y} \right) = - \psi _x .$$ Isovortical orbits in configuration space are level curves of a scalar autonomous function Ψ (x, v) satisfying a second-order, non-linear partial differential equation of the Monge-Ampere type: $$2\left( {\psi _{xx} \psi _{yy} - \psi _{xy}^2 } \right) + U_{xx} + U_{yy} = 0,$$ where U(x. y) is the autonomous potential function. The solution Soc the time variable is reduced to a quadrature following determinatio of Ψ. Self-similar solutions of the Monge-Ampere equation under Birkhoff's one-parameter transformation group are derived for homogeneous (power-law) potential functions. It is shown that Keplerian orbits belong to the class of planar isovortical flows.     

A new integration theory for conservative two degree-of-freedom systems

A two degree-of-freedom, conservative system is reduced to a single degree-of-freedom, kinematic system with Hamiltonian integral under the change of independent variable: $$dt = \zeta dt (\zeta = \upsilon _x - \upsilon _y )$$ where ζ is the curl (or vorticity) of the velocity field with cartesian inertial componentsu(x, y, t) andv(x, y, t). In the autonomous case whenu _t=v _t=0, orbits are globally represented by the level curves of an autonomous Hamiltonian functionH(x,y) satisfying a second-order quasilinear partial differential equation (Szebehely's Equation): $$2(H + U)\left( {H_{xx} H_y^2 - 2H_{xy} H_x H_y + H_{yy} H_x^2 } \right) + (H_x U_x + H_y U_y )\left( {H_x^2 + H_y^2 } \right) = 0$$ whereU(x, y) is the autonomous potential function. An inversion of dependent and independent variables reduces this equation to a second-order, ordinary differential equation for a function specifying the orbital curve. The true time variable is recovered by evaluating a quadrature. Fundamental differences exist between this approach and Hamilton-Jacobi theory.     

Periodic orbits and their bifurcations in a 3-D system

We study some simple periodic orbits and their bifurcations in the Hamiltonian . We give the forms of the orbits, the characteristics of the main families, and some existence diagrams and stability diagrams. The existence diagram of the family 1a contains regions that are stable (S), simply unstable (U), doubly unstable (DU) and complex unstable (). In the regionsS andU there are lines of equal rotation numberm/n. Along these lines we have bifurcations of families of periodic orbits of multiplicityn. When these lines reach the boundary of the complex unstable region, they are tangent to it. Inside the region there are linesm/n, along which the orbits 1a, describedn-times, are doubly unstable; however, along these lines there are no bifurcations ofn-ple periodic orbits. The families bifurcating from 1a exist only in certain regions of the parameter space (, ). The limiting lines of these regions join at particular points representing collisions of bifurcations. These collisions of bifurcations produce a nonuniqueness of the various families of periodic orbits. The complicated structure of the various bifurcations can be understood by constructing appropriate stability diagrams.     

Isovortical orbits in uniformly rotating coordinates

For the conservative, two degree-of-freedom system with autonomous potential functionV(x,y) in rotating coordinates; $$\dot u - 2n\upsilon = V_x , \dot \upsilon + 2nu = V_y $$ , vorticity (v _x -u _y ) is constant along the orbit when the relative velocity field is divergence-free such that: $$u(x,y,t) = \psi _y , \upsilon (x,y,t) = - \psi _x $$ . Unlike isoenergetic reduction using the Jacobi, integral and eliminating the time,non-singular reduction from fourth to second-order occurs when (u,v) are determined explicitly as functions of their arguments by solving for ψ (x, y, t). The orbit function ψ satisfies a second-order, non-linear partial differential equation of the Monge Ampere type: $$2(\psi _{xx} \psi _{yy} - \psi _{xy}^2 ) - 2(\psi _{xx} + \psi _{yy} ) + V_{xx} + V_{yy} = 0$$ . Isovortical orbits in the rotating frame arenot level curves of ψ because it contains time explicitly due to coriolis effects. Rather, (x, y) coordinates along the orbit are obtained, from (u, v) either by numerical integration of the kinematic equations, or by partial differentiation of the Legendre transform ? of ψ. In the latter case, ? is shown to satisfy a non-linear, second-order partial differential equation in three independent variables, derived from the Monge-Ampere Equation. Complete reduction to quadrature is possible when space-time symmetries exist, as in the case of central force motion.     

Periodic orbits around an oblate spheroid

The general properties of certain differential systems are used to prove the existence of periodic orbits for a particle around an oblate spheroid.In a fixed frame, there are periodic orbits only fori=0 andi near /2. Furthermore, the generating orbits are circles.In a rotating frame, there are three families of orbits: first a family of periodic orbits in the vicinity of the critical inclination; secondly a family of periodic orbits in the equatorial plane with 0<e<1; thirdly a family of periodic orbits for any value of the inclination ife=0.     

Three-dimensional periodic motion of three finite masses around collinear equilibrium configurations

Three-dimensional periodic motions of three bodies are shown to exist in the infinitesimal neighbourhood of their collinear equilibrium configurations. These configurations and some characteristic quantities of the emanating three-dimensional periodic orbits are given for many values of the two mass parameters, =m ₂/(m ₁+m ₂) andm ₃, of the general three-body problem, under the assumption that the straight line containing the bodies at equilibrium rotates with unit angular velocity. The analysis of the small periodic orbits near the equilibrium configurations is carried out to second-order terms in the small quantities describing the deviation from plane motion but the analytical solution obtained for the horizontal components of the state vector is valid to third-order terms in those quantities. The families of three-dimensional periodic orbits emanating from two of the collinear equilibrium configurations are continued numerically to large orbits. These families are found to terminate at large vertical-critical orbits of the familym of retrograde periodic orbits ofm ₃ around the primariesm ₁ andm ₂. The series of these termination orbits, formed when the value ofm ₃ varies, are also given. The three-dimensional orbits are computed form ₃=0.1.     

The stability parameters of three-body periodic orbits

Formulae containing the elements of the variational matrix are obtained which determine the linear iso-energetic stability parameters of periodic orbits of the general three-body problem. This requires the numerical integration of the variational equations but produces the stability parameters with the effective accuracy of the numerical integration. The procedure is applied for the determination of horizontally critical orbits among the members of sets of vertical-critical periodic orbits of the threebody problem. These critical-critical orbits have special importance as they delimit the regions in the space of initial conditions which correspond to possibly stable three-dimensional periodic motion of low inclination.     

The Third Integral in a Self-consistent Galactic Model

We apply the theory of the third integral to a self-consistent galactic model, generated by the collapse of a N-body system. The final configuration after the collapse is a stationary triaxial system, that represents an almost prolate non-rotating elliptical galaxy with its longest axis in the z-direction. This system is represented by an axisymmetric potential V plus a small triaxial perturbation V ₁. The orbits in the potential V are of three types: box orbits, tube orbits (corresponding to various resonances), and chaotic orbits.The intersections of the box and tube orbits by a Poincaré surface of section z=0 are closed invariant curves. The main tube orbits are like ellipses and form an island of stability on the (R,R) plane.We calculated the third integral I in the potential V for the general non-resonant case and for various resonant cases. The agreement between the invariant curves of the orbits and the level curves of the third integral is good for the box and tube orbits, if we truncate the third integral at an appropriate level. As expected the third integral fails in the case of chaotic orbits. The most important result is the form of the number density F on the Poincaré surface of section. This function decreases exponentially outwards for the box orbits, like Fexp(–bI), while it is constant, as expected, for the chaotic orbits. In the case of the island of the main tube orbits it has a minimum at the center of the island. This can be explained by the form of the near elliptical orbits that are elongated along R, thus they fail to support a self-consistent galaxy, which is elongated along the z-axis.     

Analytic continuation of periodic orbits from equilibria of two-degree-of-freedom systems in rotating coordinates

A new theory is formulated for the analytic continuation of periodic (and aperiodic) orbits from equilibrium solutions of a two-degree-of-freedom dynamical system in rotating coordinates:% MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiiYdd9qrFfea0dXdf9vqai-hEir8Ve% ea0de9qq-hbrpepeea0db9q8as0-LqLs-Jirpepeea0-as0Fb9pgea% 0lrP0xe9Fve9Fve9qapdbaqaaeGacaGaaiaabeqaamaabaabcaGcba% acbiGab8xDayaacaGaa8xlaiaa-jdacaWFUbGaeqyXduNaa8xpaiaa% -zfadaWgaaWcbaGccaWF4baaleqaaOGaaiilaiqbew8a1zaacaGaey% 4kaSIaaGOmaiaad6gacaWG1bGaeyypa0Jaa8NvamaaBaaaleaakiaa% -LhaaSqabaGccaGGSaGabmiEayaacaGaeyypa0JaamyDaiaacYcace% WG5bGbaiaacqGH9aqpcqaHfpqDaaa!54CD!\[\dot u - 2n\upsilon = V_x ,\dot \upsilon + 2nu = V_y ,\dot x = u,\dot y = \upsilon \]Away from resonance, a family of nonlinear, normal-mode orbits defines an autonomous velocity field u(x, y), u(x, y) represented by convergent algebraic-series expansions in the two position variables. This approach is useful for determining the global structure of solution curves and nonlinear stability of normal modes using Liapunov's direct method. At resonance, the series coefficients are time dependent because stationary modes are incompatible with the equations of motion. By eliminating small divisors, explicit time dependence provides a natural transition from non-resonance to resonance cases within the same theory.     

A theory of libration

It is shown here that many problems of libration in celestial mechanics can be reduced to a perturbation of anintermediary defined by the Hamiltonian $$F = B\left( y \right) + 2\mu ^2 A\left( y \right)f\left( x \right).$$ This generalization of the Ideal Resonance Problem, with a periodic functionf(x) replacing sin² x, is solved here toO(μ ²) by an algorithm that is essentially the same as the one used in the original formulation. The solution is of the formx=x(u), u=u(t), y=y(x), with the functionx(u) commonly involving the inversion of a hyperelliptic integralu(x), evaluated by quadrature. Libration may be simple or multiple, depending on the nature of the functionf(x) and on the initial conditions. Double libration is illustrated here by the horseshoe-shaped orbits enclosing two libration centers.     

Three-dimensional periodic oscillations generating from plane periodic ones around the collinear Lagrangian points

The three families of three-dimensional periodic oscillations which include the infinitesimal periodic oscillations about the Lagrangian equilibrium pointsL ₁,L ₂ andL ₃ are computed for the value =0.00095 (Sun-Jupiter case) of the mass parameter. From the first two vertically critical (|a _v |=1) members of the familiesa, b andc, six families of periodic orbits in three dimensions are found to bifurcate. These families are presented here together with their stability characteristics. The orbits of the nine families computed are of all types of symmetryA, B andC. Finally, examples of bifurcations between families of three-dimensional periodic solutions of different type of symmetry are given.     

Motion Around The Triangular Equilibrium Points Of The Restricted Three-Body Problem Under Angular Velocity Variation

We study numerically the asymmetric periodic orbits which emanate from the triangular equilibrium points of the restricted three-body problem under the assumption that the angular velocity ω varies and for the Sun–Jupiter mass distribution. The symmetric periodic orbits emanating from the collinear Lagrangian point L ₃, which are related to them, are also examined. The analytic determination of the initial conditions of the long- and short-period Trojan families around the equilibrium points, is given. The corresponding families were examined, for a combination of the mass ratio and the angular velocity (case of equal eigenfrequencies), and also for the critical value ω = 2

Numerical determination of families of three-dimensional double-symmetric periodic orbits in the restricted three-body problem,II

New families of three-dimensional double-symmetric periodic orbits are determined numerically in the Sun-Jupiter case of the restricted three-body problem. These families bifurcate from the vertical-critical orbits ( _v = – 1, b _v – 0) of the basic plane familiesi, g ₁, g₂, c andI. Further, the predictor-corrector procedure employed to reveal these families has been described and interesting numerical results have been pointed out. Also, computer plots of the orbits of these families have been shown in conical projections.     

Escape regions of a quartic potential

We investigate the escape regions of a quartic potential and the main types of irregular periodic orbits. Because of the symmetry of the model the zero velocity curve consists of four summetric arcs forming four open channels around the lines y = ± x through which an orbit can escape. Four unstable Lyapunov periodic orbits bridge these openings.We have found an infinite sequence of families of periodic orbits which is the outer boundary of one of the escape regions and several infinite sequences of periodic orbits inside this region that tend to homoclinic and heteroclinic orbits. Some of these sequences of periodic orbits tend to homoclinic orbits starting perpendicularly and ending asymptotically at the x-axis. The other sequences tend to heteroclinic orbits which intersect the x-axis perpendicularly for x > 0 and make infinite oscillations almost parallel to each of the two Lyapunov orbits which correspond to x > 0 or x < 0.     

Oscillatory hydromagnetic flow through a porous medium in the presence of free convection and mass transfer flow

Unsteady two-dimensional hydromagnetic free convection and mass transfer flow of an electrically-conducting viscous-incompressible fluid, through a highly porous medium bounded by a vertical plane surface of constant temperature is considered. The free-stream velocity of the fluid vibrates about a mean constant value and the surface absorbs the fluid with constant velocity. Expressions for the velocity, temperature, concentration are obtained. Effects of Gr (Grashof number), Gm (modified Grashof number),K (permeability of the porous medium), (frequency parameter), andM (magnetic parameter) upon the velocity field are discussed.     

A grid search for families of periodic orbits in the restricted problem of three bodies

A method is described for the numerical determination of families of periodic orbits in the planar restricted problem of three bodies. The families are sought in their representation as curves in a two-dimensional space of parameters. A grid search is applied to the study of the evolution of satellite motion when the mass parameter is varied. Only that part of the space of parameters is investigated for which one of them, the relative energy constant, takes values larger than that corresponding to the inner Lagrangian pointL ₂. Critical values of the mass parameter are determined for which new families of simple or double periodic orbits appear inside the closed ovals of zero velocity.     

Newtonian and relativistic periodic orbits around two fixed black holes

We compare families of simple periodic orbits of test particles in the Newtonian and relativistic problems of two fixed centers (black holes). The Newtonian problem is integrable, while the relativistic problem is highly non-integrable.The orbits are calculated on the meridian plane through the fixed centersM ₁ (atz=+1) andM ₂ (atz=–1) for energies smaller than the escape energyE=1. We use prolate spheroidal coordinates (, , =const) and also the variables =cosh and =–cos . The orbits are inside a curve of zero velocity (CZV). The Newtonian orbits are also limited by an ellipse and a hyperbola, or by two eillipses. There are 3 main types of periodic orbits (1) elliptic type (around both centers), (2) hyperbolic-type, and (3) resonant-type.The elliptic type orbits are stable in the Newtonian case and both stable and unstable in the relativistic case. From the stable orbits bifurcate double period orbits both symmetric and asymmetric with respect to thez-axis. There are also higher order bifurcations. The hyperbolic-type orbits are unstable. The Newtonian resonant orbits are defined by the ratiot _µ/t =n/m of oscillations along and during one period, and they are all marginally unstable. The corresponding relativistic orbits are stable, or unstable. The main families are figure eight orbits aroundM ₁, or aroundM ₂ (3/1 orbits); gamma, or inverse gamma orbits (4/2); higher resonant families 5/1,7/1,...,8/2,12/2,...;, more complicated orbits, like 5/3, and bifurcations from the above orbits. Satellite orbits aroundM ₁, orM ₂, and their bifurcations (e.g. double period) exist in the relativistic case but not in the Newtonian case. The characteristics of the various families are quite different in the Newtonian and the relativistic cases. The sizes of the orbits and their stabilities are also quite different in general. In the Appendix we study the various types of straight line orbits and prove that some subcases introduced by Charlier (1902) are impossible.     

Bifurcations of Plane to 3d Periodic Orbits in the Restricted Three-Body Problem with Oblateness

We consider the bifurcation of 3D periodic orbits from the plane of motion of the primaries in the restricted three-body problem with oblateness. The simplest 3D periodic orbits branch-off at the plane periodic orbits of indifferent vertical stability. We describe briefly suitable numerical techniques and apply them to produce the first few such vertical-critical orbits of the basic families of periodic orbits of the problem, for varying mass parameter and fixed oblateness coefficent A₁ = 0.005, as well as for varying A₁ and fixed = 1/2. The horizontal stability of these orbits is also determined leading to predictions about the stability of the branching 3D orbits.     

Role of magnetic field on transient forced and free-convection flow past an infinite vertical porous plate through porous medium with heat source

Effects of magnetic field and permeability of the porous medium on unsteady forced and free-convection flow past an infinite vertical porous plate in presence of temperature-dependent heat source have been analysed. The Laplace transform method is used to obtain the expression for velocity field, skin friction, and leading edge effects. During the course of discussion, the effects ofM (magnetic parameter),S (heat source parameter), (suction parameter), andK (permeability of porous medium) on velocity field, skin friction, and leading edge effect have been extensively discussed.     

Three-dimensional periodic orbits about the triangular equilibrium points of the restricted problem of three bodies

The third-order parametric expansions given by Buck in 1920 for the three-dimensional periodic solutions about the triangular equilibrium points of the restricted Problem are improved by fourthorder terms. The corresponding family of periodic orbits, which are symmetrical w.r.t. the (x, y) plane, is computed numerically for =0.00095. It is found that the family emanating from L₄ terminates at the other triangular point L₅ while it bifurcates with the family of three-dimensional periodic orbits originating at the collinear equilibrium point L₃. This family consists of stable and unstable members. A second family of nonsymmetric three-dimensional periodic orbits is found to bifurcate from the previous one. It is also determined numerically until a collision orbit is encountered with the computations.