Time domain particle tracking methods for simulating transport with retention and first‐order transformation

[1] Particle tracking in the time domain has received increasing attention as a technique for robustly simulating transport along one-dimensional subsurface pathways. Using a stochastic Lagrangian perspective, integral representations of transport including the effects of advection, longitudinal dispersion, and a broad class of retention models are derived; Monte Carlo sampling of that integral leads directly to new time domain particle tracking algorithms that represent a wide range of physical phenomena. Retention-time distributions are compiled for key retention models. An extension to accommodate linear transformations such as decay chains is also introduced. Detailed testing using first-order decay chains and four retention models (equilibrium sorption, limited diffusion, unlimited diffusion, and first-order kinetic sorption) demonstrate that the method is highly accurate. Simulations using flow fields produced by large-scale discrete-fracture network simulations, a transport problem that is difficult for conventional algorithms, demonstrate that the new algorithms are robust and highly efficient.

[1]  Jacques Bodin,et al.  Simulation of solute transport in discrete fracture networks using the time domain random walk method , 2003 .

[2]  A. Zuber,et al.  On the physical meaning of the dispersion equation and its solutions for different initial and boundary conditions , 1978 .

[3]  Peter C. Lichtner,et al.  Continuum formulation of multicomponent-multiphase reactive transport , 1996 .

[4]  G. Dagan Flow and transport in porous formations , 1989 .

[5]  B. Silverman Density estimation for statistics and data analysis , 1986 .

[6]  Scott L. Painter,et al.  Upscaling discrete fracture network simulations: An alternative to continuum transport models , 2005 .

[7]  Gedeon Dagan,et al.  Reactive transport and immiscible flow in geological media. I. General theory , 1996, Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[8]  V. Cvetkovic,et al.  Time-Domain Random-Walk Algorithms for Simulating Radionuclide Transport in Fractured Porous Rock , 2008 .

[9]  G. Dagan,et al.  Contaminant transport in aquifers with spatially variable hydraulic and sorption properties , 1998, Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[10]  Frederick Delay,et al.  Time domain random walk method to simulate transport by advection‐dispersion and matrix diffusion in fracture networks , 2001 .

[11]  Scott L. Painter,et al.  Stochastic simulation of radionuclide migration in discretely fractured rock near the Äspö Hard Rock Laboratory , 2004 .

[12]  S. James,et al.  Determining the random time step in a constant spatial step particle tracking algorithm , 2002 .

[13]  G. Dagan Stochastic modeling of groundwater flow by unconditional and conditional probabilities: 2. The solute transport , 1982 .

[14]  G. Dagan Solute transport in heterogeneous porous formations , 1984, Journal of Fluid Mechanics.

[15]  Gedeon Dagan,et al.  Transport of kinetically sorbing solute by steady random velocity in heterogeneous porous formations , 1994, Journal of Fluid Mechanics.

[16]  Vladimir Cvetkovic,et al.  Transport of reactive tracers in rock fractures , 1999, Journal of Fluid Mechanics.

[17]  G. Dagan,et al.  A first‐order analysis of solute flux statistics in aquifers: The combined effect of pore‐scale dispersion, sampling, and linear sorption kinetics , 2002 .

[18]  Ivars Neretnieks,et al.  Fluid flow and solute transport in a network of channels , 1993 .

[19]  R. Haggerty,et al.  Transport with multiple-rate exchange in disordered media. , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

[20]  Andrew F. B. Tompson,et al.  Numerical simulation of solute transport in three-dimensional, randomly heterogeneous porous media , 1990 .

[21]  H. Kimura,et al.  Particle-tracking technique for nuclide decay chain transport in fractured porous media. , 1990 .

[22]  V. Cvetkovic,et al.  Effect of Heterogeneity on Radionuclide Retardation in the Alluvial Aquifer Near Yucca Mountain, Nevada , 2001, Ground water.

[23]  C. F. Tsang,et al.  A particle‐tracking method for advective transport in fractures with diffusion into finite matrix blocks , 2001 .

[24]  R. Banks,et al.  A solution of the differential equation of longitudinal dispersion in porous media , 1961 .

[25]  Jacques Villermaux,et al.  Chemical engineering approach to dynamic modelling of linear chromatography: a flexible method for representing complex phenomena from simple concepts , 1987 .

[26]  Amvrossios C. Bagtzoglou,et al.  Projection functions for particle‐grid methods , 1992 .