广义Boussinesq方程的多辛方法

广义Boussinesq方程作为一类重要的非线性方程有着许多有趣的性质,基于Hamilton空间体系的多辛理论研究了广义Boussinesq方程的数值解法,构造了一种等价于多辛Box格式的新隐式多辛格式,该格式满足多辛守恒律、局部能量守恒律和局部动量守恒律．对广义Boussinesq方程孤子解的数值模拟结果表明,该多辛离散格式具有较好的长时间数值稳定性．     

二维非线性Schr(o)dinger方程的辛与多辛格式

对满足周期边界条件的二维非线性Schrödinger方程,运用中心差分对该方程进行空间离散, 得到一个有限维Hamilton系统,然后用隐式Euler中点格式进行时间离散得到其辛格式. 针对该方程的多辛形式, 运用有限体积法离散,得到一种直平行六面体上的中点型多辛格式. 用所构造的辛与多辛格式对二维非线性Schrödinger方程的平面波解和奇异解进行数值模拟,结果验证了所构 造格式的有效性. 最后, 根据计算结果,对两种格式进行了分析和比较.       

二维阻尼非线性sine-Gordon方程的共形多辛Fourier拟谱格式

为二维阻尼非线性sine-Gordon方程构造了一个新的共形多辛Fourier拟谱格式.基于原系统的共形多辛哈密尔顿形式,首先在时间和空间方向上分别用辛中点和Fourier拟谱方法进行离散,得到一个全离散格式.随后证明了构造的格式保持离散的共形多辛守恒律.最后数值实验验证了格式的有效性.     

2+1维sine-Gordon方程多辛格式的复合构造

以2+1维sine-Gordon方程为例,提出构造多辛算法的一种复合方法, 即: 独立地辛离散多辛偏微分方程的各个方向, 复合这些辛离散格式得到原偏微分方程的多辛算法. 通过这种复合方法, 可以构造任意阶精度的多辛算法,其中包括一些常用算法. 孤立子的数值模拟试验说明所给出的多辛算法是有效的.     

带乘性噪声的空间分数阶随机非线性Schrödinger方程的广义多辛算法

带乘性噪声的空间分数阶随机非线性Schrödinger方程是一类重要的方程,可应用于描述开放非局部量子系统的演化过程.该方程为一个无穷维分数阶随机Hamilton系统,且具有广义多辛结构和质量守恒的性质.针对该方程的广义多辛形式,在空间上采用拟谱方法离散分数阶微分算子,在时间上则采用隐式中点格式,构造出一类保持全局质量的广义多辛格式.对行波解和平面波解等进行数值模拟,结果验证了所构造格式的有效性和保结构性质,时间均方收敛阶约在0.5到1之间.     

Landau-Ginzburg-Higgs方程的多辛Runge-Kutta方法

非线性波动方程作为一类重要的数学物理方程吸引着众多的研究者,基于Hamilton空间体系的多辛理论研究了Landau-Ginzburg-Higgs方程的多辛算法,讨论了利用Runge-Kutta方法构造离散多辛格式的途径,并构造了一种典型的半隐式的多辛格式,该格式满足多辛守恒律、局部能量守恒律和局部动量守恒律.数值算例结果表明该多辛离散格式具有较好的长时间数值稳定性.     

一类强色散DGH方程的多辛普雷斯曼格式

DGH方程作为一类重要的非线性水波方程有着许多广泛的应用前景.基于Hamilton系统的多辛理论研究了一类强色散DGH方程的数值解法,利用多辛普雷斯曼方法构造了一种典型的半隐式的多辛格式.分析了该格式的局部能量和动量守恒律误差,并给出了数值算例.数值算例结果表明该多辛离散格式具有较好的长时间数值稳定性.     

一类DGH方程的多辛Fourier拟谱方法

DGH方程作为一类重要的非线性方程有着许多广泛的应用前景.通过正则变化,构造了DGH方程的多辛哈密尔顿系统.利用Fourier拟谱方法对此哈密尔顿系统进行数值离散,并构造了一种半隐式的多辛格式.数值算例结果表明该多辛离散格式具有较好的长时间数值稳定性.     

膜自由振动的多辛方法

基于Hamilton空间体系的多辛理论研究了膜自由振动问题,讨论了构造复合离散多辛格式的方法,并构造了一种典型的9×3点半隐式的多辛复合离散格式,该格式满足多辛守恒律、能量守恒律和动量守恒律．数值算例结果表明该多辛离散格式具有较好的长时间数值稳定性．     

ZK-BBM方程的多辛Preissmann格式

本文研究一类非线性ZK-BBM方程的初值问题.利用Hamilton系统的多辛Preissmann方法,获得ZK-BBM方程初值问题的数值结果,数值结果表明该多辛离散格式具有较好的长时间数值稳定性.     

非线性Pochhammer-Chree方程的多辛格式

提出非线性Pochhammer—Chree方程的多辛形式,进而得到一个等价于中心Preissmann积分的15点多辛格式．数值例子表明：多辛格式具有良好的长时间数值行为。     

A NEW MULTI-SYMPLECTIC SCHEME FOR NONLINEAR “GOOD” BOUSSINESQ EQUATION

The Hamiltonian formulations of the linear “good“ Boussinesq (L.G.B.) equation and the multi-symplectic formulation of the nonlinear “good“ Boussinesq (N.G.B.) equation are considered. For the multi-symplectic formulation, a new fifteen-point difference scheme which is equivalent to the multi-symplectic Preissmann integrator is derived. We also present numerical experiments, which show that the symplectic and multisymplectic schemes have excellent long-time numerical behavior.     

Symplectic and multi-symplectic wavelet collocation methods for two-dimensional Schrödinger equations

In this paper, we develop symplectic and multi-symplectic wavelet collocation methods to solve the two-dimensional nonlinear Schrödinger equation in wave propagation problems and the two-dimensional time-dependent linear Schrödinger equation in quantum physics. The Hamiltonian and the multi-symplectic formulations of each equation are considered. For both formulations, wavelet collocation method based on the autocorrelation function of Daubechies scaling functions is applied for spatial discretization and symplectic method is used for time integration. The conservation of energy and total norm is investigated. Combined with splitting scheme, splitting symplectic and multi-symplectic wavelet collocation methods are also constructed. Numerical experiments show the effectiveness of the proposed methods.     

Lattice Boltzmann model for two-dimensional generalized sine-Gordon equation

The nonlinear sine-Gordon equation arises in various problems in science and engineering. In this paper, we propose a numerical model based on lattice Boltmann method to obtain the numerical solutions of two-dimensional generalized sine-Gordon equation, including damped and undamped sine-Gordon equation. By choosing properly the conservation condition between the macroscopic quantity $u_t$ and the distribution functions and applying the Chapman-Enskog expansion, the governing equation is recovered correctly from the lattice Boltzmann equation. Moreover, the local equilibrium distribution function is obtained. The numerical results of the first three examples agree well with the analytic solutions, which indicates the lattice Boltzmann model is satisfactory and efficient. Numerical solutions for cases involving the most known from the bibliography line and ring solitons are given. Numerical experiments also show that the present scheme has a good long-time numerical behavior for the generalized sine-Gordon equation. Moreover, the model can also be applied to other two-dimensional nonlinear wave equations, such as nonlinear hyperbolic telegraph equation and Klein-Gordon equation.     

Explicit Multi-symplectic Method for a High Order Wave Equation of KdV Type

In this paper, we consider multi-symplectic Fourier pseudospectral method for a high order integrable equation of KdV type, which describes many important physical phenomena. The multi-symplectic structure are constructed for the equation, and the conservation laws of the continuous equation are presented. The multisymplectic discretization of each formulation is exemplified by the multi-symplectic Fourier pseudospectral scheme. The numerical experiments are given, and the results verify the efficiency of the Fourier pseudospectral method.     

Multilevel augmentation methods for solving the sine-Gordon equation

In this paper we develop the multilevel augmentation method for solving nonlinear operator equations of the second kind and apply it to solving the one-dimensional sine-Gordon equation. We first give a general setting of the multilevel augmentation method for solving the second kind nonlinear operator equations and prove that the multilevel augmentation method preserves the optimal convergence order of the projection method while reducing computational cost significantly. Then we describe the semi-discrete scheme and the fully-discrete scheme based on multiscale methods for solving the sine-Gordon equation, and apply the multilevel augmentation method to solving the discrete equation. A complete analysis for convergence order is proposed. Finally numerical experiments are presented to confirm the theoretical results and illustrate the efficiency of the method.