Birkhoffian symplectic algorithms derived from Hamiltonian symplectic algorithms

In this paper, we focus on the construction of structure preserving algorithms for Birkhoffian systems, based on existing symplectic schemes for the Hamiltonian equations. The key of the method is to seek an invertible transformation which drives the Birkhoffian equations reduce to the Hamiltonian equations. When there exists such a transformation,applying the corresponding inverse map to symplectic discretization of the Hamiltonian equations, then resulting difference schemes are verified to be Birkhoffian symplectic for the original Birkhoffian equations. To illustrate the operation process of the method, we construct several desirable algorithms for the linear damped oscillator and the single pendulum with linear dissipation respectively. All of them exhibit excellent numerical behavior, especially in preserving conserved quantities.

A Hierarchy of Integrable Nonlinear Lattice Equations and New Integrable Symplectic Map

A discrete spectral problem is discussed, and a hierarchy of integrable nonlinear lattice equations related tothis spectral problem is devised. The new integrable symplectic map and finite-dimensional integrable systems are givenby nonlinearization method. The binary Bargmann constraint gives rise to a Backlund transformation for the resultingintegrable lattice equations.

A Hierarchy of Integrable Lattice Soliton Equations and New Integrable Symplectic Map

开始从一分离光谱问题， integrable latticesoliton 方程的一个层次被导出。层次是，这被显示出完全在 theLiouville 感觉的 integrable 并且拥有分离 bi-Hamiltonian 结构。一张新 integrable symplectic 地图 andfinite 维的 integrable 系统被 nonlinearization 方法给。二进制 Bargmannconstraint 为产生 integrable 格子方程产生 Baecklund 转变。最后，层次的能量守恒定律被介绍。

A Hierarchy of Integrable Nonlinear Lattice Equations and New Integrable Symplectic Map

A discrete spectral problem is discussed, and a hierarchy of integrable nonlinear lattice equations related tothis spectral problem is devised. The new integrable symplectic map and finite-dimensional integrable systems are givenby nonlinearization method. The binary Bargmann constraint gives rise to a Backlund transformation for the resultingintegrable lattice equations.

Hamiltonian Representation of Higher Order Partial Differential Equations with Boundary Energy Flows

This paper presents a system representation that can be applied to the description of the interaction between systems connected through common boundaries. The systems consist of partial differential equations that are first order with respect to time, but spatially higher order. The representation is derived from the instantaneous multisymplectic Hamiltonian formalism;therefore, it possesses the physical consistency with respect to energy. In the interconnection, particular pairs of control inputs and observing outputs, called port variables, defined on the boundaries are used. The port variables are systematically introduced from the representation.

Structure-preserving algorithms for autonomous Birkhoffian systems

The Pfaff-Birkhoff variational principle is discretized,and based on the discrete variational principle the discrete Birkhoffian equations are obtained.Taking the discrete equations as an algorithm,the corresponding discrete flow is proved to be symplectic.That means the algorithm preserves the symplectic structure of Birkhoffian systems.Finally,simulation results of the given example indicate that structure-preserving algorithms have great advantage in stability and energy conserving.

Recent Progress in Symplectic Algorithms for Use in Quantum Systems

In this paper we survey recent progress in symplectic algorithms for use in quantum systems in the following topics:Symplectic schemes for solving Hamiltonian systems;Classical trajectories of diatomic systems,model molecule A2B,Hydrogen ion H+2 and elementary atmospheric reaction N(4S)+O2(X 3Σ−g)→NO(X 2Π)+O(3P)calculated by means of Runge-Kutta methods and symplectic methods;the classical dissociation of the HF molecule and classical dynamics of H+2 in an intense laser field;the symplectic form and symplectic-scheme shooting method for the time-independent Schr¨odinger equation;the computation of continuum eigenfunction of the Schr¨odinger equation;asymptotic boundary conditions for solving the time-dependent Schr¨odinger equation of an atom in an intense laser field;symplectic discretization based on asymptotic boundary condition and the numerical eigenfunction expansion;and applications in computing multi-photon ionization,above-threshold ionization,Rabbi oscillation and high-order harmonic generation of laser-atom interaction.

Solution of Algebraic Lyapunov Equation on Positive-Definite Hermitian Matrices by Using Extended Hamiltonian Algorithm

This communique is opted to study the approximate solution of the Algebraic Lyapunov equation on the manifold of positive-definite Hermitian matrices.We choose the geodesic distance between􀀀AHX􀀀XA and P as the cost function,and put forward the Extended Hamiltonian algorithm(EHA)and Natural gradient algorithm(NGA)for the solution.Finally,several numerical experiments give you an idea about the effectiveness of the proposed algorithms.We also show the comparison between these two algorithms EHA and NGA.Obtained results are provided and analyzed graphically.We also conclude that the extended Hamiltonian algorithm has better convergence speed than the natural gradient algorithm,whereas the trajectory of the solution matrix is optimal in case of Natural gradient algorithm(NGA)as compared to Extended Hamiltonian Algorithm(EHA).The aim of this paper is to show that the Extended Hamiltonian algorithm(EHA)has superior convergence properties as compared to Natural gradient algorithm(NGA).Upto the best of author’s knowledge,no approximate solution of the Algebraic Lyapunov equation on the manifold of positive-definite Hermitian matrices is found so far in the literature.

A quasi-dynamic model and a symplectic algorithm of super slender Kirchhoff rod

Recently the Kirchhoff rod and the methods of dynamical analogue have been widely used in modeling DNA.The features of a DNA such as its super slender and super large deformation raise new challenges in modeling and numerical simulations of a Kirchhoff rod.In this paper,Euler parameters are introduced to set up the quasi-Hamilton system of an elastic rod,then a symplectic algorithm is applied in its numerical simulations.Finally,a simplified surface model of the rod is given based on the hypothesis of rigid cross-section.

Least-squares reverse time migration in visco-acoustic media based on symplectic stereo-modeling method

The authors proposed a symplectic stereo-modeling method(SSM)in the Birkhoffian dynam-ics and apply it to the visco-acoustic least-squares reverse time migration(LSRTM).The SSM adopts ste-reo-modeling operator in space and symplectic Runge-Kutta scheme in time,resulting in great ability in suppressing numerical dispersion and long-time computing.These advantages are further confirmed by numerical dispersion analysis,long-time computation test and computational efficiency comparison.After these theoretical analyses and experiments,acoustic and visco-acoustic LSRTM are tested and compared between SSM method and the conventional symplectic method(CSM)using the fault and marmousi models.Meanwhile,dynamic source encoding and exponential decay moving average gradients method are adopted to reduce the computation cost and improve the convergence rate.The imaging results show that LSRTM based on visco-acoustic wave equations effectively takes into account the influence of viscosity can therefore compensate for the amplitude attenuation.Besides,SSM method not only has high numerical accuracy and computational efficiency,but also performs effectively in LSRTM.

一种Birkhoff形式下结构动响应问题的保辛中点格式

结构动响应预测是结构设计的基础,是结构振动控制、载荷识别的前提。本文在辛体系下针对结构动响应问题,提出了一种Birkhoff形式下的保辛中点格式。首先引入状态变量,并基于摄动方法将结构动响应方程转化为线性自治Birkhoff方程的形式,进一步利用中心差分推导出线性自治Birkhoff方程的中点格式,其证明是保辛的。该格式不要求Birkhoff方程系数矩阵非奇异,因此适用于奇数维系统。两个不同数值算例的结果充分验证了本文方法的卓越性,也凸显了相对于传统算法在计算精确度和稳定性方面的明显优势。

高阶Maggi方程的Birkhoff化及其辛算法

针对非完整系统的高阶Maggi方程,在满足一定的条件时,可以对其进行Birkhoff化.通过构造生成函数,利用Birkhoff广义辛算法对其进行数值仿真.仿真结果和传统的Runge-Kutta算法结果相比较,Birkhoff广义辛算法在长期跟踪后更加准确.

Mindlin层合圆柱壳理论的Hamiltonian结构与辛解析解

通过引进状态变量及其对偶变量。建立Ｍｉｎｄｌｉｎ层合圆柱壳的Ｈａｍｉｌｔｏｎ正则方程．在辛几何数学框架下，采用共轭辛正交归一关系给出各种复杂边界条件下的精确解．

The symplectic eigenfunction expansion theorem and its application to the plate bending equation

This paper deals with the bending problem of rectangular plates with two opposite edges simply supported.It is proved that there exists no normed symplectic orthogonal eigenfunction system for the associated infinite-dimensional Hamiltonian operator H and that the two block operators belonging to Hamiltonian operator H possess two normed symplectic orthogonal eigenfunction systems in some space.It is demonstrated by using the properties of the block operators that the above bending problem can be solved by the symplectic eigenfunction expansion theorem,thereby obtaining analytical solutions of rectangular plates with two opposite edges simply supported and the other two edges supported in any manner.

层合结构Hamiltonian元弱形式的辛方法

文章基于力学平衡方程,在柱坐标系下,导出正交各向异性层合柱壳混合方程和边界条件算子的弱形式,继而给出层合结构的Hamilton正则方程,建立半离散半解析的Hamiltonian元的微分方程,用弱形式给出的微分方程和边界条件对函数的连续性要求降低了,用于解决实际的工程问题常常比原始的微分方程更逼近真正解;针对其Hamiltonian元的矩阵结构,构建分析计算Hamiltonian元弱形式的辛方法,Hamilton结构在辛约化过程中得到充分保证,文中提出的辛方法简易可行,具有较强的有效性和稳定性。

Variational Principles and Hamiltonian Formulation for Nonlinear Water Waves

Ｖａｒｉａｔｉｏｎａｌｐｒｉｎｃｉｐｌｅｓｍａｙｓｕｃｃｉｎｃｔｌｙｌｅａｄｔｏｅｑｕａｔｉｏｎｓｏｆｍｏｔｉｏｎｆｏｒｗａｔｅｒｗａｖｅｓ，ａｌｏｗｉｎｓｉｇｈｔｉｎｔｏｔｈｅｅｆｅｃｔｏｆｐａｒａｍｅｔｅｒｓ，ａｎｄｐｒｏｖｉｄｅａｐａｔｈｆｏｒ...

Explicit Multi-Symplectic Methods for Hamiltonian Wave Equations

In this paper,based on the multi-symplecticity of concatenating symplectic Runge-Kutta-Nystrom(SRKN)methods and symplectic Runge-Kutta-type methods for numerically solving Hamiltonian PDEs,explicit multi-symplectic schemes are constructed and investigated,where the nonlinear wave equation is taken as a model problem.Numerical comparisons are made to illustrate the effectiveness of our newly derived explicit multi-symplectic integrators.

Structure-preserving algorithms for the Duffng equation

In this paper, the dissipative and the forced terms of the Duffng equation are considered as the perturbations of nonlinear Hamiltonian equations and the perturbational effect is indicated by parameter ε. Firstly, based on the gradientHamiltonian decomposition theory of vector fields, by using splitting methods, this paper constructs structure-preserving algorithms (SPAs) for the Duffng equation. Then, according to the Liouville formula, it proves that the Jacobian matrix determinants of the SPAs are equal to that of the exact flow of the Duffng equation. However, considering the explicit Runge-Kutta methods, this paper finds that there is an error term of order p+1 for the Jacobian matrix determinants. The volume evolution law of a given region in phase space is discussed for different algorithms, respectively. As a result, the sum of Lyapunov exponents is exactly invariable for the SPAs proposed in this paper. Finally, through numerical experiments, relative norm errors and absolute energy errors of phase trajectories of the SPAs and the Heun method (a second-order Runge-Kutta method) are compared. Computational results illustrate that the SPAs are evidently better than the Heun method when ε is small or equal to zero.

大气动力学方程的Hamilton算法

将大气动力学方程组写成正则算子方程的形式 ,通过引入泊松括号 ,并利用原方程组的无穷个不变量 ,深入研究其Hamilton性质 ,在忽略摩擦和外源强迫的情况下 ,证明了大气动力体系是一个Hamilton系统 ,从而构造出求解它的辛算法 ,并用数值试验检验了该算法 .

非线性Schrdinger方程的动力学行为分析

采用辛算法数值求解非线性Schr dinger方程的周期初值问题,建立不同的相空间来分析其动力学特性.首先比较分析了不同的相空间中立方非线性Schr dinger方程在不同立方非线性参数下的长时间演化的动力学特性,然后讨论了相空间中立方-五次方非线性Schr dinger方程在不同立方和五次方非线性参数下的长时间演化的动力学特性,数值结果显示,对于不同的立方非线性参数,随着五次方非线性参数的增加,动力学行为的演化路径是不一样的.