ABSTRACT
Random walk particle tracking (RWPT) is a discrete particle method that offers several advantages for simulating solute transport in heterogeneous geological systems. The formulation is a discrete solution to the advection-dispersion equation, yielding results that are not influenced by grid-related numerical dispersion. Numerical dispersion impacts the magnitude of concentrations and gradients obtained from classical grid-based solvers in advection-dominated problems with relatively large grid Péclet numbers. Accurate predictions of concentrations are crucial for reactive transport studies, and RWPT has been recognized for its potential benefits for this kind of application. This highlights the need for a solute transport program based on RWPT that can be seamlessly integrated with industry-standard groundwater flow models. This article presents a solute transport code that implements the RWPT method by extension of the particle tracking model MODPATH, which provides the base infrastructure for interacting with several variants of MODFLOW groundwater flow models. The implementation is achieved by developing a method for determining the exact cell-exit position of a particle undergoing simultaneous advection and dispersion, allowing for the sequential transfer of particles between flow model cells. The program is compatible with rectangular unstructured grids and integrates a module for the smoothed reconstruction of concentrations. In addition, the program incorporates parallel processing of particles using the OpenMP library, enabling faster simulations of solute transport in heterogeneous systems. Numerical test cases involving different applications in hydrogeology benchmark the RWPT model with well-known transport codes.
ABSTRACT
Particle tracking has several important applications for solute transport studies in aquifer systems. Travel time distribution at observation points, particle coordinates in time and streamlines are some practical results providing information of expected transport patterns and interaction with boundary conditions. However, flow model complexity and simultaneous displacement of multiple particle groups leads to rapid increase of computational requirements. MODPATH is a particle tracking engine for MODFLOW models and source code displays potential for parallel processing of particles. This article addresses the implementation of this feature with the OpenMP library. Two synthetic aquifer applications are employed for performance tests on a desktop computer with increasing number of particles. Speed up analysis shows that dynamic thread scheduling is preferable for highly heterogeneous flows, providing processing adaptivity to the presence of slow particles. In simulations writing particles position in time, thread exclusive output files lead to higher speed up factors. Results show that above a threshold number of particles, simulation runtimes become independent of flow model grid complexity and are controlled by the large number of particles, then parallel processing reduces simulation runtimes for the particle tracking model MODPATH. Article impact statement: Parallel particles processing for MODPATH using the OpenMP library.
