Open Water Flow in a Wet/Dry Multiply-Connected Channel Network: A Robust Numerical Modeling Algorithm

Our goal was to develop a robust algorithm for numerical simulation of one-dimensional shallow water flow in a complex multiply-connected channel network with arbitrary geometry and variable topography. We apply a central-upwind scheme with a novel reconstruction of the open water surface in partially flooded cells that does not require additional correction. The proposed reconstruction and an exact integration of source terms for the momentum conservation equation provide positivity preserving and well-balanced features of the scheme for various wet/dry states. We use two models based on the continuity equation and mass and momentum conservation equations integrated for a control volume around the channel junction to its treatment. These junction models permit to simulate subcritical and supercritical flows in a channel network. Numerous numerical experiments demonstrate the robustness of the proposed numerical algorithm and a good agreement of numerical results with exact solutions, experimental data, and results of the previous numerical studies. The proposed new specialized test on inundation and drying of an initially dry channel network shows the merits of the new numerical algorithm to simulate the subcritical/supercritical open water flows in the networks.

[1]  E. Tadmor,et al.  New High-Resolution Central Schemes for Nonlinear Conservation Laws and Convection—Diffusion Equations , 2000 .

[2]  N. Gouta,et al.  A finite volume solver for 1D shallow‐water equations applied to an actual river , 2002 .

[3]  Chi-Wang Shu,et al.  High-order finite volume WENO schemes for the shallow water equations with dry states , 2011 .

[4]  Christian Igel,et al.  Efficient covariance matrix update for variable metric evolution strategies , 2009, Machine Learning.

[5]  Gary W. Brunner,et al.  HEC-RAS River Analysis System. Hydraulic Reference Manual. Version 1.0. , 1995 .

[6]  B. Hodges,et al.  Conservative finite-volume forms of the Saint-Venant equations for hydrology and urban drainage , 2018, Hydrology and Earth System Sciences.

[7]  J. Zhou Shallow Water Flows , 2004 .

[8]  Notes on Galerkin-finite element methods for the Shallow Water equations with characteristic boundary conditions , 2015, 1507.08209.

[9]  Wanai Li,et al.  The Discontinuous Galerkin Method , 2014 .

[10]  G. Petrova,et al.  A SECOND-ORDER WELL-BALANCED POSITIVITY PRESERVING CENTRAL-UPWIND SCHEME FOR THE SAINT-VENANT SYSTEM ∗ , 2007 .

[11]  S. Zalesak Fully multidimensional flux-corrected transport algorithms for fluids , 1979 .

[12]  H. Yoshioka,et al.  Numerical comparison of shallow water models in multiply connected open channel networks , 2015 .

[13]  H. Yoshioka,et al.  A dual finite volume method scheme for catastrophic flash floods in channel networks , 2015 .

[14]  Wolfgang Raskob,et al.  Updated module of radionuclide hydrological dispersion of the Decision Support System RODOS , 2018 .

[15]  Bernardo Cockburn,et al.  Quasimonotone Schemes for Scalar Conservation Laws. Part II , 1990 .

[16]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme. II. Monotonicity and conservation combined in a second-order scheme , 1974 .

[17]  Gerardo Hernández-Dueñas,et al.  A POSITIVITY PRESERVING CENTRAL SCHEME FOR SHALLOW WATER FLOWS IN CHANNELS WITH WET-DRY STATES , 2014 .

[18]  Hydrological dispersion module of JRODOS: renewed chain of the emergency response models of radionuclide dispersion through watersheds and rivers , 2016 .

[19]  Alexander Kurganov,et al.  Well-Balanced Positivity Preserving Central-Upwind Scheme for the Shallow Water System with Friction Terms , 2013 .

[20]  Christian L. Müller,et al.  pCMALib: a parallel fortran 90 library for the evolution strategy with covariance matrix adaptation , 2009, GECCO '09.

[21]  Xiangxiong Zhang,et al.  Maximum-principle-satisfying and positivity-preserving high-order schemes for conservation laws: survey and new developments , 2011, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[22]  A combined computational algorithm for solving the problem of long surface waves runup on the shore , 2016 .

[23]  Raul Borsche,et al.  Junction-Generalized Riemann Problem for stiff hyperbolic balance laws in networks: An implicit solver and ADER schemes , 2016, J. Comput. Phys..

[24]  Cheng-Wei Yu Consistent Initial Conditions for the Saint-Venant Equations in River Network Modeling , 2017 .

[25]  Roger Temam,et al.  Boundary value problems for the shallow water equations with topography , 2011 .

[26]  ContarinoChristian,et al.  Junction-Generalized Riemann Problem for stiff hyperbolic balance laws in networks , 2016 .

[27]  J. Cunge,et al.  Practical aspects of computational river hydraulics , 1980 .

[28]  Charles S. Melching,et al.  Full Equations (FEQ) model for the solution of the full, dynamic equations of motion for one-dimensional unsteady flow in open channels and through control structures , 1997 .

[29]  Smadar Karni,et al.  A CENTRAL SCHEME FOR SHALLOW WATER FLOWS ALONG CHANNELS WITH IRREGULAR GEOMETRY , 2009 .

[30]  MODEL OF DAM-BREAK FLOODS FOR CHANNELS OF ARBITRARY CROSS SECTION , 1993 .

[31]  I. Bohachevsky,et al.  Finite difference method for numerical computation of discontinuous solutions of the equations of fluid dynamics , 1959 .

[32]  E. Tadmor,et al.  Non-oscillatory central differencing for hyperbolic conservation laws , 1990 .

[33]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme V. A second-order sequel to Godunov's method , 1979 .

[34]  A. V. Rodionov,et al.  Monotonic scheme of the second order of approximation for the continuous calculation of non-equilibrium flows , 1987 .

[35]  G. Hernández-Dueñas,et al.  A central-upwind scheme with artificial viscosity for shallow-water flows in channels , 2016 .

[36]  Smadar Karni,et al.  Shallow Water Flows in Channels , 2011, J. Sci. Comput..

[37]  P. Sweby High Resolution Schemes Using Flux Limiters for Hyperbolic Conservation Laws , 1984 .

[38]  Brett F. Sanders M. Iahr High-resolution and non-oscillatory solution of the St. Venant equations in non-rectangular and non-prismatic channels , 2001 .

[39]  Y. Onishi,et al.  Sediment and (137)Cs behaviors in the Ogaki Dam Reservoir during a heavy rainfall event. , 2014, Journal of environmental radioactivity.

[40]  M J Zheleznyak,et al.  Mathematical modeling of radionuclide dispersion in the Pripyat-Dnieper aquatic system after the Chernobyl accident. , 1992, The Science of the total environment.

[41]  S. Zalesak Introduction to “Flux-Corrected Transport. I. SHASTA, A Fluid Transport Algorithm That Works” , 1997 .

[42]  Bo Wang,et al.  Similarity solution of dam-break flow on horizontal frictionless channel , 2011 .

[43]  Sebastian Noelle,et al.  Finite Volume Evolution Galerkin Methods for the Shallow Water Equations with Dry Beds , 2011, 1501.03628.

[44]  Guofu Huang,et al.  Theoretical Solution of Dam-Break Shock Wave , 1999 .

[45]  S. Grant,et al.  Case Study: Modeling Tidal Transport of Urban Runoff in Channels Using the Finite-Volume Method , 2001 .

[46]  A. Chertock,et al.  Central-upwind scheme for shallow water equations with discontinuous bottom topography , 2016 .

[47]  Abdul A. Khan,et al.  Discontinuous Galerkin Method for 1D Shallow Water Flow in Nonrectangular and Nonprismatic Channels , 2012 .

[48]  Finding nonoscillatory solutions to difference schemes for the advection equation , 2008 .

[49]  Francesca Bellamoli,et al.  A numerical method for junctions in networks of shallow-water channels , 2017, Appl. Math. Comput..

[50]  Yulong Xing High order finite volume WENO schemes for the shallow water flows through channels with irregular geometry , 2016, J. Comput. Appl. Math..

[51]  Hsieh,et al.  NUMERICAL SIMULATIONS OF SCOUR AND DEPOSITION IN A CHANNEL NETWORK , 2003 .

[52]  Doron Levy,et al.  CENTRAL-UPWIND SCHEMES FOR THE SAINT-VENANT SYSTEM , 2002 .

[53]  A. Khan,et al.  Numerical solution of the Saint-Venant equations by an efficient hybrid finite-volume/finite-difference method , 2018 .