Abstract
The Riemann wave system has a fundamental role in describing waves in various nonlinear natural phenomena, for instance, tsunamis in the oceans. This paper focuses on executing the generalized exponential rational function approach and some numerical methods to obtain a distinct range of traveling wave structures and numerical results of the two-dimensional Riemann problems. The stability of obtained traveling wave solutions is analyzed by satisfying the constraint conditions of the Hamiltonian system. Numerical simulations are investigated via the finite difference method to verify the accuracy of the obtained results. To extract the approximation solutions to the underlying problem, some ODE solvers in FORTRAN software are applied, and outcomes are shown graphically. The stability and accuracy of the numerical schemes using Fourier's stabilitymethod and error analysis, respectively, to increase the reassurance are investigated. A comparison between the analytical and numerical results is obtained and graphically provided. The proposed methods are effective and practical to be applied for solving more partial differential equations (PDEs).