The multi-symplectic Fourier pseudospectral method for solving two-dimensional Hamiltonian PDEs

In this paper, the multi-symplectic Fourier pseudospectral (MSFP) method is generalized to solve two-dimensional Hamiltonian PDEs with periodic boundary conditions. Using the Fourier pseudospectral method in the space of the two-dimensional Hamiltonian PDE (2D-HPDE), the semi-discrete system obtained is proved to have semi-discrete multi-symplectic conservation laws and a global symplecticity conservation law. Then, the implicit midpoint rule is employed for time integration to obtain the MSFP method for the 2D-HPDE. The fully discrete multi-symplectic conservation laws are also obtained. In addition, the proposed method is applied to solve the Zakharov-Kuznetsov (ZK) equation and the Kadomtsev-Petviashvili (KP) equation. Numerical experiments on soliton solutions of the ZK equation and the KP equation show the high accuracy and effectiveness of the proposed method.     

Multi-Symplectic Fourier Pseudospectral Method for a Higher Order Wave Equation of KdV Type

The higher order wave equation of KdV type, which describes many important physical phenomena, has been investigated widely in last several decades. In this work, multi-symplectic formulations for the higher order wave equation of KdV type are presented, and the local conservation laws are shown to correspond to certain well-known Hamiltonian functionals. The multi-symplectic discretization of each formulation is calculated by the multi-symplectic Fourier pseudospectral scheme. Numerical experiments are carried out, which verify the efficiency of the Fourier pseudospectral method.     

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.     

Multi-symplectic methods for the Ito-type coupled KdV equation

In this paper, we find that the Ito-type coupled KdV equation can be written as a multi-symplectic Hamiltonian partial differential equation (PDE). Then, multi-symplectic Fourier pseudospectral method and multi-symlpectic wavelet collocation method are constructed for this equation. In the numerical experiments, we show the effectiveness of the proposed methods. Some comparisons between the proposed methods are also made with respect to global conservation properties.     

Zhiber-Shabat方程的多辛几何和多辛隐式格式研究

Zhiber-Shabat方程,描述许多重要的物理现象,是一类重要的非线性方程,有着许多广泛的应用前景.本文给出Zhiber-Shabat方程的多辛几何结构和多辛Fourier拟谱方法.数值算例结果表明多辛离散格式具有较好的长时间的数值稳定性.     

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.     

GSDBM方程的多辛几何和多辛Preissman格式研究（英文）

The multi-symplectic geometry for the GSDBM equation is presented in this paper. The multi-symplectic formulations for the GSDBM equation are presented and the local conservation laws are shown to correspond to certain well-known Hamiltonian functionals. The multi-symplectic discretization of each formulation is exemplified by the multisymplectic Preissmann scheme. The numerical experiments are given, and the results verify the efficiency of the Preissmann scheme.     

Multi-symplectic Geometry and Preissmann Scheme for GSDBM Equation

The multi-symplectic geometry for the GSDBM equation is presented in this paper. The multi-symplectic formulations for the GSDBM equation are presented and the local conservation laws are shown to correspond to certain well-known Hamiltonian func-tionals. The multi-symplectic discretization of each formulation is exemplified by the multi-symplectic Preissmann scheme. The numerical experiments are given, and the results verify the efficiency of the Preissmann scheme.     

Geometric numerical schemes for the KdV equation

Geometric discretizations that preserve certain Hamiltonian structures at the discrete level has been proven to enhance the accuracy of numerical schemes. In particular, numerous symplectic and multi-symplectic schemes have been proposed to solve numerically the celebrated Korteweg-de Vries equation. In this work, we show that geometrical schemes are as much robust and accurate as Fourier-type pseudospectral methods for computing the long-time KdV dynamics, and thus more suitable to model complex nonlinear wave phenomena.     

不可压饱和多孔弹性杆动力响应的多辛方法

研究了不可压饱和多孔弹性杆的一维动力响应问题.基于多孔介质理论,在流相和固相微观不可压、固相骨架小变形的假定下,建立了不可压流体饱和多孔弹性杆一维轴向动力响应的数学模型.利用Hamilton空间体系的多辛理论,构造了不可压饱和多孔弹性杆轴向振动方程的多辛形式及其多种局部守恒律.采用中点Box离散方法得到轴向振动方程的多辛离散格式和局部能量守恒律以及局部动量守恒律的离散格式;数值模拟了不可压饱和多孔弹性杆的轴向振动过程,记录了每一时间步的局部能量数值误差和局部动量数值误差.结果表明,已构造的多辛离散格式具有很高的精确性和较长时间的数值稳定性,这为解决饱和多孔介质的动力响应问题提供了新的途径.     

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

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

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

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

High Order Energy-Preserving Method of the "Good" Boussinesq Equation

The fourth order average vector field (AVF) method is applied to solve the "Good" Boussinesq equation. The semi-discrete system of the "good" Boussinesq equation obtained by the pseudo-spectral method in spatial variable, which is a classical finite dimensional Hamiltonian system, is discretized by the fourth order average vector field method. Thus, a new high order energy conservation scheme of the "good" Boussinesq equation is obtained. Numerical experiments confirm that the new high order scheme can preserve the discrete energy of the "good" Boussinesq equation exactly and simulate evolution of different solitary waves well.     

A Linearly-Implicit Energy-Preserving Algorithm for the Two-Dimensional Space-Fractional Nonlinear Schrödinger Equation Based on the SAV Approach

The main objective of this paper is to present an efficient structure-preserving scheme, which is based on the idea of the scalar auxiliary variable approach, for solving the two-dimensional space-fractional nonlinear Schrödinger equation. First, we reformulate the equation as an canonical Hamiltonian system, and obtain a new equivalent system via introducing a scalar variable. Then, we construct a semi-discrete energy-preserving scheme by using the Fourier pseudo-spectral method to discretize the equivalent system in space direction. After that, applying the Crank-Nicolson method on the temporal direction gives a linearly-implicit scheme in the fully-discrete version. As expected, the proposed scheme can preserve the energy exactly and more efficient in the sense that only decoupled equations with constant coefficients need to be solved at each time step. Finally, numerical experiments are provided to demonstrate the efficiency and conservation of the scheme.     

Fourier spectral method with an adaptive time strategy for nonlinear fractional Schrödinger equation

In this paper, a Fourier spectral method with an adaptive time step strategy is proposed to solve the fractional nonlinear Schrödinger (FNLS) equation with periodic initial value problem. First, we prove the conservation law of the mass and the energy for the semi-discrete Fourier spectral scheme. Second, the error estimation of the semi-discrete scheme is given in the relevant fractional Sobolev space. Then, an adaptive time-step strategy is designed to reduce central processing unit (CPU) time. Finally, the numerical experiments for the one-, two- and three-dimensional FNLSs, show that the adaptive strategy, compared to the constant time step, can reduce the CPU-time by almost half.     

A conservative local discontinuous Galerkin method for the solution of nonlinear Schrdinger equation in two dimensions

In this study, we present a conservative local discontinuous Galerkin(LDG) method for numerically solving the two-dimensional nonlinear Schrdinger(NLS) equation. The NLS equation is rewritten as a firstorder system and then we construct the LDG formulation with appropriate numerical flux. The mass and energy conserving laws for the semi-discrete formulation can be proved based on different choices of numerical fluxes such as the central, alternative and upwind-based flux. We will propose two kinds of time discretization methods for the semi-discrete formulation. One is based on Crank-Nicolson method and can be proved to preserve the discrete mass and energy conservation. The other one is Krylov implicit integration factor(IIF) method which demands much less computational effort. Various numerical experiments are presented to demonstrate the conservation law of mass and energy, the optimal rates of convergence, and the blow-up phenomenon.     

Quadratic invariants and multi-symplecticity of partitioned Runge-Kutta methods for Hamiltonian PDEs

In this paper, we study the preservation of quadratic conservation laws of Runge-Kutta methods and partitioned Runge-Kutta methods for Hamiltonian PDEs and establish the relation between multi-symplecticity of Runge-Kutta method and its quadratic conservation laws. For Schrödinger equations and Dirac equations, it reveals that multi-symplectic Runge-Kutta methods applied to equations with appropriate boundary conditions can preserve the global norm conservation and the global charge conservation, respectively.     

Error estimates of structure-preserving Fourier pseudospectral methods for the fractional Schrödinger equation

This paper gives a rigorous error analysis of the multisymplectic Fourier pseudospectral method for the nonlinear fractional Schrödinger equation. The method preserves some intrinsic structure properties including the generalized multisymplectic conservation law. By rewriting it in a matrix form similar to that in the finite difference method, the method is shown to be convergent in the discrete L² norm with the second-order accuracy in time and spectral accuracy in space. The key techniques in the analysis include the discrete energy method, cutoff of the nonlinearity, and a posterior bound of numerical solutions by using the inverse inequality. In a similar line, the convergence result for the symplectic Fourier pseudospectral method can also be established. Moreover, the errors in the local and global energy conservation laws of discrete systems are also investigated. Numerical tests are performed to confirm the theoretical results.     

Sine-Gordon方程的多辛Leap-frog格式

非线性发展方程由于具有多种形式的解析解而吸引着众多的研究者,借助多辛保结构理论研究了Sine-Gordon方程的多辛算法．利用Hamilton变分原理,构造出了sine-Gordon方程的多辛格式;采用显辛离散方法得到了Leap-frog多辛离散格式,该格式满足多辛守恒律;数值结果表明leap-frog多辛离散格式能够精确地模拟sine-Gordon方程的孤子解和周期解,模拟结果证实了该离散格式具有良好的数值稳定性．     

求解带刚性源项标量双曲型守恒律方程的保有界WCNS格式

带刚性源项的双曲守恒律方程是很多物理问题,特别是化学反应流的数学模型.本文考虑带刚性源项的标量双曲型守恒律方程,通过时空分离的方式,发展了一类保有界的WCNS格式.对于空间离散,我们将参数化的通量限制器推广到WCNS框架,使得方程对流项离散后满足极值原理.对于时间离散,我们将半离散的WCNS改写成指数形式,采用三阶修正指数型Runge-Kutta格式来控制方程的刚性,保持数值解的界.可以证明,本文格式对带刚性源项的一维标量守恒律方程具有保有界性和弱渐近保持性.数值试验验证了方法的有效性.