带三次非线性项的四阶Schr(o)dinger方程的分裂多辛算法

对一类带三次非线性项的四阶Schr(o)dinger方程提出分裂多辛格式.其基本思想是将多辛算法和分裂方法相结合,既具有多辛格式固有的保多辛几何结构的特性,又发挥了分裂方法在计算上灵活高效的特点.数值实验结果表明,分裂多辛格式比其它传统的多辛格式更节约计算时间和计算机的内存,从而更加优越.     

广义非线性Schrdinger方程的多辛格式与模方守恒律

通过正则变换,构造出广义非线性Schr(o)dinger方程的多辛方程组.对此多辛方程组,导出了一个新的模方守恒多辛格式.数值实验结果表明,多辛格式具有长时间的数值行为,且在保持模方守恒律方面优于蛙跳格式和辛欧拉中点格式.     

基于多步高阶差分格式求解时域Maxwell方程

提出基于无穷维哈密尔顿系统及分裂算子理论的多步高阶差分格式,求解时域Maxwell方程.在时间方向上,针对Maxwell方程采用不同阶数的辛算法进行差分离散;在空间方向上,采用四阶差分格式进行差分离散.探讨多步高阶差分格式的稳定性及数值色散性,最后给出数值计算结果.结果表明,五级四阶格式为最有效的多步高阶差分格式,具有高精度、占用较少的计算机资源等优点,适用于长时间的数值模拟.     

一类特殊非完整力学系统的辛算法计算

对一类特殊的非完整力学系统的动力学性质进行数值研究,采用当前比较优越的保结构算法进行数值计算,并和传统的Runge-Kutta方法进行比较. 通过计算结果的比较而得出辛算法在这类特殊的非完整力学系统的数值计算中的优越性. 关键词： 非完整约束 辛算法 辛差分格式     

带三次非线性项的四阶Schrdinger方程的分裂多辛算法(英文)

对一类带三次非线性项的四阶Schrdinger方程提出分裂多辛格式。其基本思想是将多辛算法和分裂方法相结合,既具有多辛格式固有的保多辛几何结构的特性,又发挥了分裂方法在计算上灵活高效的特点。数值实验结果表明,分裂多辛格式比其它传统的多辛格式更节约计算时间和计算机的内存,从而更加优越.     

SRLW方程的多辛Fourier谱格式及其守恒律

通过引进正则动量,将对称正则长波方程（简称SRLW方程）转化成多辛形式的方程组,它具有多辛守恒律;介绍了空间方向满足周期边界条件的函数的Fourier谱方法;对SRLW方程的多辛方程组在空间方向利用Fourer谱方法,时间方向上应用Euler中点格式离散,得到其多辛Fourier拟谱格式;证明此格式的一些离散守恒律．用此格式模拟了SRLW方程的单个孤立波,还模拟了多个孤立波的追赶、碰撞和分离过程．     

耦合非线性Schrdinger系统的多辛差分格式

近年来,Bridges等人在Hamiltonian力学意义下,直接把有限维Hamiltonian系统推广到无穷维,通过引入新的函数坐标,使得偏微分方程在时间和空间的各个方向上都有各自不同的有限维辛结构,这样原偏微分方程就由各个有限维辛结构以及右端的梯度函数决定,称这样的方程为多辛Hamiltonian系统.多辛Hamiltonian系统满足多辛守恒定律,满足多辛Hamiltonian系统的多辛守恒律的离散算法称为多辛算法.以耦合非线性Schr dinger方程为例,研究无穷维Hamiltonian系统的多辛算法,验证了两孤立子碰撞后会发生相互通过、反射及融合现象.     

广义Zakharov-Kuznetsov方程的多辛Preissmann格式

广义Zakharov-Kuznetsov 方程作为一类重要的非线性方程有着许广泛的应 用前景,基于Hamilton 空间体系的多辛理论研究了广义Zakharov-Kuznetsov方程的数值 解法,讨论了利用Preissmann 方法构造离散多辛格式的途径, 并构造了一种典型的半隐 式的多辛格式, 该格式满足多辛守恒律、局部能量守恒律. 数值算例结果表明该多辛离 散格式具有较好的长时间数值稳定性.     

广义Zakharov-Kuznetsov方程的多辛Preissmann格式

广义Zakharov-Kuznetsov方程作为一类重要的非线性方程有着许多广泛的应用前景,基于Hamilton空间体系的多辛理论研究了广义Zakharov-Kuznetsov方程的数值解法,讨论了利用Preissmann方法构造离散多辛格式的途径,并构造了一种典型的半隐式的多辛格式,该格式满足多辛守恒律、局部能量守恒律.数值算例结果表明该多辛离散格式具有较好的长时间数值稳定性.     

时间相关外场中量子系统时间演化的辛格式

对显含时间的可分线性哈密顿系统构造了2阶模方守恒-辛格式和显式辛格式,并计算了1维有限宽无限深势阱中的电子与模拟激光场的相互作用,结果与理论分析一致。为数值研究时间相关外场中的量子系统,特别是强激光与原子相互作用提供了合理和有效的方法。     

Explicit Multi-Symplectic Splitting Methods for the Nonlinear Dirac Equation

In this paper, we propose two new explicit multi-symplectic splitting methods for the nonlinear Dirac (NLD) equation. Based on its multi-symplectic formulation, the NLD equation is split into one linear multi-symplectic system and one nonlinear infinite Hamiltonian system. Then multi-symplectic Fourier pseudospectral method and multi-symplectic Preissmann scheme are employed to discretize the linear subproblem, respectively. And the nonlinear subsystem is solved by a symplectic scheme. Finally, a composition method is applied to obtain the final schemes for the NLD equation. We find that the two proposed schemes preserve the total symplecticity and can be solved explicitly. Numerical experiments are presented to show the effectiveness of the proposed methods.     

Multi-symplectic wavelet splitting method for the strongly coupled Schrodinger system

<正>We propose a multi-symplectic wavelet splitting method to solve the strongly coupled nonlinear Schrodinger equations.Based on its multi-symplectic formulation,the strongly coupled nonlinear Schr(o|¨)dinger equations can be split into one linear multi-symplectic subsystem and one nonlinear infinite-dimensional Hamiltonian subsystem.For the linear subsystem,the multi-symplectic wavelet collocation method and the symplectic Buler method are employed in spatial and temporal discretization,respectively.For the nonlinear subsystem,the mid-point symplectic scheme is used.Numerical simulations show the effectiveness of the proposed method during long-time numerical calculation.     

Splitting multisymplectic integrators for Maxwell’s equations

In the paper, we describe a novel kind of multisymplectic method for three-dimensional (3-D) Maxwell’s equations. Splitting the 3-D Maxwell’s equations into three local one-dimensional (LOD) equations, then applying a pair of symplectic Runge–Kutta methods to discretize each resulting LOD equation, it leads to splitting multisymplectic integrators. We say this kind of schemes to be LOD multisymplectic scheme (LOD-MS). The discrete conservation laws, convergence, dispersive relation, dissipation and stability are investigated for the schemes. Theoretical analysis shows that the schemes are unconditionally stable, non-dissipative, and of first order accuracy in time and second order accuracy in space. As a reduction, we also consider the application of LOD-MS to 2-D Maxwell’s equations. Numerical experiments match the theoretical results well. They illustrate that LOD-MS is not only efficient and simple in coding, but also has almost all the nature of multisymplectic integrators.     

Multi-symplectic Birkhoffian structure for PDEs with dissipation terms

A generalization of the multi-symplectic form for Hamiltonian systems to self-adjoint systems with dissipation terms is studied. These systems can be expressed as multi-symplectic Birkhoffian equations, which leads to a natural definition of Birkhoffian multi-symplectic structure. The concept of Birkhoffian multi-symplectic integrators for Birkhoffian PDEs is investigated. The Birkhoffian multi-symplectic structure is constructed by the continuous variational principle, and the Birkhoffian multi-symplectic integrator by the discrete variational principle. As an example, two Birkhoffian multi-symplectic integrators for the equation describing a linear damped string are given.     

Variational Formulation for the Multisymplectic Hamiltonian Systems

A variational formulation for the multisymplectic Hamiltonian systems is presented in this Letter. Using this variational formulation, we obtain multisymplectic integrators from a variational perspective. Numerical experiments are also reported.Mathematical Subject Classifications (2000). 70G50, 58Z05.     

Total Variation in Discrete Multisymplectic Field Theory and Multisymplectic-Energy-Momentum Integrators

A total variation calculus in discrete multisymplectic field theory is developed in this Letter. Using this discrete total variation calculus, we obtain multisymplectic-energy-momentum integrators. The multisymplectic discretization for the nonlinear Schrödinger equation is also presented.     

Difference Discrete Variational Principles,Euler—Lagrange Cohomology and Symplectic,Multisymplectic Structures Ⅲ：Application to Symplectic and Multisymplectic Algorithms

In the previous papers I and II,we have studied the difference discrete variational principle and the Euler-Lagrange cohomology in the framework of multi-parameter differential approach.We have gotten the difference discrete Euler-Lagrange equations and canonical ones for the difference discrete versions of classical mechanics and field theory as well as the difference discrete versions for the Euler-Lagrange cohomology and applied them to get the necessary and sufficient condition for the symplectic or multisymplectic geometry preserving properties in both the lagrangian and Hamiltonian formalisms.In this paper,we apply the difference discrete variational principle and Euler-Lagrange cohomological approach directly to the symplectic and multisymplectic algorithms.We will show that either Hamiltonian schemes of Lagrangian ones in both the symplectic and multisymplectic algorithms are variational integrators and their difference discrete symplectic structure-preserving properties can always be established not only in the solution space but also in the function space if and only if the related closed Euler-Lagrange cohomological conditions are satisfied.     

Multisymplectic Pseudospectral Discretizations for (3+1)-Dimensional Klein-Gordon Equation

We explore the multisymplectic Fourier pseudospectral discretizations for the (3+1)-dimensional Klein-Gordon equation in this paper. The corresponding multisymplectic conservation laws are derived. Two kinds of explicit symplectic integrators in time are also presented.     

Explicit multi-symplectic extended leap-frog methods for Hamiltonian wave equations

In this paper, we study the integration of Hamiltonian wave equations whose solutions have oscillatory behaviors in time and/or space. We are mainly concerned with the research for multi-symplectic extended Runge–Kutta–Nyström (ERKN) discretizations and the corresponding discrete conservation laws. We first show that the discretizations to the Hamiltonian wave equations using two symplectic ERKN methods in space and time respectively lead to an explicit multi-symplectic integrator (Eleap-frogI). Then we derive another multi-symplectic discretization using a symplectic ERKN method in time and a symplectic partitioned Runge–Kutta method, which is equivalent to the well-known Störmer–Verlet method in space (Eleap-frogII). These two new multi-symplectic schemes are extensions of the leap-frog method. The numerical stability and dispersive properties of the new schemes are analyzed. Numerical experiments with comparisons are presented, where the two new explicit multi-symplectic methods and the leap-frog method are applied to the linear wave equation and the Sine–Gordon equation. The numerical results confirm the superior performance and some significant advantages of our new integrators in the sense of structure preservation.