A Lagrangian particle model on GPU for contaminant transport in groundwater

• Published in: Computational particle mechanics Vol. 10; no. 3; pp. 587 - 601
• Main Authors: Hou, Qingzhi, Miao, Chunfu, Chen, Shaokang, Sun, Zewei, Karemat, Alireza
• Format: Journal Article
• Language: English
• Published: Cham Springer International Publishing 01.06.2023; Springer Nature B.V
• Subjects: Advection-diffusion equation; Classical and Continuum Physics; Computational Science and Engineering; Computer architecture; Conservation laws; Contaminants; Cost analysis; Engineering; Exact solutions; Graphics processing units; Groundwater; Mathematical models; Particle methods (mathematics); Pollution transport; Smooth particle hydrodynamics; Soil contamination; Theoretical and Applied Mechanics; Contaminant transport; SPH; Groundwater; GPU; Lagrangian particle model
• ISSN: 2196-4378; 2196-4386
• DOI: 10.1007/s40571-022-00495-5

To simulate contaminant transport in groundwater, this paper proposes a parallelized Lagrangian particle model using compute unified device architecture (CUDA) on graphics processing unit (GPU) based on smoothed particle hydrodynamics (SPH) method. The solved governing equation is the advection–diffusion equations (ADEs) with retardation factor for given typical flow fields. To solve the inherent particle inconsistency problem of traditional SPH method, the corrective smoothed particle method (CSPM) is applied. The speedup ratio of the parallelized SPH solver for ADEs is analyzed. The consistency and convergence of the proposed model are theoretically analyzed and numerically tested. The reduction of its computational cost and storage requirement is discussed. Numerical examples including one-dimensional (1D) and two-dimensional (2D) cases are simulated, and the results are compared with the analytical solutions and those obtained by the high-resolution monotonic upstream schemes for conservation laws (MUSCL) scheme. To further verify the practicality of the model, two engineering cases of contaminant transport through soil into groundwater are investigated. It is shown that the solutions of the developed model are in good agreement with measured data.