Non‐oscillatory relaxation methods for the shallow‐water equations in one and two space dimensions

In this paper, a new family of high‐order relaxation methods is constructed. These methods combine general higher‐order reconstruction for spatial discretization and higher order implicit‐explicit schemes or TVD Runge–Kutta schemes for time integration of relaxing systems. The new methods retain all the attractive features of classical relaxation schemes such as neither Riemann solvers nor characteristic decomposition are needed. Numerical experiments with the shallow‐water equations in both one and two space dimensions on flat and non‐flat topography demonstrate the high resolution and the ability of our relaxation schemes to better resolve the solution in the presence of shocks and dry areas without using either Riemann solvers or front tracking techniques. Copyright © 2004 John Wiley & Sons, Ltd.