Efficient numerical method for linear stability analysis of solitary waves

In this paper we use a numerical relaxation algorithm to improve and generalize the obtainment of the perturbation eigenstates of nonlinear systems. As a model problem we consider the linear stability analysis of the vortex eigenstates of the cubic–quintic nonlinear Schrödinger equation. It is shown by numerical calculations that the relaxation algorithm permits accurate tracing of complex perturbation eigenvalues.     

Orbital stability of negative solitary waves

The generalized regularized long-wave equation admits a family of negative solitary waves. The stability of these waves is investigated by numerical simulation using a spectral discretization.     

Modeling nonlinear waves in a numerical wave tank with localized meshless RBF method

This paper presents the numerical simulation of free surface waves in a 2D domain which represents a wave tank, using a localized approach of the meshless radial basis function (RBF) method. Instead of global collocation, the local approach breaks down the problem domain into subdomains, leading to a sparse global system matrix which is particularly advantageous in tackling the time consuming simulation process. Mixed Eulerian–Lagrangian approach is adopted for free surface tracking and fourth order Adams–Bashforth–Moulton scheme (ABM4) is used for time stepping. Both linear and nonlinear Stokes waves are simulated and compared to analytical solutions.     

Solitary waves,homoclinic breather waves and rogue waves of the (3+1)-dimensional Hirota bilinear equation

In this paper, the $(3+1)$-dimensional Hirota bilinear equation is investigated, which can be used to describe the nonlinear dynamic behavior in physics. By using the Bell polynomials, the bilinear form of the equation is derived in a very natural way. Based on the resulting bilinear form, its $N$-solitary waves are further obtained by using the Hirota’s bilinear theory. Finally, by using the Homoclinic test method, we obtain its rational breather wave and rogue wave solutions, respectively. In order to better understand the dynamical behaviors of the equation, some graphical analyses are discussed for these exact solutions.     

Convergence conditions for iterative methods seeking multi-component solitary waves with prescribed quadratic conserved quantities

We obtain linearized (i.e., non-global) convergence conditions for iterative methods that seek solitary waves with prescribed values of quadratic conserved quantities of multi-component Hamiltonian nonlinear wave equations. These conditions extend the ones found for single-component solitary waves in a recent publication by Yang and the present author. We also show that, and why, these convergence conditions coincide with dynamical stability conditions for ground-state solitary waves.Notably, our analysis applies regardless of whether the number of quadratic conserved quantities, s, equals or is less than the number of equations, S. To illustrate the situation when s < S, we use one of our iterative methods to find ground-state solitary waves in spin-1 Bose-Einstein condensates in a magnetic field (s = 2, S = 3).     

An automated tanh-function method for finding solitary wave solutions to non-linear evolution equations

The tanh-function method for finding explicit travelling solitary wave solutions to non-linear evolution equations is described. The method is usually extremely tedious to use by hand. We present a Mathematica package ATFM that deals with the tedious algebra and outputs directly the required solutions. The use of the package is illustrated by applying it to a variety of equations; not only are previously known solutions recovered but in some cases more general forms of solution are obtained.     

A lumped Galerkin method for the numerical solution of the modified equal-width wave equation using quadratic B-splines

The numerical solution of the one-dimensional modified equal width wave (MEW) equation is obtained by using a lumped Galerkin method based on quadratic B-spline finite elements. The motion of a single solitary wave and the interaction of two solitary waves are studied. The numerical results obtained show that the present method is a remarkably successful numerical technique for solving the MEW equation. A linear stability analysis of the scheme is also investigated.     

非均匀缓变折射率平板波导放大器中的畸形波

以缓变波导中光束传播的非线性传输方程为研究对象,研究了非均匀缓变折射率平板波导放大器中畸形波的非线性动力学性质.通过相似变换和直接假设,构建出带有自由函数的一阶精确畸形波解.在此基础上,针对不同类型的自由函数,通过数值模拟得到了不同畸形波的波形图,对于描述光纤中出现的一些物理现象具有重要的意义.     

The solitary wave solution of the two-dimensional regularized long-wave equation in fluids and plasmas

This paper investigates the solitary wave solutions of the two-dimensional regularized long-wave equation which is arising in the investigation of the Rossby waves in rotating flows and the drift waves in plasmas. The main idea behind the numerical solution is to use a combination of boundary knot method and the analog equation method. The boundary knot method is a meshless boundary-type radial basis function collocation technique. In contrast with the method of fundamental solution, the boundary knot method uses the non-singular general solution instead of the singular fundamental solution to obtain the homogeneous solution. Similar to method of fundamental solution, the radial basis function is employed to approximate the particular solution via the dual reciprocity principle. In the current paper, we applied the idea of analog equation method. According to the analog equation method, the nonlinear governing operator is replaced by an equivalent nonhomogeneous linear one with known fundamental solution and under the same boundary conditions. Furthermore, in order to show the efficiency and accuracy of the proposed method, the present work is compared with finite difference scheme. The new method is analyzed for the local truncation error and the conservation properties. The results of several numerical experiments are given for both the single and double-soliton waves.     

Solitary wave solutions for a higher order nonlinear Schrödinger equation

We consider a higher order nonlinear Schrödinger equation with third- and fourth-order dispersions, cubic–quintic nonlinearities, self steepening, and self-frequency shift effects. This model governs the propagation of femtosecond light pulses in optical fibers. In this paper, we investigate general analytic solitary wave solutions and derive explicit bright and dark solitons for the considered model. The derived analytical dark and bright wave solutions are expressed in terms of the model coefficients. These exact solutions are useful to understand the mechanism of the complicated nonlinear physical phenomena which are related to wave propagation in a higher-order nonlinear and dispersive Schrödinger system.     

三分子自催化反应扩散系统中化学波的二维格子Boltzmann模拟

建立二维区域带平方衰减反应的三分子自催化反应扩散系统的数学模型,利用浓度分布的Chapmann-Enskogz展开及多尺度技术,给出基于格子Boltzmann模型的二维三分子自催化反应扩散系统的数值求解法,得到半开放自催化反应扩散系统在反应与扩散机制同时作用所产生化学波的过程和浓度空间分布值.数值结果表明本文提供的求解化学波现象的方法是有效的.     

带强迫项的变系数KdV方程的多孤立波解及其应用

本文利用改进的齐次平衡法,首先得到了带强迫项的变系数KdV方程的多孤立波解,然后借助此解得到了强迫KdV方程的多孤立波解.最后作为应用例子,利用图形分析方法分析了Rossby孤立波的相互作用,指出了影响Rossby孤立波相对幅度、相位、传播方向及平衡位置的主要原因.     

Modelling solitary waves of a fifth-order non-linear wave equation

A numerical solution of a fifth-order non-linear dispersive wave equation is set up using collocation of seventh-order B-spline interpolation functions over finite elements. A linear stability analysis shows that this numerical scheme, based on a Crank–Nicolson approximation in time, is unconditionally stable. The method is used to model the behaviour of solitary waves.     

Existence and numerical simulation of solutions for nonlinear fractional pantograph equations

In this paper a method based on Sinc approximation is developed for the numerical solution of a nonlinear fractional pantograph equation. In order to use Sinc approximation, the problem needs to have an analytic solution. So we investigated the existence and uniqueness of analytic solutions in the proposed domain. Single and double exponential transformations are used to approximate the solution. Some test problems are given to demonstrate the efficiency and applicability of the methods. The results are compared with some other existing numerical methods to show the performance and good accuracy of the methods. The numerical orders of convergence show that the method has exponential rate of convergence.     

Visual simulation of shallow-water waves

A commodity-type graphics card (GPU) is used to simulate nonlinear water waves described by a system of balance laws called the shallow-water system. To solve this hyperbolic system we use explicit high-resolution central-upwind schemes, which are particularly well suited for exploiting the parallel processing power of the GPU. In fact, simulations on the GPU are found to run 15–30 times faster than on a CPU. The simulated cases involve dry-bed zones and nontrivial bottom topographies, which are real challenges to the robustness and accuracy of the discretization.     

数值鲁棒的非线性滤波新方法

应用多维情形的二阶插值公式构造新型非线性滤波器。该滤波器不需非线性函数的偏导计算，便能代替常规的扩展卡尔曼滤波器，并有滤波精度高、数值计算稳定和适用范围宽等优点。仿真实例表明新滤波器具有较高的性能。     

A numerical method for solving a nonlinear 2-D optimal control problem with the classical diffusion equation

This paper presents a numerical solution for solving a nonlinear 2-D optimal control problem (2DOP). The performance index of a nonlinear 2DOP is described with a state and a control function. Furthermore, dynamic constraint of the system is given by a classical diffusion equation. It is preferred to use the Ritz method for finding the numerical solution of the problem. The method is based upon the Legendre polynomial basis. By using this method, the given optimisation nonlinear 2DOP reduces to the problem of solving a system of algebraic equations. The benefit of the method is that it provides greater flexibility in which the given initial and boundary conditions of the problem are imposed. Moreover, compared with the eigenfunction method, the satisfactory results are obtained only in a small number of polynomials order. This numerical approach is applicable and effective for such a kind of nonlinear 2DOP. The convergence of the method is extensively discussed and finally two illustrative examples are included to observe the validity and applicability of the new technique developed in the current work.     

A quadrature-based method of moments for nonlinear filtering

According to the nonlinear filtering theory, optimal estimates of a general continuous-discrete nonlinear filtering problem can be obtained by solving the Fokker-Planck equation, coupled with a Bayesian update rule. This procedure does not rely on linearizations of the dynamical and/or measurement models. However, the lack of fast and efficient methods for solving the Fokker-Planck equation presents challenges in real time nonlinear filtering problems. In this paper, a direct quadrature method of moments is introduced to solve the Fokker-Planck equation efficiently and accurately. This approach involves representation of the state conditional probability density function in terms of a finite collection of Dirac delta functions. The weights and locations (abscissas) in this representation are determined by moment constraints and modified using the Bayes’ rule according to measurement updates. As demonstrated by numerical examples, this approach appears to be promising in the field of nonlinear filtering.     

Particle swarm algorithm for solving systems of nonlinear equations

Solving systems of nonlinear equations is one of the most difficult problems in all of numerical computation and in a diverse range of engineering applications. Newton’s method for solving systems of nonlinear equations can be highly sensitive to the initial guess of the solution. In this study, a new particle swarm optimization algorithm is proposed to solve systems of nonlinear equations. Some standard systems are presented to demonstrate the efficiency of this method.     

Zakai equation of nonlinear filtering with unbounded coefficients. The case of dependent noises

In this paper, we prove that the unnormalized filter associated with nonlinear filtering problems with dependent noises and a one-dimensional observations process the coefficients of which are unbounded solves a parabolic stochastic partial differential equation, the Zakai equation. The robust form of the Zakai equation is also computed.