Simulating Donnan equilibria based on the Nernst-Planck equation

Gimmi, Thomas; Alt-Epping, Peter (2018). Simulating Donnan equilibria based on the Nernst-Planck equation. Geochimica et cosmochimica acta, 232, pp. 1-13. Elsevier Science 10.1016/j.gca.2018.04.003

Preview	Text GimmiAltEpping_corr_inpressGCA_2018.pdf - Accepted Version Available under License Creative Commons: Attribution-Noncommercial-No Derivative Works (CC-BY-NC-ND). Download (2MB) | Preview
	Text Gimmi_AltEpping_GCA2018.pdf - Published Version Restricted to registered users only Available under License Publisher holds Copyright. Download (955kB) | Request a copy

Understanding ion transport through clays and clay membranes is important for many geochemical and environmental applications. Ion transport is affected by electrostatic forces exerted by charged clay surfaces. Anions are partly excluded from pore water near these surfaces, whereas cations are enriched. Such effects can be modeled by the Donnan approach. Here we introduce a new, comparatively simple way to represent Donnan equilibria in transport simulations. We include charged surfaces as immobile ions in the balance equation and calculate coupled transport of all components, including the immobile charges, with the Nernst-Planck equation. This results in an additional diffusion potential that influences ion transport, leading to Donnan ion distributions while maintaining local charge balance. The validity of our new approach was demonstrated by comparing Nernst-Planck simulations using the reactive transport code Flotran with analytical solutions available for simple Donnan systems. Attention has to be paid to the numerical evaluation of the electrochemical migration term in the Nernst-Planck equation to obtain correct results for asymmetric electrolytes. Sensitivity simulations demonstrate the influence of various Donnan model parameters on simulated anion accessible porosities. It is furthermore shown that the salt diffusion coefficient in a Donnan pore depends on local concentrations, in contrast to the aqueous salt diffusion coefficient. Our approach can be easily implemented into other transport codes. It is versatile and facilitates, for instance, assessing the implications of different activity models for the Donnan porosity.

Interest & Impact

Downloads

345 since deposited on 24 May 2018
138 in the past 12 months

Citations

5 Citations in Web of Science ®
6 Citations in Scopus

Search

in Google Scholar™

Services

Actions (login required)

Edit item

Actions (login required)

Edit item