非线性薛定谔方程的平均向量场方法

利用平均向量场方法(AVF)对非线性薛定谔方程进行求解, 在理论上得到了一个保非线性薛定谔方程描述的系统能量守恒的AVF格式, 再分别用非线性薛定谔方程的AVF格式和辛格式数值模拟孤立波的演化行为, 并比较两个格式是否保系统能量守恒特性. 数值结果表明, AVF格式也能很好地模拟孤立波的演化行为,并且比辛格式更能保持系统的能量守恒.     

Cahn-Hilliard方程的高阶保能量散逸性方法

能量散逸性是物理和力学中某些微分方程一项重要的物理特性.构造精确地保持微分方程能量散逸性的数值格式对模拟具有能量散逸性的微分方程具有重要的意义.本文利用四阶平均向量场方法和傅里叶谱方法构造了Cahn-Hilliard方程高阶保能量散逸性格式.数值结果表明高阶保能量散逸性格式能很好地模拟Cahn-Hilliard方程在不同初始条件下解的行为,并且很好地保持了Cahn-Hilliard方程的能量散逸特性.     

光伏光折变晶体中孤立波的数值模拟

高斯光束在光伏光折变晶体中孤立波的演化满足傍轴方程.傍轴方程可以看作无限维Hamil-tonian系统并可以利用辛几何算法进行计算.数值结果表明外加电场和光伏场的强弱和入射高斯光束的振辐对形成稳定的孤立波有显著的影响.傍轴方程的辛几何差分格式能很好地模拟傍轴方程中孤立波的演化行为.     

带源项浅水波方程的高分辨率熵稳定格式

提出了一种求解带源项浅水波方程的熵稳定格式.新格式利用通量限制函数将一阶熵稳定格式和高阶熵守恒格式结合,具有熵守恒格式和熵稳定格式的优点:在解的光滑区域具有高精度,在解的间断区域避免了非物理现象的产生,同时可以准确地捕捉激波,从而达到高分辨率的效果.利用新格式计算了一维和二维的经典算例,数值结果表明,新格式是模拟带源项浅水波方程的理想方法.     

三阶Airy方程满足两个守恒律的非线性差分格式

近年来,学者们对发展型偏微分方程设计了一种能保持多个守恒律的数值方法,这类方法无论在解的精度还是长时间的数值模拟方面都表现出非常好的性质.将这类思想应用到三阶Airy方程,即三阶散射方程,对其设计了满足两个守恒律的非线性差分格式.该格式不仅计算数值解,同时计算数值能量,并且保证数值解和数值能量同时守恒.从数值结果可以看出,该格式在长时间的数值模拟中具有更好的保结构性质.     

有限振幅T-S波在非平行边界层中的非线性演化研究

研究对非平行边界层稳定性有重要影响的非线性演化问题,导出与其相应的抛物化稳定性方程组,发展了求解有限振幅T-S波的非线性演化的高效数值方法。这一数值方法包括预估-校正迭代求解各模态非线性方程并避免模态间的耦合,采用高阶紧致差分格式,满足正规化条件,确定不同模态非线性项表和数值稳定地作空间推进。通过给出T-S波不同的初始幅值,研究其非线性演化。算例与全Navier-Stokes方程的直接数值模拟(DNS)的结果作了比较。     

耦合Higgs方程和Maccari系统的行波解分支

利用动力系统方法,对耦合Higgs方程和Maccari系统的定性行为和行波解进行了研究.基于这种方法,给出了系统在不同参数条件下的相图,得到了包括孤立波解和周期波解在内的行波解.运用数值模拟的方法,对方程的光滑孤立波解和周期波解进行了数值模拟.获得的结果完善了相关文献已有的研究成果.     

浅水中度振幅孤立波解的分支

利用平面动力系统分支方法研究浅水中度振幅方程的定性行为和孤立波解.给出了系统在不同参数条件下的相图.获得了光滑孤立波、cuspon解和周期波解的隐式表达式.对方程的光滑孤立波解、cuspon解和周期波解进行了数值模拟.获得的结果完善了相关文献已有的结果.     

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

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

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

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

A new high-order energy-preserving scheme for the modified Korteweg-de Vries equation

In this paper, a new high-order energy-preserving scheme is proposed for the modified Korteweg-de Vries equation. The proposed scheme is constructed by using the Hamiltonian boundary value methods in time, and Fourier pseudospectral method in space. Exploiting this method, we get second-order and fourth-order energy-preserving integrators. The proposed schemes not only have high accuracy, but also exactly conserve the total mass and energy in the discrete level. Finally, single solitary wave and the interaction of two solitary waves are presented to illustrate the effectiveness of the proposed schemes.     

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 new efficient energy-preserving finite volume element scheme for the improved Boussinesq equation

In this paper, a new linearized energy-preserving Crank-Nicolson finite volume element scheme is derived for the improved Boussinesq equation. The fully discrete scheme can be shown to conserve both mass and energy in the discrete setting. It is proved that the scheme is uniquely solvable and convergent with the rate of order two in a discrete L₂ norm. At last, a series of numerical experiments on typical improved Boussinesq and Sine–Gordon equations are provided to verify our theoretical results and to show the efficiency and accuracy of the proposed scheme.     

求解洛伦兹力系统的一个任意高阶的能量守恒格式

本文旨在提出一种求解带电粒子系统的任意高阶能量守恒格式.在使用能量不变二次化方法将原始哈密顿能量泛函转化为一个二次形式后,辛龙格-库塔方法被用来构造了一种新的能量守恒格式来求解洛伦兹力系统.新格式不仅能保持能量守恒,而且可以达到任意高阶.通过与经典的Boris方法和另一个二阶能量守恒方法对比,数值实验验证了所提算法的显著优越性.     

Numerical study of Sivashinsky equation using a splitting scheme based on Crank‐Nicolson method

The mathematical modeling of a planar solid‐liquid interface in the solidification of a dilute binary alloy is formulating by one of nonintegrable, nonlinear evolution equation known as Sivashinsky equation. In the first part of this paper, the mathematical modeling of Sivashinsky equation is briefly discussed. Since, the exact solutions of this equation is yet unknown, obtaining its numerical solution plays an important role to simulate its behavior. Therefore, in the second part, a second‐order splitting finite difference scheme, based on Crank‐Nicolson method, is investigated to approximate the solution of the Sivashinsky equation with homogeneous boundary conditions. We prove the solvability of the present scheme and establish the error estimate of the numerical scheme.     

Geometric integration of the paraxial equation

Evolution of solitary waves in photovoltaic-photorefractive crystal satisfy the paraxial equation. The paraxial equation is transformed into the symplectic structure of the infinite dimensional Hamiltonian system. The symplectic structure of the paraxial equation is discretizated by the symplectic method. The corresponding symplectic scheme preserves conservation of discrete energy which reflects conservation of energy of the paraxial equation. The symplectic scheme is applied to simulate the solitary wave behaviors of the paraxial equation. Evolution of the solitary waves with the different applied electric field and the different photovoltaic fields are investigated.     

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.