Symplectic Structure of Poisson System
-
摘要: 当Poisson系统中的Poisson矩阵是非常数时,经典的辛方法如辛Runge-Kutta方法,生成函数法一般不能保持Poisson系统的Poisson结构,利用非线性变换可把非常数Poisson结构转化成辛结构,然后任意阶的辛方法可以长时间计算Poisson系统的辛结构.自由刚体问题中Euler方程被转换成辛结构并用辛中点格式进行数值求解,数值结果给出了这种非线性变换的有效性.Abstract: When the Poisson matrix of Poisson system is non-constant,classical symplectic methods,such as symplectic Runge-Kutta method,generating function method,cannot preserve the Poisson structure.The non-constant Poisson structure was transformed into the symplectic structure by the nonlinear transform.Arbitrary order symplectic method was applied to the transformed Poisson system.The Euler equation of the free rigid body problem was transformed into the symplectic structure and computed by the mid-point scheme.Numerical results show the effectiveness of the nonlinear transform.
-
Key words:
- Poisson system /
- nonlinear transformation /
- symplectic method /
- rigid body problem
-
[1] HONG Jia-lin.A novel numerical approach to simulate nonlinear Schrdinger equations with varing coefficients[J].Applied Mathematical Letter,2003,16:759—765. doi: 10.1016/S0893-9659(03)00079-X [2] Feng K.Difference schemes for Hamiltonian formulation and symplectic geometry[J].J Comp Math,1986,4(3):279—289. [3] Feng K,Wu H M, Qin M Z,et al.Construction of canonical difference schemes for Hamiltonian formalism via generating functions[J].J Comp Math,1989,7(1):71—96. [4] QIN Meng-zhao,Li S T.A note for Lie-Poisson Hamiltonian-Jacobi equation and Lie-Poisson integrator[J].Computers Mathematical Application,1995,30(7):67—74. doi: 10.1016/0898-1221(95)00126-J [5] McLachlan R I.Explicit Lie-Poisson integration and the Euler equations[J].Phys Rev Lett,1993,71:3043—3046. doi: 10.1103/PhysRevLett.71.3043 [6] Li S T,QIN Meng-zhao.Lie-Poisson integration for rigid body dynamics[J].Computers Mathmatical Application,1995,30(9):105—118. [7] Zhu W,Qin M.Poisson schemes for Hamiltonian system on Poisson manifolds[J].Computers Mathematical Application,1994,27(12):7—16.
计量
- 文章访问数: 3087
- HTML全文浏览量: 102
- PDF下载量: 503
- 被引次数: 0