Loughborough University
A graphics processing unit based robust numerical model for solute transport driven by torrential flow condition.pdf (3.19 MB)

A graphics processing unit-based robust numerical model for solute transport driven by torrential flow condition

Download (3.19 MB)
journal contribution
posted on 2021-12-01, 10:33 authored by Jing-ming Hou, Bao-shan Shi, Qiuhua LiangQiuhua Liang, Yu Tong, Yong-de Kang, Zhao-an Zhang, Gang-gang Bai, Xu-jun Gao, Xiao Yang
Solute transport simulations are important in water pollution events. This paper introduces a finite volume Godunov-type model for solving a 4×4 matrix form of the hyperbolic conservation laws consisting of 2D shallow water equations and transport equations. The model adopts the Harten-Lax-van Leer-contact (HLLC)-approximate Riemann solution to calculate the cell interface fluxes. It can deal well with the changes in the dry and wet interfaces in an actual complex terrain, and it has a strong shock-wave capturing ability. Using monotonic upstream-centred scheme for conservation laws (MUSCL) linear reconstruction with finite slope and the Runge-Kutta time integration method can achieve second-order accuracy. At the same time, the introduction of graphics processing unit (GPU)-accelerated computing technology greatly increases the computing speed. The model is validated against multiple benchmarks, and the results are in good agreement with analytical solutions and other published numerical predictions. The third test case uses the GPU and central processing unit (CPU) calculation models which take 3.865 s and 13.865 s, respectively, indicating that the GPU calculation model can increase the calculation speed by 3.6 times. In the fourth test case, comparing the numerical model calculated by GPU with the traditional numerical model calculated by CPU, the calculation efficiencies of the numerical model calculated by GPU under different resolution grids are 9.8–44.6 times higher than those by CPU. Therefore, it has better potential than previous models for large-scale simulation of solute transport in water pollution incidents. It can provide a reliable theoretical basis and strong data support in the rapid assessment and early warning of water pollution accidents.


National Natural Science Foundation of China (Nos. 52009104 and 52079106)

Shaanxi Provincial Department of Water Resources Project (No. 2017slkj-14)

Shaanxi Provincial Department of Science and Technology Project (No. 2017JQ3043), China



  • Architecture, Building and Civil Engineering

Published in

Journal of Zhejiang University Science A






835 - 850




  • AM (Accepted Manuscript)

Rights holder

© Zhejiang University Press

Publisher statement

This paper was accepted for publication in the journal Journal of Zhejiang University Science A and the definitive published version is available at https://doi.org/10.1631/jzus.A2000585.

Acceptance date


Publication date


Copyright date







  • en


Prof Qiuhua Liang. Deposit date: 29 November 2021