Premium
Time domain particle tracking methods for simulating transport with retention and first‐order transformation
Author(s) -
Painter Scott,
Cvetkovic Vladimir,
Mancillas James,
Pensado Osvaldo
Publication year - 2008
Publication title -
water resources research
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.863
H-Index - 217
eISSN - 1944-7973
pISSN - 0043-1397
DOI - 10.1029/2007wr005944
Subject(s) - tracking (education) , monte carlo method , statistical physics , diffusion , computer science , advection , domain (mathematical analysis) , particle (ecology) , range (aeronautics) , transformation (genetics) , mathematical optimization , algorithm , mechanics , mathematics , physics , materials science , mathematical analysis , geology , chemistry , psychology , pedagogy , biochemistry , statistics , oceanography , gene , composite material , thermodynamics
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.