Balanced Meshless Method for Numerical Simulation of Pollutant Transport by Shallow Water Flow over Irregular Bed: Application in the Strait of Gibraltar

• Published in: Applied sciences Vol. 12; no. 14; p. 6849
• Main Authors: Al Nuwairan, Muneerah, Chaabelasri, Elmiloud
• Format: Journal Article
• Language: English
• Published: Basel MDPI AG 01.07.2022
• Subjects: Accuracy; Advection; Advection-diffusion equation; Approximation; Coasts; Dam failure; Diffusion; Discontinuity; Exact solutions; Finite volume method; Floods; Flow velocity; hyperviscosity stability; irregular bed; Mathematical models; meshless method; Meshless methods; Methods; MUSCL schemes; Operators (mathematics); Oscillations; Partial differential equations; pollutant transport; Pollutants; Pollution dispersion; Pollution transport; Radial basis function; Rivers; Runge-Kutta method; Shallow water; shallow water flows; Simulation; Two dimensional flow; Water depth; Water flow; Water levels; Water pollution
• ISSN: 2076-3417; 2076-3417
• DOI: 10.3390/app12146849

This paper focuses on the implementation of an accurate meshless scheme for the simulation of the advection–diffusion of non-active pollutant in a two-dimensional depth-averaged flow. The depth-averaged flow model includes the shallow water system and the pollutant propagation described by the advection–diffusion equation with diffusion tensor. The mathematical model was implemented by using a meshless method based on local radial basis functions. This method was used to numerically evaluate the spatial derivatives and complemented with the second-order Runge–Kutta method for the time evolution. To remove the non-physical oscillations, which appear at the discontinuity, a filter based on a hyperviscosity operator was applied. To investigate the effectiveness and accuracy of the proposed scheme, a number of tests are presented, including the dam break problem, and the pure transport of a pollutant in a long channel. Finally, a hypothetical example of a pollutant transport in the Strait of Gibraltar is modeled. The results obtained are compared both with analytical solutions and with simulation results obtained by a finite volume method based on the Roe-MUSCL scheme. The main advantages of the proposed method are: (i) the simplicity of implementation, (ii) the ability to handle calculations of slowly varying flows or concentrations, as well as rapidly varying flows containing shocks or discontinuities, and (iii) the ability to satisfy the C-property and guarantee positive values of both water level and pollutant concentration.