Skip to Main content Skip to Navigation
Journal articles

Solving the Nernst‐Planck Equation in Heterogeneous Porous Media with Finite Volume Methods: Averaging Approaches at Interfaces

Abstract : Molecular diffusion of dissolved species is a fundamental mass transport process affecting many environmental and technical processes. Whereas diffusive transport of single tracers can be described by Fick's law, a multicomponent approach based on the Nernst‐Planck equation is required for charge‐coupled transport of ions. The numerical solution of the Nernst‐Planck equation requires special attention with regard to properties that are required at interfaces of numerical cells when using a finite difference or finite volume method. Weighted arithmetic and harmonic averages are used in most codes that can solve the Nernst‐Planck equation. This way of averaging is correct for diffusion coefficients, but inappropriate for solute concentrations at interfaces. This averaging approach leads to charge balance problems and thus to numerical instabilities near interfaces separating grid volumes with contrasting properties. We argue that a logarithmic‐differential average should be used. Here this result is generalized, and it is demonstrated that it generally leads to improved numerical stability and accuracy of concentrations computed near material interfaces. It is particularly relevant when modeling semi‐permeable clay membranes or membranes used in water treatment processes.
Document type :
Journal articles
Complete list of metadatas

Cited literature [52 references]  Display  Hide  Download

https://hal-insu.archives-ouvertes.fr/insu-02496420
Contributor : Nathalie Pothier <>
Submitted on : Tuesday, March 17, 2020 - 8:09:07 AM
Last modification on : Wednesday, September 23, 2020 - 8:10:02 AM

File

Tournassat_et_al-2020-Water_Re...
Publisher files allowed on an open archive

Licence


Distributed under a Creative Commons Attribution - NonCommercial 4.0 International License

Identifiers

Collections

Citation

Christophe Tournassat, Carl Steefel, Thomas Gimmi. Solving the Nernst‐Planck Equation in Heterogeneous Porous Media with Finite Volume Methods: Averaging Approaches at Interfaces. Water Resources Research, American Geophysical Union, 2020, 56 (3), 10 p. ⟨10.1029/2019WR026832⟩. ⟨insu-02496420⟩

Share

Metrics

Record views

154

Files downloads

271