A Mass-transfer Particle-tracking Method for Simulating Transport with Discontinuous Diffusion Coefficients
Abstract
The problem of a spatially discontinuous diffusion coefficient () is one that may be encountered in hydrogeologic systems due to natural geological features or as a consequence of numerical discretization of flow properties. To date, mass-transfer particle-tracking (MTPT) methods, a family of Lagrangian methods in which diffusion is jointly simulated by random walk and diffusive mass transfers, have been unable to solve this problem. This manuscript presents a new mass-transfer (MT) algorithm that enables MTPT methods to accurately solve the problem of discontinuous . To achieve this, we derive a semi-analytical solution to the discontinuous problem by employing a predictor-corrector approach, and we use this semi-analytical solution as the weighting function in a reformulated MT algorithm. This semi-analytical solution is generalized for cases with multiple 1D interfaces as well as for 2D cases, including a tiling of 4 subdomains that corresponds to a numerically-generated diffusion field. The solutions generated by this new mass-transfer algorithm closely agree with an analytical 1D solution or, in more complicated cases, trusted numerical results, demonstrating the success of our proposed approach.
Cite
@article{arxiv.2001.04936,
title = {A Mass-transfer Particle-tracking Method for Simulating Transport with Discontinuous Diffusion Coefficients},
author = {Michael J. Schmidt and Nicholas B. Engdahl and Stephen D. Pankavich and Diogo Bolster},
journal= {arXiv preprint arXiv:2001.04936},
year = {2020}
}
Comments
26 pages, 11 figures