Application of the GRP scheme to open channel flow equations

Abstract

The GRP (generalized Riemann problem) scheme, originally conceived for gasdynamics, is reformulated for the numerical integration of the shallow water equations in channels of rectangular cross-section, variable width and bed profile, including a friction model for the fluid-channel shear stress. This scheme is a second-order analytic extension of the first-order Godunov-scheme, based on time-derivatives of flow variables at cell-interfaces resulting from piecewise-linear data reconstruction in cells. The second-order time-integration is based on solutions to generalized Riemann problems at cell-interfaces, thus accounting for the full governing equations, including source terms. The source term due to variable bed elevation is treated in a well-balanced way so that quiescent flow is exactly replicated; this is done by adopting the Surface Gradient Method (SGM). Several problems of steady or unsteady open channel flow are considered, including the terms corresponding to variable channel width and bed elevation, as well as to shear stress at the fluid-channel interface (using the Manning friction model). In all these examples remarkable agreement is obtained between the numerical integration and the exact or accurate solutions.