Shallow water equations (SWEs) are depth-integrated continuity and momentum equations for river flow, in which it is assumed that surface waves are long, the depth is relatively small in comparison with the horizontal scale of the flow domain, and the pressure is hydrostatic. In the version considered in the present paper, viscous terms are neglected, which restricts the applicability of the equations but is acceptable for a wide range of flows in open channel hydraulics (see, e.g., Abbott and Minns [
1] and Toro [
2]). The SWEs can be formulated as a nonlinear hyperbolic system that can be discretized using an upwind method such as the Roe [
3], Osher [
4], and AUSM [
5] schemes. Alternative higher-order schemes such as WENO [
6,
7] also work by reconstruction, using additional cells next to the interface. However, all the foregoing schemes are constrained by a time step stability condition, such as the Courant–Frederichs–Lewy number CFL < 1, where CFL =
umaxΔ
t/Δ
x, in which
umax is the maximum flow speed, Δ
t is the time step, and Δ
x is a spatial increment in the computational mesh. In the 1980s, LeVeque [
8,
9] proposed the large time step (LTS) method, which led to considerable gains in computational efficiency. In practice, LTS does not conflict with CFL theory in that the interface flux is computed across multiple cells in LTS rather than a single cell in traditional schemes, thus enabling the time step to be increased accordingly. In LTS, more cells take part directly in the computation of fluxes, unlike higher-order schemes that use reconstruction schemes to handle cells adjacent to the interface. Murillo et al. [
10] and Morales-Hernández [
11,
12,
13] then applied LTS to solve the SWEs for free surface flows, and similarly, Qian [
14] used LTS in solving the Euler equations in aerodynamics. All the foregoing LTS schemes are first order. However, as the CFL number becomes large, spurious oscillations appear in the solution, which do not occur in traditional first-order schemes. Harten introduced the total variation diminishing (TVD) concept [
15] that could suppress oscillations for second- or higher-order schemes in the traditional form. Harten [
16] then proposed a TVD-LTS scheme that was further developed by Qian and Lee [
17].
Godunov-type solvers of the SWEs can generate non-physical flow phenomena from discretized bed slope source terms that are out of balance with the flux gradient terms. Well-balanced schemes overcome such difficulties. LeVeque [
18] proposed a quasi-steady wave propagation method that incorporated source terms within high-resolution Godunov methods and obtained accurate results for the propagation of small perturbations. Hubbard and García-Navarro [
19] utilized a flux difference splitting method to ensure the flux gradient and source terms were correctly balanced. Zhou et al. [
20] suggested a surface gradient method to represent the bed slope source term. Rogers et al. [
21] presented an algebraic technique for balancing the flux gradient and source terms for single waterbodies. Liang et al. [
22] extended the algebraic technique to multiple waterbodies, with stage and discharge selected as dependent variables.
This paper applies TVD-LTS to solving the SWEs without viscous stresses and expands the scheme to non-homogeneous hyperbolic systems for the first time, to the authors’ knowledge. The results indicate that TVD-LTS gives a satisfactory representation of benchmark flows while suppressing spurious oscillations when solving homogeneous SWEs but is not as successful when applied to non-homogeneous SWEs. The paper has the following structure.
Section 2 outlines the TVD-LTS solver of the shallow water equations.
Section 3 presents results for different benchmark tests: homogeneous cases of two rarefactions; two shock waves; a single rarefaction and a shock wave; non-homogeneous cases of two rarefactions over a step; two shockwaves over a step; a single rarefaction and a bore over a step; and a transcritical flow over a bed hump.
Section 4 provides the concluding remarks.