Bedload transport in shallow water models: Why splitting (may) fail, how hyperbolicity (can) help

Abstract In this paper, we are concerned with sediment transport models consisting of a shallow water system coupled with the so called Exner equation to describe the evolution of the topography. We show that, for some bedload transport models like the well-known Meyer-Peter and Muller model, the system is hyperbolic and, thus, linearly stable, only under some constraint on the velocity. In practical situations, this condition is hopefully fulfilled. Numerical approximations of such system are often based on a splitting method, solving first shallow water equation on a time step and, updating afterwards the topography. It is shown that this strategy can create spurious/unphysical oscillations which are related to the study of hyperbolicity. Using an upper bound of the largest eigenvector may improve the results although the instabilities cannot be always avoided, e.g. in supercritical regions.

[1]  C. Parés Numerical methods for nonconservative hyperbolic systems: a theoretical framework. , 2006 .

[2]  E. Toro Shock-Capturing Methods for Free-Surface Shallow Flows , 2001 .

[3]  T. Morales de Luna,et al.  On a shallow water model for the simulation of turbidity currents , 2009 .

[4]  A. Shields,et al.  Application of similarity principles and turbulence research to bed-load movement , 1936 .

[5]  Magnus Larson,et al.  A general formula for non-cohesive bed load sediment transport , 2005 .

[6]  Eric W. Weisstein,et al.  The CRC concise encyclopedia of mathematics , 1999 .

[7]  David J. Williams,et al.  Relative celerities of mobile bed flows with finite solids concentrations , 1996 .

[8]  R. Luque Erosion And Transport Of Bed-Load Sediment , 1976 .

[9]  HYPERBOLICITY OF GRAD’S EXTENSION OF HYDRODYNAMIC MODELS FOR IONOSPHERIC PLASMA PART I: THE SINGLE SPECIES CASE , 1994 .

[10]  Olivier Delestre Simulation du ruissellement d'eau de pluie sur des surfaces agricoles. (Rain water overland flow on agricultural fields simulation) , 2010 .

[11]  Shinji Egashira,et al.  Coupled and Decoupled Numerical Modeling of Flow and Morphological Evolution in Alluvial Rivers , 2002 .

[12]  Ethan J. Kubatko,et al.  An unstructured grid morphodynamic model with a discontinuous Galerkin method for bed evolution , 2006 .

[13]  J. Ribberink Bed-load transport for steady flows and unsteady oscillatory flows , 1998 .

[14]  Dennis A. Lyn,et al.  Unsteady Sediment‐Transport Modeling , 1987 .

[15]  P. Raviart,et al.  Numerical Approximation of Hyperbolic Systems of Conservation Laws , 1996, Applied Mathematical Sciences.

[16]  M. Esteves,et al.  PSEM_2D: A physically based model of erosion processes at the plot scale , 2005 .

[17]  Tomás Morales de Luna,et al.  Mathematical Modelling and Numerical Analysis an Entropy Satisfying Scheme for Two-layer Shallow Water Equations with Uncoupled Treatment , 2022 .

[18]  P. Nielsen Coastal Bottom Boundary Layers and Sediment Transport , 1992 .

[19]  Junqing Yang Assimilation de données variationnelle pour les problèmes de transport des sédiments en rivière. (Variational data assimilation for the problems of sediment transport in rivers) , 1999 .

[20]  R. Müller,et al.  Formulas for Bed-Load transport , 1948 .

[21]  F. Bouchut Nonlinear Stability of Finite Volume Methods for Hyperbolic Conservation Laws: and Well-Balanced Schemes for Sources , 2005 .

[22]  Carlos Parés,et al.  On the well-balance property of Roe?s method for nonconservative hyperbolic systems , 2004 .

[23]  Carlos Parés,et al.  A Q-SCHEME FOR A CLASS OF SYSTEMS OF COUPLED CONSERVATION LAWS WITH SOURCE TERM. APPLICATION TO A TWO-LAYER 1-D SHALLOW WATER SYSTEM , 2001 .

[24]  I. Toumi A weak formulation of roe's approximate riemann solver , 1992 .

[25]  A. Vacca,et al.  Discussion of "Influence of Turbulence on Bed Load Sediment Transport" , 2005 .

[26]  M. Altinakar,et al.  St. Venant-Exner Equations for Near-Critical and Transcritical Flows , 2002 .

[27]  Enrique D. Fernández-Nieto,et al.  Sediment transport models in Shallow Water equations and numerical approach by high order finite volume methods , 2008 .

[28]  M. Vázquez-Cendón Improved Treatment of Source Terms in Upwind Schemes for the Shallow Water Equations in Channels with Irregular Geometry , 1999 .