Reactive Particle-tracking Solutions to a Benchmark Problem on Heavy Metal Cycling in Lake Sediments
arXiv:1908.09818 · doi:10.1016/j.jconhyd.2020.103642
Abstract
Geochemical systems are known to exhibit highly variable spatiotemporal behavior. This may be observed both in non-smooth concentration curves in space for a single sampling time and also in variability between samples taken from the same location at different times. However, most models that are designed to simulate these systems provide only single-solution smooth curves and fail to capture the noise and variability seen in the data. We apply a recently developed reactive particle-tracking method to a system that displays highly-complex geochemical behavior. When the method is made to most closely resemble a corresponding Eulerian method, in its unperturbed form, we see near-exact match between solutions of the two models. More importantly, we consider two approaches for perturbing the model and find that the spatially-perturbed condition is able to capture a greater degree of the variability present in the data. This method of perturbation is a task to which particle methods are uniquely suited and Eulerian models are not well-suited. Additionally, because of the nature of the algorithm, noisy spatial gradients can be highly resolved by a large number of mobile particles, and this incurs negligible computational cost, as compared to expensive chemistry calculations.
29 pages, 8 figures
References in corpus (12)
- The SIESTA method for ab initio order-N materials simulation
- Efficient index handling of multidimensional periodic boundary conditions
- Linear-scaling ab-initio calculations for large and complex systems
- C20: Fulleren, Bowl or Ring? New Results from Coupled-Cluster Calculations
- On the separate treatment of mixing and spreading by the reactive-particle-tracking algorithm: An example of accurate upscaling of reactive Poiseuille flow
- A Kernel-based Lagrangian Method for Imperfectly-mixed Chemical Reactions
- Accelerating and parallelizing Lagrangian simulations of mixing-limited reactive transport
- Particle Density Estimation with Grid-Projected Adaptive Kernels
- On the accuracy of simulating mixing by random-walk particle-based mass-transfer algorithms
- Numerical Equivalence Between SPH and Probabilistic Mass Transfer Methods for Lagrangian Simulation of Dispersion
- Pictures of Blockscale Transport: Effective versus Ensemble Dispersion and Its Uncertainty
- Estimates of Effective Hubbard Model Parameters for C20 isomers