Symplectic Mappings of Co-orbital Motion in the Restricted Problem of Three Bodies

The dynamics of co-orbital motion in the restricted three-body problem are investigated by symplectic mappings. Analytical and semi-numerical mappings have been developed and studied in detail. The mappings have been tested by numerical integration of the equations of motion. These mappings have been proved to be useful for a quick determination of the phase space structure reflecting the main characteristics of the dynamics of the co-orbital problem.     

Structures in the phase space of a four dimensional symplectic map

We numerically investigate the projections of non periodic orbits in a 4-dimensional (4-D) symplectic map composed of two coupled 2-dimensional (2-D) maps. We describe in detail the structures that are produced in different planes of projection and we find how the morphology of the 4-D orbits is influenced by the features of the 2-D maps as the coupling parameter increases. We give an empirical law that describes this influence.     

Frequency analysis of a dynamical system

Frequency analysis is a new method for analyzing the stability of orbits in a conservative dynamical system. It was first devised in order to study the stability of the solar system (Laskar, Icarus, 88, 1990). It is a powerful method for analyzing weakly chaotic motion in hamiltonian systems or symplectic maps. For regular motions, it yields an analytical representation of the solutions. In cases of 2 degrees of freedom system with monotonous torsion, precise numerical criterions for the destruction of KAM tori can be found. For a 4D symplectic map, plotting the frequency map in the frequency plane provides a clear representation of the global dynamics and describes the actual Arnold web of the system.     

Mass-Weighted Symplectic Forms for the N-Body Problem

Mass-weighted symplectic forms provide a unified framework for the treatment of both finite and vanishingly small masses in the N-body problem. These forms are introduced, compared to previous approaches, and their properties are discussed. Applications to symplectic mappings, the definition of action-angle variables for the Kepler problem, and Hamiltonian perturbation theory are outlined This revised version was published online in July 2006 with corrections to the Cover Date.     

An Explicit Symplectic Integrator for Cometary Orbits

A new symplectic algorithm is developed for cometary orbit integrations. The integrator can handle both high-eccentricity orbits and close encounters with planets. The method is based on time transformations for Hamiltonians separated into Keplerian and perturbation parts. The adaptive time-step of this algorithm depends on the distance from a centre and the magnitude of perturbations. The explicit leapfrog technique is simple and efficient.     

基于辛算法模拟探地雷达在复杂地电模型中的传播

近年来,探地雷达(GPR)凭借其快速、高效、无破损等特点,已经广泛应用于浅地层目标探测中.数值模拟是研究探地雷达电磁波在地下结构中传播规律的有效手段.辛算法是一种保持Hamilton系统总能量不变的时域数值计算方法.本文提出了基于一阶显式辛分块龙格库塔方法的探地雷达数值模拟方法.通过对比本文算法与时域有限差分方法计算结果可知,在同等计算精度下,本文算法可以节省25%的计算时间.并基于本文算法对两个复杂GPR模型进行正演模拟,得到模拟GPR探测wiggle图,这有助于更好的理解和分析实测雷达数据.     

地球流体力学的研究与进展

简要介绍中国科学院大气物理研究所七十多年来在理论与计算地球流体力学方面的若干研究及其新的进展.在理论地球流体力学方面,介绍了长波动力学及线性稳定性问题、弱非线性理论及行星波动力学以及用Arnold方法(能量-Casimir方法)研究大气和海洋中各种流体运动的非线性稳定性问题的成果.此外,对扰动演变、扰动和基流相互作用及热带大气动力学中的第二类不稳定条件(CISK)也作了简要的介绍.在计算地球流体力学方面,主要内容包括:用物理观点和数学分析相结合的方法阐述了造成计算紊乱和计算不稳定的机理,论证计算稳定性、算     

Recent progress in the theory and application of symplectic integrators

In this paper various aspect of symplectic integrators are reviewed. Symplectic integrators are numerical integration methods for Hamiltonian systems which are designed to conserve the symplectic structure exactly as the original flow. There are explicit symplectic schemes for systems of the formH=T(p)+V(q), and implicit schemes for general Hamiltonian systems. As a general property, symplectic integrators conserve the energy quite well and therefore an artificial damping (excitation) caused by the accumulation of the local truncation error cannot occur. Symplectic integrators have been applied to the Kepler problem, the motion of minor bodies in the solar system and the long-term evolution of outer planets.     

Existence of formal integrals of symplectic integrators

A recurrent method of solving the formal integrals of symplectic integrators is given. The special examples show that there are no long-term variations in all integrals of the Hamiltonian system in addition to the energy one when symplectic integrators are used in the numerical studies of the system. As an application of the formal integrals, the relation between them and the linear stability of symplectic integrators is discussed.     

A symplectic mapping approach to the study of the stochasticity in asteroidal resonances

In the case of the 2:1 and 3:2 resonances with Jupiter, it has not been yet possible to have a complete identification of all chaotic diffusion processes at work, mainly because the time scale of some of them are of an order still out of the reach of precise integrations. A planar Hadjidemetriou's mapping, using expansions valid for high eccentricities and scaled in order to accelerate the diffusion processes, was derived. The solutions obtained with the mapping show huge eccentricity variations in all orbits starting in the middle of the 2:I resonance, when the main short-period perturbations of Jupiter's orbit are considered. The solutions starting in the middle of the 3:2 resonance do not show any important diffusion.