Numerical Modeling of Large Scale Transport of Contminant Solutes Using the Global Random Walk Algorithm

Suciu N., Vamoş C., Hardelauf H., Vanderborght J. and Vereecken H.
Additional contact information
Suciu N.: 1. Romanian Academy, "T. Popoviciu" Institute of Numerical Analysis
Vamoş C.: 1. Romanian Academy, "T. Popoviciu" Institute of Numerical Analysis
Hardelauf H.: 2. Research Center Jülich (Germany), ICG-IV: Institut of Agrosphere
Vanderborght J.: 2. Research Center Jülich (Germany), ICG-IV: Institut of Agrosphere
Vereecken H.: 2. Research Center Jülich (Germany), ICG-IV: Institut of Agrosphere

Monte Carlo Methods and Applications, 2004, vol. 10, issue 2, 153-177

Abstract: The methods which track particles to simulate concentrations are successfully used in problems for large scale transport in groundwater, mainly when the aquifer properties are spatially heterogeneous and the process is advection dominated. These methods, sometimes called "analogical Monte Carlo methods", are not concerned with the numerical diffusion occurring in finite difference/element schemes. The limitations of classical random walk methods are due to large computation time and memory necessary to achieve statistically reliable results and accurate concentration fields. To overcome these computational limitations a "global random walk" (GRW) algorithm was developed. Unlike in the usual approach where the trajectory of each particle is simulated and stored, in GRW all the particles from a given grid node are scattered, at a given time, using a single numerical procedure, based on Bernoulli distribution, partial-deterministic, or deterministic rules. Because no restrictions are imposed for the maximum number of particles to be used in a simulation, the Monte Carlo repetitions are no longer necessary to achieve the convergence. It was proved that for simple diffusion problems GRW converges to the finite difference scheme and that for large scale transport problems in groundwater, GRW produces stable and statistically reliable simulations. A 2-dimensional transport problem was modeled by simulating local diffusion processes in realizations of a random velocity field. The behavior over 5000 days of the effective diffusion coefficients, concentration field and concentration fluctuations were investigated using 2560 realizations of the velocity field and 1010 particles in every realization. The results confirm the order of magnitude of the effective diffusion coefficients predicted by stochastic theory but the time needed to reach the asymptotic regime was found to be thousands times larger. It is also underlined that the concentration fluctuations and the dilution of contaminant solute depend essentially on local diffusion and boundary conditions.

Date: 2004
References: View complete reference list from CitEc
Citations:

Downloads: (external link)
https://doi.org/10.1515/156939604777303235 (text/html)
For access to full text, subscription to the journal or payment for the individual article is required.

Related works:
This item may be available elsewhere in EconPapers: Search for items with the same title.

Export reference: BibTeX RIS (EndNote, ProCite, RefMan) HTML/Text

Persistent link: https://EconPapers.repec.org/RePEc:bpj:mcmeap:v:10:y:2004:i:2:p:153-177:n:3

Ordering information: This journal article can be ordered from
https://www.degruyter.com/journal/key/mcma/html

DOI: 10.1515/156939604777303235

Access Statistics for this article

Monte Carlo Methods and Applications is currently edited by Karl K. Sabelfeld

More articles in Monte Carlo Methods and Applications from De Gruyter
Bibliographic data for series maintained by Peter Golla ().