A Riemann solver for unsteady computation of 2D shallow flows with variable density

A novel 2D numerical model for vertically homogeneous shallow flows with variable horizontal density is presented. Density varies according to the volumetric concentration of different components or species that can represent suspended material or dissolved solutes. The system of equations is formed by the 2D equations for mass and momentum of the mixture, supplemented by equations for the mass or volume fraction of the mixture constituents. A new formulation of the Roe-type scheme including density variation is defined to solve the system on two-dimensional meshes. By using an augmented Riemann solver, the numerical scheme is defined properly including the presence of source terms involving reaction. The numerical scheme is validated using analytical steady-state solutions of variable-density flows and exact solutions for the particular case of initial value Riemann problems with variable bed level and reaction terms. Also, a 2D case that includes interaction with obstacles illustrates the stability and robustness of the numerical scheme in presence of non-uniform bed topography and wetting/drying fronts. The obtained results point out that the new method is able to predict faithfully the overall behavior of the solution and of any type of waves.

[1]  V. Guinot.  Approximate Riemann Solvers , 2010 .

[2]  Sheng Jin,et al.  A Conservative Coupled Flow/Transport Model with Zero Mass Error , 2009 .

[3]  János Józsa,et al.  Variable density bore interaction with block obstacles , 2011 .

[4]  Javier Murillo,et al.  Generalized Roe schemes for 1D two-phase, free-surface flows over a mobile bed , 2008, J. Comput. Phys..

[5]  P ? ? ? ? ? ? ? % ? ? ? ? , 1991 .

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

[7]  Javier Murillo,et al.  Coupling between shallow water and solute flow equations: analysis and management of source terms in 2D , 2005 .

[8]  Javier Murillo,et al.  Weak solutions for partial differential equations with source terms: Application to the shallow water equations , 2010, J. Comput. Phys..

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

[10]  Javier Murillo,et al.  Time step restrictions for well‐balanced shallow water solutions in non‐zero velocity steady states , 2009 .

[11]  Javier Murillo,et al.  Improved Riemann solvers for complex transport in two-dimensional unsteady shallow flow , 2011, J. Comput. Phys..

[12]  Javier Murillo,et al.  An efficient and conservative model for solute transport in unsteady shallow water flow , 2009 .

[13]  R. LeVeque.  Finite Volume Methods for Hyperbolic Problems: Characteristics and Riemann Problems for Linear Hyperbolic Equations , 2002 .

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

[15]  Philip L. Roe,et al.  Efficient construction and utilisation of approximate riemann solutions , 1985 .

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

[17]  Ami Harten,et al.  Self adjusting grid methods for one-dimensional hyperbolic conservation laws☆ , 1983 .

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

[19]  Fayssal Benkhaldoun,et al.  Exact solutions to the Riemann problem of the shallow water equations with a bottom step , 2001 .

[20]  Javier Murillo,et al.  Conservative numerical simulation of multi-component transport in two-dimensional unsteady shallow water flow , 2009, J. Comput. Phys..

[21]  Hansong Tang,et al.  A model-coupling framework for nearshore waves, currents, sediment transport, and seabed morphology , 2009 .

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

[23]  Javier Murillo,et al.  Preserving bounded and conservative solutions of transport in one‐dimensional shallow‐water flow with upwind numerical schemes: Application to fertigation and solute transport in rivers , 2008 .

[24]  David L. George,et al.  Augmented Riemann solvers for the shallow water equations over variable topography with steady states and inundation , 2008, J. Comput. Phys..

[25]  Javier Murillo,et al.  Analysis of a second‐order upwind method for the simulation of solute transport in 2D shallow water flow , 2008 .

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

[27]  R. Courant,et al.  On the solution of nonlinear hyperbolic differential equations by finite differences , 1952 .

[28]  A. Borthwick,et al.  1‐D numerical modelling of shallow flows with variable horizontal density , 2009 .

[29]  Javier Murillo,et al.  An Exner-based coupled model for two-dimensional transient flow over erodible bed , 2010, J. Comput. Phys..

[30]  Eleuterio F. Toro,et al.  Exact solution of the Riemann problem for the shallow water equations with discontinuous bottom geometry , 2008, J. Comput. Phys..

[31]  P. Roe.  Approximate Riemann Solvers, Parameter Vectors, and Difference Schemes , 1997 .