首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 437 毫秒
1.
辛方法的校正公式   总被引:5,自引:5,他引:0  
伍歆  黄天衣  万晓生 《天文学报》2002,43(4):391-402
1996年Wisdom等提出了对辛方法进行校正的概念和实践,现在继续对辛校正进行详尽讨论和数值比较,尤其对哈密顿函数可分解为一个主要部分和多个次要部分的一般情形,用Lie级数推导任意阶的各种辛算法的一次和二次辛校正公式并对一些算法给出具体的辛校正公式。又以日、木、土三体问题为模型进行数值实验,结果表明一次辛校正能提高精度,改善数值稳定性。计算效率也比较高,因而值得推荐使用,辛方法通常用大步长数值积分,这时二次辛校正并没有显著提高结果的精度,却大大增加了计算时间,不应予以推荐。  相似文献   

2.
In the present paper we give a proof of Nekhoroshev's theorem, which is concerned with an exponential estimate for the stability times in nearly integrable Hamiltonian systems. At variance with the already published proof, which refers to the case of an unperturbed Hamiltonian having the generic property of steepness, we consider here the particular case of a convex unperturbed Hamiltonian. The corresponding simplification in the proof might be convenient for an introduction to the subject.  相似文献   

3.
We study strongly and weakly integrable 2-dimensional Hamiltonian systems with velocity dependent potentials. We determine the set of conditions which must be satisfied in order to allow the existence of an independent second invariant polynomial in the momenta. We then investigate the linear case for which a complete solution of the problem can be obtained. We recover the classical set of linear strongly integrable systems and provide several new examples of weakly integrable systems whose equations of motion can be explicitly solved at a fixed value of the energy.  相似文献   

4.
The paper deals with the problem of the existence of a normal form for a nearly-integrable real-analytic Hamiltonian with aperiodically time-dependent perturbation decaying (slowly) in time. In particular, in the case of an isochronous integrable part, the system can be cast in an exact normal form, regardless of the properties of the frequency vector. The general case is treated by a suitable adaptation of the finite order normalization techniques usually used for Nekhoroshev arguments. The key point is that the so called “geometric part” is not necessary in this case. As a consequence, no hypotheses on the integrable part are required, apart from analyticity. The work, based on two different perturbative approaches developed by Giorgilli et al., is a generalisation of the techniques used by the same authors to treat more specific aperiodically time-dependent problems.  相似文献   

5.
Differentiable Hamiltonian systems close to nondegenerate, integrable Hamiltonian systems are shown to be integrable on a Cantor set in the sense that on some Cantor set, (i) the invariant KAM-tori form a smooth foliation, (ii) there exist smooth, independent integrals in involution, and (iii) there exists a complete solution of the Hamilton Jacobi equation. The complement of the Cantor set is shown to be small in measure.Paper presented at the 1981 Oberwolfach Conference on Mathematical Methods in Celestial Mechanics.  相似文献   

6.
Most existing satellite relative motion theories utilize mean elements, and therefore cannot be used for calculating long-term bounded perturbed relative orbits. The goal of the current paper is to find an integrable approximation for the relative motion problem under the J 2 perturbation, which is adequate for long-term prediction of bounded relative orbits with arbitrary inclinations. To that end, a radial intermediary Hamiltonian is utilized. The intermediary Hamiltonian retains the original structure of the full J 2 Hamiltonian, excluding the latitude dependence. This formalism provides integrability via separation, a fact that is utilized for finding periodic relative orbits in a local-vertical local-horizontal frame and determine an initialization scheme that yields long-term boundedness of the relative distance. Numerical experiments show that the intermediary-based computation of orbits provides long-term bounded orbits in the full J 2 problem for various inclinations. In addition, a test case is shown in which the radial intermediary-based initial conditions of the chief and deputy satellites yield bounded relative distance in a high-precision orbit propagator.  相似文献   

7.
Deprit and Miller have conjectured that normalization of integrable Hamiltonians may produce normal forms exhibiting degenerate equilibria to very high order. Several examples in the class of coupled elliptic oscillators are known. In order to test the utility of normalization as a detector of integrability we normalize, to high order, a perturbed Keplerian system known to have several integrable limits; the generalized van der Waals Hamiltonian for a hydrogen atom. While the separable limits give rise to high order degeneracy we find a non-separable, integrable limit for which the normal form does not exhibit degeneracy. We conclude that normalization may, in certain cases, indicate integrability but is not guaranteed to uncover all integrable limits.  相似文献   

8.
Numerical evidence is presented which indicates that, although the third integral is tangent to the Hamiltonian (energy integral) along some periodic orbits (as has been shown by Goudas), it is not tangent to it along non-periodic orbits; therefore it is not a function of the Hamiltonian. The set of periodic orbits is probably dense in general, but a given form of the third integral is valid in the neighbourhood of a limited number of them; no form of the third integral is valid for all periodic orbits, except in integrable cases.  相似文献   

9.
We give a constructive proof of the existence of elliptic lower dimensional tori in nearly integrable Hamiltonian systems. In particular we adapt the classical Kolmogorov normalization algorithm to the case of planetary systems, for which elliptic tori may be used as replacements of elliptic Keplerian orbits in Lagrange-Laplace theory. With this paper we support with rigorous convergence estimates the semi-analytic work in our previous article (Sansottera et al., Celest Mech Dyn Astron 111:337–361, 2011), where an explicit calculation of an invariant torus for a planar model of the Sun-Jupiter-Saturn-Uranus system has been made. With respect to previous works on the same subject we exploit the characteristic of Lie series giving a precise control of all terms generated by our algorithm. This allows us to slightly relax the non-resonance conditions on the frequencies.  相似文献   

10.
We review theorems for proving non-integrability of Hamiltonian dynamical systems, which are based on properties of the variational equations in real or complex time or on the destruction of the resonant tori of an integrable system under a perturbation.  相似文献   

11.
An explicit symplectic integrator is constructed for the problem of a rotating planetary satellite on a Keplerian orbit. The spin vector is fixed perpendicularly to the orbital plane. The integrator is constructed according to the Wisdom-Holman approach: the Hamiltonian is separated in two parts so that one of them is multiplied by a small parameter. The parameter depends on the satellite’s shape or the eccentricity of its orbit. The leading part of the Hamiltonian for small eccentricity orbits is similar to the simple pendulum and hence integrable; the perturbation does not depend on angular momentum which implies a trivial ‘kick’ solution. In spite of the necessity to evaluate elliptic function at each step, the explicit symplectic integrator proves to be quite efficient. This revised version was published online in July 2006 with corrections to the Cover Date.  相似文献   

12.
In this paper, we study the existence of transversal homoclinic orbits in a planar circular restricted four-body problem, based on the perturbation theory of integrable Hamiltonian systems. We start from a planar circular restricted four-body model and regard it as a perturbation of the two-body model. Then, in order to conveniently study unbounded orbits, we transform the infinite points to finite points by a non-canonical transformation, arriving at a non-Hamiltonian system with degenerate fixed points. According to the extended Melnikov method, we finally prove that there exist transversal homoclinic orbits in this four-body model.  相似文献   

13.
The D'Alembert model for the spin/orbit problem in celestial mechanics is considered. Using a Hamiltonian formalism, it is shown that in a small neighborhood of a p:q spin/orbit resonance with (p,q) different from (1,1) and (2,1) the 'effective' D'Alembert Hamiltonian is a completely integrable system with phase space foliated by maximal invariant curves; instead, in a small neighborhood of a p:q spin/orbit resonance with (p,q) equal to (1,1) or (2,1) the 'effective' D'Alembert Hamiltonian has a phase portrait similar to that of the standard pendulum (elliptic and hyperbolic equilibria, separatrices, invariant curves of different homotopy). A fast averaging with respect to the 'mean anomaly' is also performed (by means of Nekhoroshev techniques) showing that, up to exponentially small terms, the resonant D'Alembert Hamiltonian is described by a two-degrees-of-freedom, properly degenerate Hamiltonian having the lowest order terms corresponding to the 'effective' Hamiltonian mentioned above.  相似文献   

14.
For a Hamiltonian that can be separated into N+1(N\geq 2) integrable parts, four algorithms can be built for a symplectic integrator. This research compares these algorithms for the first and second order integrators. We found that they have similar local truncation errors represented by error Hamiltonian but rather different numerical stability. When the computation of the main part of the Hamiltonian, H 0, is not expensive, we recommend to use S * type algorithm, which cuts the calculation of the H 0 system into several small time steps as Malhotra(1991) did. As to the order of the N+1 parts in one step calculation, we found that from the large to small would get a slower error accumulation. This revised version was published online in July 2006 with corrections to the Cover Date.  相似文献   

15.
The geometric approach to mechanics based on the Jacobi metric allows to easily construct natural mechanical systems which are integrable (actually separable) at a fixed value of the energy. The aim of the present paper is to investigate the dynamics of a simple prototype system outside the zero-energy hypersurface. We find that the general situation is that in which integrability is not preserved at arbitrary values of the energy. The structure of the Hamiltonian in the separating coordinates at zero energy allows a perturbation treatment of this system at energies slightly different from zero, by which we obtain an analytical proof of non-integrability.This revised version was published online in October 2005 with corrections to the Cover Date.  相似文献   

16.
In this paper, following the idea of constructing the mixed symplectic integrator (MSI) for a separable Hamiltonian system, we give a low order mixed symplectic integrator for an inseparable, but nearly integrable, Hamiltonian system, Although the difference schemes of the integrators are implicit, they not only have a small truncation error but, due to near integrability, also a faster convergence rate of iterative solution than ordinary implicit integrators, Moreover, these second order integrators are time-reversible.  相似文献   

17.
An extension ton degrees of freedom of the fact is established that forn=1 the time and the energy constant are canonically conjugate variables. This extension is useful in some cases to get action-angle variables from the general solution of a given integrable Hamiltonian system. As an example the Delaunay variables are proved to be canonical.Proceedings of the Sixth Conference on Mathematical Methods in Celestial Mechanics held at Oberwolfach (West Germany) from 14 to 19 August, 1978.  相似文献   

18.
We present a technique to control chaos in Hamiltonian systems which are close to integrable. By adding a small and simple control term to the perturbation, the system becomes more regular than the original one. We apply this technique to a forced pendulum model and show numerically that the control is able to drastically reduce chaos.  相似文献   

19.
This paper reviews various mapping techniques used in dynamical astronomy. It is mostly dealing with symplectic mappings. It is shown that used mappings can be usually interpreted as symplectic integrators. It is not necessary to introduce any functions it is just sufficient to split Hamiltonian into integrable parts. Actually it may be shown that exact mapping with function in the Hamiltonian may be non-symplectic. The application to the study of asteroid belt is emphasised but the possible use of mapping in planetary evolution studies, cometary and other problems is shortly discussed.  相似文献   

20.
In a recent paper Ballersteros and Ragnisco (1998) have proposed a new method of constructing integrable Hamiltonian systems. A new class of integrable systems may be devised using the following sequence: , where A is a Lie algebra is a Lie–Poisson structure on R 3, C is a Casimir for is a reduced Poisson bracket and (A, ▵) is a bialgebra. We study the relation between a Lie-Poisson stucture Λ and a reduced Poisson bracket , which is a key element in using the Lie algebra A to constructing this sequence. New examples of Lie algebras and their related integrable Hamiltonian systems are given. This revised version was published online in July 2006 with corrections to the Cover Date.  相似文献   

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

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