Solitary Wave Solutions of the Generalized Rosenau-KdV-RLW Equation

This paper investigates the solitary wave solutions of the generalized Rosenau–Korteweg-de Vries-regularized-long wave equation. This model is obtained by coupling the Rosenau–Korteweg-de Vries and Rosenau-regularized-long wave equations. The solution of the equation is approximated by a local meshless technique called radial basis function (RBF) and the finite-difference (FD) method. The association of the two techniques leads to a meshless algorithm that does not requires the linearization of the nonlinear terms. First, the partial differential equation is transformed into a system of ordinary differential equations (ODEs) using radial kernels. Then, the ODE system is solved by means of an ODE solver of higher-order. It is shown that the proposed method is stable. In order to illustrate the validity and the efficiency of the technique, five problems are tested and the results compared with those provided by other schemes.info:eu-repo/semantics/publishedVersio

A Space-time Smooth Artificial Viscosity Method For Nonlinear Conservation Laws

We introduce a new methodology for adding localized, space-time smooth, artificial viscosity to nonlinear systems of conservation laws which propagate shock waves, rarefactions, and contact discontinuities, which we call the $C$-method. We shall focus our attention on the compressible Euler equations in one space dimension. The novel feature of our approach involves the coupling of a linear scalar reaction-diffusion equation to our system of conservation laws, whose solution $C(x,t)$ is the coefficient to an additional (and artificial) term added to the flux, which determines the location, localization, and strength of the artificial viscosity. Near shock discontinuities, $C(x,t)$ is large and localized, and transitions smoothly in space-time to zero away from discontinuities. Our approach is a provably convergent, spacetime-regularized variant of the original idea of Richtmeyer and Von Neumann, and is provided at the level of the PDE, thus allowing a host of numerical discretization schemes to be employed. We demonstrate the effectiveness of the $C$-method with three different numerical implementations and apply these to a collection of classical problems: the Sod shock-tube, the Osher-Shu shock-tube, the Woodward-Colella blast wave and the Leblanc shock-tube. First, we use a classical continuous finite-element implementation using second-order discretization in both space and time, FEM-C. Second, we use a simplified WENO scheme within our $C$-method framework, WENO-C. Third, we use WENO with the Lax-Friedrichs flux together with the $C$-equation, and call this WENO-LF-C. All three schemes yield higher-order discretization strategies, which provide sharp shock resolution with minimal overshoot and noise, and compare well with higher-order WENO schemes that employ approximate Riemann solvers, outperforming them for the difficult Leblanc shock tube experiment.Comment: 34 pages, 27 figure

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 magnetohydro dynamic waves in plasma,nonlinear transverse waves in shallow water and phonon packets in nonlinear crystals. This paper aims to develop andanalyze a powerful numerical scheme for the nonlinear GRLWequation by Petrov–Galerkin method in which the elementshape functions are cubic and weight functions are quadratic B-splines. The proposed method is implemented to three ref-erence problems involving propagation of the single solitarywave, interaction of two solitary waves and evolution of solitons with the Maxwellian initial condition. The variational for-mulation 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 thelinearized scheme we show that the scheme is uncondition-ally stable. To verify practicality and robustness of the new scheme error norms L2, L∞ and three invariants I1, I2,and I3 are calculated. The computed numerical results are compared with other published results and confirmed to be precise and effective

Generic effective source for scalar self-force calculations

A leading approach to the modelling of extreme mass ratio inspirals involves the treatment of the smaller mass as a point particle and the computation of a regularized self-force acting on that particle. In turn, this computation requires knowledge of the regularized retarded field generated by the particle. A direct calculation of this regularized field may be achieved by replacing the point particle with an effective source and solving directly a wave equation for the regularized field. This has the advantage that all quantities are finite and require no further regularization. In this work, we present a method for computing an effective source which is finite and continuous everywhere, and which is valid for a scalar point particle in arbitrary geodesic motion in an arbitrary background spacetime. We explain in detail various technical and practical considerations that underlie its use in several numerical self-force calculations. We consider as examples the cases of a particle in a circular orbit about Schwarzschild and Kerr black holes, and also the case of a particle following a generic time-like geodesic about a highly spinning Kerr black hole. We provide numerical C code for computing an effective source for various orbital configurations about Schwarzschild and Kerr black holes.Comment: 24 pages, 7 figures, final published versio