Modeling diffusion processes in the presence of a diffuse layer at charged mineral surfaces: a benchmark exercise.
Saved in:
| Title: | Modeling diffusion processes in the presence of a diffuse layer at charged mineral surfaces: a benchmark exercise. |
|---|---|
| Authors: | Tournassat, Christophe1,2,3 (AUTHOR) c.tournassat@brgm.fr, Steefel, Carl I.1 (AUTHOR) |
| Source: | Computational Geosciences. Aug2021, Vol. 25 Issue 4, p1319-1336. 18p. |
| Subjects: | Diffusion processes, Nernst-Planck equation, Radioactive waste disposal, Electric potential, Finite differences, Isometric exercise, Radioactive wastes |
| Abstract: | The electrostatic properties of clay mineral surfaces play a significant role in their diffusion properties. The negative electrostatic potential field at clay mineral surfaces results in the presence of a diffuse layer that balances the mineral surface charge. The diffusion properties of the porosity fraction that is affected by this phenomenon are different from the diffusion properties of electroneutral bulk water. These properties have attracted growing interest from diverse communities in the past years, especially in the field of study of radioactive waste disposal. The influence of the diffuse layer can be described at the continuum scale by a set of equations that are formulated in terms of the Nernst-Planck equation. The number of codes that can handle the coupling between transport properties in clay affected by the presence of a diffuse layer in the porosity and chemical reactions is very limited, and no benchmark exercises have been published yet that make it possible to validate the numerical implementation of these equations in reactive transport codes. The present study proposes a set of benchmark exercises of increasing complexity that highlight caveats related to the finite difference (volume) treatment of the Nernst-Planck equation in the presence of a diffuse layer in heterogeneous systems. Once these problems are identified and solved, the codes PHREEQC, CrunchClay, and a new Fortran routine written for this study gave results in very good agreement for most of the benchmark exercises. When present, the differences in results were directly traceable to the differences in averaging methods at grid cell boundaries, and to the consideration or the omission of the activity gradient term in the Nernst-Planck equation. [ABSTRACT FROM AUTHOR] |
| Copyright of Computational Geosciences is the property of Springer Nature and its content may not be copied or emailed to multiple sites without the copyright holder's express written permission. Additionally, content may not be used with any artificial intelligence tools or machine learning technologies. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.) | |
| Database: | Engineering Source |
Be the first to leave a comment!
