A Lagrangian strategy for the numerical simulation of radionuclide transport problems

Abstract In this work a new algorithm for the computational implementation of the Forward Integral Tracking (FIT) method, originally introduced in Aquino et al. (2007a) , is presented and applied to the numerical solution of radionuclide transport problems in saturated heterogeneous porous media. The FIT is a semi-discrete numerical method which is virtually free of numerical diffusion, does not use Riemann solvers, and is computationally very efficient. A new, locally conservative procedure for the inclusion and removal of tracked points that represent the particles of the radionuclide that are carried by the fluid flow is introduced. Numerical results which indicate numerical convergence of the new procedure under mesh refinement are discussed.

[1]  George E. Karniadakis,et al.  Uncertainty quantification in simulation science , 2006, J. Comput. Phys..

[2]  Felipe Pereira,et al.  Crossover from Nonlinearity Controlled to Heterogeneity Controlled Mixing in Two-Phase Porous Media Flows , 2003 .

[3]  Yalchin Efendiev,et al.  Dynamic Data Driven Simulations in Stochastic Environments , 2006, Computing.

[4]  Philip John Binning,et al.  A forward particle tracking Eulerian-Lagrangian localized adjoint method for solution of the contaminant transport equation in three dimensions , 2002 .

[5]  A forward tracking scheme for solving radionuclide advective problems in unsaturated porous media , 2007 .

[6]  J. Bear.  Hydraulics of Groundwater , 1979 .

[7]  A. S. Francisco,et al.  An overview of Eulerian–Lagrangian schemes applied to radionuclide transport in unsaturated porous media , 2008 .

[8]  K. Ulrich Mayer,et al.  Reactive transport modeling in fractured rock: A state-of-the-science review , 2005 .

[9]  Felipe Pereira,et al.  A locally conservative Eulerian–Lagrangian numerical method and its application to nonlinear transport in porous media , 2000 .

[10]  J. Bahr,et al.  Direct comparison of kinetic and local equilibrium formulations for solute transport affected by surface reactions , 1987 .

[11]  A. Mantoglou,et al.  The Turning Bands Method for simulation of random fields using line generation by a spectral method , 1982 .

[12]  Numerical simulation of transient water infiltration in heterogeneous soils combining central schemes and mixed finite elements , 2006 .

[13]  Wonho Oh,et al.  Random field simulation and an application of kriging to image thresholding , 1998 .

[14]  Larry W. Lake,et al.  Flexible spectral methods for the generation of random fields with power-law semivariograms , 1997 .

[15]  Marcio R. Borges,et al.  Efficient generation of multi‐scale random fields: A hierarchical approach , 2010 .

[16]  J. Glimm,et al.  A random field model for anomalous diffusion in heterogeneous porous media , 1991 .