A mass-conservative centered finite volume model for solving two-dimensional two-layer shallow water equations for fluid mud propagation over varying topography and dry areas

In this paper we develop a finite volume model to solve the two-dimensional shallow water equations governing the propagation of two superimposed layers, with the upper water layer carrying a dilute sediment suspension, and the underlaying layer being a high concentration non-Newtonian fluid mud mixture. The model formulation contains non-conservative terms as well as source terms. We propose a scheme able to deal with varying topography and dry areas, providing well-balanced solutions when both water and fluid mud are quiescent. The model is tested against both exact solutions and numerical examples. The results show the ability of the model to deal with wetting and drying of both water and fluid mud layers, providing mass-conservative solutions. Moreover, the model solves discontinuities and steep fronts, computing accurate and oscillation-free solutions.

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

[2]  Dimitri J. Mavriplis,et al.  Revisiting the Least-squares Procedure for Gradient Reconstruction on Unstructured Meshes , 2003 .

[3]  E. Toro,et al.  Restoration of the contact surface in FORCE-type centred schemes II: Non-conservative one- and two-layer two-dimensional shallow water equations , 2012 .

[4]  M. J. Castro,et al.  FORCE schemes on unstructured meshes II: Non-conservative hyperbolic systems , 2010 .

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

[6]  P. García-Navarro,et al.  On numerical treatment of the source terms in the shallow water equations , 2000 .

[7]  Michael Dumbser,et al.  FORCE schemes on unstructured meshes I: Conservative hyperbolic systems , 2009, J. Comput. Phys..

[8]  Ken Been,et al.  Self-weight consolidation of soft soils: an experimental and theoretical study , 1981 .

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

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

[11]  P. García-Navarro,et al.  Godunov-type methods for free-surface shallow flows: A review , 2007 .

[12]  Hervé Capart,et al.  Two-layer shallow water computation of mud flow intrusions into quiescent water , 2007 .

[13]  Valerio Caleffi,et al.  High-order balanced CWENO scheme for movable bed shallow water equations , 2007 .

[14]  Randall J. LeVeque,et al.  Balancing Source Terms and Flux Gradients in High-Resolution Godunov Methods , 1998 .

[15]  E. Bertolazzi,et al.  A CELL-CENTERED SECOND-ORDER ACCURATE FINITE VOLUME METHOD FOR CONVECTION–DIFFUSION PROBLEMS ON UNSTRUCTURED MESHES , 2004 .

[16]  Mutsuto Kawahara,et al.  Finite element method for moving boundary problems in river flow , 1986 .

[17]  Jian-zhong Zhou,et al.  An unstructured finite volume model for dam‐break floods with wet/dry fronts over complex topography , 2011 .

[18]  Chongam Kim,et al.  Multi-dimensional limiting process for hyperbolic conservation laws on unstructured grids , 2005, J. Comput. Phys..

[19]  Manuel Jesús Castro Díaz,et al.  On a well-balanced high-order finite volume scheme for shallow water equations with topography and dry areas , 2007, J. Comput. Phys..

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

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

[22]  Luka Sopta,et al.  Balanced finite volume WENO and central WENO schemes for the shallow water and the open-channel flow equations , 2004 .

[23]  Alfredo Bermúdez,et al.  Upwind methods for hyperbolic conservation laws with source terms , 1994 .

[24]  Yulong Xing,et al.  High order well-balanced finite volume WENO schemes and discontinuous Galerkin methods for a class of hyperbolic systems with source terms , 2006, J. Comput. Phys..

[25]  Michael Dumbser,et al.  Well-Balanced High-Order Centred Schemes for Non-Conservative Hyperbolic Systems. Applications to Shallow Water Equations with Fixed and Mobile Bed , 2009 .

[26]  Qiuhua Liang,et al.  Adaptive quadtree simulation of shallow flows with wet-dry fronts over complex topography , 2009 .

[27]  P. Lax.  Weak solutions of nonlinear hyperbolic equations and their numerical computation , 1954 .

[28]  An efficient and conservative hybrid method for solving multidimensional conservation laws , 2009 .

[29]  A. Cooper,et al.  A two-dimensional model of the movement of fluid mud in a high energy turbid estuary , 1988 .

[30]  C. Friedrichs,et al.  Management of Fluid Mud in Estuaries, Bays, and Lakes. I: Present State of Understanding on Character and Behavior , 2007 .

[31]  Q. Liang,et al.  Numerical resolution of well-balanced shallow water equations with complex source terms , 2009 .

[32]  Manuel Jesús Castro Díaz,et al.  Available Online at Www.sciencedirect.com Mathematical and So,snos ~__d,~ot" Computer Modelling the Numerical Treatment of Wet/dry Fronts in Shallow Flows: Application to One-layer and Two-layer Systems , 2022 .

[33]  B. Flemming,et al.  Fluid mud dynamics in the Weser estuary turbidity zone tracked by high-resolution side-scan sonar and parametric sub-bottom profiler , 2006 .

[34]  G. Strang.  On the Construction and Comparison of Difference Schemes , 1968 .

[35]  David M. Farmer,et al.  Maximal two-layer exchange over a sill and through the combination of a sill and contraction with barotropic flow , 1986, Journal of Fluid Mechanics.

[36]  Carl T. Friedrichs,et al.  Effects of ambient currents and waves on gravity-driven sediment transport on continental shelves , 2001 .

[37]  J. Turner,et al.  Buoyancy Effects in Fluids , 1973 .

[38]  Manuel Jesús Castro Díaz,et al.  High order finite volume schemes based on reconstruction of states for solving hyperbolic systems with nonconservative products. Applications to shallow-water systems , 2006, Math. Comput..

[39]  Pilar García-Navarro,et al.  Flux difference splitting and the balancing of source terms and flux gradients , 2000 .

[40]  J. F. Verweij,et al.  Far-field impact of water injection dredging in the Crouch River , 2002 .

[41]  Pilar García-Navarro,et al.  Zero mass error using unsteady wetting–drying conditions in shallow flows over dry irregular topography , 2004 .

[42]  Manuel Jesús Castro Díaz,et al.  Numerical Treatment of the Loss of Hyperbolicity of the Two-Layer Shallow-Water System , 2011, J. Sci. Comput..

[43]  R. Sternberg,et al.  Fluid-mud processes on the Amazon continental shelf , 1996 .

[44]  C. Friedrichs,et al.  Observations of sediment transport on the continental shelf off the mouth of the Waiapu River, New Zealand: Evidence for current-supported gravity flows , 2008 .

[45]  I. Nikolos,et al.  An unstructured node-centered finite volume scheme for shallow water flows with wet/dry fronts over complex topography , 2009 .

[46]  Brett F. Sanders,et al.  Adaptive Godunov-Based Model for Flood Simulation , 2008 .

[47]  Emmanuel Partheniades,et al.  Erosion and Deposition of Cohesive Soils , 1965 .

[48]  E. Toro,et al.  Well-balanced high-order centered schemes on unstructured meshes for shallow water equations with fixed and mobile bed , 2010 .

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

[50]  Timothy J. Barth,et al.  The design and application of upwind schemes on unstructured meshes , 1989 .

[51]  D. Causon,et al.  The surface gradient method for the treatment of source terms in the shallow-water equations , 2001 .

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

[53]  C. Kranenburg,et al.  Erosion of Fluid Mud Layers. I: Entrainment Model , 1997 .

[54]  Pilar García-Navarro,et al.  A numerical model for the flooding and drying of irregular domains , 2002 .

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

[56]  William Roberts,et al.  Dynamics of Estuarine Muds , 2001 .