In this work, our main purpose is to develop of a sufficiently robust, accurate and efficient numerical scheme for the solution of the regularized long wave (RLW) equation, an important partial differential equation with quadratic nonlinearity, describing a large number of physical phenomena. The crucial idea is based on the discretization of the RLW equation with the aid of a combination of the discontinuous Galerkin method for the space semi-discretization and the backward difference formula for the time discretization.
Furthermore, a suitable linearization preserves a linear algebraic problem at each time level. We present error analysis of the proposed scheme for the case of nonsymmetric discretization of the dispersive term.
The appended numerical experiments confirm theoretical results and investigate the conservative properties of the RLW equation related to mass, momentum and energy. Both procedures illustrate the potency of the scheme consequently.