A Petrov-Galerkin finite element scheme for the regularized long wave equation

A Petrov-Galerkin finite element method (FEM) for the regularized long wave (RLW) equation is proposed. Finite elements are used in both the space and the time domains. Dispersion correction and a highly selective dissipation mechanism are introduced through additional streamline upwind terms in the weight functions. An implicit, conditionally stable, one-step predictor-corrector time integration scheme results. The accuracy and stability are investigated by means of local expansion by Taylor series and the resulting equivalent differential equation. An analysis based on a linear Fourier series solution and the Von Neumann's stability criterion is also performed. Based on the order of the analytical approximations and of the domain discretization it is concluded that the scheme is of third order in the nonlinear version and of fourth order in the linear version. Three numerical experiments of wave propagation are presented and their results compared with similar ones in the literature: solitary wave propagation, undular bore propagation, and cnoidal wave propagation. It is concluded that the present scheme possesses superior conservation and accuracy properties.