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

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

一维Gross-Pitaevskii方程的高阶紧致分裂步多辛格式

首先把一维Gross-Pitaevskli方程改写成多辛Hamiltonian系统的形式,把形式通过分裂变成2个子哈密尔顿系统.然后,对这些子系统用辛或者多辛算法进行离散.通过对子系统数值算法的不同组合方式,得到不同精度的具有多辛算法特征数值格式.这些格式不仅具有多辛格式、分裂步方法和高阶紧致格式的特征,而且是质量守恒的.数值实验验证了新格式的数值行为.     

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

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

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

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

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

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

两个三阶最优化力梯度辛积分器的对称组合

利用已存在的三阶最优化力梯度辛格式以对称组合方法获得两个新的四阶力梯度辛积分器.它们在求解摄动Kepler混沌问题的能量精度和一维定态Schrö,dinger方程的能量本征值精度方面比Forest-Ruth四阶非力梯度辛积分器要好得多,甚至还要明显优越于已有的四阶最优化力梯度辛积分器.     

啁啾激光场中HF分子的经典解离

应用经典理论并采用辛算法，计算了HF分子在啁啾激光作用下的经典解离.讨论了激光场中HF分子的解离概率随时间的变化以及考虑振动-转动能级跃迁对分子解离概率的影响.计算结果与理论分析相符，说明运用经典理论并采用辛格式计算双原子分子在激光作用下的经典解离是有效和可行的. 关键词： 啁啾激光场 经典解离 辛算法     

一维强场模型研究中的非齐线性正则方程的辛算法

就一维强场模型,采用对称差商代替空间变量的2阶偏导数,将含有SchrÖdinger方程的初边值问题离散成"非齐线性正则方程",它的齐方程的通解和非齐方程特解都由"辛变换生成",分别采用辛格式计算.采用这种辛算法和R-K法计算了一个数值例子,并与精确解作了比较.结果表明,经长时间计算后,辛算法保持解的固有特征,而R-K法则面目全非.     

气动声学中声传播问题的辛算法研究

引入辛算法对气动声学中的声传播问题进行了数值研究。采用Hamilton系统描述理想气体的声波方程,时间离散采用辛可分Runge-Kutta方法,空间离散采用近似解析方法,构造声波方程的保辛格式。将辛算法和有限差分算法分别在数值频散和计算效率等方面进行了对比分析,研究结果表明:辛算法能够有效地抑制数值频散,在计算效率方面具有明显的优越性。声传播特性模拟结果表明辛算法能够准确地模拟点源声辐射、声波干涉、反射及衍射现象。     

Klein-Gordon-Schrdinger方程的辛Fourier拟谱格式(英文)

主要讨论Klein-Gordon-Schrdinger方程的Fourier拟谱辛格式,包括中点公式和Strmer/Verlet格式.首先构造一个哈密尔顿方程,针对此哈密尔顿方程,在空间方向用Fourier拟谱离散得到一个有限维的哈密尔顿系统,对此有限维系统在时间方向用Strmer/Verlet方法离散得到KGS方程的完全显式的辛格式.中点格式虽然是隐式的但效率也很高,且具有质量守恒律.数值实验表明,辛格式能够在长时间内很好地模拟各类孤立波.     

Multisymplectic Geometry,Local Conservation Laws and a Multisymplectic Integrator for the Zakharov–Kuznetsov Equation

The multisymplectic geometry for the Zakharov–Kuznetsov equation is presented in this Letter. The multisymplectic form and the local energy and momentum conservation laws are derived directly from the variational principle. Based on the multisymplectic Hamiltonian formulation, we derive a 36-point multisymplectic integrator.     

Multisymplectic Integrator of the Zakharov System

A multisymplectic formulation for the Zakhaxov system is presented. The semi-explicit multisymplectic integrator of the formulation is constructed by means of the Euier-box scheme. Numerical results on simulating the propagation of one soliton and the collision of two solitons axe reported to illustrate the efficiency of the multisymplectic scheme.     

A new explicit multisymplectic integrator for the Kawahara-type equation

We derive a new multisymplectic integrator for the Kawahara-type equation which is a fully explicit scheme and thus needs less computation cost. Multisympecticity of such scheme guarantees the long-time numerical behaviors. Nu- merical experiments are presented to verify the accuracy of this scheme as well as the excellent performance on invariant preservation for three kinds of Kawahara-type equations.     

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

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

A variational approach to second-order multisymplectic field theory

This paper presents a geometric-variational approach to continuous and discrete second-order field theories following the methodology of [Marsden, Patrick, Shkoller, Comm. Math. Phys. 199 (1998) 351–395]. Staying entirely in the Lagrangian framework and letting Y denote the configuration fiber bundle, we show that both the multisymplectic structure on J³Y as well as the Noether theorem arise from the first variation of the action function. We generalize the multisymplectic form formula derived for first-order field theories in [Marsden, Patrick, Shkoller, Comm. Math. Phys. 199 (1998) 351–395], to the case of second-order field theories, and we apply our theory to the Camassa–Holm (CH) equation in both the continuous and discrete settings. Our discretization produces a multisymplectic-momentum integrator, a generalization of the Moser–Veselov rigid body algorithm to the setting of nonlinear PDEs with second-order Lagrangians.     

Variational methods, multisymplectic geometry and continuum mechanics

This paper presents a variational and multisymplectic formulation of both compressible and incompressible models of continuum mechanics on general Riemannian manifolds. A general formalism is developed for non-relativistic first-order multisymplectic field theories with constraints, such as the incompressibility constraint. The results obtained in this paper set the stage for multisymplectic reduction and for the further development of Veselov-type multisymplectic discretizations and numerical algorithms. The latter will be the subject of a companion paper.     

Total Variation and Multisymplectic Structure for CNLS System

The relation between the toal variation of classical field theory and the multisymplectic structure is shown. Then the multisymplectic structure and the corresponding multisymplectic conservation of the coupled nonlinear Schroedinger system are obtained directly from the variational principle.     

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.     

高阶非线性薛定谔方程的分步小波方法

用小波变换代替傅里叶变换解高阶非线性薛定谔方程,为高阶薛定谔方程的数值解提供了一种工具,提高了运算速度.本文分析了高阶非线性薛定谔方程分步解法的一般形式,选用Db10小波,得到了小波微分算子和色散算子对应的矩阵,得出了分步小波方法的算法公式.推导了色散算子和时域信号在小波域相乘的近似运算公式,说明了分步傅里叶方法比分步小波方法的复数乘法次数更多,同时说明了提高运算速度必须舍弃一定的运算准确度.最后以分步傅里叶方法为准,分析了分步小波方法的误差,结果表明:对于一阶孤子,分步小波方法与分步傅里叶方法间的相对误差在1.2%左右波动.