An observed error in PMPATH particle tracking algorithm for MODFLOW in case of varying porosity and a proposed correction
2012 (English)In: Computers & Geosciences, ISSN 0098-3004, E-ISSN 1873-7803, Vol. 45, 1-3 p.Article in journal (Refereed) Published
MODFLOW (Harbaugh and McDonald, 1996) is a code widely used for groundwater modeling. It has been used for numerous projects and thousands of articles have been published with results of MODFLOW, since the code was first released in 1983. Numerous codes have been developed around MODFLOW, the pre- and post processing code Processing MODFLOW (Chiang, 2005) being one example. This code includes the advective particle tracking code PMPATH (Chiang, 2005), especially convenient and commonly used for stochastic groundwater modeling analyses as it includes an easy-to-use random field generator. The particle tracking code has been regularly used to simulate the advective transport of contaminants and tracers (e.g. Stauffer, 2005; Hefting et al., 2006; Dhiman and Keshari, 2006).
A suspected error in the PMPATH code was found by us, when studying flow and particle transport in a strongly heterogeneous hydraulic conductivity fields with co-varying porosity. It emerged in the form of calculated particle flow paths converging into narrower and narrower channels. After many attempts to understand the results turned out to be unsuccessful, we decided to check into details of the code and discovered a significant but non-obvious error. This error is discussed in this paper and a procedure to obtain the correct flow paths and travel times is proposed.
Place, publisher, year, edition, pages
Elsevier, 2012. Vol. 45, 1-3 p.
PMPATH, MODFLOW, Particle tracking, Groundwater, Solute transport, Channeling
Oceanography, Hydrology, Water Resources
Research subject Hydrology
IdentifiersURN: urn:nbn:se:uu:diva-180406DOI: 10.1016/j.cageo.2012.03.001ISI: 000306771100001OAI: oai:DiVA.org:uu-180406DiVA: diva2:550042