Modeling transport of charged species in pore networks: Solution of the Nernst-Planck equations coupled with fluid flow and charge conservation equations
Abstract
A pore network modeling (PNM) framework
The proposed framework is developed by first deriving the numerical model equations (NMEs) corresponding to the partial differential equations (PDEs) based on several different time and space discretization schemes, which are compared to assess solutions accuracy. The derivation also considers various charge conservation scenarios, which also have pros and cons in terms of speed and accuracy. Ion transport problems in arbitrary pore networks were considered and solved using both PNM and finite element method (FEM) solvers. Comparisons showed an average deviation, in terms of ions concentration, between PNM and FEM below 5% with the PNM simulations being over 104 times faster than the FEM ones for a medium including about 104 pores. The improved accuracy is achieved by utilizing more accurate discretization schemes for both the advective and migrative terms, adopted from the CFD literature. The NMEs were implemented within the open-source package OpenPNM based on the iterative Gummel algorithm with relaxation. This work presents a comprehensive approach to modeling charged species transport suitable for a wide range of applications from electrochemical devices to nanoparticle movement in the subsurface.- Publication:
-
Computers and Geosciences
- Pub Date:
- July 2020
- DOI:
- 10.1016/j.cageo.2020.104505
- arXiv:
- arXiv:2001.04881
- Bibcode:
- 2020CG....14004505A
- Keywords:
-
- Porous media;
- Simulations;
- Nernst-Planck equations;
- Pore network modeling;
- openPNM;
- Physics - Computational Physics;
- Physics - Applied Physics
- E-Print:
- doi:10.1016/j.cageo.2020.104505