首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
The motion of Pluto is said to be chaotic in the sense that the maximum Lyapunov exponent is positive: the Lyapunov time (the inverse of the Lyapunov exponent) is about 20 million years. So far the longest integration up to now, over 845 million years (42 Lyapunov times), does not show any indication of a gross instability in the motion of Pluto. We carried out the numerical integration of Pluto over the age of the solar system (5.5 billion years ≈ 280 Lyapunov times). This integration also did not give any indication of chaotic evolution of Pluto. The divergences of Keplerian elements of a nearby trajectory at first grow linearly with the time and then start to increase exponentially. The exponential divergences stop at about 420 million years. The divergences in the semi-major axis and the mean anomaly ( equivalently the longitude and the distance) saturate. The divergences of the other four elements, the eccentricity, the inclination, the argument of perihelion, and the longitude of node still grow slowly after the stop of the exponential increase and finally saturate.  相似文献   

2.
The orbits of 13 Trojan asteroids have been calculated numerically in the model of the outer solar system for a time interval of 100 million years. For these asteroids Milani et al. (1997) determined Lyapunov times less than 100 000 years and introduced the notion "asteroids in stable chaotic motion". We studied the dynamical behavior of these Trojan asteroids (except the asteroid Thersites which escaped after 26 million years) within 11 time intervals - i.e. subintervals of the whole time - by means of: (1) a numerical frequency analysis (2) the root mean square (r.m.s.) of the orbital elements and (3) the proper elements. For each time interval we compared the root mean squares of the orbital elements (a, e and i) with the corresponding proper element. It turned out that the variations of the proper elements ep in the different time intervals are correlated with the corresponding r.m.s.(e); this is not the case for sin Ip with r.m.s.(i). This revised version was published online in July 2006 with corrections to the Cover Date.  相似文献   

3.
We report results from long term numerical integrations and analytical studies of particular orbits in the circular restricted three-body problem. These are mostly high-inclination trajectories in 1 : 1 resonance starting at or near the triangular Lagrangian L5 point. In some intervals of inclination these orbits have short Lyapunov times, from 100 to a few hundred periods, yet the osculating semi-major axis shows only relatively small fluctuations and there are no escapes from the 1 : 1 resonance. The eccentricity of these chaotic orbits varies in an erratic manner, some of the orbits becoming temporarily almost rectilinear. Similarly the inclination experiences large variations due to the conservation of the Jacobi constant. We studied such orbits for up to 109 periods in two cases and for 106 periods in all others, for inclinations varying from 0° to 180°. Thus our integrations extend from thousands to 10 million Lyapunov times without escapes of the massless particle. Since there are no zero-velocity curves restricting the motion this opens questions concerning the reason for the persistence of the 1 : 1 resonant motion. In the theory sections we consider the mechanism responsible for the observed behavior. We derive the averaged 1 : 1 resonance disturbing function, to second order in eccentricity, to calculate a critical inclination found in the numerical experiment, and examine motion close to this inclination. The cause of the chaos observed in the numerical experiments appears to be the emergence of saddle points in the averaged disturbing potential. We determine the location of several such saddle points in the (, ) plane, with being the mean longitude difference and the argument of pericentre. Some of the saddle points are illustrated with the aid of contour plots of the disturbing function. Motion close to these saddles is sensitive to initial conditions, thus causing chaos.This revised version was published online in October 2005 with corrections to the Cover Date.  相似文献   

4.
Qin  Zhang 《Solar physics》1998,178(2):423-431
Following the progression of nonlinear dynamical system theory, many authors have used varied methods to calculate the fractal dimension and the largest Lyapunov exponent 1 of the sunspot numbers and to evaluate the character of the chaotic attractor governing solar activity. These include the Grassberger–Procaccia algorithm, the technique provided by Wolf et al., and the nonlinear forecasting approach based on the method of distinguishing between chaos and measurement errors in time series described by Sugihara and May. In this paper, we use the Grassberger–Procaccia algorithm to estimate the other character of the chaotic attractor. This character is time scale of a transition from high-dimensional or stochastic at shorter times to a low-dimensional chaotic behavior at longer times. We find that the transitional time scale in the monthly mean sunspot numbers is about 8 yr; the low-dimensional chaotic behavior operates at time scales longer than about 8 yr and a high-dimensional or stochastic process operates at time scales shorter than about 8 yr.  相似文献   

5.
In this paper, we compare changes in the insolation at Pluto, corresponding to three epochs during the dynamical history of the planet: t = – 1, 0 and 0.5, where t is the time in millions of years A.D. The two extreme values of t coincide respectively with a maximum (126 ) and a minimum (102 ) value of the obliquity (). The other orbital elements i.e. the eccentricity (e) and the longitude of the perihelion ( p ) which affect solar radiation and which are apt to significant periodic changes are also calculated for the times under consideration. In a series of figures, the combined influence of the evolving dynamic parameters on the daily insolation and on the mean (summer, winter, annual) daily insolation is illustrated.  相似文献   

6.
The stability of the motion of a hypothetical planet in the binary system ?? Cen A?CB has been investigated. The analysis has been performed within the framework of a planar (restricted and full) three-body problem for the case of prograde orbits. Based on a representative set of initial data, we have obtained the Lyapunov spectra of the motion of a triple system with a single planet. Chaotic domains have been identified in the pericenter distance-eccentricity plane of initial conditions for the planet through a statistical analysis of the data obtained. We have studied the correspondence of these chaotic domains to the domains of initial conditions that lead to the planet??s encounter with one of the binary??s stars or to the escape of the planet from the system. We show that the stability criterion based on the maximum Lyapunov exponent gives a more clear-cut boundary of the instability domains than does the encounterescape criterion at the same integration time. The typical Lyapunov time of chaotic motion is ??500 yr for unstable outer orbits and ??60 yr for unstable inner ones. The domain of chaos expands significantly as the initial orbital eccentricity of the planet increases. The chaos-order boundary has a fractal structure due to the presence of orbital resonances.  相似文献   

7.
We investigate the potential importance of molecular cloud and stellar perturbations on the orbits of Pluto and more distant (hypothetical) planets up to 500 AU from the Sun. It is found that stellar and molecular cloud-core perturbations are of roughly equal importance. It also is found that the likelihood of substantial perturbations on Pluto is not insignificant, and that numerous substantial stellar and molecular cloud perturbations are likely to have influenced the orbits of any planets beyond 200 AU. These perturbations may contribute to a prevalence of moderate eccentricities and inclinations for planets beyond the orbit of Neptune, and may be a characteristic of distant planetary orbits in other solar systems. Given the recent discovery of chaotic behavior in Pluto's orbit (Sussman and Wisdom 1988), the effects of external perturbations on the long-term stability of Pluto's orbit warrant continued study.  相似文献   

8.
For the orbits with low to moderate inclination and eccentricity, in the asteroid main belt, the analytically computed proper elements are accurate to a level very close to the best result achievable by any analytical theory. This fundamental limitation results from the infinite web of resonances and because of the occurrence of chaotic motions. Still, there are some regions of the belt in which these proper elements are of degraded accuracy, thus preventing a reliable definition of asteroid families and detailed studies of the dynamical structure. We have used a different method to compute asteroid proper elements, following the approach introduced in the LONGSTOP project to describe the secular dynamics of the major outer planets. By applying purely numerical techniques, we produced so-called synthetic proper elements for a catalog of 10,256 asteroids with osculating semimajor axes between 2.5 and 4.0AU.The procedure consisted of simultaneous integration of asteroid and planetary orbits for 2Myr, with online filtering of the short-periodic perturbations. The output of the integration was spectrally resolved, and the principal harmonics (proper values) extracted from the time series. For each asteroid we have also tested the accuracy and stability in time of the proper elements, and estimated the maximum Lyapunov Characteristic Exponent to monitor the chaotic behaviors. This provided information on the reliability of the data for each orbit, in particular allowing to select 1,852 cases for an extended integration (10Myr) of the orbits showing instability. The results indicate that for more than half of the cases the proper elements have a time stability improved by more than a factor 3 with respect to the elements computed by the previous analytical theory. But of course there are also unstable cases for which the proper elements are less accurate and reliable, the extreme examples being 23 orbits exhibiting hyperbolic escape from the solar system. This form of escape from the asteroid belt could be responsible for a significant mass loss over the age of the solar system.  相似文献   

9.
Motions of asteroids in mean motion resonances with Jupiter are studied in three-dimensional space. Orbital changes of fictitious asteroids in the Kirkwood gaps are calculated by numerical integrations for 105 – 106 years. The main results are as follows: (1) There are various motions of resonant asteroids, and some of them are very complicated and chaotic and others are regular. (2) The eccentricity of some asteroids becomes very large, and the variation of the inclination is large while the eccentricity is large. (3) In the 3:1 resonance, there is a long periodic change in the variation of the inclination, when (7 : ) is a simple ratio (7: longitude of perihelion, : longitude of node). (4) In the 7:3 resonance, the variation of the inclination of some resonant asteroids is so large that prograde motion becomes retrograde. Some asteroids in the 7:3 resonance can collide with the Sun as well as with the inner planets.  相似文献   

10.
The chaotic orbital motion of Prometheus and Pandora, the 16th and 17th satellites of Saturn, is studied. Chaos in their orbital motion, as found by Goldreich & Rappaport and Renner & Sicardy, is due to interaction of resonances in the resonance multiplet corresponding to the 121:118 commensurability of the mean motions of the satellites. It is shown rigorously that the system moves in adiabatic regime. The Lyapunov time (the 'time horizon of predictability' of the motion) is calculated analytically and compared to the available numerical–experimental estimates. For this purpose, a method of analytical estimation of the maximum Lyapunov exponent in the perturbed pendulum model of non-linear resonance is applied. The method is based on the separatrix map theory. An analytical estimate of the width of the chaotic layer is made as well, based on the same theory. The ranges of chaotic diffusion in the mean motion are shown to be almost twice as big compared to previous estimates for both satellites.  相似文献   

11.
Three different numerical techniques are tested to determine the number of integrals of motion in dynamical systems with three degrees of freedom.
  1. The computation of the whole set of Lyapunov Characteristic Exponents (LCE).
  2. The triple sections in the configurations space.
  3. The Stine-Noid box-counting technique.
These methods are applied to a triple oscillator with coupling terms of the third order. Cases are found for which one effective integral besides the Hamiltonian subsists during a very long time. Such orbits display simultaneously chaotic and quasi-periodic motion, according to which coordinates are considered. As an application, the LCE procedure is applied to a triaxial elliptical galaxy model. Contrary to similar 2-dimensional systems, this 3-dimensional one presents noticeable zones in the phase space without any non-classical integral.  相似文献   

12.
A review is given of the determination of the long-term motion of Pluto. In particular, the discovery of the librational character of the two critical arguments is discussed. The stability of the motion of Pluto is shown to have been established when allknown gravitational forces are considered.Presented at the Symposium Star Catalogues, Positional Astronomy and Celestial Mechanics, held in honor of Paul Herget at the U.S. Naval Observatory, Washington, November 30, 1978.  相似文献   

13.
The chaotic orbital dynamics of the planet in the wide visual binary star system 16 Cyg is considered. The only planet in this system has a significant orbital eccentricity, e = 0.69. Previously, Holman et al. suggested the possibility of chaos in the orbital dynamics of the planet due to the proximity of 16 Cyg to the separatrix of the Lidov–Kozai resonance. We have calculated the Lyapunov characteristic exponents on the set of possible orbital parameters for the planet. In all cases, the dynamics of 16 Cyg is regular with a Lyapunov time of more than 30 000 yr. The dynamics is considered in detail for several possible models of the planetary orbit; the dependences of Lyapunov exponents on the time of their calculation and the time dependences of osculating orbital elements have been constructed. Phase space sections for the system dynamics near the Lidov–Kozai resonance have been constructed for all models. A chaotic behavior in the orbital motion of the planet in 16 Cyg is shown to be unlikely, because 16 Cyg in phase space is far from the separatrix of the Lidov–Kozai resonance at admissible orbital parameters, with the chaotic layer near the separatrix being very narrow.  相似文献   

14.
Five outer planets are numerically integrated over five million years in the Newtonian frame. The argument of Pluto's perihelion librates about 90 degrees with an amplitude of about 23 degrees. The period of the libration depends on the mass of Pluto: 4.0×106 years forM pluto=2.78×10–6 M sun and 3.8×106 years forM pluto=7.69×10–9 M sun, which is the newly determined mass. The motion of Neptune's perihelion is more sensitive to the mass of Pluto. ForM pluto=7.69×10–9 M sun, the perihelion of Neptune does circulate counter-clockwise and forM pluto=2.78×10–6 M sun, it does not circulate and the Neptune's eccentricity does not have a minimum. With the initial conditions which do not lie in the resonance region between Neptune and Pluto, a close approach between them takes place frequently and the orbit of Pluto becomes unstable and irregular.  相似文献   

15.
Jack Wisdom 《Icarus》1983,56(1):51-74
The sudden eccentricity increases discovered by J. Wisdom (Astron J.87, 577–593, 1982) are reproduced in numerical integrations of the planar-elliptic restricted three-body problem, verifying that this phenomenon is real. Maximum Lyapunov characteristic exponents for trajectories near the 31 commensurability are computed both with the mappings presented in Wisdom (1982) and by numerical integration of the planar-elliptic problem. In all cases the agreement is excellent, indicating that the mappings accurately reflect whether trajectories are chaotic or quasiperiodic. The mappings are used to trace out the chaotic zone near the 31 commensurability, both in the planar-elliptic problem and to a more limited extent in the three-dimensional elliptic problem. The outer boundary of the chaotic zone coincides with the boundary of the 31 Kirkwood gap in the actual distribution of asteroids within the errors of the asteroid orbital elements.  相似文献   

16.
The orbital dynamics of the single known planet in the binary star system HD 196885 has been considered. The Lyapunov characteristic exponents and Lyapunov time of the planetary system have been calculated for possible values of the planetary orbit parameters. It has been shown that the dynamics of the planetary system HD 196885 is regular with the Lyapunov time of more than 5 × 104 years (the orbital period of the planet is approximately 3.7 years), if the motion occurs at a distance from the separatrix of the Lidov–Kozai resonance. The values of the planet’s orbital inclination to the plane of the sky and longitude of the ascending node lie either within ranges 30° < i p < 90° and 30° < Ωp < 90°, or 90° < i p < 180° and 180° < Ωp < 300°.  相似文献   

17.
The dynamics of two families of minor inner solar system bodies that suffer frequent close encounters with the planets is analyzed. These families are: Jupiter family comets (JF comets) and Near Earth Asteroids (NEAs). The motion of these objects has been considered to be chaotic in a short time scale,and the close encounters are supposed to be the cause of the fast chaos. For a better understanding of the chaotic behavior we have computed Lyapunov Characteristic Exponents (LCEs) for all the observed members of both populations. LCEs are a quantitative measure of the exponential divergence of initially close orbits. We have observed that most members of the two families show a concentration of Lyapunov times (inverse of LCE) around 50–100yr. The concentration is more pronounced for JF comets than for NEAs, among which a lesser spread is observed for those that actually cross the Earth's orbit (mean perihelion distance q < 1.05 AU). It is also observed that a general correspondence exists between Lyapunov times and the time between consecutive encounters. A simple model is introduced to describe the basic characteristics of the dynamical evolution. This model considers an impulsive approach, where the particles evolve unperturbedly between encounters and suffer ‘kicks’ in semimajor axis at the encounters. It also reproduces successfully the short Lyapunov times observed in the numerical integrations and is able to estimate the dynamical lifetimes of comets during a stay in the Jupiter family in correspondence with previous estimates. It has been demonstrated with the model that the encounters with the largest effect on the exponential growth of the distance between initially nearby orbits are neither the infrequent deep encounters, nor the frequent and far ones; instead, the intermediate approaches have the most relevant contribution to the error growth. Such encounters are at a distance a few times the radius of the Hill's sphere of the planet (e.g. 3). An even simpler model allows us to get analytical estimates of the Lyapunov times in good agreement with the values coming from the model above and the numerical integrations. The predictability of the medium‐term evolution and the hazard posed to the Earth by those objects are analysed in the Discussion section. This revised version was published online in July 2006 with corrections to the Cover Date.  相似文献   

18.
It is already known (Froeschlé, Lega and Gonczi, 1997) that the Fast Lyapunov Indicator (FLI), that is the computation on a relatively short time of the largest Lyapunov indicator, allows to discriminate between ordered and weak chaotic motion. We have found that, under certain conditions, the FLI also discriminates between resonant and non-resonant orbits, not only for two-dimensional symplectic mappings but also for higher dimensional ones. Using this indicator, we present an example of the Arnold web detection for four and six-dimensional symplectic maps. We show that this method allows to detect the global transition of the system from an exponentially stable Nekhoroshevs like regime to the diffusive Chirikovs one.  相似文献   

19.
Two primary solar-activity indicators sunspot numbers(SNs)and sunspot areas(SAs)in the time interval from November 1874 to December 2012 are used to determine the chaotic and fractal properties of solar activity.The results show that(1)the long-term solar activity is governed by a low-dimensional chaotic strange attractor,and its fractal motion shows a long-term persistence on large scales;(2)both the fractal dimension and maximal Lyapunov exponent of SAs are larger than those of SNs,implying that the dynamical system of SAs is more chaotic and complex than SNs;(3)the predictions of solar activity should only be done for short-to mid-term behaviors due to its intrinsic complexity;moreover,the predictability time of SAs is obviously smaller than that of SNs and previous results.  相似文献   

20.
近地小行星(10302) 1989 ML和(4660) Nereus作为下一代深空探测的候选目标一直备受关注. 在考虑太阳系主要天体的动力学背景下, 通过计算最大Lyapunov指数(MLE)及MEGNO (Mean Exponential Growth factor of Nearby Orbits)指数讨论它们的稳定性. 同时, 对每个小行星, 在其观测误差范围内按多元正态分布各选取1000个克隆粒子, 通过统计分析显示这两个小行星在10万年内可能的运动范围, 给出半长径-偏心率空间中的出现次数分布图, 并统计小行星与地球或其他大行星之间的密近交汇及碰撞的概率. 此外还对这两个小行星的标称轨道进行长期共振、Kozai共振及平运动共振的动力学分析. 综上得出结论, 1989 ML处在平运动共振主导的区域, 发生密近交汇的概率较小, 从而其轨道相对较稳定; 而Nereus处在地球的密近交汇区域, 轨道极不稳定.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号