The 3rd USGS Modeling Conference (7-11 June 2010)

Paper No. 9
Presentation Time: 11:25 AM

CONVOLUTION-BASED PARTICLE TRACKING METHOD FOR MODELING GROUNDWATER TRANSPORT FOR TRANSIENT FLOW CONDITIONS


SRINIVANSAN, G.1, KEATING, E.2, LI, C.3, KANG, Q.4 and KWICKLIS, E.4, (1)Los Alamos National Laboratory, Theoretical Division-Mathematical Modeling and Analysis and Plasma, MS B284, PO box 1663, Los Alamos, NM 87544, (2)Los Alamos National Laboratory, Earth and Environmental Sciences, MS T003, Los Alamos, NM 87544, (3)Areva Federal Services, Albuquerque, NM 87110, (4)Los Alamos National Laboratory, Earth and Environmental Sciences, MS T003, PO Box 1663, Los Alamos, NM 87544, gowri@lanl.gov

A Convolution-Based Particle Tracking (CBPT) method is developed to calculate solute concentrations in groundwater models under transient flow conditions. This method is valid for linear transport processes such as first order decay or linear sorption. A pulse of particles is introduced at various times for every source location, depending on the nature of the flow. Numerical convolution of particle paths obtained at each release time and location with a time-varying source term is performed to yield the shape of the plume. The assumption of linearity allows for superposition of multiple solute sources. The CBPT method uses fewer particles for a single particle tracking run, from which source term variability, sorption and decay can all be simulated rapidly without additional process-model runs. We demonstrate the computational efficiency and numerical accuracy of the CBPT method for transient flow, for various flow and transport scenarios with special emphasis on diffusion, in a 3D porous medium.