A central-upwind scheme with artificial viscosity for shallow-water flows in channels

• Published in: Advances in water resources 2016-10, Vol.96, p.323-338
• Main Authors: Hernandez-Duenas, Gerardo, Beljadid, Abdelaziz
• Format: Article
• Language: English
• Subjects: Channels; Channels with irregular geometry; Computational efficiency; Hyperbolic systems of conservation and balance laws; Mathematical models; Preserves; Reconstruction; Semi-discrete central-upwind schemes; Shallow-Water equations; Topography; Viscosity; Walls
• ISSN: 0309-1708; 1872-9657
• DOI: 10.1016/j.advwatres.2016.07.021

•A scheme for shallow-water flows in channels with artificial viscosity is proposed.•The central-upwind scheme enjoys well-balance and positivity-preserving properties.•The non-oscillatory scheme uses high-order reconstructions of cell interfaces.•The computational cost is low compared to the use of nonlinear limiters.•A good agreement of numerical results and experimental data is observed. We develop a new high-resolution, non-oscillatory semi-discrete central-upwind scheme with artificial viscosity for shallow-water flows in channels with arbitrary geometry and variable topography. The artificial viscosity, proposed as an alternative to nonlinear limiters, allows us to use high-resolution reconstructions at a low computational cost. The scheme recognizes steady states at rest when a delicate balance between the source terms and flux gradients occurs. This balance in irregular geometries is more complex than that taking place in channels with vertical walls. A suitable technique is applied by properly taking into account the effects induced by the geometry. Incorporating the contributions of the artificial viscosity and an appropriate time step restriction, the scheme preserves the positivity of the water’s depth. A description of the proposed scheme, its main properties as well as the proofs of well-balance and the positivity of the scheme are provided. Our numerical experiments confirm stability, well-balance, positivity-preserving properties and high resolution of the proposed method. Comparisons of numerical solutions obtained with the proposed scheme and experimental data are conducted, showing a good agreement. This scheme can be applied to shallow-water flows in channels with complex geometry and variable bed topography.