Well-balanced finite volume schemes for pollutant transport by shallow water equations on unstructured meshes

Pollutant transport by shallow water flows on non-flat topography is presented and numerically solved using a finite volume scheme. The method uses unstructured meshes, incorporates upwinded numerical fluxes and slope limiters to provide sharp resolution of steep bathymetric gradients that may form in the approximate solution. The scheme is non-oscillatory and possesses conservation property that conserves the pollutant mass during the transport process. Numerical results are presented for three test examples which demonstrate the accuracy and robustness of the scheme and its applicability in predicting pollutant transport by shallow water flows. In this paper, we also apply the developed scheme for a pollutant transport event in the Strait of Gibraltar. The scheme is efficient, robust and may be used for practical pollutant transport phenomena.

[1]  U. C. Kothyari,et al.  Time-line cubic spline interpolation scheme for solution of advection equation , 2001 .

[2]  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..

[3]  Fayssal Benkhaldoun,et al.  A Non Homogeneous Riemann Solver for shallow water and two phase flows , 2006 .

[4]  J. G. Lafuente,et al.  Sea Level in the Strait of Gibraltar: Tides , 1990 .

[5]  F Benkhaldoun,et al.  Ignition of fuel issuing from a porous cylinder located adjacent to a heated wall: A numerical study , 2004 .

[6]  Rémi Abgrall,et al.  Efficient numerical approximation of compressible multi-material flow for unstructured meshes , 2003 .

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

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

[9]  J. Greenberg,et al.  A well-balanced scheme for the numerical processing of source terms in hyperbolic equations , 1996 .

[10]  Emmanuel Audusse,et al.  A Fast and Stable Well-Balanced Scheme with Hydrostatic Reconstruction for Shallow Water Flows , 2004, SIAM J. Sci. Comput..

[11]  F. Benkhaldoun,et al.  Positivity preserving finite volume Roe: schemes for transport-diffusion equations , 1999 .

[12]  Ming-Hseng Tseng.  The improved surface gradient method for flows simulation in variable bed topography channel using TVD‐MacCormack scheme , 2003 .

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

[14]  Luka Sopta,et al.  ENO and WENO Schemes with the Exact Conservation Property for One-Dimensional Shallow Water Equations , 2002 .

[15]  I. Babuska,et al.  A-posteriori estimation of the error in the error estimate , 1998 .

[16]  I. Kane.  A numerical investigation into the solution of the equations describing baroclinic-barotropic instability in the atmosphere , 1976 .

[17]  Thierry Gallouët,et al.  On an Approximate Godunov Scheme , 1999 .

[18]  K. Ohgushi,et al.  Refined Numerical Scheme for Advective Transport in Diffusion Simulation , 1997 .

[19]  Mohammed Seaïd,et al.  Non‐oscillatory relaxation methods for the shallow‐water equations in one and two space dimensions , 2004 .

[20]  Mourad Heniche,et al.  A two-dimensional finite element drying-wetting shallow water model for rivers and estuaries , 2000 .

[21]  Finite Element Modified Method of Characteristics for Shallow Water Flows: Application to the Strait of Gibraltar , 2006 .

[22]  G. D. van Albada,et al.  A comparative study of computational methods in cosmic gas dynamics , 1982 .

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

[24]  R Periáñez.  A particle-tracking model for simulating pollutant dispersion in the Strait of Gibraltar. , 2004, Marine pollution bulletin.

[25]  Manuel Espino Infantes,et al.  Un modelo numérico en elementos finitos para la corriente inducida por la marea. Aplicaciones al Estrecho de Gibraltar , 1995 .

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

[27]  Raytcho D. Lazarov,et al.  Modeling and simulation of pollutants transport in rivers , 2007 .

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

[29]  M. S. Horritt,et al.  Development of physically based meshes for two‐dimensional models of meandering channel flow , 2000 .

[30]  F. Gómez.  The role of the exchanges through the Strait of Gibraltar on the budget of elements in the Western Mediterranean Sea: consequences of human-induced modifications. , 2003, Marine pollution bulletin.

[31]  Abioala A. Akanbi,et al.  Model for Flood Propagation on Initially Dry Land , 1988 .

[32]  Roger Alexander Falconer,et al.  Tidal Flow and Transport Modeling Using ULTIMATE QUICKEST Scheme , 1997 .

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

[34]  Fayssal Benkhaldoun,et al.  A sign matrix based scheme for non-homogeneous PDE's with an analysis of the convergence stagnation phenomenon , 2007, J. Comput. Phys..

[35]  B. Perthame,et al.  A kinetic scheme for the Saint-Venant system¶with a source term , 2001 .

[36]  Pilar García-Navarro,et al.  A HIGH-RESOLUTION GODUNOV-TYPE SCHEME IN FINITE VOLUMES FOR THE 2D SHALLOW-WATER EQUATIONS , 1993 .