Brownian dynamics simulations were performed to examine the self-diffusion of cationic (Na+), anionic (Cl-), and water tracers in a saturated compacted clay at the micrometric scale. Clay minerals particles, representing stacks of montmorillonite layers, are modeled as cylindrical platelets with characteristics that depend on the dry bulk density of the porous medium, ρb. The tracers are allowed to diffuse between the platelets with their diffusion coefficient in bulk aqueous solution as well as inside the platelets with a diffusion coefficient that depends on the hydration level of the interlayer pores as determined from experiments or molecular dynamics simulations. The probability for a tracer to enter a platelet is related to its charge as predicted by the Poisson-Boltzmann equation. The apparent diffusion coefficients are calculated as a function of ρb for various salinities and compared with the values obtained from tracer diffusion experiments available in the literature. Despite the simplicity of the clay model, the simulations predict most features of the experimental data and keep the door open for further improvements.
All Science Journal Classification (ASJC) codes
- Geochemistry and Petrology
- Brownian dynamics
- Molecular dynamics