A meshless method for numerical solution of a linear hyperbolic equation with variable coefficients in two space dimensions

A meshless method is proposed for the numerical solution of the two space dimensional linear hyperbolic equation subject to appropriate initial and Dirichlet boundary conditions. The new developed scheme uses collocation points and approximates the solution employing thin plate splines radial basis functions. Numerical results are obtained for various cases involving variable, singular and constant coefficients, and are compared with analytical solutions to confirm the good accuracy of the presented scheme. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

A numerical method for solving the hyperbolic telegraph equation

Recently, it is found that telegraph equation is more suitable than ordinary diffusion equation in modelling reaction diffusion for such branches of sciences. In this article, we propose a numerical scheme to solve the one‐dimensional hyperbolic telegraph equation using collocation points and approximating the solution using thin plate splines radial basis function. The scheme works in a similar fashion as finite difference methods. The results of numerical experiments are presented, and are compared with analytical solutions to confirm the good accuracy of the presented scheme. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

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     

Operator splitting for numerical solutions of the RLW Equation

In this study, the numerical behavior of the one-dimensional Regularized Long Wave (RLW) equation has been sought by the Strang splitting technique with respect to time. For this purpose, cubic B-spline functions are used with the finite element collocation method. Then, single solitary wave motion, the interaction of two solitary waves and undular bore problems have been studied and the effectiveness of the method has been investigated. The new results have been compared with those of some of the previous studies available in the literature. The stability analysis has also been taken into account by the von Neumann method.     

A numerical method for one‐dimensional nonlinear Sine‐Gordon equation using collocation and radial basis functions

In this article, we propose a numerical scheme to solve the one‐dimensional undamped Sine‐Gordon equation using collocation points and approximating the solution using Thin Plate Splines (TPS) radial basis function (RBF). The scheme works in a similar fashion as finite difference methods. The results of numerical experiments are presented and are compared with analytical solutions to confirm the good accuracy of the presented scheme.© 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

A meshfree method for the numerical solution of the RLW equation

In this paper, we present a meshfree technique for the numerical solution of the regularized long wave (RLW) equation. This approach is based on a global collocation method using the radial basis functions (RBFs). Different kinds of RBFs are used for this purpose. Accuracy of the new method is tested in terms of _L2$L_2$ and _L∞$L_{\infty }$ error norms. In case of non-availability of the exact solution, performance of the new method is compared with existing methods. Stability analysis of the method is established. Propagation of single and double solitary waves, wave undulation, and conservation properties of mass, energy and momentum of the RLW equation are discussed.     

The Evolution of Internal Undular Bores over a Slope in the Presence of Rotation

The large‐amplitude internal waves commonly observed in the coastal ocean often take the form of unsteady undular bores. Hence, here, we examine the long‐time combined effect of variable topography and background rotation on the propagation of internal undular bores, using the framework of a variable‐coefficient Ostrovsky equation. Because the leading waves in an internal undular bore are close to solitary waves, we first examine the evolution of a single solitary wave. Then, we consider an internal undular bore, for which two methods of generation are used. One method is the matured undular bore developed from an initial shock box in the Korteweg–de Vries equation, that is the Ostrovsky equation with the rotational term omitted, and the other method is a modulated cnoidal wave solution of the same Korteweg–de Vries equation. It transpires that in the long‐time model simulations, the rotational effect disintegrates the nonlinear waves into inertia‐gravity waves, and then there emerge complicated interactions between these inertia‐gravity waves and the modulated periodic waves of the undular bore, especially at the rear part of the undular bore. However, near the front of the undular bore, nonlinear effects further modulate these waves, with the eventual emergence of nonlinear envelope wave packets.     

Numerical solution of the nonlinear Klein–Gordon equation using radial basis functions

The nonlinear Klein–Gordon equation is used to model many nonlinear phenomena. In this paper, we propose a numerical scheme to solve the one-dimensional nonlinear Klein–Gordon equation with quadratic and cubic nonlinearity. Our scheme uses the collocation points and approximates the solution using Thin Plate Splines (TPS) radial basis functions (RBF). The implementation of the method is simple as finite difference methods. The results of numerical experiments are presented, and are compared with analytical solutions to confirm the good accuracy of the presented scheme.     

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     

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

Both sextic and septic B‐spline collocation algorithms are presented for the numerical solutions of the RLW equation. Numerical results resolve the fine structure of the single solitary wave propagation, two and three solitary waves interaction, and evolution of solitary waves. Comparison of the numerical results is done by the results of some earlier schemes mentioned in the article. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 581–607, 2011     

Solitary waves of the EW and RLW equations

Eight finite difference methods are employed to study the solitary waves of the equal-width (EW) and regularized long–wave (RLW) equations. The methods include second-order accurate (in space) implicit and linearly implicit techniques, a three-point, fourth-order accurate, compact operator algorithm, an exponential method based on the local integration of linear, second-order ordinary differential equations, and first- and second-order accurate temporal discretizations. It is shown that the compact operator method with a Crank–Nicolson discretization is more accurate than the other seven techniques as assessed for the three invariants of the EW and RLW equations and the L₂-norm errors when the exact solution is available. It is also shown that the use of Gaussian initial conditions may result in the formation of either positive or negative secondary solitary waves for the EW equation and the formation of positive solitary waves with or without oscillating tails for the RLW equation depending on the amplitude and width of the Gaussian initial conditions. In either case, it is shown that the creation of the secondary wave may be preceded by a steepening and an narrowing of the initial condition. The creation of a secondary wave is reported to also occur in the dissipative RLW equation, whereas the effects of dissipation in the EW equation are characterized by a decrease in amplitude, an increase of the width and a curving of the trajectory of the solitary wave. The collision and divergence of solitary waves of the EW and RLW equations are also considered in terms of the wave amplitude and the invariants of these equations.     

Solitary waves of the EW and RLW equations

Eight finite difference methods are employed to study the solitary waves of the equal-width (EW) and regularized long–wave (RLW) equations. The methods include second-order accurate (in space) implicit and linearly implicit techniques, a three-point, fourth-order accurate, compact operator algorithm, an exponential method based on the local integration of linear, second-order ordinary differential equations, and first- and second-order accurate temporal discretizations. It is shown that the compact operator method with a Crank–Nicolson discretization is more accurate than the other seven techniques as assessed for the three invariants of the EW and RLW equations and the L₂-norm errors when the exact solution is available. It is also shown that the use of Gaussian initial conditions may result in the formation of either positive or negative secondary solitary waves for the EW equation and the formation of positive solitary waves with or without oscillating tails for the RLW equation depending on the amplitude and width of the Gaussian initial conditions. In either case, it is shown that the creation of the secondary wave may be preceded by a steepening and an narrowing of the initial condition. The creation of a secondary wave is reported to also occur in the dissipative RLW equation, whereas the effects of dissipation in the EW equation are characterized by a decrease in amplitude, an increase of the width and a curving of the trajectory of the solitary wave. The collision and divergence of solitary waves of the EW and RLW equations are also considered in terms of the wave amplitude and the invariants of these equations.     

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     

A Multi-Symplectic Scheme for RLW Equation

The Hamiltonian and multi-symplectic formulations for RLW equation are considered in this paper. A new twelve-point difference scheme which is equivalent to multi-symplectic Preissmann integrator is derived based on the multi-symplectic formulation of RLW equation. And the numerical experiments on solitary waves are also given. Comparing the numerical results for RLW equation with those for KdV equation, the inelastic behavior of RLW equation is shown.     

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.     

Investigation of Coriolis effect on oceanic flows and its bifurcation via geophysical Korteweg–de Vries equation

In this work, we have investigated Coriolis effect on oceanic flows in the equatorial region with the help of geophysical Korteweg–de Vries equation (GKdVE). First, Lie symmetries and conservation laws for the GKdVE have been studied. Later, we implement finite element method for numerical simulations. Propagation of nonlinear solitary structures, their interaction and advancement of solitons can be seen in the results so produced. Additionally, Gaussian initial condition and undular bore initial condition are also investigated. Results so obtained have been found in perfect agreement with the available results. Bifurcation analysis of the oceanic traveling wave of the GKdVE is presented depending on traveling wave velocity and Coriolis parameter. It is discerned that velocity of the traveling wave and Coriolis parameter affect significantly on the propagation of the nonlinear waves.     

Galerkin method for the numerical solution of the RLW equation using quintic B-splines

The regularized long wave equation (RLW) is solved numerically by using the quintic B-spline Galerkin finite element method. The same method is applied to the time-split RLW equation. Comparison is made with both analytical solutions and some previous results. Propagation of solitary waves, interaction of two solitons are studied.     

RLW方程的一个新的显式多辛格式

以多辛Euler-box格式为基础对正则长波(RLW)方程的初边值问题进行了讨论,推导了一个新的显式10点格式.模拟孤立波的数值实验表明,这个新的多辛格式是行之有效的,能很好的反映出RLW方程的非弹性性质.     

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.     

Numerical approximation of the generalized regularized long wave equation using Petrov–Galerkin finite element method

The generalized regularized long wave (GRLW) equation has been developed to model a variety of physical phenomena such as ion‐acoustic and magnetohydrodynamic waves in plasma, nonlinear transverse waves in shallow water and phonon packets in nonlinear crystals. This paper aims to develop and analyze a powerful numerical scheme for the nonlinear GRLW equation by Petrov–Galerkin method in which the element shape functions are cubic and weight functions are quadratic B‐splines. The proposed method is implemented to three reference problems involving propagation of the single solitary wave, interaction of two solitary waves and evolution of solitons with the Maxwellian initial condition. The variational formulation and semi‐discrete Galerkin scheme of the equation are firstly constituted. We estimate rate of convergence of such an approximation. Using Fourier stability analysis of the linearized scheme we show that the scheme is unconditionally stable. To verify practicality and robustness of the new scheme error norms L₂, L_∞ and three invariants I₁, I₂, and I₃ are calculated. The computed numerical results are compared with other published results and confirmed to be precise and effective.