力梯度辛方法在圆型限制性三体问题中的应用

旋转坐标系下的圆型限制性三体问 题因含非惯性系所附加的影响部分使得动能不是动量的严格二次型, 可能导致力梯度辛积分算法的应用遇到困难. 从Lie算子运算出发, 严格论证了力梯度算子在这种情形下的物理意义 仍然像质心惯性坐标系下的圆型限制性三体问题那样是引力的梯度, 而不是引力与非惯性力所得合力的梯度, 表明了力梯度辛方法适合求解旋转坐标系下的圆型限制性三体问题. 通过应用四阶力梯度辛方法、最优化四阶力梯度辛方法和Forest-Ruth 辛方法分别求解该问题, 进行了数值对比研究, 结果显示最优化型力梯度算法能够取得最好精度. 还应用最优化型算法计算两邻近轨道的Lyapunov指数和快速Lyapunov指标, 确保高精度辛方法能够贯穿于这些混沌指标计算的全过程, 以便准确刻画此系统的动力学定性性质. 关键词： 辛积分器 圆型限制性三体问题 混沌 Lyapunov指数     

Discrete and continuous integrable systems possessing the same non-dynamical r-matrix

We consider two different Lax representations with the same Lax matrix in terms of 2 × 2 traceless matrices: one produces the discrete integrable symplectic mapping resulting from the well-known Toda spectral problem under the discrete Bargmann-Garnier (BG) constraint; the other generates the continuous non-linearized integrable system for the c-KdV spectra problem under the corresponding BG constraint. We are surprised to find that the two very different (one is discrete, the other continuous) integrable systems possess the same non-dynamical r-matrix.     

The Neumann Type Systems and Algebro-Geometric Solutions of a System of Coupled Integrable Equations

A system of (1+1)-dimensional coupled integrable equations is decomposed into a pair of new Neumann type systems that separate the spatial and temporal variables for this system over a symplectic submanifold. Then, the Neumann type flows associated with the coupled integrable equations are integrated on the complex tour of a Riemann surface. Finally, the algebro-geometric solutions expressed by Riemann theta functions of the system of coupled integrable equations are obtained by means of the Jacobi inversion.     

辛算法在限制性三体问题数值研究中的应用

限制性三体问题是太阳系动力学中常采用的一种力学模型,是一哈密顿(Hamilton)系统.由于数学工具的不够,一些重要问题只能进行数值研究,但要了解系统的演化状况,必须进行长期跟踪计算.因此,对算法要求极高,应能保持运动的整体特征,而Hamilton系统的辛算法正符合这一要求,文章将利用算法合成构造旋转坐标系中圆型和椭圆型限制性三体问题(对应不可分Hamilton系统)的显式辛差分格式,并以计算实例表明方法的有效性.     

A Bäcklund Transformation of the Restricted mKdV Flow with a Rosochatius Deformation

A Bäcklund transformation of the restricted mKdV flow with a Rosochatius deformation is constructed. Its Lax representation and thus N invariants in involution are presented. Such Bäcklund transformation is a Rosochatius deformation of that of the restricted mKdV flow.     

Completely Integrable Generalized C. Neumann Systems on Several Symplectic Submanifolds

Abstract

New completely integrable generalized C. Neumann systems on several symplectic submanifolds are presented, and the relations between the generalized C. Neumann systems and the spectral problems are further discussed in this paper. In particular, a new eigenvalue problem is proposed in Part 3.3.     

The modular automorphism group of a Poisson manifold

The modular flow of Poisson manifold is a 1-parameter group of automorphisms determined by the choice of a smooth density on the manifold. When the density is changed, the generator of the group changes by a hamiltonian vector field, so one has a 1-parameter group of “outer automorphisms” intrinsically attached to any Poisson manifold. The group is trivial if and only if the manifold admits a measure which is invariant under all hamiltonian flows.

The notion of modular flow in Poisson geometry is a classical limit of the notion of modular automorphism group in the theory of von Neumann algebras. In addition, the modular flow of a Poisson manifold is related to modular cohomology classes for associated Lie algebroids and symplectic groupoids. These objects have recently turned out to be important in Poincaré duality theory for Lie algebroids.     

Probability representation of the quantum evolution and energy-level equations for optical tomograms

The von Neumann evolution equation for the density matrix and the Moyal equation for the Wigner function are mapped onto the evolution equation for the optical tomogram of the quantum state. The connection with the known evolution equation for the symplectic tomogram of the quantum state is clarified. The stationary states corresponding to quantum energy levels are associated with the probability representation of the von Neumann and Moyal equations written for optical tomograms. The classical Liouville equation for optical tomogram is obtained. An example of the parametric oscillator is considered in detail.     

Neumann-Bogolyubov-Rosochatius Oscillatory Dynamical Systems and their Integrability Via Dual Moment Maps. Part I

Abstract

The finite-dimensional invariant subspaces of the solutions of intergrable by Lax infinite-dimensional Benney-Kaup dynamical system are presented. These invariant subspaces carry the canonical symplectic structure, with relation to which the Neumann type dynamical systems are Hamiltonian and Liouville intergrable ones. For the Neumann-Bogolyubov and Neumann-Rosochatius dynamical systems, the Lax-type representations via the dual moment maps into some deformed loop algebras as well as the finite hierarchies of conservation laws are constructed.     

Finite-Difference Lattice Boltzmann Scheme for High-Speed Compressible Flow： Two-Dimensional Case

Lattice Boltzmann （LB） modeling of high-speed compressible flows has long been attempted by various authors. One common weakness of most of previous models is the instability problem when the Mach number of the flow is large. In this paper we present a finite-difference LB model, which works for flows with flexible ratios of specific heats and a wide range of Mach number, from 0 to 30 or higher. Besides the discrete-velocity-model by Watari [Physica A 382 （2007） 502], a modified Lax Wendroff finite difference scheme and an artificial viscosity are introduced. The combination of the finite-difference scheme and the adding of artificial viscosity must find a balance of numerical stability versus accuracy. The proposed model is validated by recovering results of some well-known benchmark tests： shock tubes and shock reflections. The new model may be used to track shock waves and/or to study the non-equilibrium procedure in the transition between the regular and Mach reflections of shock waves, etc.     

A Finite-Dimensional Completely Integrable System Associated with Boussinesq Hierarchy

In this paper, a new completely integrable system related to the complex spectral problem -φ xx＋（i/4）wpx＋（i/4）（wp）x＋（1/4）vφ=iλφxand the constrained flows of the Boussinesq equations axe generated. According to the viewpoint of Hamiltonian mechanics, the Euler-Lagrange equations and the Legendre transformations, a reasonable Jacobi-Ostrogradsky coordinate system is obtained. Moreover, by means of the constrained conditions between the potentiaJ u, v and the eigenfunction φ, the involutive representations of the solutions for the Boussinesq equation hieraxchy axe given.     

Integrable Rosochatius Deformations of the Restricted cKdV Flows

Three novel finite-dimensional integrable Hamiltonian systems of Rosochatius type and their Lax representations are presented. We make a deformation for the Lax matrbces of the Neumann type, the Bargmann type and the high-order symmetry type of restricted cKdV flows by adding an additional term and then prove that this kind of deformation does not change the r-matrix relations. Finally the new integrable systems are generated from these deformed Lax matrices.     

Symplectic quantum mechanics

Using the notion of symplectic structure and Weyl (or star) product of non-commutative geometry, we construct unitary representations for the Galilei group and show how to rewrite the Schrödinger equation in phase space. This approach gives rise to a new procedure to derive Wigner functions without the use of the Liouville-von Neumann equation. Applications are presented by deriving the states of linear and nonlinear oscillators in terms of amplitudes of probability in phase space. The notion of coherent states is also discussed in this context.     

Integrable systems and invariant curve flows in centro-equiaffine symplectic geometry

In this paper, the theory for curves in centro-equiaffine symplectic geometry is established. Integrable systems satisfied by the curvatures of curves under inextensible motions in centro-equiaffine symplectic geometry are identified. It is shown that certain non-stretching invariant curve flows in centro-equiaffine symplectic geometry are closely related to the matrix KdV equations and their extension.     

A note on singular lagrangian submanifolds

Classification of singular lagrangian submanifolds which appear as images of a regular one under a symplectic relation, is considered from the point of view of standard singularity theory. The classification is carried out in small dimensions and restricted to special types of symplectic objects. Normal forms for singular pullbacks and pushforwards are given using an appropriate symplectic equivalence group. It is shown that the general classification problem reduces to the classification problem for appropriate mapping diagrams. An approach to the classical theories of phase transition is given based on the geometry of singular lagrangian images. The variational open swallowtails and regularly intersecting pairs of holonomic components are resolved using an appropriate reduction relation. Examples are given of singularities encountered in physics.     

Symplectic Schemes for Birkhoffian System

A universal symplectic structure for a Newtonian system including nonconservative cases can be constructed in the framework of Birkhoffian generalization of Hamiltonian mechanics. In this paper the symplectic geometry structure of Birkhoffian system is discussed, then the symplecticity of Birkhoffian phase flow is presented. Based on these properties we give a way to construct symplectic schemes for Birkhoffian systems by using the generating function method.     

A discrete action principle for electrodynamics and the construction of explicit symplectic integrators for linear,non-dispersive media

In this work, we derive a discrete action principle for electrodynamics that can be used to construct explicit symplectic integrators for Maxwell’s equations. Different integrators are constructed depending on the choice of discrete Lagrangian used to approximate the action. By combining discrete Lagrangians in an explicit symplectic partitioned Runge–Kutta method, an integrator capable of achieving any order of accuracy is obtained. Using the von Neumann stability analysis, we show that the integrators greatly increase the numerical stability and reduce the numerical dispersion compared to other methods. For practical purposes, we demonstrate how to implement the integrators using many features of the finite-difference time-domain method. However, our approach is also applicable to other spatial discretizations, such as those used in finite element methods. Using this implementation, numerical examples are presented that demonstrate the ability of the integrators to efficiently reduce and maintain a minimal amount of numerical dispersion, particularly when the time-step is less than the stability limit. The integrators are therefore advantageous for modeling large, inhomogeneous computational domains.     

Measures in the geometric quantization of field theories

In this article, it is shown that for the standard symplectic form on the space of compactly supported sections of a symplectic fibre bundle, there is no locally-finite Borel measure which is preserved by the Hamiltonian flows of even a quite restricted set of functions on this space. As this means that some of the operators arising in geometric quantization associated to classical observables would not be Hermitean, the result suggests that one should consider quotients by gauge groups as classical phase spaces to avoid this problem.     

Numerical simulation of three-dimensional nonlinear water waves

We present an accurate and efficient numerical model for the simulation of fully nonlinear (non-breaking), three-dimensional surface water waves on infinite or finite depth. As an extension of the work of Craig and Sulem [19], the numerical method is based on the reduction of the problem to a lower-dimensional Hamiltonian system involving surface quantities alone. This is accomplished by introducing the Dirichlet–Neumann operator which is described in terms of its Taylor series expansion in homogeneous powers of the surface elevation. Each term in this Taylor series can be computed efficiently using the fast Fourier transform. An important contribution of this paper is the development and implementation of a symplectic implicit scheme for the time integration of the Hamiltonian equations of motion, as well as detailed numerical tests on the convergence of the Dirichlet–Neumann operator. The performance of the model is illustrated by simulating the long-time evolution of two-dimensional steadily progressing waves, as well as the development of three-dimensional (short-crested) nonlinear waves, both in deep and shallow water.