Coupled method for the numerical simulation of 1D shallow water and Exner transport equations in channels with variable cross-section

LIFTEC-CSIC, University of Zaragoza, Spain

^* e-mail: sermar@unizar.es

Abstract

This work is focused on the a numerical finite volume scheme for the resulting coupled shallow water-Exner system in 1D applications with arbitrary geometry. The mathematical expression modeling the the hydrodynamic and morphodynamic components of the physical phenomenon are treated to deal with cross-section shape variations and empirical solid discharge estimations. The resulting coupled system of equations can be rewritten as a nonconservative hyperbolic system with three moving waves and one stationary wave to account for the source terms discretization. But, even for the simplest solid transport models as the Grass law, to find a linearized Jacobian matrix of the system can be a challenge if one considers arbitrary shape channels. Moreover, the bottom channel slope variations depends on the erosion-deposition mechanism considered to update the channel cross-section profile. In this paper a numerical finite volume scheme is proposed, based on an augmented Roe solver (first order accurate in time and space) and dealing with solid transport flux variations caused by the channel geometry changes. Channel crosssection variations lead to the appearance of a new solid flux source term which should be discretized properly. Comparison of the numerical results for several analytical and experimental cases demonstrate the effectiveness, exact wellbalanceness and accuracy of the scheme.