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

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

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

主要讨论Klein-Gordon-Sehrodinger方程的Fourier拟谱辛格式,包括中点公式和Stormer/Vedet格式.首先构造一个哈密尔顿方程,针对此哈密尔顿方程,在空间方向用Fourier拟谱离散得到一个有限维的哈密尔顿系统,对此有限维系统在时间方向用St(o)rmer/Verlet方法离散得到KGS...     

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

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

孤立波方程的保结构算法

讨论了孤立波方程的保结构差分算法，以一些经典的孤立波方程为例，如KdV，sine-Gordon，K-P方程，给出了它们的辛和多辛结构，说明辛和多辛算法的可适用性．提出局部守恒算法和广义保结构算法的概念，它们是保结构算法的概念自然推广．还给出一种能系统构造局部守恒格式的复合方法．数值例子说明，保结构数值能很好模拟各种孤立波的演化。     

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

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

耦合非线性Schr(o)dinger系统的多辛差分格式

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

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

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

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

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

求解耦合Schrdinger方程组的数值逼近算法(英文)

对耦合Schrdinger方程组提出一个非耦合的线性化差分格式并对其进行分析.证明格式保持原方程组的守恒律,在先验估计的基础上证明格式依L2模的绝对稳定性和无条件二阶收敛性.对孤波碰撞的各种现象进行模拟.     

计算物理第23卷（2006）总目录

第1期计算中子增殖率的时间强迫碰撞方法……………………………王瑞宏邓力许海燕裴鹿成(1)Navier-Stokes方程的线性化微分求积法(英)……………………………………………李晨吴雄华(18)统一坐标系下多介质流体力学计算的一种方法…………………………………………………贾鹏彦(19)SRLW方程的多辛Fourier拟谱格式及其守恒律…………………孔令华曾文平刘儒勋孔令健(25)一类正交增强的阶谱六面体矢量单元构造……………………………………………李江海孙秦(32)SALE速度重映算法的改进……………………………………………………熊俊周海…     

A High-Order Conservative Numerical Method for Gross-Pitaevskii Equation with Time-Varying Coefficients in Modeling BEC

We propose a high-order conservative method for the nonlinear Schrodinger/Gross-Pitaevskii equation with time-varying coefficients in modeling Bose-Einstein condensation(BEC). This scheme combined with the sixth-order compact finite difference method and the fourth-order average vector field method, finely describes the condensate wave function and physical characteristics in some small potential wells. Numerical experiments are presented to demonstrate that our numerical scheme is efficient by the comparison with the Fourier pseudo-spectral method.Moreover, it preserves several conservation laws well and even exactly under some specific conditions.     

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.     

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.     

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

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

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

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

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.     

Multisymplectic Euler Box Scheme for the KdV Equation

We investigate the multisymplectic Euler box scheme for the Korteweg-de Vries （KdV） equation. A new completely explicit six-point scheme is derived. Numerical experiments of the new scheme with comparisons to the Zabusky-Kruskal scheme, the multisymplectic 12-point scheme, the narrow box scheme and the spectral method are made to show nice numerical stability and ability to preserve the integral invariant for long-time integration.     

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.     

A pseudo-spectral algorithm and test cases for the numerical solution of the two-dimensional rotating Green–Naghdi shallow water equations

A pseudo-spectral algorithm is presented for the solution of the rotating Green–Naghdi shallow water equations in two spatial dimensions. The equations are first written in vorticity–divergence form, in order to exploit the fact that time-derivatives then appear implicitly in the divergence equation only. A nonlinear equation must then be solved at each time-step in order to determine the divergence tendency. The nonlinear equation is solved by means of a simultaneous iteration in spectral space to determine each Fourier component. The key to the rapid convergence of the iteration is the use of a good initial guess for the divergence tendency, which is obtained from polynomial extrapolation of the solution obtained at previous time-levels. The algorithm is therefore best suited to be used with a standard multi-step time-stepping scheme (e.g. leap-frog).