共查询到20条相似文献，搜索用时 62 毫秒
1.
This paper gives the results of a programme attempting to exploit ‘la seule bréche’ (Poincaré, 1892, p. 82) of non-integrable systems, namely to develop an approximate general solution for the three out of its four component-solutions of the planar restricted three-body problem. This is accomplished by computing a large number of families of ‘solutions précieuses’ (periodic solutions) covering densely the space of initial conditions of this problem. More specifically, we calculated numerically and only for μ = 0.4, all families of symmetric periodic solutions (1st component of the general solution) existing in the domain D:(x
0 ∊ [−2,2],C ∊ [−2,5]) of the (x
0, C) space and consisting of symmetric solutions re-entering after 1 up to 50 revolutions (see graph in Fig. 4). Then we tested the parts of the domain D that is void of such families and established that they belong to the category of escape motions (2nd component of the general solution). The approximation of the 3rd component (asymmetric solutions) we shall present in a future publication. The 4th component of the general solution of the problem, namely the one consisting of the bounded non-periodic solutions, is considered as approximated by those of the 1st or the 2nd component on account of the `Last Geometric Theorem of Poincaré' (Birkhoff, 1913). The results obtained provoked interest to repeat the same work inside the larger closed domain D:(x
0 ∊ [−6,2], C ∊ [−5,5]) and the results are presented in Fig. 15. A test run of the programme developed led to reproduction of the results presented by Hénon (1965) with better accuracy and many additional families not included in the sited paper. Pointer directions construed from the main body of results led to the definition of useful concepts of the basic family of order
n, n = 1, 2,… and the completeness criterion of the solution inside a compact sub-domain of the (x
0, C) space. The same results inspired the ‘partition theorem’, which conjectures the possibility of partitioning an initial conditions domain D into a finite set of sub-domains D
i that fulfill the completeness criterion and allow complete approximation of the general solution of this problem by computing a relatively small number of family curves. The numerical results of this project include a large number of families that were computed in detail covering their natural termination, the morphology, and stability of their member solutions. Zooming into sub-domains of D permitted clear presentation of the families of symmetric solutions contained in them. Such zooming was made for various values of the parameter N, which defines the re-entrance revolutions number, which was selected to be from 50 to 500. The areas generating escape solutions have being investigated. In Appendix A we present families of symmetric solutions terminating at asymptotic solutions, and in Appendix B the morphology of large period symmetric solutions though examples of orbits that re-enter after from 8 to 500 revolutions. The paper concludes that approximations of the general solution of the planar restricted problem is possible and presents such approximations, only for some sub-domains that fulfill the completeness criterion, on the basis of sufficiently large number of families. 相似文献
2.
The general solution of the Henon–Heiles system is approximated inside a domain of the (x, C) of initial conditions (C is the energy constant). The method applied is that described by Poincaré as ‘the only “crack” permitting penetration into
the non-integrable problems’ and involves calculation of a dense set of families of periodic solutions that covers the solution
space of the problem. In the case of the Henon–Heiles potential we calculated the families of periodic solutions that re-enter
after 1–108 oscillations. The density of the set of such families is defined by a pre-assigned parameter ε (Poincaré parameter),
which ascertains that at least one periodic solution is computed and available within a distance ε from any point of the domain
(x, C) for which the approximate general solution computed. The approximate general solution presented here corresponds to ε =
0.07. The same solution is further improved by “zooming” into four square sub-domain of (x, C), i.e. by computing sufficient number of families that reduce the density parameter to ε = 0.003. Further zooming to reduce
the density parameter, say to ε = 10−6, or even smaller, although easily performable in both areas occupied by stable as well as unstable solutions, was found unnecessary.
The stability of all members of each and all families computed was calculated and presented in this paper for both the large
solution domain and for the sub-domains. The correspondence between areas of the approximate general solution occupied by
stable periodic solutions and Poincaré sections with well-aligned section points and also correspondence between areas occupied
by unstable solutions and Poincaré sections with randomly scattered section points is shown by calculating such sections.
All calculations were performed using the Runge-Kutta (R-K) 8th order direct integration method and the large output received,
consisting of many thousands of families is saved as “Atlas of the General Solution of the Henon–Heiles Problem,” including
their stability and is available at request. It is concluded that approximation of the general solution of this system is
straightforward and that the chaotic character of its Poincaré sections imposes no limitations or difficulties. 相似文献
3.
Two basic problems of dynamics, one of which was tackled in the extensive work of Z. Kopal (see e.g. Kopal, 1978, Dynamics of Close Binary Systems, D. Reidel Publication, Dordrecht, Holland.), are presented with their approximate general solutions. The ‘penetration’ into
the space of solution of these non-integrable autonomous and conservative systems is achieved by application of ‘The Last Geometric Theorem of Poincaré’ (Birkhoff, 1913, Am. Math. Soc. (rev. edn. 1966)) and the calculation of sub-sets of ‘solutions précieuses’ that are covering densely the spaces of all solutions
(non-periodic and periodic) of these problems. The treated problems are: 1. The two-dimensional Duffing problem, 2. The restricted
problem around the Roche limit. The approximate general solutions are developed by applying known techniques by means of which
all solutions re-entering after one, two, three, etc, revolutions are, first, located and then calculated with precision.
The properties of these general solutions, such as the morphology of their constituent periodic solutions and their stability
for both problems are discussed. Calculations of Poincaré sections verify the presence of chaos, but this does not bear on
the computability of the general solutions of the problems treated. The procedure applied seems efficient and sufficient for
developing approximate general solutions of conservative and autonomous dynamical systems that fulfil the PoincaréBirkhoff
theorems. The same procedure does not apply to the sub-set of unbounded solutions of these problems. 相似文献
4.
This paper presents the approximate general solution of the triple well, double oscillator non-linear dynamical system. This
system is non-integrable and the approximate general solution is calculated by application of the Last Geometric Theorem of
Poincaré (Birkhoff, 1913, 1925). The original problem, known as the Duffing one, is a 1 degree of freedom system that, besides
the conservative force component, includes dumping and external forcing terms (see details in the web site: http://www.uncwil.edu/people/hermanr/chaos/ted/chaos.html).
The problem considered here is a 2 degree of freedom, autonomous and conservative one, without dumping, and of axisymmetric
potential. The space of permissible motions is scanned for identification of all solutions re-entering after from one to nine
oscillations and the precise families of periodic solutions are computed, including their stability parameter, covering all
cases with periods T corresponding to 4osc/T. Seven sub-domains of the space of solutions were investigated in detail by zooming, an operation that proved the possibility
to advance the accuracy of the approximate general solution to the level permitted by the integration routine. The approximation
of the general solution, although impressive, provides clear evidence of the complexity of the problem and the need to proceed
to larger period families. Nevertheless, it allows prediction of the areas where chaos and order regions in the Poincaré surfaces
of section are to be expected. Examples of such surfaces of sections, as well as of types of closed solutions, are given.
Two peculiar points of the space of solutions were identified as crossing, or source points from which infinite families of
periodic solutions emanate. The morphology and stability of solutions of the problem are studied and discussed. 相似文献
5.
The vertical stability character of the families of short and long period solutions around the triangular equilibrium points
of the restricted three-body problem is examined. For three values of the mass parameter less than equal to the critical value
of Routh (μ
R
) i.e. for μ = 0.000953875 (Sun-Jupiter), μ = 0.01215 (Earth-Moon) and μ = μ
R
= 0.038521, it is found that all such solutions are vertically stable. For μ > (μ
R
) vertical stability is studied for a number of ‘limiting’ orbits extended to μ = 0.45. The last limiting orbit computed by
Deprit for μ = 0.044 is continued to a family of periodic orbits into which the well known families of long and short period
solutions merge. The stability characteristics of this family are also studied. 相似文献
6.
The galactic dynamical system expressed by a third-order axisymmetric polynomial potential is investigated numerically by computing periodic solutions. We define as Sthe compact set of initial conditions generating bounded motions, and as S p , with S p ? S, the countable set of all initial conditions generating periodic solutions. Then, we consider the subsets S s p and S a p of S p , where S s p ∪ S a p = S p , S s p S a p = Ø, the first of which corresponds to symmetric periodic solutions, and the second to asymmetric solutions. Then, we approximate the set S s p , leaving treatment of the set S a p of asymmetric solutions for a future publication. The set S s p is known to be dense in S (‘Last Geometric Theorem of Poincar;’, Birkhoff, 1913). Using a computer programme capable to locate all elements of the set S s p that generate symmetric periodic solutions that re-enter after intersecting the axis of symmetry from 1 to ntimes. The results of the approximation of S s p in the total domain and in the sample sub-domains of zooming, we present in graphical form as family curves in the (x, C) plane. The solutions located with the largest periods re-enter after 440 galaxy revolutions while the families calculated fully (initial conditions, period, energy, stability co-efficient) include solutions that re-enter after 340 galaxy revolutions. To advance further the approximation of the set S s p thus obtained, we applied the same procedure inside eight sub-domains of the domain Sinto which we ‘zoomed’ through selection of finer search steps and double maximum periods. The family curves thus calculated presented in the (x, C) plane do not intersect anywhere in some sub-domains and their pattern resembles that of laminar flow. In other sub-domains, however, we found family curves from which branching families emanate. The concepts of completeand non-completeapproximation of S s p in sub-domains of laminar and sub-domains with branching family curves, respectively, is introduced. Also, the concept of basic family of order1, 2, ..., n, are defined. The morphology of individual periodic solutions of all families is investigated, and the types of envelopes found are described. The approximate set S s p was also checked by computing Poincar; sections for energy values corresponding to the mean energy range of the eight sub-domains of zooming mentioned above. These sections show that most parts of the compact domain in Sgenerating non-periodic but bounded solutions correspond to with well-shaped tori that intersect the x-axis, a fact that implies that dominant to exclusive type of periodic solutions are the symmetric ones with two normal crossings of this axis. The presence of non-symmetric periodic solutions as well as of chaotic regions is encountered. All calculations reported here were performed using the variable step R-K 8th-order direct integration and setting the allowable energy variation Δ C= |C start? C end| < 10?13. The output, consisting of many thousands of families and their properties (initial conditions, morphology, stability, etc.), is stored in a directory entitled ‘Atlas of the Symmetric Periodic Solution of the Galactic Motion Problem’. 相似文献
7.
P. S. Soulis K. E. Papadakis T. Bountis 《Celestial Mechanics and Dynamical Astronomy》2008,100(4):251-266
We study the existence, linear stability and bifurcations of what we call the Sitnikov family of straight line periodic orbits
in the case of the restricted four-body problem, where the three equal mass primary bodies are rotating on a circle and the
fourth (small body) is moving in the direction vertical to the center mass of the other three. In contrast to the restricted
three-body Sitnikov problem, where the Sitnikov family has infinitely many stability intervals (hence infinitely many Sitnikov critical orbits), as the “family parameter” ż0 varies within a finite interval (while z
0 tends to infinity), in the four-body problem this family has only one stability interval and only twelve 3-dimensional (3D) families of symmetric periodic orbits exist which bifurcate from twelve
corresponding critical Sitnikov periodic orbits. We also calculate the evolution of the characteristic curves of these 3D
branch-families and determine their stability. More importantly, we study the phase space dynamics in the vicinity of these
orbits in two ways: First, we use the SALI index to investigate the extent of bounded motion of the small particle off the
z-axis along its interval of stable Sitnikov orbits, and secondly, through suitably chosen Poincaré maps, we chart the motion
near one of the 3D families of plane-symmetric periodic orbits. Our study reveals in both cases a fascinating structure of
ordered motion surrounded by “sticky” and chaotic orbits as well as orbits which rapidly escape to infinity. 相似文献
8.
For the equation describing plane oscillations and rotations of a satellite, we consider families of symmetric generalized
periodic solutions with integral rotation number p. We give new confirmations of the hypothesis: there are only four classes of these families with topologically different
structures, namely, the classes of families of periodic solutions with p≥ 1, p= 0, p=−1, and p≤−2. Besides, we demonstrate that the vertices of cusps of these families are placed on some analytical curves, and the same
is true for the multiple intersections of these families with other families. 相似文献
9.
This paper presents the procedure of a computational scheme leading to approximate general solution of the axi-symmetric,2-degrees
of freedom dynamical systems. Also the results of application of this scheme in two such systems of the non-linear double
oscillator with third and fifth order potentials in position variables. Their approximate general solution is constructed
by computing a dense set of families of periodic solutions and their presentation is made through plots of initial conditions.
The accuracy of the approximate general solution is defined by two error parameters, one giving a measure of the accuracy
of the integration and calculation of periodic solutions procedure, and the second the density in the initial conditions space
of the periodic solutions calculated. Due to the need to compute families of periodic solutions of large periods the numerical
integrations were carried out using the eighth order, variable step, R-K algorithm, which secured for almost all results presented
here conservation of the energy constant between 10-9 and 10-12 for single runs of any and all solutions. The accuracy of the approximate general solution is controlled by increasing the
number of family curves and also by `zooming' into parts of the space of initial conditions. All families of periodic solutions
were checked for their stability. The computation of such families within areas of `deterministic chaos' did not encounter
any difficulty other than poorer precision. Furthermore, on the basis of the stability study of the computed families, the
boundaries of areas of `order' and `chaos' were approximately defined. On the basis of these results it is concluded that
investigations in thePoincaré sections have to disclose 3 distinct types of areas of `order' and 2 distinct types of areas
of `chaos'. Verification of the `order'/`chaos' boundary calculation was made by working out several Poincaré surfaces of
sections.
This revised version was published online in July 2006 with corrections to the Cover Date. 相似文献
10.
The predictor-corrector method is described for numerically extending with respect to the parameters of the periodic solutions of a Lagrangian system, including recurrent solutions. The orbital stability in linear approximation is investigated simultaneously with its construction.The method is applied to the investigation of periodic motions, generated from Lagrangian solutions of the circular restricted three body problem. Small short-period motions are extended in the plane problem with respect to the parameters h, µ (h = energy constant, µ = mass ratio of the two doninant gravitators); small vertical oscillations are extended in the three-dimensional problem with respect to the parameters h, µ. For both problems in parameter's plane h, µ domaines of existince and stability of derived periodic motions are constructed, resonance curves of third and fourth orders are distinguished. 相似文献
11.
Ramadan M. Ali 《Astrophysics and Space Science》2007,310(3-4):201-204
Cylindrically symmetric perfect fluid solutions are derived for the Levi-Civita metric. The pressure P is finite. The matter density is greater than the stresses in the material. The solutions are inside cylinders of bounded
radius at which the pressure vanishes. The range of σ, for which the sources have been matched to the Levi-Civeta metric is ∞>σ>0. The solutions are regular and satisfy energy conditions 相似文献
12.
13.
Mercè Ollé Joan R. Pacha Jordi Villanueva 《Celestial Mechanics and Dynamical Astronomy》2004,90(1-2):87-107
The paper deals with different kinds of invariant motions (periodic orbits, 2D and 3D invariant tori and invariant manifolds of periodic orbits) in order to analyze the Hamiltonian direct Hopf bifurcation that
takes place close to the Lyapunov vertical family of periodic orbits of the triangular equilibrium point L4 in the 3D restricted three-body problem (RTBP) for the mass parameter, μ greater than (and close to) μR (Routh’s mass parameter). Consequences of such bifurcation, concerning the confinement of the motion close to the hyperbolic
orbits and the 3D nearby tori are also described. 相似文献
14.
For the n-centre problem of one particle moving in the potential of attracting centres of small mass fixed in an arbitrary smooth potential
and magnetic field, we prove the existence of periodic and chaotic trajectories shadowing sequences of collision orbits. In
particular, we obtain large subshifts of solutions of this type for the circular restricted 3-body problem of celestial mechanics.
Poincaré had conjectured existence of the periodic ones and given them the name ‘second species solutions’.
This revised version was published online in July 2006 with corrections to the Cover Date. 相似文献
15.
Poincaré surface of section technique is used to study the evolution of a family ‘f’ of simply symmetric retrograde periodic orbits around the smaller primary in the framework of restricted three-body problem for a number of systems, actual and hypothetical, with mass ratio varying from 10−7 to 0.015. It is found that as the mass ratio decreases the region of phase space containing the two separatrices shrinks in size and moves closer to the smaller primary. Also the corresponding value of Jacobi constant tends towards 3. 相似文献
16.
E. A. Perdios 《Celestial Mechanics and Dynamical Astronomy》2007,99(2):85-104
This paper deals with the Sitnikov family of straight-line motions of the circular restricted three-body problem, viewed as
generator of families of three-dimensional periodic orbits. We study the linear stability of the family, determine several
new critical orbits at which families of three dimensional periodic orbits of the same or double period bifurcate and present
an extensive numerical exploration of the bifurcating families. In the case of the same period bifurcations, 44 families are
determined. All these families are computed for equal as well as for nearly equal primaries (μ = 0.5, μ = 0.4995). Some of the bifurcating families are determined for all values of the mass parameter μ for which they exist. Examples of families of three dimensional periodic orbits bifurcating from the Sitnikov family at double
period bifurcations are also given. These are the only families of three-dimensional periodic orbits presented in the paper
which do not terminate with coplanar orbits and some of them contain stable parts. By contrast, all families bifurcating at
single-period bifurcations consist entirely of unstable orbits and terminate with coplanar orbits. 相似文献
17.
It is shown (1) that the coefficients Ai of the limb darkening functions I(μ)/Icenter = P5 (μ) = ∑Ai μi (i = 0... 5; μ = cos ϑ), which had been published by Neckel and Labs (Solar Phys. 153, 91, 1994), can well be approximated by analytical functions of wavelength λ, and (2) that at first sight purely formal extrapolation
of the functions P5(μ) to the very limb (μ = 0.0) is not meaningless: in combination with absolute intensities for the disk center these functions
yield ‘limb intensities’ which all correspond to almost the same ‘limb temperature’, Tlimb≈4746 K. Together these results lead to ‘reference functions’ which can quickly yield rather reliable values of the Sun's
continuum intensities, for any values of μ and λ. 相似文献
18.
D. Viswanath 《Celestial Mechanics and Dynamical Astronomy》2006,94(2):213-235
The restricted three-body problem describes the motion of a massless particle under the influence of two primaries of masses
1− μ and μ that circle each other with period equal to 2π. For small μ, a resonant periodic motion of the massless particle
in the rotating frame can be described by relatively prime integers p and q, if its period around the heavier primary is approximately 2π p/q, and by its approximate eccentricity e. We give a method for the formal development of the stable and unstable manifolds associated with these resonant motions.
We prove the validity of this formal development and the existence of homoclinic points in the resonant region. In the study
of the Kirkwood gaps in the asteroid belt, the separatrices of the averaged equations of the restricted three-body problem
are commonly used to derive analytical approximations to the boundaries of the resonances. We use the unaveraged equations
to find values of asteroid eccentricity below which these approximations will not hold for the Kirkwood gaps with q/p equal to 2/1, 7/3, 5/2, 3/1, and 4/1. Another application is to the existence of asymmetric librations in the exterior resonances.
We give values of asteroid eccentricity below which asymmetric librations will not exist for the 1/7, 1/6, 1/5, 1/4, 1/3,
and 1/2 resonances for any μ however small. But if the eccentricity exceeds these thresholds, asymmetric librations will exist
for μ small enough in the unaveraged restricted three-body problem. 相似文献
19.
A. P. Markeev 《Astronomy Letters》2001,27(7):475-479
The centers of the gaps observed in the asteroid belt are displaced toward Jupiter from their positions that correspond to the exact commensurability between the mean motions of an asteroid and Jupiter. Using the current theory of stability and nonlinear oscillations of Hamiltonian systems, we point out the dynamical causes of this asymmetry. Our analysis is performed in terms of the plane circular restricted three-body problem. The orbits that correspond to Poincaré periodic solutions of the first kind are taken as unperturbed asteroid orbits. 相似文献
20.
We consider the problem of 4 bodies of equal masses in R
3 for the Newtonian r−1 potential. We address the question of the absolute minima of the action integral among (anti)symmetric loops of class H
1 whose period is fixed. It is the simplest case for which the results of [4] (corrected in [5]) do not apply: the minima cannot
be the relative equilibria whose configuration is an absolute minimum of the potential among the configurations having a given
moment of inertia with respect to their center of mass. This is because the regular tetrahedron cannot have a relative equilibrium
motion in R
3 (see [2]). We show that the absolute minima of the action are not homographic motions. We also show that if we force the
configuration to admit a certain type of symmetry of order 4, the absolute minimum is a collisionless orbit whose configuration
‘hesitates’ between the central configuration of the square and the one of the tetrahedron. We call these orbits ‘hip-hop’.
A similar result holds in case of a symmetry of order 3 where the central configuration of the equilateral triangle with a
body at the center of mass replaces the square.
This revised version was published online in July 2006 with corrections to the Cover Date. 相似文献