A New Depth-Averaged Eulerian SPH Model for Passive Pollutant Transport in Open Channel Flows

Various nature-based solutions (NbS)—such as constructed wetlands, drainage ditches, and vegetated buffer strips—have recently demonstrated strong potential for mitigating pollutant transport in open channels and river systems. Numerical modeling is a widely adopted and effective approach for assessing the performance of these interventions. This study presents the first development of a two-dimensional (2D) meshless advection–diffusion model based on an Eulerian smoothed particle hydrodynamics (SPH) framework, specifically designed to simulate passive pollutant transport in open channel flows. The proposed model marks a pioneering application of the ESPH technique to environmental pollutant transport problems. It couples the 2D depth-averaged shallow water equations with an advection–diffusion equation to represent both fluid motion and pollutant concentration dynamics. A uniform particle arrangement ensures that each fluid particle interacts symmetrically with eight neighboring particles for flux computation. To represent the pollutant transport process, the dispersion coefficient is defined as the sum of molecular and turbulent diffusion components. The turbulent diffusion coefficient is calculated using a prescribed turbulent Schmidt number and the eddy viscosity obtained from a Smagorinsky-type mixing-length turbulence model. Three analytical case studies, including one-dimensional transcritical open channel flow, 2D isotropic and anisotropic diffusion in still water, and advection–diffusion in a 2D uniform flow, are employed to verify the model’s accuracy and convergence. The model demonstrates first-order convergence, with relative root mean square errors (RRMSEs) of approximately 0.2% for water depth and velocity, and 0.1–0.5% for concentration. Additionally, the model is applied to a laboratory experiment involving 2D pollutant dispersion in a 90° junction channel. The simulated results show good agreement with measured velocity and concentration distributions. These findings indicate that the developed model is a reliable and effective tool for evaluating the performance of NbS in mitigating pollutant transport in open channels and river systems.

Authors

Institutions

Publication Details

Journal
Water
Published
2025-07-24
DOI
https://doi.org/10.3390/w17152205
Citations
1
Primary Topic
Fluid Dynamics Simulations and Interactions
Type
article
Field-Weighted Citation Impact
1.33

Funders

Controls
|||
ALL TIME
JAN
FEB
MAR
APR
MAY
JUN
JUL
AUG
SEP
OCT
article

A New Depth-Averaged Eulerian SPH Model for Passive Pollutant Transport in Open Channel Flows

Yung‐Chieh Wang, Kaimin Shih, Kao-Hua Chang
1 citations
Water
Fluid Dynamics Simulations and Interactions
1.33
article

A New Depth-Averaged Eulerian SPH Model for Passive Pollutant Transport in Open Channel Flows

Yung‐Chieh Wang, Kaimin Shih, Kao-Hua Chang
article en
1 citations

Abstract

Various nature-based solutions (NbS)—such as constructed wetlands, drainage ditches, and vegetated buffer strips—have recently demonstrated strong potential for mitigating pollutant transport in open channels and river systems. Numerical modeling is a widely adopted and effective approach for assessing the performance of these interventions. This study presents the first development of a two-dimensional (2D) meshless advection–diffusion model based on an Eulerian smoothed particle hydrodynamics (SPH) framework, specifically designed to simulate passive pollutant transport in open channel flows. The proposed model marks a pioneering application of the ESPH technique to environmental pollutant transport problems. It couples the 2D depth-averaged shallow water equations with an advection–diffusion equation to represent both fluid motion and pollutant concentration dynamics. A uniform particle arrangement ensures that each fluid particle interacts symmetrically with eight neighboring particles for flux computation. To represent the pollutant transport process, the dispersion coefficient is defined as the sum of molecular and turbulent diffusion components. The turbulent diffusion coefficient is calculated using a prescribed turbulent Schmidt number and the eddy viscosity obtained from a Smagorinsky-type mixing-length turbulence model. Three analytical case studies, including one-dimensional transcritical open channel flow, 2D isotropic and anisotropic diffusion in still water, and advection–diffusion in a 2D uniform flow, are employed to verify the model’s accuracy and convergence. The model demonstrates first-order convergence, with relative root mean square errors (RRMSEs) of approximately 0.2% for water depth and velocity, and 0.1–0.5% for concentration. Additionally, the model is applied to a laboratory experiment involving 2D pollutant dispersion in a 90° junction channel. The simulated results show good agreement with measured velocity and concentration distributions. These findings indicate that the developed model is a reliable and effective tool for evaluating the performance of NbS in mitigating pollutant transport in open channels and river systems.

WaterVol. 17(15)
National Chung Hsing University (TW)
National Science and Technology Council
Openalex Percentile: Top 21%
Fluid Dynamics Simulations and Interactions
1.33
AI Navigator

Ask Laika to Summarize, Analyze, and Connect papers live on the map.

Summarize Papers & Methodologies

Extract key findings, datasets, and comparative methods across publications.

Benchmark Rankings & Visual Analytics

Rank top research institutions, authors, funders, topics, and journals by Field-Weighted Citation Impact (FWCI) and paper volume with instant charts.

Connect Distant Disciplines

Bridge topological clusters on the map to find hidden collaborative intersections.