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     

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.     

Solitary wave solutions of the MRLW equation using radial basis functions

In this study, traveling wave solutions of the modified regularized long wave (MRLW) equation are simulated by using the meshless method based on collocation with well‐known radial basis functions. The method is tested for three test problems which are single solitary wave motion, interaction of two solitary waves and interaction of three solitary waves. Invariant values for all test problems are calculated, also L₂, L_∞ norms and values of the absolute error for single solitary wave motion are calculated. Numerical results by using the meshless method with different radial basis functions are presented. Figures of wave motions for all test problems are shown. Altogether, meshless methods with radial basis functions solve the MRLW equation very satisfactorily.© 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 28: 235–247, 2012     

Study of magnetosonic solitary waves for the electron magnetohydrodynamics equations

Electron magnetohydrodynamics equations are derived with allowance for nonlinearity, dispersion, and dissipation caused by friction between the ions and electrons. These equations are transformed into a form convenient for the construction of a numerical scheme. The interaction of codirectional and oppositely directed magnetosonic solitary waves with no dissipation is computed. In the first case, the solitary waves are found to behave as solitons (i.e., their amplitudes after the interaction remain the same), while, in the second case, waves are emitted that lead to decreased amplitudes. The decay of a solitary wave due to dissipation is computed. In the case of weak dissipation, the solution is similar to that of the Riemann problem with a structure combining a discontinuity and a solitary wave. The decay of a solitary wave due to dispersion is also computed, in which case the solution can also be interpreted as one with a discontinuity. The decay of a solitary wave caused by the combined effect of dissipation and dispersion is analyzed.     

Nonlinear Waves in a Rod: Results for Incompressible Elastic Materials

A nonlinear intrinsic theory is used to describe the motions of a straight round elastic rod including the influence of radial shear and inertia. Consideration of steady wave motions reduces the two coupled partial differential equations to ordinary differential equations for which two integrals of the motion may be found. For incompressible elastic materials with the restriction of small strain gradients, but arbitrary finite strains, a large variety of exact solutions may be found by quadrature. These include large amplitude periodic waves (which may contain shocks), solitary waves, and in some cases waves that are transitional from one stress level to another. Such solutions may be found for uniform stress strain curves that are concave up or down or that contain inflections, and even for nonmontonic curves, which have been used to represent phase transitions.     

Quartic B‐spline collocation algorithms for numerical solution of the RLW equation

The collocation method based on quartic B‐spline interpolation is studied for numerical solution of the regularized long wave (RLW) equation. The time‐split RLW equation is also solved with the quartic B‐spline collocation method. Numerical accuracy is tested by obtaining the single solitary wave solution. Then, interaction, undulation and evolution of solitary waves are studied. Solutions are compared with available results. Conservation quantities are computed for all test experiments. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2007     

Solitary waves and the structures of discontinuities in non-dissipative models with complex dispersion

Stationary solutions of reversible evolutionary equations of mechanics with higher derivatives are analysed. A two-dimensional graphical method for investigating the solutions of systems of ordinary differential equations is described, which enables one to find special types of solutions: periodic waves, solitary waves and the structures of discontinuities. At the same time, solitary waves can be obtained by taking the limit of sequences of periodic waves and the structures of discontinuities obtained by taking the limit of sequences of solitary waves. This general approach has enabled the existence of all earlier predicted structures to be verified has enabled new types of structures (three-wave structures) to be revealed and has enabled all the necessary conditions at the discontinuities to be found. All the previously known types of solitary waves are found and new types of solitary waves are revealed (generalized ordinary and 1:1 multisolitons). Methods of finding generalized solitary waves, including those with a finite amplitude of the periodic component, are determined. Examples of the solution of the following problems are given for a fourth-order system: generalized solitary waves as the limiting solutions of two-wave resonance solutions, generalized solitary waves and the structure of a discontinuity with three waves, a 1:1 soliton and the structure of a discontinuity with a single radiated wave, a solitary wave with fixed propagation velocity, and the structure of a discontinuity in the form of a kink with radiation. A generalized 1:1 soliton and the structure of a discontinuity with two radiated waves is considered in the case of sixth-order systems. The discussion is mainly based on the example of travelling waves described by the generalized Korteweg-de Vries equations. Other models with complex dispersion (a plasma and a stratified fluid) are also considered.     

Convective Linear Stability of Solitary Waves for Boussinesq Equations

Boussinesq was the first to explain the existence of Scott Russell's solitary wave mathematically. He employed a variety of asymptotically equivalent equations to describe water waves in the small-amplitude, long-wave regime. We study the linearized stability of solitary waves for three linearly well-posed Boussinesq models. These are problems for which well-developed Lyapunov methods of stability analysis appear to fail. However, we are able to analyze the eigenvalue problem for small-amplitude solitary waves, by comparison to the equation that Boussinesq himself used to describe the solitary wave, which is now called the Korteweg–de Vries equation. With respect to a weighted norm designed to diminish as perturbations convect away from the wave profile, we prove that nonzero eigenvalues are absent in a half-plane of the form R λ>− b for some b >0, for all three Boussinesq models. This result is used to prove the decay of solutions of the evolution equations linearized about the solitary wave, in two of the models. This "convective linear stability" property has played a central role in the proof of nonlinear asymptotic stability of solitary-wave-like solutions in other systems.     

A Model Equation for Wavepacket Solitary Waves Arising from Capillary-Gravity Flows

A model equation governing the primitive dynamics of wave packets near an extremum of the linear dispersion relation at finite wavenumber is derived. In two spatial dimensions, we include the effects of weak variation of the wave in the direction transverse to the direction of propagation. The resulting equation is contrasted with the Kadomtsev–Petviashvilli and Nonlinear Schrödinger (NLS) equations. The model is derived as an approximation to the equations for deep water gravity-capillary waves, but has wider applications. Both line solitary waves and solitary waves which decay in both the transverse and propagating directions—lump solitary waves—are computed. The stability of these waves is investigated and their dynamics are studied via numerical time evolution of the equation.     

Classification of Solitary Wave Bifurcations in Generalized Nonlinear Schrödinger Equations

Bifurcations of solitary waves are classified for the generalized nonlinear Schrödinger equations with arbitrary nonlinearities and external potentials in arbitrary spatial dimensions. Analytical conditions are derived for three major types of solitary wave bifurcations, namely, saddle‐node, pitchfork, and transcritical bifurcations. Shapes of power diagrams near these bifurcations are also obtained. It is shown that for pitchfork and transcritical bifurcations, their power diagrams look differently from their familiar solution‐bifurcation diagrams. Numerical examples for these three types of bifurcations are given as well. Of these numerical examples, one shows a transcritical bifurcation, which is the first report of transcritical bifurcations in the generalized nonlinear Schrödinger equations. Another shows a power loop phenomenon which contains several saddle‐node bifurcations, and a third example shows double pitchfork bifurcations. These numerical examples are in good agreement with the analytical results.     

Accelerated Imaginary-time Evolution Methods for the Computation of Solitary Waves

Two accelerated imaginary-time evolution methods are proposed for the computation of solitary waves in arbitrary spatial dimensions. For the first method (with traditional power normalization), the convergence conditions as well as conditions for optimal accelerations are derived. In addition, it is shown that for nodeless solitary waves, this method converges if and only if the solitary wave is linearly stable. The second method is similar to the first method except that it uses a novel amplitude normalization. The performance of these methods is illustrated on various examples. It is found that while the first method is competitive with the Petviashvili method, the second method delivers much better performance than the first method and the Petviashvili method.     

An iterative method with error estimators

Iterative methods for the solution of linear systems of equations produce a sequence of approximate solutions. In many applications it is desirable to be able to compute estimates of the norm of the error in the approximate solutions generated and terminate the iterations when the estimates are sufficiently small. This paper presents a new iterative method based on the Lanczos process for the solution of linear systems of equations with a symmetric matrix. The method is designed to allow the computation of estimates of the Euclidean norm of the error in the computed approximate solutions. These estimates are determined by evaluating certain Gauss, anti-Gauss, or Gauss–Radau quadrature rules.     

The error norm of Gaussian quadrature formulae for weight functions of Bernstein-Szegö type

Summary We consider the Gaussian quadrature formulae for the Bernstein-Szegö weight functions consisting of any one of the four Chebyshev weights divided by an arbitrary quadratic polynomial that remains positive on [–1, 1]. Using the method in Akrivis (1985), we compute the norm of the error functional of these quadrature formulae. The quality of the bounds for the error functional, that can be obtained in this way, is demonstrated by two numerical examples.     

Do peaked solitary water waves indeed exist?

Many models of shallow water waves, such as the famous Camassa–Holm equation, admit peaked solitary waves. However, it is an open question whether or not the widely accepted peaked solitary waves can be derived from the fully nonlinear wave equations. In this paper, a unified wave model (UWM) based on the symmetry and the fully nonlinear wave equations is put forward for progressive waves with permanent form in finite water depth. Different from traditional wave models, the flows described by the UWM are not necessarily irrotational at crest, so that it is more general. The unified wave model admits not only the traditional progressive waves with smooth crest, but also a new kind of solitary waves with peaked crest that include the famous peaked solitary waves given by the Camassa–Holm equation. Besides, it is proved that Kelvin’s theorem still holds everywhere for the newly found peaked solitary waves. Thus, the UWM unifies, for the first time, both of the traditional smooth waves and the peaked solitary waves. In other words, the peaked solitary waves are consistent with the traditional smooth ones. So, in the frame of inviscid fluid, the peaked solitary waves are as acceptable and reasonable as the traditional smooth ones. It is found that the peaked solitary waves have some unusual and unique characteristics. First of all, they have a peaked crest with a discontinuous vertical velocity at crest. Especially, unlike the traditional smooth waves that are dispersive with wave height, the phase speed of the peaked solitary waves has nothing to do with wave height, but depends (for a fixed wave height) on its decay length, i.e., the actual wavelength: in fact, the peaked solitary waves are dispersive with the actual wavelength when wave height is fixed. In addition, unlike traditional smooth waves whose kinetic energy decays exponentially from free surface to bottom, the kinetic energy of the peaked solitary waves either increases or almost keeps the same. All of these unusual properties show the novelty of the peaked solitary waves, although it is still an open question whether or not they are reasonable in physics if the viscosity of fluid and surface tension are considered.     

Universally-Convergent Squared-Operator Iteration Methods for Solitary Waves in General Nonlinear Wave Equations

Three new iteration methods, namely the squared-operator method, the modified squared-operator method, and the power-conserving squared-operator method, for solitary waves in general scalar and vector nonlinear wave equations are proposed. These methods are based on iterating new differential equations whose linearization operators are squares of those for the original equations, together with acceleration techniques. The first two methods keep the propagation constants fixed, while the third method keeps the powers (or other arbitrary functionals) of the solution fixed. It is proved that all these methods are guaranteed to converge to any solitary wave (either ground state or not) as long as the initial condition is sufficiently close to the corresponding exact solution, and the time step in the iteration schemes is below a certain threshold value. Furthermore, these schemes are fast-converging, highly accurate, and easy to implement. If the solitary wave exists only at isolated propagation constant values, the corresponding squared-operator methods are developed as well. These methods are applied to various solitary wave problems of physical interest, such as higher-gap vortex solitons in the two-dimensional nonlinear Schrödinger equations with periodic potentials, and isolated solitons in Ginzburg–Landau equations, and some new types of solitary wave solutions are obtained. It is also demonstrated that the modified squared-operator method delivers the best performance among the methods proposed in this article.     

Solitary waves for Korteweg-de Vries equation with small delay

This paper deals with the existence of solitary waves for Korteweg-de Vries equation with time delay. Based upon the inertial manifold theory and differential manifold geometric theory, the existence of solitary wave solution is proved when the delay is small enough. Up to now, studies on solitary wave for such delay differential equation are not available, so the results of this paper are new.     

Alternating Direction Finite Volume Element Methods for Three-Dimensional Parabolic Equations

<正>This paper presents alternating direction finite volume element methods for three-dimensional parabolic partial differential equations and gives four computational schemes,one is analogous to Douglas finite difference scheme with second-order splitting error,the other two schemes have third-order splitting error,and the last one is an extended LOD scheme.The L~2 norm and H~1 semi-norm error estimates are obtained for the first scheme and second one,respectively.Finally,two numerical examples are provided to illustrate the efficiency and accuracy of the methods.     

Quadrature error of the load vector in the finite element method

We analyze the error in thep version of the of the finite element method when the effect of the quadrature error is taken in the load vector. We briefly study some results on theH ¹ norm error and present some new results for the error in theL ₂ norm. We investigate the quadrature error due to the numerical integration of the right hand side We present theoretical and computational examples showing the sharpness of our results.     

Structural Stability of Discontinuous Galerkin Schemes

The goal of this work is to determine classes of traveling solitary wave solutions for a differential approximation of a discontinuous Galerkin finite difference scheme by means of an hyperbolic ansatz. It is shown that spurious solitary waves can occur in finite-difference solutions of nonlinear wave equation. The occurence of such a spurious solitary wave, which exhibits a very long life time, results in a non-vanishing numerical error for arbitrary time in unbounded numerical domain. Such a behavior is referred here to have a structural instability of the scheme, since the space of solutions spanned by the numerical scheme encompasses types of solutions (solitary waves in the present case) that are not solutions of the original continuous equations. This paper extends our previous work about classical schemes to discontinuous Galerkin schemes (David and Sagaut in Chaos Solitons Fractals 41(4):2193?C2199, 2009; Chaos Solitons Fractals 41(2):655?C660, 2009).     

Regular and embedded solitons in a generalized Pochammer PDE

Variational methods are employed to generate families of both regular and embedded solitary wave solutions for a generalized Pochammer PDE that is currently of great interest. The technique for obtaining the embedded solitons incorporates several recent generalizations of the usual variational technique and is thus topical in itself. One unusual feature of the solitary waves derived here is that we are able to obtain them in analytical form (within the family of the trial functions). Thus, a direct error analysis is performed, showing the accuracy of the resulting solitary waves. Given the importance of solitary wave solutions in wave dynamics and information propagation in nonlinear PDEs, as well as the fact that only the parameter regimes for the existence of solitary waves had previously been analyzed for the microstructure PDE considered here, the results obtained here are both new and timely.