Multispecies diffusion models: A study of uranyl species diffusion

Abstract

[1] Rigorous numerical description of multispecies 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 in imposing charge conservation. In this study, several diffusion models with variable complexity in charge and species coupling were formulated and compared to describe reactive multispecies diffusion in groundwater. Diffusion of uranyl [U(VI)] species was used as an example in demonstrating the effectiveness of the models in describing multispecies diffusion. Numerical simulations found that a diffusion model with a single, common diffusion coefficient for all species was sufficient to describe multispecies U(VI) diffusion under a steady state condition of major chemical composition, but not under transient chemical conditions. Simulations revealed that for multispecies U(VI) diffusion under transient chemical conditions, a fully coupled diffusion model could be well approximated by a component-based diffusion model when the diffusion coefficient for each chemical component was properly selected. The component-based diffusion model considers the 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 enforced, if necessary, by adding a secondary migration term resulting from model simplification. The effect of ion activity coefficient gradients on multispecies diffusion is also discussed. The diffusion models were applied to describe U(VI) diffusive mass transfer in intragranular domains in two sediments collected from U.S. Department of Energy's Hanford 300A, where intragranular 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 had been previously described using a semiempirical, multirate model. Compared with the multirate model, the diffusion models have the advantage to provide spatiotemporal speciation evolution within the diffusion domains.Citation: Liu, C., J. Shang, and J. M. Zachara (2011), Multispecies diffusion models: A study of uranyl species diffusion, Water Resour.

Pore-Scale Process Coupling and Effective Surface Reaction Rates in Heterogeneous Subsurface Materials

“…Various numerical approaches exist for solving the multi-species diffusion problem described by Equation (14) using implicit or explicit numerical scheme, or their combination (Yeh and Tripathi 1991;Press et al 1992;Liu et al 2011). …”

Lattice Boltzmann-Based Approaches for Pore-Scale Reactive Transport

“…The multispecies diffusive fl uxes including the electrochemical migration can be important in dilute solutions (Van Cappellen and Gaillard 1996;Li et al 2007;Lichtner and Kang 2007;Liu et al 2011;Molins et al 2012;Zhang and Klapper 2014). In this article, we will briefl y introduce the lattice Poisson-Boltzmann method to account for the electroosmotic fl ows as a special case.…”

Benchmark reactive transport simulations of a column experiment in compacted bentonite with multispecies diffusion and explicit treatment of electrostatic effects

“…The simplification of using a single diffusion coefficient offers the advantage that local charge balance is guaranteed and charge-induced electrostatic interaction is avoided which simplifies numerical execution. The disadvantage is that species properties such as ion size, charge, and hydration energy which may affect species-specific diffusivities are ignored [17]. Rigorous modelling of diffusion in systems containing charged species requires electrochemical migration to be included as a transport process in such a way that positively and negatively charged species move in a coordinated manner to maintain local charge balance.…”