PNNL

Abstract

Rigorous numerical description of multi-species diffusion requires coupling of species, charge, and aqueous and surface complexation reactions that collectively affect diffusive fluxes. The applicability of a fully coupled diffusion model is, however, often constrained by the availability of species self-diffusion coefficients, as well as by computational complication for imposing charge conservation. In this study, several diffusion models with variable complexity in charge and species coupling were formulated and compared to describe reactive multi-species diffusion in groundwater. Diffusion of uranyl [U(VI)] species was used as an example in demonstrating the effectiveness of the models in describing multi-species diffusion. Numerical simulations found that a diffusion model with a single, common diffusion coefficient for all species was sufficient to describe multi-species U(VI) diffusion under steady-state condition of major chemical composition, but not under transient chemical conditions. Simulations revealed that a fully coupled diffusion model can be well approximated by a component-based diffusion model, which considers difference in diffusion coefficients between chemical components, but not between the species within each chemical component. This treatment significantly enhanced computational efficiency at the expense of minor charge conservation. The charge balance in the component-based diffusion model can be rigorously enforced, if necessary, by adding an artificial kinetic reaction term induced by the charge separation. The diffusion models were applied to describe U(VI) diffusive mass transfer in intragranular domains in two sediments collected from US Department of Energy’s Hanford 300A where intragrain diffusion is a rate-limiting process controlling U(VI) adsorption and desorption. The grain-scale reactive diffusion model was able to describe U(VI) adsorption/desorption kinetics that has been described using a semi-empirical, multi-rate model. Compared with the multi-rate model, the diffusion models have the advantage to provide spatiotemporal speciation evolution within the diffusion domains.