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

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

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.     

Hamiltonian spatial structure for three-dimensional water waves in a moving frame of reference

Summary The governing equations for three-dimensional time-dependent water waves in a moving frame of reference are reformulated in terms of the energy and momentum flux. The novelty of this approach is that time-independent motions of the system—that is, motions that are steady in a moving frame of reference—satisfy a partial differential equation, which is shown to be Hamiltonian. The theory of Hamiltonian evolution equations (canonical variables, Poisson brackets, symplectic form, conservation laws) is applied to the spatial Hamiltonian system derived for pure gravity waves. The addition of surface tension changes the spatial Hamiltonian structure in such a way that the symplectic operator becomes degenerate, and the properties of this generalized Hamiltonian system are also studied. Hamiltonian bifurcation theory is applied to the linear spatial Hamiltonian system for capillary-gravity waves, showing how new waves can be found in this framework.     

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.     

A collocation method with cubic B-splines for solving the MRLW equation

The modified regularized long wave (MRLW) equation is solved numerically by collocation method using cubic B-splines finite element. A linear stability analysis of the scheme is shown to be marginally stable. Three invariants of motion are evaluated to determine the conservation properties of the algorithm, also the numerical scheme leads to accurate and efficient results. Moreover, interaction of two and three solitary waves are studied through computer simulation and the development of the Maxwellian initial condition into solitary waves is also shown.     

Symplectic structure-preserving integrators for the two-dimensional Gross-Pitaevskii equation for BEC

Symplectic integrators have been developed for solving the two-dimensional Gross-Pitaevskii equation. The equation is transformed into a Hamiltonian form with symplectic structure. Then, symplectic integrators, including the midpoint rule, and a splitting symplectic scheme are developed for treating this equation. It is shown that the proposed codes fulfill the discrete charge conservation law. Furthermore, the global error of the numerical solution is theoretically estimated. The theoretical analysis is supported by some numerical simulations.     

A new lattice hierarchy: Hamiltonian structures,symplectic map and N-fold Darboux transformation

A new lattice hierarchy is constructed from a discrete matrix spectral problem. By the Tu scheme technique, the associated Hamiltonian structures and infinitely many conservation laws of this hierarchy are derived. Then a symplectic map is proposed based on the Lax pair and the adjoint Lax pair. Furthermore, the N-fold Darboux transformation and explicitly exact solutions of the first two equations in the hierarchy are investigated. Finally, the density profiles of these exact solutions are presented to illustrate the overtaking collisions of solitary waves.     

A new conservative fourth‐order accurate difference scheme for solving a model of nonlinear dispersive equations

This article is devoted to the study of a nonlinear conservative fourth‐order difference scheme for a model of nonlinear dispersive equations that is governed by the RLW‐KdV equation. The existence of the approximate solution and the convergence of the difference scheme are proved, by using the energy method. In addition, the convergent order in maximum norm is 2 in temporal direction and 4 in spatial direction. The unconditional stability as well as uniqueness of the difference scheme is also derived. An application on the RLW and MRLW equations is discussed numerically in details. Furthermore, interaction of solitary waves with different amplitudes are shown. The 3 invariants of the motion are evaluated to determine the conservation proprieties of the system. The temporal evaluation of a Maxwellian initial pulse is then studied. Some numerical examples are given to validate the theoretical results.     

Cosine expansion‐based differential quadrature algorithm for numerical solution of the RLW equation

The differential quadrature method based on cosine expansion is applied to obtain numerical solutions of the RLW equation. The propagation of single solitary wave is studied to validate the efficiency of the algorithm. Then, test problems including interaction of two and three solitary waves, undulation, and evolution of solitary waves are implemented. Solutions are compared with earlier results. Discrete conservation quantities are computed for test experiments. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

Poisson系统的辛结构

当Poisson系统中的Poisson矩阵是非常数时,经典的辛方法如辛Runge-Kutta方法,生成函数法一般不能保持Poisson系统的Poisson结构,利用非线性变换可把非常数Poisson结构转化成辛结构,然后任意阶的辛方法可以长时间计算Poisson系统的辛结构．自由刚体问题中Euler方程被转换成辛结构并用辛中点格式进行数值求解,数值结果给出了这种非线性变换的有效性．     

A meshless method using the radial basis functions for numerical solution of the regularized long wave equation

This article discusses on the solution of the regularized long wave (RLW) equation, which is introduced to describe the development of the undular bore, has been used for modeling in many branches of science and engineering. A numerical method is presented to solve the RLW equation. The main idea behind this numerical simulation is to use the collocation and approximating the solution by radial basis functions (RBFs). To avoid solving the nonlinear system, a predictor‐corrector scheme is proposed. Several test problems are given to validate the new technique. The numerical simulation, includes the propagation of a solitary wave, interaction of two positive solitary waves, interaction of a positive and a negative solitary wave, the evaluation of Maxwellian pulse into stable solitary waves and the development of an undular bore. The three invariants of the motion are calculated to determine the conservation properties of the algorithm. The results of numerical experiments are compared with analytical solution and with those of other recently published methods to confirm the accuracy and efficiency of the presented scheme.© 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010     

Quadratic trigonometric b-spline galerkin methods for the regularized long wave equation

In this study, a numerical solution of the Regularized Long Wave (RLW) equation is obtained using Galerkin finite element method, based on two and three steps Adams Moulton method for the time integration and quadratic trigonometric B-spline functions for the space integration. After two different linearization techniques are applied, the proposed algorithms are tested on the problems of propagation of a solitary wave and interaction of two solitary waves. For the first test problem, the rate of convergence and the running time of the proposed algorithms are computed and the error norm $L_{\infty }$ is used to measure the differences between exact and numerical solutions. The three conservation quantities of the motion are calculated to determine the conservation properties of the proposed algorithms for both of the test problems.     

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.     

Navier Stokes model of solitary wave collision

Wave collision and its interaction characteristics is one of the important challenges in coastal engineering. This article concerns the collision of solitary waves over a horizontal bottom considering unsteady, incompressible viscous flow with free surface. The method solves the two dimensional Naiver–Stokes equations for conservation of momentum, continuity equation, and full nonlinear kinematic free-surface equation for Newtonian fluids, as the governing equations in a vertical plan. A mapping was developed to trace the deformed free surface encountered during wave propagation, transforms and interaction by transferring the governing equations from the physical domain to a computational domain. Also a numerical scheme is developed using finite element modeling technique in order to predict the solitary wave collision. Consequently results compared with other researches and show the inelastic behavior of solitary wave collision.     

Numerical study using ADM for the modified regularized long wave equation

The modified regularized long wave (MRLW) equation is solved numerically by Adomian decomposition method (ADM) with some initial conditions. The method leads to high accurate and efficient results. Three polynomial invariant conditions are evaluated to determine the conservation properties of the problem. The convergence of Adomian decomposition method applied to the MRLW equation is proved. Moreover, the interaction of solitons and the development of the Maxwellian initial condition into solitary waves are considered.     

一维非线性周期结构中弹性波传播的辛数学方法

利用辛数学方法分析了质量-弹簧非线性周期结构链中弹性波的传播问题．首先利用能量方法得到频域动力方程,随后通过小量变换将非线性动力方程线性化,得到辛矩阵,进而通过求解辛矩阵的本征值问题来研究波的传播性能．质量-弹簧模型中的弹簧刚度非线性对结构链的传播特性影响很大,研究发现非线性明显改变了周期结构的传播性能,而且不同于线性结构,非线性结构的传播特性与入射波强度有关．数值算例表明随着非线性强度及入射波强度的增大,传播通带宽度逐渐减小,禁带宽度逐渐增大．当入射波强度增大到一定值时,弹性波无法在结构中进行传播．与一般递归方法的比较分析,验证了辛数学方法在非线性周期结构波传播问题中的有效性与优越性．     

Scattering of solitons for Dirac equation coupled to a particle

We establish soliton-like asymptotics for finite energy solutions to the Dirac equation coupled to a relativistic particle. Any solution with initial state close to the solitary manifold, converges in long time limit to a sum of traveling wave and outgoing free wave. The convergence holds in global energy norm. The proof uses spectral theory and symplectic projection onto solitary manifold in the Hilbert phase space.     

Two Novel Classes of Arbitrary High-Order Structure-Preserving Algorithms for Canonical Hamiltonian Systems

In this paper, we systematically construct two classes of structure-preserving schemes with arbitrary order of accuracy for canonical Hamiltonian systems. The one class is the symplectic scheme, which contains two new families of parameterized symplectic schemes that are derived by basing on the generating function method and the symmetric composition method, respectively. Each member in these schemes is symplectic for any fixed parameter. A more general form of generating functions is introduced, which generalizes the three classical generating functions that are widely used to construct symplectic algorithms. The other class is a novel family of energy and quadratic invariants preserving schemes, which is devised by adjusting the parameter in parameterized symplectic schemes to guarantee energy conservation at each time step. The existence of the solutions of these schemes is verified. Numerical experiments demonstrate the theoretical analysis and conservation of the proposed schemes.     

A Note on Conservation Laws of Symplectic Difference Schemes for Hamiltonian Systems

In this paper we consider the necessary conditions of conservation laws of symplectic difference schemes for Hamiltonian systems and give an example which shows that there does not exist any centered symplectic difference scheme which preserves all Hamiltonian energy.