A HLLC scheme for nonconservative hyperbolic problems. Application to turbidity currents with sediment transport

The goal of this paper is to obtain a well-balanced, stable, fast, and robust HLLC-type approximate Riemann solver for a hyperbolic nonconservative PDE system arising in a turbidity current model. The main difficulties come from the nonconservative nature of the system. A general strategy to derive simple approximate Riemann solvers for nonconservative systems is introduced, which is applied to the turbidity current model to obtain two different HLLC solvers. Some results concerning the non- negativity preserving property of the corresponding numerical methods are presented. The numerical results provided by the two HLLC solvers are compared between them and also with those obtained with a Roe-type method in a number of 1d and 2d test problems. This comparison shows that, while the quality of the numerical solutions is comparable, the computational cost of the HLLC solvers is lower, as only some partial information of the eigenstructure of the matrix system is needed.

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

[2]  Manuel Jesús Castro Díaz,et al.  High Order Extensions of Roe Schemes for Two-Dimensional Nonconservative Hyperbolic Systems , 2009, J. Sci. Comput..

[3]  Van Rijn,et al.  Sediment transport; Part I, Bed load transport , 1984 .

[4]  W. H. Graf,et al.  Flow structure in turbidity currents , 1996 .

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

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

[7]  Enrique D. Fernández-Nieto,et al.  A consistent intermediate wave speed for a well-balanced HLLC solver , 2008 .

[8]  Michael C. Quick,et al.  Sediment transport by waves and currents , 1983 .

[9]  J. F. A. Sleath,et al.  Sediment transport by waves and currents , 1995 .

[10]  E. Toro,et al.  Restoration of the contact surface in the HLL-Riemann solver , 1994 .

[11]  Stéphane Cordier,et al.  Bedload transport in shallow water models: Why splitting (may) fail, how hyperbolicity (can) help , 2011 .

[12]  T. Hou,et al.  Why nonconservative schemes converge to wrong solutions: error analysis , 1994 .

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

[14]  Y. Kubo,et al.  Experimental and numerical study of topographic effects on deposition from two-dimensional, particle-driven density currents , 2004 .

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

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

[17]  Carlos Parés Madroñal,et al.  On some difficulties of the numerical approximation of nonconservative hyperbolic systems , 2009 .

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

[19]  Manuel Jesús Castro Díaz,et al.  A Duality Method for Sediment Transport Based on a Modified Meyer-Peter & Müller Model , 2011, J. Sci. Comput..

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

[21]  Nikolaos D. Katopodes,et al.  Hydrodynamics of Turbid Underflows. I: Formulation and Numerical Analysis , 1999 .

[22]  P. Lax,et al.  On Upstream Differencing and Godunov-Type Schemes for Hyperbolic Conservation Laws , 1983 .

[23]  E. Toro.  Riemann Solvers and Numerical Methods for Fluid Dynamics , 1997 .

[24]  Enrique D. Fernández-Nieto.  Modelling and numerical simulation of submarine sediment shallow flows: transport and avalanches , 2009 .

[25]  Tomás Morales de Luna,et al.  Semi-discrete Entropy Satisfying Approximate Riemann Solvers. The Case of the Suliciu Relaxation Approximation , 2009, J. Sci. Comput..

[26]  Y. Kubo,et al.  Laboratory experiments and numerical simulation of sediment-wave formation by turbidity currents , 2002 .

[27]  Yusuke Fukushima,et al.  Self-accelerating turbidity currents , 1986, Journal of Fluid Mechanics.

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

[29]  Carlos Parés,et al.  Godunov method for nonconservative hyperbolic systems , 2007 .

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

[31]  Scott F. Bradford,et al.  Numerical modeling of hyperpycnal plume , 2005 .

[32]  Manuel J. Castro,et al.  WELL-BALANCED NUMERICAL SCHEMES BASED ON A GENERALIZED HYDROSTATIC RECONSTRUCTION TECHNIQUE , 2007 .

[33]  Natalia Kopteva,et al.  The Formation of River Channels , 2007, SIAM J. Appl. Math..

[34]  Rémi Abgrall,et al.  A comment on the computation of non-conservative products , 2010, J. Comput. Phys..

[35]  Carlos Parés,et al.  Two-dimensional sediment transport models in shallow water equations. A second order finite volume approach on unstructured meshes , 2009 .

[36]  Carlos Parés Madroñal,et al.  On the Convergence and Well-Balanced Property of Path-Conservative Numerical Schemes for Systems of Balance Laws , 2011, J. Sci. Comput..

[37]  José Miguel Mantas,et al.  Two-Dimensional Compact Third-Order Polynomial Reconstructions. Solving Nonconservative Hyperbolic Systems Using GPUs , 2011, J. Sci. Comput..

[38]  Manuel Jesús Castro Díaz,et al.  Why many theories of shock waves are necessary: Convergence error in formally path-consistent schemes , 2008, J. Comput. Phys..